j <- [0..ymax],
k <- [0..zmax],
let tet_vol = (1/24)*(delta^(3::Int)),
+ let fvs' = make_values fvs i j k,
let cube_ijk =
- Cube delta i j k (make_values fvs i j k) tet_vol]
+ Cube delta i j k fvs' tet_vol]
where
xmax = xsize - 1
ymax = ysize - 1
o' = (fromIntegral o) / (fromIntegral sfz) - offset
p = (m', n', o') :: Point
cube = find_containing_cube g p
- -- Figure out i,j,k without importing those functions.
t = find_containing_tetrahedron cube p
f = polynomial t