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