Skip to content
Snippets Groups Projects
Select Git revision
  • 7b1ed83f18dc9505ed9557a30edb296c089ee775
  • master default protected
  • dof-renumbering
  • gdemesy-master-patch-30528
  • eval-space-time
  • oscillating_multiharm
  • MH_movement
  • axisqu
  • write_vtu_and_ensight_formats
  • movingband
  • CP_1972_add_vtu_file_writing
  • mortar
  • fast_freq_sweep_Resolution
  • applyresolvent_again
  • marteaua-master-patch-54323
  • patch-1
  • binde-master-patch-08072
  • binde-master-patch-52461
  • BCGSL
  • resolvent
  • TreeElementsOf
  • getdp_3_5_0
  • getdp_3_4_0
  • getdp_3_3_0
  • getdp_3_2_0
  • getdp_3_1_0
  • getdp_3_0_4
  • getdp_3_0_3
  • getdp_3_0_2
  • getdp_3_0_1
  • getdp_3_0_0
  • onelab_mobile_2.1.0
  • getdp_2_11_3 protected
  • getdp_2_11_2 protected
  • getdp_2_11_1 protected
  • getdp_2_11_0 protected
  • getdp_2_10_0 protected
  • getdp_2_9_2 protected
  • getdp_2_9_1 protected
  • getdp_2_9_0 protected
  • getdp_2_8_0 protected
41 results

Bessel.h

Blame
  • Bessel.h 1.41 KiB
    // GetDP - Copyright (C) 1997-2015 P. Dular and C. Geuzaine, University of Liege
    //
    // See the LICENSE.txt file for license information. Please report all
    // bugs and problems to the public mailing list <getdp@geuz.org>.
    
    #ifndef _BESSEL_H_
    #define _BESSEL_H_
    
    #include <math.h>
    #include <complex>
    
    #define BESSEL_ERROR_INPUT     1
    #define BESSEL_OVERFLOW        2
    #define BESSEL_HALF_ACCURACY   3
    #define BESSEL_NO_ACCURACY     4
    #define BESSEL_NO_CONVERGENCE  5
    
    // These routines provide a C++ interface to the Fortran Bessel
    // functions from Donald E. Amos (Sandia National Laboratories)
    
    int BesselJn(double n, int num, double x, double *val);
    int BesselSphericalJn(double n, int num, double x, double *val);
    int BesselAltSphericalJn(double n, int num, double x, double *val);
    
    int BesselYn(double n, int num, double x, double *val);
    int BesselSphericalYn(double n, int num, double x, double *val);
    int BesselAltSphericalYn(double n, int num, double x, double *val);
    
    int BesselHn(int type, double n, int num, double x, std::complex<double> *val);
    int BesselSphericalHn(int type, double n, int num, double x, std::complex<double> *val);
    int BesselAltSphericalHn(int type, double n, int num, double x, std::complex<double> *val);
    
    // Utilities for backward compatibility
    
    double Spherical_j_n(int n, double x);
    double AltSpherical_j_n(int n, double x);
    double Spherical_y_n(int n, double x);
    double AltSpherical_y_n(int n, double x);
    
    #endif