-
Notifications
You must be signed in to change notification settings - Fork 0
/
AsWrap.cpp
39 lines (28 loc) · 938 Bytes
/
AsWrap.cpp
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
#include <Rcpp.h>
#include <boost/numeric/ublas/vector.hpp>
namespace Rcpp {
namespace traits {
template <typename T>
SEXP wrap(const boost::numeric::ublas::vector<T>& obj) {
const int RTYPE = Rcpp::traits::r_sexptype_traits<T>::rtype;
// Rcpp::traits::storage_type<RTYPE>::type
return Rcpp::Vector<RTYPE>(obj.begin(), obj.end());
}
template <typename T>
class Exporter<boost::numeric::ublas::vector<T>> {
typedef typename boost::numeric::ublas::vector<T> out_type;
const static int RTYPE = Rcpp::traits::r_sexptype_traits<T>::rtype;
Rcpp::Vector<RTYPE> vec;
public:
Exporter(SEXP x): vec(x) {
if (TYPEOF(x) != RTYPE)
throw std::invalid_argument("Wrong R type");
}
out_type get() {
out_type x(vec.size());
std::copy(vec.begin(), vec.end(), x.begin());
return x;
}
};
}
}