Displaying 1 result from an estimated 1 matches for "rpostlamb".
2005 Nov 30
1
RNG stuck via Fortran call
...ng a Markow chain in Fortran interfaced with R.
Each loop of my Fortran calls various functions of the R RNG through
the wrapper given below.
In a run of 100 iterations of the Markov kernel,
after 20 iterations, the RNG seems like frozen.
For example, the first call to the RNG in my loop is:
rpostlamb = ggrgam(dble(m),1.d0)
write(*,*) 'rpostlamb=',rpostlamb
after the 20th iteration, it will return the same value
rpostlamb= 1.24634557
for all the following interations.
Is there any suggestion of explanation for this strange fact ?
Thanks in advance
Gilles
R Version 2.2.0 comp...