curvature_DG_operation.hpp

Go to the documentation of this file.

#pragma once
 
#include <deal.II/base/vectorization.h>
 
#include <deal.II/lac/la_parallel_vector.h>
#include <deal.II/lac/trilinos_sparse_matrix.h>
 
#include <deal.II/matrix_free/matrix_free.h>
 
#include <meltpooldg/core/scratch_data.hpp>
#include <meltpooldg/level_set/curvature_data.hpp>
#include <meltpooldg/level_set/curvature_operation_base.hpp>
#include <meltpooldg/level_set/helmholtz_DG_operator.hpp>
#include <meltpooldg/time_integration/solution_history.hpp>
#include <meltpooldg/utilities/utility_functions.hpp>
#include <meltpooldg/utilities/vector_tools.hpp>
 
namespace MeltPoolDG::LevelSet
{
  template <int dim, typename number>
  class CurvatureDGOperation
  {
  private:
    using VectorType       = dealii::LinearAlgebra::distributed::Vector<number>;
    using BlockVectorType  = dealii::LinearAlgebra::distributed::BlockVector<number>;
    using SparseMatrixType = dealii::TrilinosWrappers::SparseMatrix;
 
  public:
    CurvatureDGOperation(const ScratchData<dim, dim, number> &scratch_data_in,
                         const unsigned int                   curvature_dof_idx_in,
                         const unsigned int                   curvature_quad_idx_in,
                         const BlockVectorType               &solution_normal_vector_in,
                         const CurvatureData<number>         &curvature_data);
 
    void
    solve();
 
    const dealii::LinearAlgebra::distributed::Vector<number> &
    get_curvature() const;
 
    dealii::LinearAlgebra::distributed::Vector<number> &
    get_curvature();
 
    void
    reinit();
 
    void
    attach_vectors(std::vector<dealii::LinearAlgebra::distributed::Vector<number> *> &vectors);
 
  private:
    const ScratchData<dim, dim, number> &scratch_data;
    const BlockVectorType               &solution_normal_vector;
    const CurvatureData<number>          curvature_data;
 
    TimeIntegration::SolutionHistory<VectorType> solution_history;
 
    /*
     *  Based on the following indices the correct DoFHandler or quadrature rule from
     *  ScratchData<dim,dim,number> object is selected. This is important when
     * ScratchData<dim,dim,number> holds multiple DoFHandlers, quadrature rules, etc.
     */
    const unsigned int curvature_dof_idx;
    const unsigned int curvature_quad_idx;
 
    const HelmholtzDGOperator<dim, number> helmholtz_operator;
 
    void
    right_hand_side_domain(const dealii::MatrixFree<dim, number>       &data,
                           VectorType                                  &dst,
                           const BlockVectorType                       &src,
                           const std::pair<unsigned int, unsigned int> &cell_range) const;
 
    template <uint direction>
    void
    right_hand_side_inner_face(
      [[maybe_unused]] const dealii::MatrixFree<dim, number>       &data,
      [[maybe_unused]] VectorType                                  &dst,
      [[maybe_unused]] const BlockVectorType                       &src,
      [[maybe_unused]] const std::pair<unsigned int, unsigned int> &face_range) const
    {}
 
    template <uint direction>
    void
    right_hand_side_boundary_face(
      [[maybe_unused]] const dealii::MatrixFree<dim, number>       &data,
      [[maybe_unused]] VectorType                                  &dst,
      [[maybe_unused]] const BlockVectorType                       &src,
      [[maybe_unused]] const std::pair<unsigned int, unsigned int> &face_range) const
    {}
  };
} // namespace MeltPoolDG::LevelSet