8   PetscInt mx;    /* grid points in each direction */  member
109   user.mx = 8;  in main()
111 …PetscCall(PetscOptionsInt("-mx", "Number of grid points in each direction", "", user.mx, &user.mx,…  in main()
130   user.m       = user.ns * user.mx * user.mx * user.mx; /* number of constraints */  in main()
132   user.ndesign = user.mx * user.mx * user.mx;  in main()
655   PetscCall(PetscMalloc1(user->mx, &x));  in EllipticInitialize()
656   PetscCall(PetscMalloc1(user->mx, &y));  in EllipticInitialize()
657   PetscCall(PetscMalloc1(user->mx, &z));  in EllipticInitialize()
662   n         = user->mx * user->mx * user->mx;  in EllipticInitialize()
663   m         = 3 * user->mx * user->mx * (user->mx - 1);  in EllipticInitialize()
688   h        = 1.0 / user->mx;  in EllipticInitialize()
689   hinv     = user->mx;  in EllipticInitialize()
694     i  = linear_index % user->mx;  in EllipticInitialize()
695     j  = ((linear_index - i) / user->mx) % user->mx;  in EllipticInitialize()
696     k  = ((linear_index - i) / user->mx - j) / user->mx;  in EllipticInitialize()
802       iblock = i / (user->mx - 1);  in EllipticInitialize()
803       j      = iblock * user->mx + (i % (user->mx - 1));  in EllipticInitialize()
809       iblock = (i - m / 3) / (user->mx * (user->mx - 1));  in EllipticInitialize()
810       j      = iblock * user->mx * user->mx + ((i - m / 3) % (user->mx * (user->mx - 1)));  in EllipticInitialize()
812       j = j + user->mx;  in EllipticInitialize()
818       j = j + user->mx * user->mx;  in EllipticInitialize()
836       iblock = i / (user->mx - 1);  in EllipticInitialize()
837       j      = iblock * user->mx + (i % (user->mx - 1));  in EllipticInitialize()
843       iblock = (i - m / 3) / (user->mx * (user->mx - 1));  in EllipticInitialize()
844       j      = iblock * user->mx * user->mx + ((i - m / 3) % (user->mx * (user->mx - 1)));  in EllipticInitialize()
846       j = j + user->mx;  in EllipticInitialize()
852       j = j + user->mx * user->mx;  in EllipticInitialize()
869       iblock = i / (user->mx - 1);  in EllipticInitialize()
870       j      = iblock * user->mx + (i % (user->mx - 1));  in EllipticInitialize()
876       iblock = (i - m / 3) / (user->mx * (user->mx - 1));  in EllipticInitialize()
877       j      = iblock * user->mx * user->mx + ((i - m / 3) % (user->mx * (user->mx - 1)));  in EllipticInitialize()
879       j = j + user->mx;  in EllipticInitialize()
885       j = j + user->mx * user->mx;  in EllipticInitialize()
909         iblock = i / (user->mx - 1);  in EllipticInitialize()
910         j      = iblock * user->mx + (i % (user->mx - 1));  in EllipticInitialize()
916         iblock = (i - m / 3) / (user->mx * (user->mx - 1));  in EllipticInitialize()
917         j      = iblock * user->mx * user->mx + ((i - m / 3) % (user->mx * (user->mx - 1)));  in EllipticInitialize()
919         j = j + user->mx;  in EllipticInitialize()
925         j = j + user->mx * user->mx;  in EllipticInitialize()
947   PetscCall(VecSetSizes(user->lwork, PETSC_DECIDE, m + user->mx * user->mx * user->mx));  in EllipticInitialize()
1061   for (i = 0; i < user->mx; i++) {  in EllipticInitialize()
1068   nx = user->mx;  in EllipticInitialize()
1069   ny = user->mx;  in EllipticInitialize()
1070   nz = user->mx;  in EllipticInitialize()