0

我在 Eigen 文档中关注此页面,试图了解 Eigen 参数的使用

https://eigen.tuxfamily.org/dox/TopicFunctionTakingEigenTypes.html

以下代码按预期工作

#include <iostream>
#include <Eigen/Dense>

// a - nonwriteable, b - writeable
void eigen_reference_class(
    const Eigen::Ref<const Eigen::Array<double, 
    Eigen::Dynamic, 3> >& a, Eigen::Ref<Eigen::Array<double, Eigen::Dynamic, 3> > b) {
    b = 2 * a;
}

void eigen_reference(const Eigen::Array<double, 1, 3>& a, Eigen::Array<double, 1, 3>& b) {
    b = 2*a;
}

template<typename Derived> 
void eigen_template(const Eigen::PlainObjectBase<Derived>& a, Eigen::PlainObjectBase<Derived>& b) {
    b = 2*a;
}

int main()
{
    Eigen::Array<double, 1, 3> a, b, c, d;
    a << 1, 2, 3;
    eigen_reference_class(a, b);
    eigen_reference(a, c);
    eigen_template(a,  d);
    std::cout << "a : \n" << a << std::endl;
    std::cout << "b: \n" << b << std::endl;
    std::cout << "c: \n" << c << std::endl;
    std::cout << "d: \n" << d << std::endl;

}

但是,如果数组的初始声明更改为

Eigen::Array<double, Eigen::Dynamic, 3> a, b, c, d;

然后程序将无法通过以下方式编译:

error: invalid initialization of non-const reference of type ‘Eigen::Array<double, 1, 3>&’ from an rvalue of type ‘Eigen::Array<double, 1, 3>’       
     eigen_reference(a, c);

或者即使只a保留定义,也会因分段错误而失败。

4

1 回答 1

1

我喜欢看看int在这种情况下会发生什么,因为它可以帮助我更容易地理解在这种情况下发生的更复杂的事情。

想象一下,你有一个struct可以int隐式转换为 an 的int东西,这样你就可以用它来代替 an int,就像这样:

struct Integer {
    int m_int;
    operator int() { return m_int; }
};

现在定义两个函数,一个取一个int const &,一个取一个int &

void const_ref(int const &) {}
void ref(int &) {}

如果我们在这些函数中使用常规整数,则不会有任何意外。

int j = 3;
const_ref(j);
ref(j);

但是如果我们Integer改用它,它就不再编译了。

Integer i{2};
const_ref(i);
ref(i);

错误类似于

error: no matching function for call to 'ref'
    ref(i);
    ^~~
note: candidate function not viable: no known conversion from 'Integer' to 'int &' for 1st argument
void ref(int &) {}
     ^
1 error generated.

现在的问题是,为什么const_ref有效但ref无效?

在调用中const_ref发生的事情是调用转换运算符 to ,它返回一个我们可以绑定常量引用int的临时值。int我们不能将可变引用绑定到这个临时对象,因为这会导致奇怪的效果。想象一下,该函数ref修改了参数,但在这种情况下,它会修改一个临时的,并且不会将任何内容写入Integer您传入的原始参数......一个保证的错误。


在你的代码中是一样的。 Eigen::Array<double, Eigen::Dynamic, 3>可以隐式转换为,Eigen::Array<double, 1, 3>但不能将可变引用绑定到临时引用。

于 2018-03-14T05:21:08.440 回答