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 165f80ce2aSJacob 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; 195f80ce2aSJacob 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++) { 235f80ce2aSJacob 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 } 295f80ce2aSJacob 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 41*b122ec5aSJacob Faibussowitsch CHKERRQ(PetscInitialize(&argc,&argv,(char*)0,help)); 425f80ce2aSJacob Faibussowitsch CHKERRQ(PetscOptionsGetInt(NULL,0,"-dof",&dof,0)); 43c4762a1bSJed Brown /* Create distributed array and get vectors */ 445f80ce2aSJacob Faibussowitsch CHKERRQ(DMDACreate2d(PETSC_COMM_WORLD, DM_BOUNDARY_NONE, DM_BOUNDARY_NONE,DMDA_STENCIL_BOX,M,N,m,n,dof,1,NULL,NULL,&da)); 455f80ce2aSJacob Faibussowitsch CHKERRQ(DMSetFromOptions(da)); 465f80ce2aSJacob Faibussowitsch CHKERRQ(DMSetUp(da)); 475f80ce2aSJacob Faibussowitsch CHKERRQ(DMCreateGlobalVector(da,&global)); 485f80ce2aSJacob Faibussowitsch CHKERRQ(DMCreateLocalVector(da,&local)); 49c4762a1bSJed Brown 505f80ce2aSJacob Faibussowitsch CHKERRQ(doit(da,global)); 51c4762a1bSJed Brown 525f80ce2aSJacob Faibussowitsch CHKERRQ(VecView(global,0)); 53c4762a1bSJed Brown 54c4762a1bSJed Brown /* Free memory */ 555f80ce2aSJacob Faibussowitsch CHKERRQ(VecDestroy(&local)); 565f80ce2aSJacob Faibussowitsch CHKERRQ(VecDestroy(&global)); 575f80ce2aSJacob Faibussowitsch CHKERRQ(DMDestroy(&da)); 58*b122ec5aSJacob Faibussowitsch CHKERRQ(PetscFinalize()); 59*b122ec5aSJacob Faibussowitsch return 0; 60c4762a1bSJed Brown } 61