我在C 中写了一个类class_a,并使用RCPP模块框架将其暴露于R。
#include "Rcpp.h"
using namespace Rcpp;
class class_A {
public:
class_A(double num){this->num = num;};
double get_num(){return this->num;};
private:
double num;
};
RCPP_MODULE(class_A_module) {
class_<class_A>("class_A")
.constructor<double>()
.property("num", &class_A::get_num)
;
}
RCPP_EXPOSED_CLASS(class_A);
现在,我想拥有一个可以存储class_a的多个实例(以及其他属性(的类class_b。我的想法是为此目的使用std :: vector:
class class_B {
public:
class_B(std::vector<class_A> num_list){this->num_list = num_list;};
std::vector<class_A> get_num_list(){return this->num_list;};
private:
std::vector<class_A> num_list;
};
RCPP_MODULE(class_B_module) {
class_<class_B>("class_B")
.constructor<std::vector<class_A>>()
.property("num_list", &class_B::get_num_list)
;
}
RCPP_EXPOSED_CLASS(class_B);
不编译:
no matching function for call to 'export_range__dispatch(SEXPREC*&,
__gnu_cxx::__normal_iterator<class_A*, std::vector<class_A>
>&, Rcpp::traits::r_type_traits<class_A>::r_category)'
这并不奇怪:RCPP应该如何将std ::向量展示到r?
从这个高度相关的RCPP邮件邮件列表线程中,我知道,尽管如此,可以以某种方式创建此简单的存储class_b。我想我必须编写自己的wrap
和as
实现。
1。到目前为止这是正确的吗?
2。在我的简单代码示例的情况下,我该如何确切地执行此操作?我发现所有的例子都比我预期的最低实施要复杂得多。
编辑:我修改了示例代码,并试图查明问题。
标题和实现现在分开:
ex.h
class class_A {
public:
class_A(double num);
double get_num();
private:
double num;
};
class class_B {
public:
class_B(std::vector<class_A> num_list);
std::vector<class_A> get_num_list();
private:
std::vector<class_A> num_list;
};
ex.cpp
#include "Rcpp.h"
#include "ex.h"
using namespace Rcpp;
class_A::class_A(double num){this->num = num;}
double class_A::get_num(){return this->num;}
RCPP_MODULE(class_A_module) {
class_<class_A>("class_A")
.constructor<double>()
.property("num", &class_A::get_num)
;
}
RCPP_EXPOSED_CLASS(class_A);
class_B::class_B(std::vector<class_A> num_list){this->num_list = num_list;}
std::vector<class_A> class_B::get_num_list(){return this->num_list;}
RCPP_MODULE(class_B_module) {
class_<class_B>("class_B")
.constructor<std::vector<class_A>>()
.property("num_list", &class_B::get_num_list)
;
}
RCPP_EXPOSED_CLASS(class_B);
这是我尝试实现"包装"one_answers" AS"的尝试。
fix.cpp
#include <RcppCommon.h>
#include "ex.h"
// non-intrusive extension via template specialisation
namespace Rcpp {
template <> class_A as(SEXP aa);
template <> SEXP wrap(const class_A &a);
template <> class_B as(SEXP bb);
template <> SEXP wrap(const class_B &b);
}
#include <Rcpp.h>
// define template specialisations for as and wrap
namespace Rcpp {
template <> class_A as(SEXP aasexp) {
double aa = as<double>(aasexp);
return class_A(aa);
}
template <> SEXP wrap(const class_A &a) {
return Rcpp::wrap(a);
}
template <> class_B as(SEXP bbsexp) {
// ?
}
template <> SEXP wrap(const class_B &b) {
// ?
}
}
我猜有趣的部分是我无法转过头的?
一种方法是让构造函数服用Rcpp::List
或工厂,如果您不想弄乱class_B
的语义:
class_B* class_B_ctor( List data){
std::vector<class_A> v ;
for( int i=0; i<data.size(); i++){
v.push_back( as<class_A>(data[i]) );
}
return new class_B(v) ;
}
您可以与模块中的.factory
方法一起使用:
RCPP_MODULE(class_A_module) {
class_<class_A>("class_A")
.constructor<double>()
.property("num", &class_A::get_num)
;
class_<class_B>("class_B")
.factory< List >( &class_B_ctor )
.property("num_list", &class_B::get_num_list)
;
}
这样:
> a1 <- new( class_A )
> a2 <- new( class_B )
> a1
C++ object <0x103008b78> of class 'class_A' <0x10e4bcd20>
> a2
C++ object <0x103008b78> of class 'class_B' <0x10e4bce20>
> new( class_B, list(a1,a2) )
C++ object <0x10e4de420> of class 'class_B' <0x10e4bce20>