PetscErrorCode DMCoarsen_ADDA(DM dm, MPI_Comm comm,DM *dmc)
{
PetscErrorCode ierr;
PetscInt *nodesc;
PetscInt dofc;
PetscInt i;
DM_ADDA *dd = (DM_ADDA*)dm->data;
PetscFunctionBegin;
PetscValidHeaderSpecific(dm, DM_CLASSID, 1);
PetscValidPointer(dmc, 3);
ierr = PetscMalloc(dd->dim*sizeof(PetscInt), &nodesc);CHKERRQ(ierr);
for (i=0; i<dd->dim; i++) {
nodesc[i] = (dd->nodes[i] % dd->refine[i]) ? dd->nodes[i] / dd->refine[i] + 1 : dd->nodes[i] / dd->refine[i];
}
dofc = (dd->dof % dd->dofrefine) ? dd->dof / dd->dofrefine + 1 : dd->dof / dd->dofrefine;
ierr = DMADDACreate(PetscObjectComm((PetscObject)dm), dd->dim, nodesc, dd->procs, dofc, dd->periodic, dmc);CHKERRQ(ierr);
ierr = PetscFree(nodesc);CHKERRQ(ierr);
/* copy refinement factors */
ierr = DMADDASetRefinement(*dmc, dd->refine, dd->dofrefine);CHKERRQ(ierr);
PetscFunctionReturn(0);
}
void PETSC_STDCALL dmaddasetrefinement_(DM dm,PetscInt *refine,PetscInt *dofrefine, int *__ierr ){
*__ierr = DMADDASetRefinement(
(DM)PetscToPointer((dm) ),refine,*dofrefine);
}
