profilebuilder.h

Go to the documentation of this file.

 
#ifndef PROFILE_BUILDER_H
#define PROFILE_BUILDER_H
 
#include <QDateTime>
#include <QPointF>
#include <QString>
#include <QVector>
 
#include <limits>
 
namespace ProfileBuilder
{
 
inline constexpr double kGravityFps2 = 32.174;   
inline constexpr double kGravityMps2 = 9.80665;  
enum class NodeKind
{
    Junction = 0,
    Outfall,
    Storage,
    Divider,
    VirtualJunction,
};
 
enum class LinkKind
{
    Conduit = 0,
    Pump,
    Orifice,
    Weir,
    Outlet,
};
 
enum class OutputKind
{
    HGL = 0,
    EGL,
    WaterSurface,
    MaxHGL,
    MaxEGL,
    MinHGL,
    MinEGL,
    MaxWaterSurface,
    MinWaterSurface,
};
 
inline constexpr double kNoBearing = 1.0e9;
 
struct BranchLink
{
    QString  name;
    LinkKind kind       = LinkKind::Conduit;
    double   invertElev = 0.0;   
    double   maxDepth   = 0.0;   
    bool     openTop    = false; 
    double   bearingRad = 0.0;   
    bool     intoNode   = false; 
    bool     onPath     = false; 
};
 
struct NodeStatic
{
    QString  name;
    double   invertElev      = 0.0;
    double   maxDepth        = 0.0;
    double   surchargeDepth  = 0.0;   
    NodeKind kind            = NodeKind::Junction;
    bool     isInlet         = false;
    QVector<BranchLink> branches;
};
 
struct LinkStatic
{
    QString  name;
    double   length      = 0.0;
    double   offset1     = 0.0;   
    double   offset2     = 0.0;   
    double   maxDepth    = 0.0;   
    double   crestHeight = 0.0;   
    bool     openTop     = false; 
    bool     isStreet    = false; 
    LinkKind kind        = LinkKind::Conduit;
    bool     reversed    = false; 
};
 
struct PathStatic
{
    QVector<NodeStatic> nodes;
    QVector<LinkStatic> links;
    QVector<double>     chainage;
    QVector<QPointF>    terrainSamples;
};
 
struct SourceSeries
{
    QString             sourceId;
    QDateTime           startTime;
    int                 reportStepSec = 0;
    int                 periodCount   = 0;
 
    QVector<QVector<float>>  nodeHead;
    QVector<QVector<float>>  nodeDepth;
    QVector<QVector<float>>  linkVelocity;
};
 
struct SourceDerived
{
    QVector<QVector<double>> hglByPeriod;
    QVector<QVector<double>> eglByPeriod;
    QVector<QVector<double>> waterSurfaceByPeriod;
    QVector<double>          minHgl;
    QVector<double>          maxHgl;
    QVector<double>          minEgl;
    QVector<double>          maxEgl;
    QVector<double>          minWaterSurface;
    QVector<double>          maxWaterSurface;
};
 
struct Diagnostic
{
    QString error;
};
 
// --------------------------------------------------------------------------
// Static-topology utilities
// --------------------------------------------------------------------------
 
[[nodiscard]] QVector<double> computeChainage(const QVector<LinkStatic> &links);
 
[[nodiscard]] double groundElev(const NodeStatic &n);
 
[[nodiscard]] double conduitInvertUpstream(const PathStatic &path, int linkIdx);
 
[[nodiscard]] double conduitInvertDownstream(const PathStatic &path, int linkIdx);
 
// --------------------------------------------------------------------------
// Per-period EGL helper (exposed for callers that want a point query rather
// than a full envelope compute, e.g. tooltip / hover read-outs).
// --------------------------------------------------------------------------
 
[[nodiscard]] double velocityHead(const PathStatic &path,
                                  const SourceSeries &src,
                                  int nodeIdx,
                                  int period,
                                  double gravity);
 
[[nodiscard]] double eglAtPeriod(const PathStatic &path,
                                 const SourceSeries &src,
                                 int nodeIdx,
                                 int period,
                                 double gravity);
 
// --------------------------------------------------------------------------
// Validation + full compute
// --------------------------------------------------------------------------
 
[[nodiscard]] Diagnostic validate(const PathStatic &path,
                                  const SourceSeries &src);
 
[[nodiscard]] SourceDerived compute(const PathStatic &path,
                                    const SourceSeries &src,
                                    double gravity);
 
bool appendPeriods(const PathStatic &path,
                   const SourceSeries &src,
                   double gravity,
                   SourceDerived &derived);
 
void accumulatePeriod(const PathStatic &path,
                      const SourceSeries &src,
                      int p,
                      double gravity,
                      SourceDerived &out);
 
} // namespace ProfileBuilder
 
#endif // PROFILE_BUILDER_H