[petsc-users] Gmsh 8-noded quadrilateral

Matthew Knepley knepley at gmail.com
Thu Feb 10 09:23:44 CST 2022

On Thu, Feb 10, 2022 at 10:12 AM Susanne Claus <susanne.claus at onera.fr>

> Hello,
> I am using DMPlex for the mesh structure of a solid mechanics finite
> element code. I mainly use gmsh as input file format. When I try to read in
> 8-noded Quadrilaterals (Element type 16 in gmsh) DMPlex tells me that this
> element type is unknown. However a 9-noded Quadrilateral can be read
> without problem. On inspecting the plexgmsh.c source code I can see that
> 8-noded quadrilaterals are deactivated:
> #if 0146:   {20, GMSH_TRI, 2, 3, 3,  9, NULL},147:   {16, GMSH_QUA, 2, 2, 4,  8, NULL},
> For our application these 8-noded quadrilateral are very important.
> Is there any reason why they have not been implemented/deactivated in the dmplex gmsh reader?
> No, we can handle them in the same way I think. Let me look at it.
Hopefully it is easy.



> Thank you for all the great work you are doing. PETSc is amazing.
> Best wishes,
> Susanne Claus

What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener

https://www.cse.buffalo.edu/~knepley/ <http://www.cse.buffalo.edu/~knepley/>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20220210/f6d14980/attachment-0001.html>

More information about the petsc-users mailing list