<div dir="ltr"><div dir="ltr">On Mon, Nov 22, 2021 at 12:10 PM Bhargav Subramanya <<a href="mailto:bhargav.subramanya@kaust.edu.sa" target="_blank">bhargav.subramanya@kaust.edu.sa</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 dir="ltr">Dear All,<div><br></div><div>I have a prismatic mesh generated by extruding a base icosahedron spherical surface mesh. The mesh is distributed on 2 processes after extrusion. The dof for my problem are defined only on the simplices of the spherical layers.</div><div><br></div><div>1. For a given spherical layer, I want to get points, which are basically simplices, lying on that layer only. DMPlexGetHeightStratum returns the points on all the spherical Layers. I can probably use DMPlexFindVertices (based on the radius of the spherical layer) and DMPlexGetSimplexOrBoxCells. Could you suggest if there is a better way to retrieve the points (which are simplices) on a given spherical layer in the extruded mesh?</div></div></blockquote><div><br></div><div>DMPlexGetHeightStratum() refers to height in the Hasse Diagram, which is a DAG, not in the mesh. For example, height 0 points are the cells, height 1 are the faces, etc.<br></div><div><br></div><div>I believe the default numbering for extrusion (in the main branch), is that all vertices produced from a given vertex be in order. This would mean that if v were the vertex point number, then</div><div><br></div><div>  (v - vStart) % Nl == l</div><div><br></div><div>where Nl is the number of layers and l is the layer of that vertex. It is also the same for triangles. So if you want to segregate each shell, after extrusion, make a label that gives triangles this marker, and then use DMPlexLabelComplete(). Then after refinement you will still have the shells labeled correctly.</div><div><br></div><div>I would be happy to help you make an example that does this. It seems cool.</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 dir="ltr"><div>2. I am able to visualize the entire mesh using dm_view. Is there a way to get the mesh file for the local dm from a specific process? </div></div></blockquote><div><br></div><div>You can use -dm_partition_view which outputs a field with the process number. Then use Clip inside Paraview and clip to the process number you want,</div><div>or just view that field so each process has a different color.</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 dir="ltr"><div>3. One of my previous emails might have got lost as I mistakenly attached a large mesh file and sent it. So I am repeating the question here. DMPlexExtrude gives the following error after distributing the base 2D spherical mesh. Both refinement or/and extrusion after dm distribution does not work for me. In fact, I tried with src/dm/impls/plex/tutorials/ex10.c also. Although the mesh is distributed after the extrusion in ex10.c (which is working fine for me), I tried to distribute before extrusion, and I get the following error. Could you please suggest where I might be making a mistake?</div></div></blockquote><div><br></div><div>So you want to distribute the mesh before extruding. For that small example (using the main branch), I run</div><div><br></div><div>  mpiexec -n 3 ./meshtest -dm_plex_shape sphere -dm_refine_pre 2 -dm_distribute -dm_partition_view -dm_view hdf5:mesh.h5 -dm_extrude 3</div><div><br></div><div>and I get the attached picture.</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 dir="ltr"><div>[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------<br>[0]PETSC ERROR: Object is in wrong state<br>[0]PETSC ERROR: This DMPlex is distributed but its PointSF has no graph set<br>[0]PETSC ERROR: See <a href="https://petsc.org/release/faq/" target="_blank">https://petsc.org/release/faq/</a> for trouble shooting.<br>[0]PETSC ERROR: Petsc Release Version 3.16.1, unknown<br>[0]PETSC ERROR: ./cavity_flow.out on a arch-darwin-c-debug named kl-21859 by subrambm Mon Nov 22 19:47:14 2021<br>[0]PETSC ERROR: Configure options --with-cc=gcc --with-cxx=g++ --with-fc=gfortran --download-mpich --download-fblaslapack --download-superlu_dist --download-hypre --download-fiat --download-generator --download-triangle --download-tetgen --download-chaco --download-make -download-boost --download-cmake --download-ml --download-mumps=<a href="https://bitbucket.org/petsc/pkg-mumps.git" target="_blank">https://bitbucket.org/petsc/pkg-mumps.git</a> --download-mumps-commit=v5.4.1-p1 --download-scalapack --download-ptscotch --download-hdf5 --force<br>[0]PETSC ERROR: #1 DMPlexCheckPointSF() at /Users/subrambm/petsc/src/dm/impls/plex/plex.c:8554<br>[0]PETSC ERROR: #2 DMPlexOrientInterface_Internal() at /Users/subrambm/petsc/src/dm/impls/plex/plexinterpolate.c:595<br>[0]PETSC ERROR: #3 DMPlexInterpolate() at /Users/subrambm/petsc/src/dm/impls/plex/plexinterpolate.c:1357<br>[0]PETSC ERROR: #4 DMPlexExtrude() at /Users/subrambm/petsc/src/dm/impls/plex/plexcreate.c:1543<br>[0]PETSC ERROR: #5 CreateMesh() at ex10.c:161<br>[0]PETSC ERROR: #6 main() at ex10.c:180<br>[0]PETSC ERROR: PETSc Option Table entries:<br>[0]PETSC ERROR: -dm_plex_extrude_layers 3<br>[0]PETSC ERROR: -dm_view vtk:mesh.vtk<br>[0]PETSC ERROR: -init_dm_plex_dim 2<br>[0]PETSC ERROR: -petscpartitioner_type simple<br>[0]PETSC ERROR: -srf_dm_distribute<br>[0]PETSC ERROR: -srf_dm_refine 0<br>[0]PETSC ERROR: ----------------End of Error Message -------send entire error message to petsc-maint@mcs.anl.gov----------</div><div><br></div><div>Thanks,</div><div>Bhargav</div><div><br></div></div>

<br>
<div><hr></div><font face="Arial" size="1">This message and its contents, including attachments are intended solely for the original recipient. If you are not the intended recipient or have received this message in error, please notify me immediately and delete this message from your computer system. Any unauthorized use or distribution is prohibited. Please consider the environment before printing this email.</font></blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr"><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>