[petsc-users] Inquiry about reading the P2 tetrahedron mesh from GMSH

[Matthew Knepley]

On Fri, Jun 30, 2023 at 4:40 PM neil liu <liufield at gmail.com> wrote:

> Dear Petsc developers,
>
> I am reading P2 mesh from GMSH. And used DMFieldGetClosure_Internal to
> check the coordinates for each tetrahedron, It seems reasonable.
> But when I tried DMGetCoordinates (dm, &global),  it seems the vector
> global is not consistent with the node number, Then what is global here?
>

Reading high order coordinates works for quads/hexes, but not simplices
right now because I need to figure out what order GMsh stores the
coordinates in. Thus we are defaulting to affine coordinates here.

  Thanks,

     Matt


> Thanks,
>
> Xiaodong
>


-- 
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/20230717/c7ddab37/attachment.html>