202   PetscInt   xs, xm, gxs, gxm, ys, ym, gys, gym;  in FormFunctionGradient()  local
214   PetscCall(DMDAGetCorners(user->dm, &xs, &ys, NULL, &xm, &ym, NULL));  in FormFunctionGradient()
233   for (j = ys; j < ys + ym; j++) {  in FormFunctionGradient()
241         xl  = left[j - ys + 1];  in FormFunctionGradient()
242         xlt = left[j - ys + 2];  in FormFunctionGradient()
255         xr  = right[j - ys + 1];  in FormFunctionGradient()
256         xrb = right[j - ys];  in FormFunctionGradient()
318     for (j = ys; j < ys + ym; j++) {  in FormFunctionGradient()
319       d3 = (left[j - ys + 1] - left[j - ys + 2]) * rhy;  in FormFunctionGradient()
320       d2 = (left[j - ys + 1] - x[(j - gys) * gxm]) * rhx;  in FormFunctionGradient()
324   if (ys == 0) { /* bottom side */  in FormFunctionGradient()
333     for (j = ys; j < ys + ym; j++) {  in FormFunctionGradient()
334       d1 = (x[(j + 1 - gys) * gxm - 1] - right[j - ys + 1]) * rhx;  in FormFunctionGradient()
335       d4 = (right[j - ys] - right[j - ys + 1]) * rhy;  in FormFunctionGradient()
339   if (ys + ym == my) { /* top side */  in FormFunctionGradient()
347   if (ys == 0 && xs == 0) {  in FormFunctionGradient()
352   if (ys + ym == my && xs + xm == mx) {  in FormFunctionGradient()
420   PetscInt   xs, xm, gxs, gxm, ys, ym, gys, gym, col[7];  in FormHessian()  local
439   PetscCall(DMDAGetCorners(user->dm, &xs, &ys, NULL, &xm, &ym, NULL));  in FormHessian()
456     for (j = ys; j < ys + ym; j++) {  in FormHessian()
464         xl  = left[j - ys + 1];  in FormHessian()
465         xlt = left[j - ys + 2];  in FormHessian()
478         xr  = right[j - ys + 1];  in FormHessian()
479         xrb = right[j - ys];  in FormHessian()
606   PetscInt   xs, ys, xm, ym, gxs, gys, gxm, gym;  in MSA_BoundaryConditions()  local
619   PetscCall(DMDAGetCorners(user->dm, &xs, &ys, NULL, &xm, &ym, NULL));  in MSA_BoundaryConditions()
652       yt    = b + hy * ys;  in MSA_BoundaryConditions()
657       yt    = b + hy * ys;  in MSA_BoundaryConditions()
726   PetscInt   xs, ys, xm, ym;  in MSA_Plate()  local
740   PetscCall(DMDAGetCorners(user->dm, &xs, &ys, NULL, &xm, &ym, NULL));  in MSA_Plate()
751       for (j = ys; j < ys + ym; j++) {  in MSA_Plate()
752         row = (j - ys) * xm + (i - xs);  in MSA_Plate()
762       for (j = ys; j < ys + ym; j++) {  in MSA_Plate()
763         row = (j - ys) * xm + (i - xs);  in MSA_Plate()
804     PetscInt   row, xs, xm, gxs, gxm, ys, ym, gys, gym;  in MSA_InitialPoint()  local
810     PetscCall(DMDAGetCorners(user->dm, &xs, &ys, NULL, &xm, &ym, NULL));  in MSA_InitialPoint()
821     for (j = ys; j < ys + ym; j++) {  in MSA_InitialPoint()
824 …+ 1) * top[i - xs + 1] / (my + 2) + (i + 1) * left[j - ys + 1] / mx + (mx - i + 1) * right[j - ys …  in MSA_InitialPoint()