22   PetscInt   mx, my;                      /* discretization in x, y directions */  member
53   user.mx = 4;  in main()
57   PetscCall(PetscOptionsGetInt(NULL, NULL, "-mx", &user.mx, &flg));  in main()
61 …intf(PETSC_COMM_SELF, "mx: %" PetscInt_FMT "     my: %" PetscInt_FMT "   \n\n", user.mx, user.my));  in main()
64   N = user.mx * user.my;  in main()
94 …if (ksp) PetscCall(KSPSetTolerances(ksp, PETSC_CURRENT, PETSC_CURRENT, PETSC_CURRENT, user.mx * us…  in main()
128   PetscInt           mx = user->mx, my = user->my;  in FormFunctionGradient()  local
129   PetscReal          rhx = mx + 1, rhy = my + 1;  in FormFunctionGradient()
130 …PetscReal          hx = 1.0 / (mx + 1), hy = 1.0 / (my + 1), hydhx = hy / hx, hxdhy = hx / hy, are…  in FormFunctionGradient()
145     for (i = 0; i < mx; i++) {  in FormFunctionGradient()
146       row = (j)*mx + (i);  in FormFunctionGradient()
160         xb = x[row - mx];  in FormFunctionGradient()
163       if (i + 1 == mx) { /* right side */  in FormFunctionGradient()
174         xt = x[row + mx];  in FormFunctionGradient()
177       if (i > 0 && j + 1 < my) xlt = x[row - 1 + mx];  in FormFunctionGradient()
178       if (j > 0 && i + 1 < mx) xrb = x[row + 1 - mx];  in FormFunctionGradient()
227     d2 = (user->left[j + 1] - x[j * mx]) * rhx;  in FormFunctionGradient()
231   for (i = 0; i < mx; i++) { /* bottom */  in FormFunctionGradient()
238     d1 = (x[(j + 1) * mx - 1] - user->right[j + 1]) * rhx;  in FormFunctionGradient()
243   for (i = 0; i < mx; i++) { /* top side */  in FormFunctionGradient()
244     d1 = (x[(my - 1) * mx + i] - user->top[i + 1]) * rhy;  in FormFunctionGradient()
256   d2 = (user->top[mx + 1] - user->top[mx]) * rhx;  in FormFunctionGradient()
264   PetscCall(PetscLogFlops(67.0 * mx * my));  in FormFunctionGradient()
306   PetscInt           mx = user->mx, my = user->my;  in QuadraticH()  local
308   PetscReal          hx = 1.0 / (mx + 1), hy = 1.0 / (my + 1), hydhx = hy / hx, hxdhy = hx / hy;  in QuadraticH()
309   PetscReal          rhx = mx + 1, rhy = my + 1;  in QuadraticH()
326   for (i = 0; i < mx; i++) {  in QuadraticH()
328       row = (j)*mx + (i);  in QuadraticH()
345         xb = x[row - mx];  in QuadraticH()
348       if (i + 1 == mx) {  in QuadraticH()
359         xt = x[row + mx];  in QuadraticH()
362       if (i > 0 && j + 1 < my) xlt = x[row - 1 + mx];  in QuadraticH()
363       if (j > 0 && i + 1 < mx) xrb = x[row + 1 - mx];  in QuadraticH()
402         col[k] = row - mx;  in QuadraticH()
406       if (j > 0 && i < mx - 1) {  in QuadraticH()
408         col[k] = row - mx + 1;  in QuadraticH()
422       if (i < mx - 1) {  in QuadraticH()
430         col[k] = row + mx - 1;  in QuadraticH()
436         col[k] = row + mx;  in QuadraticH()
455   PetscCall(PetscLogFlops(199.0 * mx * my));  in QuadraticH()
474   PetscInt   mx = user->mx, my = user->my;  in MSA_BoundaryConditions()  local
483   bsize = mx + 2;  in MSA_BoundaryConditions()
486   tsize = mx + 2;  in MSA_BoundaryConditions()
493   hx = (r - l) / (mx + 1);  in MSA_BoundaryConditions()
572     PetscInt   mx = user->mx, my = user->my;  in MSA_InitialPoint()  local
579       for (i = 0; i < mx; i++) {  in MSA_InitialPoint()
580         row    = (j)*mx + (i);  in MSA_InitialPoint()
581 …top[i + 1]) / (my + 2) + ((i + 1) * user->left[j + 1] + (mx - i + 1) * user->right[j + 1]) / (mx +…  in MSA_InitialPoint()