202   PetscInt   xs, xm, gxs, gxm, ys, ym, gys, gym;  in FormFunctionGradient()  local
215   PetscCall(DMDAGetGhostCorners(user->dm, &gxs, &gys, NULL, &gxm, &gym, NULL));  in FormFunctionGradient()
235       row = (j - gys) * gxm + (i - gxs);  in FormFunctionGradient()
261       if (j + 1 == gys + gym) { /* top side */  in FormFunctionGradient()
268       if (i > gxs && j + 1 < gys + gym) xlt = x[row - 1 + gxm];  in FormFunctionGradient()
269       if (j > gys && i + 1 < gxs + gxm) xrb = x[row + 1 - gxm];  in FormFunctionGradient()
320       d2 = (left[j - ys + 1] - x[(j - gys) * gxm]) * rhx;  in FormFunctionGradient()
334       d1 = (x[(j + 1 - gys) * gxm - 1] - right[j - ys + 1]) * rhx;  in FormFunctionGradient()
420   PetscInt   xs, xm, gxs, gxm, ys, ym, gys, gym, col[7];  in FormHessian()  local
440   PetscCall(DMDAGetGhostCorners(user->dm, &gxs, &gys, NULL, &gxm, &gym, NULL));  in FormHessian()
457       row = (j - gys) * gxm + (i - gxs);  in FormHessian()
470       if (j == gys) {  in FormHessian()
484       if (j + 1 == gys + gym) {  in FormHessian()
491       if (i > gxs && j + 1 < gys + gym) xlt = x[row - 1 + gxm];  in FormHessian()
492       if (j > gys && i + 1 < gxs + gxm) xrb = x[row + 1 - gxm];  in FormHessian()
606   PetscInt   xs, ys, xm, ym, gxs, gys, gxm, gym;  in MSA_BoundaryConditions()  local
620   PetscCall(DMDAGetGhostCorners(user->dm, &gxs, &gys, NULL, &gxm, &gym, NULL));  in MSA_BoundaryConditions()
804     PetscInt   row, xs, xm, gxs, gxm, ys, ym, gys, gym;  in MSA_InitialPoint()  local
811     PetscCall(DMDAGetGhostCorners(user->dm, &gxs, &gys, NULL, &gxm, &gym, NULL));  in MSA_InitialPoint()
823         row    = (j - gys) * gxm + (i - gxs);  in MSA_InitialPoint()