GMSHAdaptiveMeshDensity.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#pragma once
 
#ifndef NDEBUG
#include <string>
#endif
 
#include "GMSHMeshDensityStrategy.h"
 
namespace GeoLib
{
class Point;
template <typename T>
class QuadTree;
#ifndef NDEBUG
class GEOObjects;
#endif
}  // namespace GeoLib
 
namespace FileIO
{
namespace GMSH
{
 
class GMSHAdaptiveMeshDensity : public GMSHMeshDensityStrategy
{
public:
    GMSHAdaptiveMeshDensity(double pnt_density,
                            double station_density,
                            std::size_t max_pnts_per_leaf);
    ~GMSHAdaptiveMeshDensity() override;
    void initialize(std::vector<GeoLib::Point const*> const& pnts) override;
    double getMeshDensityAtPoint(GeoLib::Point const* const pnt) const override;
    void addPoints(std::vector<GeoLib::Point const*> const& pnts);
    double getMeshDensityAtStation(
        GeoLib::Point const* const /*unused*/) const override;
    void getSteinerPoints(std::vector<GeoLib::Point*>& pnts,
                          std::size_t additional_levels = 0) const;
#ifndef NDEBUG
    std::string getQuadTreeGeometry(GeoLib::GEOObjects& geo_objs) const;
#endif
 
private:
    double _pnt_density;
    double _station_density;
    std::size_t _max_pnts_per_leaf;
    GeoLib::QuadTree<GeoLib::Point>* _quad_tree = nullptr;
};
 
}  // namespace GMSH
}  // end namespace FileIO