74   DMDALocalInfo info;  in main()  local
105   PetscCall(DMDAGetLocalInfo(da_after, &info));  in main()
107   PetscCall(FormExactSolution(&info, u_exact));  in main()
110   error1 /= (PetscReal)info.mx * (PetscReal)info.my; /* average error */  in main()
112 …tscInt_FMT " grid:  av |u-uexact|  = %.3e,  |u-uexact|_inf = %.3e\n", info.mx, info.my, (double)er…  in main()
120 PetscErrorCode FormExactSolution(DMDALocalInfo *info, Vec u)  in FormExactSolution()  argument
126   dx = 4.0 / (PetscReal)(info->mx - 1);  in FormExactSolution()
127   dy = 4.0 / (PetscReal)(info->my - 1);  in FormExactSolution()
128   PetscCall(DMDAVecGetArray(info->da, u, &au));  in FormExactSolution()
129   for (j = info->ys; j < info->ys + info->ym; j++) {  in FormExactSolution()
131     for (i = info->xs; i < info->xs + info->xm; i++) {  in FormExactSolution()
136   PetscCall(DMDAVecRestoreArray(info->da, u, &au));  in FormExactSolution()
143   DMDALocalInfo info;  in FormBounds()  local
149   PetscCall(DMDAGetLocalInfo(da, &info));  in FormBounds()
150   dx = 4.0 / (PetscReal)(info.mx - 1);  in FormBounds()
151   dy = 4.0 / (PetscReal)(info.my - 1);  in FormBounds()
153   for (j = info.ys; j < info.ys + info.ym; j++) {  in FormBounds()
155     for (i = info.xs; i < info.xs + info.xm; i++) {  in FormBounds()
165 PetscErrorCode FormFunctionLocal(DMDALocalInfo *info, PetscScalar **au, PetscScalar **af, void *use…  in FormFunctionLocal()  argument
171   dx = 4.0 / (PetscReal)(info->mx - 1);  in FormFunctionLocal()
172   dy = 4.0 / (PetscReal)(info->my - 1);  in FormFunctionLocal()
173   for (j = info->ys; j < info->ys + info->ym; j++) {  in FormFunctionLocal()
175     for (i = info->xs; i < info->xs + info->xm; i++) {  in FormFunctionLocal()
177       if (i == 0 || j == 0 || i == info->mx - 1 || j == info->my - 1) {  in FormFunctionLocal()
181         ue       = (i + 1 == info->mx - 1) ? u_exact(x + dx, y) : au[j][i + 1];  in FormFunctionLocal()
183         un       = (j + 1 == info->my - 1) ? u_exact(x, y + dy) : au[j + 1][i];  in FormFunctionLocal()
188   PetscCall(PetscLogFlops(12.0 * info->ym * info->xm));  in FormFunctionLocal()
192 PetscErrorCode FormJacobianLocal(DMDALocalInfo *info, PetscScalar **au, Mat A, Mat jac, void *user)  in FormJacobianLocal()  argument
199   dx  = 4.0 / (PetscReal)(info->mx - 1);  in FormJacobianLocal()
200   dy  = 4.0 / (PetscReal)(info->my - 1);  in FormJacobianLocal()
203   for (j = info->ys; j < info->ys + info->ym; j++) {  in FormJacobianLocal()
204     for (i = info->xs; i < info->xs + info->xm; i++) {  in FormJacobianLocal()
207       if (i == 0 || j == 0 || i == info->mx - 1 || j == info->my - 1) { /* boundary */  in FormJacobianLocal()
220         if (i + 1 < info->mx - 1) {  in FormJacobianLocal()
230         if (j + 1 < info->my - 1) {  in FormJacobianLocal()
247   PetscCall(PetscLogFlops(2.0 * info->ym * info->xm));  in FormJacobianLocal()