c++ - Eigen - 将 vector reshape 为矩阵

Question

我正在尝试将 vector reshape 为矩阵，但出现以下错误

unsigned int Nx     = 8;
unsigned int Ny     = 7;
Eigen::VectorXi G_temp = Eigen::VectorXi::LinSpaced((Nx + 2) * (Ny + 2),0,(Nx + 2) * (Ny + 2)-1);
Eigen::MatrixXd G = Eigen::Map<Eigen::MatrixXd>(G_temp.data(),Nx+2, Ny+2); // error: no matching constructor for initialization of 'Eigen::Map<Eigen::MatrixXd>'

我按照写的here进行操作，但我不明白我做错的方式。

Answer 1

Eigen 中没有从整数值表达式到双值表达式的隐式转换。要么只使用 VectorXd对于 G_temp (和 LinSpaced 表达式):

Eigen::VectorXd G_temp = Eigen::VectorXd::LinSpaced((Nx + 2) * (Ny + 2),0,(Nx + 2) * (Ny + 2)-1);

或者使用 MatrixXi - map 和.cast<double>()将其分配给 G 之前的结果.

Eigen::MatrixXd G = Eigen::Map<Eigen::MatrixXi>(G_temp.data(),Nx+2, Ny+2).cast<double>();

为了避免任何临时的情况，您还可以分配 MatrixXd并直接在里面分配适当的值:

Eigen::MatrixXd G(Nx+2, Ny+2); // allocate matrix
// set values in-place:
Eigen::VectorXd::Map(G.data(), (Nx + 2) * (Ny + 2)).setLinSpaced(0,(Nx + 2) * (Ny + 2)-1);

或者使用 master/3.4 分支:

G.reshaped().setLinSpaced(0,(Nx + 2) * (Ny + 2)-1);