meshcellstats.h

Go to the documentation of this file.

 
#ifndef OPENSWMMVIS_MESH_MESHCELLSTATS_H
#define OPENSWMMVIS_MESH_MESHCELLSTATS_H
 
#include "meshresult.h"
 
#include <QHash>
#include <QPair>
 
namespace mesh {
 
struct CellAreaStats
{
    double min    = 0.0;
    double max    = 0.0;
    double mean   = 0.0;
    double median = 0.0;
    int    count  = 0;    
};
 
[[nodiscard]] double triangleArea(const MeshResult &mesh, int tri);
 
[[nodiscard]] QPointF cellCentroid(const MeshResult &mesh, int cell);
 
struct QuadStats
{
    int    count           = 0;    
    double minAngleDeg     = 360.0;
    double maxAngleDeg     = 0.0;  
    double maxNonPlanarity = 0.0;
    // QUAD_MESHING_REDESIGN_PLAN_2026-09-06.md §5 — shape metrics
    // (mesh::quadQuality). All 0 when count == 0.
    double minScaledJacobian    = 1.0;  
    double medianRectangularity = 0.0;  
    double maxAspect            = 1.0;  
    int    nonConvex            = 0;    
    int    irregularVertices    = 0;    
    int    sjHistogram[10] = {0,0,0,0,0,0,0,0,0,0};
};
 
[[nodiscard]] QuadStats computeQuadStats(const MeshResult &mesh);
 
[[nodiscard]] CellAreaStats computeCellAreaStats(const MeshResult &mesh);
 
using EdgeVertexPair = QPair<int, int>;
 
[[nodiscard]] QHash<EdgeVertexPair, QVector<int>>
buildEdgeTriangles(const MeshResult &mesh);
 
struct CflStats
{
    double min    = 0.0;   
    double p10    = 0.0;
    double p50    = 0.0;
    double p90    = 0.0;
    double max    = 0.0;
    double ratio  = 0.0;   
 
    int impliedTiers = 0;
 
    double ltsWork = 0.0;
 
    int count          = 0;  
    int interiorFaces  = 0;
    int isolatedCells  = 0;  
};
 
[[nodiscard]] QVector<double> computeCellLchar(const MeshResult &mesh);
 
[[nodiscard]] CflStats computeCflStats(const MeshResult &mesh, int ltsTiers = 4);
 
} // namespace mesh
 
#endif // OPENSWMMVIS_MESH_MESHCELLSTATS_H