lagrangecoefficients.hh

// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
// vi: set et ts=4 sw=2 sts=2:
#ifndef DUNE_LAGRANGECOEFFICIENTS_HH
#define DUNE_LAGRANGECOEFFICIENTS_HH
 
#include <vector>
 
#include <dune/common/fvector.hh>
 
#include <dune/geometry/type.hh>
 
#include <dune/localfunctions/utility/field.hh>
#include <dune/localfunctions/common/localkey.hh>
 
namespace Dune
{
 
  template< template <class,unsigned int> class LP,
      unsigned int dim, class F>
  struct LagrangeCoefficientsFactory
  {
    static const unsigned int dimension = dim;
    const typedef LP<F,dim> Object;
    typedef std::size_t Key;
 
    template< GeometryType::Id geometryId >
    static Object *create ( const Key &order )
    {
      if (order == 0 || !Object::template supports<geometryId>(order))
        return 0;
      typedef typename std::remove_const<Object>::type LagrangeCoefficients;
      LagrangeCoefficients *object = new LagrangeCoefficients(order);
      if ( !object->template build<geometryId>() )
      {
        delete object;
        object = nullptr;
      }
      return object;
    }
    static void release( Object *object ) { delete object; }
  };
 
}
 
#endif // DUNE_LAGRANGECOEFFICIENTS_HH