Eigen_LUSolver.hpp

Go to the documentation of this file.

// _LICENSE_HEADER_
//
// Copyright (C) 2019 - 2025.
// Terms register on the GPL-3.0 license.
//
// This file can be redistributed and/or modified under the license terms.
//
// See top level LICENSE file for more details.
//
// This file can be used citing references in CITATION.cff file.
 
#ifndef __EIGENLUSOLVER_H
#define __EIGENLUSOLVER_H
 
#include "Eigen/Eigen"
#include "ILinearSolver.hpp"
 
namespace Gedim
{
template <typename Eigen_ArrayType = Eigen::VectorXd, typename Eigen_SparseArrayType = Eigen::SparseMatrix<double>>
class Eigen_LUSolver final : public ILinearSolver
{
  private:
    Eigen::SparseLU<Eigen_SparseArrayType> linearSolver; 
    const IArray *_rightHandSide;                        
    IArray *_solution;                                   
 
  public:
    Eigen_LUSolver()
    {
        _rightHandSide = nullptr;
        _solution = nullptr;
    }
    ~Eigen_LUSolver()
    {
        _rightHandSide = nullptr;
        _solution = nullptr;
    }
 
    void Initialize(const ISparseArray &matrix, const IArray &rightHandSide, IArray &solution, const Configuration &config = Configuration());
 
    void Initialize(const IArray &rightHandSide, IArray &solution, const Configuration &config = Configuration());
 
    SolutionInfo Solve() const;
 
    void Initialize(const ISparseArray &matrix, const Configuration &config = Configuration());
 
    SolutionInfo Solve(const IArray &rightHandSide, IArray &solution) const;
};
} // namespace Gedim
 
#endif // __EIGENLUSOLVER_H