我将 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*=
但我有点意识到它仍然需要一个中间矩阵来进行计算。
最佳答案
result = left * right
分配一个临时矩阵来保存乘法的结果,将乘积计算到临时矩阵中,然后将结果从临时矩阵复制到result
。这是为了处理像 A = A * B
这样需要临时矩阵的语句。
如果您知道结果与产品中的矩阵不同,那么您可以编写result.noalias() = left * right
。在这种情况下,Eigen 不会使用临时矩阵。
关于 Eigen 中别名的更多解释在 http://eigen.tuxfamily.org/dox/group__TopicAliasing.html
关于c++ - Eigen 中没有复制乘法,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/24259998/