<div dir="ltr"><div dir="ltr">On Tue, Sep 10, 2019 at 12:20 PM Amir <<a href="mailto:mhbaghaei@mail.sjtu.edu.cn">mhbaghaei@mail.sjtu.edu.cn</a>> wrote:<br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><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></blockquote><div><br></div><div>How about first using DMPlexCreateBoxMesh(), and seeing if you can visualize it. Then if that works, we can talk about</div><div>inputting your mesh from scratch.</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:1px solid rgb(204,204,204);padding-left:1ex"><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="gmail-m_3292043112138582813mailspring-open" alt="Sent from Mailspring" width="0" height="0" style="border: 0px; width: 0px; height: 0px;" src="https://link.getmailspring.com/open/CE107FD7-8BE8-492B-98A7-90F4EEAA18F2@getmailspring.com?me=97a6ec29&recipient=a25lcGxleUBnbWFpbC5jb20%3D"><div class="gmail_quote"><br><div>---------- Forwarded Message ---------</div><br><div>From: Matthew Knepley <<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>></div><div>Subject: Re: [petsc-users] View 3D DMPlex</div><div>Date: Sep 10 2019, at 9:08 pm</div><div>To: Amir <<a href="mailto:mhbaghaei@mail.sjtu.edu.cn" target="_blank">mhbaghaei@mail.sjtu.edu.cn</a>></div><div>Cc: PETSc <<a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a>></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" target="_blank">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-m_3292043112138582813gmail_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/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a></div></div></div></div></div></div></div></div></div></div></div><div class="gmail-m_3292043112138582813gmail_quote_attribution">On Sep 10 2019, at 9:08 pm, Matthew Knepley <<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>> 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=a25lcGxleUBnbWFpbC5jb20%3D" title="mailto:petsc-users@mcs.anl.gov" target="_blank">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-m_3292043112138582813gmail_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=a25lcGxleUBnbWFpbC5jb20%3D" title="http://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a></div></div></div></div></div></div></div></div></div></div></blockquote></blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>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><br></div><div><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>