sizefield.hpp

Go to the documentation of this file.

 
#include "meshkit/detail/planmodel.hpp"
#ifndef MESHKIT_NATIVE_SIZEFIELD_HPP
#define MESHKIT_NATIVE_SIZEFIELD_HPP
 
namespace meshkit::plan
{
 
struct SizeFieldOptions {
    bool footprintBounds = false; 
    double nearSize = 0.0;
 
    double gradation = 0.25;
 
    double areaFloor = 0.0;
 
    std::int64_t maxGridCells = 4'000'000;
 
    // ── Overhaul additions (MESH_OVERHAUL_PLAN_2026-09-29.md §3 Stage 2) ──
    double maxSize = 0.0;
    std::function<double(double, double)> terrainSizeAt;
    struct Region {
        Polygon ring;
        double h = 0.0;
    };
    Vector<Region> regions;
    Vector<Region>
        terrainRegions; 
 
    Vector<Vector<Point2>> steps;
};
 
class SizeField
{
public:
    // Conservative bound on every bilinear target over this rectangle. Zero
    // means unbounded/invalid field. Includes all interpolation support nodes.
    double minimumSizeInBox(const Box2 &) const;
    double minimumTerrainSizeInBox(const Box2 &)
        const; 
 
    bool build(const Box2 &bbox, const Vector<ConstraintSegment> &segs,
               const Vector<Vector<Point2>> &rings, const Vector<SteinerPoint> &pts,
               const SizeFieldOptions &opt);
 
    [[nodiscard]] double targetAreaAt(double x, double y) const;
    [[nodiscard]] double sizeAt(double x, double y) const;
    [[nodiscard]] bool isValid() const { return m_cols > 0 && m_rows > 0; }
    [[nodiscard]] bool hasSizeBound() const { return isValid() && sizeAt(m_x0, m_y0) > 0.0; }
    [[nodiscard]] double pitch() const { return m_pitch; }
    [[nodiscard]] int cols() const { return m_cols; }
    [[nodiscard]] int rows() const { return m_rows; }
 
    [[nodiscard]] double distanceAt(double x, double y) const;
 
private:
    [[nodiscard]] double cellDist(int cx, int cy) const;
    [[nodiscard]] double cellSize(int cx, int cy) const;
    void buildSizeGrid(const SizeFieldOptions &opt, const Vector<float> &stepDist);
    void stampSeedPoint(const Point2 &p, Vector<float> &into);
    void stampSeedSegment(const Point2 &a, const Point2 &b, Vector<float> &into);
 
    Vector<float> m_dist; 
    Vector<float> m_h;    
    struct MinLevel {
        int columns, rows;
        Vector<float> values;
    };
    Vector<MinLevel> m_minLevels;
    double minimumSizeInBoxImpl(const Box2 &, bool includeTerrain) const;
    struct RegionBound {
        Box2 box;
        double size = 0;
        int end = 0;
        bool terrainOnly = false;
    };
    Vector<RegionBound> m_regionBounds; 
    int m_cols = 0, m_rows = 0;
    double m_x0 = 0.0, m_y0 = 0.0; 
    double m_pitch = 0.0;
    double m_near = 0.0, m_g = 0.0, m_floor = 0.0;
};
 
} // namespace meshkit::plan
 
#endif // MESHKIT_NATIVE_SIZEFIELD_HPP