Eigen动态矩阵按列更新编译报错，求正确实现方案

你的编译错误主要来自两个核心问题：

• VectorXd存储的是double类型元素，order(i)返回的浮点数无法直接作为Eigen矩阵的列索引（Eigen要求索引为整数类型Eigen::Index）。
• 你使用all的方式不符合Eigen语法规范，直接调用列访问接口col()会更简洁可靠。

修复手动循环版本 #

将order(i)转换为整数索引，同时改用col()接口操作整列：

#include "eigen.h"

using namespace Eigen;

MatrixXd sort_cols(const MatrixXd& mat, const VectorXd& order)
{
    const Index rows = mat.rows();
    const Index cols = mat.cols();
    MatrixXd ordered_mat(rows, cols);

    for (Index i = 0; i < cols; ++i)
    {
        // 将double类型的order值转为Eigen标准整数索引
        const Index target_col = static_cast<Index>(order(i));
        ordered_mat.col(target_col) = mat.col(i);
    }

    return ordered_mat;
}

更高效的Eigen原生实现 #

Eigen内置了PermutationMatrix类，专门用于矩阵的行列重排，性能比手动循环更优，代码也更简洁：

#include "eigen.h"

using namespace Eigen;

MatrixXd sort_cols(const MatrixXd& mat, const VectorXd& order)
{
    // 将double类型的order向量转为整数索引向量
    VectorXi perm_vec = order.cast<Index>();
    // 根据索引向量创建排列矩阵
    PermutationMatrix<Dynamic, Dynamic> perm(perm_vec);
    // 矩阵右乘排列矩阵完成列重排
    return mat * perm;
}

注意：如果你的order向量存储的是1-based索引，需要先对每个元素减1再转换，比如perm_vec = (order.array() - 1).cast<Index>();。

关键细节 #

• 用Eigen::Index替代int作为行列数和索引类型，符合Eigen的跨平台类型规范。
• PermutationMatrix是Eigen针对排列操作优化的工具，内部实现更高效，避免手动循环的冗余操作。

内容的提问来源于stack exchange，提问作者P_B