/* "Unintrusive" multi-physics DM */
#if !defined(PETSCDMCOMPOSITE_H)
#define PETSCDMCOMPOSITE_H

#include <petscdm.h>

PETSC_EXTERN PetscErrorCode DMCompositeCreate(MPI_Comm, DM *);
PETSC_EXTERN PetscErrorCode DMCompositeAddDM(DM, DM);
PETSC_EXTERN PetscErrorCode DMCompositeSetCoupling(DM, PetscErrorCode (*)(DM, Mat, PetscInt *, PetscInt *, PetscInt, PetscInt, PetscInt, PetscInt));
PETSC_EXTERN PetscErrorCode DMCompositeAddVecScatter(DM, VecScatter);
PETSC_EXTERN PetscErrorCode DMCompositeScatter(DM, Vec, ...);
PETSC_EXTERN PetscErrorCode DMCompositeScatterArray(DM, Vec, Vec *);
PETSC_EXTERN PetscErrorCode DMCompositeGather(DM, InsertMode, Vec, ...);
PETSC_EXTERN PetscErrorCode DMCompositeGatherArray(DM, InsertMode, Vec, Vec *);
PETSC_EXTERN PetscErrorCode DMCompositeGetNumberDM(DM, PetscInt *);
PETSC_EXTERN PetscErrorCode DMCompositeGetAccess(DM, Vec, ...);
PETSC_EXTERN PetscErrorCode DMCompositeRestoreAccess(DM, Vec, ...);
PETSC_EXTERN PetscErrorCode DMCompositeGetAccessArray(DM, Vec, PetscInt, const PetscInt *, Vec *);
PETSC_EXTERN PetscErrorCode DMCompositeRestoreAccessArray(DM, Vec, PetscInt, const PetscInt *, Vec *);
PETSC_EXTERN PetscErrorCode DMCompositeGetLocalAccessArray(DM, Vec, PetscInt, const PetscInt *, Vec *);
PETSC_EXTERN PetscErrorCode DMCompositeRestoreLocalAccessArray(DM, Vec, PetscInt, const PetscInt *, Vec *);
PETSC_EXTERN PetscErrorCode DMCompositeGetLocalVectors(DM, ...);
PETSC_EXTERN PetscErrorCode DMCompositeGetEntries(DM, ...);
PETSC_EXTERN PetscErrorCode DMCompositeGetEntriesArray(DM, DM[]);
PETSC_EXTERN PetscErrorCode DMCompositeRestoreLocalVectors(DM, ...);
PETSC_EXTERN PetscErrorCode DMCompositeGetGlobalISs(DM, IS *[]);
PETSC_EXTERN PetscErrorCode DMCompositeGetLocalISs(DM, IS **);
PETSC_EXTERN PetscErrorCode DMCompositeGetISLocalToGlobalMappings(DM, ISLocalToGlobalMapping **);

#endif