diff options
Diffstat (limited to 'sysdeps/ieee754/dbl-64/e_j1.c')
| -rw-r--r-- | sysdeps/ieee754/dbl-64/e_j1.c | 8 |
1 files changed, 4 insertions, 4 deletions
diff --git a/sysdeps/ieee754/dbl-64/e_j1.c b/sysdeps/ieee754/dbl-64/e_j1.c index daf025fdb7..f5f5c28830 100644 --- a/sysdeps/ieee754/dbl-64/e_j1.c +++ b/sysdeps/ieee754/dbl-64/e_j1.c @@ -125,10 +125,10 @@ static double zero = 0.0; * j1(x) = 1/sqrt(pi) * (P(1,x)*cc - Q(1,x)*ss) / sqrt(x) * y1(x) = 1/sqrt(pi) * (P(1,x)*ss + Q(1,x)*cc) / sqrt(x) */ - if(ix>0x48000000) z = (invsqrtpi*cc)/__sqrt(y); + if(ix>0x48000000) z = (invsqrtpi*cc)/__ieee754_sqrt(y); else { u = pone(y); v = qone(y); - z = invsqrtpi*(u*cc-v*ss)/__sqrt(y); + z = invsqrtpi*(u*cc-v*ss)/__ieee754_sqrt(y); } if(hx<0) return -z; else return z; @@ -214,10 +214,10 @@ static double V0[5] = { * sin(x) +- cos(x) = -cos(2x)/(sin(x) -+ cos(x)) * to compute the worse one. */ - if(ix>0x48000000) z = (invsqrtpi*ss)/__sqrt(x); + if(ix>0x48000000) z = (invsqrtpi*ss)/__ieee754_sqrt(x); else { u = pone(x); v = qone(x); - z = invsqrtpi*(u*ss+v*cc)/__sqrt(x); + z = invsqrtpi*(u*ss+v*cc)/__ieee754_sqrt(x); } return z; } |