topology.hh

#ifndef DUNE_ALU3DGRIDTOPOLOGY_HH
#define DUNE_ALU3DGRIDTOPOLOGY_HH
 
#include <utility>
 
#include <dune/alugrid/common/alugrid_assert.hh>
 
namespace Dune {
 
  // types of the elementes,
  // i.e . tetra or hexa, mixed is not implemeneted
  enum ALU3dGridElementType { tetra = 4, hexa = 7, mixed, error };
 
  template <ALU3dGridElementType type>
  struct EntityCount {};
 
  template <>
  struct EntityCount<tetra> {
    enum {numFaces = 4};
    enum {numVertices = 4};
    enum {numEdges = 6};
    enum {numVerticesPerFace = 3};
    enum {numEdgesPerFace = 3};
  };
 
  template <>
  struct EntityCount<hexa> {
    enum {numFaces = 6};
    enum {numVertices = 8};
    enum {numEdges = 12};
    enum {numVerticesPerFace = 4};
    enum {numEdgesPerFace = 4};
  };
 
 
  template <ALU3dGridElementType type>
  class ElementTopologyMapping
  {
  public:
    enum { numFaces = EntityCount<type>::numFaces };
    enum { numVertices = EntityCount<type>::numVertices };
    enum { numEdges = EntityCount<type>::numEdges };
    enum { numVerticesPerFace = EntityCount<type>::numVerticesPerFace };
 
    static int dune2aluFace(int index);
    static int alu2duneFace(int index);
 
    static int dune2aluEdge(int index);
    static int alu2duneEdge(int index);
 
    static int dune2aluVertex(int index);
    static int alu2duneVertex(int index);
 
    static int generic2aluFace ( const int index );
    static int alu2genericFace ( const int index );
 
    static int generic2aluVertex ( const int index );
    static int alu2genericVertex ( const int index );
 
    static int faceOrientation(int index);
 
    static int dune2aluFaceVertex(int face, int localVertex);
    static int alu2duneFaceVertex(int face, int localVertex);
 
    static int duneFaceTwist ( int face );
 
    static std::pair< int, int > duneEdgeMap ( int edge );
 
    static int faceVertex ( int face, int local );
 
  private:
    const static int dune2aluFace_[numFaces];
    const static int alu2duneFace_[numFaces];
 
    const static int dune2aluEdge_[numEdges];
    const static int alu2duneEdge_[numEdges];
 
    const static int dune2aluVertex_[numVertices];
    const static int alu2duneVertex_[numVertices];
 
    static const int generic2aluFace_[ numFaces ];
    static const int alu2genericFace_[ numFaces ];
 
    static const int generic2aluVertex_[ numVertices ];
    static const int alu2genericVertex_[ numVertices ];
 
    const static int faceOrientation_[numFaces];
 
    const static int dune2aluFaceVertex_[numFaces][numVerticesPerFace];
    const static int alu2duneFaceVertex_[numFaces][numVerticesPerFace];
 
    static const int duneFaceTwist_[ numFaces ];
 
    static const int duneEdgeMap_[ numEdges ][ 2 ];
 
    static const int faceVertex_[ numFaces ][ numVerticesPerFace ];
  };
 
  template <ALU3dGridElementType type>
  class FaceTopologyMapping {
  public:
    static int dune2aluVertex(int index);
    static int dune2aluVertex(int index, int twist);
    static int alu2duneVertex(int index);
    static int alu2duneVertex(int index, int twist);
    static int dune2aluEdge(int index);
    static int alu2duneEdge(int index);
    //  private:
    static int twist(int index, int faceTwist);
    static int invTwist(int index, int faceTwist);
 
    static int twistedDuneIndex( const int idx, const int twist );
 
    // for each aluTwist apply additional mapping
    static int aluTwistMap(const int aluTwist);
  private:
    const static int dune2aluVertex_[EntityCount<type>::numVerticesPerFace];
    const static int alu2duneVertex_[EntityCount<type>::numVerticesPerFace];
 
    const static int dune2aluEdge_[EntityCount<type>::numEdgesPerFace];
    const static int alu2duneEdge_[EntityCount<type>::numEdgesPerFace];
 
    const static int alu2duneTwist_[ 2 * EntityCount<type>::numVerticesPerFace ];
    const static int aluTwistMap_[ 2 * EntityCount<type>::numVerticesPerFace ];
  };
 
} // end namespace Dune
#endif