On Tue, Dec 12, 2023 at 12:22 PM onur.notonur via petsc-users <
petsc-users@mcs.anl.gov> wrote:
> Hi,
>
> I hope this email finds you well. I am currently working on importing an
> OpenFOAM PolyMesh into DMPlex, and I've encountered an issue. The PolyMesh
> format includes face owner cells/neighbor cells and face-to-vertex
> connectivity. I was using the "DMPlexCreateFromCellListPetsc()" function,
> which required cell-to-vertex connectivity. However, when attempting to
> create the cell connectivity using an edge loop [p_0, p_1, ..., p_7] (p_n
> and p_(n+1) are valid edges in my mesh), I encountered an error stating,
> "Could not find orientation for quadrilateral."
>
> (Actually at first, I generated the connectivity list by simply creating a
> cell-to-face list and then using that to create a cell-to-vertex list.
> (just map over the list and remove duplicates) This created a DMPlex
> successfully, however, resulted in a mesh that was incorrect when looking
> with ParaView. I think that was because of I stated wrong edge loop to
> create cells)
>
> I understand that I may need to follow a different format for
> connectivity, but I'm not sure what that format is. My current mesh is
> hexahedral, consisting of 8 corner elements(if important). I would
> appreciate any guidance on a general approach to address this issue.
>
Can you start by giving the PolyMesh format, or some URL with it documented?
Thanks,
Matt
> Thank you for your time and assistance.
> Best,
> Onur
> Sent with Proton Mail secure email.
>
--
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
https://www.cse.buffalo.edu/~knepley/ <http://www.cse.buffalo.edu/~knepley/>
