Eigen和std :: vector

Question

我有一个矩阵，其名称为：

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

我想将它映射到Eigen线性代数库，如下所示：

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

但代码失败了

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

是否有转换矢量的矢量以便Eigen可以使用它？

Answer 1

Eigen使用连续内存，std::vector也是如此。但是，外部std::vector包含一组连续的std::vector<std::complex<double> >，每个都指向一组不同的复数（并且可以是不同的长度）。因此，std＆＃34;矩阵＆＃34;不是连续的。你可以做的是将数据复制到Eigen矩阵，有多种方法可以做到这一点。最简单的是循环遍历i和j，更好的选择是

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