Tag: 犰狳

RcppArmadillo传递用户定义的函数

考虑下面的R代码， ## ———– R version ———– caller <- function(x=1:3, fun = "identity", …){ ## do some other stuff ## … ## then call the function eval(call(fun, x)) } fun1 <- function(x, …){ x + x } fun2 <- function(x, a = 10) a * x caller(fun = "fun1") caller(fun = "fun2") 用户可以传递函数名称“fun”，这是由caller使用的。 我希望用RcppArmadillo对象执行相同的任务（显然，作为更复杂任务的一部分）。 函数将在C++定义，用户通过引用其名称在R级select它： caller_cpp(1:3, […]

是犰狳解决（）线程安全吗？

在我的代码中，我有一个循环，在其中我构build和确定的线性系统，并试图解决它： #pragma omp parallel for for (int i = 0; i < n[0]+1; i++) { for (int j = 0; j < n[1]+1; j++) { for (int k = 0; k < n[2]+1; k++) { arma::mat A(max_points, 2); arma::mat y(max_points, 1); // initialize A and y arma::vec solution = solve(A,y); } } } 有时候，程序会非常随机地挂起来，或者解决scheme向量中的结果是NaN。 如果我这样做： […]