gfortran - random number generator in fortran -
gfortran - random number generator in fortran -
i testing rng code mentioned in link (fortran 77 code): https://www.cisl.ucar.edu/zine/96/spring/articles/3.random-6.html applies park & miller algorithm.
the function phone call of programme is
call srand(seed) x=rand()
however programme not seem respond seeding srand(), i.e., x-value unaffected initial seeding, , equal ~0.218. suspect has definition of mutual block , info block, value of x equivalent putting seed 123456789, initialization value defined in datablock.
any ideas?
i compiling gfortran.
ok, problem original srand
, rand
may called. renamed functions srand2
, rand2
, start work expected.
for seed = 1
result 7.82636926e-06
, seed = 2
result 1.56527385e-05
. checked in gfortran , in intel's fortran.
fortran gfortran prng fortran-common-block
Comments
Post a Comment