/home66/gary/public_html/cloudy/c10_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 F1(x, y, z)   (z ^ (x & (y ^ z)))
#define F2(x, y, z)   F1(z, x, y)
#define F3(x, y, z)   (x ^ y ^ z)
#define F4(x, y, z)   (y ^ (x | ~z))
#define MD5STEP(f, w, x, y, z, data, s)   w = rotlFixed(w + f(x, y, z) + data, s) + x

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, const double xorg[], const double yorg[], 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)
unsigned long MIXBITS (unsigned long u, unsigned long v)
unsigned long TWIST (unsigned long u, unsigned long v)
void init_genrand (unsigned long s)
void init_by_array (unsigned long init_key[], int key_length)
static void next_state ()
unsigned long genrand_int32 ()
long genrand_int31 ()
double genrand_real1 ()
double genrand_real2 ()
double genrand_real3 ()
double genrand_res53 ()
void humlik (int n, realnum x[], realnum y, realnum k[])
STATIC uint32 MD5swap (uint32 word)
STATIC void MD5_Transform (uint32 *digest, const uint32 *in)
string MD5file (const char *fnam, access_scheme scheme)
string MD5datafile (const char *fnam, access_scheme scheme)
string MD5string (const string &str)
uint32 rotlFixed (uint32 x, unsigned int y)

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 const int N = 624
static const int M = 397
static const unsigned long MATRIX_A = 0x9908b0dfUL
static const unsigned long UMASK = 0x80000000UL
static const unsigned long LMASK = 0x7fffffffUL
static unsigned long state [N]
static int nleft = 1
static int initf = 0
static unsigned long * nexxt

Define Documentation

#define F1 ( x,
y,
z   )     (z ^ (x & (y ^ z)))

Referenced by MD5_Transform(), and phymir_state< X, Y, NP, NSTR >::p_setup_next_hyperblock().

#define F2 ( x,
y,
z   )     F1(z, x, y)

Referenced by MD5_Transform(), and phymir_state< X, Y, NP, NSTR >::p_setup_next_hyperblock().

#define F3 ( x,
y,
z   )     (x ^ y ^ z)

Referenced by MD5_Transform().

#define F4 ( x,
y,
z   )     (y ^ (x | ~z))

Referenced by MD5_Transform().

#define MD5STEP ( f,
w,
x,
y,
z,
data,
s   )     w = rotlFixed(w + f(x, y, z) + data, s) + x

Referenced by MD5_Transform().

Function Documentation

double bessel_i0 ( double  x  ) 

Definition at line 1725 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 1742 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 1907 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 1928 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 707 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 938 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 1041 of file thirdparty.cpp.

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

Here is the call graph for this function:

double bessel_k0 ( double  x  ) 

Definition at line 1358 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 1381 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 1534 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 1556 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 745 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 965 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 1176 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 431 of file thirdparty.cpp.

References DEBUG_ENTRY.

Referenced by DoSutherland(), and Hypergeometric2F1().

double chbevl ( double  x,
const double  array[],
int  n 
) [inline]

Definition at line 2402 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 2040 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 2120 of file thirdparty.cpp.

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

Referenced by ParseDrive(), and 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 355 of file thirdparty.cpp.

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

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

long genrand_int31 ( void   ) 

Definition at line 2586 of file thirdparty.cpp.

References next_state(), nexxt, and nleft.

Here is the call graph for this function:

unsigned long genrand_int32 ( void   ) 

Definition at line 2569 of file thirdparty.cpp.

References next_state(), nexxt, and nleft.

Referenced by genrand_res53().

Here is the call graph for this function:

double genrand_real1 ( void   ) 

Definition at line 2603 of file thirdparty.cpp.

References next_state(), nexxt, and nleft.

Here is the call graph for this function:

double genrand_real2 ( void   ) 

Definition at line 2621 of file thirdparty.cpp.

References next_state(), nexxt, and nleft.

Here is the call graph for this function:

double genrand_real3 ( void   ) 

Definition at line 2639 of file thirdparty.cpp.

References next_state(), nexxt, and nleft.

Referenced by RandGauss().

Here is the call graph for this function:

double genrand_res53 ( void   ) 

Definition at line 2657 of file thirdparty.cpp.

References genrand_int32().

Here is the call graph for this function:

void humlik ( int  n,
realnum  x[],
realnum  y,
realnum  k[] 
)

Definition at line 2679 of file thirdparty.cpp.

References a0, e2(), and h2.

Referenced by getVoigt().

Here is the call graph for this function:

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

Definition at line 2521 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 2502 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 398 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,
const double  xorg[],
const double  yorg[],
double &  a,
double &  siga,
double &  b,
double &  sigb 
)

Definition at line 45 of file thirdparty.cpp.

References ASSERT, DEBUG_ENTRY, and pow2().

Referenced by iter_track::deriv(), and iter_track::zero_fit().

Here is the call graph for this function:

STATIC void MD5_Transform ( uint32 *  digest,
const uint32 *  in 
)

Definition at line 3060 of file thirdparty.cpp.

References DEBUG_ENTRY, F1, F2, F3, F4, and MD5STEP.

Referenced by MD5string().

string MD5datafile ( const char *  fnam,
access_scheme  scheme = AS_DATA_ONLY 
)

non-standard MD5 algorithm that skips eol characters and comments lines

Definition at line 2880 of file thirdparty.cpp.

References DEBUG_ENTRY, MD5string(), mode_r, and open_data().

Referenced by t_continuum::t_continuum().

Here is the call graph for this function:

string MD5file ( const char *  fnam,
access_scheme  scheme = AS_DATA_ONLY 
)

calculate the MD5 sum of a file

Definition at line 2861 of file thirdparty.cpp.

References DEBUG_ENTRY, MD5string(), mode_r, and open_data().

Here is the call graph for this function:

string MD5string ( const string &  str  ) 

calculate the MD5 sum of a string

Definition at line 2897 of file thirdparty.cpp.

References ASSERT, cpu, DEBUG_ENTRY, t_cpu::little_endian(), MD5_Transform(), and state.

Referenced by MD5datafile(), and MD5file().

Here is the call graph for this function:

STATIC uint32 MD5swap ( uint32  word  ) 

Definition at line 2964 of file thirdparty.cpp.

References DEBUG_ENTRY.

unsigned long MIXBITS ( unsigned long  u,
unsigned long  v 
) [inline]

Definition at line 2487 of file thirdparty.cpp.

References LMASK, and UMASK.

Referenced by TWIST().

static void next_state (  )  [static]

Definition at line 2547 of file thirdparty.cpp.

References init_genrand(), initf, M, N, nexxt, 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 2328 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 2306 of file thirdparty.cpp.

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

uint32 rotlFixed ( uint32  x,
unsigned int  y 
) [inline]

Definition at line 3055 of file thirdparty.cpp.

unsigned long TWIST ( unsigned long  u,
unsigned long  v 
) [inline]

Definition at line 2491 of file thirdparty.cpp.

References MATRIX_A, and MIXBITS().

Referenced by next_state().

Here is the call graph for this function:

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 615 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 625 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 635 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 646 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 684 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 691 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 657 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 668 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 870 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 880 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 890 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 901 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 851 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 858 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 912 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 921 of file thirdparty.cpp.

Referenced by bessel_y1().

const double BIG = 1.44115188075855872E+17 [static]

Definition at line 2117 of file thirdparty.cpp.

Referenced by expn(), and ffun1().

const double C1 = 1.3862943611198906188e0 [static]

Definition at line 2038 of file thirdparty.cpp.

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

const double DR1 = 5.78318596294678452118e0 [static]

Definition at line 680 of file thirdparty.cpp.

Referenced by bessel_j0().

const double DR2 = 3.04712623436620863991e1 [static]

Definition at line 682 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 2008 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 2023 of file thirdparty.cpp.

Referenced by ellpk().

const double i0_A[] [static]

Definition at line 1656 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 1696 of file thirdparty.cpp.

Referenced by bessel_i0(), and bessel_i0_scaled().

double i1_A[] [static]

Definition at line 1839 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 1878 of file thirdparty.cpp.

Referenced by bessel_i1(), and bessel_i1_scaled().

int initf = 0 [static]

Definition at line 2498 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 1310 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 1330 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 1484 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 1505 of file thirdparty.cpp.

Referenced by bessel_k1(), and bessel_k1_scaled().

const unsigned long LMASK = 0x7fffffffUL [static]

Definition at line 2486 of file thirdparty.cpp.

Referenced by MIXBITS().

const int M = 397 [static]

Definition at line 2483 of file thirdparty.cpp.

Referenced by iso_state_lifetime(), next_state(), and pah2_fun().

const unsigned long MATRIX_A = 0x9908b0dfUL [static]

Definition at line 2484 of file thirdparty.cpp.

Referenced by TWIST().

const double MAXLOG = log(DBL_MAX) [static]

Definition at line 2116 of file thirdparty.cpp.

Referenced by expn().

const int N = 624 [static]

Definition at line 2482 of file thirdparty.cpp.

Referenced by Badnell_rec_init(), Fe_10_11_13_cs(), HyperfineCreate(), init_by_array(), init_genrand(), next_state(), const_n_pointer< T, N, C_TYPE, true >::operator[](), const_n_pointer< T, N, ARPA_TYPE, true >::operator[](), const_n_pointer< T, N, C_TYPE, false >::operator[](), const_n_pointer< T, N, ARPA_TYPE, false >::operator[](), n_pointer< T, N, C_TYPE, true >::operator[](), n_pointer< T, N, ARPA_TYPE, true >::operator[](), n_pointer< T, N, C_TYPE, false >::operator[](), n_pointer< T, N, ARPA_TYPE, false >::operator[](), and read_SH98_He1_cross_sections().

unsigned long* nexxt [static]

Definition at line 2499 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 2497 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 705 of file thirdparty.cpp.

Referenced by bessel_j0(), and bessel_y0().

const double pre_factorial[NPRE_FACTORIAL] [static]

Definition at line 180 of file thirdparty.cpp.

Referenced by factorial().

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

Definition at line 704 of file thirdparty.cpp.

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

unsigned long state[N] [static]

Definition at line 2496 of file thirdparty.cpp.

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

Definition at line 936 of file thirdparty.cpp.

Referenced by bessel_j1(), and bessel_y1().

const double TWOOPI = 2./PI [static]

Definition at line 703 of file thirdparty.cpp.

Referenced by bessel_y0(), and bessel_y1().

const unsigned long UMASK = 0x80000000UL [static]

Definition at line 2485 of file thirdparty.cpp.

Referenced by MIXBITS().

const double Z1 = 1.46819706421238932572E1 [static]

Definition at line 933 of file thirdparty.cpp.

Referenced by bessel_j1().

const double Z2 = 4.92184563216946036703E1 [static]

Definition at line 934 of file thirdparty.cpp.

Referenced by bessel_j1().

 All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Defines
Generated on 15 Nov 2012 for cloudy by  doxygen 1.6.1