Eigen::Isometry3d中rotation()与linear()函数性能差异问询

核心原因：返回类型与表达式模板的开销差异 #

在Eigen 3.3.4中，Isometry3d的rotation()和linear()虽然语义上都指向变换的线性（旋转）部分，但返回类型完全不同，这直接导致了性能差异：

• linear()直接返回Matrix3d&（const版本为const Matrix3d&），是对变换内部矩阵块的直接引用，没有额外包装。
• rotation()返回的是RotationMatrix3d（Eigen::RotationBase的派生类），这是一个表达式模板对象，它对linear()返回的矩阵做了一层包装，用于保留“旋转矩阵”的语义（比如正交性）。

不同优化级别下的表现 #

1. -O0 无优化场景
   无优化时，编译器不会消除表达式模板带来的临时对象构造、函数调用开销。rotation()每次调用都会创建一个RotationMatrix3d临时对象，两次rotation()调用加上乘法操作，会触发更多的函数跳转和临时内存操作，这就是测试中耗时差两个数量级的原因。而linear()直接操作矩阵引用，乘法逻辑更直接，开销极小。

2. -O3 最高优化场景
   高优化级别下，编译器会尝试折叠表达式模板的额外开销，但Eigen 3.3.4版本对RotationMatrix3d乘法的优化可能不如直接的Matrix3d乘法彻底。部分额外的语义检查或包装逻辑可能无法完全被优化掉，导致rotation()仍略慢于linear()。升级到Eigen 3.4及以上版本后，这类表达式模板的性能通常会有明显改善。

代码层面的佐证 #

查看Eigen 3.3.4的源码实现：

• linear()的实现是直接返回内部矩阵的块引用：

  EIGEN_DEVICE_FUNC LinearMatrixType& linear() { return m_matrix.block<Dim,Dim>(0,0); }
  

• rotation()则是通过包装linear()返回的矩阵构造一个RotationMatrixType临时对象：

  EIGEN_DEVICE_FUNC RotationMatrixType rotation() { return RotationMatrixType(linear()); }
  

这种包装在无优化时会产生显著的额外开销，即使优化后也可能残留细微差异。

总结 #

如果只需要执行矩阵乘法这类纯数值操作，直接使用linear()可以避免表达式模板带来的额外开销；如果需要利用旋转矩阵的特殊语义（比如正交化、转换为四元数等），再使用rotation()会更合适。

内容的提问来源于stack exchange，提问作者DarioZ