Re: BZDEV: Re: RNGs

From: Luc Bourhis (luc.bourhis@durham.ac.uk)
Date: Wed Jan 27 1999 - 10:31:00 EST


Allan Stokes [cbi] wrote:
>
> I'll do it if no-one else volunteers.
>
> Luc Bourhis wrote:
> >I did some researchs on the web and it appears that Makoto Matsumoto and
> >Takuji Nishimura have created a new generator even better than their
> >previous TT800 currently implemented in Blitz++ : see
> >http://www.math.keio.ac.jp/~matumoto/emt.html
>
> The "cokus.c" version looks promising for incorporating into
> Blitz. If anyone wants to repackage the C code as a blitz-style
> RNG, I'd be happy to fold it into the sources. Otherwise it'll
> go on the TODO.
>
> http://www.math.keio.ac.jp/~matumoto/cokus.c
>
> Cheers,
> Todd
>

warning : the main difference between the "cokus" version and the
original version is the use of "*a" coding instead of "a[i]". We know
that the interest of this trick is platform and compiler dependent. So,
the implementer of this class should use the preprocessor flag
introduced by Todd or Blitz++ array classes and their iterators. For
example on my MacOS powerPC box with the compiler Codewarrior Pro 4, the
speed difference between these two versions is very small, with a small
advantage to the original one.

-- 
                     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