/home66/gary/public_html/cloudy/c08_branch/source/thirdparty.cpp File Reference

#include "cddefines.h"
#include "thirdparty.h"
#include "physconst.h"

Include dependency graph for thirdparty.cpp:

Go to the source code of this file.

Data Structures

class  t_lfact

Defines

#define N   624
#define M   397
#define MATRIX_A   0x9908b0dfUL
#define UMASK   0x80000000UL
#define LMASK   0x7fffffffUL
#define MIXBITS(u, v)   ( ((u) & UMASK) | ((v) & LMASK) )
#define TWIST(u, v)   ((MIXBITS(u,v) >> 1) ^ ((v)&1UL ? MATRIX_A : 0UL))

Functions

double polevl (double x, const double coef[], int N)
double p1evl (double x, const double coef[], int N)
double chbevl (double, const double[], int)
bool linfit (long n, double x[], double y[], double &a, double &siga, double &b, double &sigb)
double factorial (long n)
double lfactorial (long n)
complex< double > cdgamma (complex< double > x)
double bessel_j0 (double x)
double bessel_y0 (double x)
double bessel_j1 (double x)
double bessel_y1 (double x)
double bessel_jn (int n, double x)
double bessel_yn (int n, double x)
double bessel_k0 (double x)
double bessel_k0_scaled (double x)
double bessel_k1 (double x)
double bessel_k1_scaled (double x)
double bessel_i0 (double x)
double bessel_i0_scaled (double x)
double bessel_i1 (double x)
double bessel_i1_scaled (double x)
double ellpk (double x)
double expn (int n, double x)
void init_genrand (unsigned long s)
void init_by_array (unsigned long init_key[], int key_length)
static void next_state (void)
unsigned long genrand_int32 (void)
long genrand_int31 (void)
double genrand_real1 (void)
double genrand_real2 (void)
double genrand_real3 (void)
double genrand_res53 (void)

Variables

static const double pre_factorial [NPRE_FACTORIAL]
static const double b0_PP [7]
static const double b0_PQ [7]
static const double b0_QP [8]
static const double b0_QQ [7]
static const double b0_YP [8]
static const double b0_YQ [7]
static const double DR1 = 5.78318596294678452118e0
static const double DR2 = 3.04712623436620863991e1
static double b0_RP [4]
static double b0_RQ [8]
static const double TWOOPI = 2./PI
static const double SQ2OPI = sqrt(2./PI)
static const double PIO4 = PI/4.
static const double b1_RP [4]
static const double b1_RQ [8]
static const double b1_PP [7]
static const double b1_PQ [7]
static const double b1_QP [8]
static const double b1_QQ [7]
static const double b1_YP [6]
static const double b1_YQ [8]
static const double Z1 = 1.46819706421238932572E1
static const double Z2 = 4.92184563216946036703E1
static const double THPIO4 = 3.*PI/4.
static const double k0_A []
static const double k0_B []
static const double k1_A []
static const double k1_B []
static const double i0_A []
static const double i0_B []
static double i1_A []
static double i1_B []
static const double elk_P []
static const double elk_Q []
static const double C1 = 1.3862943611198906188e0
static const double MAXLOG = log(DBL_MAX)
static const double BIG = 1.44115188075855872E+17
static unsigned long state [N]
static int nleft = 1
static int initf = 0
static unsigned long * next


Define Documentation

#define LMASK   0x7fffffffUL

Definition at line 2455 of file thirdparty.cpp.

#define M   397

Definition at line 2452 of file thirdparty.cpp.

Referenced by iso_state_lifetime(), and next_state().

#define MATRIX_A   0x9908b0dfUL

Definition at line 2453 of file thirdparty.cpp.

#define MIXBITS ( u,
 )     ( ((u) & UMASK) | ((v) & LMASK) )

Definition at line 2456 of file thirdparty.cpp.

#define N   624

Definition at line 2451 of file thirdparty.cpp.

#define TWIST ( u,
 )     ((MIXBITS(u,v) >> 1) ^ ((v)&1UL ? MATRIX_A : 0UL))

Definition at line 2457 of file thirdparty.cpp.

Referenced by next_state().

#define UMASK   0x80000000UL

Definition at line 2454 of file thirdparty.cpp.


Function Documentation

double bessel_i0 ( double  x  ) 

Definition at line 1700 of file thirdparty.cpp.

References chbevl(), DEBUG_ENTRY, i0_A, and i0_B.

Referenced by bessel_k0(), and bessel_k0_scaled().

Here is the call graph for this function:

double bessel_i0_scaled ( double  x  ) 

Definition at line 1717 of file thirdparty.cpp.

References chbevl(), DEBUG_ENTRY, i0_A, and i0_B.

Here is the call graph for this function:

double bessel_i1 ( double  x  ) 

Definition at line 1882 of file thirdparty.cpp.

References chbevl(), DEBUG_ENTRY, i1_A, and i1_B.

Referenced by bessel_k1(), and bessel_k1_scaled().

Here is the call graph for this function:

double bessel_i1_scaled ( double  x  ) 

Definition at line 1903 of file thirdparty.cpp.

References chbevl(), DEBUG_ENTRY, i1_A, and i1_B.

Here is the call graph for this function:

double bessel_j0 ( double  x  ) 

Definition at line 683 of file thirdparty.cpp.

References b0_PP, b0_PQ, b0_QP, b0_QQ, b0_RP, b0_RQ, DEBUG_ENTRY, DR1, DR2, p1evl(), PIO4, polevl(), and SQ2OPI.

Referenced by bessel_jn(), and bessel_y0().

Here is the call graph for this function:

double bessel_j1 ( double  x  ) 

Definition at line 914 of file thirdparty.cpp.

References b1_PP, b1_PQ, b1_QP, b1_QQ, b1_RP, b1_RQ, DEBUG_ENTRY, p1evl(), polevl(), SQ2OPI, THPIO4, Z1, and Z2.

Referenced by bessel_jn(), and bessel_y1().

Here is the call graph for this function:

double bessel_jn ( int  n,
double  x 
)

Definition at line 1017 of file thirdparty.cpp.

References bessel_j0(), bessel_j1(), DEBUG_ENTRY, and sign().

Here is the call graph for this function:

double bessel_k0 ( double  x  ) 

Definition at line 1333 of file thirdparty.cpp.

References bessel_i0(), cdEXIT, chbevl(), DEBUG_ENTRY, ioQQQ, k0_A, and k0_B.

Referenced by S62_Therm_ave_coll_str().

Here is the call graph for this function:

double bessel_k0_scaled ( double  x  ) 

Definition at line 1356 of file thirdparty.cpp.

References bessel_i0(), cdEXIT, chbevl(), DEBUG_ENTRY, ioQQQ, k0_A, and k0_B.

Here is the call graph for this function:

double bessel_k1 ( double  x  ) 

Definition at line 1509 of file thirdparty.cpp.

References bessel_i1(), cdEXIT, chbevl(), DEBUG_ENTRY, ioQQQ, k1_A, and k1_B.

Referenced by S62_Therm_ave_coll_str().

Here is the call graph for this function:

double bessel_k1_scaled ( double  x  ) 

Definition at line 1531 of file thirdparty.cpp.

References bessel_i1(), cdEXIT, chbevl(), DEBUG_ENTRY, ioQQQ, k1_A, and k1_B.

Here is the call graph for this function:

double bessel_y0 ( double  x  ) 

Definition at line 721 of file thirdparty.cpp.

References b0_PP, b0_PQ, b0_QP, b0_QQ, b0_YP, b0_YQ, bessel_j0(), cdEXIT, DEBUG_ENTRY, ioQQQ, p1evl(), PIO4, polevl(), SQ2OPI, and TWOOPI.

Referenced by bessel_yn().

Here is the call graph for this function:

double bessel_y1 ( double  x  ) 

Definition at line 941 of file thirdparty.cpp.

References b1_PP, b1_PQ, b1_QP, b1_QQ, b1_YP, b1_YQ, bessel_j1(), cdEXIT, DEBUG_ENTRY, ioQQQ, p1evl(), polevl(), SQ2OPI, THPIO4, and TWOOPI.

Referenced by bessel_yn().

Here is the call graph for this function:

double bessel_yn ( int  n,
double  x 
)

Definition at line 1151 of file thirdparty.cpp.

References bessel_y0(), bessel_y1(), cdEXIT, DEBUG_ENTRY, ioQQQ, and sign().

Here is the call graph for this function:

complex<double> cdgamma ( complex< double >  x  ) 

Definition at line 411 of file thirdparty.cpp.

References DEBUG_ENTRY.

Referenced by DoSutherland(), and Hypergeometric2F1().

double chbevl ( double  ,
const   double[],
int   
) [inline]

Definition at line 2376 of file thirdparty.cpp.

References b0, b1, and b2.

Referenced by bessel_i0(), bessel_i0_scaled(), bessel_i1(), bessel_i1_scaled(), bessel_k0(), bessel_k0_scaled(), bessel_k1(), and bessel_k1_scaled().

double ellpk ( double  x  ) 

Definition at line 2015 of file thirdparty.cpp.

References C1, cdEXIT, DEBUG_ENTRY, elk_P, elk_Q, ioQQQ, and polevl().

Referenced by StarkCollTransProb_VF01().

Here is the call graph for this function:

double expn ( int  n,
double  x 
)

expn, returns exponential integral,

Parameters:
n is order, 1 for first integral integral
x is argument, must be positive

Definition at line 2094 of file thirdparty.cpp.

References BIG, cdEXIT, DEBUG_ENTRY, EULER, factorial(), ioQQQ, is_odd(), MAXLOG, and powi().

Here is the call graph for this function:

double factorial ( long  n  ) 

factorial: compute n! by lookup in table of predefined factorials

Definition at line 340 of file thirdparty.cpp.

References cdEXIT, DEBUG_ENTRY, ioQQQ, NPRE_FACTORIAL, and pre_factorial.

Referenced by bhg(), bhG(), DoSutherland(), expn(), fsff(), and hrii().

long genrand_int31 ( void   ) 

Definition at line 2549 of file thirdparty.cpp.

References next, next_state(), and nleft.

Here is the call graph for this function:

unsigned long genrand_int32 ( void   ) 

Definition at line 2532 of file thirdparty.cpp.

References next, next_state(), and nleft.

Referenced by genrand_res53().

Here is the call graph for this function:

double genrand_real1 ( void   ) 

Definition at line 2566 of file thirdparty.cpp.

References next, next_state(), and nleft.

Here is the call graph for this function:

double genrand_real2 ( void   ) 

Definition at line 2584 of file thirdparty.cpp.

References next, next_state(), and nleft.

Here is the call graph for this function:

double genrand_real3 ( void   ) 

Definition at line 2602 of file thirdparty.cpp.

References next, next_state(), and nleft.

Referenced by RandGauss().

Here is the call graph for this function:

double genrand_res53 ( void   ) 

Definition at line 2620 of file thirdparty.cpp.

References genrand_int32().

Here is the call graph for this function:

void init_by_array ( unsigned long  init_key[],
int  key_length 
)

Definition at line 2484 of file thirdparty.cpp.

References init_genrand(), initf, N, nleft, and state.

Here is the call graph for this function:

void init_genrand ( unsigned long  s  ) 

Definition at line 2465 of file thirdparty.cpp.

References initf, N, nleft, and state.

Referenced by init_by_array(), next_state(), ParseAtomH2(), ParseAtomISO(), and ParseSet().

double lfactorial ( long  n  ) 

lfactorial: compute log10(n!), this sroutine cahes its results for efficiency

Definition at line 383 of file thirdparty.cpp.

References cdEXIT, DEBUG_ENTRY, Singleton< T >::Inst(), and ioQQQ.

Referenced by bhg_log(), bhG_mx(), hrii_log(), and log10_fsff().

Here is the call graph for this function:

bool linfit ( long  n,
double  x[],
double  y[],
double &  a,
double &  siga,
double &  b,
double &  sigb 
)

Definition at line 45 of file thirdparty.cpp.

References DEBUG_ENTRY, and pow2().

Referenced by ConvEdenIoniz().

Here is the call graph for this function:

static void next_state ( void   )  [static]

Definition at line 2510 of file thirdparty.cpp.

References init_genrand(), initf, M, N, next, nleft, state, and TWIST.

Referenced by genrand_int31(), genrand_int32(), genrand_real1(), genrand_real2(), and genrand_real3().

Here is the call graph for this function:

double p1evl ( double  x,
const double  coef[],
int  N 
) [inline]

Definition at line 2302 of file thirdparty.cpp.

Referenced by bessel_j0(), bessel_j1(), bessel_y0(), and bessel_y1().

double polevl ( double  x,
const double  coef[],
int  N 
) [inline]

Definition at line 2280 of file thirdparty.cpp.

Referenced by bessel_j0(), bessel_j1(), bessel_y0(), bessel_y1(), and ellpk().


Variable Documentation

const double b0_PP[7] [static]

Initial value:

 {
        7.96936729297347051624e-4,
        8.28352392107440799803e-2,
        1.23953371646414299388e0,
        5.44725003058768775090e0,
        8.74716500199817011941e0,
        5.30324038235394892183e0,
        9.99999999999999997821e-1,
}

Definition at line 591 of file thirdparty.cpp.

Referenced by bessel_j0(), and bessel_y0().

const double b0_PQ[7] [static]

Initial value:

 {
        9.24408810558863637013e-4,
        8.56288474354474431428e-2,
        1.25352743901058953537e0,
        5.47097740330417105182e0,
        8.76190883237069594232e0,
        5.30605288235394617618e0,
        1.00000000000000000218e0,
}

Definition at line 601 of file thirdparty.cpp.

Referenced by bessel_j0(), and bessel_y0().

const double b0_QP[8] [static]

Initial value:

 {
        -1.13663838898469149931e-2,
        -1.28252718670509318512e0,
        -1.95539544257735972385e1,
        -9.32060152123768231369e1,
        -1.77681167980488050595e2,
        -1.47077505154951170175e2,
        -5.14105326766599330220e1,
        -6.05014350600728481186e0,
}

Definition at line 611 of file thirdparty.cpp.

Referenced by bessel_j0(), and bessel_y0().

const double b0_QQ[7] [static]

Initial value:

 {
        
        6.43178256118178023184e1,
        8.56430025976980587198e2,
        3.88240183605401609683e3,
        7.24046774195652478189e3,
        5.93072701187316984827e3,
        2.06209331660327847417e3,
        2.42005740240291393179e2,
}

Definition at line 622 of file thirdparty.cpp.

Referenced by bessel_j0(), and bessel_y0().

double b0_RP[4] [static]

Initial value:

 {
        -4.79443220978201773821e9,
        1.95617491946556577543e12,
        -2.49248344360967716204e14,
        9.70862251047306323952e15,
}

Definition at line 660 of file thirdparty.cpp.

Referenced by bessel_j0().

double b0_RQ[8] [static]

Initial value:

 {
        
        4.99563147152651017219e2,
        1.73785401676374683123e5,
        4.84409658339962045305e7,
        1.11855537045356834862e10,
        2.11277520115489217587e12,
        3.10518229857422583814e14,
        3.18121955943204943306e16,
        1.71086294081043136091e18,
}

Definition at line 667 of file thirdparty.cpp.

Referenced by bessel_j0().

const double b0_YP[8] [static]

Initial value:

 {
        1.55924367855235737965e4,
        -1.46639295903971606143e7,
        5.43526477051876500413e9,
        -9.82136065717911466409e11,
        8.75906394395366999549e13,
        -3.46628303384729719441e15,
        4.42733268572569800351e16,
        -1.84950800436986690637e16,
}

Definition at line 633 of file thirdparty.cpp.

Referenced by bessel_y0().

const double b0_YQ[7] [static]

Initial value:

 {
        
        1.04128353664259848412e3,
        6.26107330137134956842e5,
        2.68919633393814121987e8,
        8.64002487103935000337e10,
        2.02979612750105546709e13,
        3.17157752842975028269e15,
        2.50596256172653059228e17,
}

Definition at line 644 of file thirdparty.cpp.

Referenced by bessel_y0().

const double b1_PP[7] [static]

Initial value:

 {
        7.62125616208173112003e-4,
        7.31397056940917570436e-2,
        1.12719608129684925192e0,
        5.11207951146807644818e0,
        8.42404590141772420927e0,
        5.21451598682361504063e0,
        1.00000000000000000254e0,
}

Definition at line 846 of file thirdparty.cpp.

Referenced by bessel_j1(), and bessel_y1().

const double b1_PQ[7] [static]

Initial value:

 {
        5.71323128072548699714e-4,
        6.88455908754495404082e-2,
        1.10514232634061696926e0,
        5.07386386128601488557e0,
        8.39985554327604159757e0,
        5.20982848682361821619e0,
        9.99999999999999997461e-1,
}

Definition at line 856 of file thirdparty.cpp.

Referenced by bessel_j1(), and bessel_y1().

const double b1_QP[8] [static]

Initial value:

 {
        5.10862594750176621635e-2,
        4.98213872951233449420e0,
        7.58238284132545283818e1,
        3.66779609360150777800e2,
        7.10856304998926107277e2,
        5.97489612400613639965e2,
        2.11688757100572135698e2,
        2.52070205858023719784e1,
}

Definition at line 866 of file thirdparty.cpp.

Referenced by bessel_j1(), and bessel_y1().

const double b1_QQ[7] [static]

Initial value:

 {
        
        7.42373277035675149943e1,
        1.05644886038262816351e3,
        4.98641058337653607651e3,
        9.56231892404756170795e3,
        7.99704160447350683650e3,
        2.82619278517639096600e3,
        3.36093607810698293419e2,
}

Definition at line 877 of file thirdparty.cpp.

Referenced by bessel_j1(), and bessel_y1().

const double b1_RP[4] [static]

Initial value:

 {
        -8.99971225705559398224e8,
        4.52228297998194034323e11,
        -7.27494245221818276015e13,
        3.68295732863852883286e15,
}

Definition at line 827 of file thirdparty.cpp.

Referenced by bessel_j1().

const double b1_RQ[8] [static]

Initial value:

 {
        
        6.20836478118054335476e2,
        2.56987256757748830383e5,
        8.35146791431949253037e7,
        2.21511595479792499675e10,
        4.74914122079991414898e12,
        7.84369607876235854894e14,
        8.95222336184627338078e16,
        5.32278620332680085395e18,
}

Definition at line 834 of file thirdparty.cpp.

Referenced by bessel_j1().

const double b1_YP[6] [static]

Initial value:

 {
        1.26320474790178026440e9,
        -6.47355876379160291031e11,
        1.14509511541823727583e14,
        -8.12770255501325109621e15,
        2.02439475713594898196e17,
        -7.78877196265950026825e17,
}

Definition at line 888 of file thirdparty.cpp.

Referenced by bessel_y1().

const double b1_YQ[8] [static]

Initial value:

 {
        
        5.94301592346128195359E2,
        2.35564092943068577943E5,
        7.34811944459721705660E7,
        1.87601316108706159478E10,
        3.88231277496238566008E12,
        6.20557727146953693363E14,
        6.87141087355300489866E16,
        3.97270608116560655612E18,
}

Definition at line 897 of file thirdparty.cpp.

Referenced by bessel_y1().

const double BIG = 1.44115188075855872E+17 [static]

Definition at line 2092 of file thirdparty.cpp.

Referenced by expn(), and ffun1().

const double C1 = 1.3862943611198906188e0 [static]

Definition at line 2013 of file thirdparty.cpp.

Referenced by ellpk(), Hydcs123(), and ILAENV().

const double DR1 = 5.78318596294678452118e0 [static]

Definition at line 656 of file thirdparty.cpp.

Referenced by bessel_j0().

const double DR2 = 3.04712623436620863991e1 [static]

Definition at line 658 of file thirdparty.cpp.

Referenced by bessel_j0().

const double elk_P[] [static]

Initial value:

{
        1.37982864606273237150e-4,
        2.28025724005875567385e-3,
        7.97404013220415179367e-3,
        9.85821379021226008714e-3,
        6.87489687449949877925e-3,
        6.18901033637687613229e-3,
        8.79078273952743772254e-3,
        1.49380448916805252718e-2,
        3.08851465246711995998e-2,
        9.65735902811690126535e-2,
        1.38629436111989062502e0
}

Definition at line 1983 of file thirdparty.cpp.

Referenced by ellpk().

const double elk_Q[] [static]

Initial value:

{
        2.94078955048598507511e-5,
        9.14184723865917226571e-4,
        5.94058303753167793257e-3,
        1.54850516649762399335e-2,
        2.39089602715924892727e-2,
        3.01204715227604046988e-2,
        3.73774314173823228969e-2,
        4.88280347570998239232e-2,
        7.03124996963957469739e-2,
        1.24999999999870820058e-1,
        4.99999999999999999821e-1
}

Definition at line 1998 of file thirdparty.cpp.

Referenced by ellpk().

const double i0_A[] [static]

Definition at line 1631 of file thirdparty.cpp.

Referenced by bessel_i0(), and bessel_i0_scaled().

const double i0_B[] [static]

Initial value:

{
        -7.23318048787475395456e-18,
        -4.83050448594418207126e-18,
        4.46562142029675999901e-17,
        3.46122286769746109310e-17,
        -2.82762398051658348494e-16,
        -3.42548561967721913462e-16,
        1.77256013305652638360e-15,
        3.81168066935262242075e-15,
        -9.55484669882830764870e-15,
        -4.15056934728722208663e-14,
        1.54008621752140982691e-14,
        3.85277838274214270114e-13,
        7.18012445138366623367e-13,
        -1.79417853150680611778e-12,
        -1.32158118404477131188e-11,
        -3.14991652796324136454e-11,
        1.18891471078464383424e-11,
        4.94060238822496958910e-10,
        3.39623202570838634515e-9,
        2.26666899049817806459e-8,
        2.04891858946906374183e-7,
        2.89137052083475648297e-6,
        6.88975834691682398426e-5,
        3.36911647825569408990e-3,
        8.04490411014108831608e-1
}

Definition at line 1671 of file thirdparty.cpp.

Referenced by bessel_i0(), and bessel_i0_scaled().

double i1_A[] [static]

Definition at line 1814 of file thirdparty.cpp.

Referenced by bessel_i1(), and bessel_i1_scaled().

double i1_B[] [static]

Initial value:

{
        7.51729631084210481353e-18,
        4.41434832307170791151e-18,
        -4.65030536848935832153e-17,
        -3.20952592199342395980e-17,
        2.96262899764595013876e-16,
        3.30820231092092828324e-16,
        -1.88035477551078244854e-15,
        -3.81440307243700780478e-15,
        1.04202769841288027642e-14,
        4.27244001671195135429e-14,
        -2.10154184277266431302e-14,
        -4.08355111109219731823e-13,
        -7.19855177624590851209e-13,
        2.03562854414708950722e-12,
        1.41258074366137813316e-11,
        3.25260358301548823856e-11,
        -1.89749581235054123450e-11,
        -5.58974346219658380687e-10,
        -3.83538038596423702205e-9,
        -2.63146884688951950684e-8,
        -2.51223623787020892529e-7,
        -3.88256480887769039346e-6,
        -1.10588938762623716291e-4,
        -9.76109749136146840777e-3,
        7.78576235018280120474e-1
}

Definition at line 1853 of file thirdparty.cpp.

Referenced by bessel_i1(), and bessel_i1_scaled().

int initf = 0 [static]

Definition at line 2461 of file thirdparty.cpp.

Referenced by init_by_array(), init_genrand(), and next_state().

const double k0_A[] [static]

Initial value:

{
        1.37446543561352307156e-16,
        4.25981614279661018399e-14,
        1.03496952576338420167e-11,
        1.90451637722020886025e-9,
        2.53479107902614945675e-7,
        2.28621210311945178607e-5,
        1.26461541144692592338e-3,
        3.59799365153615016266e-2,
        3.44289899924628486886e-1,
        -5.35327393233902768720e-1
}

Definition at line 1285 of file thirdparty.cpp.

Referenced by bessel_k0(), and bessel_k0_scaled().

const double k0_B[] [static]

Initial value:

 {
        5.30043377268626276149e-18,
        -1.64758043015242134646e-17,
        5.21039150503902756861e-17,
        -1.67823109680541210385e-16,
        5.51205597852431940784e-16,
        -1.84859337734377901440e-15,
        6.34007647740507060557e-15,
        -2.22751332699166985548e-14,
        8.03289077536357521100e-14,
        -2.98009692317273043925e-13,
        1.14034058820847496303e-12,
        -4.51459788337394416547e-12,
        1.85594911495471785253e-11,
        -7.95748924447710747776e-11,
        3.57739728140030116597e-10,
        -1.69753450938905987466e-9,
        8.57403401741422608519e-9,
        -4.66048989768794782956e-8,
        2.76681363944501510342e-7,
        -1.83175552271911948767e-6,
        1.39498137188764993662e-5,
        -1.28495495816278026384e-4,
        1.56988388573005337491e-3,
        -3.14481013119645005427e-2,
        2.44030308206595545468e0
}

Definition at line 1305 of file thirdparty.cpp.

Referenced by bessel_k0(), and bessel_k0_scaled().

const double k1_A[] [static]

Initial value:

{
        -7.02386347938628759343e-18,
        -2.42744985051936593393e-15,
        -6.66690169419932900609e-13,
        -1.41148839263352776110e-10,
        -2.21338763073472585583e-8,
        -2.43340614156596823496e-6,
        -1.73028895751305206302e-4,
        -6.97572385963986435018e-3,
        -1.22611180822657148235e-1,
        -3.53155960776544875667e-1,
        1.52530022733894777053e0
}

Definition at line 1459 of file thirdparty.cpp.

Referenced by bessel_k1(), and bessel_k1_scaled().

const double k1_B[] [static]

Initial value:

{
        -5.75674448366501715755e-18,
        1.79405087314755922667e-17,
        -5.68946255844285935196e-17,
        1.83809354436663880070e-16,
        -6.05704724837331885336e-16,
        2.03870316562433424052e-15,
        -7.01983709041831346144e-15,
        2.47715442448130437068e-14,
        -8.97670518232499435011e-14,
        3.34841966607842919884e-13,
        -1.28917396095102890680e-12,
        5.13963967348173025100e-12,
        -2.12996783842756842877e-11,
        9.21831518760500529508e-11,
        -4.19035475934189648750e-10,
        2.01504975519703286596e-9,
        -1.03457624656780970260e-8,
        5.74108412545004946722e-8,
        -3.50196060308781257119e-7,
        2.40648494783721712015e-6,
        -1.93619797416608296024e-5,
        1.95215518471351631108e-4,
        -2.85781685962277938680e-3,
        1.03923736576817238437e-1,
        2.72062619048444266945e0
}

Definition at line 1480 of file thirdparty.cpp.

Referenced by bessel_k1(), and bessel_k1_scaled().

const double MAXLOG = log(DBL_MAX) [static]

Definition at line 2091 of file thirdparty.cpp.

Referenced by expn().

unsigned long* next [static]

Definition at line 2462 of file thirdparty.cpp.

Referenced by genrand_int31(), genrand_int32(), genrand_real1(), genrand_real2(), genrand_real3(), and next_state().

int nleft = 1 [static]

Definition at line 2460 of file thirdparty.cpp.

Referenced by genrand_int31(), genrand_int32(), genrand_real1(), genrand_real2(), genrand_real3(), init_by_array(), init_genrand(), next_state(), and partx().

const double PIO4 = PI/4. [static]

Definition at line 681 of file thirdparty.cpp.

Referenced by bessel_j0(), and bessel_y0().

const double pre_factorial[NPRE_FACTORIAL] [static]

Definition at line 165 of file thirdparty.cpp.

Referenced by factorial().

const double SQ2OPI = sqrt(2./PI) [static]

Definition at line 680 of file thirdparty.cpp.

Referenced by bessel_j0(), bessel_j1(), bessel_y0(), and bessel_y1().

unsigned long state[N] [static]

Definition at line 2459 of file thirdparty.cpp.

const double THPIO4 = 3.*PI/4. [static]

Definition at line 912 of file thirdparty.cpp.

Referenced by bessel_j1(), and bessel_y1().

const double TWOOPI = 2./PI [static]

Definition at line 679 of file thirdparty.cpp.

Referenced by bessel_y0(), and bessel_y1().

const double Z1 = 1.46819706421238932572E1 [static]

Definition at line 909 of file thirdparty.cpp.

Referenced by bessel_j1().

const double Z2 = 4.92184563216946036703E1 [static]

Definition at line 910 of file thirdparty.cpp.

Referenced by bessel_j1().


Generated on Mon Feb 16 12:13:43 2009 for cloudy by  doxygen 1.4.7