ublas矩阵表达教程/示例

Question

我正在尝试实现某些矩阵操作,但我迷失在ublas库的内部.是否有一个资源,如教程或如何实现新的ublas矩阵表达式的例子？

谢谢

Answer 1

我的建议是仅对新函数进行模板化，这样您就不必担心矩阵表达式或 ublas 内部结构。例如，如果您想编写自己的反函数，请将其编写为模板：

template<typename MATRIX_IN, typename MATRIX_OUT>
bool invert_matrix(const MATRIX& A, MATRIX_OUT& B)
{
    // now use A(i,j), B(i,j), etc
    ...
}