RadiusEdgeRatioMetric.cpp

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#include "RadiusEdgeRatioMetric.h"
 
#include "GeoLib/MinimalBoundingSphere.h"
#include "MeshLib/Node.h"
 
namespace MeshToolsLib
{
void RadiusEdgeRatioMetric::calculateQuality()
{
    std::vector<MeshLib::Element*> const& elements(_mesh.getElements());
    std::size_t const nElements(_mesh.getNumberOfElements());
    for (std::size_t k(0); k < nElements; k++)
    {
        MeshLib::Element const& elem(*elements[k]);
        std::size_t const n_nodes(elem.getNumberOfBaseNodes());
        std::vector<MathLib::Point3d*> pnts(n_nodes);
        std::copy_n(elem.getNodes(), n_nodes, pnts.begin());
        GeoLib::MinimalBoundingSphere const s(pnts);
        auto const& [min, max] = computeSqrEdgeLengthRange(elem);
        _element_quality_metric[k] = std::sqrt(min) / (2 * s.getRadius());
    }
}
 
}  // namespace MeshToolsLib