其他分享
首页 > 其他分享> > c – Rc包中的ARMA_NO_DEBUG与RcppArmadillo

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