Skip to content
Snippets Groups Projects
Select Git revision
  • 6eedc659d9b7cd595d6b7aa513db8e7fb6c95de8
  • develop default protected
  • feature/gmsh-reader
  • origin/stage/bouguettaia
  • feature/kinetic-schemes
  • feature/reconstruction
  • feature/local-dt-fsi
  • feature/composite-scheme-sources
  • feature/composite-scheme-other-fluxes
  • feature/serraille
  • feature/variational-hydro
  • feature/composite-scheme
  • hyperplastic
  • feature/polynomials
  • feature/gks
  • feature/implicit-solver-o2
  • feature/coupling_module
  • feature/implicit-solver
  • feature/merge-local-dt-fsi
  • master protected
  • feature/escobar-smoother
  • v0.5.0 protected
  • v0.4.1 protected
  • v0.4.0 protected
  • v0.3.0 protected
  • v0.2.0 protected
  • v0.1.0 protected
  • Kidder
  • v0.0.4 protected
  • v0.0.3 protected
  • v0.0.2 protected
  • v0 protected
  • v0.0.1 protected
33 results

PolynomialBasis.hpp

Blame
  • PolynomialBasis.hpp 1.97 KiB
    #ifndef POLYNOMIALBASIS_HPP
    #define POLYNOMIALBASIS_HPP
    
    #include <algebra/TinyVector.hpp>
    #include <analysis/Polynomial.hpp>
    #include <utils/Messenger.hpp>
    
    enum class BasisType
    {
      canonical
    };
    
    template <size_t N>
    class PolynomialBasis
    {
     private:
      static_assert((N >= 0), "Number of elements in the basis must be non-negative");
      TinyVector<N + 1, Polynomial<N>> m_elements;
      BasisType m_basis_type;
      PUGS_INLINE
      constexpr PolynomialBasis<N>&
      _buildCanonicalBasis()
      {
        for (size_t i = 0; i <= N; i++) {
          TinyVector<N + 1> coefficients(zero);
          coefficients[i] = 1;
          elements()[i]   = Polynomial<N>(coefficients);
        }
        return *this;
      }
    
     public:
      PUGS_INLINE
      constexpr size_t
      size() const
      {
        return N + 1;
      }
    
      PUGS_INLINE
      constexpr size_t
      degree() const
      {
        return N;
      }
    
      PUGS_INLINE
      constexpr BasisType&
      type()
      {
        return m_basis_type;
      }
    
      PUGS_INLINE
      std::string_view
      displayType()
      {
        switch (m_basis_type) {
        case BasisType::canonical:
          return "canonical";
        default:
          return "unknown basis type";
        }
      }
    
      PUGS_INLINE
      constexpr const TinyVector<N + 1, Polynomial<N>>&
      elements() const
      {
        return m_elements;
      }
    
      PUGS_INLINE
      constexpr TinyVector<N + 1, Polynomial<N>>&
      elements()
      {
        return m_elements;
      }
    
      PUGS_INLINE
      constexpr PolynomialBasis<N>&
      build(BasisType basis_type)
      {
        type() = basis_type;
        switch (basis_type) {
        case BasisType::canonical: {
          return _buildCanonicalBasis();
          break;
        }
        default:
          throw UnexpectedError("unknown basis type");
        }
      }
      PUGS_INLINE constexpr PolynomialBasis(const TinyVector<N + 1, Polynomial<N>>& elements) noexcept
        : m_elements{elements}
      {}
    
      PUGS_INLINE
      constexpr PolynomialBasis(TinyVector<N + 1, Polynomial<N>>&& elements) noexcept : m_elements{elements} {}
    
      PUGS_INLINE
      constexpr PolynomialBasis() noexcept = default;
      ~PolynomialBasis()                   = default;
    };
    #endif   // POLYNOMIAL_HPP