cloudy  trunk
 All Data Structures Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
Data Structures | Macros | Functions | Variables
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
 

Macros

#define N   624
 
#define M   397
 
#define MATRIX_A   0x9908b0dfUL /* constant vector a */
 
#define UMASK   0x80000000UL /* most significant w-r bits */
 
#define LMASK   0x7fffffffUL /* least significant r bits */
 
#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
 

Macro Definition Documentation

#define LMASK   0x7fffffffUL /* least significant r bits */

Definition at line 2455 of file thirdparty.cpp.

#define M   397

Definition at line 2452 of file thirdparty.cpp.

Referenced by DGEMM(), DGER(), DTRSM(), iso_state_lifetime(), and next_state().

#define MATRIX_A   0x9908b0dfUL /* constant vector a */

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.

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

#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 /* most significant w-r bits */

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  x,
const double  array[],
int  n 
)
inline
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
nis order, 1 for first integral integral
xis 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().

Referenced by SanityCheckBegin().

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, t_lfact::get_lfact(), Singleton< t_lfact >::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.

References N.

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
int nleft = 1
static
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 for cloudy by doxygen 1.8.3.1