STK_BSplineCoefficients.h

Go to the documentation of this file.

/*--------------------------------------------------------------------*/
/*     Copyright (C) 2004-2010  Serge Iovleff

 This program is free software; you can redistribute it and/or modify
 it under the terms of the GNU Lesser General Public License as
 published by the Free Software Foundation; either version 2 of the
 License, or (at your option) any later version.

 This program is distributed in the hope that it will be useful,
 but WITHOUT ANY WARRANTY; without even the implied warranty of
 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 GNU Lesser General Public License for more details.

 You should have received a copy of the GNU Lesser General Public
 License along with this program; if not, write to the
 Free Software Foundation, Inc.,
 59 Temple Place,
 Suite 330,
 Boston, MA 02111-1307
 USA

 Contact : Serge.Iovleff@stkpp.org
 */

/*
 * Project:  stkpp::Regress
 * created on: 25 juin 2010
 * Purpose:  Compute the coefficient of a B-Spline curves.
 * Author:   iovleff, serge.iovleff@stkpp.org
 **/

#ifndef STK_BSPLINECOEFFICIENTS_H
#define STK_BSPLINECOEFFICIENTS_H

#include "../../Arrays/include/STK_Matrix.h"

namespace STK
{


class BSplineCoefficients
{
  public:
    enum KnotsPosition
    {
      uniform_  
    , periodic_ 
    , density_  
    , unknown_  
    };

    static KnotsPosition StringToKnotsPosition( String const& type);

    static String KnotsPositionToString( KnotsPosition const& type);

    BSplineCoefficients( Vector const* p_data
                       , Integer const& nbControlPoints
                       , Integer const& degree = 3
                       , KnotsPosition const& position = uniform_
                       );

    BSplineCoefficients( Vector const& data
                       , Integer const& nbControlPoints
                       , Integer const& degree = 3
                       , KnotsPosition const& position = uniform_
                       );

    virtual ~BSplineCoefficients();

    void run();

    void setData( Vector const* p_data
                , Integer const& nbControlPoints
                , Integer const& degree = 3
                , KnotsPosition const& position = uniform_
                );

    inline Integer const& degree() const { return degree_;}
    inline Integer const& nbKnots() const { return nbKnots_;}
    inline Integer const& nbControlPoints() const { return nbControlPoints_;}
    inline Vector const& knots() const { return knots_;}
    inline Matrix const& coefficients() const { return coefficients_;}

  protected:
    Vector const* p_data_;
    Integer nbKnots_;
    Integer lastKnot_;
    Integer nbControlPoints_;
    Integer lastControlPoint_;
    Integer degree_;
    KnotsPosition position_;
    Vector knots_;
    Matrix coefficients_;

  private:
    Real minValue_;
    Real maxValue_;
    void computeKnots();
    void computeUniformKnots();
    void computePeriodicKnots();
    void computeDensityKnots();
    void computeCoefficients();
    void computeCoefficientsRow(Integer const& irow);
};

}

#endif /* STK_BSPLINECOEFFICIENTS_H */