cloudy trunk
Loading...
Searching...
No Matches
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 F1(x, y, z)
#define F2(x, y, z)
#define F3(x, y, z)
#define F4(x, y, z)
#define MD5STEP(f, w, x, y, z, data, s)

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)
double erfce (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 []
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]

Macro Definition Documentation

◆ F1

#define F1 ( x,
y,
z )

◆ F2

#define F2 ( x,
y,
z )

◆ F3

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

Referenced by t_lfact::Singleton< t_lfact >.

◆ F4

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

Referenced by t_lfact::Singleton< t_lfact >.

◆ MD5STEP

#define MD5STEP ( f,
w,
x,
y,
z,
data,
s )
Value:
w = rotlFixed(w + f(x, y, z) + data, s) + x
uint32 rotlFixed(uint32 x, unsigned int y)

Referenced by t_lfact::Singleton< t_lfact >.

Function Documentation

◆ bessel_i0()

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:

◆ bessel_i0_scaled()

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:

◆ bessel_i1()

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:

◆ bessel_i1_scaled()

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:

◆ bessel_j0()

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:

◆ bessel_j1()

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:

◆ bessel_jn()

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:

◆ bessel_k0()

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:

◆ bessel_k0_scaled()

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:

◆ bessel_k1()

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:

◆ bessel_k1_scaled()

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:

◆ bessel_y0()

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:

◆ bessel_y1()

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:

◆ bessel_yn()

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:

◆ cdgamma()

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

Definition at line 432 of file thirdparty.cpp.

References DEBUG_ENTRY.

Referenced by DoSutherland(), and Hypergeometric2F1().

◆ chbevl()

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

◆ dawson()

double dawson ( double x,
int order )
inline

Definition at line 3137 of file thirdparty.cpp.

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

Referenced by FastVoigtH().

Here is the call graph for this function:

◆ ellpk()

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:

◆ erfce()

double erfce ( double x)

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:

◆ expn()

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 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:

◆ factorial()

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().

◆ FastVoigtH()

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:

◆ genrand_int31()

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:

◆ genrand_int32()

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:

◆ genrand_real1()

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:

◆ genrand_real2()

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:

◆ genrand_real3()

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:

◆ genrand_res53()

double genrand_res53 ( void )

Definition at line 2989 of file thirdparty.cpp.

References genrand_int32().

Here is the call graph for this function:

◆ humlik()

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:

◆ init_by_array()

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:

◆ init_genrand()

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().

◆ lfactorial()

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:

◆ linfit()

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:

◆ MD5_Transform()

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

Definition at line 3624 of file thirdparty.cpp.

Referenced by MD5string().

◆ MD5datafile()

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:

◆ MD5file()

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:

◆ MD5string()

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, MD5_Transform(), MD5swap(), and state.

Referenced by MD5datafile(), and MD5file().

Here is the call graph for this function:

◆ MD5swap()

STATIC uint32 MD5swap ( uint32 word)

Definition at line 3528 of file thirdparty.cpp.

References DEBUG_ENTRY, and STATIC.

Referenced by MD5string().

◆ MIXBITS()

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().

◆ next_state()

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:

◆ p1evl()

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

Definition at line 2660 of file thirdparty.cpp.

References N.

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

◆ polevl()

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

Definition at line 2638 of file thirdparty.cpp.

References N.

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

◆ rotlFixed()

uint32 rotlFixed ( uint32 x,
unsigned int y )
inline

Definition at line 3619 of file thirdparty.cpp.

◆ TWIST()

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

◆ b0_PP

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().

◆ b0_PQ

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().

◆ b0_QP

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().

◆ b0_QQ

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().

◆ b0_RP

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().

◆ b0_RQ

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().

◆ b0_YP

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().

◆ b0_YQ

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().

◆ b1_PP

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().

◆ b1_PQ

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().

◆ b1_QP

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().

◆ b1_QQ

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().

◆ b1_RP

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().

◆ b1_RQ

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().

◆ b1_YP

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().

◆ b1_YQ

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().

◆ BIG

const double BIG = 1.44115188075855872E+17
static

Definition at line 2118 of file thirdparty.cpp.

Referenced by expn(), and ffun1().

◆ C1

const double C1 = 1.3862943611198906188e0
static

Definition at line 2039 of file thirdparty.cpp.

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

◆ DR1

const double DR1 = 5.78318596294678452118e0
static

Definition at line 681 of file thirdparty.cpp.

Referenced by bessel_j0().

◆ DR2

const double DR2 = 3.04712623436620863991e1
static

Definition at line 683 of file thirdparty.cpp.

Referenced by bessel_j0().

◆ elk_P

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().

◆ elk_Q

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().

◆ erf_P

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 erfce().

◆ erf_Q

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 erfce().

◆ erf_R

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 erfce().

◆ erf_S

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 erfce().

◆ i0_A

const double i0_A[]
static

Definition at line 1657 of file thirdparty.cpp.

Referenced by bessel_i0(), and bessel_i0_scaled().

◆ i0_B

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().

◆ i1_A

double i1_A[]
static

Definition at line 1840 of file thirdparty.cpp.

Referenced by bessel_i1(), and bessel_i1_scaled().

◆ i1_B

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().

◆ initf

int initf = 0
static

Definition at line 2830 of file thirdparty.cpp.

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

◆ k0_A

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().

◆ k0_B

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().

◆ k1_A

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().

◆ k1_B

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().

◆ LMASK

const unsigned long LMASK = 0x7fffffffUL
static

Definition at line 2818 of file thirdparty.cpp.

Referenced by MIXBITS().

◆ M

const int M = 397
static

◆ MATRIX_A

const unsigned long MATRIX_A = 0x9908b0dfUL
static

Definition at line 2816 of file thirdparty.cpp.

Referenced by TWIST().

◆ MAXLOG

const double MAXLOG = log(DBL_MAX)
static

Definition at line 2117 of file thirdparty.cpp.

Referenced by expn().

◆ N

◆ N_DAWSON

const int N_DAWSON = 100

Definition at line 3002 of file thirdparty.cpp.

Referenced by dawson().

◆ nexxt

unsigned long* nexxt
static

◆ nleft

◆ PIO4

const double PIO4 = PI/4.
static

Definition at line 706 of file thirdparty.cpp.

Referenced by bessel_j0(), and bessel_y0().

◆ pre_factorial

const double pre_factorial[NPRE_FACTORIAL]
static

Definition at line 181 of file thirdparty.cpp.

Referenced by factorial().

◆ SQ2OPI

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().

◆ state

unsigned long state[N]
static

Definition at line 2828 of file thirdparty.cpp.

◆ tbl_dawson

const double tbl_dawson[N_DAWSON+1]
static

Definition at line 3006 of file thirdparty.cpp.

Referenced by dawson().

◆ THPIO4

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

Definition at line 937 of file thirdparty.cpp.

Referenced by bessel_j1(), and bessel_y1().

◆ TWOOPI

const double TWOOPI = 2./PI
static

Definition at line 704 of file thirdparty.cpp.

Referenced by bessel_y0(), and bessel_y1().

◆ UMASK

const unsigned long UMASK = 0x80000000UL
static

Definition at line 2817 of file thirdparty.cpp.

Referenced by MIXBITS().

◆ Z1

const double Z1 = 1.46819706421238932572E1
static

Definition at line 934 of file thirdparty.cpp.

Referenced by bessel_j1().

◆ Z2

const double Z2 = 4.92184563216946036703E1
static

Definition at line 935 of file thirdparty.cpp.

Referenced by bessel_j1().