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“矩阵”不是连续的。 你可以做的是将数据复制到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);