import System.Environment (getArgs)
-import Cube (tetrahedron)
-import Grid (cube_at, make_grid, zoom)
-import PolynomialArray (make_polynomial_array)
-import Tetrahedron (polynomial)
+import Grid (make_grid, zoom)
import Values (read_values_3d, write_values_1d)
mri_shape :: DIM3
mridata <- read_values_3d mri_shape in_file
let g = make_grid 1 mridata
- let polynomials = make_polynomial_array (255,255,0,23)
- [ ((i,j,k,tet), polynomial t) | i <- [0..255],
- j <- [0..255],
- k <- [0],
- tet <- [0..23],
- let c = cube_at g i j k,
- let t = tetrahedron c tet ]
-
- let output = zoom g polynomials (8,8,1)
+ let output = zoom g (4,4,1)
write_values_1d output out_file