MathLib::details Namespace Reference

Classes
class	EigenDirectLinearSolver
	Template class for Eigen direct linear solvers. More...
class	EigenIterativeLinearSolver
	Template class for Eigen iterative linear solvers. More...

Typedefs
template<typename Mat , typename Precon >
using	EigenCGSolver = Eigen::ConjugateGradient< Mat, Eigen::Lower, Precon >

Functions
template<template< typename, typename > class Solver, typename Precon >
std::unique_ptr< EigenLinearSolverBase >	createIterativeSolver ()
template<template< typename, typename > class Solver>
std::unique_ptr< EigenLinearSolverBase >	createIterativeSolver (EigenOption::PreconType precon_type)
std::unique_ptr< EigenLinearSolverBase >	createIterativeSolver (EigenOption::SolverType solver_type, EigenOption::PreconType precon_type)

Typedef Documentation

EigenCGSolver

template<typename Mat , typename Precon >

using MathLib::details::EigenCGSolver = typedef Eigen::ConjugateGradient<Mat, Eigen::Lower, Precon>

Definition at line 190 of file EigenLinearSolver.cpp.

Function Documentation

createIterativeSolver() [1/3]

template<template< typename, typename > class Solver, typename Precon >

std::unique_ptr<EigenLinearSolverBase> MathLib::details::createIterativeSolver

(

)

Definition at line 159 of file EigenLinearSolver.cpp.

{
    using Slv =
        EigenIterativeLinearSolver<Solver<EigenMatrix::RawMatrixType, Precon>>;
    return std::make_unique<Slv>();
}

Referenced by MathLib::EigenLinearSolver::EigenLinearSolver(), and createIterativeSolver().

createIterativeSolver() [2/3]

template<template< typename, typename > class Solver>

std::unique_ptr<EigenLinearSolverBase> MathLib::details::createIterativeSolver

(

EigenOption::PreconType

precon_type

)

Definition at line 167 of file EigenLinearSolver.cpp.

{
    switch (precon_type)
    {
        case EigenOption::PreconType::NONE:
            return createIterativeSolver<Solver,
                                         Eigen::IdentityPreconditioner>();
        case EigenOption::PreconType::DIAGONAL:
            return createIterativeSolver<
                Solver, Eigen::DiagonalPreconditioner<double>>();
        case EigenOption::PreconType::ILUT:
            // TODO for this preconditioner further options can be passed.
            // see
            // https://eigen.tuxfamily.org/dox/classEigen_1_1IncompleteLUT.html
            return createIterativeSolver<Solver,
                                         Eigen::IncompleteLUT<double>>();
        default:
            OGS_FATAL("Invalid Eigen preconditioner type.");
    }
}

References createIterativeSolver(), MathLib::EigenOption::DIAGONAL, MathLib::EigenOption::ILUT, MathLib::EigenOption::NONE, and OGS_FATAL.

createIterativeSolver() [3/3]

std::unique_ptr<EigenLinearSolverBase> MathLib::details::createIterativeSolver	(	EigenOption::SolverType	solver_type,
		EigenOption::PreconType	precon_type
	)

Definition at line 192 of file EigenLinearSolver.cpp.

{
    switch (solver_type)
    {
        case EigenOption::SolverType::BiCGSTAB:
        {
            return createIterativeSolver<Eigen::BiCGSTAB>(precon_type);
        }
        case EigenOption::SolverType::CG:
        {
            return createIterativeSolver<EigenCGSolver>(precon_type);
        }
        case EigenOption::SolverType::GMRES:
        {
#ifdef USE_EIGEN_UNSUPPORTED
            return createIterativeSolver<Eigen::GMRES>(precon_type);
#else
            OGS_FATAL(
                "The code is not compiled with the Eigen unsupported modules. "
                "Linear solver type GMRES is not available.");
#endif
        }
        default:
            OGS_FATAL("Invalid Eigen iterative linear solver type. Aborting.");
    }
}

References MathLib::EigenOption::BiCGSTAB, MathLib::EigenOption::CG, MathLib::EigenOption::GMRES, and OGS_FATAL.