1: #include <petsc/private/dmdaimpl.h>
3: extern PetscErrorCode DMSetUp_DA_1D(DM);
4: extern PetscErrorCode DMSetUp_DA_2D(DM);
5: extern PetscErrorCode DMSetUp_DA_3D(DM);
7: PetscErrorCode DMSetUp_DA(DM da)
8: {
9: DM_DA *dd = (DM_DA *)da->data;
15: PetscCalloc1(dd->w + 1, &dd->fieldname);
16: PetscCalloc1(da->dim, &dd->coordinatename);
17: if (da->dim == 1) {
18: DMSetUp_DA_1D(da);
19: } else if (da->dim == 2) {
20: DMSetUp_DA_2D(da);
21: } else if (da->dim == 3) {
22: DMSetUp_DA_3D(da);
23: } else SETERRQ(PetscObjectComm((PetscObject)da), PETSC_ERR_SUP, "DMs only supported for 1, 2, and 3d");
24: DMViewFromOptions(da, NULL, "-dm_view");
25: return 0;
26: }