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 PetscInt i,j,k,M,N,dof; 14c4762a1bSJed Brown 159566063dSJacob Faibussowitsch PetscCall(DMDAGetInfo(da,0,&M,&N,0,0,0,0,&dof,0,0,0,0,0)); 16c4762a1bSJed Brown { 17c4762a1bSJed Brown struct {PetscScalar inside[dof];} **mystruct; 189566063dSJacob Faibussowitsch PetscCall(DMDAVecGetArrayRead(da,global,(void*) &mystruct)); 19c4762a1bSJed Brown for (i=0; i<N; i++) { 20c4762a1bSJed Brown for (j=0; j<M; j++) { 21c4762a1bSJed Brown for (k=0; k<dof; k++) { 22*63a3b9bcSJacob Faibussowitsch PetscCall(PetscPrintf(PETSC_COMM_WORLD,"%" PetscInt_FMT " %" PetscInt_FMT " %g\n",i,j,(double)mystruct[i][j].inside[0])); 23c4762a1bSJed Brown 24c4762a1bSJed Brown mystruct[i][j].inside[1] = 2.1; 25c4762a1bSJed Brown } 26c4762a1bSJed Brown } 27c4762a1bSJed Brown } 289566063dSJacob Faibussowitsch PetscCall(DMDAVecRestoreArrayRead(da,global,(void*) &mystruct)); 29c4762a1bSJed Brown } 30c4762a1bSJed Brown PetscFunctionReturn(0); 31c4762a1bSJed Brown } 32c4762a1bSJed Brown 33c4762a1bSJed Brown int main(int argc,char **argv) 34c4762a1bSJed Brown { 35c4762a1bSJed Brown PetscInt dof = 2,M = 3,N = 3,m = PETSC_DECIDE,n = PETSC_DECIDE; 36c4762a1bSJed Brown DM da; 37c4762a1bSJed Brown Vec global,local; 38c4762a1bSJed Brown 399566063dSJacob Faibussowitsch PetscCall(PetscInitialize(&argc,&argv,(char*)0,help)); 409566063dSJacob Faibussowitsch PetscCall(PetscOptionsGetInt(NULL,0,"-dof",&dof,0)); 41c4762a1bSJed Brown /* Create distributed array and get vectors */ 429566063dSJacob Faibussowitsch PetscCall(DMDACreate2d(PETSC_COMM_WORLD, DM_BOUNDARY_NONE, DM_BOUNDARY_NONE,DMDA_STENCIL_BOX,M,N,m,n,dof,1,NULL,NULL,&da)); 439566063dSJacob Faibussowitsch PetscCall(DMSetFromOptions(da)); 449566063dSJacob Faibussowitsch PetscCall(DMSetUp(da)); 459566063dSJacob Faibussowitsch PetscCall(DMCreateGlobalVector(da,&global)); 469566063dSJacob Faibussowitsch PetscCall(DMCreateLocalVector(da,&local)); 47c4762a1bSJed Brown 489566063dSJacob Faibussowitsch PetscCall(doit(da,global)); 49c4762a1bSJed Brown 509566063dSJacob Faibussowitsch PetscCall(VecView(global,0)); 51c4762a1bSJed Brown 52c4762a1bSJed Brown /* Free memory */ 539566063dSJacob Faibussowitsch PetscCall(VecDestroy(&local)); 549566063dSJacob Faibussowitsch PetscCall(VecDestroy(&global)); 559566063dSJacob Faibussowitsch PetscCall(DMDestroy(&da)); 569566063dSJacob Faibussowitsch PetscCall(PetscFinalize()); 57b122ec5aSJacob Faibussowitsch return 0; 58c4762a1bSJed Brown } 59