time_integrator_base.hpp

Go to the documentation of this file.

 
#pragma once
#include <deal.II/lac/la_parallel_vector.h>
 
#include <meltpooldg/time_integration/solution_history.hpp>
#include <meltpooldg/time_integration/time_integrator_data.hpp>
 
#include <functional>
 
namespace MeltPoolDG::TimeIntegration
{
  template <typename number>
  class TimeIntegratorBase
  {
    using VectorType = dealii::LinearAlgebra::distributed::Vector<number>;
 
  public:
    explicit TimeIntegratorBase(const TimeIntegratorData<number> &time_integrator_data_in)
      : time_integrator_data(time_integrator_data_in)
    {}
 
    virtual ~TimeIntegratorBase() = default;
 
    virtual unsigned int
    required_solution_history_size() const = 0;
 
    virtual void
    reinit(const VectorType &vector_template) = 0;
 
    virtual void
    reinit(const SolutionHistory<VectorType> &solution_history)
    {
      reinit(solution_history.get_current_solution());
    }
 
    virtual void
    perform_time_step(number                       current_time,
                      number                       time_step,
                      SolutionHistory<VectorType> &solution_history,
                      const std::function<void(number, number, VectorType &, const VectorType &)>
                        &pre_processing = {},
                      const std::function<void(number, number, VectorType &, const VectorType &)>
                        &post_processing = {}) = 0;
 
    TimeIntegratorSchemes
    get_integrator_type() const
    {
      return time_integrator_data.integrator_type;
    }
 
    void
    set_monitoring_vector(VectorType &monitoring_vector_in)
    {
      monitoring_vector = &monitoring_vector_in;
    }
 
    void
    reset_monitoring_vector()
    {
      monitoring_vector = nullptr;
    }
 
  protected:
    const TimeIntegratorData<number> time_integrator_data;
    VectorType                      *monitoring_vector = nullptr;
  };
} // namespace MeltPoolDG::TimeIntegration