module Main
where

import Data.Array.Repa (
  DIM3,
  Z(..),
  (:.)(..),
  )

import Values
import Grid(make_grid, zoom)

mri_shape :: DIM3
mri_shape = (Z :. 256 :. 256 :. 109)

main :: IO ()
main = do
  mridata <- read_values_3d mri_shape "./data/mridata.txt"
  let g = make_grid 1 mridata
  let output = zoom g 1
  write_values_1d output "output.txt"