Question

我有兴趣创建一个函数Derivative，它返回一个函数，该函数是某些函数的衍生物，在某些时候传递给它。但是，我希望能够对此进行专门化，以便对于特定的功能，我可以返回分析解决方案。

所以，我正在寻找这样的东西：

auto Derivate(alias Function)(x)
{ return (Function(x+h) - Function(x-h))/(2h);}

auto Derivate(BSpline!(k)(x))(x)
{ return k * BSpline!(k-1)(x) + x * BSpline!(k-1)(x); }

但是，我目前以这种方式定义了BSpline：

pure Real BSpline(int k : 0, Real)(scope Real x, scope const(Real)[] t)
{
    if (t[0] <= x && x < t[k+1])
        return 1;
    else
        return 0;
}

pure Real BSpline(int k, Real)(scope Real x, scope const(Real)[] t)
{
    if (t[0] <= x && x < t[k+1])
    {
        Real a = (x - t[0]) / (t[k] - t[0]);
        Real b = (t[k+1] - x) / (t[k+1] - t[1]);
        Real c = BSpline!(k-1,Real)(x, t[0..k+1]);
        Real d = BSpline!(k-1,Real)(x, t[1..k+2]);
        Real rv = (c?c*a:c) + (d?d*b:d);
        return rv;
    }
    else
        return 0;
}

因此BSpline上的类型签名将是Real函数（Real，Real），这与任何其他类型的函数都不可区分。解决这个问题的方法是创建一个定义了opCall的“BSpline”类吗？或者我可以使用某种typedef来识别这个函数吗？

谢谢！

Answer 1

要专门化模板，您必须使用:表示法：

auto foo(alias F_, X_)(X_ x) {
    /* code here ... */
}

auto foo(alias F_ : BSpline, X_)(X_ x) {
    /* specialized version here */
}

D中的函数类型

1 个答案: