160 …scErrorCode DMCreateSubDM_DA(DM dm, PetscInt numFields, const PetscInt fields[], IS *is, DM *subdm)  in DMCreateSubDM_DA()  argument
165   if (subdm) {  in DMCreateSubDM_DA()
171     PetscCall(DMCreate(PetscObjectComm((PetscObject)dm), subdm));  in DMCreateSubDM_DA()
173     PetscCall(DMSetPointSF(*subdm, sf));  in DMCreateSubDM_DA()
175     PetscCall(DMSetApplicationContext(*subdm, ctx));  in DMCreateSubDM_DA()
178       PetscCall(DMSetCoordinatesLocal(*subdm, coords));  in DMCreateSubDM_DA()
181       if (coords) PetscCall(DMSetCoordinates(*subdm, coords));  in DMCreateSubDM_DA()
184     PetscCall(DMSetType(*subdm, DMDA));  in DMCreateSubDM_DA()
185     PetscCall(DMSetDimension(*subdm, dm->dim));  in DMCreateSubDM_DA()
186     PetscCall(DMDASetSizes(*subdm, da->M, da->N, da->P));  in DMCreateSubDM_DA()
187     PetscCall(DMDASetNumProcs(*subdm, da->m, da->n, da->p));  in DMCreateSubDM_DA()
188     PetscCall(DMDASetBoundaryType(*subdm, da->bx, da->by, da->bz));  in DMCreateSubDM_DA()
189     PetscCall(DMDASetDof(*subdm, numFields));  in DMCreateSubDM_DA()
190     PetscCall(DMDASetStencilType(*subdm, da->stencil_type));  in DMCreateSubDM_DA()
191     PetscCall(DMDASetStencilWidth(*subdm, da->s));  in DMCreateSubDM_DA()
192     PetscCall(DMDASetOwnershipRanges(*subdm, da->lx, da->ly, da->lz));  in DMCreateSubDM_DA()