159   PetscInt    i, j, Mx, My, xs, ys, xm, ym;  in FormIFunction()  local
166 …PetscCall(DMDAGetInfo(da, PETSC_IGNORE, &Mx, &My, PETSC_IGNORE, PETSC_IGNORE, PETSC_IGNORE, PETSC_…  in FormIFunction()
170   hy = 1.0 / (PetscReal)(My - 1);  in FormIFunction()
195       if (i == 0 || j == 0 || i == Mx - 1 || j == My - 1) {  in FormIFunction()
203           } else if (i == 0 && j == My - 1) { /* NW corner */  in FormIFunction()
205           } else if (i == Mx - 1 && j == My - 1) { /* NE corner */  in FormIFunction()
213           } else if (j == My - 1) { /* Top */  in FormIFunction()
248   PetscInt    i, j, Mx, My, xs, ys, xm, ym, nc;  in FormIJacobian()  local
255 …PetscCall(DMDAGetInfo(da, PETSC_IGNORE, &Mx, &My, PETSC_IGNORE, PETSC_IGNORE, PETSC_IGNORE, PETSC_…  in FormIJacobian()
260   hy = 1.0 / (PetscReal)(My - 1);  in FormIJacobian()
268       if (user->boundary == 0 && (i == 0 || i == Mx - 1 || j == 0 || j == My - 1)) {  in FormIJacobian()
294       } else if (user->boundary > 0 && j == My - 1) { /* Top Neumann */  in FormIJacobian()
341   PetscInt      i, j, xs, ys, xm, ym, Mx, My;  in FormInitialSolution()  local
346 …PetscCall(DMDAGetInfo(da, PETSC_IGNORE, &Mx, &My, PETSC_IGNORE, PETSC_IGNORE, PETSC_IGNORE, PETSC_…  in FormInitialSolution()
349   hy = 1.0 / (PetscReal)(My - 1);  in FormInitialSolution()