AMReX_EB2_Graph.H

Go to the documentation of this file.

#ifndef AMREX_EB2_GRAPH_H_
#define AMREX_EB2_GRAPH_H_
#include <AMReX_Config.H>
 
#include <AMReX_Array.H>
#include <AMReX_Vector.H>
#include <AMReX_IntVect.H>
 
#include <unordered_map>
#include <utility>
#include <cstdint>
#include <limits>
 
namespace amrex::EB2 {
 
struct Cell
{
    RealArray m_volfrac;
    // and more such as boundary centroid and normal
};
 
struct Face
{
    RealArray m_facecent;
    Real m_areafrac;
};
 
using VertexID = std::pair<IntVect,int>;
 
struct Edge
{
    Face m_face;
    VertexID m_vid;
};
 
using Type_t = uint32_t;
 
namespace Type {
    static constexpr Type_t regular=0x0;
    static constexpr Type_t covered=0x1;
    static constexpr Type_t irregular=0x2;
    static constexpr Type_t unknown=0x3;
}
 
static constexpr Real no_intercept = std::numeric_limits<Real>::max();
 
struct Neighbor
{
    void setXLo (Type_t t) { flag &= zero_xlo_mask; flag |= t; }
    void setXHi (Type_t t) { flag &= zero_xhi_mask; flag |= (t << 2); }
    void setYLo (Type_t t) { flag &= zero_ylo_mask; flag |= (t << 4); }
    void setYHi (Type_t t) { flag &= zero_yhi_mask; flag |= (t << 6); }
    void setZLo (Type_t t) { flag &= zero_zlo_mask; flag |= (t << 8); }
    void setZHi (Type_t t) { flag &= zero_zhi_mask; flag |= (t << 10); }
 
    [[nodiscard]] Type_t getXLo () const { return (flag & one_xlo_mask); }
    [[nodiscard]] Type_t getXHi () const { return (flag & one_xhi_mask) >> 2; }
    [[nodiscard]] Type_t getYLo () const { return (flag & one_ylo_mask) >> 4; }
    [[nodiscard]] Type_t getYHi () const { return (flag & one_yhi_mask) >> 6; }
    [[nodiscard]] Type_t getZLo () const { return (flag & one_zlo_mask) >> 8; }
    [[nodiscard]] Type_t getZHi () const { return (flag & one_zhi_mask) >> 10; }
 
private:
 
    Type_t flag = 0;
 
    static constexpr Type_t  one_xlo_mask =  0x3;
    static constexpr Type_t zero_xlo_mask = ~0x3;
    static constexpr Type_t  one_xhi_mask =  0xC;
    static constexpr Type_t zero_xhi_mask = ~0xC;
    static constexpr Type_t  one_ylo_mask =  0x30;
    static constexpr Type_t zero_ylo_mask = ~0x30;
    static constexpr Type_t  one_yhi_mask =  0xC0;
    static constexpr Type_t zero_yhi_mask = ~0xC0;
    static constexpr Type_t  one_zlo_mask =  0x300;
    static constexpr Type_t zero_zlo_mask = ~0x300;
    static constexpr Type_t  one_zhi_mask =  0xC00;
    static constexpr Type_t zero_zhi_mask = ~0xC00;
};
 
struct Vertex
{
    Cell m_cell;
    Neighbor m_nghb;
    Vector<Edge> m_edges;
};
 
using Graph = std::unordered_map<IntVect, Vector<Vertex>, IntVect::hasher>;
 
}
 
#endif