c++ - boost::uBLAS 如何处理矩阵的嵌套积？

Question

我读了一篇关于使用动态规划优化矩阵嵌套乘积的文章，我想看看它是如何在 boost::uBLAS 中实现的。

我不确定我是否理解 the documentation (他们在页面的最后谈论它)，但似乎他们处理了这个问题。当您编写 R = prod(A, prod(B,C)); 时，库会计算 A x (B x C) 或 (A x B ) x C 取决于 A、B 和 C 的大小。

如何实现？图书馆如何“动括号”？在编写这样的代码行时，我认为 prod 的参数将被评估，然后该函数将被运行。

FAQ提到表达式模板的概念。有关联吗？

Answer 1

是的，涉及表达式模板(您可以在名称“表达式树”下找到它)。

基本上，prod 不返回结果，而是一个包含操作(矩阵乘法)和指向两个输入的指针的包装对象。如果使用这样的包装器作为输入调用 prod，它可以应用重新排序优化。

但是，根据我对该页面的阅读，没有应用此类优化(它说它支持括号结构)。