<br><br><div>The mesh contains a cube. I tried to change the ordering in cone of dm. In some ordering, in Paraview, I also noticed too many interior edges and saw the interior nodes. I have not yet been able to see the interior edge correctly placed. Do you suggest to output in other format? I do not really know where this misplacing of edges comes from.</div><div>PetscInt    numPoints[2]         = {27, 8};</div><div>PetscInt    coneSize[35]        = {8,8,8,8,8,8,8,8, 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0};</div><div>PetscInt    cones[64]                = { 8, 9, 16, 15,   17, 24, 25, 18,</div><div>  9, 10,11, 16,   18, 25, 20, 19,</div><div>                 16,11, 12,13,   25, 22, 21, 20,</div><div>15,16, 13,14,   24, 23, 22, 25,    /////////// Not see the interior node and edge</div><div>17,18,25, 24,   26, 33, 34, 27,</div><div>18,19,20, 25,   27, 34, 29, 28,</div><div>        25,20,21, 22,   34, 31, 30, 29,</div><div>24,25,22, 23,   33, 32, 31, 34} ;</div><div>PetscInt    cones2[64]            = { 8, 15, 16, 9,   17, 18, 25, 24,</div><div>             9, 16,11, 10,   18, 19, 20, 25,         /////////// See the interior node and edge</div><div>     16,13, 12,11,   25, 20, 21, 22,</div><div>    15,14, 13,16,   24, 25, 22, 23,</div><div>    17,18,25, 24,   26, 27, 34, 33,</div><div>    18,19,20, 25,   27, 28, 29, 34,</div><div>    25,20,21, 22,   34, 29, 30, 31,</div><div>    24,25,22, 23,   33, 34, 31, 32} ;</div><div>PetscInt    coneOrientations[64] =  {0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0};</div><div>PetscScalar vertexCoords[81]     = {0.0,0.0,0.0,  1.0,0.0, 0.0,  2.0,0.0,0.0,  2.0,1.0,0.0,  2.0,2.0, 0.0,  1.0,2.0,0.0,</div><div>      0.0,2.0,0.0,  0.0,1.0, 0.0, 1.0,1.0,0.0,</div><div>      0.0,0.0,1.0,  1.0,0.0, 1.0,  2.0,0.0,1.0,  2.0,1.0,1.0,  2.0,2.0, 1.0,  1.0,2.0,1.0,</div><div>      0.0,2.0,1.0,  0.0,1.0, 1.0, 1.0,1.0,1.0,</div><div>      0.0,0.0,2.0,  1.0,0.0, 2.0,  2.0,0.0,2.0,  2.0,1.0,2.0,  2.0,2.0, 2.0,  1.0,2.0,2.0,</div><div>      0.0,2.0,2.0,  0.0,1.0, 2.0, 1.0,1.0,2.0};</div><div>Thanks</div><div>Amir</div><img class="mailspring-open" alt="Sent from Mailspring" width="0" height="0" style="border:0; width:0; height:0;" src="https://link.getmailspring.com/open/CE107FD7-8BE8-492B-98A7-90F4EEAA18F2@getmailspring.com?me=97a6ec29&recipient=cGV0c2MtdXNlcnNAbWNzLmFubC5nb3Y%3D"><div class="gmail_quote"><br><div>---------- Forwarded Message ---------</div><br><div>From: Matthew Knepley <knepley@gmail.com></div><div>Subject: Re: [petsc-users] View 3D DMPlex</div><div>Date: Sep 10 2019, at 9:08 pm</div><div>To: Amir <mhbaghaei@mail.sjtu.edu.cn></div><div>Cc: PETSc <petsc-users@mcs.anl.gov></div><br><div><div><div><div>On Tue, Sep 10, 2019 at 9:00 AM Amir via petsc-users <<a href="mailto:petsc-users@mcs.anl.gov" title="mailto:petsc-users@mcs.anl.gov">petsc-users@mcs.anl.gov</a>> wrote:</div></div><div class="gmail_quote"><blockquote><div>Hi</div><div>I am trying to view a cubic mesh constructed by DMPlex. I noticed that the interior point is not seen in the output VTU file.</div><div>It means that I do not see an edge inside the cube. I tried to check some detail of DM using --in_dm_view. The detail does not show any problem. Do you think there is a problem in my vtk output or dm setup.</div></blockquote><div><br></div><div>By default, Paraview does not show interior edges. You have to use a filter, like "Extract Edges".</div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><div> </div><blockquote><div>Thanks</div><div>Amir</div><div>DM_0x84000000_0 in 3 dimensions:</div><div>  0-cells: 27</div><div>  1-cells: 54</div><div>  2-cells: 36</div><div>  3-cells: 8</div><div>Labels:</div><div>  depth: 4 strata with value/size (0 (27), 1 (54), 2 (36), 3 (8))</div><br></blockquote></div><br><div><br></div><div>--</div><div class="gmail_signature"><div><div><div><div><div><div><div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.</div><div>-- Norbert Wiener</div></div><div><br></div><div><div><a href="http://www.cse.buffalo.edu/~knepley/" title="http://www.cse.buffalo.edu/~knepley/">https://www.cse.buffalo.edu/~knepley/</a></div></div></div></div></div></div></div></div></div></div></div><div class="gmail_quote_attribution">On Sep 10 2019, at 9:08 pm, Matthew Knepley <knepley@gmail.com> wrote:</div><blockquote><div><div><div><div>On Tue, Sep 10, 2019 at 9:00 AM Amir via petsc-users <<a href="https://link.getmailspring.com/link/CE107FD7-8BE8-492B-98A7-90F4EEAA18F2@getmailspring.com/0?redirect=mailto%3Apetsc-users%40mcs.anl.gov&recipient=cGV0c2MtdXNlcnNAbWNzLmFubC5nb3Y%3D" title="mailto:petsc-users@mcs.anl.gov">petsc-users@mcs.anl.gov</a>> wrote:</div></div><div class="gmail_quote"><blockquote><div>Hi</div><div>I am trying to view a cubic mesh constructed by DMPlex. I noticed that the interior point is not seen in the output VTU file.</div><div>It means that I do not see an edge inside the cube. I tried to check some detail of DM using --in_dm_view. The detail does not show any problem. Do you think there is a problem in my vtk output or dm setup.</div></blockquote><div><br></div><div>By default, Paraview does not show interior edges. You have to use a filter, like "Extract Edges".</div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><div> </div><blockquote><div>Thanks</div><div>Amir</div><div>DM_0x84000000_0 in 3 dimensions:</div><div>  0-cells: 27</div><div>  1-cells: 54</div><div>  2-cells: 36</div><div>  3-cells: 8</div><div>Labels:</div><div>  depth: 4 strata with value/size (0 (27), 1 (54), 2 (36), 3 (8))</div><br></blockquote></div><br><div><br></div><div>--</div><div class="gmail_signature"><div><div><div><div><div><div><div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.</div><div>-- Norbert Wiener</div></div><div><br></div><div><div><a href="https://link.getmailspring.com/link/CE107FD7-8BE8-492B-98A7-90F4EEAA18F2@getmailspring.com/1?redirect=http%3A%2F%2Fwww.cse.buffalo.edu%2F~knepley%2F&recipient=cGV0c2MtdXNlcnNAbWNzLmFubC5nb3Y%3D" title="http://www.cse.buffalo.edu/~knepley/">https://www.cse.buffalo.edu/~knepley/</a></div></div></div></div></div></div></div></div></div></div></blockquote>