limiters.hpp

Go to the documentation of this file.

#pragma once
 
#include <deal.II/base/parameter_handler.h>
#include <deal.II/base/tensor.h>
#include <deal.II/base/vectorization.h>
 
#include <deal.II/matrix_free/matrix_free.h>
 
#include <meltpooldg/utilities/better_enum.hpp>
#include <meltpooldg/utilities/dealii_tensor.hpp>
#include <meltpooldg/utilities/matrix_free_util.hpp>
 
#include <array>
#include <vector>
 
namespace MeltPoolDG::Utilities
{
  BETTER_ENUM(LimiterType, char, tvd_minmod, tvb_minmod)
 
  
  template <typename number>
  struct LimiterData
  {
    bool apply_limiter = false;
 
    std::vector<number> tvb_constant;
 
    LimiterType type = LimiterType::tvd_minmod;
 
    void
    add_parameters(dealii::ParameterHandler &prm);
  };
 
  template <int n_components, typename TensorType, typename Container>
  TensorType
  tvd_minmod(const Container &values);
 
  template <typename number, int n_components, typename TensorType, typename Container>
  TensorType
  tvb_minmod(const Container           &values,
             const std::vector<number> &tvb_constant,
             const number               cell_size);
 
  template <int dim, int n_components, typename number>
  std::array<dealii::Tensor<1, n_components, number>, dim>
  compute_minmod_type_limited_slopes(
    const std::vector<std::pair<dealii::Tensor<1, n_components, number>,
                                dealii::Tensor<1, n_components, dealii::Tensor<1, dim, number>>>>
                                                                    &cell_average_values,
    const typename dealii::Triangulation<dim>::active_cell_iterator &cell,
    const LimiterData<number>                                       &limiter_data);
 
  template <int dim,
            int n_components,
            typename number,
            typename VectorizedArrayType = dealii::VectorizedArray<number>,
            typename VectorType          = dealii::LinearAlgebra::distributed::Vector<number>>
  void
  apply_minmod_type_limiter(const MatrixFreeContext<dim, number> &mf_context,
                            VectorType                           &dst,
                            const VectorType                     &src,
                            const LimiterData<number>            &limiter_data);
 
} // namespace MeltPoolDG::Utilities