Displaying 2 results from an estimated 2 matches for "ggrexp".
Did you mean:
ggrep
2005 Nov 30
1
RNG stuck via Fortran call
...<R.h>
#include <Rmath.h>
/******************/
/* random numbers */
/******************/
void F77_SUB(rndstart)(void)
{ GetRNGstate(); }
void F77_SUB(rndend)(void)
{ PutRNGstate(); }
double F77_SUB(ggrnorm)(double *mu, double *sigma)
{ return rnorm(*mu, *sigma); }
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 *lamb...
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