将 v.cwiseAbs() 传递给接受 Ref<VectorXd> 的函数

Question

以下似乎不起作用。

void Foo(Ref<VectorXd> v) {
  // modifies v
}

Eigen::VectorXd v;
Foo(v.cwiseAbs());

出现以下错误信息

错误：无法用 Derived = Eigen::ArrayWrapper, -1, 1, true>, -1, 1, false> 转换 'Eigen::ArrayBase::cwiseAbs() const >; Eigen::ArrayBase::CwiseAbsReturnType = Eigen::CwiseUnaryOp, const Eigen::ArrayWrapper, -1, 1, true>, -1, 1, false> > >; typename Eigen::internal::traits::Scalar = double' from 'const CwiseAbsReturnType {aka const Eigen::CwiseUnaryOp, const Eigen::ArrayWrapper, -1, 1, true>, -1, 1, false> > >} ' 到 'Eigen::Ref >'

任何建议为什么以及如何解决？

Answer 1

这不起作用，因为需要两个额外的const 限定符才能匹配.cwiseAbs() 返回的表达式。这是有道理的，因为不可能通过以引用形式接受此参数的函数来修改v.cwiseAbs() 的结果。以下代码编译：

void Foo(const Ref<const VectorXd>& v) {
  std::cout << v << std::endl;
}

int main() {
  Eigen::VectorXd v(3);
  v << 1,-2,3;
  Foo(v.cwiseAbs());
}

但是，通过此修改，不允许在Foo() 内修改v。

最简单的解决方案可能是删除Ref<> 并使用

 Foo(VectorXd v) {...}

相反。这会生成一个本地副本，但这在性能方面应该不是问题，此外，如果要在Foo() 中修改v，无论如何都无法避免任何副本。如果Ref 保留在Foo 的签名中，则可以复制v.cwiseAbs() 并使用该副本调用Foo()：

void Foo(Ref<VectorXd> v) {...}
...
Eigen::VectorXd w = v.cwiseAbs();
Foo(w);

正如@chtz 所建议的，C++11 提供的替代方案是

void Foo(VectorXd&& v) {...}

虽然这允许在Foo() 中修改v，但使用它可能有些危险，因为在Foo() 中只更改了一个临时值。在Foo() 中对v 的修改不会改变main() 中的v。