如何解决调用`segment`之后未分配的`Eigen :: VectorXd`对象会怎样?
我有一个类Foo
,它存储匹配尺寸的本征MatrixXd
和VectorXd
,并提供了执行矩阵矢量乘法的简单函数.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 举报,一经查实,本站将立刻删除。