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
    };
}
 
    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"