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.