213   PetscInt       *bases, *ldims, base, x_t, y_t, z_t, s_t, count, s_x, s_y, s_z;  in DMSetUp_DA_3D()  local
456   s_y = s;  in DMSetUp_DA_3D()
813     for (i = 1; i <= s_y; i++) {  in DMSetUp_DA_3D()
818         s_t = bases[n0] + x_t * y_t * z_t - (s_y - i) * x_t - s_x - (s_z - k - 1) * x_t * y_t;  in DMSetUp_DA_3D()
819 …    if (twod && (s_t < 0)) s_t = bases[n0] + x_t * y_t * z_t - (s_y - i) * x_t - s_x; /* 2D case */  in DMSetUp_DA_3D()
826         s_t = bases[n1] + x_t * y_t * z_t - (s_y + 1 - i) * x_t - (s_z - k - 1) * x_t * y_t;  in DMSetUp_DA_3D()
827 …      if (twod && (s_t < 0)) s_t = bases[n1] + x_t * y_t * z_t - (s_y + 1 - i) * x_t; /* 2D case */  in DMSetUp_DA_3D()
834         s_t = bases[n2] + x_t * y_t * z_t - (s_y + 1 - i) * x_t - (s_z - k - 1) * x_t * y_t;  in DMSetUp_DA_3D()
835 …      if (twod && (s_t < 0)) s_t = bases[n2] + x_t * y_t * z_t - (s_y + 1 - i) * x_t; /* 2D case */  in DMSetUp_DA_3D()
871     for (i = 1; i <= s_y; i++) {  in DMSetUp_DA_3D()
901     for (i = 1; i <= s_y; i++) {  in DMSetUp_DA_3D()
906         s_t = bases[n9] - (s_y - i) * x_t - s_x + (k + 1) * x_t * y_t;  in DMSetUp_DA_3D()
913         s_t = bases[n10] - (s_y + 1 - i) * x_t + (k + 1) * x_t * y_t;  in DMSetUp_DA_3D()
916         for (j = 0; j < x; j++) idx[nn++] = bases[rank] + j + k * x * y + (s_y - i) * x;  in DMSetUp_DA_3D()
922         s_t = bases[n11] - (s_y + 1 - i) * x_t + (k + 1) * x_t * y_t;  in DMSetUp_DA_3D()
953     for (i = 1; i <= s_y; i++) {  in DMSetUp_DA_3D()
982     for (i = 1; i <= s_y; i++) {  in DMSetUp_DA_3D()
987         s_t = bases[n18] - (s_y - i) * x_t - s_x + (k + 1) * x_t * y_t;  in DMSetUp_DA_3D()
988 …if (twod && (s_t >= M * N * P)) s_t = bases[n18] - (s_y - i) * x_t - s_x + x_t * y_t; /* 2d case */  in DMSetUp_DA_3D()
995         s_t = bases[n19] - (s_y + 1 - i) * x_t + (k + 1) * x_t * y_t;  in DMSetUp_DA_3D()
996 …  if (twod && (s_t >= M * N * P)) s_t = bases[n19] - (s_y + 1 - i) * x_t + x_t * y_t; /* 2d case */  in DMSetUp_DA_3D()
1003         s_t = bases[n20] - (s_y + 1 - i) * x_t + (k + 1) * x_t * y_t;  in DMSetUp_DA_3D()
1004 …  if (twod && (s_t >= M * N * P)) s_t = bases[n20] - (s_y + 1 - i) * x_t + x_t * y_t; /* 2d case */  in DMSetUp_DA_3D()
1040     for (i = 1; i <= s_y; i++) {  in DMSetUp_DA_3D()
1104       for (i = 1; i <= s_y; i++) {  in DMSetUp_DA_3D()
1109           s_t = bases[n0] + x_t * y_t * z_t - (s_y - i) * x_t - s_x - (s_z - k - 1) * x_t * y_t;  in DMSetUp_DA_3D()
1118           s_t = bases[n1] + x_t * y_t * z_t - (s_y + 1 - i) * x_t - (s_z - k - 1) * x_t * y_t;  in DMSetUp_DA_3D()
1127           s_t = bases[n2] + x_t * y_t * z_t - (s_y + 1 - i) * x_t - (s_z - k - 1) * x_t * y_t;  in DMSetUp_DA_3D()
1170       for (i = 1; i <= s_y; i++) {  in DMSetUp_DA_3D()
1203       for (i = 1; i <= s_y; i++) {  in DMSetUp_DA_3D()
1208           s_t = bases[n9] - (s_y - i) * x_t - s_x + (k + 1) * x_t * y_t;  in DMSetUp_DA_3D()
1217           s_t = bases[n10] - (s_y + 1 - i) * x_t + (k + 1) * x_t * y_t;  in DMSetUp_DA_3D()
1221             for (j = 0; j < x; j++) idx[nn++] = bases[rank] + j + k * x * y + (s_y - i + 1) * x;  in DMSetUp_DA_3D()
1230           s_t = bases[n11] - (s_y + 1 - i) * x_t + (k + 1) * x_t * y_t;  in DMSetUp_DA_3D()
1271       for (i = 1; i <= s_y; i++) {  in DMSetUp_DA_3D()
1308       for (i = 1; i <= s_y; i++) {  in DMSetUp_DA_3D()
1313           s_t = bases[n18] - (s_y - i) * x_t - s_x + (k + 1) * x_t * y_t;  in DMSetUp_DA_3D()
1322           s_t = bases[n19] - (s_y + 1 - i) * x_t + (k + 1) * x_t * y_t;  in DMSetUp_DA_3D()
1331           s_t = bases[n20] - (s_y + 1 - i) * x_t + (k + 1) * x_t * y_t;  in DMSetUp_DA_3D()
1374       for (i = 1; i <= s_y; i++) {  in DMSetUp_DA_3D()