src/Main.hs

   1 module Main
   2 where
   3
   4 import Data.Array.Repa (
   5   DIM3,
   6   Z(..),
   7   (:.)(..),
   8   )
   9
  10 import Values
  11 import Grid(make_grid, zoom)
  12
  13 mri_shape :: DIM3
  14 mri_shape = (Z :. 256 :. 256 :. 109)
  15
  16 main :: IO ()
  17 main = do
  18   mridata <- read_values_3d mri_shape "./data/mridata.txt"
  19   let g = make_grid 1 mridata
  20   let output = zoom g 1
  21   write_values_1d output "output.txt"