InertialEdges.hpp

Go to the documentation of this file.

 
#ifndef OPENSWMM_ENGINE_2D_INERTIAL_EDGES_HPP
#define OPENSWMM_ENGINE_2D_INERTIAL_EDGES_HPP
 
#include <vector>
 
namespace openswmm::twoD {
 
struct MeshData;

struct InertialEdges {
    int ne = 0;                       
 
    // Per-edge arrays (size ne), oriented cL→cR (cL = min(t,nbr), cR = max).
    std::vector<int>    cL, cR;       
    std::vector<double> xi;           
    std::vector<double> inv_dx;       
    std::vector<double> zface;        
    std::vector<int>    slotL, slotR; 
 
    // Per-cell CSR incidence for the conservative continuity gather. For cell i,
    // the incident edges are cell_edge[cell_ptr[i] .. cell_ptr[i+1]) with sign
    // cell_sign (+1 if i==cL, −1 if i==cR). Then
    //   dV_i/dt (flux part) = − Σ cell_sign · q[edge] · ξ[edge].
    std::vector<int>    cell_ptr;     
    std::vector<int>    cell_edge;    
    std::vector<double> cell_sign;    

    void build(const MeshData& mesh);
 
    bool empty() const noexcept { return ne == 0; }
};
 
} // namespace openswmm::twoD
 
#endif // OPENSWMM_ENGINE_2D_INERTIAL_EDGES_HPP