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
457   s_z = s;  in DMSetUp_DA_3D()
812   for (k = 0; k < s_z; k++) {  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()
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()
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()
845         s_t = bases[n3] + (i + 1) * x_t - s_x + x_t * y_t * z_t - (s_z - k) * x_t * y_t;  in DMSetUp_DA_3D()
854         s_t = bases[n4] + i * x_t + x_t * y_t * z_t - (s_z - k) * x_t * y_t;  in DMSetUp_DA_3D()
858         for (j = 0; j < x; j++) idx[nn++] = bases[rank] + j + x * i + (s_z - k - 1) * x * y;  in DMSetUp_DA_3D()
865         s_t = bases[n5] + i * x_t + x_t * y_t * z_t - (s_z - k) * x_t * y_t;  in DMSetUp_DA_3D()
876         s_t = bases[n6] + i * x_t - s_x + x_t * y_t * z_t - (s_z - k) * x_t * y_t;  in DMSetUp_DA_3D()
884         s_t = bases[n7] + (i - 1) * x_t + x_t * y_t * z_t - (s_z - k) * x_t * y_t;  in DMSetUp_DA_3D()
892         s_t = bases[n8] + (i - 1) * x_t + x_t * y_t * z_t - (s_z - k) * x_t * y_t;  in DMSetUp_DA_3D()
981   for (k = 0; k < s_z; k++) {  in DMSetUp_DA_3D()
1103     for (k = 0; k < s_z; k++) {  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()
1139           s_t = bases[n3] + (i + 1) * x_t - s_x + x_t * y_t * z_t - (s_z - k) * x_t * y_t;  in DMSetUp_DA_3D()
1149           s_t = bases[n4] + i * x_t + x_t * y_t * z_t - (s_z - k) * x_t * y_t;  in DMSetUp_DA_3D()
1153             for (j = 0; j < x; j++) idx[nn++] = bases[rank] + j + x * i + (s_z - k) * x * y;  in DMSetUp_DA_3D()
1163           s_t = bases[n5] + i * x_t + x_t * y_t * z_t - (s_z - k) * x_t * y_t;  in DMSetUp_DA_3D()
1175           s_t = bases[n6] + i * x_t - s_x + x_t * y_t * z_t - (s_z - k) * x_t * y_t;  in DMSetUp_DA_3D()
1184           s_t = bases[n7] + (i - 1) * x_t + x_t * y_t * z_t - (s_z - k) * x_t * y_t;  in DMSetUp_DA_3D()
1193           s_t = bases[n8] + (i - 1) * x_t + x_t * y_t * z_t - (s_z - k) * x_t * y_t;  in DMSetUp_DA_3D()
1307     for (k = 0; k < s_z; k++) {  in DMSetUp_DA_3D()