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
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
|
// -*- mode: C++; c-indent-level: 4; c-basic-offset: 4; indent-tabs-mode: nil; -*-
//
// SugarBlock.h: Rcpp R/C++ interface class library -- sugar functions
//
// Copyright (C) 2010 - 2011 Dirk Eddelbuettel and Romain Francois
//
// This file is part of Rcpp.
//
// Rcpp is free software: you can redistribute it and/or modify it
// under the terms of the GNU General Public License as published by
// the Free Software Foundation, either version 2 of the License, or
// (at your option) any later version.
//
// Rcpp is distributed in the hope that it will be useful, but
// WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
// GNU General Public License for more details.
//
// You should have received a copy of the GNU General Public License
// along with Rcpp. If not, see <http://www.gnu.org/licenses/>.
#ifndef RCPP_SUGAR_BLOCK_2_H
#define RCPP_SUGAR_BLOCK_2_H
namespace Rcpp{
namespace sugar{
template <bool NA, typename RESULT_TYPE, typename U1, typename T1, typename U2, typename T2>
class SugarBlock_2 : public Rcpp::VectorBase< Rcpp::traits::r_sexptype_traits<RESULT_TYPE>::rtype , NA, SugarBlock_2<NA,RESULT_TYPE,U1,T1,U2,T2> > {
public:
typedef RESULT_TYPE (*FunPtr)(U1,U2) ;
SugarBlock_2( FunPtr ptr_, const T1 & x_, const T2& y_ ) :
ptr(ptr_), x(x_), y(y_){
// TODO: check that x and y have same size
}
inline RESULT_TYPE operator[]( R_xlen_t i) const {
return ptr( x[i], y[i] ) ;
}
inline R_xlen_t size() const { return x.size() ; }
private:
FunPtr ptr ;
const T1& x ;
const T2& y ;
};
template <bool NA, typename RESULT_TYPE, typename U1, typename T1, typename U2>
class SugarBlock_2__VP : public Rcpp::VectorBase< Rcpp::traits::r_sexptype_traits<RESULT_TYPE>::rtype , NA, SugarBlock_2__VP<NA,RESULT_TYPE,U1,T1,U2> > {
public:
typedef RESULT_TYPE (*FunPtr)(U1,U2) ;
SugarBlock_2__VP( FunPtr ptr_, const T1 & x_, U2 u2 ) :
ptr(ptr_), x(x_), y(u2){}
inline RESULT_TYPE operator[]( R_xlen_t i) const {
return ptr( x[i], y ) ;
}
inline R_xlen_t size() const { return x.size() ; }
private:
FunPtr ptr ;
const T1& x ;
U2 y ;
};
template <bool NA, typename RESULT_TYPE, typename U1, typename U2, typename T2>
class SugarBlock_2__PV : public Rcpp::VectorBase< Rcpp::traits::r_sexptype_traits<RESULT_TYPE>::rtype , NA, SugarBlock_2__PV<NA,RESULT_TYPE,U1,U2,T2> > {
public:
typedef RESULT_TYPE (*FunPtr)(U1,U2) ;
SugarBlock_2__PV( FunPtr ptr_, U1 u1, const T2& y_ ) :
ptr(ptr_), x(u1), y(y_){}
inline RESULT_TYPE operator[]( R_xlen_t i) const {
return ptr( x, y[i] ) ;
}
inline R_xlen_t size() const { return y.size() ; }
private:
FunPtr ptr ;
U1 x ;
const T2& y ;
};
} // sugar
} // Rcpp
#define SB2_LHT VectorBase<REALSXP,LHS_NA,LHS_T>
#define SB2_RHT VectorBase<REALSXP,RHS_NA,RHS_T>
#define SUGAR_BLOCK_2(__NAME__,__SYMBOL__) \
namespace Rcpp{ \
template <bool LHS_NA, typename LHS_T, bool RHS_NA, typename RHS_T > \
inline sugar::SugarBlock_2< (LHS_NA||RHS_NA) ,double,double,SB2_LHT,double,SB2_RHT> \
__NAME__( \
const SB2_LHT& lhs, \
const SB2_RHT& rhs \
){ \
return sugar::SugarBlock_2< (LHS_NA||RHS_NA) ,double,double,SB2_LHT,double,SB2_RHT >(\
__SYMBOL__ , lhs, rhs \
) ; \
} \
template <bool LHS_NA, typename LHS_T> \
inline sugar::SugarBlock_2__VP<LHS_NA,double,double,SB2_LHT,double> \
__NAME__( \
const SB2_LHT& lhs, \
double rhs \
){ \
return sugar::SugarBlock_2__VP<LHS_NA,double,double,SB2_LHT,double>( \
__SYMBOL__ , lhs, rhs \
) ; \
} \
template <bool RHS_NA, typename RHS_T> \
inline sugar::SugarBlock_2__PV<RHS_NA,double,double,double,SB2_RHT> \
__NAME__( \
double lhs, \
const SB2_RHT& rhs \
){ \
return sugar::SugarBlock_2__PV<RHS_NA,double,double,double,SB2_RHT >( \
__SYMBOL__ , lhs, rhs \
) ; \
} \
}
#define SUGAR_BLOCK_2_NA(__NAME__,__SYMBOL__,__NA__) \
namespace Rcpp{ \
template <bool LHS_NA, typename LHS_T, bool RHS_NA, typename RHS_T > \
inline sugar::SugarBlock_2< __NA__ ,double,double,SB2_LHT,double,SB2_RHT> \
__NAME__( \
const SB2_LHT& lhs, \
const SB2_RHT& rhs \
){ \
return sugar::SugarBlock_2< __NA__ ,double,double,SB2_LHT,double,SB2_RHT \
>( \
__SYMBOL__ , lhs, rhs \
) ; \
} \
template <bool LHS_NA, typename LHS_T> \
inline sugar::SugarBlock_2__VP<__NA__,double,double,SB2_LHT,double> \
__NAME__( \
const SB2_LHT& lhs, \
double rhs \
){ \
return sugar::SugarBlock_2__VP<__NA__,double,double,SB2_LHT,double>( \
__SYMBOL__ , lhs, rhs \
) ; \
} \
template <bool RHS_NA, typename RHS_T> \
inline sugar::SugarBlock_2__PV<__NA__,double,double,double,SB2_RHT> \
__NAME__( \
double lhs, \
const SB2_RHT& rhs \
){ \
return sugar::SugarBlock_2__PV<__NA__,double,double,double,SB2_RHT \
>( \
__SYMBOL__ , lhs, rhs \
) ; \
} \
}
#endif
|
