Pardiso_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 __Pardiso_LUSolver_H
#define __Pardiso_LUSolver_H
 
#include "Eigen/Eigen"
#include "Gedim_Macro.hpp"
#include "ILinearSolver.hpp"
 
#if ENABLE_MKL
#include <Eigen/PardisoSupport>
#endif
 
namespace Gedim
{
template <typename Eigen_ArrayType = Eigen::VectorXd, typename Eigen_SparseArrayType = Eigen::SparseMatrix<double>>
class Pardiso_LUSolver final : public ILinearSolver
{
  private:
#if ENABLE_MKL
    Eigen::PardisoLU<Eigen_SparseArrayType> linearSolver; 
#endif
    const IArray *_rightHandSide; 
    IArray *_solution;            
 
  public:
    Pardiso_LUSolver();
    ~Pardiso_LUSolver();
 
    void Initialize(const ISparseArray &matrix,
                    const IArray &rightHandSide,
                    IArray &solution,
                    const ILinearSolver::Configuration &config = Configuration());
 
    ILinearSolver::SolutionInfo Solve() const;
 
    void Initialize(const ISparseArray &matrix, const ILinearSolver::Configuration &config = Configuration());
 
    ILinearSolver::SolutionInfo Solve(const IArray &rightHandSide, IArray &solution) const;
};
} // namespace Gedim
 
#endif // __Pardiso_LUSolver_H