-{-# LANGUAGE TypeSynonymInstances #-}
+{-# LANGUAGE FlexibleInstances #-}
module Values
where
fromList,
index,
reshape,
+ size
)
-import Data.Array.Repa.IO.Vector (readVectorFromTextFile)
+import Data.Array.Repa.IO.Vector (readVectorFromTextFile,
+ writeVectorToTextFile)
import System.FilePath ()
import Test.QuickCheck (Arbitrary(..), Gen)
+import ScaleFactor
+
+
type Values1D = Array DIM1 Double
type Values2D = Array DIM2 Double
type Values3D = Array DIM3 Double
read_values_1d :: FilePath -> IO Values1D
-read_values_1d path = readVectorFromTextFile path
+read_values_1d = readVectorFromTextFile
read_values_3d :: DIM3 -> FilePath -> IO Values3D
one_d <- read_values_1d path
return $ reshape sh one_d
+write_values_1d :: Values3D -> FilePath -> IO ()
+write_values_1d v3d path = do
+ let size3d = size $ extent v3d
+ let shape1d = (Z :. size3d)
+ let v1d = reshape shape1d v3d
+ writeVectorToTextFile v1d path
empty3d :: Values3D
empty3d = Data.Array.Repa.fromList (Z :. 0 :. 0 :. 0) []
where
shape :: DIM3
shape = (Z :. i :. j :. k)
+
+
+zoom_shape :: ScaleFactor -> DIM3 -> DIM3
+zoom_shape (sfx, sfy, sfz) sh =
+ let (Z :. x :. y :. z) = sh
+ x' = x * sfx
+ y' = y * sfy
+ z' = z * sfz
+ in
+ (Z :. x' :. y' :. z')