我正在尝试实现某些矩阵运算,但我迷失在 ublas 库的内部。是否有关于如何实现新的 ublas 矩阵表达式的教程或示例等资源?
谢谢
不知道它是否有帮助,但这里有一个关于扩展 uBlas 的 wiki页面。那个表达模板的东西真的让我大吃一惊。:)
我的建议是对你的新函数进行模板化,这样你就不必担心矩阵表达式或 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
...
}
这里也有一些例子,如果它们有帮助的话。