sf_bessel.hpp

Go to the documentation of this file.

/*
 * $Id: sf_bessel.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_BESSEL_HPP
#define CCGSL_SF_BESSEL_HPP
 
#include<gsl/gsl_sf_bessel.h>
#include"mode.hpp"
#include"sf_result.hpp"
 
namespace gsl {
  namespace sf {
    namespace bessel {
      inline int J0_e( double const x, result& result ){ return gsl_sf_bessel_J0_e( x, &result ); }  
      inline double J0( double const x ){ return gsl_sf_bessel_J0( x ); } 
      inline int J1_e( double const x, result& result ){ return gsl_sf_bessel_J1_e( x, &result ); } 
      inline double J1( double const x ){ return gsl_sf_bessel_J1( x ); } 
      inline int Jn_e( int n, double x, result& result ){ return gsl_sf_bessel_Jn_e( n, x, &result ); } 
      inline double Jn( int const n, double const x ){ return gsl_sf_bessel_Jn( n, x ); } 
#ifndef DOXYGEN_SKIP
      inline int Jn_array( int nmin, int nmax, double x, double* result_array ){
    return gsl_sf_bessel_Jn_array( nmin, nmax, x, result_array ); } 
#endif // DOXYGEN_SKIP
      template<typename DATA>
      inline int Jn_array( int nmin, int nmax, double x, DATA& result_array ){
    if(result_array.size() < static_cast<size_t>(nmax - nmin + 1)){
      return GSL_EBADLEN;
    }
    return gsl_sf_bessel_Jn_array( nmin, nmax, x, result_array.data() ); } 
      inline int Y0_e( double const x, result& result ){ return gsl_sf_bessel_Y0_e( x, &result ); } 
      inline double Y0( double const x ){ return gsl_sf_bessel_Y0( x ); } 
      inline int Y1_e( double const x, result& result ){ return gsl_sf_bessel_Y1_e( x, &result ); } 
      inline double Y1( double const x ){ return gsl_sf_bessel_Y1( x ); } 
      inline int Yn_e( int n, double const x, result& result ){ return gsl_sf_bessel_Yn_e( n, x, &result ); } 
      inline double Yn( int const n, double const x ){ return gsl_sf_bessel_Yn( n, x ); } 
#ifndef DOXYGEN_SKIP
      inline int Yn_array( int const nmin, int const nmax, double const x,
               double* result_array ){
    return gsl_sf_bessel_Yn_array( nmin, nmax, x, result_array ); } 
#endif // DOXYGEN_SKIP
      template<typename DATA>
      inline int Yn_array( int const nmin, int const nmax, double const x,
               DATA& result_array ){
    if(result_array.size() < static_cast<size_t>(nmax - nmin + 1)){
      return GSL_EBADLEN;
    }
    return gsl_sf_bessel_Yn_array( nmin, nmax, x, result_array.data() ); } 
      inline int I0_e( double const x, result& result ){ return gsl_sf_bessel_I0_e( x, &result ); } 
      inline double I0( double const x ){ return gsl_sf_bessel_I0( x ); } 
      inline int I1_e( double const x, result& result ){ return gsl_sf_bessel_I1_e( x, &result ); } 
      inline double I1( double const x ){ return gsl_sf_bessel_I1( x ); } 
      inline int In_e( int const n, double const x, result& result ){ return gsl_sf_bessel_In_e( n, x, &result ); } 
      inline double In( int const n, double const x ){ return gsl_sf_bessel_In( n, x ); } 
#ifndef DOXYGEN_SKIP
      inline int In_array( int const nmin, int const nmax, double const x,
               double* result_array ){
    return gsl_sf_bessel_In_array( nmin, nmax, x, result_array ); } 
#endif // DOXYGEN_SKIP
      template<typename DATA>
      inline int In_array( int const nmin, int const nmax, double const x,
               DATA& result_array ){
    if(result_array.size() < static_cast<size_t>(nmax - nmin + 1)){
      return GSL_EBADLEN;
    }
    return gsl_sf_bessel_In_array( nmin, nmax, x, result_array.data() ); } 
      inline int I0_scaled_e( double const x, result& result ){ return gsl_sf_bessel_I0_scaled_e( x, &result ); } 
      inline double I0_scaled( double const x ){ return gsl_sf_bessel_I0_scaled( x ); } 
      inline int I1_scaled_e( double const x, result& result ){ return gsl_sf_bessel_I1_scaled_e( x, &result ); } 
      inline double I1_scaled( double const x ){ return gsl_sf_bessel_I1_scaled( x ); } 
      inline int In_scaled_e( int n, double const x, result& result ){
    return gsl_sf_bessel_In_scaled_e( n, x, &result ); } 
      inline double In_scaled( int const n, double const x ){ return gsl_sf_bessel_In_scaled( n, x ); } 
#ifndef DOXYGEN_SKIP
      inline int In_scaled_array( int const nmin, int const nmax, double const x,
                  double* result_array ){
    return gsl_sf_bessel_In_scaled_array( nmin, nmax, x, result_array ); } 
#endif // DOXYGEN_SKIP
      template<typename DATA>
      inline int In_scaled_array( int const nmin, int const nmax, double const x,
                  DATA& result_array ){
    if(result_array.size() < static_cast<size_t>(nmax - nmin + 1)){
      return GSL_EBADLEN;
    }
    return gsl_sf_bessel_In_scaled_array( nmin, nmax, x, result_array.data() ); } 
      inline int K0_e( double const x, result& result ){ return gsl_sf_bessel_K0_e( x, &result ); } 
      inline double K0( double const x ){ return gsl_sf_bessel_K0( x ); } 
      inline int K1_e( double const x, result& result ){ return gsl_sf_bessel_K1_e( x, &result ); } 
      inline double K1( double const x ){ return gsl_sf_bessel_K1( x ); } 
      inline int Kn_e( int const n, double const x, result& result ){
    return gsl_sf_bessel_Kn_e( n, x, &result ); } 
      inline double Kn( int const n, double const x ){ return gsl_sf_bessel_Kn( n, x ); } 
#ifndef DOXYGEN_SKIP
      inline int Kn_array( int const nmin, int const nmax, double const x,
               double* result_array ){
    return gsl_sf_bessel_Kn_array( nmin, nmax, x, result_array ); } 
#endif // DOXYGEN_SKIP
      template<typename DATA>
      inline int Kn_array( int const nmin, int const nmax, double const x,
               DATA& result_array ){
    if(result_array.size() < static_cast<size_t>(nmax - nmin + 1)){
      return GSL_EBADLEN;
    }
    return gsl_sf_bessel_Kn_array( nmin, nmax, x, result_array.data() ); } 
      inline int K0_scaled_e( double const x, result& result ){
    return gsl_sf_bessel_K0_scaled_e( x, &result ); } 
      inline double K0_scaled( double const x ){ return gsl_sf_bessel_K0_scaled( x ); } 
      inline int K1_scaled_e( double const x, result& result ){
    return gsl_sf_bessel_K1_scaled_e( x, &result ); } 
      inline double K1_scaled( double const x ){ return gsl_sf_bessel_K1_scaled( x ); } 
      inline int Kn_scaled_e( int n, double const x, result& result ){
    return gsl_sf_bessel_Kn_scaled_e( n, x, &result ); } 
      inline double Kn_scaled( int const n, double const x ){ return gsl_sf_bessel_Kn_scaled( n, x ); } 
      template<typename DATA>
      inline int Kn_scaled_array( int const nmin, int const nmax, double const x,
                  DATA& result_array ){
    if(result_array.size() < static_cast<size_t>(nmax - nmin + 1)){
      return GSL_EBADLEN;
    }
    return gsl_sf_bessel_Kn_scaled_array( nmin, nmax, x, result_array.data() ); } 
      inline int j0_e( double const x, result& result ){ return gsl_sf_bessel_j0_e( x, &result ); } 
      inline double j0( double const x ){ return gsl_sf_bessel_j0( x ); } 
      inline int j1_e( double const x, result& result ){ return gsl_sf_bessel_j1_e( x, &result ); } 
      inline double j1( double const x ){ return gsl_sf_bessel_j1( x ); } 
      inline int j2_e( double const x, result& result ){ return gsl_sf_bessel_j2_e( x, &result ); } 
      inline double j2( double const x ){ return gsl_sf_bessel_j2( x ); } 
      inline int jl_e( int const l, double const x, result& result ){
    return gsl_sf_bessel_jl_e( l, x, &result ); } 
      inline double jl( int const l, double const x ){ return gsl_sf_bessel_jl( l, x ); } 
#ifndef DOXYGEN_SKIP
      inline int jl_array( int const lmax, double const x, double* result_array ){
    return gsl_sf_bessel_jl_array( lmax, x, result_array ); } 
#endif // DOXYGEN_SKIP
      template<typename DATA>
      inline int jl_array( int const lmax, double const x, DATA& result_array ){
    if(result_array.size() < static_cast<size_t>(lmax + 1)){
      return GSL_EBADLEN;
    }
    return gsl_sf_bessel_jl_array( lmax, x, result_array.data() ); } 
#ifndef DOXYGEN_SKIP
      inline int jl_steed_array( int const lmax, double const x, double* jl_x_array ){
    return gsl_sf_bessel_jl_steed_array( lmax, x, jl_x_array ); } 
#endif // DOXYGEN_SKIP
      template<typename DATA>
      inline int jl_steed_array( int const lmax, double const x, DATA& jl_x_array ){
    if(jl_x_array.size() < static_cast<size_t>(lmax + 1)){
      return GSL_EBADLEN;
    }
    return gsl_sf_bessel_jl_steed_array( lmax, x, jl_x_array.data() ); } 
      inline int y0_e( double const x, result& result ){
    return gsl_sf_bessel_y0_e( x, &result ); } 
      inline double y0( double const x ){ return gsl_sf_bessel_y0( x ); } 
      inline int y1_e( double const x, result& result ){ return gsl_sf_bessel_y1_e( x, &result ); } 
      inline double y1( double const x ){ return gsl_sf_bessel_y1( x ); } 
      inline int y2_e( double const x, result& result ){ return gsl_sf_bessel_y2_e( x, &result ); } 
      inline double y2( double const x ){ return gsl_sf_bessel_y2( x ); } 
      inline int yl_e( int l, double const x, result& result ){ return gsl_sf_bessel_yl_e( l, x, &result ); } 
      inline double yl( int const l, double const x ){ return gsl_sf_bessel_yl( l, x ); } 
#ifndef DOXYGEN_SKIP
      inline int yl_array( int const lmax, double const x, double* result_array ){
    return gsl_sf_bessel_yl_array( lmax, x, result_array ); } 
#endif // DOXYGEN_SKIP
      template<typename DATA>
      inline int yl_array( int const lmax, double const x, DATA& result_array ){
    if(result_array.size() < static_cast<size_t>(lmax + 1)){
      return GSL_EBADLEN;
    }
    return gsl_sf_bessel_yl_array( lmax, x, result_array.data() ); } 
      inline int i0_scaled_e( double const x, result& result ){
    return gsl_sf_bessel_i0_scaled_e( x, &result ); } 
      inline double i0_scaled( double const x ){ return gsl_sf_bessel_i0_scaled( x ); } 
      inline int i1_scaled_e( double const x, result& result ){ return gsl_sf_bessel_i1_scaled_e( x, &result ); } 
      inline double i1_scaled( double const x ){ return gsl_sf_bessel_i1_scaled( x ); } 
      inline int i2_scaled_e( double const x, result& result ){ return gsl_sf_bessel_i2_scaled_e( x, &result ); } 
      inline double i2_scaled( double const x ){ return gsl_sf_bessel_i2_scaled( x ); } 
      inline int il_scaled_e( int const l, double x, result& result ){
    return gsl_sf_bessel_il_scaled_e( l, x, &result ); } 
      inline double il_scaled( int const l, double const x ){ return gsl_sf_bessel_il_scaled( l, x ); } 
#ifndef DOXYGEN_SKIP
      inline int il_scaled_array( int const lmax, double const x, double* result_array ){
    return gsl_sf_bessel_il_scaled_array( lmax, x, result_array ); } 
#endif // DOXYGEN_SKIP
      template<typename DATA>
      inline int il_scaled_array( int const lmax, double const x, DATA& result_array ){
    if( result_array.size() < static_cast<size_t>(lmax + 1)) {
      return GSL_EBADLEN;
    }
    return gsl_sf_bessel_il_scaled_array( lmax, x, result_array.data() ); } 
      inline int k0_scaled_e( double const x, result& result ){
    return gsl_sf_bessel_k0_scaled_e( x, &result ); } 
      inline double k0_scaled( double const x ){ return gsl_sf_bessel_k0_scaled( x ); } 
      inline int k1_scaled_e( double const x, result& result ){
    return gsl_sf_bessel_k1_scaled_e( x, &result ); } 
      inline double k1_scaled( double const x ){ return gsl_sf_bessel_k1_scaled( x ); } 
      inline int k2_scaled_e( double const x, result& result ){
    return gsl_sf_bessel_k2_scaled_e( x, &result ); } 
      inline double k2_scaled( double const x ){ return gsl_sf_bessel_k2_scaled( x ); } 
      inline int kl_scaled_e( int l, double const x, result& result ){
    return gsl_sf_bessel_kl_scaled_e( l, x, &result ); } 
      inline double kl_scaled( int const l, double const x ){ return gsl_sf_bessel_kl_scaled( l, x ); } 
#ifndef DOXYGEN_SKIP
      inline int kl_scaled_array( int const lmax, double const x, double* result_array ){
    return gsl_sf_bessel_kl_scaled_array( lmax, x, result_array ); } 
#endif // DOXYGEN_SKIP
      template<typename DATA>
      inline int kl_scaled_array( int const lmax, double const x, DATA& result_array ){
    if(result_array.size() < static_cast<size_t>(lmax + 1)){
      return GSL_EBADLEN;
    }
    return gsl_sf_bessel_kl_scaled_array( lmax, x, result_array.data() ); } 
      inline int Jnu_e( double const nu, double const x, result& result ){
    return gsl_sf_bessel_Jnu_e( nu, x, &result ); } 
      inline double Jnu( double const nu, double const x ){ return gsl_sf_bessel_Jnu( nu, x ); } 
      inline int Ynu_e( double nu, double x, result& result ){
    return gsl_sf_bessel_Ynu_e( nu, x, &result ); } 
      inline double Ynu( double const nu, double const x ){ return gsl_sf_bessel_Ynu( nu, x ); } 
#ifndef DOXYGEN_SKIP
      inline int sequence_Jnu_e( double nu, mode_t mode, size_t size, double* v ){
    return gsl_sf_bessel_sequence_Jnu_e( nu, mode, size, v ); } 
#endif // DOXYGEN_SKIP
      template<typename DATA>
      inline int sequence_Jnu_e( double nu, mode_t mode, DATA& v ){
    return gsl_sf_bessel_sequence_Jnu_e( nu, mode, v.size(), v.data() ); } 
      inline int Inu_scaled_e( double nu, double x, result& result ){
    return gsl_sf_bessel_Inu_scaled_e( nu, x, &result ); } 
      inline double Inu_scaled( double nu, double x ){ return gsl_sf_bessel_Inu_scaled( nu, x ); } 
      inline int Inu_e( double nu, double x, result& result ){
    return gsl_sf_bessel_Inu_e( nu, x, &result ); } 
      inline double Inu( double nu, double x ){ return gsl_sf_bessel_Inu( nu, x ); } 
      inline int Knu_scaled_e( double const nu, double const x, result& result ){
    return gsl_sf_bessel_Knu_scaled_e( nu, x, &result ); } 
      inline double Knu_scaled( double const nu, double const x ){
    return gsl_sf_bessel_Knu_scaled( nu, x ); } 
      inline int Knu_e( double const nu, double const x, result& result ){
    return gsl_sf_bessel_Knu_e( nu, x, &result ); } 
      inline double Knu( double const nu, double const x ){ return gsl_sf_bessel_Knu( nu, x ); } 
      inline int lnKnu_e( double const nu, double const x, result& result ){
    return gsl_sf_bessel_lnKnu_e( nu, x, &result ); } 
      inline double lnKnu( double const nu, double const x ){
    return gsl_sf_bessel_lnKnu( nu, x ); } 
      inline int zero_J0_e( unsigned int s, result& result ){
    return gsl_sf_bessel_zero_J0_e( s, &result ); } 
      inline double zero_J0( unsigned int s ){ return gsl_sf_bessel_zero_J0( s ); } 
      inline int zero_J1_e( unsigned int s, result& result ){
    return gsl_sf_bessel_zero_J1_e( s, &result ); } 
      inline double zero_J1( unsigned int s ){ return gsl_sf_bessel_zero_J1( s ); } 
      inline int zero_Jnu_e( double nu, unsigned int s, result& result ){
    return gsl_sf_bessel_zero_Jnu_e( nu, s, &result ); } 
      inline double zero_Jnu( double nu, unsigned int s ){
    return gsl_sf_bessel_zero_Jnu( nu, s ); }
    }
  }
}
 
#endif