eigen：如何只计算矩阵内积的下/上部分

Question

我需要使用Eigen 计算像“A'*A”这样的公式，其中 A 是一个 m x n 矩阵。这样做的直观方法是，

result = A.transpose()*A;

但是由于结果是对称的，是否可以只计算结果的下半部分或上半部分？

Answer 1

是的，使用 selfadjointView 和 rankUpdate:

result.setZero();
result.selfadjointView<Lower>().rankUpdate(A.transpose());

这只会更新result的下部。