Displaying 2 results from an estimated 2 matches for "ggrpois".
Did you mean:
ggros
2005 Nov 30
1
RNG stuck via Fortran call
...double F77_SUB(ggrexp)(double *scale)
{return rexp(*scale);}
double F77_SUB(ggrgam)(double *a, double *scale)
{return rgamma(*a, *scale);}
double F77_SUB(ggrunif)(double *a, double *b)
{return runif(*a, *b);}
double F77_SUB(ggrbinom)(double *n, double *p)
{return rbinom(*n, *p);}
double F77_SUB(ggrpois)(double *lambda)
{return rpois(*lambda);}
2005 Nov 30
2
R binomial RNG stuck at 1 via Fortran call
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