static double erfc1(double x)
{
- double s,P,Q;
+ double_t s,P,Q;
s = fabs(x) - 1;
P = pa0+s*(pa1+s*(pa2+s*(pa3+s*(pa4+s*(pa5+s*pa6)))));
static double erfc2(uint32_t ix, double x)
{
- double s,z,R,S;
+ double_t s,R,S;
+ double z;
if (ix < 0x3ff40000) /* |x| < 1.25 */
return erfc1(x);