principal_curvature.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_MESH_COMPONENTS_PRINCIPAL_CURVATURE_H
#define VCL_MESH_COMPONENTS_PRINCIPAL_CURVATURE_H
 
#include "base/component.h"
#include "base/predicates.h"
 
#include <vclib/space/core.h>
 
namespace vcl::comp {
 
template<typename Scalar, typename ParentElemType = void, bool OPT = false>
class PrincipalCurvature :
        public Component<
            PrincipalCurvature<Scalar, ParentElemType, OPT>,
            CompId::PRINCIPAL_CURVATURE,
            vcl::PrincipalCurvature<Scalar>,
            ParentElemType,
            !std::is_same_v<ParentElemType, void>,
            OPT>
{
    using Base = Component<
        PrincipalCurvature<Scalar, ParentElemType, OPT>,
        CompId::PRINCIPAL_CURVATURE,
        vcl::PrincipalCurvature<Scalar>,
        ParentElemType,
        !std::is_same_v<ParentElemType, void>,
        OPT>;
 
public:
    using PrincipalCurvatureType = vcl::PrincipalCurvature<Scalar>;
 
    /* Constructors */
 
    PrincipalCurvature() = default;
 
    /* Member functions */
 
    const PrincipalCurvatureType& principalCurvature() const
    {
        return Base::data();
    }
 
    PrincipalCurvatureType& principalCurvature() { return Base::data(); }
 
protected:
    // Component interface functions
    template<typename Element>
    void importFrom(const Element& e, bool = true);
 
    void serialize(std::ostream& os) const
    {
        principalCurvature().serialize(os);
    }
 
    void deserialize(std::istream& is) { principalCurvature().deserialize(is); }
};
 
/* concepts */
 
template<typename T>
concept HasPrincipalCurvature =
    TTB::IsDerivedFromSpecializationOfV<T, PrincipalCurvature>;
 
template<typename T>
concept HasOptionalPrincipalCurvature =
    HasPrincipalCurvature<T> &&
    IsOptionalComponent<typename RemoveRef<T>::PrincipalCurvature>;
 
/* importFrom function */
 
template<typename Scalar, typename ParentElemType, bool OPT>
template<typename Element>
void PrincipalCurvature<Scalar, ParentElemType, OPT>::importFrom(
    const Element& e,
    bool)
{
    if constexpr (HasPrincipalCurvature<Element>) {
        if (isPrincipalCurvatureAvailableOn(e)) {
            principalCurvature() =
                e.principalCurvature().template cast<Scalar>();
        }
    }
}
 
/* Detector function to check if a class has PrincipalCurvature available */
 
bool isPrincipalCurvatureAvailableOn(const auto& element)
{
    return isComponentAvailableOn<CompId::PRINCIPAL_CURVATURE>(element);
}
 
/* Specialization Aliases */
 
template<typename ElementType = void, bool OPT = false>
using PrincipalCurvaturef = PrincipalCurvature<float, ElementType, OPT>;
 
template<typename ElementType = void, bool OPT = false>
using PrincipalCurvatured = PrincipalCurvature<double, ElementType, OPT>;
 
} // namespace vcl::comp
 
#endif // VCL_MESH_COMPONENTS_PRINCIPAL_CURVATURE_H