Eigen và std :: vector

Question

Tôi có một ma trận, được đưa ra như:

std::vector<std::vector<std::complex<double>>> A; 

Và tôi muốn để lập bản đồ đó đến Eigen thư viện đại số tuyến tính như thế này:

Eigen::Map<Eigen::MatrixXcd, Eigen::RowMajor> mat(A.data(),51,51); 

Nhưng mã không thành công với

error: no matching function for call to   
‘Eigen::Map<Eigen::Matrix<std::complex<double>, -1, -1>, 1>:: 

Có cách nào để chuyển đổi vectơ của vectơ sao cho Eigen có thể sử dụng nó?

Answer 1

Eigen sử dụng bộ nhớ liền kề, cũng như std::vector. Tuy nhiên, bên ngoài std::vector chứa một tập liền kề std::vector<std::complex<double> >, mỗi trỏ đến một tập các số phức khác (và có thể có độ dài khác nhau). Vì vậy, std "ma trận" không tiếp giáp. Những gì bạn có thể làm là sao chép dữ liệu vào ma trận Eigen, có nhiều cách để thực hiện điều đó. Cách đơn giản nhất sẽ là để lặp qua i và j, với một lựa chọn tốt hơn là một cái gì đó giống như

Eigen::MatrixXcd mat(rows, cols); 
for(int i = 0; i < cols; i++) 
    mat.col(i) = Eigen::Map<Eigen::VectorXcd> (A[i].data(), rows);