MathLib::EigenLinearSolver Class Reference

Detailed Description

Definition at line 26 of file EigenLinearSolver.h.

#include <EigenLinearSolver.h>

Collaboration diagram for MathLib::EigenLinearSolver:

Public Member Functions
	EigenLinearSolver (std::string const &solver_name, EigenOption const &option)
	~EigenLinearSolver ()
void	setOption (const EigenOption &option)
EigenOption &	getOption ()
bool	compute (EigenMatrix &A, MathLib::LinearSolverBehaviour const linear_solver_behaviour)
bool	solve (EigenVector &b, EigenVector &x)
bool	solve (EigenMatrix &A, EigenVector &b, EigenVector &x, MathLib::LinearSolverBehaviour const linear_solver_behaviour=MathLib::LinearSolverBehaviour::RECOMPUTE)
	Computes and solves in a single call.
bool	canSolveRectangular () const
	Get, if the solver can handle rectangular equation systems.

Protected Member Functions
void	setRestart ()
void	setL ()
void	setS ()
void	setSmoothing ()
void	setAngle ()
void	setResidualUpdate ()

Protected Attributes
EigenOption	option_
std::unique_ptr< EigenLinearSolverBase >	solver_
bool	can_solve_rectangular_ = false

Constructor & Destructor Documentation

EigenLinearSolver()

MathLib::EigenLinearSolver::EigenLinearSolver ( std::string const & solver_name, EigenOption const & option )

explicit

Constructor

Parameters

solver_name	A name used as a prefix for command line options if there are such options available.
option	Eigen linear solver options.

Definition at line 501 of file EigenLinearSolver.cpp.

    : option_(option)
{
    using Matrix = EigenMatrix::RawMatrixType;
 
    switch (option_.solver_type)
    {
        case EigenOption::SolverType::SparseLU:
        {
            using SolverType =
                Eigen::SparseLU<Matrix, Eigen::COLAMDOrdering<int>>;
            solver_ = std::make_unique<
                details::EigenDirectLinearSolver<SolverType>>();
            can_solve_rectangular_ = false;
            return;
        }
        case EigenOption::SolverType::BiCGSTAB:
        case EigenOption::SolverType::BiCGSTABL:
        case EigenOption::SolverType::CG:
        case EigenOption::SolverType::GMRES:
        case EigenOption::SolverType::IDRS:
        case EigenOption::SolverType::IDRSTABL:
            solver_ =
                details::createIterativeSolver(option_.solver_type,
                                               option_.precon_type,
                                               option_.triangular_matrix_type);
            can_solve_rectangular_ = false;
            return;
        case EigenOption::SolverType::LeastSquareCG:
            solver_ =
                details::createIterativeSolver(option_.solver_type,
                                               option_.precon_type,
                                               option_.triangular_matrix_type);
            can_solve_rectangular_ = true;
            return;
        case EigenOption::SolverType::PardisoLU:
        {
#ifdef USE_MKL
            using SolverType = Eigen::PardisoLU<EigenMatrix::RawMatrixType>;
            solver_.reset(new details::EigenDirectLinearSolver<SolverType>);
            can_solve_rectangular_ = false;
            return;
#else
            OGS_FATAL(
                "The code is not compiled with Intel MKL. Linear solver type "
                "PardisoLU is not available.");
#endif
        }
    }
 
    OGS_FATAL("Invalid Eigen linear solver type. Aborting.");
}

References MathLib::EigenOption::BiCGSTAB, MathLib::EigenOption::BiCGSTABL, can_solve_rectangular_, MathLib::EigenOption::CG, MathLib::details::createIterativeSolver(), MathLib::EigenOption::GMRES, MathLib::EigenOption::IDRS, MathLib::EigenOption::IDRSTABL, MathLib::EigenOption::LeastSquareCG, option_, MathLib::EigenOption::PardisoLU, MathLib::EigenOption::precon_type, solver_, MathLib::EigenOption::solver_type, and MathLib::EigenOption::triangular_matrix_type.

~EigenLinearSolver()

MathLib::EigenLinearSolver::~EigenLinearSolver ( )

default

Member Function Documentation

canSolveRectangular()

bool MathLib::EigenLinearSolver::canSolveRectangular ( ) const

inline

Get, if the solver can handle rectangular equation systems.

Definition at line 75 of file EigenLinearSolver.h.

{ return can_solve_rectangular_; }

References can_solve_rectangular_.

compute()

bool MathLib::EigenLinearSolver::compute	(	EigenMatrix &	A,
		MathLib::LinearSolverBehaviour const	linear_solver_behaviour )

Performs the compute() step of the Eigen linear solver.

I.e., computes the (LU) decomposition in case of a direct solver, or computes the preconditioner of an iterative solver.

Definition at line 557 of file EigenLinearSolver.cpp.

{
    INFO("------------------------------------------------------------------");
    INFO("*** Eigen solver compute()");
 
    return solver_->compute(A.getRawMatrix(), option_, linear_solver_behaviour);
}

References MathLib::EigenMatrix::getRawMatrix(), and INFO().

getOption()

EigenOption & MathLib::EigenLinearSolver::getOption ( )

inline

get linear solver options

Definition at line 48 of file EigenLinearSolver.h.

{ return option_; }

References option_.

setAngle()

void MathLib::EigenLinearSolver::setAngle ( )

protected

setL()

void MathLib::EigenLinearSolver::setL ( )

protected

setOption()

void MathLib::EigenLinearSolver::setOption ( const EigenOption & option )

inline

copy linear solvers options

Definition at line 43 of file EigenLinearSolver.h.

{ option_ = option; }

References option_.

setResidualUpdate()

void MathLib::EigenLinearSolver::setResidualUpdate ( )

protected

setRestart()

void MathLib::EigenLinearSolver::setRestart ( )

protected

setS()

void MathLib::EigenLinearSolver::setS ( )

protected

setSmoothing()

void MathLib::EigenLinearSolver::setSmoothing ( )

protected

solve() [1/2]

bool MathLib::EigenLinearSolver::solve	(	EigenMatrix &	A,
		EigenVector &	b,
		EigenVector &	x,
		MathLib::LinearSolverBehaviour const	linear_solver_behaviour = MathLib::LinearSolverBehaviour::RECOMPUTE )

Computes and solves in a single call.

Definition at line 575 of file EigenLinearSolver.cpp.

{
    return solver_->compute(A.getRawMatrix(), option_,
                            linear_solver_behaviour) &&
           solver_->solve(b.getRawVector(), x.getRawVector(), option_);
}

References MathLib::EigenMatrix::getRawMatrix(), and MathLib::EigenVector::getRawVector().

solve() [2/2]

bool MathLib::EigenLinearSolver::solve	(	EigenVector &	b,
		EigenVector &	x )

Solves the linear system for the given right-hand side b and initial guess x.

Precondition

compute() must have been called before. (Not necessarily for every x and b separately, but for every new/changed matrix A.

Definition at line 567 of file EigenLinearSolver.cpp.

{
    INFO("------------------------------------------------------------------");
    INFO("*** Eigen solver solve()");
 
    return solver_->solve(b.getRawVector(), x.getRawVector(), option_);
}

References MathLib::EigenVector::getRawVector(), and INFO().

Member Data Documentation

can_solve_rectangular_

bool MathLib::EigenLinearSolver::can_solve_rectangular_ = false

protected

Definition at line 80 of file EigenLinearSolver.h.

Referenced by EigenLinearSolver(), and canSolveRectangular().

option_

EigenOption MathLib::EigenLinearSolver::option_

protected

Definition at line 78 of file EigenLinearSolver.h.

Referenced by EigenLinearSolver(), getOption(), and setOption().

solver_

std::unique_ptr<EigenLinearSolverBase> MathLib::EigenLinearSolver::solver_

protected

Definition at line 79 of file EigenLinearSolver.h.

Referenced by EigenLinearSolver().

---

The documentation for this class was generated from the following files:

• MathLib/LinAlg/Eigen/EigenLinearSolver.h
• MathLib/LinAlg/Eigen/EigenLinearSolver.cpp