PollutantData.hpp

Go to the documentation of this file.

 
#ifndef OPENSWMM_ENGINE_POLLUTANT_DATA_HPP
#define OPENSWMM_ENGINE_POLLUTANT_DATA_HPP
 
#include <vector>
#include <string>
#include <cstdint>
 
namespace openswmm {
 
// ============================================================================
// Concentration units
// ============================================================================
 
enum class MassUnits : int8_t {
    MG_PER_L  = 0,  
    UG_PER_L  = 1,  
    COUNTS_PER_L = 2 
};
 
// ============================================================================
// PollutantData — pollutant definitions
// ============================================================================
 
struct PollutantData {
 
    // -----------------------------------------------------------------------
    // Pollutant definition arrays (indexed by pollutant index)
    // -----------------------------------------------------------------------
 
    std::vector<MassUnits>  units;
 
    std::vector<double>     mwt;
 
    std::vector<double>     k_decay;
 
    std::vector<double>     c_rain;
 
    std::vector<double>     c_gw;
 
    std::vector<double>     c_rdii;
 
    std::vector<double>     init_conc;
 
    std::vector<int>        co_pollut;
 
    std::vector<double>     co_frac;
 
    std::vector<bool>       snow_only;
 
    // -----------------------------------------------------------------------
    // Per-node quality state — flat 2D: [node * n_pollutants + pollutant]
    // -----------------------------------------------------------------------
 
    std::vector<double>     node_conc;
 
    std::vector<double>     node_conc_old;
 
    // -----------------------------------------------------------------------
    // Per-link quality state — flat 2D: [link * n_pollutants + pollutant]
    // -----------------------------------------------------------------------
 
    std::vector<double>     link_conc;
 
    std::vector<double>     link_conc_old;
 
    // -----------------------------------------------------------------------
    // Per-subcatch quality state — flat 2D: [subcatch * n_pollutants + pollutant]
    // -----------------------------------------------------------------------
 
    std::vector<double>     subcatch_conc;
 
    std::vector<double>     subcatch_conc_old;
 
    // -----------------------------------------------------------------------
    // Capacity management
    // -----------------------------------------------------------------------
 
    int n_pollutants() const noexcept { return static_cast<int>(units.size()); }
 
    void resize_pollutants(int n) {
        const auto un = static_cast<std::size_t>(n);
        units.assign(un, MassUnits::MG_PER_L);
        mwt.assign(un, 1.0);
        k_decay.assign(un, 0.0);
        c_rain.assign(un, 0.0);
        c_gw.assign(un, 0.0);
        c_rdii.assign(un, 0.0);
        init_conc.assign(un, 0.0);
        co_pollut.assign(un, -1);
        co_frac.assign(un, 0.0);
        snow_only.assign(un, false);
    }
 
    void resize_quality(int n_nodes, int n_links, int n_subcatches) {
        const auto np = static_cast<std::size_t>(n_pollutants());
        node_conc.assign(static_cast<std::size_t>(n_nodes) * np, 0.0);
        node_conc_old.assign(static_cast<std::size_t>(n_nodes) * np, 0.0);
        link_conc.assign(static_cast<std::size_t>(n_links) * np, 0.0);
        link_conc_old.assign(static_cast<std::size_t>(n_links) * np, 0.0);
        subcatch_conc.assign(static_cast<std::size_t>(n_subcatches) * np, 0.0);
        subcatch_conc_old.assign(static_cast<std::size_t>(n_subcatches) * np, 0.0);
    }
 
    void save_quality_state() noexcept {
        node_conc_old = node_conc;
        link_conc_old = link_conc;
        subcatch_conc_old = subcatch_conc;
    }
 
    void reset_quality_state() noexcept {
        std::fill(node_conc.begin(),    node_conc.end(),    0.0);
        std::fill(node_conc_old.begin(),node_conc_old.end(),0.0);
        std::fill(link_conc.begin(),    link_conc.end(),    0.0);
        std::fill(link_conc_old.begin(),link_conc_old.end(),0.0);
        std::fill(subcatch_conc.begin(),    subcatch_conc.end(),    0.0);
        std::fill(subcatch_conc_old.begin(),subcatch_conc_old.end(),0.0);
    }
};
 
} /* namespace openswmm */
 
#endif /* OPENSWMM_ENGINE_POLLUTANT_DATA_HPP */