power.hh

// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
// vi: set et ts=4 sw=2 sts=2:
 
#ifndef DUNE_LOCALFUNCTIONS_META_POWER_HH
#define DUNE_LOCALFUNCTIONS_META_POWER_HH
 
#include <cstddef>
#include <memory>
 
#include <dune/geometry/type.hh>
 
#include <dune/localfunctions/meta/power/basis.hh>
#include <dune/localfunctions/meta/power/coefficients.hh>
#include <dune/localfunctions/meta/power/interpolation.hh>
 
namespace Dune {
 
 
  template<class Backend, std::size_t dimR>
  class PowerFiniteElement {
  public:
    struct Traits {
      typedef PowerBasis<typename Backend::Traits::Basis, dimR> Basis;
      typedef PowerCoefficients Coefficients;
      typedef PowerInterpolation<typename Backend::Traits::Interpolation,
          typename Basis::Traits> Interpolation;
    };
  private:
    std::shared_ptr<const Backend> backend;
    typename Traits::Basis basis_;
    typename Traits::Coefficients coefficients_;
    typename Traits::Interpolation interpolation_;
 
  public:
 
    PowerFiniteElement(const Backend &backend_) :
      backend(new Backend(backend_)),
      basis_(backend->basis()),
      coefficients_(backend->coefficients(), dimR),
      interpolation_(backend->interpolation())
    { }
 
 
    PowerFiniteElement(const std::shared_ptr<const Backend> &backendSPtr) :
      backend(backendSPtr),
      basis_(backend->basis()),
      coefficients_(backend->coefficients(), dimR),
      interpolation_(backend->interpolation())
    { }
 
 
    const typename Traits::Basis& basis() const { return basis_; }
 
    const typename Traits::Coefficients& coefficients() const
    { return coefficients_; }
 
    const typename Traits::Interpolation& interpolation() const
    { return interpolation_; }
    GeometryType type() const { return backend->type(); }
  };
 
 
  template<class BackendFiniteElement, std::size_t dimR>
  class PowerFiniteElementFactory
  {
  public:
    typedef PowerFiniteElement<BackendFiniteElement, dimR> FiniteElement;
 
 
    const FiniteElement make(const BackendFiniteElement &backend) const
    { return FiniteElement(backend); }
 
    const FiniteElement
    make(const std::shared_ptr<const BackendFiniteElement> &backendSPtr) const
    { return FiniteElement(backendSPtr); }
 
  };
 
} // namespace Dune
 
#endif // DUNE_LOCALFUNCTIONS_META_POWER_HH