FuncGrid_Bspline.h

Go to the documentation of this file.

/*******************************************************************
 *  Copyright (C) 2010 by Federico Marulli and Alfonso Veropalumbo *
 *  federico.marulli3@unibo.it                                     *
 *                                                                 *
 *  This program is free software; you can redistribute it and/or  *
 *  modify it under the terms of the GNU 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 General Public License for more details.                   *
 *                                                                 *
 *  You should have received a copy of the GNU 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.      *
 *******************************************************************/
 
#ifndef __FUNCGRIDBSPL__
#define __FUNCGRIDBSPL__ 
 
#include "GSLwrapper.h"
 
// =====================================================================================
 
 
namespace cbl {
 
  namespace glob {
    
    class FuncGrid_Bspline
    {
      private:
 
    std::vector<double> m_x;
 
    std::vector<double> m_fx;
 
    int m_nbreakpoints;
 
    int m_ncoefficients;
 
    int m_order;
 
    std::shared_ptr<gsl_bspline_workspace> m_bspline;
 
    std::shared_ptr<gsl_vector> m_Bcoeff; 
    
    std::shared_ptr<gsl_vector> m_weights; 
    
    std::shared_ptr<gsl_matrix> m_covariance;
 
    double m_integral;
 
 
    void m_set_bspline(const std::vector<double> x, const std::vector<double> fx, const int nbreakpoints, const int order);
 
    void m_set_knots (const double xmin=cbl::par::defaultDouble, const double xmax=cbl::par::defaultDouble);
 
    void m_set_knots (const std::vector<double> breakpoints);
 
    void m_linear_fit (const double frac=0.1);
 
    void m_compute_func_integral ();
    
 
      public:
 
    FuncGrid_Bspline () {}
 
    FuncGrid_Bspline (const std::vector<double> x, const std::vector<double> fx, const int nbreakpoints, const int order=4, const double frac=0.1, const double xmin=cbl::par::defaultDouble, const double xmax=cbl::par::defaultDouble);
 
    FuncGrid_Bspline (const std::vector<double> x, const std::vector<double> fx, const std::vector<double> breakpoints, const int order=4, const double frac=0.1);
 
    ~FuncGrid_Bspline () = default;
 
 
 
    void set (const std::vector<double> x, const std::vector<double> fx, const int nbreakpoints, const int order=4, const double frac=0.1, const double xmin=cbl::par::defaultDouble, const double xmax=cbl::par::defaultDouble);
 
    void set (const std::vector<double> x, const std::vector<double> fx, const std::vector<double> breakpoints, const int order=4, const double frac=0.1);
        
 
 
        double operator () (const double xx, const double integral=cbl::par::defaultDouble) const;
 
        std::vector<double> eval_func (const std::vector<double> xx, const double integral=cbl::par::defaultDouble) const;
        
    };
 
  }
}
#endif