Actual source code: ex31.c
1: static char help[] = "Tests MAIJ matrix for large DOF\n\n";
3: #include petscda.h
7: int main(int argc,char *argv[])
8: {
9: Mat M;
10: Vec x,y;
12: DA da,daf;
14: PetscInitialize(&argc,&argv,0,help);
15: DACreate2d(PETSC_COMM_WORLD,DA_NONPERIODIC,DA_STENCIL_STAR,4,5,PETSC_DECIDE,PETSC_DECIDE,41,1,0,0,&da);
16: DARefine(da,PETSC_COMM_WORLD,&daf);
17: DAGetInterpolation(da,daf,&M,PETSC_NULL);
18: DACreateGlobalVector(da,&x);
19: DACreateGlobalVector(daf,&y);
21: MatMult(M,x,y);
22: MatMultTranspose(M,y,x);
23: DADestroy(da);
24: DADestroy(daf);
25: VecDestroy(x);
26: VecDestroy(y);
27: MatDestroy(M);
28: PetscFinalize();
29: return 0;
30: }