1 #if !defined(__PETSCMOAB_H) 2 #define __PETSCMOAB_H 3 4 #include <petscvec.h> 5 #include <petscmat.h> 6 #include <petscdm.h> 7 8 #include <moab/Core.hpp> 9 #include <moab/ParallelComm.hpp> 10 11 /* The MBERR macro is used to save typing. It checks a MOAB error code 12 * (rval) and calls SETERRQ if not MB_SUCCESS. A message (msg) can 13 * also be passed in. */ 14 #define MBERR(msg,rval) do{if(rval != moab::MB_SUCCESS) SETERRQ2(PETSC_COMM_SELF,PETSC_ERR_LIB,"MOAB ERROR (%i): %s",rval,msg);} while(0) 15 #define MBERRNM(rval) do{if(rval != moab::MB_SUCCESS) SETERRQ1(PETSC_COMM_SELF,PETSC_ERR_LIB,"MOAB ERROR (%i)",rval);} while(0) 16 17 PETSC_EXTERN PetscErrorCode DMMoabCreate(MPI_Comm comm, DM *moab); 18 PETSC_EXTERN PetscErrorCode DMMoabSetParallelComm(DM dm,moab::ParallelComm *pcomm); 19 PETSC_EXTERN PetscErrorCode DMMoabGetParallelComm(DM dm,moab::ParallelComm **pcomm); 20 PETSC_EXTERN PetscErrorCode DMMoabSetInterface(DM dm,moab::Interface *iface); 21 PETSC_EXTERN PetscErrorCode DMMoabGetInterface(DM dm,moab::Interface **mbint); 22 PETSC_EXTERN PetscErrorCode DMMoabSetRange(DM dm,moab::Range range); 23 PETSC_EXTERN PetscErrorCode DMMoabGetRange(DM dm,moab::Range *range); 24 PETSC_EXTERN PetscErrorCode DMMoabSetLocalToGlobalTag(DM dm,moab::Tag ltogtag); 25 PETSC_EXTERN PetscErrorCode DMMoabGetLocalToGlobalTag(DM dm,moab::Tag *ltog_tag); 26 27 PETSC_EXTERN PetscErrorCode DMMoabCreateFromInstance(MPI_Comm comm,moab::Interface *iface,moab::ParallelComm *pcomm,moab::Tag ltog_tag,moab::Range *range,DM *moab); 28 PETSC_EXTERN PetscErrorCode DMMoabCreateDMAndInstance(MPI_Comm comm,DM *moab); 29 PETSC_EXTERN PetscErrorCode DMMoabCreateVector(DM dm,PetscInt tag_size,moab::Range range,PetscBool serial,PetscBool destroy_tag,Vec *vec); 30 PETSC_EXTERN PetscErrorCode DMMoabCreateVectorFromTag(DM dm, moab::Tag tag,moab::Range range,PetscBool serial, PetscBool destroy_tag,Vec *X); 31 PETSC_EXTERN PetscErrorCode DMMoabGetVecTag(Vec vec,moab::Tag *tag); 32 PETSC_EXTERN PetscErrorCode DMMoabGetVecRange(Vec vec,moab::Range *range); 33 PETSC_EXTERN PetscErrorCode DMDestroy_Moab(DM dm); 34 PETSC_EXTERN PetscErrorCode DMInitialize_Moab(DM dm); 35 36 PETSC_EXTERN PetscErrorCode DMCreateGlobalVector_Moab(DM dm,Vec *gvec); 37 PETSC_EXTERN PetscErrorCode DMCreateLocalVector_Moab(DM dm,Vec *gvec); 38 PETSC_EXTERN PetscErrorCode DMCreateGlobalVectorFromTag_Moab(DM dm,moab::Tag tag,Vec *gvec); 39 #endif 40