Eigen

Question

我正在尝试将向量整形为矩阵，但出现以下错误

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);

或者在将结果分配给 G 之前使用 MatrixXi-Map 和 .cast<double>() 结果。

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);

Eigen - 将向量重塑为矩阵

Eigen - reshape a vector to matrix

c++

reshape