gls.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 "./curvature.h"
 
namespace Ponca
{
 
template < class DataPoint, class _WFunctor, typename T>
class GLSParam : public T
{
    PONCA_FITTING_DECLARE_DEFAULT_TYPES
 
protected:
    enum
    {
        Check = Base::PROVIDES_ALGEBRAIC_SPHERE,
        PROVIDES_GLS_PARAMETRIZATION
    };
 
protected:
    Scalar m_fitness {0};   
public:
    PONCA_EXPLICIT_CAST_OPERATORS(GLSParam,glsParam)
    PONCA_FITTING_DECLARE_FINALIZE
 
    /**************************************************************************/
    /* Use results                                                            */
    /**************************************************************************/
    PONCA_MULTIARCH inline Scalar tau() const
    {
        return Base::isNormalized() ? Base::m_uc : Base::m_uc / Base::prattNorm();
    }
 
    PONCA_MULTIARCH inline VectorType eta() const { return Base::primitiveGradient(); }
 
    PONCA_MULTIARCH inline Scalar kappa() const
    {
        return Scalar(2.) * (Base::isNormalized() ? Base::m_uq : Base::m_uq / Base::prattNorm());
    }
 
    PONCA_MULTIARCH inline Scalar tau_normalized() const { return tau() / Base::m_w.evalScale(); }
 
    PONCA_MULTIARCH inline VectorType eta_normalized() const { return eta(); }
 
    PONCA_MULTIARCH inline Scalar kappa_normalized() const { return kappa() * Base::m_w.evalScale(); }
 
    PONCA_MULTIARCH inline Scalar fitness() const { return m_fitness; }
 
    PONCA_MULTIARCH inline Scalar compareTo (const GLSParam<DataPoint, _WFunctor, T>& _other,
                                        bool _useFitness = true) const
    {
        Scalar nTau     = this->tau_normalized()   - _other.tau_normalized();
        Scalar nKappa   = this->kappa_normalized() - _other.kappa_normalized();
        Scalar nFitness = _useFitness ? this->fitness() - _other.fitness() : Scalar(0.);
 
        return nTau * nTau + nKappa * nKappa + nFitness * nFitness;
    }
 
}; //class GLSParam
 
 
template < class DataPoint, class _WFunctor, int DiffType, typename T>
class GLSDer : public T
{
PONCA_FITTING_DECLARE_DEFAULT_TYPES
PONCA_FITTING_DECLARE_DEFAULT_DER_TYPES
 
protected:
    enum
    {
        Check = Base::PROVIDES_GLS_PARAMETRIZATION &
                Base::PROVIDES_PRIMITIVE_DERIVATIVE &
                Base::PROVIDES_ALGEBRAIC_SPHERE_DERIVATIVE,
        PROVIDES_GLS_DERIVATIVE,
        PROVIDES_GLS_GEOM_VAR
    };
 
public:
    PONCA_EXPLICIT_CAST_OPERATORS_DER(GLSDer,glsDer)
 
    PONCA_MULTIARCH inline ScalarArray dtau()   const; 
    PONCA_MULTIARCH inline VectorArray deta()   const; 
    PONCA_MULTIARCH inline ScalarArray dkappa() const; 
    PONCA_MULTIARCH inline ScalarArray dtau_normalized()   const; 
    PONCA_MULTIARCH inline VectorArray deta_normalized()   const; 
    PONCA_MULTIARCH inline ScalarArray dkappa_normalized() const; 
    PONCA_MULTIARCH inline Scalar geomVar(Scalar wtau   = Scalar(1),
                                          Scalar weta   = Scalar(1),
                                          Scalar wkappa = Scalar(1)) const;
}; //class GLSScaleDer
 
 
#include "gls.hpp"
 
} //namespace Ponca