Lines Matching refs:s_t
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
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()
820 for (j = 0; j < s_x; j++) idx[nn++] = s_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()
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()
828 for (j = 0; j < x_t; j++) idx[nn++] = s_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()
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()
836 for (j = 0; j < s_x; j++) idx[nn++] = s_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()
846 …if (twod && (s_t < 0)) s_t = bases[n3] + (i + 1) * x_t - s_x + x_t * y_t * z_t - x_t * y_t; /* 2D … in DMSetUp_DA_3D()
847 for (j = 0; j < s_x; j++) idx[nn++] = s_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()
855 … if (twod && (s_t < 0)) s_t = bases[n4] + i * x_t + x_t * y_t * z_t - x_t * y_t; /* 2D case */ in DMSetUp_DA_3D()
856 for (j = 0; j < x_t; j++) idx[nn++] = s_t++; 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()
866 … if (twod && (s_t < 0)) s_t = bases[n5] + i * x_t + x_t * y_t * z_t - x_t * y_t; /* 2D case */ in DMSetUp_DA_3D()
867 for (j = 0; j < s_x; j++) idx[nn++] = s_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()
877 …if (twod && (s_t < 0)) s_t = bases[n6] + i * x_t - s_x + x_t * y_t * z_t - x_t * y_t; /* 2D case */ in DMSetUp_DA_3D()
878 for (j = 0; j < s_x; j++) idx[nn++] = s_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()
885 …if (twod && (s_t < 0)) s_t = bases[n7] + (i - 1) * x_t + x_t * y_t * z_t - x_t * y_t; /* 2D case */ in DMSetUp_DA_3D()
886 for (j = 0; j < x_t; j++) idx[nn++] = s_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()
893 …if (twod && (s_t < 0)) s_t = bases[n8] + (i - 1) * x_t + x_t * y_t * z_t - x_t * y_t; /* 2D case */ in DMSetUp_DA_3D()
894 for (j = 0; j < s_x; j++) idx[nn++] = s_t++; 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()
907 for (j = 0; j < s_x; j++) idx[nn++] = s_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()
914 for (j = 0; j < x_t; j++) idx[nn++] = s_t++; 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()
923 for (j = 0; j < s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
932 s_t = bases[n12] + (i + 1) * x_t - s_x + k * x_t * y_t; in DMSetUp_DA_3D()
933 for (j = 0; j < s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
939 s_t = bases[rank] + i * x + k * x * y; in DMSetUp_DA_3D()
940 for (j = 0; j < x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
946 s_t = bases[n14] + i * x_t + k * x_t * y_t; in DMSetUp_DA_3D()
947 for (j = 0; j < s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
958 s_t = bases[n15] + i * x_t - s_x + k * x_t * y_t; in DMSetUp_DA_3D()
959 for (j = 0; j < s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
965 s_t = bases[n16] + (i - 1) * x_t + k * x_t * y_t; in DMSetUp_DA_3D()
966 for (j = 0; j < x_t; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
974 s_t = bases[n17] + (i - 1) * x_t + k * x_t * y_t; in DMSetUp_DA_3D()
975 for (j = 0; j < s_x; j++) idx[nn++] = s_t++; 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()
989 for (j = 0; j < s_x; j++) idx[nn++] = s_t++; 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()
997 for (j = 0; j < x_t; j++) idx[nn++] = s_t++; 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()
1005 for (j = 0; j < s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1014 s_t = bases[n21] + (i + 1) * x_t - s_x + k * x_t * y_t; in DMSetUp_DA_3D()
1015 if (twod && (s_t >= M * N * P)) s_t = bases[n21] + (i + 1) * x_t - s_x; /* 2d case */ in DMSetUp_DA_3D()
1016 for (j = 0; j < s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1023 s_t = bases[n22] + i * x_t + k * x_t * y_t; in DMSetUp_DA_3D()
1024 if (twod && (s_t >= M * N * P)) s_t = bases[n22] + i * x_t; /* 2d case */ in DMSetUp_DA_3D()
1025 for (j = 0; j < x_t; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1034 s_t = bases[n23] + i * x_t + k * x_t * y_t; in DMSetUp_DA_3D()
1035 if (twod && (s_t >= M * N * P)) s_t = bases[n23] + i * x_t; /* 2d case */ in DMSetUp_DA_3D()
1036 for (j = 0; j < s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1045 s_t = bases[n24] + i * x_t - s_x + k * x_t * y_t; in DMSetUp_DA_3D()
1046 if (twod && (s_t >= M * N * P)) s_t = bases[n24] + i * x_t - s_x; /* 2d case */ in DMSetUp_DA_3D()
1047 for (j = 0; j < s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1053 s_t = bases[n25] + (i - 1) * x_t + k * x_t * y_t; in DMSetUp_DA_3D()
1054 if (twod && (s_t >= M * N * P)) s_t = bases[n25] + (i - 1) * x_t; /* 2d case */ in DMSetUp_DA_3D()
1055 for (j = 0; j < x_t; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1061 s_t = bases[n26] + (i - 1) * x_t + k * x_t * y_t; in DMSetUp_DA_3D()
1062 if (twod && (s_t >= M * N * P)) s_t = bases[n26] + (i - 1) * x_t; /* 2d case */ in DMSetUp_DA_3D()
1063 for (j = 0; j < s_x; j++) idx[nn++] = s_t++; 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()
1110 for (j = 0; j < s_x; j++) idx[nn++] = s_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()
1119 for (j = 0; j < x_t; j++) idx[nn++] = s_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()
1128 for (j = 0; j < s_x; j++) idx[nn++] = s_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()
1140 for (j = 0; j < s_x; j++) idx[nn++] = s_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()
1150 for (j = 0; j < x_t; j++) idx[nn++] = s_t++; 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()
1164 for (j = 0; j < s_x; j++) idx[nn++] = s_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()
1176 for (j = 0; j < s_x; j++) idx[nn++] = s_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()
1185 for (j = 0; j < x_t; j++) idx[nn++] = s_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()
1194 for (j = 0; j < s_x; j++) idx[nn++] = s_t++; 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()
1209 for (j = 0; j < s_x; j++) idx[nn++] = s_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()
1218 for (j = 0; j < x_t; j++) idx[nn++] = s_t++; 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()
1231 for (j = 0; j < s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1242 s_t = bases[n12] + (i + 1) * x_t - s_x + k * x_t * y_t; in DMSetUp_DA_3D()
1243 for (j = 0; j < s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1253 s_t = bases[rank] + i * x + k * x * y; in DMSetUp_DA_3D()
1254 for (j = 0; j < x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1260 s_t = bases[n14] + i * x_t + k * x_t * y_t; in DMSetUp_DA_3D()
1261 for (j = 0; j < s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1276 s_t = bases[n15] + i * x_t - s_x + k * x_t * y_t; in DMSetUp_DA_3D()
1277 for (j = 0; j < s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1285 s_t = bases[n16] + (i - 1) * x_t + k * x_t * y_t; in DMSetUp_DA_3D()
1286 for (j = 0; j < x_t; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1298 s_t = bases[n17] + (i - 1) * x_t + k * x_t * y_t; in DMSetUp_DA_3D()
1299 for (j = 0; j < s_x; j++) idx[nn++] = s_t++; 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()
1314 for (j = 0; j < s_x; j++) idx[nn++] = s_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()
1323 for (j = 0; j < x_t; j++) idx[nn++] = s_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()
1332 for (j = 0; j < s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1343 s_t = bases[n21] + (i + 1) * x_t - s_x + k * x_t * y_t; in DMSetUp_DA_3D()
1344 for (j = 0; j < s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1353 s_t = bases[n22] + i * x_t + k * x_t * y_t; in DMSetUp_DA_3D()
1354 for (j = 0; j < x_t; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1367 s_t = bases[n23] + i * x_t + k * x_t * y_t; in DMSetUp_DA_3D()
1368 for (j = 0; j < s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1379 s_t = bases[n24] + i * x_t - s_x + k * x_t * y_t; in DMSetUp_DA_3D()
1380 for (j = 0; j < s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1388 s_t = bases[n25] + (i - 1) * x_t + k * x_t * y_t; in DMSetUp_DA_3D()
1389 for (j = 0; j < x_t; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1397 s_t = bases[n26] + (i - 1) * x_t + k * x_t * y_t; in DMSetUp_DA_3D()
1398 for (j = 0; j < s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()