diff options
Diffstat (limited to 'sysdeps/ieee754/dbl-64/e_sqrt.c')
-rw-r--r-- | sysdeps/ieee754/dbl-64/e_sqrt.c | 5 |
1 files changed, 2 insertions, 3 deletions
diff --git a/sysdeps/ieee754/dbl-64/e_sqrt.c b/sysdeps/ieee754/dbl-64/e_sqrt.c index fff6d148fe..c1fed7d97e 100644 --- a/sysdeps/ieee754/dbl-64/e_sqrt.c +++ b/sysdeps/ieee754/dbl-64/e_sqrt.c @@ -21,7 +21,7 @@ /* */ /* FUNCTION: usqrt */ /* */ -/* FILES NEEDED: dla.h endian.h mydefs.h uroot.h */ +/* FILES NEEDED: dla.h endian.h mydefs.h */ /* uroot.tbl */ /* */ /* An ultimate sqrt routine. Given an IEEE double machine number x */ @@ -47,7 +47,6 @@ double __ieee754_sqrt (double x) { -#include "uroot.h" static const double rt0 = 9.99999999859990725855365213134618E-01, rt1 = 4.99999999495955425917856814202739E-01, @@ -134,7 +133,7 @@ __ieee754_sqrt (double x) return x; /* sqrt(+0)=+0, sqrt(-0)=-0 */ if (k < 0) return (x - x) / (x - x); /* sqrt(-ve)=sNaN */ - return tm256.x * __ieee754_sqrt (x * t512.x); + return 0x1p-256 * __ieee754_sqrt (x * 0x1p512); } } strong_alias (__ieee754_sqrt, __sqrt_finite) |