25   PetscInt   mx, my;  member
54   user.mx = 4;  in main()
58   PetscCall(PetscOptionsGetInt(NULL, NULL, "-mx", &user.mx, &flg));  in main()
62   N = user.mx * user.my;  in main()
65 …PetscCall(PetscPrintf(PETSC_COMM_SELF, "mx:%" PetscInt_FMT ", my:%" PetscInt_FMT "\n", user.mx, us…  in main()
140   PetscInt     mx = user->mx, my = user->my;  in FormConstraints()  local
141   PetscReal    hx = 1.0 / (mx + 1), hy = 1.0 / (my + 1), hydhx = hy / hx, hxdhy = hx / hy;  in FormConstraints()
157     for (i = 0; i < mx; i++) {  in FormConstraints()
158       row = j * mx + i;  in FormConstraints()
174         xb = x[row - mx];  in FormConstraints()
177       if (i + 1 == mx) { /* right side */  in FormConstraints()
188         xt = x[row + mx];  in FormConstraints()
191       if (i > 0 && j + 1 < my) xlt = x[row - 1 + mx];  in FormConstraints()
192       if (j > 0 && i + 1 < mx) xrb = x[row + 1 - mx];  in FormConstraints()
240   PetscCall(PetscLogFlops(67 * mx * my));  in FormConstraints()
261   PetscInt           mx = user->mx, my = user->my;  in FormJacobian()  local
263   PetscReal          hx = 1.0 / (mx + 1), hy = 1.0 / (my + 1), hydhx = hy / hx, hxdhy = hx / hy;  in FormJacobian()
280   for (i = 0; i < mx; i++) {  in FormJacobian()
282       row = j * mx + i;  in FormJacobian()
299         xb = x[row - mx];  in FormJacobian()
302       if (i + 1 == mx) {  in FormJacobian()
313         xt = x[row + mx];  in FormJacobian()
316       if (i > 0 && j + 1 < my) xlt = x[row - 1 + mx];  in FormJacobian()
317       if (j > 0 && i + 1 < mx) xrb = x[row + 1 - mx];  in FormJacobian()
356         col[k] = row - mx;  in FormJacobian()
360       if (j > 0 && i < mx - 1) {  in FormJacobian()
362         col[k] = row - mx + 1;  in FormJacobian()
376       if (i < mx - 1) {  in FormJacobian()
384         col[k] = row + mx - 1;  in FormJacobian()
390         col[k] = row + mx;  in FormJacobian()
408   PetscCall(PetscLogFlops(199 * mx * my));  in FormJacobian()
426   PetscInt   mx = user->mx, my = user->my;  in MSA_BoundaryConditions()  local
435   bsize = mx + 2;  in MSA_BoundaryConditions()
438   tsize = mx + 2;  in MSA_BoundaryConditions()
445   hx = (r - l) / (mx + 1);  in MSA_BoundaryConditions()
523     PetscInt     mx = user->mx, my = user->my;  in MSA_InitialPoint()  local
531       for (i = 0; i < mx; i++) {  in MSA_InitialPoint()
532         row    = (j)*mx + (i);  in MSA_InitialPoint()
533 …top[i + 1]) / (my + 2) + ((i + 1) * user->left[j + 1] + (mx - i + 1) * user->right[j + 1]) / (mx +…  in MSA_InitialPoint()