reinitialization_geometric_operation.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/reinitialization_data.hpp>
#include <meltpooldg/level_set/reinitialization_operation_base.hpp>
#include <meltpooldg/level_set/signed_distance_solver.hpp>
#include <meltpooldg/post_processing/generic_data_out.hpp>
#include <meltpooldg/utilities/journal.hpp>
 
namespace MeltPoolDG::LevelSet
{
  using namespace dealii;
 
  template <int dim, typename number>
  class ReinitializationGeometricOperation : public ReinitializationOperationBase<dim, number>
  {
  public:
    using VectorType = dealii::LinearAlgebra::distributed::Vector<number>;
 
    using BlockVectorType = dealii::LinearAlgebra::distributed::BlockVector<number>;
 
    ReinitializationGeometricOperation(const ScratchData<dim, dim, number> &scratch_data_in,
                                       const ReinitializationData<number>  &reinit_data_in,
                                       const unsigned int                   ls_dof_idx_in,
                                       const unsigned int                   ls_quad_idx_in);
    void
    reinit() override;
 
    void
    solve() override;
 
    void
    set_initial_condition(const VectorType &solution_level_set_in) override;
 
    void
    set_initial_condition(const Function<dim> &initial_field_function) override;
 
    const VectorType &
    get_level_set() const override;
 
    VectorType &
    get_level_set() override;
 
    void
    attach_vectors([[maybe_unused]] std::vector<VectorType *> &vectors) override;
 
    void
    attach_output_vectors(GenericDataOut<dim, number> &data_out) const override;
 
  private:
    void
    create_solver();
 
    void
    transform_signed_distance_to_tanh_level_set();
 
  private:
    std::unique_ptr<SignedDistanceSolver<dim, VectorType>> signed_distance_solver;
 
    const ScratchData<dim, dim, number> &scratch_data;
 
    const ReinitializationData<number> &reinit_data;
 
    unsigned int ls_dof_idx;
 
    unsigned int ls_quad_idx;
 
    VectorType level_set;
  };
 
} // namespace MeltPoolDG::LevelSet