<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Wed, Feb 26, 2014 at 11:43 PM, Dharmendar Reddy <span dir="ltr"><<a href="mailto:dharmareddy84@gmail.com" target="_blank">dharmareddy84@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">Hello,<br>
        I need to load a mesh file into DMPlex object. The mesh is<br>
generated by a software which outputs the mesh in this format:<br>
<br>
List of Vertex Coordinates (i=1 to N_V) :  P_i = X_i Y_i Z_i<br>
<br>
List of Edges (i= 1 to N_E):  E_i = P_m P_n where m,n <= N_V<br>
<br>
List of Faces (i=1 to N_F):  F_i = E_m E_n E_o where m,n,o < N_E<br>
(edge ids can be negative or positive depending on the orientation)<br>
<br>
List of Cells (i = 1 to N_C): C_i = F_m F_n F_o F_p .... for say Tetrahedron<br>
                                                  or<br>
                                                  = E_m E_n E_o for<br>
triangles on the boundary regions.<br></blockquote><div><br></div><div>We were cool until here. This sounds like Gmsh. This aliasing of things in</div><div>the domain is really unacceptable. The boundary faces have already been defined</div>
<div>above in the faces section. They should just be tagged.</div><div><br></div><div>If this is Gmsh, there is trial support in next, DMPlexCreateGmsh(). If not, I see a few</div><div>choices:</div><div><br></div><div>  a) Use DMPlexCreateFromCellList(). This means tracing through the info to get just</div>
<div>       the vertices on cells. However, it would mean that you get consistent ordering,</div><div>       which is not true for the faces in Gmsh files. You could renumber faces afterwards to</div><div>       conform to the original file, once all orientations have been calculated.</div>
<div><br></div><div>  b) Use DMPlexCreateFromDAG(). This gets the numbering from above, but it may</div><div>       have orientation inconsistencies in 3D.</div><div><br></div><div>With either one, you would ignore the boundary faces, and use DMPlexJoin() to locate</div>
<div>each one and tag it in a label after the mesh was built.</div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">

Any suggestions on how i can read this mesh file into DMPlex.<br>
<br>
Thanks<br>
<span class=""><font color="#888888">Reddy<br>
</font></span></blockquote></div><br><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener
</div></div>