1 static char help[] = "Tests DMDACreate3d() memory usage\n\n"; 2 3 #include <petscdmda.h> 4 5 6 int main(int argc,char **argv) 7 { 8 DM dm; 9 Vec X,Y; 10 PetscErrorCode ierr; 11 PetscInt dof = 10; 12 13 ierr = PetscInitialize(&argc,&argv,(char*)0,help);if (ierr) return ierr; 14 ierr = PetscOptionsGetInt(NULL,NULL,"-dof",&dof,NULL);CHKERRQ(ierr); 15 ierr = DMDACreate3d(PETSC_COMM_WORLD,DM_BOUNDARY_NONE,DM_BOUNDARY_NONE,DM_BOUNDARY_NONE,DMDA_STENCIL_STAR,-128,-128,-128,PETSC_DECIDE,PETSC_DECIDE,PETSC_DECIDE,dof,1,NULL,NULL,NULL,&dm);CHKERRQ(ierr); 16 ierr = DMSetFromOptions(dm);CHKERRQ(ierr); 17 ierr = DMSetUp(dm);CHKERRQ(ierr); 18 ierr = PetscMemoryTrace("DMDACreate3d ");CHKERRQ(ierr); 19 20 ierr = DMCreateGlobalVector(dm,&X);CHKERRQ(ierr); 21 ierr = PetscMemoryTrace("DMCreateGlobalVector"); 22 ierr = DMCreateGlobalVector(dm,&Y);CHKERRQ(ierr); 23 ierr = PetscMemoryTrace("DMCreateGlobalVector");CHKERRQ(ierr); 24 25 ierr = VecDestroy(&X);CHKERRQ(ierr); 26 ierr = VecDestroy(&Y);CHKERRQ(ierr); 27 ierr = DMDestroy(&dm);CHKERRQ(ierr); 28 ierr = PetscFinalize(); 29 return ierr; 30 } 31 32