C++和R接口,获得输出

Question

我刚开始学习R,因此陷入了琐碎的问题.我试图找出如何在R中输出值.在C++中我们只使用返回变量,但这似乎不是R的情况.假设我有一个带有4个输入参数的函数,将这些参数传递给c ++函数执行所需的计算,现在如果我想在R中加载这个myfun并从c ++函数获取输出,我需要做什么？以下是我尝试使用的模板.

extern "C" {
   SEXP myfun(SEXP S, SEXP A, SEXP B, SEXP C) { //will call this function from R.
     SEXP rate, dir, list, list_names; //declare variables

     PROTECT( rate = allocMatrix(REALSXP, 10, 2) ); //allocate 10x2 matrix of double type?
     PROTECT( dir = allocVector(INTSXP, 10) ); //allocated vector(10) of int type?

     double* p_rate = REAL(rate); //why do I need pointers? 
     int* p_dir = INTEGER(dir);

//这里我调用一个C++函数来计算vector<vector<double> > someVal和vector<int> someVal2.

现在我想将这些值传递给rate和dir.

 for(int i =0; i < 10; i++){
      rate[i][0] = someVal1[i][0];
      rate[i][1] = someVal1[i][1];
      dir[i] = someVal2[i]; 
  }

但编译器不喜欢它.在这种情况下,我不确定如何使用p_rate和p_dir指针.

另外,我不太明白为什么下面的代码是这样的.它似乎将字符串,rate和dir传递给list_names,并将值传递给rate和dir.为什么我们不能简单地做cout并打印所需的值？

   char *tag[2] = {"rate","dir"};
   PROTECT(list_names = allocVector(STRSXP ,2));
   SET_STRING_ELT(list_names, 0, mkChar(tag[0]));
   SET_STRING_ELT(list_names, 1, mkChar(tag[1]));
   PROTECT(list = allocVector(VECSXP ,2));
   SET_VECTOR_ELT(list, 0, rate);
   SET_VECTOR_ELT(list, 1, dir);
   setAttrib(list, R_NamesSymbol , list_names);

   UNPROTECT(8);
   return (list);
}

我没有得到在R中返回变量的基本想法,这看起来很混乱.如果有人能指出适当的资源解释这些东西,我真的很感激.

感谢您的帮助!

Answer 1

如果您已经熟悉C++,您可能会发现我们用于无缝R和C++ 集成的R​​cpp项目非常有用.它包含大量文档和示例,而且CRAN现在有超过80个使用它的包,它提供了另一个大型的示例集.

您可以从与我们的JSS论文相对应的"简介"插图开始,并在[rcpp]此处仔细阅读标记.

否则,你回到R API,这是更简单的基于C的.

Answer 2

正如我们在回答您之前提出的问题时所说的那样,给Rcpp一个机会,然后拿走代码中的所有kludge.

#include <Rcpp.h>
using namespace Rcpp ;

void yourCode( std::vector< std::vector<double> >& someVal, 
               std::vector<double>& someVal2 ){ ... }

extern "C" SEXP myfun() { 
    NumericMatrix rate(10, 2); 
    NumericVector dir(10) ;

    std::vector< std::vector<double> > someVal ;
    std::vector< double > someVal2 ;

    // call whatever function that fills someVal and someVal2
    yourCode( someVal, someVal2 ) ;

    // fills rate and dir
    for(int i =0; i < 10; i++){
      rate(i, 0) = someVal1[i][0];
      rate(i, 1) = someVal1[i][1];
      dir[i] = someVal2[i]; 
    }

    // structure the output
    List result = List::create( 
        _["rate"] = rate,
        _["dir"] = dir
    ) ;

    return result ;

}