sphereFit.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"
 
#include "../primitive/algebraicSphere.h" // used to define SphereFit
 
#define SPHERE_FIT_REQUIREMENTS ProvidesAlgebraicSphere<T>
 
namespace Ponca
{
 
    template <class DataPoint, class _NFilter, typename T>
        requires SPHERE_FIT_REQUIREMENTS
    class SphereFitImpl : public T
    {
        PONCA_FITTING_DECLARE_DEFAULT_TYPES
    protected:
        using VectorA = Eigen::Matrix<Scalar, DataPoint::Dim + 2, 1>;
        using MatrixA = Eigen::Matrix<Scalar, DataPoint::Dim + 2, DataPoint::Dim + 2>;
 
    public:
        using Solver = Eigen::EigenSolver<MatrixA>;
 
    protected:
        // computation data
        MatrixA m_matA{MatrixA::Zero()}; 
        Solver m_solver;
 
    public:
        PONCA_EXPLICIT_CAST_OPERATORS(SphereFitImpl, sphereFit)
        PONCA_FITTING_DECLARE_INIT_ADD_FINALIZE
        PONCA_FITTING_IS_SIGNED(false)
 
        PONCA_MULTIARCH inline const Solver& solver() const { return m_solver; }
 
    }; // class SphereFit
 
    template <class DataPoint, class _NFilter, typename T>
    using SphereFit = SphereFitImpl<DataPoint, _NFilter, AlgebraicSphere<DataPoint, _NFilter, T>>;
 
#include "sphereFit.hpp"
 
} // namespace Ponca