pk.hh

// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
// vi: set et ts=4 sw=2 sts=2:
/* vim: set ai expandtab sw=4 ts=4: */
#ifndef DUNE_PK_LOCALFINITEELEMENT_HH
#define DUNE_PK_LOCALFINITEELEMENT_HH
 
#include "p0.hh"
#include "p1.hh"
#include "pk1d.hh"
#include "pk2d.hh"
#include "pk3d.hh"
 
namespace Dune
{
 
  template<class D, class R, int d, int k>
  class PkLocalFiniteElement
  {
  public:
    PkLocalFiniteElement()
    {}
 
    PkLocalFiniteElement(const unsigned int vertexmap[k+1])
    {}
  };
#if 0
  template<class D, class R>
  class PkLocalFiniteElement<D, R, 1, 1>
    : public P1LocalFiniteElement<D, R, 1>
  {
  public:
    PkLocalFiniteElement()
    {}
 
    PkLocalFiniteElement(const unsigned int vertexmap[2])
    {}
  };
 
  template<class D, class R>
  class PkLocalFiniteElement<D, R, 1, 0>
    : public P0LocalFiniteElement<D, R, 1>
  {
  public:
    PkLocalFiniteElement()
      : P0LocalFiniteElement<D,R,1>(GeometryType(0,1))
    {}
 
    PkLocalFiniteElement(const unsigned int vertexmap[2])
    {}
  };
#endif
  template<class D, class R, int k>
  class PkLocalFiniteElement<D, R, 1, k>
    : public Pk1DLocalFiniteElement<D, R, k>
  {
  public:
    PkLocalFiniteElement()
    {}
 
    PkLocalFiniteElement(const unsigned int vertexmap[2]) :
      Pk1DLocalFiniteElement<D, R, k>(vertexmap)
    {}
  };
 
  template<class D, class R, int k>
  class PkLocalFiniteElement<D, R, 2, k>
    : public Pk2DLocalFiniteElement<D, R, k>
  {
  public:
    PkLocalFiniteElement()
    {}
 
    PkLocalFiniteElement(const unsigned int vertexmap[3]) :
      Pk2DLocalFiniteElement<D, R, k>(vertexmap)
    {}
  };
 
  template<class D, class R, int k>
  class PkLocalFiniteElement<D, R, 3, k>
    : public Pk3DLocalFiniteElement<D, R, k>
  {
  public:
    PkLocalFiniteElement()
    {}
 
    PkLocalFiniteElement(const unsigned int vertexmap[4]) :
      Pk3DLocalFiniteElement<D, R, k>(vertexmap)
    {}
  };
 
}
 
#endif