topologyutility.hh

// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
// vi: set et ts=4 sw=2 sts=2:
#ifndef DUNE_PDELAB_COMMON_TOPOLOGYUTILITY_HH
#define DUNE_PDELAB_COMMON_TOPOLOGYUTILITY_HH
 
#include <dune/geometry/type.hh>
 
namespace Dune {
 
  namespace PDELab {
 
 
 
    template<int dimension, unsigned int topologyId>
    struct BasicTypeFromDimensionAndTopologyId
    {
      static const bool isCube =
        ((topologyId ^ ((1 << dimension)-1)) >> 1 == 0);
 
      static const bool isSimplex =
        (topologyId | 1) == 1;
 
      static const GeometryType::BasicType value =
        isSimplex ? GeometryType::simplex
        : (
          isCube ? GeometryType::cube
          : GeometryType::none
          );
    };
 
 
  } // namespace PDELab
} // namespace Dune
 
 
#endif // DUNE_PDELAB_COMMON_TOPOLOGYUTILITY_HH