Actual source code: ex13.c
2: static char help[] = "Tests loading DA vector from file.\n\n";
4: #include petscda.h
8: int main(int argc,char **argv)
9: {
11: PetscInt M = PETSC_DECIDE,N = PETSC_DECIDE;
12: DA da;
13: Vec global;
14: PetscViewer bviewer;
16: PetscInitialize(&argc,&argv,(char*)0,help);
18: /* Read options */
19: PetscOptionsGetInt(PETSC_NULL,"-M",&M,PETSC_NULL);
20: PetscOptionsGetInt(PETSC_NULL,"-N",&N,PETSC_NULL);
22: PetscViewerBinaryOpen(PETSC_COMM_WORLD,"daoutput",FILE_MODE_READ,&bviewer);
23: DALoad(bviewer,M,N,PETSC_DECIDE,&da);
24: DACreateGlobalVector(da,&global);
25: VecLoadIntoVector(bviewer,global);
26: PetscViewerDestroy(bviewer);
29: VecView(global,PETSC_VIEWER_DRAW_WORLD);
32: /* Free memory */
33: VecDestroy(global);
34: DADestroy(da);
35: PetscFinalize();
36: return 0;
37: }
38: