reinitialization_DG_grad_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/utilities/fe_integrator.hpp>
 
namespace MeltPoolDG::LevelSet
{
  template <int dim, typename number>
  class RIGradOperator
  {
  public:
    using VectorType      = dealii::LinearAlgebra::distributed::Vector<number>;
    using BlockVectorType = dealii::LinearAlgebra::distributed::BlockVector<number>;
 
 
    RIGradOperator(const MeltPoolDG::ScratchData<dim, dim, number> &scratch_data_in,
                   const unsigned int                               reinit_dof_idx_in,
                   const unsigned int                               reinit_quad_idx_in);
 
 
    template <bool is_right, uint component>
    void
    apply(const VectorType &src, VectorType &dst);
 
  private:
    const MeltPoolDG::ScratchData<dim, dim, number> &scratch_data;
 
    const unsigned int reinit_dof_idx;
    const unsigned int reinit_quad_idx;
 
    template <uint component>
    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;
 
    template <bool is_right, uint component>
    void
    local_apply_inner_face(const dealii::MatrixFree<dim, number>       &data,
                           VectorType                                  &dst,
                           const VectorType                            &src,
                           const std::pair<unsigned int, unsigned int> &face_range) const;
 
    template <bool is_right, uint component>
    void
    local_apply_boundary_face(const dealii::MatrixFree<dim, number>       &data,
                              VectorType                                  &dst,
                              const VectorType                            &src,
                              const std::pair<unsigned int, unsigned int> &face_range) const;
  };
} // namespace MeltPoolDG::LevelSet