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 = Eigen::ConjugateGradient<Mat, Eigen::Lower, Precon>

Definition at line 403 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 372 of file EigenLinearSolver.cpp.

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

Referenced by MathLib::EigenLinearSolver::EigenLinearSolver(), createIterativeSolver(), 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 380 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://libeigen.gitlab.io/docs/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 405 of file EigenLinearSolver.cpp.

{
    switch (solver_type)
    {
        case EigenOption::SolverType::BiCGSTAB:
        {
            return createIterativeSolver<Eigen::BiCGSTAB>(precon_type);
        }
        case EigenOption::SolverType::BiCGSTABL:
        {
#ifdef USE_EIGEN_UNSUPPORTED
#if EIGEN_VERSION_AT_LEAST(3, 4, 90)
            return createIterativeSolver<Eigen::BiCGSTABL>(precon_type);
#else
            OGS_FATAL("BiCGSTABL requires at least Eigen version 3.4.90");
#endif
#else
            OGS_FATAL(
                "The code is not compiled with the Eigen unsupported modules. "
                "Linear solver type BiCGSTABL is not available.");
#endif
        }
        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
        }
        case EigenOption::SolverType::IDRS:
        {
#ifdef USE_EIGEN_UNSUPPORTED
#if EIGEN_VERSION_AT_LEAST(3, 4, 90)
            return createIterativeSolver<Eigen::IDRS>(precon_type);
#else
            OGS_FATAL("IDRS requires at least Eigen version 3.4.90");
#endif
#else
            OGS_FATAL(
                "The code is not compiled with the Eigen unsupported modules. "
                "Linear solver type IDRS is not available.");
#endif
        }
        case EigenOption::SolverType::IDRSTABL:
        {
#ifdef USE_EIGEN_UNSUPPORTED
#if EIGEN_VERSION_AT_LEAST(3, 4, 90)
            return createIterativeSolver<Eigen::IDRSTABL>(precon_type);
#else
            OGS_FATAL("IDRSTABL requires at least Eigen version 3.4.90");
#endif
#else
            OGS_FATAL(
                "The code is not compiled with the Eigen unsupported modules. "
                "Linear solver type IDRSTABL is not available.");
#endif
        }
        default:
            OGS_FATAL("Invalid Eigen iterative linear solver type. Aborting.");
    }
}

References MathLib::EigenOption::BiCGSTAB, MathLib::EigenOption::BiCGSTABL, MathLib::EigenOption::CG, createIterativeSolver(), MathLib::EigenOption::GMRES, MathLib::EigenOption::IDRS, MathLib::EigenOption::IDRSTABL, and OGS_FATAL.