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()
254       if (i + 1 == gxs + gxm) { /* right 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()
327       d3 = (bottom[i - xs + 1] - x[i - gxs]) * rhy;  in FormFunctionGradient()
341       d1 = (x[(gym - 1) * gxm + i - gxs] - top[i - xs + 1]) * rhy;  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()
463       if (i == gxs) {  in FormHessian()
477       if (i + 1 == gxs + gxm) {  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()