c++ - Eigen SuiteSparse (Eigen::SPQR) 不匹配 num 行

Question

我正在使用 Eigen::SPQR 来解决稀疏矩阵系统这是会给出运行时错误的最小示例

# include <Eigen/Core>
# include <Eigen/Sparse>
# include <Eigen/SparseLU>
# include <Eigen/SparseQR>
///  test eigen SPQR
    Eigen::SparseMatrix < double > A(500, 400);
    Eigen::VectorXd b(500);
    std::vector < Eigen::Triplet < double > > triplet;
    for( int i = 0; i < 500; i++ )
    {
        int r = rand() % 500;
        int c = rand() % 400;
        triplet.push_back( Eigen::Triplet < double > ( r, c, static_cast<double>(rand())/(RAND_MAX+1.0) ) );

        b(i) = static_cast<double>(rand())/(RAND_MAX+1.0);
    }
    A.setFromTriplets(triplet.begin(), triplet.end());
    Eigen::SPQR< Eigen::SparseMatrix < double > > spqrsolver(A);

    std::cout << A.rows() << " " << A.cols() << " "
              << b.rows() << " " << b.cols() << " "
              << spqrsolver.rows() << " " << spqrsolver.cols() << std::endl;
    Eigen::VectorXd X = spqrsolver.solve(b);

输出为 500 400 500 1 400 400。错误消息为“断言 `this->rows()==B.rows() && “SPQR::solve(): 右侧矩阵 B 的行数无效” ' 失败的。”