typefromvertexcount.hh

// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
// vi: set et ts=4 sw=2 sts=2:
#ifndef DUNE_GEOMETRY_TYPE_FROM_VERTEX_COUNT_HH
#define DUNE_GEOMETRY_TYPE_FROM_VERTEX_COUNT_HH
 
#include <dune/geometry/type.hh>
 
namespace Dune {
 
  inline
  GeometryType geometryTypeFromVertexCount(unsigned int dim, unsigned int vertices)
  {
    switch (dim)
    {
      case 0 :
        return GeometryTypes::vertex;
      case 1 :
        return GeometryTypes::line;
      case 2 :
        switch (vertices) {
          case 3 :
            return GeometryTypes::triangle;
          case 4 :
            return GeometryTypes::quadrilateral;
          default :
            DUNE_THROW(NotImplemented, "2d elements with " << vertices << " corners are not supported!");
        }
      case 3 :
        switch (vertices) {
          case 4 :
            return GeometryTypes::tetrahedron;
          case 5 :
            return GeometryTypes::pyramid;
          case 6 :
            return GeometryTypes::prism;
          case 8 :
            return GeometryTypes::hexahedron;
          default :
            DUNE_THROW(NotImplemented, "3d elements with " << vertices << " corners are not supported!");
        }
      default :
        DUNE_THROW(NotImplemented, "geometryTypeFromVertexCount works only up to dim=3");
    }
  }
 
}
 
#endif // DUNE_GEOMETRY_TYPE_FROM_VERTEX_COUNT_HH