Printing the solution via X Window

29: #include petscda.h

251: PetscOptionsHasName(PETSC_NULL,"–view_sol_x",&flg);

252: if (flg) {

253: /* view solution grid in an X window */

254: PetscScalar *xx; DA da;

int size; Vec x_da;

258: MPI_Comm_size(PETSC_COMM_WORLD,&size);

259: /* force 1-dimensional distribution, same as with x !!! */

260: DACreate2d(PETSC_COMM_WORLD,DA_NONPERIODIC,DA_STENCIL_STAR,

261: n,m,1,size,1,0,PETSC_NULL,PETSC_NULL,&da);

262: DACreateGlobalVector(da, &x_da);

263: VecGetOwnershipRange(x, &Istart, &Iend);

264: VecGetArray(x, &xx);

265: for (I=Istart; I<Iend; I++)

266: VecSetValues(x_da, 1, &I, &xx[I-Istart], INSERT_VALUES);

267: VecRestoreArray(x, &xx);

268: VecAssemblyBegin(x_da); VecAssemblyEnd(x_da);

269: PetscOptionsHasName(PETSC_NULL,"-view_sol_x_da",&flg);

270: if (flg) VecView(x_da,PETSC_VIEWER_STDOUT_WORLD);

271: VecView(x_da,PETSC_VIEWER_DRAW_(PETSC_COMM_WORLD));

272: VecDestroy(x_da);  DADestroy(da);

273: }

Create a 2-dimensional vector x_da

Copy values of x into x_da

View the solution as a 2-dimensional plot

Controling thedistribution of x_da

Vorherige Folie

Nächste Folie

Zurück zur ersten Folie

Graphik-Version anzeigen