1 static char help[] = "Tests MAIJ matrix for large DOF\n\n"; 2 3 #include <petscdm.h> 4 #include <petscdmda.h> 5 6 int main(int argc,char *argv[]) 7 { 8 Mat M; 9 Vec x,y; 10 PetscErrorCode ierr; 11 DM da,daf; 12 13 ierr = PetscInitialize(&argc,&argv,0,help);if (ierr) return ierr; 14 ierr = DMDACreate2d(PETSC_COMM_WORLD, DM_BOUNDARY_NONE, DM_BOUNDARY_NONE,DMDA_STENCIL_STAR,4,5,PETSC_DECIDE,PETSC_DECIDE,41,1,0,0,&da);CHKERRQ(ierr); 15 ierr = DMSetFromOptions(da);CHKERRQ(ierr); 16 ierr = DMSetUp(da);CHKERRQ(ierr); 17 ierr = DMRefine(da,PETSC_COMM_WORLD,&daf);CHKERRQ(ierr); 18 ierr = DMCreateInterpolation(da,daf,&M,NULL);CHKERRQ(ierr); 19 ierr = DMCreateGlobalVector(da,&x);CHKERRQ(ierr); 20 ierr = DMCreateGlobalVector(daf,&y);CHKERRQ(ierr); 21 22 ierr = MatMult(M,x,y);CHKERRQ(ierr); 23 ierr = MatMultTranspose(M,y,x);CHKERRQ(ierr); 24 ierr = DMDestroy(&da);CHKERRQ(ierr); 25 ierr = DMDestroy(&daf);CHKERRQ(ierr); 26 ierr = VecDestroy(&x);CHKERRQ(ierr); 27 ierr = VecDestroy(&y);CHKERRQ(ierr); 28 ierr = MatDestroy(&M);CHKERRQ(ierr); 29 ierr = PetscFinalize(); 30 return ierr; 31 } 32 33 34 /*TEST 35 36 test: 37 nsize: 2 38 39 TEST*/ 40