#include "petscdm.h" #include "petscdmlabel.h" #include "petscds.h" PetscErrorCode DMGetVecType(DM da,VecType *ctype)