GMSHPolygonTree.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#pragma once
 
#include <vector>
#include <string>
 
#include "GMSHMeshDensityStrategy.h"
#include "GMSHPoint.h"
#include "GMSHLine.h"
 
#include "GeoLib/SimplePolygonTree.h"
 
namespace GeoLib
{
class GEOObjects;
class Polygon;
class PolylineWithSegmentMarker;
class PolygonWithSegmentMarker;
}
 
namespace FileIO
{
namespace GMSH
{
 
class GMSHPolygonTree: public GeoLib::SimplePolygonTree {
public:
    GMSHPolygonTree(GeoLib::PolygonWithSegmentMarker* polygon,
                    GMSHPolygonTree* parent, GeoLib::GEOObjects& geo_objs,
                    std::string const& geo_name,
                    GMSHMeshDensityStrategy& mesh_density_strategy);
    ~GMSHPolygonTree() override;

    void markSharedSegments();

    bool insertStation(GeoLib::Point const* station);
    void insertPolyline(GeoLib::PolylineWithSegmentMarker * ply);

    void initMeshDensityStrategy();

    void createGMSHPoints(std::vector<GMSHPoint*> & gmsh_pnts) const;
 
    virtual void writeLineLoop(std::size_t& line_offset,
                               std::size_t& sfc_offset, std::ostream& out,
                               bool const write_physical) const;
    void writeSubPolygonsAsLineConstraints(std::size_t &line_offset, std::size_t sfc_number, std::ostream& out) const;
    virtual void writeLineConstraints(std::size_t &line_offset, std::size_t sfc_number, std::ostream& out) const;
    void writeStations(std::size_t & pnt_id_offset, std::size_t sfc_number, std::ostream& out) const;
    void writeAdditionalPointData(std::size_t & pnt_id_offset, std::size_t sfc_number, std::ostream& out) const;
 
private:
    void getPointsFromSubPolygons(
        std::vector<GeoLib::Point const*>& pnts) const;
    void getStationsInsideSubPolygons(
        std::vector<GeoLib::Point const*>& stations) const;
    void checkIntersectionsSegmentExistingPolylines(
        GeoLib::PolylineWithSegmentMarker* ply,
        GeoLib::Polyline::SegmentIterator const& seg_it);
 
    GeoLib::GEOObjects & _geo_objs;
    std::string const& _geo_name;
    std::vector<GeoLib::Point const*> _stations;
    std::vector<GeoLib::PolylineWithSegmentMarker*> _plys;
    std::vector<GMSHLine*> _gmsh_lines_for_constraints;
 
    GMSHMeshDensityStrategy& _mesh_density_strategy;
};
 
}  // end namespace GMSH
}  // end namespace FileIO