LegendrePolynomials.h

Go to the documentation of this file.

/*******************************************************************
 *  Copyright (C) 2010 by Federico Marulli                         *
 *  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 __LEGPOLS__
#define __LEGPOLS__
 
#include "Func.h"
 
namespace cbl {
 
  namespace glob {
 
    class LegendrePolynomials {
 
    protected:
 
      size_t m_nOrders;
 
      Eigen::MatrixXd m_coefficients;
 
      void m_set_coefficients (const int lMax);
    
    public:
 
 
      LegendrePolynomials () { set(0); }
 
      LegendrePolynomials (const int lMax, const bool safe=false);
 
      ~LegendrePolynomials () {}
 
    
      void set (const int lMax);
 
      double operator () (const double x, const int ell);
 
      std::vector<double> operator () (const double x);
 
      double integral (const double x_min, const double x_max, const int ell);
 
      std::vector<double> integral (const double x_min, const double x_max);
 
      std::vector<std::vector<double>> triangle (const double r12, const double r13, const double r23);
 
      std::vector<double> triangle_integral (const double r12_min, const double r12_max, const double r13_min, const double r13_max, const double r23_min, const double r23_max, const double rel_err=1.e-4, const int nevals=10000);
    };
  }
}
 
#endif