1af0996ceSBarry Smith #include <petsc/private/dmdaimpl.h> /*I "petscdmda.h" I*/ 247c6ae99SBarry Smith 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