I have a Monte Carlo library written in FORTRAN that I want to port in
C++. Therefore, I need a generator of random doubles uniformly
distributed in [0,1[. Blitz++ provides two of them (classes Uniform ant
TT800) but there are several "NEEDS WORK" comments : so what are their
status ?
The random generator used by my Monte Carlo generates uniformly two
integers x and y and then builds a floating-point number z=(x+y/m1)/m0
where m0 = 1664557 and m1 = 1048583. If z is a double and assuming that
integers have 32 bits, does this algorithm gives 53 significant bits for
the significand of z ? (Sorry to ask a question not directly related to
Blitz++)
--
Luc Bourhis
Center for Particle Theory
University of Durham, UK
--------------------- blitz-dev list --------------------------------
* To subscribe/unsubscribe: mail to majordomo@oonumerics.org, with
"subscribe blitz-dev" or "unsubscribe blitz-dev" in the body of the message
* Blitz++ web page: http://oonumerics.org/blitz/
This archive was generated by hypermail 2b29 : Wed Feb 20 2002 - 04:30:08 EST