quad.h

#ifndef QUAD_H
#define QUAD_H
 
#include <vector>
#include <array>
#include <iostream>
#include <cmath>
#include <unordered_set>
#include <unordered_map>
#include <fstream>
#include <sstream>
#include <string>
#include <algorithm>
 
 
#include "geometry/point.h"
#include "geometry/edge.h"
 
namespace Enmesh {

struct Quad {
 
    static constexpr size_t elementType = 3; 
    static constexpr size_t numVertices = 4; 
 
    std::array<size_t, 4> v;   
 

    bool containsEdge(size_t v1, size_t v2) const;
 

    bool operator==(const Quad& other) const {
        std::array<size_t, 4> sortedV1 = v;
        std::array<size_t, 4> sortedV2 = other.v;
        std::sort(sortedV1.begin(), sortedV1.end());
        std::sort(sortedV2.begin(), sortedV2.end());
        return sortedV1 == sortedV2;
    }
 
};
 

struct QuadHash {
    size_t operator()(const Quad& quad) const {
        // Combine the hashes of the four vertex indices to create a unique hash for the quad
        size_t h = std::hash<size_t>{}(quad.v[0]);
        h ^= std::hash<size_t>{}(quad.v[1]) + 0x9e3779b9 + (h << 6) + (h >> 2);  // Combine with the second vertex index
        h ^= std::hash<size_t>{}(quad.v[2]) + 0x9e3779b9 + (h << 6) + (h >> 2);  // Combine with the third vertex index
        h ^= std::hash<size_t>{}(quad.v[3]) + 0x9e3779b9 + (h << 6) + (h >> 2);  // Combine with the fourth vertex index
        return h;
    }
};
 
} // namespace Enmesh
 
#endif