pointmapper.hh

#ifndef DUNE_FEM_POINTMAPPER_HH
#define DUNE_FEM_POINTMAPPER_HH
 
//- system includes
#include <vector>
 
//- Dune includes
#include <dune/common/version.hh>
#include <dune/fem/quadrature/quadrature.hh>
#include <dune/common/hash.hh>
 
namespace Dune
{
 
  namespace Fem
  {
 
    struct QuadratureKey
    {
      QuadratureKey ( const GeometryType &geoType, const size_t id )
      : id_( ((topologyId( geoType ) >> 1) << 16) + id )
      {
        assert( id < (1 << 16) );
      }
 
      bool operator< ( const QuadratureKey &other ) const
      {
        return (id_ < other.id_);
      }
 
      bool operator> ( const QuadratureKey &other ) const
      {
        return (id_ > other.id_);
      }
 
      bool operator== ( const QuadratureKey &other ) const
      {
        return (id_ == other.id_);
      }
 
      friend std::ostream &operator<< ( std::ostream &out, const QuadratureKey &key )
      {
        return out << "(topologyId " << ((key.id_ >> 16) << 1) << ", quadId " << (key.id_ & ((1u << 16)-1)) << ")";
      }
 
      inline friend std::size_t hash_value(const QuadratureKey& arg)
      {
        std::size_t seed = 0;
        hash_combine(seed,arg.id_);
        return seed;
      }
 
    protected:
      static unsigned int topologyId ( const GeometryType &type )
      {
        return type.id();
      }
 
      const size_t id_;
    };
 
 
 
    template <class ct, int dim>
    struct CachingTraits {
      typedef IntegrationPointListImp<ct, dim> QuadratureType;
      typedef typename QuadratureType::CoordinateType PointType;
      typedef std::vector<PointType>    PointVectorType;
      typedef std::vector<size_t>       MapperType;
      typedef std::pair< MapperType, MapperType > MapperPairType;
      typedef std::vector<MapperType>   MapperVectorType;
 
      typedef QuadratureKey QuadratureKeyType;
 
      // minimal twist is -4 for hexahedrons
      // so we add 4 to start from zero
      enum { twistOffset_ = 4 };
 
    };
 
  } // namespace Fem
 
} // namespace Dune
 
DUNE_DEFINE_HASH(DUNE_HASH_TEMPLATE_ARGS(),DUNE_HASH_TYPE(Dune::Fem::QuadratureKey))
 
#endif // #ifndef DUNE_FEM_POINTMAPPER_HH