c = 0.003379167095512573896158903121545171688;
80 result= z * 1.125 + z * c;
90 static const doubleY = 1.044948577880859375;
91 static const double P[] = {0.0834305892146531832907,
92-0.338165134459360935041,
93-0.0509990735146777432841,
94-0.00772758345802133288487,
95-0.000322780120964605683831};
97 static const doubleQ[] = {1.0,
980.455004033050794024546,
990.0875222600142252549554,
1000.00858571925074406212772,
1010.000370900071787748000569};
105 doublep = (((
P[4] * zz +
P[3]) * zz +
P[2]) * zz +
P[1]) * zz +
P[0];
106 doubleq = (((Q[4] * zz + Q[3]) * zz + Q[2]) * zz + Q[1]) * zz + Q[0];
111 else if(
invert? (z < 28.0) : (z < 5.8)) {
122 static const doubleY = 0.405935764312744140625;
123 static const double P[] = {-0.098090592216281240205,
1240.178114665841120341155,
1250.191003695796775433986,
1260.0888900368967884466578,
1270.0195049001251218801359,
1280.00180424538297014223957};
130 static const doubleQ[] = {1.0,
1311.84759070983002217845,
1321.42628004845511324508,
1330.578052804889902404909,
1340.12385097467900864233,
1350.0113385233577001411017,
1360.337511472483094676155e-5};
139 doublep = ((((
P[5] * x +
P[4]) * x +
P[3]) * x +
P[2]) * x +
P[1]) * x +
P[0];
140 doubleq = ((((Q[5] * x + Q[4]) * x + Q[3]) * x + Q[2]) * x + Q[1]) * x + Q[0];
143 result*= expl(-z * z) / z;
151 static const doubleY = 0.50672817230224609375;
152 static const double P[] = {-0.0243500476207698441272,
1530.0386540375035707201728,
1540.04394818964209516296,
1550.0175679436311802092299,
1560.00323962406290842133584,
1570.000235839115596880717416};
159 static const doubleQ[] = {1.0,
1601.53991494948552447182,
1610.982403709157920235114,
1620.325732924782444448493,
1630.0563921837420478160373,
1640.00410369723978904575884};
167 doublep = ((((
P[5] * x +
P[4]) * x +
P[3]) * x +
P[2]) * x +
P[1]) * x +
P[0];
168 doubleq = ((((Q[5] * x + Q[4]) * x + Q[3]) * x + Q[2]) * x + Q[1]) * x + Q[0];
171 result*= expl(-z * z) / z;
179 static const doubleY = 0.5405750274658203125;
180 static const double P[] = {0.00295276716530971662634,
1810.0137384425896355332126,
1820.00840807615555585383007,
1830.00212825620914618649141,
1840.000250269961544794627958,
1850.113212406648847561139e-4};
187 static const doubleQ[] = {1.0,
1881.04217814166938418171,
1890.442597659481563127003,
1900.0958492726301061423444,
1910.0105982906484876531489,
1920.000479411269521714493907};
195 doublep = ((((
P[5] * x +
P[4]) * x +
P[3]) * x +
P[2]) * x +
P[1]) * x +
P[0];
196 doubleq = ((((Q[5] * x + Q[4]) * x + Q[3]) * x + Q[2]) * x + Q[1]) * x + Q[0];
199 result*= expl(-z * z) / z;
207 static const doubleY = 0.5579090118408203125;
208 static const double P[] = {0.00628057170626964891937,
2090.0175389834052493308818,
210-0.212652252872804219852,
211-0.687717681153649930619,
212-2.5518551727311523996,
213-3.22729451764143718517,
214-2.8175401114513378771};
216 static const doubleQ[] = {1.0,
2172.79257750980575282228,
21811.0567237927800161565,
21915.930646027911794143,
22022.9367376522880577224,
22113.5064170191802889145,
2225.48409182238641741584};
225 doublep = (((((
P[6] * x +
P[5]) * x +
P[4]) * x +
P[3]) * x +
P[2]) * x +
P[1]) * x +
P[0];
226 doubleq = (((((Q[6] * x + Q[5]) * x + Q[4]) * x + Q[3]) * x + Q[2]) * x + Q[1]) * x + Q[0];
229 result*= expl(-z * z) / z;
double ErfC(double z)
Complementary error function.
double Erf(double z)
Error function.
static double ErfImpl(double z, int invert)
RetroSearch is an open source project built by @garambo | Open a GitHub Issue
Search and Browse the WWW like it's 1997 | Search results from DuckDuckGo
HTML:
3.2
| Encoding:
UTF-8
| Version:
0.7.4