MeshInformation.cpp

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#include "MeshInformation.h"
 
#include <range/v3/action/sort.hpp>
#include <range/v3/action/unique.hpp>
#include <range/v3/range/conversion.hpp>
#include <range/v3/view/unique.hpp>
 
#include "BaseLib/DemangleTypeInfo.h"
#include "MeshLib/Elements/Element.h"
#include "MeshQuality/MeshValidation.h"
 
namespace
{
template <typename T>
void printBounds(MeshLib::PropertyVector<T> const& property)
{
    auto const bounds = MeshToolsLib::MeshInformation::getValueBounds(property);
    if (!bounds.has_value())
    {
        INFO("\t{:s}: Could not get value bounds for property vector.",
             property.getPropertyName());
        return;
    }
    INFO("\t{:s}: ({:d} values) [{}, {}]", property.getPropertyName(),
         property.size(), bounds->first, bounds->second);
}
}  // namespace
 
namespace MeshToolsLib
{
GeoLib::AABB MeshInformation::getBoundingBox(const MeshLib::Mesh& mesh)
{
    const std::vector<MeshLib::Node*>& nodes(mesh.getNodes());
    return GeoLib::AABB(nodes.begin(), nodes.end());
}
 
std::map<MeshLib::MeshElemType, unsigned>
MeshInformation::getNumberOfElementTypes(const MeshLib::Mesh& mesh)
{
    std::map<MeshLib::MeshElemType, unsigned> n_element_types;
    const std::vector<MeshLib::Element*>& elements(mesh.getElements());
    for (auto element : elements)
    {
        MeshLib::MeshElemType t = element->getGeomType();
        n_element_types[t]++;
    }
    return n_element_types;
}
 
void MeshInformation::writeAllNumbersOfElementTypes(const MeshLib::Mesh& mesh)
{
    auto const& nr_ele_types =
        MeshToolsLib::MeshInformation::getNumberOfElementTypes(mesh);
 
    INFO("Number of elements in the mesh:");
    for (auto entry : nr_ele_types)
    {
        INFO("\t{:s}s: {:d}",
             MeshLib::MeshElemType2String(
                 static_cast<MeshLib::MeshElemType>(entry.first)),
             entry.second);
    }
}
 
void MeshInformation::writePropertyVectorInformation(const MeshLib::Mesh& mesh)
{
    auto const& properties = mesh.getProperties();
    INFO("There are {:d} properties in the mesh:", properties.size());
 
    for (auto [name, property] : properties)
    {
        if (auto const* p =
                dynamic_cast<MeshLib::PropertyVector<double>*>(property))
        {
            printBounds(*p);
        }
        else if (auto const* p =
                     dynamic_cast<MeshLib::PropertyVector<float>*>(property))
        {
            printBounds(*p);
        }
        else if (auto const* p =
                     dynamic_cast<MeshLib::PropertyVector<int>*>(property))
        {
            printBounds(*p);
        }
        else if (auto const* p =
                     dynamic_cast<MeshLib::PropertyVector<unsigned>*>(property))
        {
            printBounds(*p);
        }
        else if (auto const* p =
                     dynamic_cast<MeshLib::PropertyVector<long>*>(property))
        {
            printBounds(*p);
        }
        else if (auto const* p =
                     dynamic_cast<MeshLib::PropertyVector<long long>*>(
                         property))
        {
            printBounds(*p);
        }
        else if (auto const* p =
                     dynamic_cast<MeshLib::PropertyVector<unsigned long>*>(
                         property))
        {
            printBounds(*p);
        }
        else if (auto const* p =
                     dynamic_cast<MeshLib::PropertyVector<unsigned long long>*>(
                         property))
        {
            printBounds(*p);
        }
        else if (auto const* p =
                     dynamic_cast<MeshLib::PropertyVector<std::size_t>*>(
                         property))
        {
            printBounds(*p);
        }
        else if (auto const* p =
                     dynamic_cast<MeshLib::PropertyVector<char>*>(property))
        {
            printBounds(*p);
        }
        else
        {
            INFO(
                "\t{:s}: Could not get value bounds for property vector of "
                "type '{:s}'.",
                name, BaseLib::typeToString(*property));
        }
    }
}
 
void MeshInformation::writeMeshValidationResults(MeshLib::Mesh& mesh)
{
    INFO("Mesh Quality Control:");
    MeshToolsLib::MeshValidation::allNodesUsed(mesh);
    MeshToolsLib::MeshValidation::existCollapsibleNodes(mesh);
    MeshToolsLib::MeshValidation::evaluateElementGeometry(mesh);
 
    unsigned const n_holes(MeshToolsLib::MeshValidation::detectHoles(mesh));
    if (n_holes > 0)
    {
        INFO("{:d} hole(s) detected within the mesh", n_holes);
    }
    else
    {
        INFO("No holes found within the mesh.");
    }
}
 
std::vector<int> MeshInformation::getMaterialIDs(const MeshLib::Mesh& mesh)
{
    auto const* matids = MeshLib::materialIDs(mesh);
    if (!matids)
    {
        INFO("No MaterialIDs were found.");
        return {};
    }
 
    return *matids | ranges::views::unique | ranges::to<std::vector> |
           ranges::actions::sort | ranges::actions::unique |
           ranges::to<std::vector>;
}
 
}  // namespace MeshToolsLib