#include "petscdm.h" #include "petscdmlabel.h" #include "petscds.h" PetscErrorCode DMCreateMassMatrix(DM dac, DM daf, Mat *mat)