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. Thank you for your time and assistance. Best, Onur Sent with Proton Mail secure email.
