819 const V3d &vsP)
const
834 V3i c1(
static_cast<int>(floor(p.x)),
835 static_cast<int>(floor(p.y)),
836 static_cast<int>(floor(p.z)));
847 c1.x = std::min(dataWindow.max.x + 1, std::max(dataWindow.min.x, c1.x));
848 c1.y = std::min(dataWindow.max.y, std::max(dataWindow.min.y, c1.y));
849 c1.z = std::min(dataWindow.max.z, std::max(dataWindow.min.z, c1.z));
850 c2.x = std::min(dataWindow.max.x + 1, std::max(dataWindow.min.x, c2.x));
851 c2.y = std::min(dataWindow.max.y, std::max(dataWindow.min.y, c2.y));
852 c2.z = std::min(dataWindow.max.z, std::max(dataWindow.min.z, c2.z));
854 ret.x = (f1.x * (f1.y * (f1.z * data.
u(c1.x, c1.y, c1.z) +
855 f2.z * data.
u(c1.x, c1.y, c2.z)) +
856 f2.y * (f1.z * data.
u(c1.x, c2.y, c1.z) +
857 f2.z * data.
u(c1.x, c2.y, c2.z))) +
858 f2.x * (f1.y * (f1.z * data.
u(c2.x, c1.y, c1.z) +
859 f2.z * data.
u(c2.x, c1.y, c2.z)) +
860 f2.y * (f1.z * data.
u(c2.x, c2.y, c1.z) +
861 f2.z * data.
u(c2.x, c2.y, c2.z))));
865 p.setValue(vsP.x - 0.5, vsP.y , vsP.z - 0.5);
868 c1.x =
static_cast<int>(floor(p.x ));
869 c1.y =
static_cast<int>(floor(p.y ));
870 c1.z =
static_cast<int>(floor(p.z ));
883 c1.x = std::min(dataWindow.max.x, std::max(dataWindow.min.x, c1.x));
884 c1.y = std::min(dataWindow.max.y + 1, std::max(dataWindow.min.y, c1.y));
885 c1.z = std::min(dataWindow.max.z, std::max(dataWindow.min.z, c1.z));
886 c2.x = std::min(dataWindow.max.x, std::max(dataWindow.min.x, c2.x));
887 c2.y = std::min(dataWindow.max.y + 1, std::max(dataWindow.min.y, c2.y));
888 c2.z = std::min(dataWindow.max.z, std::max(dataWindow.min.z, c2.z));
890 ret.y = (f1.x * (f1.y * (f1.z * data.
v(c1.x, c1.y, c1.z) +
891 f2.z * data.
v(c1.x, c1.y, c2.z)) +
892 f2.y * (f1.z * data.
v(c1.x, c2.y, c1.z) +
893 f2.z * data.
v(c1.x, c2.y, c2.z))) +
894 f2.x * (f1.y * (f1.z * data.
v(c2.x, c1.y, c1.z) +
895 f2.z * data.
v(c2.x, c1.y, c2.z)) +
896 f2.y * (f1.z * data.
v(c2.x, c2.y, c1.z) +
897 f2.z * data.
v(c2.x, c2.y, c2.z))));
901 p.setValue(vsP.x - 0.5 , vsP.y - 0.5, vsP.z);
904 c1.x =
static_cast<int>(floor(p.x ));
905 c1.y =
static_cast<int>(floor(p.y ));
906 c1.z =
static_cast<int>(floor(p.z ));
919 c1.x = std::min(dataWindow.max.x, std::max(dataWindow.min.x, c1.x));
920 c1.y = std::min(dataWindow.max.y, std::max(dataWindow.min.y, c1.y));
921 c1.z = std::min(dataWindow.max.z + 1, std::max(dataWindow.min.z, c1.z));
922 c2.x = std::min(dataWindow.max.x, std::max(dataWindow.min.x, c2.x));
923 c2.y = std::min(dataWindow.max.y, std::max(dataWindow.min.y, c2.y));
924 c2.z = std::min(dataWindow.max.z + 1, std::max(dataWindow.min.z, c2.z));
926 ret.z = (f1.x * (f1.y * (f1.z * data.
w(c1.x, c1.y, c1.z) +
927 f2.z * data.
w(c1.x, c1.y, c2.z)) +
928 f2.y * (f1.z * data.
w(c1.x, c2.y, c1.z) +
929 f2.z * data.
w(c1.x, c2.y, c2.z))) +
930 f2.x * (f1.y * (f1.z * data.
w(c2.x, c1.y, c1.z) +
931 f2.z * data.
w(c2.x, c1.y, c2.z)) +
932 f2.y * (f1.z * data.
w(c2.x, c2.y, c1.z) +
933 f2.z * data.
w(c2.x, c2.y, c2.z))));
943 const V3d &vsP)
const
961 p.setValue<>(vsP.x, vsP.y-0.5, vsP.z-0.5);
964 c1.x =
static_cast<int>(floor(p.x));
965 c1.y =
static_cast<int>(floor(p.y));
966 c1.z =
static_cast<int>(floor(p.z));
979 c1.x = std::min(dataWindow.max.x + 1, std::max(dataWindow.min.x, c1.x));
980 c1.y = std::min(dataWindow.max.y, std::max(dataWindow.min.y, c1.y));
981 c1.z = std::min(dataWindow.max.z, std::max(dataWindow.min.z, c1.z));
982 c2.x = std::min(dataWindow.max.x + 1, std::max(dataWindow.min.x, c2.x));
983 c2.y = std::min(dataWindow.max.y, std::max(dataWindow.min.y, c2.y));
984 c2.z = std::min(dataWindow.max.z, std::max(dataWindow.min.z, c2.z));
986 ret = (f1.x * (f1.y * (f1.z * data.
u(c1.x, c1.y, c1.z) +
987 f2.z * data.
u(c1.x, c1.y, c2.z)) +
988 f2.y * (f1.z * data.
u(c1.x, c2.y, c1.z) +
989 f2.z * data.
u(c1.x, c2.y, c2.z))) +
990 f2.x * (f1.y * (f1.z * data.
u(c2.x, c1.y, c1.z) +
991 f2.z * data.
u(c2.x, c1.y, c2.z)) +
992 f2.y * (f1.z * data.
u(c2.x, c2.y, c1.z) +
993 f2.z * data.
u(c2.x, c2.y, c2.z))));
999 p.setValue(vsP.x-0.5, vsP.y, vsP.z-0.5);
1002 c1.x =
static_cast<int>(floor(p.x ));
1003 c1.y =
static_cast<int>(floor(p.y ));
1004 c1.z =
static_cast<int>(floor(p.z ));
1017 c1.x = std::min(dataWindow.max.x, std::max(dataWindow.min.x, c1.x));
1018 c1.y = std::min(dataWindow.max.y + 1, std::max(dataWindow.min.y, c1.y));
1019 c1.z = std::min(dataWindow.max.z, std::max(dataWindow.min.z, c1.z));
1020 c2.x = std::min(dataWindow.max.x, std::max(dataWindow.min.x, c2.x));
1021 c2.y = std::min(dataWindow.max.y + 1, std::max(dataWindow.min.y, c2.y));
1022 c2.z = std::min(dataWindow.max.z, std::max(dataWindow.min.z, c2.z));
1024 ret = (f1.x * (f1.y * (f1.z * data.
v(c1.x, c1.y, c1.z) +
1025 f2.z * data.
v(c1.x, c1.y, c2.z)) +
1026 f2.y * (f1.z * data.
v(c1.x, c2.y, c1.z) +
1027 f2.z * data.
v(c1.x, c2.y, c2.z))) +
1028 f2.x * (f1.y * (f1.z * data.
v(c2.x, c1.y, c1.z) +
1029 f2.z * data.
v(c2.x, c1.y, c2.z)) +
1030 f2.y * (f1.z * data.
v(c2.x, c2.y, c1.z) +
1031 f2.z * data.
v(c2.x, c2.y, c2.z))));
1037 p.setValue(vsP.x-0.5, vsP.y-0.5, vsP.z);
1040 c1.x =
static_cast<int>(floor(p.x ));
1041 c1.y =
static_cast<int>(floor(p.y ));
1042 c1.z =
static_cast<int>(floor(p.z ));
1055 c1.x = std::min(dataWindow.max.x, std::max(dataWindow.min.x, c1.x));
1056 c1.y = std::min(dataWindow.max.y, std::max(dataWindow.min.y, c1.y));
1057 c1.z = std::min(dataWindow.max.z + 1, std::max(dataWindow.min.z, c1.z));
1058 c2.x = std::min(dataWindow.max.x, std::max(dataWindow.min.x, c2.x));
1059 c2.y = std::min(dataWindow.max.y, std::max(dataWindow.min.y, c2.y));
1060 c2.z = std::min(dataWindow.max.z + 1, std::max(dataWindow.min.z, c2.z));
1062 ret = (f1.x * (f1.y * (f1.z * data.
w(c1.x, c1.y, c1.z) +
1063 f2.z * data.
w(c1.x, c1.y, c2.z)) +
1064 f2.y * (f1.z * data.
w(c1.x, c2.y, c1.z) +
1065 f2.z * data.
w(c1.x, c2.y, c2.z))) +
1066 f2.x * (f1.y * (f1.z * data.
w(c2.x, c1.y, c1.z) +
1067 f2.z * data.
w(c2.x, c1.y, c2.z)) +
1068 f2.y * (f1.z * data.
w(c2.x, c2.y, c1.z) +
1069 f2.z * data.
w(c2.x, c2.y, c2.z))));
1084 const V3d &vsP)
const
1086 typedef typename Field_T::value_type Data_T;
1091 V3d clampedVsP(std::max(0.5, vsP.x),
1092 std::max(0.5, vsP.y),
1093 std::max(0.5, vsP.z));
1096 const Box3i &dataWindow = data.dataWindow();
1099 V3i c(
static_cast<int>(floor(p.x)),
1100 static_cast<int>(floor(p.y)),
1101 static_cast<int>(floor(p.z)));
1108 im = std::max(dataWindow.min.x, std::min(c.x, dataWindow.max.x));
1109 jm = std::max(dataWindow.min.y, std::min(c.y, dataWindow.max.y));
1110 km = std::max(dataWindow.min.z, std::min(c.z, dataWindow.max.z));
1111 int im_1, jm_1, km_1;
1112 im_1 = std::max(dataWindow.min.x, std::min(im - 1, dataWindow.max.x));
1113 jm_1 = std::max(dataWindow.min.y, std::min(jm - 1, dataWindow.max.y));
1114 km_1 = std::max(dataWindow.min.z, std::min(km - 1, dataWindow.max.z));
1116 im1 = std::max(dataWindow.min.x, std::min(im + 1, dataWindow.max.x));
1117 jm1 = std::max(dataWindow.min.y, std::min(jm + 1, dataWindow.max.y));
1118 km1 = std::max(dataWindow.min.z, std::min(km + 1, dataWindow.max.z));
1120 im2 = std::max(dataWindow.min.x, std::min(im + 2, dataWindow.max.x));
1121 jm2 = std::max(dataWindow.min.y, std::min(jm + 2, dataWindow.max.y));
1122 km2 = std::max(dataWindow.min.z, std::min(km + 2, dataWindow.max.z));
1125 data.fastValue(im_1, jm_1, km),
1126 data.fastValue(im_1, jm_1, km1),
1127 data.fastValue(im_1, jm_1, km2), t.z);
1129 data.fastValue(im_1, jm, km),
1130 data.fastValue(im_1, jm, km1),
1131 data.fastValue(im_1, jm, km2), t.z);
1133 data.fastValue(im_1, jm1, km),
1134 data.fastValue(im_1, jm1, km1),
1135 data.fastValue(im_1, jm1, km2), t.z);
1137 data.fastValue(im_1, jm2, km),
1138 data.fastValue(im_1, jm2, km1),
1139 data.fastValue(im_1, jm2, km2), t.z);
1142 data.fastValue(im, jm_1, km),
1143 data.fastValue(im, jm_1, km1),
1144 data.fastValue(im, jm_1, km2), t.z);
1146 data.fastValue(im, jm, km),
1147 data.fastValue(im, jm, km1),
1148 data.fastValue(im, jm, km2), t.z);
1150 data.fastValue(im, jm1, km),
1151 data.fastValue(im, jm1, km1),
1152 data.fastValue(im, jm1, km2), t.z);
1154 data.fastValue(im, jm2, km),
1155 data.fastValue(im, jm2, km1),
1156 data.fastValue(im, jm2, km2), t.z);
1159 data.fastValue(im1, jm_1, km),
1160 data.fastValue(im1, jm_1, km1),
1161 data.fastValue(im1, jm_1, km2), t.z);
1163 data.fastValue(im1, jm, km),
1164 data.fastValue(im1, jm, km1),
1165 data.fastValue(im1, jm, km2), t.z);
1167 data.fastValue(im1, jm1, km),
1168 data.fastValue(im1, jm1, km1),
1169 data.fastValue(im1, jm1, km2), t.z);
1171 data.fastValue(im1, jm2, km),
1172 data.fastValue(im1, jm2, km1),
1173 data.fastValue(im1, jm2, km2), t.z);
1176 data.fastValue(im2, jm_1, km),
1177 data.fastValue(im2, jm_1, km1),
1178 data.fastValue(im2, jm_1, km2), t.z);
1180 data.fastValue(im2, jm, km),
1181 data.fastValue(im2, jm, km1),
1182 data.fastValue(im2, jm, km2), t.z);
1184 data.fastValue(im2, jm1, km),
1185 data.fastValue(im2, jm1, km1),
1186 data.fastValue(im2, jm1, km2), t.z);
1188 data.fastValue(im2, jm2, km),
1189 data.fastValue(im2, jm2, km1),
1190 data.fastValue(im2, jm2, km2), t.z);
1206 const V3d &vsP)
const
1208 typedef typename Data_T::BaseType T;
1220 V3d clampedVsP(std::max(0.5, vsP.x),
1221 std::max(0.5, vsP.y),
1222 std::max(0.5, vsP.z));
1225 clampedVsP.z - 0.5);
1228 V3i c(
static_cast<int>(floor(p.x)),
1229 static_cast<int>(floor(p.y)),
1230 static_cast<int>(floor(p.z)));
1237 im = std::max(dataWindow.min.x, std::min(c.x, dataWindow.max.x + 1));
1238 jm = std::max(dataWindow.min.y, std::min(c.y, dataWindow.max.y));
1239 km = std::max(dataWindow.min.z, std::min(c.z, dataWindow.max.z));
1240 int im_1, jm_1, km_1;
1241 im_1 = std::max(dataWindow.min.x, std::min(im - 1, dataWindow.max.x + 1));
1242 jm_1 = std::max(dataWindow.min.y, std::min(jm - 1, dataWindow.max.y));
1243 km_1 = std::max(dataWindow.min.z, std::min(km - 1, dataWindow.max.z));
1245 im1 = std::max(dataWindow.min.x, std::min(im + 1, dataWindow.max.x + 1));
1246 jm1 = std::max(dataWindow.min.y, std::min(jm + 1, dataWindow.max.y));
1247 km1 = std::max(dataWindow.min.z, std::min(km + 1, dataWindow.max.z));
1249 im2 = std::max(dataWindow.min.x, std::min(im + 2, dataWindow.max.x + 1));
1250 jm2 = std::max(dataWindow.min.y, std::min(jm + 2, dataWindow.max.y));
1251 km2 = std::max(dataWindow.min.z, std::min(km + 2, dataWindow.max.z));
1254 data.
u(im_1, jm_1, km),
1255 data.
u(im_1, jm_1, km1),
1256 data.
u(im_1, jm_1, km2), t.z);
1258 data.
u(im_1, jm, km),
1259 data.
u(im_1, jm, km1),
1260 data.
u(im_1, jm, km2), t.z);
1262 data.
u(im_1, jm1, km),
1263 data.
u(im_1, jm1, km1),
1264 data.
u(im_1, jm1, km2), t.z);
1266 data.
u(im_1, jm2, km),
1267 data.
u(im_1, jm2, km1),
1268 data.
u(im_1, jm2, km2), t.z);
1271 data.
u(im, jm_1, km),
1272 data.
u(im, jm_1, km1),
1273 data.
u(im, jm_1, km2), t.z);
1276 data.
u(im, jm, km1),
1277 data.
u(im, jm, km2), t.z);
1279 data.
u(im, jm1, km),
1280 data.
u(im, jm1, km1),
1281 data.
u(im, jm1, km2), t.z);
1283 data.
u(im, jm2, km),
1284 data.
u(im, jm2, km1),
1285 data.
u(im, jm2, km2), t.z);
1288 data.
u(im1, jm_1, km),
1289 data.
u(im1, jm_1, km1),
1290 data.
u(im1, jm_1, km2), t.z);
1292 data.
u(im1, jm, km),
1293 data.
u(im1, jm, km1),
1294 data.
u(im1, jm, km2), t.z);
1296 data.
u(im1, jm1, km),
1297 data.
u(im1, jm1, km1),
1298 data.
u(im1, jm1, km2), t.z);
1300 data.
u(im1, jm2, km),
1301 data.
u(im1, jm2, km1),
1302 data.
u(im1, jm2, km2), t.z);
1305 data.
u(im2, jm_1, km),
1306 data.
u(im2, jm_1, km1),
1307 data.
u(im2, jm_1, km2), t.z);
1309 data.
u(im2, jm, km),
1310 data.
u(im2, jm, km1),
1311 data.
u(im2, jm, km2), t.z);
1313 data.
u(im2, jm1, km),
1314 data.
u(im2, jm1, km1),
1315 data.
u(im2, jm1, km2), t.z);
1317 data.
u(im2, jm2, km),
1318 data.
u(im2, jm2, km1),
1319 data.
u(im2, jm2, km2), t.z);
1332 p.setValue(clampedVsP.x - 0.5, vsP.y , clampedVsP.z - 0.5);
1335 c.x =
static_cast<int>(floor(p.x));
1336 c.y =
static_cast<int>(floor(p.y));
1337 c.z =
static_cast<int>(floor(p.z));
1343 im = std::max(dataWindow.min.x, std::min(c.x, dataWindow.max.x));
1344 jm = std::max(dataWindow.min.y, std::min(c.y, dataWindow.max.y + 1));
1345 km = std::max(dataWindow.min.z, std::min(c.z, dataWindow.max.z));
1346 int im_1, jm_1, km_1;
1347 im_1 = std::max(dataWindow.min.x, std::min(im - 1, dataWindow.max.x));
1348 jm_1 = std::max(dataWindow.min.y, std::min(jm - 1, dataWindow.max.y + 1));
1349 km_1 = std::max(dataWindow.min.z, std::min(km - 1, dataWindow.max.z));
1351 im1 = std::max(dataWindow.min.x, std::min(im + 1, dataWindow.max.x));
1352 jm1 = std::max(dataWindow.min.y, std::min(jm + 1, dataWindow.max.y + 1));
1353 km1 = std::max(dataWindow.min.z, std::min(km + 1, dataWindow.max.z));
1355 im2 = std::max(dataWindow.min.x, std::min(im + 2, dataWindow.max.x));
1356 jm2 = std::max(dataWindow.min.y, std::min(jm + 2, dataWindow.max.y + 1));
1357 km2 = std::max(dataWindow.min.z, std::min(km + 2, dataWindow.max.z));
1360 data.
v(im_1, jm_1, km),
1361 data.
v(im_1, jm_1, km1),
1362 data.
v(im_1, jm_1, km2), t.z);
1364 data.
v(im_1, jm, km),
1365 data.
v(im_1, jm, km1),
1366 data.
v(im_1, jm, km2), t.z);
1368 data.
v(im_1, jm1, km),
1369 data.
v(im_1, jm1, km1),
1370 data.
v(im_1, jm1, km2), t.z);
1372 data.
v(im_1, jm2, km),
1373 data.
v(im_1, jm2, km1),
1374 data.
v(im_1, jm2, km2), t.z);
1377 data.
v(im, jm_1, km),
1378 data.
v(im, jm_1, km1),
1379 data.
v(im, jm_1, km2), t.z);
1382 data.
v(im, jm, km1),
1383 data.
v(im, jm, km2), t.z);
1385 data.
v(im, jm1, km),
1386 data.
v(im, jm1, km1),
1387 data.
v(im, jm1, km2), t.z);
1389 data.
v(im, jm2, km),
1390 data.
v(im, jm2, km1),
1391 data.
v(im, jm2, km2), t.z);
1394 data.
v(im1, jm_1, km),
1395 data.
v(im1, jm_1, km1),
1396 data.
v(im1, jm_1, km2), t.z);
1398 data.
v(im1, jm, km),
1399 data.
v(im1, jm, km1),
1400 data.
v(im1, jm, km2), t.z);
1402 data.
v(im1, jm1, km),
1403 data.
v(im1, jm1, km1),
1404 data.
v(im1, jm1, km2), t.z);
1406 data.
v(im1, jm2, km),
1407 data.
v(im1, jm2, km1),
1408 data.
v(im1, jm2, km2), t.z);
1411 data.
v(im2, jm_1, km),
1412 data.
v(im2, jm_1, km1),
1413 data.
v(im2, jm_1, km2), t.z);
1415 data.
v(im2, jm, km),
1416 data.
v(im2, jm, km1),
1417 data.
v(im2, jm, km2), t.z);
1419 data.
v(im2, jm1, km),
1420 data.
v(im2, jm1, km1),
1421 data.
v(im2, jm1, km2), t.z);
1423 data.
v(im2, jm2, km),
1424 data.
v(im2, jm2, km1),
1425 data.
v(im2, jm2, km2), t.z);
1437 p.setValue(clampedVsP.x - 0.5 , clampedVsP.y - 0.5, vsP.z);
1440 c.x =
static_cast<int>(floor(p.x));
1441 c.y =
static_cast<int>(floor(p.y));
1442 c.z =
static_cast<int>(floor(p.z));
1448 im = std::max(dataWindow.min.x, std::min(c.x, dataWindow.max.x));
1449 jm = std::max(dataWindow.min.y, std::min(c.y, dataWindow.max.y));
1450 km = std::max(dataWindow.min.z, std::min(c.z, dataWindow.max.z + 1));
1451 int im_1, jm_1, km_1;
1452 im_1 = std::max(dataWindow.min.x, std::min(im - 1, dataWindow.max.x));
1453 jm_1 = std::max(dataWindow.min.y, std::min(jm - 1, dataWindow.max.y));
1454 km_1 = std::max(dataWindow.min.z, std::min(km - 1, dataWindow.max.z + 1));
1456 im1 = std::max(dataWindow.min.x, std::min(im + 1, dataWindow.max.x));
1457 jm1 = std::max(dataWindow.min.y, std::min(jm + 1, dataWindow.max.y));
1458 km1 = std::max(dataWindow.min.z, std::min(km + 1, dataWindow.max.z + 1));
1460 im2 = std::max(dataWindow.min.x, std::min(im + 2, dataWindow.max.x));
1461 jm2 = std::max(dataWindow.min.y, std::min(jm + 2, dataWindow.max.y));
1462 km2 = std::max(dataWindow.min.z, std::min(km + 2, dataWindow.max.z + 1));
1465 data.
w(im_1, jm_1, km),
1466 data.
w(im_1, jm_1, km1),
1467 data.
w(im_1, jm_1, km2), t.z);
1469 data.
w(im_1, jm, km),
1470 data.
w(im_1, jm, km1),
1471 data.
w(im_1, jm, km2), t.z);
1473 data.
w(im_1, jm1, km),
1474 data.
w(im_1, jm1, km1),
1475 data.
w(im_1, jm1, km2), t.z);
1477 data.
w(im_1, jm2, km),
1478 data.
w(im_1, jm2, km1),
1479 data.
w(im_1, jm2, km2), t.z);
1482 data.
w(im, jm_1, km),
1483 data.
w(im, jm_1, km1),
1484 data.
w(im, jm_1, km2), t.z);
1487 data.
w(im, jm, km1),
1488 data.
w(im, jm, km2), t.z);
1490 data.
w(im, jm1, km),
1491 data.
w(im, jm1, km1),
1492 data.
w(im, jm1, km2), t.z);
1494 data.
w(im, jm2, km),
1495 data.
w(im, jm2, km1),
1496 data.
w(im, jm2, km2), t.z);
1499 data.
w(im1, jm_1, km),
1500 data.
w(im1, jm_1, km1),
1501 data.
w(im1, jm_1, km2), t.z);
1503 data.
w(im1, jm, km),
1504 data.
w(im1, jm, km1),
1505 data.
w(im1, jm, km2), t.z);
1507 data.
w(im1, jm1, km),
1508 data.
w(im1, jm1, km1),
1509 data.
w(im1, jm1, km2), t.z);
1511 data.
w(im1, jm2, km),
1512 data.
w(im1, jm2, km1),
1513 data.
w(im1, jm2, km2), t.z);
1516 data.
w(im2, jm_1, km),
1517 data.
w(im2, jm_1, km1),
1518 data.
w(im2, jm_1, km2), t.z);
1520 data.
w(im2, jm, km),
1521 data.
w(im2, jm, km1),
1522 data.
w(im2, jm, km2), t.z);
1524 data.
w(im2, jm1, km),
1525 data.
w(im2, jm1, km1),
1526 data.
w(im2, jm1, km2), t.z);
1528 data.
w(im2, jm2, km),
1529 data.
w(im2, jm2, km1),
1530 data.
w(im2, jm2, km2), t.z);