Hello,

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

vdrnguniform

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....