[petsc-users] DMPlex H27 elements

[Yann Jobic]

Dear Petsc Users,

I'm using DMPlexCreateFromCellList to create my DM. I would like to have 
an order 2 geometry.

It's working fine in 2D for elements of type Q9.

I checked that it's working correctly by using 
DMPlexComputeCellGeometryFEM, and compute the value of the determinant 
(CheckMeshGeometry of dm/impls/plex/examples/tutorials/ex2.c)

I can import H8 elements, it's working fine.

But for H27 element, it's not working. I really don't know how to order 
my vertex in the cells array. So far, the determinant is zero...

I don't know where to look in order to find this information. I tried 
the function DMPlexGetRawFaces_Internal of the file

https://www.mcs.anl.gov/petsc/petsc-current/src/dm/impls/plex/plexinterpolate.c.html#DMPlexGetRawFaces_Internal

but it didn't help me.

Could you please point me where to look ?

Thanks, and happy new year !

Yann