使用 .C() 从 R 调用 C 代码：随机值每次都相同

Question

我对这个提问者有相反的问题：Passing seed/Setting seed/ C within R code。我的问题是我想使用 R 随机数生成器在 C 代码中生成一个随机数，但是每次运行它时，我都会得到相同的“随机”值。这是我的代码（包括我碰巧正在使用的头文件，这些对于代码 sn-p 来说并非都是必需的）：

#include <stdio.h>
#include <stdlib.h>
#include <malloc.h>
#include <math.h>   // C header files

#include <R.h>
#include <R_ext/Utils.h>
#include <Rmath.h>

void main(){
GetRNGstate();
double u = unif_rand();
PutRNGstate();
Rprintf("%f\n",u);
}

现在我在使用 R CMD SHLIB 编译并使用 dyn.load() 后从 R Windows GUI 调用它：

> .C("main")
0.000889
list()

现在如果我关闭 R，重新启动，然后 dyn.load() 再次（但没有第二次编译）：

> .C("main")
0.000889
list()

根据我在互联网上找到的各种代码和“编写 R 扩展”手册，这不应该发生。但我可能只是错过了一步？

Answer 1

我无法复制这个。

这是一个完整且可重现的示例，仅依靠 inline 包来帮助构建、编译、链接和加载最小函数：

R> library(inline)
R> src <- '
+   GetRNGstate();
+   double u = unif_rand();
+   Rprintf("%f\\n",u);
+   PutRNGstate();
+   return R_NilValue;
+ '
R> fun <- cxxfunction(signature(), body=src)
R> set.seed(42)
R> invisible(fun())   # invisible to suppress the NULL return in foo()
0.914806
R> invisible(fun())
0.937075
R> invisible(fun())
0.286140
R> invisible(fun())
0.830448
R> invisible(fun())
0.641746
R> 

如果我重置种子，同样的序列开始：

R> set.seed(42)
R> invisible(fun())
0.914806
R> invisible(fun())
0.937075
R> invisible(fun())
0.286140
R> 

您可能正在加载/运行具有与设置 RNG 相关的副作用的不同功能。此外，以防万一，您的示例无法加载，因为您无法将带有 main() 的函数加载到 R 本身中。