* Add all the new generators from Knuth's 3rd Edition 

* Add SWNS generator Phys.Rev.E 50 (2) p. 1607-1615 (1994), Phys.Rev.E
60 (6), p.7626-7628 (1999)

* Add get_array type methods which can provide optimized versions of
each function (?). We should offer the possibility of eliminating
function call overhead -- there are various possible ways to implement
it.

* Add ISAAC generator (??)

* Add (A)RC4 ? (and maybe MD5)

* Check that RANLXS will work on machines with non-standard width of
float/dbl (original has checks for DBL_MANT_DIG ..)

* mention more clearly why not all Cokus used (or recheck MT pages for
improvements)

* run the DIEHARD tests on all the generators, especially the ones we
are listing as "Simulation Quality" -- some of those are a bit old and
might fail one or two diehard tests.

* Add  NAG,       missing, gave up! 
       CDC 48-bit missing 

* Check out the bug fix to mrand48 that was made in glibc2, pr757

* Check out  the following paper,

    On the anomaly of ran1() in Monte Carlo pricing of financial
    derivatives; Akira Tajima , Syoiti Ninomiya , and Shu Tezuka ; Winter
    simulation , 1996, Page 360, from ACM

* The following papers have been published, I think we refer to them
(or could do),

    Pierre L'Ecuyer, "Tables of Linear Congruential Generators of different
    size and good lattice structure", Mathematics of Computation, Vol 68,
    No 225, Jan 1999, p249-260

    Pierre L'Ecuyer, "Tables of Maximally equidistributed combined LSFR
    generators", ibid, p261-270

* Look at this paper: I. Vattulainen, "Framework for testing random numbers
  in parallel calculations", Phys. Rev. E, 59, 6, June 1999, p7200

----------------------------------------------------------------------
DONE

x1. Improve the seeding for routines that use the LCG seed generator.
It can only generate 130,000 different initial states. We should
change it to provide 2^31 different initial states (this will also
prevent the high bits being zero). DONE - we now use a 32-bit
generator.

x8. Get the macros from the faster MT19937 generator and use them. We
need to make MT be the fastest of the simulation quality generators if
it is the default. DONE. It didn't improve the speed on other
platforms, so I just used the tricks which also worked on the pentium
(e.g. changing mag[x&1] to x&1 ? mag[1] : mag[0])
