On Fri, Jun 30, 2023 at 4:40 PM neil liu <liufi...@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/>
