传递指针中的 Eigen::Ref

Question

类似于提问Pointer vs Reference difference when passing Eigen objects as arguments

假设我们有 foo1 和矩阵 mat2by2：

void foo1(MatrixXd& container){
    //...container matrix is modified here
}

和

Matrix33d mat2by2;
mat2by2 << 1,2,
           3,4;

在阅读http://eigen.tuxfamily.org/dox/TopicFunctionTakingEigenTypes.html 之后，似乎 foo1 的更好替代方案可能是：

void foo2(Ref<MatrixXd> container){
    //...container matrix is modified here
}

如果 foo2 的参数作为引用传递， 使用 Eigen::Ref 类的指针传递等效于什么？？

void foo(Eigen::MatrixXd* container){
    //...container matrix is modified here
}

Answer 1

我认为基本思想是不直接使用指针或引用。 Eigen 使用模板表达式来表示计算。这意味着类型会根据用于计算矩阵的表达式而变化，并且表达式可能会在未计算的情况下进行。

如有必要，Ref 会将模板表达式评估为与您请求作为参数传递的内存布局匹配的临时对象。如果参数的内存布局与参数所需的内存布局相匹配，则 Ref 将充当透明引用。

直接从文档中借用：您的输入参数应该声明为常量，而非 const 参数可以用作输出参数。

void cov(const Ref<const MatrixXf> x, Ref<MatrixXf> C)
{
   ...
   C = ...; // Your return value here
}

如果你读写一个矩阵，参数显然也应该是非常量的。

对于可选参数，您可以使用指向 Ref 的指针。

编辑：文档确实注意到您可以直接使用常量引用来传递参数。这只是因为编译器很乐意将临时对象转换为常量引用。它不适用于指针。