0

我正在尝试在 Rcpp 中实现粒子过滤器并使用 OpenMP 来并行化转换步骤。我正在使用这里的分发函数dqrng来创建线程安全的 RNG 。boost

R API 的代码可以在这里找到并在此处介绍dqrng

我遇到的问题是，使用 R API，我获得了正确的结果，并针对替代实现进行了验证，估计器的密度按预期大致呈正态分布。但是，对于该dqrng版本，估计器的密度似乎不正确，并获得了不同的结果。密度图如下所示。

dqrng 估计器的密度

使用 R API 的估计器密度

有谁知道为什么会这样？

4

0 回答 0