Kết quả sai khi sử dụng tự động với Eigen

Question

Tôi nhận được kết quả khác nhau bằng tự động và sử dụng Vector khi tổng hợp hai vectơ.

Mã của tôi:

#include "stdafx.h" 
#include <iostream> 
#include "D:\externals\eigen_3_1_2\include\Eigen\Geometry" 

typedef Eigen::Matrix<double, 3, 1>  Vector3; 

void foo(const Vector3& Ha, volatile int j) 
{ 
    const auto resAuto = Ha + Vector3(0.,0.,j * 2.567); 
    const Vector3 resVector3 = Ha + Vector3(0.,0.,j * 2.567); 

    std::cout << "resAuto = " << resAuto <<std::endl; 
    std::cout << "resVector3 = " << resVector3 <<std::endl; 
} 

int main(int argc, _TCHAR* argv[]) 
{ 
    Vector3 Ha(-24.9536,-29.3876,65.801); 
    Vector3 z(0.,0.,2.567); 

    int j = 7; 

    foo(Ha,j); 
    return 0; 
} 

Kết quả:

resAuto = -24,9536, -29.3876,65.801

resVector3 = -24,9536, -29.3876,83.77

Press bất kỳ chìa khóa nào để tiếp tục. . .

Tôi hiểu rằng Eigen tối ưu hóa nội bộ tạo ra các kết quả khác nhau. Nhưng nó trông giống như một lỗi trong Eigen và C++ 11.

Answer 1

Từ khóa auto yêu cầu trình biên dịch "đoán" đối tượng tốt nhất dựa trên bên phải của =. Bạn có thể kiểm tra kết quả bằng cách thêm

std::cout << typeid(resAuto).name() <<std::endl; 
std::cout << typeid(resVector3).name() <<std::endl; 

để foo (đừng quên để bao gồm <typeinfo>).

Trong trường hợp này, sau khi xây dựng tạm thời Vector3, phương thức operator+ được gọi, tạo đối tượng CwiseBinaryOp. Đối tượng này là một phần của Eigens lazy evaluation (có thể tăng hiệu suất). Nếu bạn muốn buộc đánh giá háo hức (và do đó hãy xác định loại), bạn có thể sử dụng

const auto resAuto = (Ha + Vector3(0.,0.,j * 2.567)).eval(); 

thay vì dòng của bạn trong foo.

Một vài lưu ý phụ:

• Vector3 là giống hệt với lớp Vector3d quy định tại Eigen
• Bạn có thể sử dụng #include <Eigen/Core> thay vì #include <Eigen/Geometry> để bao gồm hầu hết các tiêu đề Eigen, cộng với những điều nào đó được định nghĩa đó mà nên là.