这个问题有很多很好的答案。还应该提到模板支持开放设计。在面向对象编程语言的当前状态下,在处理此类问题时必须使用访问者模式,而真正的OOP应该支持多个动态绑定。参见c++的开放多方法,P. Pirkelbauer等。非常有趣的阅读。
模板的另一个有趣之处在于,它们也可以用于运行时多态性。例如
template<class Value,class T>
Value euler_fwd(size_t N,double t_0,double t_end,Value y_0,const T& func)
{
auto dt=(t_end-t_0)/N;
for(size_t k=0;k<N;++k)
{y_0+=func(t_0 + k*dt,y_0)*dt;}
return y_0;
}
注意,如果Value是某种类型的向量(不是std::vector,应该称为std::dynamic_array以避免混淆),则此函数也可以工作。
如果func很小,这个函数将从内联中获得很多。示例使用
auto result=euler_fwd(10000,0.0,1.0,1.0,[](double x,double y)
{return y;});
在这种情况下,你应该知道确切的答案(2.718…),但是很容易构造一个没有初等解的简单ODE(提示:在y中使用多项式)。
现在,您在func中有一个大表达式,并且在许多地方使用ODE求解器,因此您的可执行文件到处都受到模板实例化的污染。怎么办呢?首先要注意的是,常规函数指针可以工作。然后,您希望添加curry,以便编写接口和显式实例化
class OdeFunction
{
public:
virtual double operator()(double t,double y) const=0;
};
template
double euler_fwd(size_t N,double t_0,double t_end,double y_0,const OdeFunction& func);
但是上面的实例化只适用于double,为什么不把接口写成模板呢:
template<class Value=double>
class OdeFunction
{
public:
virtual Value operator()(double t,const Value& y) const=0;
};
并专门化一些常见的值类型:
template double euler_fwd(size_t N,double t_0,double t_end,double y_0,const OdeFunction<double>& func);
template vec4_t<double> euler_fwd(size_t N,double t_0,double t_end,vec4_t<double> y_0,const OdeFunction< vec4_t<double> >& func); // (Native AVX vector with four components)
template vec8_t<float> euler_fwd(size_t N,double t_0,double t_end,vec8_t<float> y_0,const OdeFunction< vec8_t<float> >& func); // (Native AVX vector with 8 components)
template Vector<double> euler_fwd(size_t N,double t_0,double t_end,Vector<double> y_0,const OdeFunction< Vector<double> >& func); // (A N-dimensional real vector, *not* `std::vector`, see above)
如果函数首先是围绕接口设计的,那么您将被迫继承ABC。现在您有了这个选项,还有函数指针、lambda或任何其他函数对象。这里的关键是我们必须有operator()(),我们必须能够在它的返回类型上使用一些算术运算符。因此,在这种情况下,如果c++没有操作符重载,模板机制就会中断。