linear_solver_data.hpp

Go to the documentation of this file.

#pragma once
 
#include <meltpooldg/utilities/enum.hpp>
 
namespace MeltPoolDG
{
  BETTER_ENUM(PreconditionerType,
              char,
              // No preconditioner.
              Identity,
              // Algebraic multigrid preconditioner from the Trilinos package ...
              AMG,
              // Incomplete LU factorization preconditioner from the Trilinos package ...
              ILU,
              // Use the inverse diagonal of the system matrix as preconditioner ...
              Diagonal)
 
  BETTER_ENUM(LinearSolverType,
              char,
              // conjugate-gradient
              CG,
              // generalized minimal residual
              GMRES)
 
  BETTER_ENUM(LinearSolverMonitorType,
              char,
              // do not monitor history (production run)
              none,
              // print first and last residual
              reduced,
              // print full history
              all)
 
  template <typename number = double>
  struct LinearSolverData
  {
    bool               do_matrix_free      = true;
    PreconditionerType preconditioner_type = PreconditionerType::Identity;
    LinearSolverType   solver_type         = LinearSolverType::GMRES;
    unsigned int       max_iterations      = 10000;
    number             rel_tolerance       = 1e-12;
    number             abs_tolerance       = 1e-20;
 
    LinearSolverMonitorType monitor_type = LinearSolverMonitorType::none;
 
    void
    add_parameters(dealii::ParameterHandler &prm);
  };
 
} // namespace MeltPoolDG