27   PetscCallMPI(MPI_Comm_rank(PETSC_COMM_WORLD, &rank));  in main()
28   PetscCallMPI(MPI_Comm_size(PETSC_COMM_WORLD, &size));  in main()
56 …PetscCall(DMDACreate2d(PETSC_COMM_WORLD, bx, by, stencil_type, M, N, PETSC_DECIDE, PETSC_DECIDE, d…  in main()
58 …PetscCall(DMDACreate3d(PETSC_COMM_WORLD, bx, by, bz, stencil_type, M, N, P, PETSC_DECIDE, PETSC_DE…  in main()
60     PetscCall(DMDACreate1d(PETSC_COMM_WORLD, bx, M, dof, stencil_width, NULL, &da));  in main()
67   PetscCall(PetscRandomCreate(PETSC_COMM_WORLD, &rdm));  in main()
72     PetscCall(PetscViewerBinaryOpen(PETSC_COMM_WORLD, "temp", FILE_MODE_WRITE, &viewer));  in main()
75     PetscCall(PetscViewerHDF5Open(PETSC_COMM_WORLD, "temp", FILE_MODE_WRITE, &viewer));  in main()
77 …} else SETERRQ(PETSC_COMM_WORLD, PETSC_ERR_SUP, "Invalid Viewer : Run with -binary or -hdf5 option…  in main()
81   PetscCall(PetscPrintf(PETSC_COMM_WORLD, "Global vector written to temp file is \n"));  in main()
85 …PetscCall(DMDACreate2d(PETSC_COMM_WORLD, bx, by, stencil_type, M, N, PETSC_DECIDE, PETSC_DECIDE, d…  in main()
87 …PetscCall(DMDACreate3d(PETSC_COMM_WORLD, bx, by, bz, stencil_type, M, N, P, PETSC_DECIDE, PETSC_DE…  in main()
89     PetscCall(DMDACreate1d(PETSC_COMM_WORLD, bx, M, dof, stencil_width, NULL, &da2));  in main()
95     PetscCall(PetscViewerBinaryOpen(PETSC_COMM_WORLD, "temp", FILE_MODE_READ, &viewer));  in main()
98     PetscCall(PetscViewerHDF5Open(PETSC_COMM_WORLD, "temp", FILE_MODE_READ, &viewer));  in main()
100 …} else SETERRQ(PETSC_COMM_WORLD, PETSC_ERR_SUP, "Invalid Viewer : Run with -binary or -hdf5 option…  in main()
107   PetscCall(PetscPrintf(PETSC_COMM_WORLD, "Global vector read from temp file is \n"));  in main()
112 …PetscCall(PetscPrintf(PETSC_COMM_WORLD, "ex23: Norm of difference %g should be zero\n", (double)no…  in main()
113     PetscCall(PetscPrintf(PETSC_COMM_WORLD, "  Number of processors %d\n", size));  in main()
114 …PetscCall(PetscPrintf(PETSC_COMM_WORLD, "  M,N,P,dof %" PetscInt_FMT " %" PetscInt_FMT " %" PetscI…  in main()
115 …PetscCall(PetscPrintf(PETSC_COMM_WORLD, "  stencil_width %" PetscInt_FMT " stencil_type %d periodi…  in main()
116     PetscCall(PetscPrintf(PETSC_COMM_WORLD, "  dimension %d\n", 1 + (int)flg2 + (int)flg3));  in main()