Actual source code: dareg.c

  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: }