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: }