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

在Eigen中没有复制乘法

1 个答案: