Question

尝试使用sourceCpp包中的Rcpp进行编译时出现以下错误：

`my path to R/.../Rcpp/internal/Exporter.h`
no matching function for call to 'arma::Cube<double>::Cube(SEXPREC*&)'

对象cube是armadillo中与array相当的R。

编辑：请注意，问题似乎是函数无法接受arma::cube对象作为参数。如果我们将arma::cube B更改为arma::mat B，则确实有效：

#include <RcppArmadillo.h>
// [[Rcpp::depends(RcppArmadillo)]]

using namespace arma;

//  [[Rcpp::export]]
arma::cube ssmooth(arma::mat  A, 
                   arma::cube B) {

int ns = A.n_rows;    
int nk = A.n_cols;    
int np = B.n_rows;    

arma::mat C = zeros<mat>(nk, ns);
arma::cube D = zeros<cube>(nk, nk, ns);

return D;

}

我很感激任何提示。

Answer 1

基本示例有效：

R> cppFunction("arma::cube getCube(int n) { arma::cube a(n,n,n);\
                    a.zeros(); return a; }", depends="RcppArmadillo")
R> getCube(2)
, , 1

     [,1] [,2]
[1,]    0    0
[2,]    0    0

, , 2

     [,1] [,2]
[1,]    0    0
[2,]    0    0

R>

所以你要么做错了，要么你的装置已关闭。

Answer 2

我有同样的问题。问题似乎与组合“Rcpp :: export”和cube作为导出函数的参数有关。我的猜测是从sexp到cube的转换器可能还没有实现（没有双关语;-)）。（或者我们都错过了一些东西......）。

当你想在Rcpp :: export函数中有一个arma :: cube参数时解决方法：先将它作为NumericVector获取，然后再创建一个立方体......

#include <RcppArmadillo.h>
// [[Rcpp::depends(RcppArmadillo)]]

using namespace arma;

//  [[Rcpp::export]]
arma::cube ssmooth(arma::mat  A, 
                   NumericVector B_) {
  IntegerVector dimB=B_.attr("dim");
  arma::cube B(B_.begin(), dimB[0], dimB[1], dimB[2]);

  //(rest of your code unchanged...)
  int ns = A.n_rows;    
  int nk = A.n_cols;    
  int np = B.n_rows;    

  arma::mat C = zeros<mat>(nk, ns);
  arma::cube D = zeros<cube>(nk, nk, ns);

  return D;

}

Answer 3

我认为你的代码失败了，因为它隐式地尝试像这样进行转换：

#include<RcppArmadillo.h>
using namespace Rcpp;
using namespace arma;
// [[Rcpp::depends(RcppArmadillo)]]
// [[Rcpp::export]]
cube return_thing(SEXP thing1){
  cube thing2 = as<cube>(thing1);
  return thing2;
}


/***R
thing <- 1:8
dim(thing) <- c(2, 2, 2)
return_thing(thing)
*/

哪个不起作用，而它适用于矩阵：

#include<RcppArmadillo.h>
using namespace Rcpp;
using namespace arma;
// [[Rcpp::depends(RcppArmadillo)]]

//[[Rcpp::export]]
mat return_thing(SEXP thing1){
  mat thing2 = as<mat>(thing1);
  return thing2;
}


/***R
thing <- 1:4
dim(thing) <- c(2, 2)
return_thing(thing)
*/

Answer 4

我能够读取并返回具有以下功能的 arma 立方体：

#include <RcppArmadillo.h>

using namespace Rcpp;
using namespace arma;

// [[Rcpp::export]]
arma::cube return_cube(arma::cube X)
{
  return(X);
}

例如，当我在 R 中运行以下命令时，我得到以下结果：

my_cube <- array(data = rnorm(5 * 3 * 2), dim = c(5,3, 2))
return_cube(my_cube)
, , 1

          [,1]       [,2]       [,3]
[1,] 0.4815994  1.0863765  0.3278728
[2,] 1.4138699 -0.7809922  0.8341867
[3,] 0.6555752 -0.2708001  0.7701501
[4,] 1.1447104 -1.4064894 -0.2653888
[5,] 1.5972670  1.8368235 -2.2814959

, , 2

             [,1]       [,2]       [,3]
[1,] -0.485091067  1.1826162 -0.3524851
[2,]  0.227652584  0.3005968 -0.6079604
[3,] -0.147653664  1.3463318 -1.2238623
[4,]  0.067090580 -0.8982740 -0.8903684
[5,]  0.006421618 -1.7156955 -1.2813880

使用RcppArmadillo将`arma :: cube`argument传递给函数时出错

4 个答案: