<div dir="ltr">Dear all,<div><br></div><div>I am working on a vertex-centered Finite Volume <span class="" id=":4j7.1" tabindex="-1">CFD</span> solver (with only tetrahedral or triangular cells). I'd like to be able to use <span class="" id=":4j7.2" tabindex="-1">DMPlex</span>, but I have a couple of questions:</div><div><br></div><div>For me, the most convenient way to store the mesh is to have, in addition to the standard <span class="" id=":4j7.3" tabindex="-1">DMPlex</span> built from a cell to node connectivity is to also have 2 types of fictitious edges:</div><div>   1- for periodic surfaces, I build ghost cells and ghost nodes prior to building the <span class="" id=":4j7.5" tabindex="-1">DMPlex</span>, and I'd like to add "edges" that will give the matching between the ghost nodes and the corresponding nodes. These edges will not belong to any face or cell --> what would be the most convenient way to add such edges ?</div><div>    2- for each node on each boundary surface, I'd like to add a fictitious edge to easily compute the Finite Volume operators. --> can I have edges that link a node to itself (possible several times if the node belongs to several surfaces)? alternatively, should I add a fictitious node for each surface, and link all the nodes in the surface to this fictitious node?<br></div><div><br></div><div>For the mesh partition, which option should I use in <br></div><div><br></div><div>    <span class="" id=":4j7.8" tabindex="-1">ierr</span> = <span class="" id=":4j7.9" tabindex="-1">DMPlexSetAdjacencyUseCone</span>(_<span class="" id=":4j7.10" tabindex="-1">dm</span>,<span class="" id=":4j7.11" tabindex="-1">PETSC</span>_TRUE);</div><div>    <span class="" id=":4j7.12" tabindex="-1">ierr</span> = <span class="" id=":4j7.13" tabindex="-1">DMPlexSetAdjacencyUseClosure</span>(_<span class="" id=":4j7.14" tabindex="-1">dm</span>,<span class="" id=":4j7.15" tabindex="-1">PETSC</span>_TRUE);</div><div><br></div><div>knowing that the unknowns are stores at the <span class="" id=":4j7.16" tabindex="-1">vertices</span> and that two <span class="" id=":4j7.17" tabindex="-1">vertices</span> are connected if there is an edge between them?</div><div><br></div><div>Thank you very much for your library, and for your help.</div><div><br></div><div>Best regards,</div><div><br></div><div>Xavier</div></div>