1af0996ceSBarry Smith #include <petsc/private/dmdaimpl.h> /*I "petscdmda.h" I*/ 247c6ae99SBarry Smith 37087cfbeSBarry Smith extern PetscErrorCode DMSetUp_DA_1D(DM); 47087cfbeSBarry Smith extern PetscErrorCode DMSetUp_DA_2D(DM); 57087cfbeSBarry Smith extern PetscErrorCode DMSetUp_DA_3D(DM); 647c6ae99SBarry Smith 77087cfbeSBarry Smith PetscErrorCode DMSetUp_DA(DM da) 847c6ae99SBarry Smith { 947c6ae99SBarry Smith PetscErrorCode ierr; 10d7bf68aeSBarry Smith DM_DA *dd = (DM_DA*)da->data; 1147c6ae99SBarry Smith 1247c6ae99SBarry Smith PetscFunctionBegin; 1347c6ae99SBarry Smith PetscValidHeaderSpecific(da, DM_CLASSID,1); 14*2c71b3e2SJacob Faibussowitsch PetscCheckFalse(dd->w < 1,PetscObjectComm((PetscObject)da),PETSC_ERR_ARG_OUTOFRANGE,"Must have 1 or more degrees of freedom per node: %D",dd->w); 15*2c71b3e2SJacob Faibussowitsch PetscCheckFalse(dd->s < 0,PetscObjectComm((PetscObject)da),PETSC_ERR_ARG_OUTOFRANGE,"Stencil width cannot be negative: %D",dd->s); 16109c9344SBarry Smith 17c629b14aSBarry Smith ierr = PetscCalloc1(dd->w+1,&dd->fieldname);CHKERRQ(ierr); 18c73cfb54SMatthew G. Knepley ierr = PetscCalloc1(da->dim,&dd->coordinatename);CHKERRQ(ierr); 19c73cfb54SMatthew G. Knepley if (da->dim == 1) { 209a42bb27SBarry Smith ierr = DMSetUp_DA_1D(da);CHKERRQ(ierr); 21c73cfb54SMatthew G. Knepley } else if (da->dim == 2) { 229a42bb27SBarry Smith ierr = DMSetUp_DA_2D(da);CHKERRQ(ierr); 23c73cfb54SMatthew G. Knepley } else if (da->dim == 3) { 249a42bb27SBarry Smith ierr = DMSetUp_DA_3D(da);CHKERRQ(ierr); 25ce94432eSBarry Smith } else SETERRQ(PetscObjectComm((PetscObject)da),PETSC_ERR_SUP,"DMs only supported for 1, 2, and 3d"); 269200b928SMatthew G. Knepley ierr = DMViewFromOptions(da,NULL,"-dm_view");CHKERRQ(ierr); 2747c6ae99SBarry Smith PetscFunctionReturn(0); 2847c6ae99SBarry Smith } 29