<div dir="ltr"><div dir="ltr">On Wed, Sep 11, 2019 at 3:05 AM 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"><div>I tried DMPlexCreateBoxMesh(). It shows points and nodes on the boundary surface. Plus, I do  see easily and correctly the interior edge/node inside by clipping the volume mesh. Do you think my dm setup is wrong that I cant see the same using my DAG settings in Paraview.</div><div></div></blockquote><div><br></div><div>Yes. You can see the order of the BoxMesh cells by using</div><div><br></div><div>  -dm_view ::ascii_info_detail</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>Thanks</div><div>Amir</div><div> </div><br><div class="gmail-m_8695782344422025668gmail_quote_attribution">On Sep 11 2019, at 2:27 am, 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 12:20 PM Amir <<a href="https://link.getmailspring.com/link/D5CB7218-13F0-46E6-9D32-1BEE866C43CA@getmailspring.com/0?redirect=mailto%3Amhbaghaei%40mail.sjtu.edu.cn&recipient=a25lcGxleUBnbWFpbC5jb20%3D" title="mailto:mhbaghaei@mail.sjtu.edu.cn" target="_blank">mhbaghaei@mail.sjtu.edu.cn</a>> wrote:</div></div><img class="gmail-m_8695782344422025668mailspring-open" alt="Sent from Mailspring" width="0" height="0" style="border: 0px; width: 0px; height: 0px;" src="https://link.getmailspring.com/open/D5CB7218-13F0-46E6-9D32-1BEE866C43CA@getmailspring.com?me=97a6ec29&recipient=a25lcGxleUBnbWFpbC5jb20%3D"><div class="gmail_quote"><blockquote><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><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><div class="gmail_quote"><br><div>---------- Forwarded Message ---------</div><br><div>From: Matthew Knepley <<a href="mailto:knepley@gmail.com" title="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" title="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" title="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_8695782344422025668gmail-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_8695782344422025668gmail-m_3292043112138582813gmail_quote_attribution">On Sep 10 2019, at 9:08 pm, Matthew Knepley <<a href="mailto:knepley@gmail.com" title="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="https://link.getmailspring.com/link/CE107FD7-8BE8-492B-98A7-90F4EEAA18F2@getmailspring.com/0?redirect=mailto%3Apetsc-users%40mcs.anl.gov&recipient=a25lcGxleUBnbWFpbC5jb20%3D" 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_8695782344422025668gmail-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="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" target="_blank">https://www.cse.buffalo.edu/~knepley/</a></div></div></div></div></div></div></div></div></div></div></blockquote></blockquote></div><br><div><br></div><div>--</div><div class="gmail-m_8695782344422025668gmail_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/D5CB7218-13F0-46E6-9D32-1BEE866C43CA@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>