pslgminsize.h

Go to the documentation of this file.

 
#ifndef OPENSWMMVIS_MESH_PSLGMINSIZE_H
#define OPENSWMMVIS_MESH_PSLGMINSIZE_H
 
#include "mesh/meshgenerator.h"
 
#include <QPointF>
#include <QPolygonF>
#include <QString>
#include <QVector>
 
#include <functional>
 
namespace mesh {
namespace pslg {
 
struct MinSizePolicy
{
    double minCellSize = 0.0;
 
    double weldRadius    = -1.0;  
    double minSegmentLen = -1.0;  
    double maxDeviation  = -1.0;  
 
    double trimAngleDeg = 30.0;
 
    bool trimAtTaggedNodes = false;
 
    bool dropSubScaleHoles = true;
 
    bool collapseSubScalePaths = true;
 
    bool ringsReadOnly = false;
 
    bool allowIdentityMerge = false;
 
    double identityMergeRadius = -1.0;
 
    int maxIterations = 3;        
    int maxResiduals  = 200;      
 
    void resolveDefaults();
 
    [[nodiscard]] bool enabled() const { return minCellSize > 0.0; }
 
    [[nodiscard]] double minTriangleArea() const;
 
    [[nodiscard]] double refinementAreaCap(double uniformCap) const;
};
 
enum class ViolationCause
{
    ShortSegment,   
    CloseFeatures,  
    SmallAngle,     
    SubScaleRing,   
    IdentityMerged  
};
 
[[nodiscard]] QString violationCauseName(ViolationCause c);
 
struct Violation
{
    QPointF        xy;
    double         lfs = 0.0;   
    ViolationCause cause = ViolationCause::CloseFeatures;
    QString        tagA;        
    QString        tagB;        
};
 
struct ConditionReport
{
    int verticesWelded    = 0;
    int segmentsSplit     = 0;   
    int pathsCollapsed    = 0;
    int pathsShortened    = 0;   
    int shortEdgesKept    = 0;   
    int holesDropped      = 0;
    int cornersTrimmed    = 0;
    int cornersSkipped    = 0;   
    int crossingsBefore   = 0;
    int crossingsAfter    = 0;
    int crossingsRepaired = 0;
    int iterationsUsed    = 0;
 
    int duplicateSegments = 0;
 
    int identitiesMerged = 0;
 
    bool   weldRadiusDominatesDomain = false;
    double domainNarrowestExtent     = 0.0;
 
    double maxDisplacement  = 0.0;
    double domainAreaBefore = 0.0;
    double domainAreaAfter  = 0.0;
    double predictedMinLfs  = 0.0;  
 
    bool conditioningAbandoned = false;
 
    QString abandonReason;
 
    QVector<Violation> residuals;   
 
    [[nodiscard]] QString summary() const;
};
 
[[nodiscard]] QVector<Violation> analyseLocalFeatureSize(
    const QVector<QPolygonF>          &domains,
    const QVector<QVector<QPointF>>   &holeRings,
    const QVector<ConstraintSegment>  &segs,
    const QVector<SteinerPoint>       &pts,
    double h,
    int maxReported = 200);
 
bool conditionMinSize(QVector<QPolygonF>         *domains,
                      QVector<QVector<QPointF>>  *holeRings,
                      QVector<ConstraintSegment> *segs,
                      QVector<SteinerPoint>      *pts,
                      const MinSizePolicy        &policy,
                      ConditionReport            *report,
                      const std::function<bool()> &isCancelled = {});
 
} // namespace pslg
} // namespace mesh
 
#endif // OPENSWMMVIS_MESH_PSLGMINSIZE_H