FileIO::GMSH::GMSHAdaptiveMeshDensity Class Reference

Detailed Description

Definition at line 34 of file GMSHAdaptiveMeshDensity.h.

#include <GMSHAdaptiveMeshDensity.h>

Inheritance diagram for FileIO::GMSH::GMSHAdaptiveMeshDensity:

Collaboration diagram for FileIO::GMSH::GMSHAdaptiveMeshDensity:

Public Member Functions
	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) const override
void	getSteinerPoints (std::vector< GeoLib::Point * > &pnts, std::size_t additional_levels=0) const
std::string	getQuadTreeGeometry (GeoLib::GEOObjects &geo_objs) const
Public Member Functions inherited from FileIO::GMSH::GMSHMeshDensityStrategy
virtual	~GMSHMeshDensityStrategy ()=default

Private Attributes
double	_pnt_density
double	_station_density
std::size_t	_max_pnts_per_leaf
GeoLib::QuadTree< GeoLib::Point > *	_quad_tree

Constructor & Destructor Documentation

GMSHAdaptiveMeshDensity()

FileIO::GMSH::GMSHAdaptiveMeshDensity::GMSHAdaptiveMeshDensity	(	double	pnt_density,
		double	station_density,
		std::size_t	max_pnts_per_leaf )

Definition at line 32 of file GMSHAdaptiveMeshDensity.cpp.

    : _pnt_density(pnt_density),
      _station_density(station_density),
      _max_pnts_per_leaf(max_pnts_per_leaf),
      _quad_tree(nullptr)
{
}

~GMSHAdaptiveMeshDensity()

FileIO::GMSH::GMSHAdaptiveMeshDensity::~GMSHAdaptiveMeshDensity ( )

override

Definition at line 42 of file GMSHAdaptiveMeshDensity.cpp.

{
    delete _quad_tree;
}

References _quad_tree.

Member Function Documentation

addPoints()

void FileIO::GMSH::GMSHAdaptiveMeshDensity::addPoints

(

std::vector< GeoLib::Point const * > const &

pnts

)

Definition at line 89 of file GMSHAdaptiveMeshDensity.cpp.

{
    // *** QuadTree - insert points
    const std::size_t n_pnts(pnts.size());
    DBUG(
        "GMSHAdaptiveMeshDensity::addPoints(): Inserting {:d} points into "
        "quadtree.",
        n_pnts);
    for (std::size_t k(0); k < n_pnts; k++)
    {
        _quad_tree->addPoint(pnts[k]);
    }
    DBUG("GMSHAdaptiveMeshDensity::addPoints(): \tok.");
    _quad_tree->balance();
}

References _quad_tree, GeoLib::QuadTree< POINT >::addPoint(), GeoLib::QuadTree< POINT >::balance(), and DBUG().

Referenced by initialize().

getMeshDensityAtPoint()

double FileIO::GMSH::GMSHAdaptiveMeshDensity::getMeshDensityAtPoint ( GeoLib::Point const *const pnt ) const

overridevirtual

Implements FileIO::GMSH::GMSHMeshDensityStrategy.

Definition at line 106 of file GMSHAdaptiveMeshDensity.cpp.

{
    GeoLib::Point ll;
    GeoLib::Point ur;
    _quad_tree->getLeaf(*pnt, ll, ur);
    return _pnt_density * (ur[0] - ll[0]);
}

References _pnt_density, _quad_tree, and GeoLib::QuadTree< POINT >::getLeaf().

getMeshDensityAtStation()

double FileIO::GMSH::GMSHAdaptiveMeshDensity::getMeshDensityAtStation ( GeoLib::Point const * const pnt ) const

overridevirtual

Implements FileIO::GMSH::GMSHMeshDensityStrategy.

Definition at line 115 of file GMSHAdaptiveMeshDensity.cpp.

{
    GeoLib::Point ll;
    GeoLib::Point ur;
    _quad_tree->getLeaf(*pnt, ll, ur);
    return _station_density * (ur[0] - ll[0]);
}

References _quad_tree, _station_density, and GeoLib::QuadTree< POINT >::getLeaf().

getQuadTreeGeometry()

std::string FileIO::GMSH::GMSHAdaptiveMeshDensity::getQuadTreeGeometry

(

GeoLib::GEOObjects &

geo_objs

)

const

Definition at line 166 of file GMSHAdaptiveMeshDensity.cpp.

{
    std::list<GeoLib::QuadTree<GeoLib::Point>*> leaf_list;
    _quad_tree->getLeafs(leaf_list);
 
    std::string quad_tree_geo("QuadTree");
    {
        std::vector<GeoLib::Point*> points{};
        for (auto const leaf : leaf_list)
        {
            // fetch corner points from leaf
            GeoLib::Point ll;
            GeoLib::Point ur;
            leaf->getSquarePoints(ll, ur);
            std::size_t const pnt_offset(points.size());
            points.push_back(new GeoLib::Point(ll, pnt_offset));
            points.push_back(
                new GeoLib::Point(ur[0], ll[1], 0.0, pnt_offset + 1));
            points.push_back(new GeoLib::Point(ur, pnt_offset + 2));
            points.push_back(
                new GeoLib::Point(ll[0], ur[1], 0.0, pnt_offset + 3));
        }
        geo_objs.addPointVec(std::move(points), quad_tree_geo,
                             GeoLib::PointVec::NameIdMap{});
    }
    auto& points = geo_objs.getPointVecObj(quad_tree_geo)->getVector();
 
    std::vector<GeoLib::Polyline*> polylines{};
    for (std::size_t l = 0; l < leaf_list.size(); ++l)
    {
        auto* polyline = new GeoLib::Polyline(points);
        for (std::size_t p = 0; p < 4; ++p)
        {
            polyline->addPoint(4 * l + p);
        }
        polyline->closePolyline();
        polylines.push_back(polyline);
    }
    geo_objs.addPolylineVec(std::move(polylines), quad_tree_geo,
                            GeoLib::PolylineVec::NameIdMap{});
 
    return quad_tree_geo;
}

References _quad_tree, GeoLib::GEOObjects::addPointVec(), GeoLib::GEOObjects::addPolylineVec(), GeoLib::QuadTree< POINT >::getLeafs(), GeoLib::GEOObjects::getPointVecObj(), and GeoLib::TemplateVec< T >::getVector().

getSteinerPoints()

void FileIO::GMSH::GMSHAdaptiveMeshDensity::getSteinerPoints	(	std::vector< GeoLib::Point * > &	pnts,
		std::size_t	additional_levels = 0 ) const

Definition at line 124 of file GMSHAdaptiveMeshDensity.cpp.

{
    // get Steiner points
    std::size_t max_depth(0);
    _quad_tree->getMaxDepth(max_depth);
 
    std::list<GeoLib::QuadTree<GeoLib::Point>*> leaf_list;
    _quad_tree->getLeafs(leaf_list);
 
    for (std::list<GeoLib::QuadTree<GeoLib::Point>*>::const_iterator it(
             leaf_list.begin());
         it != leaf_list.end();
         ++it)
    {
        if ((*it)->getPoints().empty())
        {
            // compute point from square
            GeoLib::Point ll;
            GeoLib::Point ur;
            (*it)->getSquarePoints(ll, ur);
            if ((*it)->getDepth() + additional_levels > max_depth)
            {
                additional_levels = max_depth - (*it)->getDepth();
            }
            const std::size_t n_pnts_per_quad_dim = static_cast<std::size_t>(1)
                                                    << additional_levels;
            const double delta((ur[0] - ll[0]) / (2 * n_pnts_per_quad_dim));
            for (std::size_t i(0); i < n_pnts_per_quad_dim; i++)
            {
                for (std::size_t j(0); j < n_pnts_per_quad_dim; j++)
                {
                    pnts.push_back(new GeoLib::Point(
                        ll[0] + (2 * i + 1) * delta,
                        ll[1] + (2 * j + 1) * delta, 0.0, pnts.size()));
                }
            }
        }
    }
}

References _quad_tree, GeoLib::QuadTree< POINT >::getLeafs(), and GeoLib::QuadTree< POINT >::getMaxDepth().

initialize()

void FileIO::GMSH::GMSHAdaptiveMeshDensity::initialize ( std::vector< GeoLib::Point const * > const & pnts )

overridevirtual

Implements FileIO::GMSH::GMSHMeshDensityStrategy.

Definition at line 47 of file GMSHAdaptiveMeshDensity.cpp.

{
    // *** QuadTree - determining bounding box
    DBUG(
        "GMSHAdaptiveMeshDensity::init(): computing axis aligned bounding box "
        "(2D) for quadtree.");
 
    GeoLib::Point min(*pnts[0]);
    GeoLib::Point max(*pnts[0]);
    std::size_t n_pnts(pnts.size());
    for (std::size_t k(1); k < n_pnts; k++)
    {
        for (std::size_t j(0); j < 2; j++)
        {
            if ((*(pnts[k]))[j] < min[j])
            {
                min[j] = (*(pnts[k]))[j];
            }
        }
        for (std::size_t j(0); j < 2; j++)
        {
            if ((*(pnts[k]))[j] > max[j])
            {
                max[j] = (*(pnts[k]))[j];
            }
        }
    }
    min[2] = 0.0;
    max[2] = 0.0;
    DBUG("GMSHAdaptiveMeshDensity::init(): \tok");
 
    // *** QuadTree - create object
    DBUG("GMSHAdaptiveMeshDensity::init(): Creating quadtree.");
    _quad_tree =
        new GeoLib::QuadTree<GeoLib::Point>(min, max, _max_pnts_per_leaf);
    DBUG("GMSHAdaptiveMeshDensity::init(): \tok.");
 
    // *** QuadTree - insert points
    addPoints(pnts);
}

References _max_pnts_per_leaf, _quad_tree, addPoints(), and DBUG().

Member Data Documentation

_max_pnts_per_leaf

std::size_t FileIO::GMSH::GMSHAdaptiveMeshDensity::_max_pnts_per_leaf

private

Definition at line 55 of file GMSHAdaptiveMeshDensity.h.

Referenced by initialize().

_pnt_density

double FileIO::GMSH::GMSHAdaptiveMeshDensity::_pnt_density

private

Definition at line 53 of file GMSHAdaptiveMeshDensity.h.

Referenced by getMeshDensityAtPoint().

_quad_tree

GeoLib::QuadTree<GeoLib::Point>* FileIO::GMSH::GMSHAdaptiveMeshDensity::_quad_tree

private

Definition at line 56 of file GMSHAdaptiveMeshDensity.h.

Referenced by ~GMSHAdaptiveMeshDensity(), addPoints(), getMeshDensityAtPoint(), getMeshDensityAtStation(), getQuadTreeGeometry(), getSteinerPoints(), and initialize().

_station_density

double FileIO::GMSH::GMSHAdaptiveMeshDensity::_station_density

private

Definition at line 54 of file GMSHAdaptiveMeshDensity.h.

Referenced by getMeshDensityAtStation().

---

The documentation for this class was generated from the following files:

• Applications/FileIO/Gmsh/GMSHAdaptiveMeshDensity.h
• Applications/FileIO/Gmsh/GMSHAdaptiveMeshDensity.cpp