How to advance the RNG state when sampling random distributions using hmatrix-gsl-stats?

Question

How do I advance the RNG in the hmatrix-gsl-stats when I'm sampling from a random distribution: For example:

λ> CN.random_1p Exponential 235252 1.0
0.8742160459331277
λ> CN.random_1p Exponential 235252 1.0
0.8742160459331277

I don't see an obvious interface to maintain the state of the RNG in https://hackage.haskell.org/package/hmatrix-gsl-stats-0.2.1/docs/Numeric-GSL-Distribution-Continuous.html

score 3 · Accepted Answer · answered May 06 '15 at 22:43

3

The API in that package doesn't seem to offer any operations that return updated seeds (which is obviously a big no-no when doing randomness in a pure language). You will have to patch the package to return updated seeds or use a different package.

answered May 06 '15 at 22:43

Daniel Wagner

145,880
9
220
380

1

that's a bummer, doesn't anyone actually _use_ these libraries!? – daj May 07 '15 at 03:44

How to advance the RNG state when sampling random distributions using hmatrix-gsl-stats?

1 Answers1