22   PetscInt   mx, my;                      /* discretization in x, y directions */  member
54   user.my = 4;  in main()
58   PetscCall(PetscOptionsGetInt(NULL, NULL, "-my", &user.my, &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 …) PetscCall(KSPSetTolerances(ksp, PETSC_CURRENT, PETSC_CURRENT, PETSC_CURRENT, user.mx * user.my));  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()
144   for (j = 0; j < my; j++) {  in FormFunctionGradient()
170       if (j + 1 == 0 + my) { /* top side */  in FormFunctionGradient()
177       if (i > 0 && j + 1 < my) xlt = x[row - 1 + mx];  in FormFunctionGradient()
225   for (j = 0; j < my; j++) { /* left side */  in FormFunctionGradient()
237   for (j = 0; j < my; j++) { /* right side */  in FormFunctionGradient()
244     d1 = (x[(my - 1) * mx + i] - user->top[i + 1]) * rhy;  in FormFunctionGradient()
255   d1 = (user->right[my + 1] - user->right[my]) * rhy;  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()
327     for (j = 0; j < my; j++) {  in QuadraticH()
355       if (j + 1 == my) {  in QuadraticH()
362       if (i > 0 && j + 1 < my) xlt = x[row - 1 + mx];  in QuadraticH()
428       if (i > 0 && j < my - 1) {  in QuadraticH()
434       if (j < my - 1) {  in QuadraticH()
455   PetscCall(PetscLogFlops(199.0 * mx * my));  in QuadraticH()
474   PetscInt   mx = user->mx, my = user->my;  in MSA_BoundaryConditions()  local
484   lsize = my + 2;  in MSA_BoundaryConditions()
485   rsize = my + 2;  in MSA_BoundaryConditions()
494   hy = (t - b) / (my + 1);  in MSA_BoundaryConditions()
572     PetscInt   mx = user->mx, my = user->my;  in MSA_InitialPoint()  local
578     for (j = 0; j < my; j++) {  in MSA_InitialPoint()
581 …x[row] = (((j + 1) * user->bottom[i + 1] + (my - j + 1) * user->top[i + 1]) / (my + 2) + ((i + 1) …  in MSA_InitialPoint()