正定稀疏矩阵求系数解X

#include <Eigen/Sparse>
#include <Eigen/IterativeLinearSolvers>
#include <iostream>

int main() {
    int num = 100; // 假设的尺寸
    Eigen::SparseMatrix<float> A(num, num);
    Eigen::VectorXf B(num);

    // 填充 A 和 B
    // ...

    // 创建共轭梯度求解器实例
    Eigen::ConjugateGradient<Eigen::SparseMatrix<float>, Eigen::Lower|Eigen::Upper> cg;

    // 计算 A
    cg.compute(A);

    // 检查是否成功分解
    if (cg.info() != Eigen::Success) {
        std::cerr << "Matrix decomposition failed" << std::endl;
        return -1;
    }

    // 求解线性系统
    Eigen::VectorXf X = cg.solve(B);

    // 检查求解是否成功
    if (cg.info() != Eigen::Success) {
        std::cerr << "Solving failed" << std::endl;
        return -1;
    }

    // X 现在包含解向量
    // ...

    return 0;
}