-{-# LANGUAGE TypeSynonymInstances #-}
+{-# LANGUAGE FlexibleInstances #-}
module Values
where
import Data.Array.Repa.IO.Vector (readVectorFromTextFile,
writeVectorToTextFile)
import System.FilePath ()
-import Test.QuickCheck (Arbitrary(..), Gen)
+import Test.QuickCheck (Arbitrary(..), Gen, choose, vectorOf)
+
+
+import ScaleFactor
type Values1D = Array DIM1 Double
instance Arbitrary Values3D where
arbitrary = do
- x_dim <- arbitrary :: Gen Int
- y_dim <- arbitrary :: Gen Int
- z_dim <- arbitrary :: Gen Int
- one_d <- arbitrary :: Gen Values1D
+ -- I declare not to care about empty lists.
+ x_dim <- choose (1, 27)
+ y_dim <- choose (1, 27)
+ z_dim <- choose (1, 27)
+ elements <- vectorOf (x_dim * y_dim * z_dim) (arbitrary :: Gen Double)
let new_shape = (Z :. x_dim :. y_dim :. z_dim)
- let three_d = reshape new_shape one_d
+ let three_d = Data.Array.Repa.fromList new_shape elements
return three_d
-instance Arbitrary Values1D where
- arbitrary = do
- x <- arbitrary :: Gen [Double]
- let shape = (Z :. (length x))
- let one_d = Data.Array.Repa.fromList shape x
- return one_d
-
-
read_values_1d :: FilePath -> IO Values1D
-read_values_1d path = readVectorFromTextFile path
+read_values_1d = readVectorFromTextFile
read_values_3d :: DIM3 -> FilePath -> IO Values3D
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')