coordinateChange.hh

Go to the documentation of this file.

#ifndef COORDINATECHANGE_HH
#define COORDINATECHANGE_HH
 
#include <array>
#include <unordered_map>
#include <utility>
 
#include "basics.hh"
 
namespace concepts {
    
// ******************************* Coordinate Transformation Isomorphisms **
 
    template<uint pdim>
    class QuadCoordinateChange;
 
    template<uint pdim>
    class EdgeCoordinateChange;
 
 
    // ***************************************** Boundary elements of Hexahedron **
 
    template<>
    class QuadCoordinateChange<3> {
    public:
        static std::pair<uint,uint> TopToGeo(uint faceIndex, std::pair<Z2,Z4> topologicalOrientation) {
            return TopToGeo_[faceIndex].at(topologicalOrientation);
        }
        static std::pair<uint,uint> TopToCon(uint faceIndex, std::pair<Z2,Z4> topologicalOrientation) {
            return TopToCon_[faceIndex].at(topologicalOrientation);
        }
        static std::pair<uint,uint> GeoToCon(uint faceIndex) {
            return GeoToCon_[faceIndex];
        }
 
 
    private:
      
      struct QuadHash {
        std::size_t operator()(const std::pair<Z2,Z4>& key) const {
          return key.first*4+key.second;
        }
      };
 
      const static std::array< std::unordered_map< std::pair<Z2,Z4>, std::pair<uint,uint>, QuadHash >,6  > TopToGeo_;
      const static std::array< std::unordered_map< std::pair<Z2,Z4>, std::pair<uint,uint>, QuadHash >,6  > TopToCon_ ;
      const static std::array< std::pair<uint,uint>,6 > GeoToCon_;
 
    };
 
    template<>
    class EdgeCoordinateChange<3> {
    public:
        static uint GeoToCon(uint edgeIndex) {
            return GeoToCon_[edgeIndex];
        }
 
    private:
      const static std::array<uint,12> GeoToCon_;
    };
 
 
    // ***************************************** Boundary elements of Quad **
 
    template<>
    class EdgeCoordinateChange<2> {
    public:
        static uint TopToGeo(uint edgeIndex, Z2 topologicalOrientation) {
            return TopToGeo_[edgeIndex].at(topologicalOrientation);
        }
        static uint TopToCon(uint edgeIndex, Z2 topologicalOrientation) {
            return TopToCon_[edgeIndex].at(topologicalOrientation);
        }
        static uint GeoToCon(uint edgeIndex) {
            return GeoToCon_[edgeIndex];
        }
 
    private:
      
      struct EdgeHash {
        std::size_t operator()(const Z2& key) const {
          return key;
        }
      };
 
      const static std::array< std::unordered_map< Z2, uint, EdgeHash >, 4 > TopToGeo_;
      const static std::array< std::unordered_map< Z2, uint, EdgeHash >, 4 > TopToCon_;
      const static std::array<uint,4> GeoToCon_;
    };
 
}
 
#endif // COORDINATECHANGE_HH