<div dir="ltr"><div><div><div><div><div>Hello all,<br><br></div>Recently I've been trying to run a fem problem using DMPlex, PetscFE, etc. using a mesh created with gmsh. I used DMPlexReverseCell to get rid of the negative jacobian determinants. However, I noticed another problem: when I projected the coordinates on a quadratic fem-field and interpolated them they didn't match the interpolated values of the original coordinate field anymore.<br><br></div>I attached a simple test to reproduce this. With linear interpolation I see no problem, it's only when I start using quadratic interpolation.<br><br></div>I have no idea why this is giving problems. Am I missing something simple here or is this a bug?<br><br></div>Thanks,<br></div>Sander<br></div>