reinitialization_DG_diffusion_operator.hpp

Go to the documentation of this file.

#pragma once
 
#include <deal.II/matrix_free/operators.h>
 
#include <deal.II/numerics/vector_tools.h>
 
#include <meltpooldg/core/scratch_data.hpp>
#include <meltpooldg/level_set/reinitialization_data.hpp>
#include <meltpooldg/utilities/fe_integrator.hpp>
 
namespace MeltPoolDG::LevelSet
{
  template <int dim, typename number>
  class ReinitializationDGDiffusionOperator
  {
  public:
    using VectorType      = dealii::LinearAlgebra::distributed::Vector<number>;
    using BlockVectorType = dealii::LinearAlgebra::distributed::BlockVector<number>;
 
    ReinitializationDGDiffusionOperator(
      const MeltPoolDG::ScratchData<dim, dim, number> &scratch_datain,
      const ReinitializationData<number>              &reinit_data_in,
      const unsigned int                               reinit_dof_idx_in,
      const unsigned int                               reinit_quad_idx_in,
      const VectorType                                &curvature_in,
      const BlockVectorType                           &normal_vector_in,
      const VectorType                                &smooth_signum_in);
 
    void
    compute_diffusitivity_value();
 
    void
    compute_penalty_parameter();
 
    void
    set_field_functions([[maybe_unused]] const number time) const {};
 
    void
    apply_operator([[maybe_unused]] const number                          time,
                   VectorType                                            &dst,
                   const VectorType                                      &src,
                   const std::function<void(unsigned int, unsigned int)> &func = {}) const;
 
    void
    apply_dirichlet_boundary_operator([[maybe_unused]] const number      time,
                                      [[maybe_unused]] VectorType       &dst,
                                      [[maybe_unused]] const VectorType &src) const {};
 
    number
    get_max_diffusitivity() const;
 
    bool update_field_functions = false;
 
  private:
    const MeltPoolDG::ScratchData<dim, dim, number> &scratch_data;
    const ReinitializationData<number>              &reinit_data;
 
    const unsigned int reinit_dof_idx;
    const unsigned int reinit_quad_idx;
 
    mutable VectorType                                             diffusitivity;
    mutable dealii::AlignedVector<dealii::VectorizedArray<number>> array_penalty_parameter;
 
    const VectorType      &curvature;
    const BlockVectorType &normal_vector;
    const VectorType      &smooth_signum;
 
    void
    local_apply_domain(const dealii::MatrixFree<dim, number>       &data,
                       VectorType                                  &dst,
                       const VectorType                            &src,
                       const std::pair<unsigned int, unsigned int> &cell_range) const;
 
    void
    local_apply_inner_face(const dealii::MatrixFree<dim, number>       &data,
                           VectorType                                  &dst,
                           const VectorType                            &src,
                           const std::pair<unsigned int, unsigned int> &cell_range) const;
 
    void
    local_apply_boundary_face(const dealii::MatrixFree<dim, number>       &data,
                              VectorType                                  &dst,
                              const VectorType                            &src,
                              const std::pair<unsigned int, unsigned int> &cell_range) const;
  };
} // namespace MeltPoolDG::LevelSet