X-Git-Url: http://gitweb.michael.orlitzky.com/?a=blobdiff_plain;f=src%2FValues.hs;h=070abbbf2554ea5eec5276e5496fe1554fcaf3d3;hb=638f544b522157f83b075e07a70b2f106e01ec68;hp=2efd80b6f8b167e58dde3efe1250cf48f9df79db;hpb=c954154b379c5bd444d527298c33142fb150711b;p=spline3.git diff --git a/src/Values.hs b/src/Values.hs index 2efd80b..070abbb 100644 --- a/src/Values.hs +++ b/src/Values.hs @@ -12,7 +12,7 @@ import Data.Array.Repa ( DIM3, extent, fromList, - index, + unsafeIndex, reshape, size ) @@ -20,7 +20,7 @@ import Data.Array.Repa ( 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 @@ -33,23 +33,16 @@ type Values3D = Array DIM3 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 = readVectorFromTextFile @@ -79,7 +72,7 @@ dims v3d = 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)