Spline8interp.h

#pragma once

/*interpoliert Punktmenge mit 8 Basisfunktionen folgender Lagrangebasis:
Zu n+1 Knoten n  intervallweise 
Splines mit Funktionswerten [0 0 ...1 ...0] (sog. Lagrange-Basis). Periodizitätwird gewaehrleistet, indem eine ungerade Knotenzahl verwendet wird und der
Basisspline fuer den mittleren knoten (bei Null-RB an die erste Abletung, siehe
octave help spline, "complete spline") nur verschoben wird.
*/
class Spline8interp :
    public H1interp
{
protected:
    static double baseCoeffs[8*4*8];
    //sonst kann die abgeleitete Klasse nicht auf die eigenen Koeffs zugreifen:
    double* baseCoeffsPtr;
    double evalPP    (double x, const  double  coeffs[4] ) const;
    double findInterval( unsigned int *i, double x ) const;
    void    evalMonoms( double baseVals[4], const double x ) const ;
    void    evalMonomsD( double baseVals[4], const double x ) const;

public:
    double evalSpline    (double x, const double  coeffs[8*4] ) const;
        virtual double  operator()(double * coeffs, const cv::Mat pts, const cv::Mat& firstDerivs) { return -1; }
    int  operator()( double& residual,
            double&     maxDist,
            cv::Mat&     approxCoords,
            size_t&      sysSize,
            const std::vector<cv::Mat> arg,
            const std::vector<cv::Mat> val,
            const std::vector<cv::Mat> argD,
            const std::vector<cv::Mat> valD,
            const std::vector<bool>& combineFlag);
    /*Dieser Konstruktor kannn von abgeleiteten Klassen verwendet werden, die hoehere Zahl
    von splines verwenden. Genaugenommen  ist das natuerlich keine Ist-ein-Beziehung, 
    man muesste eine gemeinsame Basis machen. Andererseits ist das nicht schlimm.*/
    Spline8interp(const  int nrBasisFcts_In = 8);
    ~Spline8interp(void);
    private:
/*    const double* diffx_i;
    double p(double x,int i);
    double p_(double x,int i);
    double q(double x,int i);
    double q_(double x,int i);
*/
};


wird gewaehrleistet, indem eine ungerade Knotenzahl verwendet wird und der
Basisspline fuer den mittleren knoten (bei Null-RB an die erste Abletung, siehe
octave help spline, "complete spline") nur verschoben wird.
*/
class Spline8interp :
    public H1interp
{
protected:
    static double baseCoeffs[8*4*8];
    //sonst kann die abgeleitete Klasse nicht auf die eigenen Koeffs zugreifen:
    double* baseCoeffsPtr;
    double evalPP    (double x, const  double  coeffs[4] ) const;
    double findInterval( unsigned int *i, double x ) const;
    void    evalMonoms( double baseVals[4], const double x ) const ;
    void    evalMonomsD( double baseVals[4], const double x ) const;

public:
    double evalSpline    (double x, const double  coeffs[8*4] ) const;
        virtual double  operator()(double * coeffs, const cv::Mat pts, const cv::Mat& firstDerivs) { return -1; }
    int  operator()( double& residual,
            double&     maxDist,
            cv::Mat&     approxCoords,
            size_t&      sysSize,
            const std::vector<cv::Mat> arg,
            const std::vector<cv::Mat> val,
            const std::vector<cv::Mat> argD,
            const std::vector<cv::Mat> valD,
            const std::vector<bool>& combineFlag);
    /*Dieser Konstruktor kannn von abgeleiteten Klassen verwendet werden, die hoehere Zahl
    von splines verwenden. Genaugenommen  ist das natuerlich keine Ist-ein-Beziehung, 
    man muesste eine gemeinsame Basis machen. Andererseits ist das nicht schlimm.*/
    Spline8interp(const  int nrBasisFcts_In = 8);
    ~Spline8interp(void);
    private:
/*    const double* diffx_i;
    double p(double x,int i);
    double p_(double x,int i);
    double q(double x,int i);
    double q_(double x,int i);
*/
};