curvature.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
{
    namespace internal
    {
        template <class DataPoint, class _NFilter, typename T>
        class CurvatureEstimatorBase : public T
        {
            PONCA_FITTING_DECLARE_DEFAULT_TYPES
            PONCA_FITTING_DECLARE_MATRIX_TYPE
 
        protected:
            enum
            {
                PROVIDES_PRINCIPAL_CURVATURES
            };
 
        private:
            Scalar m_kmin{0},
                m_kmax{0};
            VectorType m_vmin{VectorType::Zero()},
                m_vmax{VectorType::Zero()};
 
            bool m_isValid{false};
 
            static_assert(DataPoint::Dim == 3, "CurvatureEstimatorBase is only valid in 3D");
 
        public:
            PONCA_FITTING_DECLARE_INIT
 
            PONCA_MULTIARCH [[nodiscard]] inline bool isValid() const { return m_isValid; }
 
            /**************************************************************************/
            /* Use results                                                            */
            /**************************************************************************/
            PONCA_MULTIARCH [[nodiscard]] inline Scalar kmin() const { return m_kmin; }
 
            PONCA_MULTIARCH [[nodiscard]] inline Scalar kmax() const { return m_kmax; }
 
            PONCA_MULTIARCH [[nodiscard]] inline VectorType kminDirection() const { return m_vmin; }
 
            PONCA_MULTIARCH [[nodiscard]] inline VectorType kmaxDirection() const { return m_vmax; }
 
            PONCA_MULTIARCH [[nodiscard]] inline Scalar kMean() const { return (m_kmin + m_kmax) / Scalar(2); }
 
            PONCA_MULTIARCH [[nodiscard]] inline Scalar GaussianCurvature() const { return m_kmin * m_kmax; }
 
        protected:
            PONCA_MULTIARCH inline void setCurvatureValues(Scalar kmin, Scalar kmax, const VectorType& vmin,
                                                           const VectorType& vmax);
        };
    } // namespace internal
 
    template <class DataPoint, class _NFilter, typename T>
    class CurvatureEstimator : public internal::CurvatureEstimatorBase<DataPoint, _NFilter, T>
    {
    public:
        PONCA_EXPLICIT_CAST_OPERATORS(CurvatureEstimator, curvatureEstimator)
    };
 
    template <class DataPoint, class _NFilter, int DiffType, typename T>
    class CurvatureEstimatorDer : public internal::CurvatureEstimatorBase<DataPoint, _NFilter, T>
    {
    public:
        PONCA_EXPLICIT_CAST_OPERATORS_DER(CurvatureEstimatorDer, curvatureEstimator)
    };
 
} // namespace Ponca
 
#include "curvature.hpp"