sf_gegenbauer.hpp

Go to the documentation of this file.

/*
 * $Id: sf_gegenbauer.hpp 9 2010-06-13 14:02:43Z jdl3 $
 * Copyright (C) 2010, 2020 John D Lamb
 *
 * 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., 675 Mass Ave, Cambridge, MA 02139, USA.
 */
 
#ifndef CCGSL_SF_GEGENBAUER_HPP
#define CCGSL_SF_GEGENBAUER_HPP
 
#include<gsl/gsl_sf_gegenbauer.h>
#include"mode.hpp"
#include"sf_result.hpp"
 
namespace gsl {
  namespace sf {
    inline int gegenpoly_1_e( double lambda, double x, result& result ){
      return gsl_sf_gegenpoly_1_e( lambda, x, &result ); } 
    inline int gegenpoly_2_e( double lambda, double x, result& result ){
      return gsl_sf_gegenpoly_2_e( lambda, x, &result ); } 
    inline int gegenpoly_3_e( double lambda, double x, result& result ){
      return gsl_sf_gegenpoly_3_e( lambda, x, &result ); } 
    inline double gegenpoly_1( double lambda, double x ){
      return gsl_sf_gegenpoly_1( lambda, x ); } 
    inline double gegenpoly_2( double lambda, double x ){
      return gsl_sf_gegenpoly_2( lambda, x ); } 
    inline double gegenpoly_3( double lambda, double x ){
      return gsl_sf_gegenpoly_3( lambda, x ); } 
    inline int gegenpoly_n_e( int n, double lambda, double x, result& result ){
      return gsl_sf_gegenpoly_n_e( n, lambda, x, &result ); } 
    inline double gegenpoly_n( int n, double lambda, double x ){
      return gsl_sf_gegenpoly_n( n, lambda, x ); } 
#ifndef DOXYGEN_SKIP
    inline int gegenpoly_array( int nmax, double lambda, double x, double* result_array ){
      return gsl_sf_gegenpoly_array( nmax, lambda, x, result_array );
    }
#endif // DOXYGEN_SKIP
    template<typename DATA>
    inline int gegenpoly_array( int nmax, double lambda, double x, DATA& result_array ){
      if(result_array.size() < static_cast<size_t>(nmax + 1)){ return GSL_EBADLEN;}
      return gsl_sf_gegenpoly_array( nmax, lambda, x, result_array.data() );
    }
  }
}
 
#endif