diff options
| author | Vineet Gupta <vgupta@synopsys.com> | 2020-05-29 15:21:39 -0700 |
|---|---|---|
| committer | Vineet Gupta <vgupta@synopsys.com> | 2020-06-03 10:23:22 -0700 |
| commit | 3374868668e708581ca06e256f6122518c89a8ad (patch) | |
| tree | 597e6a659674c190f2f388fa9f5e58b5c7cc0dfb /sysdeps/ieee754/flt-32 | |
| parent | ba9f6ee9bb8a894c9e2fb715edf693dd157b420a (diff) | |
| download | glibc-3374868668e708581ca06e256f6122518c89a8ad.zip glibc-3374868668e708581ca06e256f6122518c89a8ad.tar.gz glibc-3374868668e708581ca06e256f6122518c89a8ad.tar.bz2 | |
ieee754: provide gcc builtins based generic sqrt functions
Reviewed-by: Adhemerval Zanella <adhemerval.zanella@linaro.org>
Diffstat (limited to 'sysdeps/ieee754/flt-32')
| -rw-r--r-- | sysdeps/ieee754/flt-32/e_sqrtf.c | 16 |
1 files changed, 10 insertions, 6 deletions
diff --git a/sysdeps/ieee754/flt-32/e_sqrtf.c b/sysdeps/ieee754/flt-32/e_sqrtf.c index b339444..177bf5d 100644 --- a/sysdeps/ieee754/flt-32/e_sqrtf.c +++ b/sysdeps/ieee754/flt-32/e_sqrtf.c @@ -16,12 +16,15 @@ #include <math.h> #include <math_private.h> #include <libm-alias-finite.h> - -static const float one = 1.0, tiny=1.0e-30; +#include <math-use-builtins.h> float __ieee754_sqrtf(float x) { +#if USE_SQRTF_BUILTIN + return __builtin_sqrtf (x); +#else + /* Use generic implementation. */ float z; int32_t sign = (int)0x80000000; int32_t ix,s,q,m,t,i; @@ -70,10 +73,10 @@ __ieee754_sqrtf(float x) /* use floating add to find out rounding direction */ if(ix!=0) { - z = one-tiny; /* trigger inexact flag */ - if (z>=one) { - z = one+tiny; - if (z>one) + z = 0x1p0 - 0x1.4484cp-100; /* trigger inexact flag. */ + if (z >= 0x1p0) { + z = 0x1p0 + 0x1.4484cp-100; + if (z > 0x1p0) q += 2; else q += (q&1); @@ -83,6 +86,7 @@ __ieee754_sqrtf(float x) ix += (m <<23); SET_FLOAT_WORD(z,ix); return z; +#endif /* ! USE_SQRTF_BUILTIN */ } #ifndef __ieee754_sqrtf libm_alias_finite (__ieee754_sqrtf, __sqrtf) |