I want to generate random number with the SFMT19937 generator. Everything works fine as long as I store the random numbers to a vector, e.g
.... brng=VSL_BRNG_SFMT19937 method=VSL_RNG_METHOD_UNIFORM_STD seed=1777 errcode=vslnewstream( stream, brng, seed ) do j = 1,1000 errcode=vdrnguniform( method, stream, n, r, lb, ub ) write(1,*) r(j) end do
with this I fill r(1000) with double real random numbers, since the argument r in
is a vector.
Is there a way just to generate them without storing them to a vector?
I want to generate several trillions of random numbers for a monte-carlo simulation.
Storing them to an array doesn't seem like a good idea....