Lines Matching refs:s_t
221 PetscInt *bases,*ldims,base,x_t,y_t,z_t,s_t,count,s_x,s_y,s_z; in DMSetUp_DA_3D() local
778 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()
779 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()
780 for (j=0; j<s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
786 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()
787 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()
788 for (j=0; j<x_t; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
794 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()
795 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()
796 for (j=0; j<s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
805 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()
806 … 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 case */ in DMSetUp_DA_3D()
807 for (j=0; j<s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
814 s_t = bases[n4] + i*x_t + x_t*y_t*z_t - (s_z-k)*x_t*y_t; in DMSetUp_DA_3D()
815 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()
816 for (j=0; j<x_t; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
825 s_t = bases[n5] + i*x_t + x_t*y_t*z_t - (s_z-k)*x_t*y_t; in DMSetUp_DA_3D()
826 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()
827 for (j=0; j<s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
836 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()
837 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()
838 for (j=0; j<s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
844 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()
845 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()
846 for (j=0; j<x_t; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
852 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()
853 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()
854 for (j=0; j<s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
866 s_t = bases[n9] - (s_y-i)*x_t -s_x + (k+1)*x_t*y_t; in DMSetUp_DA_3D()
867 for (j=0; j<s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
873 s_t = bases[n10] - (s_y+1-i)*x_t + (k+1)*x_t*y_t; in DMSetUp_DA_3D()
874 for (j=0; j<x_t; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
882 s_t = bases[n11] - (s_y+1-i)*x_t + (k+1)*x_t*y_t; in DMSetUp_DA_3D()
883 for (j=0; j<s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
892 s_t = bases[n12] + (i+1)*x_t - s_x + k*x_t*y_t; in DMSetUp_DA_3D()
893 for (j=0; j<s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
899 s_t = bases[rank] + i*x + k*x*y; in DMSetUp_DA_3D()
900 for (j=0; j<x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
906 s_t = bases[n14] + i*x_t + k*x_t*y_t; in DMSetUp_DA_3D()
907 for (j=0; j<s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
918 s_t = bases[n15] + i*x_t - s_x + k*x_t*y_t; in DMSetUp_DA_3D()
919 for (j=0; j<s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
925 s_t = bases[n16] + (i-1)*x_t + k*x_t*y_t; in DMSetUp_DA_3D()
926 for (j=0; j<x_t; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
934 s_t = bases[n17] + (i-1)*x_t + k*x_t*y_t; in DMSetUp_DA_3D()
935 for (j=0; j<s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
947 s_t = bases[n18] - (s_y-i)*x_t -s_x + (k+1)*x_t*y_t; in DMSetUp_DA_3D()
948 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()
949 for (j=0; j<s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
955 s_t = bases[n19] - (s_y+1-i)*x_t + (k+1)*x_t*y_t; in DMSetUp_DA_3D()
956 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()
957 for (j=0; j<x_t; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
963 s_t = bases[n20] - (s_y+1-i)*x_t + (k+1)*x_t*y_t; in DMSetUp_DA_3D()
964 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()
965 for (j=0; j<s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
974 s_t = bases[n21] + (i+1)*x_t - s_x + k*x_t*y_t; in DMSetUp_DA_3D()
975 if (twod && (s_t >= M*N*P)) s_t = bases[n21] + (i+1)*x_t - s_x; /* 2d case */ in DMSetUp_DA_3D()
976 for (j=0; j<s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
983 s_t = bases[n22] + i*x_t + k*x_t*y_t; in DMSetUp_DA_3D()
984 if (twod && (s_t >= M*N*P)) s_t = bases[n22] + i*x_t; /* 2d case */ in DMSetUp_DA_3D()
985 for (j=0; j<x_t; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
994 s_t = bases[n23] + i*x_t + k*x_t*y_t; in DMSetUp_DA_3D()
995 if (twod && (s_t >= M*N*P)) s_t = bases[n23] + i*x_t; /* 2d case */ in DMSetUp_DA_3D()
996 for (j=0; j<s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1005 s_t = bases[n24] + i*x_t - s_x + k*x_t*y_t; in DMSetUp_DA_3D()
1006 if (twod && (s_t >= M*N*P)) s_t = bases[n24] + i*x_t - s_x; /* 2d case */ in DMSetUp_DA_3D()
1007 for (j=0; j<s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1013 s_t = bases[n25] + (i-1)*x_t + k*x_t*y_t; in DMSetUp_DA_3D()
1014 if (twod && (s_t >= M*N*P)) s_t = bases[n25] + (i-1)*x_t; /* 2d case */ in DMSetUp_DA_3D()
1015 for (j=0; j<x_t; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1021 s_t = bases[n26] + (i-1)*x_t + k*x_t*y_t; in DMSetUp_DA_3D()
1022 if (twod && (s_t >= M*N*P)) s_t = bases[n26] + (i-1)*x_t; /* 2d case */ in DMSetUp_DA_3D()
1023 for (j=0; j<s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1054 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()
1055 for (j=0; j<s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1063 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()
1064 for (j=0; j<x_t; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1072 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()
1073 for (j=0; j<s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1084 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()
1085 for (j=0; j<s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1094 s_t = bases[n4] + i*x_t + x_t*y_t*z_t - (s_z-k)*x_t*y_t; in DMSetUp_DA_3D()
1095 for (j=0; j<x_t; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1108 s_t = bases[n5] + i*x_t + x_t*y_t*z_t - (s_z-k)*x_t*y_t; in DMSetUp_DA_3D()
1109 for (j=0; j<s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1120 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()
1121 for (j=0; j<s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1129 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()
1130 for (j=0; j<x_t; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1138 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()
1139 for (j=0; j<s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1153 s_t = bases[n9] - (s_y-i)*x_t -s_x + (k+1)*x_t*y_t; in DMSetUp_DA_3D()
1154 for (j=0; j<s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1162 s_t = bases[n10] - (s_y+1-i)*x_t + (k+1)*x_t*y_t; in DMSetUp_DA_3D()
1163 for (j=0; j<x_t; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1175 s_t = bases[n11] - (s_y+1-i)*x_t + (k+1)*x_t*y_t; in DMSetUp_DA_3D()
1176 for (j=0; j<s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1187 s_t = bases[n12] + (i+1)*x_t - s_x + k*x_t*y_t; in DMSetUp_DA_3D()
1188 for (j=0; j<s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1198 s_t = bases[rank] + i*x + k*x*y; in DMSetUp_DA_3D()
1199 for (j=0; j<x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1205 s_t = bases[n14] + i*x_t + k*x_t*y_t; in DMSetUp_DA_3D()
1206 for (j=0; j<s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1221 s_t = bases[n15] + i*x_t - s_x + k*x_t*y_t; in DMSetUp_DA_3D()
1222 for (j=0; j<s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1230 s_t = bases[n16] + (i-1)*x_t + k*x_t*y_t; in DMSetUp_DA_3D()
1231 for (j=0; j<x_t; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1243 s_t = bases[n17] + (i-1)*x_t + k*x_t*y_t; in DMSetUp_DA_3D()
1244 for (j=0; j<s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1258 s_t = bases[n18] - (s_y-i)*x_t -s_x + (k+1)*x_t*y_t; in DMSetUp_DA_3D()
1259 for (j=0; j<s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1267 s_t = bases[n19] - (s_y+1-i)*x_t + (k+1)*x_t*y_t; in DMSetUp_DA_3D()
1268 for (j=0; j<x_t; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1276 s_t = bases[n20] - (s_y+1-i)*x_t + (k+1)*x_t*y_t; in DMSetUp_DA_3D()
1277 for (j=0; j<s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1288 s_t = bases[n21] + (i+1)*x_t - s_x + k*x_t*y_t; in DMSetUp_DA_3D()
1289 for (j=0; j<s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1298 s_t = bases[n22] + i*x_t + k*x_t*y_t; in DMSetUp_DA_3D()
1299 for (j=0; j<x_t; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1312 s_t = bases[n23] + i*x_t + k*x_t*y_t; in DMSetUp_DA_3D()
1313 for (j=0; j<s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1324 s_t = bases[n24] + i*x_t - s_x + k*x_t*y_t; in DMSetUp_DA_3D()
1325 for (j=0; j<s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1333 s_t = bases[n25] + (i-1)*x_t + k*x_t*y_t; in DMSetUp_DA_3D()
1334 for (j=0; j<x_t; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()
1342 s_t = bases[n26] + (i-1)*x_t + k*x_t*y_t; in DMSetUp_DA_3D()
1343 for (j=0; j<s_x; j++) idx[nn++] = s_t++; in DMSetUp_DA_3D()