alpinist
03-21-2007, 09:15 PM
Hi,dear everyone, if I'd like to get random numbers of double precision, can I just change the "real" to "double precision" in subroutine ran0,or ran1,etc.
Hi,dear everyone, if I'd like to get random numbers of double precision, can I just change the "real" to "double precision" in subroutine ran0,or ran1,etc.
For the Fortran 90 recipes take a look at
http://www.numerical-recipes.com/forum/showthread.php?t=300
I suppose you will need to make a similar change for amm in subroutine ran_init in module ran_state as well as changing real to double precision as you have mentioned.