R binomial RNG stuck at 1 via Fortran call

Wed, Nov 30, 2005 1:51 AM

Hi, 

I have some trouble with the result of a fortran function calling the R
binomial RNG:

The C function rbinom is wrapped as in the file attached.

My main fortran program starts as
call rndstart()

and ends as
call rndend()  

I happen to call the binomial RNG within a loop as 

b = ggrbinom(1.d0,0.5d0)
write(*,*) 'b=',b 

In certain cases, after a few iterations in the loop,
b get stuck at 1

Any hint to explain that would help.

Gilles

???????????????????????????????????????????????????????????????????
Gilles Guillot 
INRA - MIA Paris 

Currently working from G?teborg Stochastic Centre
Eklandagatan 86 - Rum1439 
Chalmers University of Technology
SE 412-96 G?teborg Sweden 
Phone +46 31 772 3514 / Fax +46 31772 3508

www.inapg.inra.fr/ens_rech/mathinfo/personnel/guillot/welcome.html
???????????????????????????????????????????????????????????????????