FuncGrid.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 __FUNCGRID__
#define __FUNCGRID__ 
 
#include "GSLwrapper.h"
#include "CUBAwrapper.h"
 
 
// =====================================================================================
 
 
namespace cbl {
 
  namespace glob {
    
    class FuncGrid
    {
      
    private:
 
      std::vector<double> m_x;
 
      std::vector<double> m_y;
      
      size_t m_size;
      
      std::string m_interpType;
 
      std::shared_ptr<gsl_spline> m_spline;
 
      const gsl_interp_type *m_type;
 
      std::shared_ptr<gsl_interp_accel> m_acc;
 
      double m_xmin;
 
      double m_xmax;
 
      BinType m_binType;
 
      
    public:
 
      
      FuncGrid () = default;
 
      FuncGrid (const std::vector<double> x, const std::vector<double> y, const std::string interpType, const BinType bin_type=BinType::_linear_);
 
      ~FuncGrid () = default;
 
 
      
 
      std::vector<double> x () const { return m_x; }
      
      double x (const int i) const { return m_x[i]; }
      
      std::vector<double> y () const { return m_y; }
      
      double y (const int i) const { return m_y[i]; }
 
      double size () const { return m_size; }
 
      double xmin () const { return m_xmin; }
 
      double xmax () const { return m_xmax; }
             
 
      
 
      void free ();
 
      double operator () (const double xx) const;
 
      std::vector<double> eval_func (const std::vector<double> xx) const;
 
      double D1v (const double xx) const;
      
      double D2v (const double xx) const;
   
      double integrate_qag (const double a, const double b, const double rel_err=1.e-2, const double abs_err=1.e-6, const int limit_size=1000, const int rule=6);
 
      double integrate_cquad (const double a, const double b, const double rel_err=1.e-2, const double abs_err=1.e-6, const int nevals=100);
 
      double integrate_qaws (const double a, const double b, const double alpha=0, const double beta=0, const int mu=0, const int nu=0, const double rel_err=1.e-2, const double abs_err=1.e-6, const int limit_size=1000);
      
      double root (const double x_low, const double x_up, const double fx0=0, const double rel_err=1.e-2, const double abs_err=1.e-6);
 
      double root_D1v (const double x_low, const double x_up, const double fx0=0, const double rel_err=1.e-2, const double abs_err=1.e-6); 
 
      double root_D2v (const double x_low, const double x_up, const double fx0=0, const double rel_err=1.e-2, const double abs_err=1.e-6);
    };
 
    class FuncGrid2D
    {
      
    private:
 
      std::vector<double> m_x;
      
      std::vector<double> m_y;
 
      std::shared_ptr<double> m_fxy;
      
      size_t m_size_x;
 
      size_t m_size_y; 
      
      std::string m_interpType;
 
      const gsl_interp2d_type *m_type;
 
      std::shared_ptr<gsl_interp_accel> m_acc_x;
 
      std::shared_ptr<gsl_interp_accel> m_acc_y;
 
      std::shared_ptr<gsl_spline2d> m_spline;
 
      double m_xmin;
 
      double m_xmax;
 
      double m_ymin;
 
      double m_ymax;
      
    public:
 
      
      FuncGrid2D () = default;
 
      FuncGrid2D (const std::vector<double> x, const std::vector<double> y, const std::vector<std::vector<double>> fxy, const std::string interpType);
 
      ~FuncGrid2D () = default;
 
 
      
 
      std::vector<double> x () const { return m_x; }
 
      double x (const int i) const { return m_x[i]; }
      
      std::vector<double> y () const { return m_y; }
      
      double y (const int i) const { return m_y[i]; }
 
      double size_x () const { return m_size_x; }
 
      double xmin () const { return m_xmin; }
 
      double xmax () const { return m_xmax; }
        
      double size_y () const { return m_size_y; }
 
      double ymin () const { return m_ymin; }
 
      double ymax () const { return m_ymax; }    
 
 
      
 
      void free ();
 
      double operator () (const double xx, const double yy) const;
 
      std::vector<double> eval_func (const std::vector<std::vector<double>> xx) const;
 
      double IntegrateVegas (const double xmin, const double xmax, const double ymin, const double ymax) const;
 
      double IntegrateSuave (const double xmin, const double xmax, const double ymin, const double ymax) const;
  
      double IntegrateDivonne (const double xmin, const double xmax, const double ymin, const double ymax) const;
  
      double IntegrateCuhre (const double xmin, const double xmax, const double ymin, const double ymax) const;
 
    };
 
 
  }
}
    
#endif