mlsSphereFitDer.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"
#include "../concepts.h"
 
#define MLS_SPHERE_FIT_DER_REQUIREMENTS \
    ProvidesAlgebraicSphere<T>&& ProvidesImplicitPrimitiveDerivative<T>&& ProvidesNormal<DataPoint>
 
namespace Ponca
{


    template <class DataPoint, class _NFilter, int DiffType, typename T>
        requires MLS_SPHERE_FIT_DER_REQUIREMENTS
    class MlsSphereFitDer : public T
    {
        PONCA_FITTING_DECLARE_DEFAULT_TYPES
        PONCA_FITTING_DECLARE_DEFAULT_DER_TYPES
    protected:
        enum
        {
            Dim    = DataPoint::Dim,     
            DerDim = Base::NbDerivatives 
        };
 
    public:
        using Matrix = Eigen::Matrix<Scalar, DerDim, DerDim>;


        using MatrixArray = Eigen::Matrix<Scalar, DerDim, Dim * DerDim>;
 
    protected:
        // computation data
        Matrix m_d2SumDotPN, 
            m_d2SumDotPP,    
            m_d2SumW;        
 
        MatrixArray m_d2SumP, 
            m_d2SumN;         
 
    public:
        // results
        Matrix m_d2Uc,      
            m_d2Uq;         
        MatrixArray m_d2Ul; 
 
    public:
        PONCA_EXPLICIT_BASKET_CAST_OPERATORS_DER(MlsSphereFitDer, mlsSphereFitDer)
        PONCA_EXPLICIT_BASKET_CAST_OPERATORS_DER(MlsSphereFitDer, implicitPrimitiveDer)
 
        PONCA_FITTING_DECLARE_INIT_ADDDER_FINALIZE

        PONCA_MULTIARCH [[nodiscard]] inline ScalarArray dPotential() const;

        PONCA_MULTIARCH [[nodiscard]] inline VectorType primitiveGradient() const;

        PONCA_MULTIARCH [[nodiscard]] inline VectorArray dNormal() const;
 
    }; // class MlsSphereFitDer
 
#include "mlsSphereFitDer.hpp"
 
} // namespace Ponca