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