sf_coulomb.hpp

Go to the documentation of this file.

/*
 * $Id: sf_coulomb.hpp 9 2010-06-13 14:02:43Z jdl3 $
 * Copyright (C) 2010 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_COULOMB_HPP
#define CCGSL_SF_COULOMB_HPP
 
#include<gsl/gsl_sf_coulomb.h>
#include"mode.hpp"
#include"sf_result.hpp"
 
namespace gsl {
  namespace sf {
    inline int hydrogenicR_1_e( double const Z, double const r, result& result ){
      return gsl_sf_hydrogenicR_1_e( Z, r, &result ); }
    inline double hydrogenicR_1( double const Z, double const r ){
      return gsl_sf_hydrogenicR_1( Z, r ); } 
    inline int hydrogenicR_e( int const n, int const l, double const Z, double const r, result& result ){
      return gsl_sf_hydrogenicR_e( n, l, Z, r, &result ); } 
    inline double hydrogenicR( int const n, int const l, double const Z, double const r ){
      return gsl_sf_hydrogenicR( n, l, Z, r ); } 
    namespace coulomb {
#ifndef DOXYGEN_SKIP
      inline int wave_FG_e( double const eta, double const x, double const lam_F, int const k_lam_G,
                result& F, result& Fp, result& G, result& Gp, double* exp_F, double* exp_G ){
    return gsl_sf_coulomb_wave_FG_e( eta, x, lam_F, k_lam_G, &F, &Fp, &G, &Gp, exp_F, exp_G ); } 
#endif // DOXYGEN_SKIP
      inline int wave_FG_e( double const eta, double const x, double const lam_F, int const k_lam_G,
                result& F, result& Fp, result& G, result& Gp, double& exp_F, double& exp_G ){
    return gsl_sf_coulomb_wave_FG_e( eta, x, lam_F, k_lam_G, &F, &Fp, &G, &Gp, &exp_F, &exp_G ); } 
#ifndef DOXYGEN_SKIP
      inline int wave_F_array( double lam_min, int kmax, double eta, double x,
                   double* fc_array, double* F_exponent ){
    return gsl_sf_coulomb_wave_F_array( lam_min, kmax, eta, x, fc_array, F_exponent ); } 
#endif // DOXYGEN_SKIP
      template<typename DATA>
      inline int wave_F_array( double lam_min, int kmax, double eta, double x,
                   DATA& fc_array, double& F_exponent ){
    if(fc_array.size() < static_cast<size_t>(kmax)){ return GSL_EBADLEN;}
    return gsl_sf_coulomb_wave_F_array( lam_min, kmax, eta, x, fc_array.data(),
                        &F_exponent ); } 
#ifndef DOXYGEN_SKIP
      inline int wave_FG_array( double lam_min, int kmax, double eta, double x, double* fc_array,
                double* gc_array, double* F_exponent, double* G_exponent ){
    return gsl_sf_coulomb_wave_FG_array( lam_min, kmax, eta, x, fc_array, gc_array,
                         F_exponent, G_exponent ); } 
#endif // DOXYGEN_SKIP
      template<typename DATA>
      inline int wave_FG_array( double lam_min, int kmax, double eta, double x, DATA& fc_array,
                DATA& gc_array, double& F_exponent, double& G_exponent ){
    if(fc_array.size() < static_cast<size_t>(kmax)){ return GSL_EBADLEN;}
    if(gc_array.size() < static_cast<size_t>(kmax)){ return GSL_EBADLEN;}
    return gsl_sf_coulomb_wave_FG_array( lam_min, kmax, eta, x, fc_array.data(),
                         gc_array.data(), &F_exponent, &G_exponent ); } 
#ifndef DOXYGEN_SKIP
      inline int wave_FGp_array( double lam_min, int kmax, double eta, double x,
                 double* fc_array, double* fcp_array, double* gc_array,
                 double* gcp_array, double* F_exponent,
                 double* G_exponent ){
    return gsl_sf_coulomb_wave_FGp_array( lam_min, kmax, eta, x, fc_array, fcp_array,
                          gc_array, gcp_array, F_exponent, G_exponent ); } 
#endif // DOXYGEN_SKIP
      template<typename DATA>
      inline int wave_FGp_array( double lam_min, int kmax, double eta, double x,
                 DATA& fc_array, DATA& fcp_array, DATA& gc_array,
                 DATA& gcp_array, double& F_exponent,
                 double& G_exponent ){
    if(fc_array.size() < static_cast<size_t>(kmax)){ return GSL_EBADLEN;}
    if(fcp_array.size() < static_cast<size_t>(kmax)){ return GSL_EBADLEN;}
    if(gc_array.size() < static_cast<size_t>(kmax)){ return GSL_EBADLEN;}
    if(gcp_array.size() < static_cast<size_t>(kmax)){ return GSL_EBADLEN;}
    return gsl_sf_coulomb_wave_FGp_array( lam_min, kmax, eta, x, fc_array.data(),
                          fcp_array.data(), gc_array.data(),
                          gcp_array.data(), &F_exponent, &G_exponent ); } 
#ifndef DOXYGEN_SKIP
      inline int wave_sphF_array( double lam_min, int kmax, double eta, double x,
                  double* fc_array, double* F_exponent ){
    return gsl_sf_coulomb_wave_sphF_array( lam_min, kmax, eta, x, fc_array, F_exponent ); } 
#endif // DOXYGEN_SKIP
      template<typename DATA>
      inline int wave_sphF_array( double lam_min, int kmax, double eta, double x,
                  DATA& fc_array, double& F_exponent ){
    if(fc_array.size() < static_cast<size_t>(kmax)){ return GSL_EBADLEN;}
    return gsl_sf_coulomb_wave_sphF_array( lam_min, kmax, eta, x, fc_array.data(),
                           &F_exponent ); } 
      inline int CL_e( double L, double eta, result& result ){
    return gsl_sf_coulomb_CL_e( L, eta, &result ); } 
#ifndef DOXYGEN_SKIP
      inline int CL_array( double Lmin, int kmax, double eta, double* cl ){
    return gsl_sf_coulomb_CL_array( Lmin, kmax, eta, cl ); }
#endif // DOXYGEN_SKIP
      template<typename DATA>
      inline int CL_array( double Lmin, int kmax, double eta, DATA& cl ){
    if(cl.size() < static_cast<size_t>(kmax)){ return GSL_EBADLEN;}
    return gsl_sf_coulomb_CL_array( Lmin, kmax, eta, cl.data() ); }
    }
  }
}
 
#endif