c++ - Eigen 中没有复制乘法

Question

我将 Eigen 与大矩阵一起使用，我正在考虑优化代码的方法，重点是减少动态内存分配。

我正在尝试将两个矩阵相乘。这些矩阵时不时地发生一点变化，但它们的大小保持不变。

我希望看到乘法的输出进入一个预定义的矩阵(已经为它分配了内存)。

下面是我正在做的一个例子:

Eigen::MatrixXd   left, right,result;
// ... 
result = left * right;
// ... left and right values change a little
result = left * right;

我正在寻找这样的解决方案:

  void Multiply(
    Eigen::MatrixXd const&  left,
    Eigen::MatrixXd const&  right,
    Eigen::MatrixXd&        result);

  void Example()
  {
    Eigen::MatrixXd   left, right, result;
    // ... 
    Multiply(left, right, result);
    // ... 
    Multiply(left, right, result);
  }

主要目的是重用result 矩阵内存，因为理论上它不应该改变维度。我在考虑使用 operator*= 但我有点意识到它仍然需要一个中间矩阵来进行计算。

Answer 1

result = left * right 分配一个临时矩阵来保存乘法的结果，将乘积计算到临时矩阵中，然后将结果从临时矩阵复制到result。这是为了处理像 A = A * B 这样需要临时矩阵的语句。

如果您知道结果与产品中的矩阵不同，那么您可以编写result.noalias() = left * right。在这种情况下，Eigen 不会使用临时矩阵。

关于 Eigen 中别名的更多解释在 http://eigen.tuxfamily.org/dox/group__TopicAliasing.html