PetscErrorCode DMPatchCreateGrid(MPI_Comm comm, PetscInt dim, MatStencil patchSize, MatStencil commSize, MatStencil gridSize, DM *dm)
{
DM_Patch *mesh;
DM da;
PetscInt dof = 1, width = 1;
PetscErrorCode ierr;
PetscFunctionBegin;
ierr = DMPatchCreate(comm, dm);CHKERRQ(ierr);
mesh = (DM_Patch*) (*dm)->data;
if (dim < 2) {
gridSize.j = 1;
patchSize.j = 1;
}
if (dim < 3) {
gridSize.k = 1;
patchSize.k = 1;
}
ierr = DMCreate(comm, &da);CHKERRQ(ierr);
ierr = DMSetType(da, DMDA);CHKERRQ(ierr);
ierr = DMSetDimension(da, dim);CHKERRQ(ierr);
ierr = DMDASetSizes(da, gridSize.i, gridSize.j, gridSize.k);CHKERRQ(ierr);
ierr = DMDASetBoundaryType(da, DM_BOUNDARY_NONE, DM_BOUNDARY_NONE, DM_BOUNDARY_NONE);CHKERRQ(ierr);
ierr = DMDASetDof(da, dof);CHKERRQ(ierr);
ierr = DMDASetStencilType(da, DMDA_STENCIL_BOX);CHKERRQ(ierr);
ierr = DMDASetStencilWidth(da, width);CHKERRQ(ierr);
mesh->dmCoarse = da;
ierr = DMPatchSetPatchSize(*dm, patchSize);CHKERRQ(ierr);
ierr = DMPatchSetCommSize(*dm, commSize);CHKERRQ(ierr);
PetscFunctionReturn(0);
}
PETSC_EXTERN void PETSC_STDCALL dmpatchcreate_(MPI_Fint * comm,DM *mesh, int *__ierr ){
*__ierr = DMPatchCreate(
MPI_Comm_f2c( *(comm) ),mesh);
}
