spacetime_fe_values.hh

// ---------------------------------------------------------------------
//
// Copyright (C) 2022 - 2023 by the ideal.II authors
//
// This file is part of the ideal.II library.
//
// The ideal.II library is free software; you can use it, redistribute
// it, and/or modify it under the terms of the GNU Lesser General
// Public License as published by the Free Software Foundation; either
// version 3.0 of the License, or (at your option) any later version.
// The full text of the license can be found in the file LICENSE.md at
// the top level directory of ideal.II.
//
// ---------------------------------------------------------------------
 
#ifndef INCLUDE_IDEAL_II_FE_SPACETIME_FE_VALUES_HH_
#define INCLUDE_IDEAL_II_FE_SPACETIME_FE_VALUES_HH_
 
#include <ideal.II/base/quadrature_lib.hh>
 
#include <ideal.II/fe/fe_dg.hh>
 
#include <deal.II/fe/fe_values.h>
 
namespace idealii::spacetime
{
 
  template <int dim>
  class FEValues
  {
  public:
    FEValues(DG_FiniteElement<dim>    &fe,
             Quadrature<dim>          &quad,
             const dealii::UpdateFlags uflags);
 
    void
    reinit_space(const typename dealii::TriaIterator<
                 dealii::DoFCellAccessor<dim, dim, false>> &cell_space);
    void
    reinit_time(
      const typename dealii::TriaIterator<dealii::DoFCellAccessor<1, 1, false>>
        &cell_time);
 
    double
    shape_value(unsigned int function_no, unsigned int point_no);
 
    double
    shape_dt(unsigned int function_no, unsigned int point_no);
 
    dealii::Tensor<1, dim>
    shape_space_grad(unsigned int function_no, unsigned int point_no);
 
    template <class InputVector>
    void
    get_function_values(
      const InputVector &fe_function,
      std::vector<dealii::Vector<typename InputVector::value_type>> &values)
      const;
 
    template <class InputVector>
    void
    get_function_dt(
      const InputVector &fe_function,
      std::vector<dealii::Vector<typename InputVector::value_type>> &values)
      const;
 
    template <class InputVector>
    void
    get_function_space_gradients(
      const InputVector &fe_function,
      std::vector<
        std::vector<dealii::Tensor<1, dim, typename InputVector::value_type>>>
        &gradients) const;
 
 
    typename dealii::FEValuesViews::Scalar<dim>::value_type
    scalar_value(const typename dealii::FEValuesExtractors::Scalar &extractor,
                 unsigned int                                       function_no,
                 unsigned int                                       point_no);
 
    typename dealii::FEValuesViews::Scalar<dim>::value_type
    scalar_dt(const typename dealii::FEValuesExtractors::Scalar &extractor,
              unsigned int                                       function_no,
              unsigned int                                       point_no);
    typename dealii::FEValuesViews::Scalar<dim>::gradient_type
    scalar_space_grad(
      const typename dealii::FEValuesExtractors::Scalar &extractor,
      unsigned int                                       function_no,
      unsigned int                                       point_no);
 
    typename dealii::FEValuesViews::Vector<dim>::value_type
    vector_value(const typename dealii::FEValuesExtractors::Vector &extractor,
                 unsigned int                                       function_no,
                 unsigned int                                       point_no);
 
    typename dealii::FEValuesViews::Vector<dim>::value_type
    vector_dt(const typename dealii::FEValuesExtractors::Vector &extractor,
              unsigned int                                       function_no,
              unsigned int                                       point_no);
 
    typename dealii::FEValuesViews::Vector<dim>::divergence_type
    vector_divergence(
      const typename dealii::FEValuesExtractors::Vector &extractor,
      unsigned int                                       function_no,
      unsigned int                                       point_no);
 
    typename dealii::FEValuesViews::Vector<dim>::gradient_type
    vector_space_grad(
      const typename dealii::FEValuesExtractors::Vector &extractor,
      unsigned int                                       function_no,
      unsigned int                                       point_no);
 
    typename dealii::FEValuesViews::Vector<dim>::curl_type
    vector_space_curl(
      const typename dealii::FEValuesExtractors::Vector &extractor,
      unsigned int                                       function_no,
      unsigned int                                       point_no);
 
    double
    time_quadrature_point(unsigned int quadrature_point);
 
    dealii::Point<dim>
    space_quadrature_point(unsigned int quadrature_point);
 
    double
    JxW(const unsigned int quadrature_point);
    void
    get_local_dof_indices(
      std::vector<dealii::types::global_dof_index> &indices);
 
    std::shared_ptr<dealii::FEValues<dim>>
    spatial();
    std::shared_ptr<dealii::FEValues<1>>
    temporal();
 
  private:
    DG_FiniteElement<dim> &_fe;
    Quadrature<dim>       &_quad;
 
    std::shared_ptr<dealii::FEValues<dim>> _fev_space;
    std::shared_ptr<dealii::FEValues<1>>   _fev_time;
 
    std::vector<dealii::types::global_dof_index> local_space_dof_index;
    std::vector<dealii::types::global_dof_index> local_time_dof_index;
 
    unsigned int       n_dofs_space;
    unsigned int       time_cell_index;
    const unsigned int n_dofs_space_cell;
    const unsigned int n_quads_space;
 
  public:
    unsigned int n_quadrature_points;
  };
 
  template <int dim>
  class FEJumpValues
  {
  public:
    FEJumpValues(DG_FiniteElement<dim>    &fe,
                 Quadrature<dim>          &quad,
                 const dealii::UpdateFlags uflags);
 
    void
    reinit_space(const typename dealii::TriaIterator<
                 dealii::DoFCellAccessor<dim, dim, false>> &cell_space);
 
    void
    reinit_time(
      const typename dealii::TriaIterator<dealii::DoFCellAccessor<1, 1, false>>
        &cell_time);
 
    double
    shape_value_plus(unsigned int function_no, unsigned int point_no);
 
    double
    shape_value_minus(unsigned int function_no, unsigned int point_no);
 
    template <class InputVector>
    void
    get_function_values_minus(
      const InputVector &fe_function,
      std::vector<dealii::Vector<typename InputVector::value_type>> &values)
      const;
 
    template <class InputVector>
    void
    get_function_values_plus(
      const InputVector &fe_function,
      std::vector<dealii::Vector<typename InputVector::value_type>> &values)
      const;
 
    typename dealii::FEValuesViews::Scalar<dim>::value_type
    scalar_value_plus(
      const typename dealii::FEValuesExtractors::Scalar &extractor,
      unsigned int                                       function_no,
      unsigned int                                       point_no);
 
    typename dealii::FEValuesViews::Scalar<dim>::value_type
    scalar_value_minus(
      const typename dealii::FEValuesExtractors::Scalar &extractor,
      unsigned int                                       function_no,
      unsigned int                                       point_no);
 
    typename dealii::FEValuesViews::Vector<dim>::value_type
    vector_value_plus(
      const typename dealii::FEValuesExtractors::Vector &extractor,
      unsigned int                                       function_no,
      unsigned int                                       points_no);
 
    typename dealii::FEValuesViews::Vector<dim>::value_type
    vector_value_minus(
      const typename dealii::FEValuesExtractors::Vector &extractor,
      unsigned int                                       function_no,
      unsigned int                                       point_no);
    std::shared_ptr<dealii::FEValues<dim>>
    spatial();
 
    std::shared_ptr<dealii::FEValues<1>>
    temporal();
 
    double
    JxW(const unsigned int quadrature_point);
 
    unsigned int n_quadrature_points;
 
  private:
    unsigned int           n_dofs_space;
    DG_FiniteElement<dim> &_fe;
    Quadrature<dim>       &_quad;
 
    std::shared_ptr<dealii::FEValues<dim>> _fev_space;
    std::shared_ptr<dealii::FEValues<1>>   _fev_time;
 
    std::vector<dealii::types::global_dof_index> local_space_dof_index;
    std::vector<dealii::types::global_dof_index> local_time_dof_index;
  };
 
  template <int dim>
  class FEFaceValues
  {
  public:
    FEFaceValues(DG_FiniteElement<dim>    &fe,
                 Quadrature<dim - 1>      &quad,
                 const dealii::UpdateFlags uflags,
                 const dealii::UpdateFlags additional_flags);
 
    void
    reinit_space(const typename dealii::TriaIterator<
                   dealii::DoFCellAccessor<dim, dim, false>> &cell_space,
                 const unsigned int                           face_no);
    void
    reinit_time(
      const typename dealii::TriaIterator<dealii::DoFCellAccessor<1, 1, false>>
        &cell_time);
 
    double
    shape_value(unsigned int function_no, unsigned int point_no);
 
    typename dealii::FEValuesViews::Scalar<dim>::value_type
    scalar_value(const typename dealii::FEValuesExtractors::Scalar &extractor,
                 unsigned int                                       function_no,
                 unsigned int                                       point_no);
 
    typename dealii::FEValuesViews::Vector<dim>::value_type
    vector_value(const typename dealii::FEValuesExtractors::Vector &extractor,
                 unsigned int                                       function_no,
                 unsigned int                                       point_no);
 
    double
    time_quadrature_point(unsigned int quadrature_point);
 
    dealii::Point<dim>
    space_quadrature_point(unsigned int quadrature_point);
 
    const dealii::Tensor<1, dim> &
    space_normal_vector(unsigned int i);
 
    double
    JxW(const unsigned int quadrature_point);
    void
    get_local_dof_indices(
      std::vector<dealii::types::global_dof_index> &indices);
 
    std::shared_ptr<dealii::FEFaceValues<dim>>
    spatial();
    std::shared_ptr<dealii::FEValues<1>>
    temporal();
 
  private:
    DG_FiniteElement<dim> &_fe;
    Quadrature<dim - 1>   &_quad;
 
    std::shared_ptr<dealii::FEFaceValues<dim>> _fev_space;
    std::shared_ptr<dealii::FEValues<1>>       _fev_time;
 
    std::vector<dealii::types::global_dof_index> local_space_dof_index;
    std::vector<dealii::types::global_dof_index> local_time_dof_index;
 
    unsigned int       n_dofs_space;
    unsigned int       time_cell_index;
    const unsigned int n_dofs_space_cell;
    const unsigned int n_quads_space;
 
  public:
    unsigned int n_quadrature_points;
  };
} // namespace idealii::spacetime
 
#endif /* INCLUDE_IDEAL_II_BASE_SPACETIME_QUADRATURE_HH_ */