#include "petscdmplex.h" #include "petscfe.h" PetscErrorCode DMPlexComputeGeometryFEM(DM dm, Vec *cellgeom)