特征：：SelfAdjointView：：rankUpdate更新慢于+= w*w.transpose()

Question

本征的测试速度：：SelfAdjointView：：用本征：：Matrix4d更新

与朴素的+= w*w.transpose()相比，它慢了2倍。

1. 我做错什么了？
2. 我能加快计算速度吗？

Answer 1

对于小的固定大小的表达式，您不能用SelfAdjointView::rankUpdate保存任何东西，相反，它会增加开销，因为它需要确保只修改一半的元素。在你的例子中一个简单的

A.noalias() += w*w.adjoint();

应该给出接近最优的代码(添加.noalias()可以避免将副本复制到临时代码中)。