<div dir="ltr"><div dir="ltr">On Thu, Jul 14, 2022 at 6:33 PM Mike Michell <<a href="mailto:mi.mike1021@gmail.com">mi.mike1021@gmail.com</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"><div dir="ltr">VTK format itself supports pyramids for sure. But I do not think PETSc output part through VTK format does not. Is there any possibility to debug this issue?<br></div></div></blockquote><div><br></div><div>VTK output now works with pyramids in</div><div><br></div><div> <a href="https://gitlab.com/petsc/petsc/-/merge_requests/5422">https://gitlab.com/petsc/petsc/-/merge_requests/5422</a></div><div><br></div><div>Could you send me a much smaller GMsh that has pyramids and tets or hexes in it? I would like to put in a test</div><div>before we merge this.</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 dir="ltr">Thanks,<br></div><br><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"><div dir="ltr"><br>On Wed, Jul 13, 2022 at 4:18 PM Mike Michell <<a href="mailto:mi.mike1021@gmail.com" target="_blank">mi.mike1021@gmail.com</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"><div dir="ltr">Before your pull request merged to main, I copied your change to my local /src/dm/impls/plex/plexgmsh.c to test pyramids.<br><br>Loading DMPlex with the mesh including pyramids works okay, but printing out the solution field to .vtk or .vtu format causes a problem with "Unknown Cell Type Error" from post-processing tools (such as tecplot or paraview).<br><br>I had a look at /src/dm/impls/plex/plexvtk.c the function, DMPlexVTKGetCellType_Internal() seems not explicitly including a cellType marker for prism and pyramids. Probably, "VTK_WEDGE" is prism with triangular base, but still, there seems no pyramids marker.<br></div></div></blockquote><div><br></div><div>Yes, I am not sure it supports pyramids. Will have to look.</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 dir="ltr">Thanks,<br></div><br><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"><div dir="ltr"><br>On Wed, Jul 13, 2022 at 11:18 AM Mike Michell <<a href="mailto:mi.mike1021@gmail.com" target="_blank">mi.mike1021@gmail.com</a>> wrote:<br></div><div dir="ltr"><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"><div dir="ltr">Attached is a mixed mesh file that I am testing. I cannot see any special marker for pyramid cells. Version of gmsh is 4.9.0. </div></div></blockquote><div><br></div><div>You are correct. Pyramids were disabled. I have activated them here:</div><div><br></div><div> <a href="https://gitlab.com/petsc/petsc/-/merge_requests/5422" target="_blank">https://gitlab.com/petsc/petsc/-/merge_requests/5422</a></div><div><br></div><div>You mesh runs fine after this for me. I am not sure things like geometry will work for pyramids. However, if you find something</div><div>broken, just let me know.</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 dir="ltr"><div>Thanks,</div></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"><div dir="ltr"><br>On Wed, Jul 13, 2022 at 10:41 AM Mike Michell <<a href="mailto:mi.mike1021@gmail.com" target="_blank">mi.mike1021@gmail.com</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"><div dir="ltr">Thank you for the quick response. Below is the full error message I get. <div><br></div><div><br></div><div>[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------<br>[0]PETSC ERROR: Argument out of range<br>[0]PETSC ERROR: No face description for cell type unknown<br></div></div></div></blockquote><div><br></div><div>Here is the problem. The pyramid (or some other cell) was classified as "unknown". First, make sure the Gmsh file is version 4.1.</div><div>If that fails, send the Gmsh file and I will try to figure out why the cell is not coming up as a pyramid.</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 dir="ltr"><div>[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.17.0, unknown <br>[0]PETSC ERROR: /home/Mike/Workspace/test on a named Mike Wed Jul 13 10:38:31 2022<br>[0]PETSC ERROR: Configure options --prefix=/home/Mike/Library/petsc_partition/install-intel PETSC_ARCH=linux-gnu-intel --with-cc=mpiicc --with-cxx=mpiicpc --with-fc=mpiifort --download-fblaslapack --download-metis --download-parmetis --download-eigen --download-pragmatic --download-hdf5 --download-triangle --with-debugging=1 COPTFLAGS="-O3 -mtune=native" CXXOPTFLAGS=-O3 FOPTFLAGS=-O3<br>[0]PETSC ERROR: #1 DMPlexGetRawFaces_Internal() at /home/Mike/Library/petsc_partition/src/dm/impls/plex/plexinterpolate.c:312<br>[0]PETSC ERROR: #2 DMPlexInterpolateFaces_Internal() at /home/Mike/Library/petsc_partition/src/dm/impls/plex/plexinterpolate.c:350<br>[0]PETSC ERROR: #3 DMPlexInterpolate() at /home/Mike/Library/petsc_partition/src/dm/impls/plex/plexinterpolate.c:1327<br>[0]PETSC ERROR: #4 DMPlexCreateGmsh() at /home/Mike/Library/petsc_partition/src/dm/impls/plex/plexgmsh.c:1634<br>[0]PETSC ERROR: #5 DMPlexCreateGmshFromFile() at /home/Mike/Library/petsc_partition/src/dm/impls/plex/plexgmsh.c:1418<br>[0]PETSC ERROR: #6 DMPlexCreateFromFile() at /home/Mike/Library/petsc_partition/src/dm/impls/plex/plexcreate.c:4721<br>[0]PETSC ERROR: #7 DMPlexCreateFromOptions_Internal() at /home/Mike/Library/petsc_partition/src/dm/impls/plex/plexcreate.c:3212<br>[0]PETSC ERROR: #8 DMSetFromOptions_Plex() at /home/Mike/Library/petsc_partition/src/dm/impls/plex/plexcreate.c:3433<br>[0]PETSC ERROR: #9 DMSetFromOptions() at /home/Mike/Library/petsc_partition/src/dm/interface/dm.c:887<br>[0]PETSC ERROR: #10 User provided function() at User file:0<br>Abort(63) on node 0 (rank 0 in comm 16): application called MPI_Abort(MPI_COMM_SELF, 63) - process 0<br></div><div><br></div><div>Thanks,</div><div>Mike</div></div><br><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"><div dir="ltr"><br>On Wed, Jul 13, 2022 at 10:30 AM Mike Michell <<a href="mailto:mi.mike1021@gmail.com" target="_blank">mi.mike1021@gmail.com</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">Hi, DMCreate() is used to load/distribute grid built from gmsh, and the function crashes when a mixed mesh of tetra and pyramids in 3D. It looks PETSc can handle the cell types of Tetra/Hexa/Prism/Pyramids/Polygon/Polyhedra in 3D. Thus I am unsure why it crashes when pyramids included in gmsh file. <br><br>Below functions are used to mesh distribution.<br>call DMCreate(PETSC_COMM_WORLD, dm_g, ierr);CHKERRA(ierr)<br>call DMSetType(dm_g, DMPLEX, ierr);CHKERRA(ierr)<br>call DMSetFromOptions(dm_g, ierr);CHKERRA(ierr)<br><br>Below error messages got from run.<br></div></blockquote><div><br></div><div>I need the complete error message. This is only the stack. I cannot see the error message, or the version of PETSc you are using.</div><div><br></div><div>This should work, so I do not immediately know what is wrong.</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">[0]PETSC ERROR: #1 DMPlexGetRawFaces_Internal() at /home/Mike/Library/petsc_partition/src/dm/impls/plex/plexinterpolate.c:312<br>[0]PETSC ERROR: #2 DMPlexInterpolateFaces_Internal() at /home/Mike/Library/petsc_partition/src/dm/impls/plex/plexinterpolate.c:350<br>[0]PETSC ERROR: #3 DMPlexInterpolate() at /home/Mike/Library/petsc_partition/src/dm/impls/plex/plexinterpolate.c:1327<br>[0]PETSC ERROR: #4 DMPlexCreateGmsh() at /home/Mike/Library/petsc_partition/src/dm/impls/plex/plexgmsh.c:1634<br>[0]PETSC ERROR: #5 DMPlexCreateGmshFromFile() at /home/Mike/Library/petsc_partition/src/dm/impls/plex/plexgmsh.c:1418<br>[0]PETSC ERROR: #6 DMPlexCreateFromFile() at /home/Mike/Library/petsc_partition/src/dm/impls/plex/plexcreate.c:4721<br>[0]PETSC ERROR: #7 DMPlexCreateFromOptions_Internal() at /home/Mike/Library/petsc_partition/src/dm/impls/plex/plexcreate.c:3212<br>[0]PETSC ERROR: #8 DMSetFromOptions_Plex() at /home/Mike/Library/petsc_partition/src/dm/impls/plex/plexcreate.c:3433<br>[0]PETSC ERROR: #9 DMSetFromOptions() at /home/Mike/Library/petsc_partition/src/dm/interface/dm.c:887<br>[0]PETSC ERROR: #10 User provided function() at User file:0<br>Abort(63) on node 0 (rank 0 in comm 16): application called MPI_Abort(MPI_COMM_SELF, 63) - process 0<br><br>Can I get any comments on that? <br><br>Thanks,<br><br><br></div>
</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>
</blockquote></div></div>
</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>
</blockquote></div></div>
</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>
</div>
</blockquote></div></div>
</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>
</blockquote></div></div>
</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>