Question

在我正在RStudio中编辑的.cpp文件中，以下代码

// [[Rcpp::depends(RcppArmadillo)]]

#include <RcppArmadillo.h>
using namespace Rcpp;
using namespace arma;

// [[Rcpp::export]]
mat mkGramMatrix (mat X, int n, double (*kernel) (mat, mat, double), double k_param=0){
  double (*func_kernel) (mat, mat, double);
  func_kernel = kernel;

  mat K(n, n);

  for (int i = 0; i < n; i++){
    for (int j = 0; j < n; j++){
      K(i, j) = func_kernel(X.row(i), X.row(j), k_param);
    }
  }
  return K;
}

当我尝试在RStudio中输入时，

产生以下错误：

KFDA_Rcpp.cpp:138:63: error: cannot combine with previous 'type-name' declaration specifier
RcppExport SEXP sourceCpp_32_mkGramMatrix(SEXP XSEXP, SEXP nSEXP, SEXP double)SEXP, SEXP k_paramSEXP) {
                                                                       ^
KFDA_Rcpp.cpp:138:70: error: expected ';' after top level declarator
RcppExport SEXP sourceCpp_32_mkGramMatrix(SEXP XSEXP, SEXP nSEXP, SEXP double)SEXP, SEXP k_paramSEXP) {
                                                                              ^
                                                                              ;
KFDA_Rcpp.cpp:138:74: error: expected unqualified-id
RcppExport SEXP sourceCpp_32_mkGramMatrix(SEXP XSEXP, SEXP nSEXP, SEXP double)SEXP, SEXP k_paramSEXP) {
                                                                                  ^
3 errors generated.

为什么在

中放置右括号

SEXP double)SEXP

当我相信我正确使用函数指针时？我已成功编译并使用其他RcppArmadillo函数没有问题，但这是我正在制作的第一个利用函数指针的函数，所以我可以在调用此函数时指定内核。

Answer 1

您不能在界面中放置函数指针，因为它不是映射到SEXP的类型。

我们只能从R调用映射到需要.Call()参数的SEXP接口的函数---所以我们可以使用我们有来自SEXP的映射器的类型。

这里简单的方法是使用Rcpp::XPtr。它提供了一个外部指针接口＆＃39;通过它你可以传递你的内核函数。在mkGramMatrix()函数中，您可以将XPtr取消引用到要调用的函数。这里和Rcpp Gallery上有几个帖子可以帮助您入门。它是一个伟大的机制：挂在那里并阅读它。

RcppExport中的错误原因

1 个答案: