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/>

Reply via email to