Eigen 稀疏矩阵LU分解解方程组

下载Eigen、配置目录

这里用LU分解举例

 #include <Eigen/Dense>
   #include <Eigen/Sparse>
   #include <Eigen/SparseLU>
   using namespace Eigen;
   //处理矩阵部分
   vector<Triplet<float>> triple;//用三元组赋值
   triple.push_back(Triplet<float>(i, j,A[i][j]));//ij处的值是A[i][j]
   SparseMatrix<float> A(n, n);//定义n*n的系数矩阵
  A.setFromTriplets(triple.begin(), triple.end());//三元数给矩阵赋值
  SparseLU<SparseMatrix<float>> solver
   solver.compute(A);//对A进行预分解
  if (solver.info()!=Success)
    {
        cout << "Compute Matrix is error" << endl;
        return;
    }
   //右端项
  VectorXf b(n);//n行的方程组
  //求解
  x = solver.solve(b);