1c4762a1bSJed Brown 2c4762a1bSJed Brown static char help[] = "Tests DMDA with variable multiple degrees of freedom per node.\n\n"; 3c4762a1bSJed Brown 4c4762a1bSJed Brown /* 5c4762a1bSJed Brown This code only compiles with gcc, since it is not ANSI C 6c4762a1bSJed Brown */ 7c4762a1bSJed Brown 8c4762a1bSJed Brown #include <petscdm.h> 9c4762a1bSJed Brown #include <petscdmda.h> 10c4762a1bSJed Brown 11c4762a1bSJed Brown PetscErrorCode doit(DM da,Vec global) 12c4762a1bSJed Brown { 13c4762a1bSJed Brown PetscErrorCode ierr; 14c4762a1bSJed Brown PetscInt i,j,k,M,N,dof; 15c4762a1bSJed Brown 16*5f80ce2aSJacob Faibussowitsch CHKERRQ(DMDAGetInfo(da,0,&M,&N,0,0,0,0,&dof,0,0,0,0,0)); 17c4762a1bSJed Brown { 18c4762a1bSJed Brown struct {PetscScalar inside[dof];} **mystruct; 19*5f80ce2aSJacob Faibussowitsch CHKERRQ(DMDAVecGetArrayRead(da,global,(void*) &mystruct)); 20c4762a1bSJed Brown for (i=0; i<N; i++) { 21c4762a1bSJed Brown for (j=0; j<M; j++) { 22c4762a1bSJed Brown for (k=0; k<dof; k++) { 23*5f80ce2aSJacob Faibussowitsch CHKERRQ(PetscPrintf(PETSC_COMM_WORLD,"%D %D %g\n",i,j,(double)mystruct[i][j].inside[0])); 24c4762a1bSJed Brown 25c4762a1bSJed Brown mystruct[i][j].inside[1] = 2.1; 26c4762a1bSJed Brown } 27c4762a1bSJed Brown } 28c4762a1bSJed Brown } 29*5f80ce2aSJacob Faibussowitsch CHKERRQ(DMDAVecRestoreArrayRead(da,global,(void*) &mystruct)); 30c4762a1bSJed Brown } 31c4762a1bSJed Brown PetscFunctionReturn(0); 32c4762a1bSJed Brown } 33c4762a1bSJed Brown 34c4762a1bSJed Brown int main(int argc,char **argv) 35c4762a1bSJed Brown { 36c4762a1bSJed Brown PetscInt dof = 2,M = 3,N = 3,m = PETSC_DECIDE,n = PETSC_DECIDE; 37c4762a1bSJed Brown PetscErrorCode ierr; 38c4762a1bSJed Brown DM da; 39c4762a1bSJed Brown Vec global,local; 40c4762a1bSJed Brown 41c4762a1bSJed Brown ierr = PetscInitialize(&argc,&argv,(char*)0,help);if (ierr) return ierr; 42*5f80ce2aSJacob Faibussowitsch CHKERRQ(PetscOptionsGetInt(NULL,0,"-dof",&dof,0)); 43c4762a1bSJed Brown /* Create distributed array and get vectors */ 44*5f80ce2aSJacob Faibussowitsch CHKERRQ(DMDACreate2d(PETSC_COMM_WORLD, DM_BOUNDARY_NONE, DM_BOUNDARY_NONE,DMDA_STENCIL_BOX,M,N,m,n,dof,1,NULL,NULL,&da)); 45*5f80ce2aSJacob Faibussowitsch CHKERRQ(DMSetFromOptions(da)); 46*5f80ce2aSJacob Faibussowitsch CHKERRQ(DMSetUp(da)); 47*5f80ce2aSJacob Faibussowitsch CHKERRQ(DMCreateGlobalVector(da,&global)); 48*5f80ce2aSJacob Faibussowitsch CHKERRQ(DMCreateLocalVector(da,&local)); 49c4762a1bSJed Brown 50*5f80ce2aSJacob Faibussowitsch CHKERRQ(doit(da,global)); 51c4762a1bSJed Brown 52*5f80ce2aSJacob Faibussowitsch CHKERRQ(VecView(global,0)); 53c4762a1bSJed Brown 54c4762a1bSJed Brown /* Free memory */ 55*5f80ce2aSJacob Faibussowitsch CHKERRQ(VecDestroy(&local)); 56*5f80ce2aSJacob Faibussowitsch CHKERRQ(VecDestroy(&global)); 57*5f80ce2aSJacob Faibussowitsch CHKERRQ(DMDestroy(&da)); 58c4762a1bSJed Brown ierr = PetscFinalize(); 59c4762a1bSJed Brown return ierr; 60c4762a1bSJed Brown } 61