quality.h

/*****************************************************************************
 * VCLib                                                                     *
 * Visual Computing Library                                                  *
 *                                                                           *
 * Copyright(C) 2021-2025                                                    *
 * Visual Computing Lab                                                      *
 * ISTI - Italian National Research Council                                  *
 *                                                                           *
 * All rights reserved.                                                      *
 *                                                                           *
 * This program is free software; you can redistribute it and/or modify      *
 * it under the terms of the Mozilla Public License Version 2.0 as published *
 * by the Mozilla Foundation; either version 2 of the License, or            *
 * (at your option) any later version.                                       *
 *                                                                           *
 * This program is distributed in the hope that it will be useful,           *
 * but WITHOUT ANY WARRANTY; without even the implied warranty of            *
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the              *
 * Mozilla Public License Version 2.0                                        *
 * (https://www.mozilla.org/en-US/MPL/2.0/) for more details.                *
 ****************************************************************************/
 
#ifndef VCL_ALGORITHMS_MESH_UPDATE_QUALITY_H
#define VCL_ALGORITHMS_MESH_UPDATE_QUALITY_H
 
#include <vclib/algorithms/core/polygon.h>
#include <vclib/algorithms/mesh/stat.h>
#include <vclib/mesh/requirements.h>
 
namespace vcl {
 
template<MeshConcept MeshType>
void setPerVertexQuality(
    MeshType&                                  m,
    typename MeshType::VertexType::QualityType s)
{
    requirePerVertexQuality(m);
 
    using VertexType = MeshType::VertexType;
 
    for (VertexType& v : m.vertices()) {
        v.quality() = s;
    }
}
 
template<FaceMeshConcept MeshType>
void setPerFaceQuality(MeshType& m, typename MeshType::FaceType::QualityType s)
{
    requirePerFaceQuality(m);
 
    using FaceType = MeshType::FaceType;
 
    for (FaceType& f : m.faces()) {
        f.quality() = s;
    }
}
 
template<MeshConcept MeshType>
void clampPerVertexQuality(
    MeshType&                                  m,
    typename MeshType::VertexType::QualityType minS,
    typename MeshType::VertexType::QualityType maxS)
{
    requirePerVertexQuality(m);
 
    using VertexType = MeshType::VertexType;
 
    for (VertexType& v : m.vertices()) {
        v.quality() = std::min(maxS, std::max(minS, v.quality()));
    }
}
 
template<FaceMeshConcept MeshType>
void clampPerFaceQuality(
    MeshType&                                m,
    typename MeshType::FaceType::QualityType minS,
    typename MeshType::FaceType::QualityType maxS)
{
    requirePerFaceQuality(m);
 
    using FaceType = MeshType::FaceType;
 
    for (FaceType& f : m.faces()) {
        f.quality() = std::min(maxS, std::max(minS, f.quality()));
    }
}
 
template<MeshConcept MeshType>
void normalizePerVertexQuality(
    MeshType&                                  m,
    typename MeshType::VertexType::QualityType minS = 0,
    typename MeshType::VertexType::QualityType maxS = 1)
{
    requirePerVertexQuality(m);
 
    using VertexType  = MeshType::VertexType;
    using QualityType = VertexType::QualityType;
 
    QualityType                         range = maxS - minS;
    std::pair<QualityType, QualityType> p     = vertexQualityMinMax(m);
 
    for (VertexType& v : m.vertices()) {
        v.quality() =
            minS + range * ((v.quality() - p.first) / (p.second - p.first));
    }
}
 
template<FaceMeshConcept MeshType>
void normalizePerFaceQuality(
    MeshType&                                m,
    typename MeshType::FaceType::QualityType minS = 0,
    typename MeshType::FaceType::QualityType maxS = 1)
{
    requirePerFaceQuality(m);
 
    using FaceType    = MeshType::FaceType;
    using QualityType = FaceType::QualityType;
 
    QualityType                         range = maxS - minS;
    std::pair<QualityType, QualityType> p     = faceQualityMinMax(m);
 
    for (FaceType& f : m.faces()) {
        f.quality() =
            minS + range * ((f.quality() - p.first) / (p.second - p.first));
    }
}
 
template<FaceMeshConcept MeshType>
void setPerVertexQualityFromVertexValence(MeshType& m)
{
    requirePerVertexQuality(m);
 
    using VertexType = MeshType::VertexType;
    using FaceType   = MeshType::FaceType;
 
    setPerVertexQuality(m, 0);
 
    for (FaceType& f : m.faces()) {
        for (VertexType* v : f.vertices()) {
            v->quality() += 1;
        }
    }
}
 
template<FaceMeshConcept MeshType>
void setPerFaceQualityFromFaceArea(MeshType& m)
{
    requirePerFaceQuality(m);
 
    using FaceType = MeshType::FaceType;
 
    for (FaceType& f : m.faces()) {
        f.quality() = faceArea(f);
    }
}
 
template<MeshConcept MeshType>
void setPerVertexQualityFromPrincipalCurvatureGaussian(MeshType& m)
{
    requirePerVertexQuality(m);
    requirePerVertexPrincipalCurvature(m);
 
    using VertexType = MeshType::VertexType;
 
    for (VertexType& v : m.vertices()) {
        v.quality() = v.principalCurvature().maxValue() *
                      v.principalCurvature().minValue();
    }
}
 
template<MeshConcept MeshType>
void setPerVertexQualityFromPrincipalCurvatureMean(MeshType& m)
{
    requirePerVertexQuality(m);
    requirePerVertexPrincipalCurvature(m);
 
    using VertexType = MeshType::VertexType;
 
    for (VertexType& v : m.vertices()) {
        v.quality() = (v.principalCurvature().maxValue() +
                       v.principalCurvature().minValue()) /
                      2;
    }
}
 
template<MeshConcept MeshType>
void setPerVertexQualityFromPrincipalCurvatureMinValue(MeshType& m)
{
    requirePerVertexQuality(m);
    requirePerVertexPrincipalCurvature(m);
 
    using VertexType = MeshType::VertexType;
 
    for (VertexType& v : m.vertices()) {
        v.quality() = v.principalCurvature().minValue();
    }
}
 
template<MeshConcept MeshType>
void setPerVertexQualityFromPrincipalCurvatureMaxValue(MeshType& m)
{
    requirePerVertexQuality(m);
    requirePerVertexPrincipalCurvature(m);
 
    using VertexType = MeshType::VertexType;
 
    for (VertexType& v : m.vertices()) {
        v.quality() = v.principalCurvature().maxValue();
    }
}
 
template<MeshConcept MeshType>
void setPerVertexQualityFromPrincipalCurvatureShapeIndex(MeshType& m)
{
    requirePerVertexQuality(m);
    requirePerVertexPrincipalCurvature(m);
 
    using VertexType = MeshType::VertexType;
    using ScalarType = VertexType::PrincipalCurvatureType::ScalarType;
 
    for (VertexType& v : m.vertices()) {
        ScalarType k1 = v.principalCurvature().maxValue();
        ScalarType k2 = v.principalCurvature().minValue();
        if (k1 < k2)
            std::swap(k1, k2);
        v.quality() = (2.0 / M_PI) * std::atan2(k1 + k2, k1 - k2);
    }
}
 
template<MeshConcept MeshType>
void setPerVertexQualityFromPrincipalCurvatureCurvedness(MeshType& m)
{
    requirePerVertexQuality(m);
    requirePerVertexPrincipalCurvature(m);
 
    using VertexType = MeshType::VertexType;
    using ScalarType = VertexType::PrincipalCurvatureType::ScalarType;
 
    for (VertexType& v : m.vertices()) {
        ScalarType k1 = v.principalCurvature().maxValue();
        ScalarType k2 = v.principalCurvature().minValue();
 
        v.quality() = std::sqrt((k1 * k1 + k2 * k2) / 2.0);
    }
}
 
} // namespace vcl
 
#endif // VCL_ALGORITHMS_MESH_UPDATE_QUALITY_H