DIM3,
extent,
fromList,
- index,
+ unsafeIndex,
reshape,
size
)
import Data.Array.Repa.IO.Vector (readVectorFromTextFile,
writeVectorToTextFile)
import System.FilePath ()
-import Test.QuickCheck (Arbitrary(..), Gen, Positive(..), vectorOf)
+import Test.QuickCheck (Arbitrary(..), Gen, choose, vectorOf)
import ScaleFactor
-type Values1D = Array DIM1 Double
-type Values2D = Array DIM2 Double
-type Values3D = Array DIM3 Double
+type Values sh = Array sh Double
+type Values1D = Values DIM1
+type Values2D = Values DIM2
+type Values3D = Values DIM3
instance Arbitrary Values3D where
arbitrary = do
- (Positive x_dim) <- arbitrary :: Gen (Positive Int)
- (Positive y_dim) <- arbitrary :: Gen (Positive Int)
- (Positive z_dim) <- arbitrary :: Gen (Positive Int)
+ -- 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.fromList new_shape elements
idx :: Values3D -> Int -> Int -> Int -> Double
idx v3d i j k =
- index v3d shape
+ unsafeIndex v3d shape
where
shape :: DIM3
shape = (Z :. i :. j :. k)
z' = z * sfz
in
(Z :. x' :. y' :. z')
+
+
+drop_z :: DIM3 -> DIM2
+drop_z (Z :. 1 :. y :. x) = (Z :. y :. x)
+drop_z _ = error "can't drop the z-dimension unless its size is 1"