表达式必须具有常数值本征矩阵

Question

我一直在尝试通过函数来​​给出运行时的矩阵大小，例如

Eigen::MatrixXd FileReader::load_csv(const std::string & path, const int rows_csv, const int columns_csv) {
    std::ifstream indata;
    indata.open(path);
    std::string line;
    std::vector<double> values;
    int rows = 0;
    while (getline(indata, line)) {
        std::stringstream lineStream(line);
        std::string cell;
        while (std::getline(lineStream, cell, ',')) {
            values.push_back(std::stod(cell));
        }
        ++rows;
    }
    std::cout << "loading";
    return Eigen::Map<const Eigen::Matrix<double, columns_csv, rows_csv, Eigen::RowMajor>>(values.data(), rows, values.size() / rows);
}

（以上代码来自https://stackoverflow.com/a/39146048/3782963 ），我无法将矩阵大小的常量发送给函数，我得到Expression must have constant value错误。 有什么办法可以做这样的事情：

Eigen::MatrixXd mat = load_csv("some_path", 20, 30);

任何想法？

Answer 1

这行是你的问题：

Eigen::Map<const Eigen::Matrix<double, columns_csv, rows_csv, Eigen::RowMajor>>
                                          (values.data(), rows, values.size() / rows);

的

Eigen::Map<const Eigen::Matrix<double, columns_csv, rows_csv, Eigen::RowMajor>>
                                           ^^^         ^^^

告诉Map ，它应该期望尺寸不变的矩阵，并且

(values.data(), rows, values.size() / rows);
                ^^^              ^^^

是一个动态大小的矩阵。 如果由于某种原因，您想要忽略文件中矩阵的形状并在程序中定义它，请尝试：

Eigen::Map<Eigen::Matrix<double, Eigen::Dynamic, Eigen::Dynamic, Eigen::RowMajor>>
                                                 (values.data(), rows_csv, columns_csv);

Answer 2

阅读有关矩阵大小分配的文档也有帮助，该文档指出：

固定尺寸与动态尺寸：

固定大小意味着已知的行数和列数是编译时的。 在这种情况下，Eigen将系数数组分配为固定大小的数组，作为类成员。 对于非常小的矩阵（通常最大为4x4，有时最大为16x16），这是有意义的。 即使较大的矩阵在编译时就知道其大小，也应将其声明为动态大小。 动态大小意味着在编译时不一定知道行数或列数。 在这种情况下，它们是运行时变量，并且系数数组在堆上动态分配。 请注意，稠密矩阵（无论是固定大小还是动态大小）在std :: map的意义上不会动态扩展。 如果您想要这种行为，请参阅稀疏模块。