使用表达式模板时,如何创建专业化?从Wikipedia example,我可以像这样制作一个Vector sum模板类:
template <typename E1, typename E2>
class VecSum : public VecExpression<VecSum<E1, E2> > {
E1 const& _u;
E2 const& _v;
public:
VecSum(E1 const& u, E2 const& v) : _u(u), _v(v) {
assert(u.size() == v.size());
}
double operator[](size_t i) const { return _u[i] + _v[i]; }
size_t size() const { return _v.size(); }
};
template <typename E1, typename E2>
VecSum<E1,E2> const
operator+(E1 const& u, E2 const& v) {
return VecSum<E1, E2>(u, v);
}
根据维基百科的说法,如果我有一个扩展Vector
的{{1}}类和一个使用VecExpression<Vector>
运算符和循环的VecExpression
类的构造函数,允许循环合并,所以如下所述的语句只使用一个循环:
[]
我知道为什么这样可行,但我不确定如何将其扩展为标量。我希望能够向整个Vector添加标量(int,float或double),但我不知道如何执行此操作。我最好的猜测是为VecSum类创建专门化,如:
Vector a = ...;
Vector b = ...;
Vector c = ...;
Vector d = a+b+c;
但这似乎比必要的工作更多,还有另一种解决方案吗?
答案 0 :(得分:3)
只需使用SFINAE检查类型是否为算术类型,并根据需要进行专门化。
示例:
template <typename E1, typename E2, typename Enable = void > class VecSum;
template <typename E1, typename E2>
class VecSum< E1, E2,
typename std::enable_if_t<!std::is_arithmetic<E1>::value && !std::is_arithmetic<E2>::value>
> : public VecExpression<VecSum<E1, E2> >
{
E1 const& _u;
E2 const& _v;
public:
VecSum(E1 const& u, E2 const& v) : _u(u), _v(v)
{
assert(u.size() == v.size());
}
double operator[](size_t i) const { return _u[i] + _v[i]; }
size_t size() const { return _v.size(); }
};
template <typename E1, typename E2>
class VecSum < E1, E2,
typename std::enable_if_t< std::is_arithmetic<E1>::value && !std::is_arithmetic<E2>::value>
> : public VecExpression<VecSum<E1, E2> >
{
E1 const& _u;
E2 const& _v;
public:
VecSum(E1 const& u, E2 const& v) : _u(u), _v(v)
{
}
double operator[](size_t i) const { return _u + _v[i]; }
size_t size() const { return _v.size(); }
};
template <typename E1, typename E2>
class VecSum < E1, E2,
typename std::enable_if_t< !std::is_arithmetic<E1>::value && std::is_arithmetic<E2>::value>
> : public VecExpression<VecSum<E1, E2> >
{
E1 const& _u;
E2 const& _v;
public:
VecSum(E1 const& u, E2 const& v) : _u(u), _v(v)
{
}
double operator[](size_t i) const { return _u[i] + _v; }
size_t size() const { return _u.size(); }
};
int main(){
Vec v0 = { 1, 2, 3 ,4 };
Vec v1 = {10, 20,30,40 };
Vec v2 = {100,200,300,400 };
{
Vec sum = v0+v1+v2;
Vec v3(4);
for(int i=0;i<4;++i)
v3[i]=sum[i];
for(unsigned int i=0;i<v3.size();++i)
std::cout << v3[i] << std::endl;
}
std::cout << "with lhs skalar" << std::endl;
{
Vec sum = 5 + 50 + v1;
Vec v3(4);
for(int i=0;i<4;++i)
v3[i]=sum[i];
for(unsigned int i=0;i<v3.size();++i)
std::cout << v3[i] << std::endl;
}
std::cout << "with rhs skalar" << std::endl;
{
Vec sum = v1 + 5 + 50 ;
Vec v3(4);
for(int i=0;i<4;++i)
v3[i]=sum[i];
for(unsigned int i=0;i<v3.size();++i)
std::cout << v3[i] << std::endl;
}
}
答案 1 :(得分:1)
另一种选择:
class AnySize { };
bool operator==(AnySize, AnySize) { return true; }
bool operator==(size_t, AnySize) { return true; }
bool operator==(AnySize, size_t) { return true; }
template <typename T>
class Scalar {
T _value;
public:
explicit Scalar(T value) : _value(value) {}
T operator[](int) { return _value; }
AnySize size() { return {}; }
}
您还需要更改VecExpression::size
,以提供合适的尺寸,而不是始终在左侧
AnySize size_for(AnySize, AnySize) { return {}; }
size_t size_for(size_t lhs, AnySize) { return lhs; }
size_t size_for(AnySize, size_t rhs) { return rhs; }
size_t size_for(size_t lhs, size_t rhs) { return lhs; }
auto VecExpression<E1,E2>::size() { return size_for(_u.size(), _v.size()); }