SparseCholesky.hpp

Go to the documentation of this file.

#pragma once

#include <memory>
#include <Eigen/Sparse>
#include <solver/interface/SolverSetting.hpp>

namespace linalg {

  class SparseCholesky
  {
    public:
      SparseCholesky(){}
      ~SparseCholesky(){}

      void analyzePattern(const Eigen::SparseMatrix<double>& mat, const solver::SolverSetting& stgs);
      int  factorize(const Eigen::SparseMatrix<double>& mat, const Eigen::Ref<const Eigen::VectorXd>& sign);
      void solve(const Eigen::Ref<const Eigen::VectorXd>& b, double* x);
      Eigen::VectorXd& solve(const Eigen::VectorXd& b);

    private:
      inline const solver::SolverSetting& getSetting() const { return *setting_; }

      int n_;
      double eps_, delta_;
      Eigen::VectorXd D_, Y_, X_;
      Eigen::SparseMatrix<double> L_;
      Eigen::VectorXi Parent_, Pattern_, Flag_, Lnnz_;
      std::shared_ptr<const solver::SolverSetting> setting_;

  };

}