time_integrator_data.hpp

Go to the documentation of this file.

#pragma once
 
#include <deal.II/base/parameter_handler.h>
 
#include <meltpooldg/linear_algebra/linear_solver_data.hpp>
#include <meltpooldg/linear_algebra/nonlinear_solver_data.hpp>
#include <meltpooldg/utilities/better_enum.hpp>
#include <meltpooldg/utilities/enum.hpp>
 
#include <string>
 
namespace MeltPoolDG::TimeIntegration
{
  BETTER_ENUM(TimeIntegratorSchemes,
              int,
              not_initialized,
              LSRK_stage_1_order_1,
              LSRK_stage_3_order_3, /* Kennedy, Carpenter, Lewis, 2000 */
              LSRK_stage_5_order_4, /* Kennedy, Carpenter, Lewis, 2000 */
              LSRK_stage_7_order_4, /* Tselios, Simos, 2007 */
              LSRK_stage_9_order_5, /* Kennedy, Carpenter, Lewis, 2000 */
              RKC_n_stages,         /* Sommeijer, Shampine, Verwer, 1997*/
              implicit_euler,
              explicit_euler,
              crank_nicolson,
              bdf_1,
              bdf_2,
              bdf_3,
              bdf_4,
              bdf_5,
              bdf_6,
              imex /* first order */)
 
  
  template <typename number = double>
  struct TimeIntegratorData
  {
    TimeIntegratorData() = default;
 
    explicit TimeIntegratorData(const TimeIntegratorSchemes scheme)
      : integrator_type(scheme)
    {}
 
 
    TimeIntegratorSchemes integrator_type = TimeIntegratorSchemes::not_initialized;
 
    unsigned int rkc_n_stages = 5;
    number       epsilon      = 2. / 13.;
 
    unsigned int                preconditioner_update_frequency = 100;
    NonlinearSolverData<number> nlsolver_data;
    LinearSolverData<number>    linear_solver_data;
 
    void
    add_parameters(dealii::ParameterHandler &prm)
    {
      prm.enter_subsection("time integration");
      {
        prm.add_parameter("type", integrator_type, "Name of the time integration scheme.");
        prm.add_parameter("preconditioner update frequency",
                          preconditioner_update_frequency,
                          "Frequency at which the preconditioner gets updated.");
        prm.add_parameter("RKC n stages", rkc_n_stages, "Number of stages for the RKC scheme.");
        prm.add_parameter("epsilon", epsilon, "Damping parameter for the RKC scheme.");
        nlsolver_data.add_parameters(prm);
        linear_solver_data.add_parameters(prm);
      }
      prm.leave_subsection();
    }
  };
 
} // namespace MeltPoolDG::TimeIntegration