weingarten.h

/*
This Source Code Form is subject to the terms of the Mozilla Public
 License, v. 2.0. If a copy of the MPL was not distributed with this
 file, You can obtain one at http://mozilla.org/MPL/2.0/.
*/
 
#pragma once
 
#include "./defines.h"
 
namespace Ponca
{
    template <class DataPoint, class _NFilter, typename T>
    class FundamentalFormWeingartenEstimator : public T
    {
        PONCA_FITTING_DECLARE_DEFAULT_TYPES
        using Matrix2 = Eigen::Matrix<Scalar, 2, 2>;
        static_assert(DataPoint::Dim == 3, "FundamentalFormWeingartenEstimator is only valid in 3D");
 
    protected:
        enum
        {
            Check =
                Base::PROVIDES_FIRST_FUNDAMENTAL_FORM_COMPONENTS && Base::PROVIDES_SECOND_FUNDAMENTAL_FORM_COMPONENTS,
            PROVIDES_WEINGARTEN_MAP
        };
 
    public:
        PONCA_EXPLICIT_CAST_OPERATORS(FundamentalFormWeingartenEstimator, fundamentalFormWeingartenEstimator)
 
        
        PONCA_MULTIARCH [[nodiscard]] inline Matrix2 firstFundamentalForm() const;
 
        template <typename Matrix2Derived>
        PONCA_MULTIARCH inline void firstFundamentalForm(Matrix2Derived& first) const;
 
        PONCA_MULTIARCH [[nodiscard]] inline Matrix2 secondFundamentalForm() const;
 
        template <typename Matrix2Derived>
        PONCA_MULTIARCH inline void secondFundamentalForm(Matrix2Derived& second) const;
 
        PONCA_MULTIARCH [[nodiscard]] inline Matrix2 weingartenMap() const;
 
        template <typename Matrix2Derived>
        PONCA_MULTIARCH inline void weingartenMap(Matrix2Derived& w) const;
 
        PONCA_MULTIARCH inline Scalar kMean() const;
 
        PONCA_MULTIARCH inline Scalar GaussianCurvature() const;
    };
 
    template <class DataPoint, class _NFilter, int DiffType, typename T>
    class NormalDerivativeWeingartenEstimator : public T
    {
        PONCA_FITTING_DECLARE_DEFAULT_TYPES
        PONCA_FITTING_DECLARE_MATRIX_TYPE
        PONCA_FITTING_DECLARE_DEFAULT_DER_TYPES
        using Matrix2 = Eigen::Matrix<Scalar, 2, 2>;
        static_assert(DataPoint::Dim == 3, "NormalDerivativeWeingartenEstimator is only valid in 3D");
        static_assert(Base::isSpaceDer(), "NormalDerivativeWeingartenEstimator requires spatial derivation");
 
    protected:
        enum
        {
            Check = Base::PROVIDES_NORMAL_DERIVATIVE,
            PROVIDES_WEINGARTEN_MAP,
            PROVIDES_TANGENT_PLANE_BASIS
        };
 
    private:
        MatrixType m_tangentBasis{MatrixType::Zero()};
 
    public:
        PONCA_EXPLICIT_CAST_OPERATORS_DER(NormalDerivativeWeingartenEstimator, normalDerivativeWeingartenEstimator)
        PONCA_FITTING_DECLARE_FINALIZE
 
        PONCA_MULTIARCH [[nodiscard]] inline Matrix2 weingartenMap() const;
 
        template <typename Matrix2Derived>
        PONCA_MULTIARCH inline void weingartenMap(Matrix2Derived& w) const;
 
        PONCA_MULTIARCH inline VectorType worldToTangentPlane(const VectorType& _q,
                                                              bool _isPositionVector = true) const;
 
        PONCA_MULTIARCH inline VectorType tangentPlaneToWorld(const VectorType& _q,
                                                              bool _isPositionVector = true) const;
    };
 
    namespace internal
    {
        template <class DataPoint, class _NFilter, typename T>
        class WeingartenCurvatureEstimatorBase : public T
        {
            PONCA_FITTING_DECLARE_DEFAULT_TYPES
            using Matrix2 = Eigen::Matrix<Scalar, 2, 2>;
            static_assert(DataPoint::Dim == 3, "WeingartenCurvatureEstimator is only valid in 3D");
 
        protected:
            enum
            {
                Check = Base::PROVIDES_TANGENT_PLANE_BASIS &&  // required for tangentPlaneToWorld
                        Base::PROVIDES_PRINCIPAL_CURVATURES && // required curvature storage
                        Base::PROVIDES_WEINGARTEN_MAP
            };
 
        public:
            PONCA_FITTING_DECLARE_FINALIZE
        };
    } // namespace internal
 
    template <class DataPoint, class _NFilter, typename T>
    struct WeingartenCurvatureEstimator : public internal::WeingartenCurvatureEstimatorBase<DataPoint, _NFilter, T>
    {
        PONCA_FITTING_DECLARE_DEFAULT_TYPES
        PONCA_EXPLICIT_CAST_OPERATORS(WeingartenCurvatureEstimator, weingartenCurvatureEstimator)
    };
 
    template <class DataPoint, class _NFilter, int DiffType, typename T>
    struct WeingartenCurvatureEstimatorDer : public internal::WeingartenCurvatureEstimatorBase<DataPoint, _NFilter, T>
    {
        PONCA_FITTING_DECLARE_DEFAULT_TYPES
        PONCA_FITTING_DECLARE_DEFAULT_DER_TYPES
        PONCA_EXPLICIT_CAST_OPERATORS_DER(WeingartenCurvatureEstimatorDer, weingartenCurvatureEstimator)
    };
 
} // namespace Ponca
 
#include "weingarten.hpp"