c – Rc包中的ARMA_NO_DEBUG与RcppArmadillo
作者:互联网
我想在访问RcppArmadillo中的矩阵元素时禁用绑定检查.
犰狳的文件说
Armadillo can be configured via editing the file
include/armadillo_bits/config.hpp. Specific functionality can be
enabled or disabled by uncommenting or commenting out a particular
#define, listed below.
但是在R包的上下文中,我该如何激活该指令?
我试图用.创建一个config.h文件
#ifndef CONFIG_LOADED
#define CONFIG_LOADED
#define ARMA_NO_DEBUG
#endif
然后将它包含在my / src文件夹的每个.cpp文件中,但我不确定它是否正常工作,或者除了在每个.cpp文件中添加#include“config.h”之外还有其他方法.
目前我有一个.cpp(包含主算法的那个),它以:
#include "configs.h"
#include <RcppArmadillo.h>
using namespace Rcpp;
using namespace arma;
// [[Rcpp::export]]
SEXP sample_gibbs_cpp(const arma::vec& v_n, const arma::mat& W,
arma::vec h_n, double alpha = 1, double beta = 1, int iter=100,
double burnin = 0.5){
... code ...
}
然后是其他人
#include <RcppArmadillo.h>
using namespace Rcpp;
using namespace arma;
... code ...
我的描述文件:
Package: mypackage
Title: What the Package Does (one line, title case)
Version: 0.0.0.9000
Authors@R: person("First", "Last", email = "first.last@example.com", role = c("aut", "cre"))
Description: What the package does (one paragraph).
Depends:
R (>= 3.2.3)
License: What license is it under?
Encoding: UTF-8
LazyData: true
RoxygenNote: 5.0.1
Imports:
ggplot2,
dplyr,
tidyr,
rstan
LinkingTo: Rcpp, RcppArmadillo, RcppEigen
SystemRequirements: C++11
我编译我的包:
devtools::load_all()
解决方法:
订单在这里很重要在包含#include< RcppArmadillo.h>之前,必须包含此#define语句.
一个例子:
custom_config.h
#ifndef CONFIG_LOADED
#define CONFIG_LOADED
#define ARMA_NO_DEBUG
#endif
example_compiled_file.cpp
#include "custom_config.h"
#include <RcppArmadillo.h>
// [[Rcpp::export]]
void test_pkg(const arma::vec& x) {
// Should not trigger error bound checking with debug flag on.
double my_val_protected = x(0);
// Never triggers error bound checking
double my_val = x.at(0);
}
注意:由于这是一个包,因此不需要使用// [[Rcpp :: depends(RcppArmadillo)]].相反,您必须在DESCRIPTION文件的LinkingTo:字段中指定RcppArmadillo和Rcpp,并在Imports:字段中包含Rcpp.您必须从Rcpp最小程度地导入一个函数(最好是:evalCpp).
例如说明必须具备:
Imports: Rcpp (>= 0.12.15)
LinkingTo: Rcpp, RcppArmadillo
标签:c,r,rcpp,armadillo 来源: https://codeday.me/bug/20190828/1748640.html