我对Rcpp
或者更具体地说是RcppEigen
非常新,正在努力学习如何在我的代码中使用pi
作为一个常数。代码在MCMC算法中运行多次,因此任何速度的提升都是完美的。目前,每当我调用函数时,都手动传递pi,就像下面的代码一样:
require(RcppEigen)
require(inline)
I.Cpp <- "
using Eigen::Map;
using Eigen::MatrixXd;
using Eigen::VectorXd;
using Rcpp::NumericVector;
const Map<MatrixXd> delta(as<Map<MatrixXd> >(delta0));
const Map<VectorXd> d(as<Map<VectorXd> >(DD));
const Rcpp::NumericVector tpi(pie);
double pi = tpi[0];
const MatrixXd I = delta.transpose() * d.asDiagonal() * pi * pi;
return wrap(I);
"
I.cpp <- cxxfunction(signature(delta0 = "matrix", DD = "numeric", pie = "numeric"), I.Cpp, plugin = "RcppEigen")
delta0 <- matrix(rnorm(25), 5)
DD <- rnorm(5)
I.cpp(delta0, DD, pi) # this piece of code gets called multiple times?
我的问题是:如何在不传递的情况下在
RcppEigen
中使用常数 pi
?不需要在每次调用 I.cpp
时都传递它。
M_PI
替换pi
,它就可以工作了。(我猜我仍然没有完全理解我可以从Rcpp中访问什么,但显然这是没有特定包含语句的) - HenrikM_PI
是基本的C语言常量用法,大约在20世纪70年代末。关于“我从哪里获取pi值”的问题与Rcpp或RcppEigen无关。首先要做的事情是 :) - Dirk Eddelbuettel