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()
234     for (i = xs; i < xs + xm; i++) {  in FormFunctionGradient()
248         xb  = bottom[i - xs + 1];  in FormFunctionGradient()
249         xrb = bottom[i - xs + 2];  in FormFunctionGradient()
262         xt  = top[i - xs + 1];  in FormFunctionGradient()
263         xlt = top[i - xs];  in FormFunctionGradient()
317   if (xs == 0) { /* left side */  in FormFunctionGradient()
325     for (i = xs; i < xs + xm; i++) {  in FormFunctionGradient()
326       d2 = (bottom[i + 1 - xs] - bottom[i - xs + 2]) * rhx;  in FormFunctionGradient()
327       d3 = (bottom[i - xs + 1] - x[i - gxs]) * rhy;  in FormFunctionGradient()
332   if (xs + xm == mx) { /* right side */  in FormFunctionGradient()
340     for (i = xs; i < xs + xm; i++) {  in FormFunctionGradient()
341       d1 = (x[(gym - 1) * gxm + i - gxs] - top[i - xs + 1]) * rhy;  in FormFunctionGradient()
342       d4 = (top[i - xs + 1] - top[i - xs]) * rhx;  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()
455   for (i = xs; i < xs + xm; i++) {  in FormHessian()
471         xb  = bottom[i - xs + 1];  in FormHessian()
472         xrb = bottom[i - xs + 2];  in FormHessian()
485         xt  = top[i - xs + 1];  in FormHessian()
486         xlt = top[i - xs];  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()
643       xt    = l + hx * xs;  in MSA_BoundaryConditions()
648       xt    = l + hx * xs;  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()
750     for (i = xs; i < xs + xm; i++) {  in MSA_Plate()
752         row = (j - ys) * xm + (i - xs);  in MSA_Plate()
761     for (i = xs; i < xs + xm; i++) {  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()
822       for (i = xs; i < xs + xm; i++) {  in MSA_InitialPoint()
824 …x[row] = ((j + 1) * bottom[i - xs + 1] / my + (my - j + 1) * top[i - xs + 1] / (my + 2) + (i + 1) …  in MSA_InitialPoint()