c++ - 使用模板允许使用任意 LinAlg 库

Question

我有一个执行一些线性代数计算的函数，我目前使用 the Eigen library .

函数看起来像这样:

void Predict(const MatrixXd& F,
        VectorXd& x,
        VectorXd& u,
        MatrixXd& P,
        MatrixXd& Q)
    {
        x = F*x + u;
        P = F*P*F.transpose() + Q;
    }

现在我想使用模板概括这个函数，这样客户端就可以使用任意的 LinAlg 库。我的第一种方法如下所示:

template<typename Matrix, typename Vector>
    void Predict(const Matrix& F,
        Vector& x,
        Vector& u,
        Matrix& P,
        Matrix& Q)
    {
        x = F*x + u;
        P = F*P*F.transpose() + Q;
    }

然而，明显的问题是，该函数依赖于 transpose() 方法。尽管底层功能应由任何 LinAlg 实现实现，但相应方法的签名可能而且很可能会有所不同。一种解决方案是将转置矩阵作为参数传递，这在这种情况下有效，但在无法事先计算的情况下无效。那么像这样抽象出约束的最佳方法是什么？

Answer 1

提供包装模板函数:

template<typename Matrix>
Matrix transpose(const Matrix& m
{
    return m.transpose();
}

并使用它:

template<typename Matrix, typename Vector>
void Predict(const Matrix& F,
    Vector& x,
    Vector& u,
    Matrix& P,
    Matrix& Q)
{
    x = F*x + u;
    P = F*P*transpose(F) + Q;
}

如果您曾经有一个具有一些不同语法的库，请专门化或重载模板:

SpecialMatrixImplementation transpose(SpecialMatrixImplementation const& m)
{
    m.calculateTransposition();
}

这与 e 的想法相同。 G。 std::begin()和 std::end()函数...