1af0996ceSBarry Smith #include <petsc/private/dmdaimpl.h> /*I "petscdmda.h" I*/
247c6ae99SBarry Smith
DMSetUp_DA(DM da)3d71ae5a4SJacob Faibussowitsch PetscErrorCode DMSetUp_DA(DM da)
4d71ae5a4SJacob Faibussowitsch {
5d7bf68aeSBarry Smith DM_DA *dd = (DM_DA *)da->data;
647c6ae99SBarry Smith
747c6ae99SBarry Smith PetscFunctionBegin;
847c6ae99SBarry Smith PetscValidHeaderSpecific(da, DM_CLASSID, 1);
963a3b9bcSJacob Faibussowitsch PetscCheck(dd->w >= 1, PetscObjectComm((PetscObject)da), PETSC_ERR_ARG_OUTOFRANGE, "Must have 1 or more degrees of freedom per node: %" PetscInt_FMT, dd->w);
1063a3b9bcSJacob Faibussowitsch PetscCheck(dd->s >= 0, PetscObjectComm((PetscObject)da), PETSC_ERR_ARG_OUTOFRANGE, "Stencil width cannot be negative: %" PetscInt_FMT, dd->s);
11109c9344SBarry Smith
129566063dSJacob Faibussowitsch PetscCall(PetscCalloc1(dd->w + 1, &dd->fieldname));
139566063dSJacob Faibussowitsch PetscCall(PetscCalloc1(da->dim, &dd->coordinatename));
14c73cfb54SMatthew G. Knepley if (da->dim == 1) {
159566063dSJacob Faibussowitsch PetscCall(DMSetUp_DA_1D(da));
16c73cfb54SMatthew G. Knepley } else if (da->dim == 2) {
179566063dSJacob Faibussowitsch PetscCall(DMSetUp_DA_2D(da));
18c73cfb54SMatthew G. Knepley } else if (da->dim == 3) {
199566063dSJacob Faibussowitsch PetscCall(DMSetUp_DA_3D(da));
20ce94432eSBarry Smith } else SETERRQ(PetscObjectComm((PetscObject)da), PETSC_ERR_SUP, "DMs only supported for 1, 2, and 3d");
21*4c4edb6fSBarry Smith da->Nf = dd->w;
229566063dSJacob Faibussowitsch PetscCall(DMViewFromOptions(da, NULL, "-dm_view"));
233ba16761SJacob Faibussowitsch PetscFunctionReturn(PETSC_SUCCESS);
2447c6ae99SBarry Smith }
25