fspecial.cpp File Reference

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <cmath>
#include <float.h>
#include "compat.h"
#include "constants.h"
#include "fspecial.h"

Go to the source code of this file.

Data Structures
struct	cheb_series_t

Typedefs
typedef struct cheb_series_t	cheb_series

Functions
static nr_double_t	cheb_eval (const cheb_series *cs, const nr_double_t x)
static nr_double_t	erfc8 (nr_double_t x)
static nr_double_t	erfseries (nr_double_t x)

Variables
const nr_double_t	SN_ACC = 1e-5
const nr_double_t	K_ERR = 1e-8
static nr_double_t	erfc_xlt1_data [20]
static cheb_series	erfc_xlt1_cs
static nr_double_t	erfc_x15_data [25]
static cheb_series	erfc_x15_cs
static nr_double_t	erfc_x510_data [20]
static cheb_series	erfc_x510_cs
static nr_double_t	bi0_data [12]
static cheb_series	bi0_cs
static nr_double_t	ai0_data [21]
static cheb_series	ai0_cs
static nr_double_t	ai02_data [22]
static cheb_series	ai02_cs

Typedef Documentation

typedef struct cheb_series_t cheb_series

Definition at line 191 of file fspecial.cpp.

Function Documentation

static nr_double_t cheb_eval ( const cheb_series *  cs, const nr_double_t  x  )

static

Definition at line 193 of file fspecial.cpp.

static nr_double_t erfc8 ( nr_double_t  x )

static

Definition at line 298 of file fspecial.cpp.

static nr_double_t erfseries ( nr_double_t  x )

static

Definition at line 362 of file fspecial.cpp.

Variable Documentation

cheb_series ai02_cs

static

Initial value:

= {
  ai02_data, 21, -1, 1
}

Definition at line 496 of file fspecial.cpp.

nr_double_t ai02_data[22]

static

Initial value:

= {
   .05449041101410882,
   .00336911647825569,
   .00006889758346918,
   .00000289137052082,
   .00000020489185893,
   .00000002266668991,
   .00000000339623203,
   .00000000049406022,
   .00000000001188914,
  -.00000000003149915,
  -.00000000001321580,
  -.00000000000179419,
   .00000000000071801,
   .00000000000038529,
   .00000000000001539,
  -.00000000000004151,
  -.00000000000000954,
   .00000000000000382,
   .00000000000000176,
  -.00000000000000034,
  -.00000000000000027,
   .00000000000000003
}

Definition at line 472 of file fspecial.cpp.

cheb_series ai0_cs

static

Initial value:

= {
  ai0_data, 20, -1, 1
}

Definition at line 468 of file fspecial.cpp.

nr_double_t ai0_data[21]

static

Initial value:

= {
   .07575994494023796,
   .00759138081082334,
   .00041531313389237,
   .00001070076463439,
  -.00000790117997921,
  -.00000078261435014,
   .00000027838499429,
   .00000000825247260,
  -.00000001204463945,
   .00000000155964859,
   .00000000022925563,
  -.00000000011916228,
   .00000000001757854,
   .00000000000112822,
  -.00000000000114684,
   .00000000000027155,
  -.00000000000002415,
  -.00000000000000608,
   .00000000000000314,
  -.00000000000000071,
   .00000000000000007
}

Definition at line 445 of file fspecial.cpp.

cheb_series bi0_cs

static

Initial value:

= {
  bi0_data, 11, -1, 1
}

Definition at line 441 of file fspecial.cpp.

nr_double_t bi0_data[12]

static

Initial value:

= {
  -.07660547252839144951,
  1.92733795399380827000,
   .22826445869203013390,
   .01304891466707290428,
   .00043442709008164874,
   .00000942265768600193,
   .00000014340062895106,
   .00000000161384906966,
   .00000000001396650044,
   .00000000000009579451,
   .00000000000000053339,
   .00000000000000000245
}

Definition at line 427 of file fspecial.cpp.

cheb_series erfc_x15_cs

static

Initial value:

= {
  erfc_x15_data, 24, -1, 1
}

Definition at line 264 of file fspecial.cpp.

nr_double_t erfc_x15_data[25]

static

Initial value:

= {
  0.44045832024338111077637466616,
 -0.143958836762168335790826895326,
  0.044786499817939267247056666937,
 -0.013343124200271211203618353102,
  0.003824682739750469767692372556,
 -0.001058699227195126547306482530,
  0.000283859419210073742736310108,
 -0.000073906170662206760483959432,
  0.000018725312521489179015872934,
 -4.62530981164919445131297264430e-6,
  1.11558657244432857487884006422e-6,
 -2.63098662650834130067808832725e-7,
  6.07462122724551777372119408710e-8,
 -1.37460865539865444777251011793e-8,
  3.05157051905475145520096717210e-9,
 -6.65174789720310713757307724790e-10,
  1.42483346273207784489792999706e-10,
 -3.00141127395323902092018744545e-11,
  6.22171792645348091472914001250e-12,
 -1.26994639225668496876152836555e-12,
  2.55385883033257575402681845385e-13,
 -5.06258237507038698392265499770e-14,
  9.89705409478327321641264227110e-15,
 -1.90685978789192181051961024995e-15,
  3.50826648032737849245113757340e-16
}

Definition at line 237 of file fspecial.cpp.

cheb_series erfc_x510_cs

static

Initial value:

= {
  erfc_x510_data, 19, -1, 1
}

Definition at line 292 of file fspecial.cpp.

nr_double_t erfc_x510_data[20]

static

Initial value:

= {
  1.11684990123545698684297865808,
  0.003736240359381998520654927536,
 -0.000916623948045470238763619870,
  0.000199094325044940833965078819,
 -0.000040276384918650072591781859,
  7.76515264697061049477127605790e-6,
 -1.44464794206689070402099225301e-6,
  2.61311930343463958393485241947e-7,
 -4.61833026634844152345304095560e-8,
  8.00253111512943601598732144340e-9,
 -1.36291114862793031395712122089e-9,
  2.28570483090160869607683087722e-10,
 -3.78022521563251805044056974560e-11,
  6.17253683874528285729910462130e-12,
 -9.96019290955316888445830597430e-13,
  1.58953143706980770269506726000e-13,
 -2.51045971047162509999527428316e-14,
  3.92607828989125810013581287560e-15,
 -6.07970619384160374392535453420e-16,
  9.12600607264794717315507477670e-17
}

Definition at line 270 of file fspecial.cpp.

cheb_series erfc_xlt1_cs

static

Initial value:

= {
  erfc_xlt1_data, 19, -1, 1
}

Definition at line 232 of file fspecial.cpp.

nr_double_t erfc_xlt1_data[20]

static

Initial value:

= {
  1.06073416421769980345174155056,
 -0.42582445804381043569204735291,
  0.04955262679620434040357683080,
  0.00449293488768382749558001242,
 -0.00129194104658496953494224761,
 -0.00001836389292149396270416979,
  0.00002211114704099526291538556,
 -5.23337485234257134673693179020e-7,
 -2.78184788833537885382530989578e-7,
  1.41158092748813114560316684249e-8,
  2.72571296330561699984539141865e-9,
 -2.06343904872070629406401492476e-10,
 -2.14273991996785367924201401812e-11,
  2.22990255539358204580285098119e-12,
  1.36250074650698280575807934155e-13,
 -1.95144010922293091898995913038e-14,
 -6.85627169231704599442806370690e-16,
  1.44506492869699938239521607493e-16,
  2.45935306460536488037576200030e-18,
 -9.29599561220523396007359328540e-19
}

Definition at line 210 of file fspecial.cpp.

const nr_double_t K_ERR = 1e-8

Definition at line 89 of file fspecial.cpp.

const nr_double_t SN_ACC = 1e-5

Definition at line 88 of file fspecial.cpp.