-module Values
+{-# LANGUAGE FlexibleInstances #-}
+
+module Values (
+ Values,
+ Values2D,
+ Values3D,
+ dims,
+ empty3d,
+ idx,
+ zoom_shape
+ )
where
import Data.Array.Repa (
- Array,
- DIM1,
+ Array, U,
+ Z(..),
+ (:.)(..),
+ DIM2,
DIM3,
- reshape,
+ extent,
+ fromListUnboxed,
+ unsafeIndex,
)
-import Data.Array.Repa.IO.Vector (readVectorFromTextFile)
import System.FilePath ()
+import Test.QuickCheck (Arbitrary(..), Gen, choose, vectorOf)
+
+
+import ScaleFactor (ScaleFactor)
+
+
+type Values sh = Array U sh Double
+type Values2D = Values DIM2
+type Values3D = Values DIM3
+
+
+instance Arbitrary Values3D where
+ arbitrary = do
+ -- 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 = Data.Array.Repa.fromListUnboxed new_shape elements
+ return three_d
+
+
+empty3d :: Values3D
+empty3d = Data.Array.Repa.fromListUnboxed (Z :. 0 :. 0 :. 0) []
+
+
+dims :: Values3D -> (Int, Int, Int)
+dims v3d =
+ let (Z :. x :. y :. z) = extent v3d
+ in
+ (x,y,z)
-type Values1D = Array DIM1 Double
-type Values3D = Array DIM3 Double
+idx :: Values3D -> Int -> Int -> Int -> Double
+idx v3d i j k =
+ unsafeIndex v3d shape
+ where
+ shape :: DIM3
+ shape = (Z :. i :. j :. k)
-read_values_1d :: FilePath -> IO Values1D
-read_values_1d path = readVectorFromTextFile path
-read_values_3d :: DIM3 -> FilePath -> IO Values3D
-read_values_3d sh path = do
- one_d <- read_values_1d path
- return $ reshape sh one_d
+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')