此代碼是爲賦值給出的。但是,蝙蝠在計算pi蒙特卡洛函數中存在錯誤。從int到unsigned int的轉換錯誤統計無效。林不知道這是否是有意的,但我一直無法糾正這個問題。如果你想知道這個函數是線程示例的一部分。任何建議是非常感謝 void *compute_pi(void *s)
{
int seed;
int ii;
int *hit_pointer;
int local_hits;
前一段時間(我不記得在哪裏),我發現這條巨蟒片斷這implents插口算法計算圓周率的數字: def pi_digits():
"""generator for digits of pi"""
q,r,t,k,n,l = 1,0,1,1,3,3
while True:
if 4*q+r-t < n*t:
yield n
q,r,
我對Rcpp或更具體地說是RcppEigen非常新,並且在如何使用pi作爲我的代碼中的常量方面苦苦掙扎。代碼在MCMC算法中運行很多次,所以速度的改進將是完美的。目前,我交出我每次調用函數時PI,如下面的代碼: require(RcppEigen)
require(inline)
I.Cpp <- "
using Eigen::Map;
using Eigen::MatrixXd;
u