c++ - 如何将 ProductReturnType 转换为矩阵？

Question

我正在使用 C++ 线性代数库 eigen。我尝试将 2 个矩阵相乘:

static void do_stuff_with_matrix(Eigen::MatrixXf& mat) {
  return;
}

Eigen::MatrixXf a(3, 4);
Eigen::MatrixXf b(4, 5);

Eigen::MatrixXf c = a * b;
do_stuff_with_matrix(c);

不幸的是，我收到一个编译器错误，指出 ProductReturnType(c 是)无法转换为 Eigen::MatrixXf&。我如何执行此转换？

Answer 1

Eigen 使用惰性求值来防止不必要的临时值和其他东西。结果 c 本质上是一个 ProductReturnType , 矩阵乘积的优化结构:

template<typename Lhs, typename Rhs, int ProductType>
class Eigen::ProductReturnType< Lhs, Rhs, ProductType >

Helper class to get the correct and optimized returned type of operator*. [see also 2]

为了从 A * B 形式的表达式创建一个真正的矩阵，你需要直接计算它:

Eigen::MatrixXf c = (a * b).eval();
do_stuff_with_matrix(c);

参见 this page有关 Eigen 的惰性求值和别名的更多信息。