heat_operation_base.hpp

Go to the documentation of this file.

#pragma once
 
#include <deal.II/base/function.h>
 
#include <deal.II/lac/la_parallel_block_vector.h>
#include <deal.II/lac/la_parallel_vector.h>
 
#include <meltpooldg/core/scratch_data.hpp>
#include <meltpooldg/level_set/nearest_point_data.hpp>
#include <meltpooldg/post_processing/generic_data_out.hpp>
 
#include <vector>
 
namespace MeltPoolDG::Heat
{
  template <int dim, typename number>
  class HeatOperationBase
  {
  public:
    using VectorType      = dealii::LinearAlgebra::distributed::Vector<number>;
    using BlockVectorType = dealii::LinearAlgebra::distributed::BlockVector<number>;
 
    virtual ~HeatOperationBase() = default;
 
    virtual void
    reinit() = 0;
 
    virtual void
    distribute_dofs(ScratchData<dim, dim, number> &scratch_data) const = 0;
 
    virtual void
    setup_constraints(ScratchData<dim, dim, number> &scratch_data) const = 0;
 
    virtual void
    set_initial_condition(const dealii::Function<dim> &initial_temperature) = 0;
 
    virtual void
    distribute_constraints() = 0;
 
    virtual void
    init_time_advance() = 0;
 
    virtual void
    solve() = 0;
 
    virtual void
    register_interface_projection_data(
      const VectorType                         &distance,
      const BlockVectorType                    &normal_vector,
      const LevelSet::NearestPointData<double> &nearest_point_data) = 0;
 
    virtual void
    compute_interface_temperature() = 0;
 
    virtual void
    attach_vectors(std::vector<VectorType *> &vectors) = 0;
 
    virtual void
    attach_output_vectors(GenericDataOut<dim, number> &data_out) const = 0;
 
    virtual void
    attach_output_vectors_failed_step(GenericDataOut<dim, number> &data_out) const = 0;
 
    /*
     * getters
     */
    virtual const VectorType &
    get_temperature() const = 0;
 
    virtual VectorType &
    get_temperature() = 0;
 
    virtual const VectorType &
    get_interface_temperature() const = 0;
 
    virtual VectorType &
    get_interface_temperature() = 0;
 
    virtual const VectorType &
    get_heat_source() const = 0;
 
    virtual VectorType &
    get_heat_source() = 0;
  };
} // namespace MeltPoolDG::Heat