[Rd] Problem with rnorm ?

[Petr Savicky]

> #include <stdio.h>
> #include <Rmath.h>
> 
> int main(void){
>   double x[3];
>   for(int i=0;i<3;i++)  x[i]=rnorm(0,1);
>   printf("%lf,%lf,%lf",x[0],x[1],x[2]);
>   return 0;
> }
> 
> output :   -8.773321,-8.773321,-8.773321

You probably have to call GetRNGstate() before rnorm and PutRNGstate() after.
At least for extension packages, this is necessary, see section 6.3 in R-exts.pdf.

Petr.