/home66/gary/public_html/cloudy/c13_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)
double dawson (double x, int order)
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)
STATIC double expx2 (double x, int sign)
double erfc (double a)
double erfce (double x)
double erf (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 ()
realnum FastVoigtH (realnum a, realnum v)
void humlik (int n, const 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 double erf_P []
static double erf_Q []
static double erf_R []
static double erf_S []
const bool lgUSE_EXPXSQ = true
static double erf_T []
static double erf_U []
const double exp_M = 128.0
const double exp_MINV = .0078125
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
const int N_DAWSON = 100
static const double tbl_dawson [N_DAWSON+1]

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 1726 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 1743 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 1908 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 1929 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 708 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 939 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 1042 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 1359 of file thirdparty.cpp.

References bessel_i0(), cdEXIT, chbevl(), DEBUG_ENTRY, EXIT_FAILURE, 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 1382 of file thirdparty.cpp.

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

Here is the call graph for this function:

double bessel_k1 ( double  x  ) 

Definition at line 1535 of file thirdparty.cpp.

References bessel_i1(), cdEXIT, chbevl(), DEBUG_ENTRY, EXIT_FAILURE, 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 1557 of file thirdparty.cpp.

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

Here is the call graph for this function:

double bessel_y0 ( double  x  ) 

Definition at line 746 of file thirdparty.cpp.

References b0_PP, b0_PQ, b0_QP, b0_QQ, b0_YP, b0_YQ, bessel_j0(), cdEXIT, DEBUG_ENTRY, EXIT_FAILURE, 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 966 of file thirdparty.cpp.

References b1_PP, b1_PQ, b1_QP, b1_QQ, b1_YP, b1_YQ, bessel_j1(), cdEXIT, DEBUG_ENTRY, EXIT_FAILURE, 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 1177 of file thirdparty.cpp.

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

Here is the call graph for this function:

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

Definition at line 432 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 2734 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 dawson ( double  x,
int  order 
) [inline]

Definition at line 3137 of file thirdparty.cpp.

References max(), min(), N_DAWSON, tbl_dawson, and TotalInsanity().

Referenced by FastVoigtH().

Here is the call graph for this function:

double ellpk ( double  x  ) 

Definition at line 2041 of file thirdparty.cpp.

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

Referenced by StarkCollTransProb_VF01().

Here is the call graph for this function:

double erf ( double   ) 

the standard error functions

Definition at line 2492 of file thirdparty.cpp.

References DEBUG_ENTRY, erf_T, erf_U, erfc(), p1evl(), and polevl().

Referenced by erfc().

Here is the call graph for this function:

double erfc ( double  a  ) 

Definition at line 2404 of file thirdparty.cpp.

References DEBUG_ENTRY, erf(), erf_P, erf_Q, erf_R, erf_S, expx2(), lgUSE_EXPXSQ, MAXLOG, p1evl(), and polevl().

Referenced by erf().

Here is the call graph for this function:

double erfce ( double   ) 

erfce(a) returns erfc(a)*exp(a^2)

Definition at line 2454 of file thirdparty.cpp.

References DEBUG_ENTRY, erf_P, erf_Q, erf_R, erf_S, p1evl(), and polevl().

Referenced by VoigtH0(), and VoigtU0().

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 2121 of file thirdparty.cpp.

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

Referenced by k2DampArg::operator()(), ParseDrive(), and SanityCheckBegin().

Here is the call graph for this function:

STATIC double expx2 ( double  x,
int  sign 
)

Definition at line 2551 of file thirdparty.cpp.

References DEBUG_ENTRY, exp_M, exp_MINV, and MAXLOG.

Referenced by erfc().

double factorial ( long  n  ) 

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

Definition at line 356 of file thirdparty.cpp.

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

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

realnum FastVoigtH ( realnum  a,
realnum  v 
)

Definition at line 3171 of file thirdparty.cpp.

References ASSERT, dawson(), DEBUG_ENTRY, dsexp(), order(), pow2(), and SQRTPI.

Referenced by VoigtH().

Here is the call graph for this function:

long genrand_int31 ( void   ) 

Definition at line 2918 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 2901 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 2935 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 2953 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 2971 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 2989 of file thirdparty.cpp.

References genrand_int32().

Here is the call graph for this function:

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

Definition at line 3257 of file thirdparty.cpp.

References a0, DEBUG_ENTRY, e2(), and h2.

Referenced by VoigtH().

Here is the call graph for this function:

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

Definition at line 2853 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 2834 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 399 of file thirdparty.cpp.

References cdEXIT, DEBUG_ENTRY, EXIT_FAILURE, 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 46 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 3624 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 3444 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 3425 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 3461 of file thirdparty.cpp.

References ASSERT, cpu, DEBUG_ENTRY, t_cpu::i(), t_cpu_i::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 3528 of file thirdparty.cpp.

References DEBUG_ENTRY.

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

Definition at line 2819 of file thirdparty.cpp.

References LMASK, and UMASK.

Referenced by TWIST().

static void next_state (  )  [static]

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

Referenced by bessel_j0(), bessel_j1(), bessel_y0(), bessel_y1(), erf(), erfc(), and erfce().

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

Definition at line 2638 of file thirdparty.cpp.

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

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

Definition at line 3619 of file thirdparty.cpp.

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

Definition at line 2823 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 616 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 626 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 636 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 647 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 685 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 692 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 658 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 669 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 871 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 881 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 891 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 902 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 852 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 859 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 913 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 922 of file thirdparty.cpp.

Referenced by bessel_y1().

const double BIG = 1.44115188075855872E+17 [static]

Definition at line 2118 of file thirdparty.cpp.

Referenced by expn(), and ffun1().

const double C1 = 1.3862943611198906188e0 [static]

Definition at line 2039 of file thirdparty.cpp.

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

const double DR1 = 5.78318596294678452118e0 [static]

Definition at line 681 of file thirdparty.cpp.

Referenced by bessel_j0().

const double DR2 = 3.04712623436620863991e1 [static]

Definition at line 683 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 2009 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 2024 of file thirdparty.cpp.

Referenced by ellpk().

double erf_P[] [static]
Initial value:
 {
 2.46196981473530512524e-10,
 5.64189564831068821977e-1,
 7.46321056442269912687e0,
 4.86371970985681366614e1,
 1.96520832956077098242e2,
 5.26445194995477358631e2,
 9.34528527171957607540e2,
 1.02755188689515710272e3,
 5.57535335369399327526e2
}

Definition at line 2354 of file thirdparty.cpp.

Referenced by erfc(), and erfce().

double erf_Q[] [static]
Initial value:
 {

 1.32281951154744992508e1,
 8.67072140885989742329e1,
 3.54937778887819891062e2,
 9.75708501743205489753e2,
 1.82390916687909736289e3,
 2.24633760818710981792e3,
 1.65666309194161350182e3,
 5.57535340817727675546e2
}

Definition at line 2365 of file thirdparty.cpp.

Referenced by erfc(), and erfce().

double erf_R[] [static]
Initial value:
 {
 5.64189583547755073984e-1,
 1.27536670759978104416e0,
 5.01905042251180477414e0,
 6.16021097993053585195e0,
 7.40974269950448939160e0,
 2.97886665372100240670e0
}

Definition at line 2376 of file thirdparty.cpp.

Referenced by erfc(), and erfce().

double erf_S[] [static]
Initial value:
 {

 2.26052863220117276590e0,
 9.39603524938001434673e0,
 1.20489539808096656605e1,
 1.70814450747565897222e1,
 9.60896809063285878198e0,
 3.36907645100081516050e0
}

Definition at line 2384 of file thirdparty.cpp.

Referenced by erfc(), and erfce().

double erf_T[] [static]
Initial value:
 {
 9.60497373987051638749e0,
 9.00260197203842689217e1,
 2.23200534594684319226e3,
 7.00332514112805075473e3,
 5.55923013010394962768e4
}

Definition at line 2476 of file thirdparty.cpp.

Referenced by erf().

double erf_U[] [static]
Initial value:
 {

 3.35617141647503099647e1,
 5.21357949780152679795e2,
 4.59432382970980127987e3,
 2.26290000613890934246e4,
 4.92673942608635921086e4
}

Definition at line 2483 of file thirdparty.cpp.

Referenced by erf().

const double exp_M = 128.0

Definition at line 2548 of file thirdparty.cpp.

Referenced by expx2().

const double exp_MINV = .0078125

Definition at line 2549 of file thirdparty.cpp.

Referenced by expx2().

const double i0_A[] [static]

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

Referenced by bessel_i0(), and bessel_i0_scaled().

double i1_A[] [static]

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

Referenced by bessel_i1(), and bessel_i1_scaled().

int initf = 0 [static]

Definition at line 2830 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 1311 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 1331 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 1485 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 1506 of file thirdparty.cpp.

Referenced by bessel_k1(), and bessel_k1_scaled().

const bool lgUSE_EXPXSQ = true

Definition at line 2402 of file thirdparty.cpp.

Referenced by erfc().

const unsigned long LMASK = 0x7fffffffUL [static]

Definition at line 2818 of file thirdparty.cpp.

Referenced by MIXBITS().

const int M = 397 [static]

Definition at line 2815 of file thirdparty.cpp.

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

const unsigned long MATRIX_A = 0x9908b0dfUL [static]

Definition at line 2816 of file thirdparty.cpp.

Referenced by TWIST().

const double MAXLOG = log(DBL_MAX) [static]

Definition at line 2117 of file thirdparty.cpp.

Referenced by erfc(), expn(), and expx2().

const int N = 624 [static]

Definition at line 2814 of file thirdparty.cpp.

Referenced by Badnell_rec_init(), HyperfineCreate(), init_by_array(), init_genrand(), next_state(), and read_SH98_He1_cross_sections().

const int N_DAWSON = 100

Definition at line 3002 of file thirdparty.cpp.

Referenced by dawson().

unsigned long* nexxt [static]

Definition at line 2831 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 2829 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 706 of file thirdparty.cpp.

Referenced by bessel_j0(), and bessel_y0().

const double pre_factorial[NPRE_FACTORIAL] [static]

Definition at line 181 of file thirdparty.cpp.

Referenced by factorial().

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

Definition at line 705 of file thirdparty.cpp.

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

unsigned long state[N] [static]

Definition at line 2828 of file thirdparty.cpp.

const double tbl_dawson[N_DAWSON+1] [static]

Definition at line 3006 of file thirdparty.cpp.

Referenced by dawson().

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

Definition at line 937 of file thirdparty.cpp.

Referenced by bessel_j1(), and bessel_y1().

const double TWOOPI = 2./PI [static]

Definition at line 704 of file thirdparty.cpp.

Referenced by bessel_y0(), and bessel_y1().

const unsigned long UMASK = 0x80000000UL [static]

Definition at line 2817 of file thirdparty.cpp.

Referenced by MIXBITS().

const double Z1 = 1.46819706421238932572E1 [static]

Definition at line 934 of file thirdparty.cpp.

Referenced by bessel_j1().

const double Z2 = 4.92184563216946036703E1 [static]

Definition at line 935 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