166   PetscInt    xs, xm, gxs, gxm, ys, ym, gys, gym;  in FormFunction()  local
177   PetscCall(DMDAGetCorners(user->dm, &xs, &ys, NULL, &xm, &ym, NULL));  in FormFunction()
188   for (j = ys; j < ys + ym; j++) {  in FormFunction()
193         xl = user->left[j - ys + 1];  in FormFunction()
205         xr = user->right[j - ys + 1];  in FormFunction()
235     for (j = ys; j < ys + ym; j++) {  in FormFunction()
236       d3 = (user->left[j - ys + 1] - user->left[j - ys + 2]) * rhy;  in FormFunction()
237       d2 = (user->left[j - ys + 1] - x[j][0]) * rhx;  in FormFunction()
241   if (ys == 0) { /* bottom side */  in FormFunction()
249     for (j = ys; j < ys + ym; j++) {  in FormFunction()
250       d1 = (x[j][mx - 1] - user->right[j - ys + 1]) * rhx;  in FormFunction()
251       d4 = (user->right[j - ys] - user->right[j - ys + 1]) * rhy;  in FormFunction()
255   if (ys + ym == my) { /* top side */  in FormFunction()
262   if (ys == 0 && xs == 0) {  in FormFunction()
267   if (ys + ym == my && xs + xm == mx) {  in FormFunction()
300   PetscInt    xs, xm, gxs, gxm, ys, ym, gys, gym;  in FormFunctionGradient()  local
312   PetscCall(DMDAGetCorners(user->dm, &xs, &ys, NULL, &xm, &ym, NULL));  in FormFunctionGradient()
324   for (j = ys; j < ys + ym; j++) {  in FormFunctionGradient()
330         xl  = user->left[j - ys + 1];  in FormFunctionGradient()
331         xlt = user->left[j - ys + 2];  in FormFunctionGradient()
344         xr  = user->right[j - ys + 1];  in FormFunctionGradient()
345         xrb = user->right[j - ys];  in FormFunctionGradient()
407     for (j = ys; j < ys + ym; j++) {  in FormFunctionGradient()
408       d3 = (user->left[j - ys + 1] - user->left[j - ys + 2]) * rhy;  in FormFunctionGradient()
409       d2 = (user->left[j - ys + 1] - x[j][0]) * rhx;  in FormFunctionGradient()
413   if (ys == 0) { /* bottom side */  in FormFunctionGradient()
422     for (j = ys; j < ys + ym; j++) {  in FormFunctionGradient()
423       d1 = (x[j][mx - 1] - user->right[j - ys + 1]) * rhx;  in FormFunctionGradient()
424       d4 = (user->right[j - ys] - user->right[j - ys + 1]) * rhy;  in FormFunctionGradient()
428   if (ys + ym == my) { /* top side */  in FormFunctionGradient()
436   if (ys == 0 && xs == 0) {  in FormFunctionGradient()
441   if (ys + ym == my && xs + xm == mx) {  in FormFunctionGradient()
506   PetscInt    xs, xm, gxs, gxm, ys, ym, gys, gym;  in QuadraticH()  local
519   PetscCall(DMDAGetCorners(user->dm, &xs, &ys, NULL, &xm, &ym, NULL));  in QuadraticH()
537   for (j = ys; j < ys + ym; j++) {  in QuadraticH()
544         xl  = user->left[j - ys + 1];  in QuadraticH()
545         xlt = user->left[j - ys + 2];  in QuadraticH()
558         xr  = user->right[j - ys + 1];  in QuadraticH()
559         xrb = user->right[j - ys];  in QuadraticH()
690   PetscInt   xs, ys, xm, ym, gxs, gys, gxm, gym;  in MSA_BoundaryConditions()  local
702   PetscCall(DMDAGetCorners(user->dm, &xs, &ys, NULL, &xm, &ym, NULL));  in MSA_BoundaryConditions()
730       yt       = b + hy * ys;  in MSA_BoundaryConditions()
735       yt       = b + hy * ys;  in MSA_BoundaryConditions()
826     PetscInt    xs, xm, ys, ym;  in MSA_InitialPoint()  local
831     PetscCall(DMDAGetCorners(user->dm, &xs, &ys, NULL, &xm, &ym, NULL));  in MSA_InitialPoint()
837     for (j = ys; j < ys + ym; j++) {  in MSA_InitialPoint()
838 …op[i - xs + 1]) / (my + 2) + ((i + 1) * user->left[j - ys + 1] + (mx - i + 1) * user->right[j - ys…  in MSA_InitialPoint()