fe_dg.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_FE_DG_HH_
#define INCLUDE_IDEAL_II_FE_FE_DG_HH_
 
#include <deal.II/fe/fe.h>
#include <deal.II/fe/fe_dgq.h>
 
#include <memory>
 
namespace idealii::spacetime
{
  template <int dim>
  class DG_FiniteElement
  {
  public:
    // Info/Todo: With proper Quadrature class could be extended to RaudauLeft
    // and RadauRight Would need Radau based on Legendre polynomial roots
    enum support_type
    {
      Legendre,
      Lobatto,
      RadauLeft,
      RadauRight
    };
    DG_FiniteElement(std::shared_ptr<dealii::FiniteElement<dim>> fe_space,
                     const unsigned int                          r,
                     support_type type = support_type::Lobatto);
    std::shared_ptr<dealii::FiniteElement<dim>>
    spatial();
    std::shared_ptr<dealii::FiniteElement<1>>
    temporal();
    const unsigned int dofs_per_cell;
    support_type
    type();
 
  private:
    std::shared_ptr<dealii::FiniteElement<dim>> _fe_space;
    std::shared_ptr<dealii::FiniteElement<1>>   _fe_time;
    support_type                                _type;
  };
} // namespace idealii::spacetime
 
#endif /* INCLUDE_IDEAL_II_FE_FE_DG_HH_ */