微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

调用`segment`之后未分配的`Eigen :: VectorXd`对象会怎样?

如何解决调用`segment`之后未分配的`Eigen :: VectorXd`对象会怎样?

我有一个Foo,它存储匹配尺寸的本征MatrixXdVectorXd,并提供了执行矩阵矢量乘法的简单函数.mult()我有一个函数调用mult()获取返回向量的一部分,并进行一些进一步的计算。这是代码的简化版本:

#include <Eigen/Dense>

class Foo {
  Eigen::MatrixXd X;
  Eigen::VectorXd v;
public:
  Foo() {
    X.resize(2,2);
    v.resize(2);
    X << 1,4,2,4;
    v << -1,1;
  }
  Eigen::VectorXd mult() {
    return X * v;
  }
};

Eigen::VectorXd bar() {
  Foo f;
  Eigen::VectorXd val = f.mult().segment(1,1);
// Use val to do some computation.
  Eigen::VectorXd ret = val.array() + 2;
  return ret;
}

代码审查期间,有人提出一条评论,认为这可能导致对存储f.mult()结果的未命名/临时变量的引用悬挂。对于未分配名称的临时变量的生存期,人们提出了一个普遍的问题。

如果我正确理解RVO,则会在堆栈上创建一个临时对象,其地址将传递到f.mult(),然后将X * v的结果存储到传递的内存地址中。当bar()终止时,此临时对象将超出范围。对于临时和未分配变量的范围,不应有任何悬而未决的参考或关注。但是我不太确定。

Eigen是一个优雅而复杂的库,我希望确保尽可能安全地使用它。预先感谢!

解决方法

您正在做的是节省(但效率低下)。

f.mult()按值返回Eigen::VectorXd对象,该对象存在于语句的其余部分(“直到下一个;”)。因此,只要将结果存储到segment()(而不是存储到VectorXd对象中),就可以对其进行调用。{p>

要获得矩阵向量乘积的一个分段,将矩阵的相应子块乘以它会更有效,即:

auto

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。