From 27ebc2cd08a0a65522bb571c4866aa7f41875ac3 Mon Sep 17 00:00:00 2001
From: Alberto Ruiz <aruiz@um.es>
Date: Fri, 9 Oct 2015 18:18:44 +0200
Subject: mwc-random example (wip)

---
 examples/random.hs | 23 +++++++++++++++++++++++
 1 file changed, 23 insertions(+)
 create mode 100644 examples/random.hs

(limited to 'examples/random.hs')

diff --git a/examples/random.hs b/examples/random.hs
new file mode 100644
index 0000000..9be8eca
--- /dev/null
+++ b/examples/random.hs
@@ -0,0 +1,23 @@
+import System.Random.MWC
+import qualified System.Random.MWC.Distributions as D
+import Numeric.LinearAlgebra
+import Numeric.LinearAlgebra.Devel
+  
+rvec :: Vector Double
+rvec = runSTVector $ do
+    v <- newUndefinedVector 10
+    g <- initialize (fromList [4, 8, 15, 16, 23, 42])
+    mapM_ (\k -> writeVector v k =<< D.standard g) [0..9]
+    return v
+
+
+main = do
+    v <- withSystemRandom . asGenST $ \gen -> uniformVector gen 20
+    print (v :: Vector Double)
+
+    g <- initialize (fromList [4, 8, 15, 16, 23, 42])
+    x <- uniform g :: IO Double
+    print x
+    print =<< (uniform g :: IO Double)
+    print =<< (uniformVector g 20 :: IO (Vector Double))
+
-- 
cgit v1.2.3