Printing the solution
229: /* Draw solution grid */
233: PetscOptionsHasName(PETSC_NULL,"–view_sol_serial",&flg);
234: if (flg) {VecView( x, PETSC_VIEWER_STDOUT_WORLD); }
236: PetscOptionsHasName(PETSC_NULL,"–view_sol",&flg);
241: PetscPrintf( PETSC_COMM_WORLD, "Solution Grid (without boundary conditions):\n" );
242: for (I=Istart; I<Iend; I++) {
243: i = I/n; j = I – i*n;