bug-gsl
[Top][All Lists]
Advanced

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

gsl_ran_gaussian_ziggurat doesn't work with gsl_rng_slatec and gsl_rng_u


From: camel-cdr
Subject: gsl_ran_gaussian_ziggurat doesn't work with gsl_rng_slatec and gsl_rng_uni
Date: Fri, 05 Feb 2021 07:52:54 +0000

I've been developing a ziggurat implementation and used 
https://www.seehuhn.de/pages/ziggurat to benchmark mine against the one used in 
gsl.
I ran[timegauss.c,](https://m.seehuhn.de/data/ziggurat/timegauss.c) but it 
fails with gsl_rng_slatec and gsl_rng_uni.
Consequently, I checked randist/gausszig.c and gsl_ran_gaussian_ziggurat 
doesn't support generators with smaller ranges than 16777216.

On a related note: I'd also be interested in improving the performance of 
gsl_ran_gaussian_ziggurat outside of this bug, but I'm new to this project, and 
I'm not sure what the procedure would be. (In my current benchmarks my 
implementation is 1.3 times faster and uses 2/3 off the memory)

reply via email to

[Prev in Thread] Current Thread [Next in Thread]