Use R random number generator inside Fortran

[cjvanlissa]

Here is some code that I found, that might be helpful:

myrandom.c:

#include <R.h>
#include <Rmath.h>

void F77_SUB(fseedi)(void)
{
   GetRNGstate();
}


void F77_SUB(fseedo)(void)
{
   PutRNGstate();
}


void F77_SUB(myrbeta)(double* px, double* pa, double* pb)
{
        *px = rbeta(*pa, *pb);
}


and example.f

       subroutine example(x, a, b)
       implicit double precision (a-h,o-z)
       call fseedi()
       call myrbeta(x, a, b)
       call fseedo()
       end