reinitialization_elliptic_operator_CG_newton.hpp

Go to the documentation of this file.

#pragma once
 
#include <deal.II/fe/fe_dgq.h>
#include <deal.II/fe/fe_system.h>
 
#include <deal.II/lac/generic_linear_algebra.h>
 
#include <deal.II/matrix_free/fe_evaluation.h>
#include <deal.II/matrix_free/fe_point_evaluation.h>
#include <deal.II/matrix_free/matrix_free.h>
#include <deal.II/matrix_free/operators.h>
 
#include <meltpooldg/core/operator_base.hpp>
#include <meltpooldg/core/scratch_data.hpp>
#include <meltpooldg/level_set/reinitialization_data.hpp>
#include <meltpooldg/utilities/fe_integrator.hpp>
#include <meltpooldg/utilities/vector_tools.hpp>
 
 
namespace MeltPoolDG::LevelSet
{
  template <int dim, typename number>
  class ReinitializationEllipticOperatorNewton : public OperatorMatrixFree<dim, number>
  {
    using OperatorMatrixFree<dim, number>::vmult;
    using OperatorMatrixFree<dim, number>::create_rhs;
    using OperatorMatrixFree<dim, number>::compute_inverse_diagonal_from_matrixfree;
 
  public:
    using VectorType          = dealii::LinearAlgebra::distributed::Vector<number>;
    using VectorizedArrayType = dealii::VectorizedArray<number>;
    using MappingInfoType     = CutUtil::MappingInfoType<dim, number>;
    using PointEvaluationType = dealii::FEPointEvaluation<1, dim, dim, VectorizedArrayType>;
 
 
    ReinitializationEllipticOperatorNewton(
      const MeltPoolDG::ScratchData<dim, dim, number> &scratch_data_in,
      const ReinitializationData<number>              &reinit_data_in,
      const unsigned int                               reinit_dof_idx_in,
      const unsigned int                               reinit_quad_idx_in,
      const MappingInfoType                           &mapping_info_surface_in,
      const unsigned int                               ls_dof_idx_in);
 
    void
    vmult(VectorType &dst, const VectorType &src) const final;
 
    void
    create_residual(VectorType &rhs, const VectorType &level_set) const;
 
    void
    reinit() final;
 
    VectorType solution_old;
 
    void
    compute_system_matrix_from_matrixfree(
      dealii::TrilinosWrappers::SparseMatrix &system_matrix) const final;
 
    void
    compute_inverse_diagonal_from_matrixfree(VectorType &diagonal) const final;
 
  private:
    template <int n_components>
    VectorizedArrayType
    evaluate_coefficient(const FECellIntegrator<dim, n_components, number> &phi_old,
                         const unsigned int                                 q_index) const;
 
    template <int n_components>
    VectorizedArrayType
    evaluate_coefficient_derivative(const FECellIntegrator<dim, n_components, number> &phi_old,
                                    const unsigned int q_index) const;
 
    template <int n_components>
    VectorizedArrayType
    evaluate_coefficient_derivative(const FEFaceIntegrator<dim, n_components, number> &phi_old,
                                    const unsigned int q_index) const;
 
    void
    residual_laplace_operation(FECellIntegrator<dim, 1, number> &cell_eval) const;
 
    void
    tangent_laplace_operation(FECellIntegrator<dim, 1, number> &cell_eval,
                              FECellIntegrator<dim, 1, number> &phi_old) const;
 
    void
    cell_operation(FECellIntegrator<dim, 1, number> &interface_penalty,
                   FECellIntegrator<dim, 1, number> &cell_eval,
                   PointEvaluationType              &interface_penalty_surface) const;
    void
    interface_penalty_cell_operation(PointEvaluationType              &interface_penalty_surface,
                                     FECellIntegrator<dim, 1, number> &interface_penalty,
                                     const unsigned int                lane,
                                     const number                      penalty_coefficient) const;
 
    const ScratchData<dim, dim, number> &scratch_data;
 
    const ReinitializationData<number> &reinit_data;
 
    const unsigned int reinit_quad_idx;
 
    const MappingInfoType &mapping_info_surface;
 
    dealii::FE_Q<dim> fe_point_level_set;
 
    const unsigned int n_dofs_per_cell;
 
    const unsigned int ls_dof_idx;
 
    VectorType zero_interface;
  };
} // namespace MeltPoolDG::LevelSet