diff options
author | Jeff Johnston <jjohnstn@redhat.com> | 2007-09-28 18:58:10 +0000 |
---|---|---|
committer | Jeff Johnston <jjohnstn@redhat.com> | 2007-09-28 18:58:10 +0000 |
commit | 1671fbe1ca0825a77cb9137a08cf4ad8a48d0f95 (patch) | |
tree | 79e8a71cd7363e82effb317acda5cd10c04050c7 /newlib/libm/machine/spu/headers/sqrt.h | |
parent | e30a7b84aa1527df63f20fe0cb3605afab652d94 (diff) | |
download | cygnal-1671fbe1ca0825a77cb9137a08cf4ad8a48d0f95.tar.gz cygnal-1671fbe1ca0825a77cb9137a08cf4ad8a48d0f95.tar.bz2 cygnal-1671fbe1ca0825a77cb9137a08cf4ad8a48d0f95.zip |
2007-09-28 Patrick Mansfield <patmans@us.ibm.com>
* libm/machine/spu/headers/dom_chkd_less_than.h: New file.
* libm/machine/spu/headers/dom_chkd_negone_one.h: Ditto.
* libm/machine/spu/headers/dom_chkf_less_than.h: Ditto.
* libm/machine/spu/headers/dom_chkf_negone_one.h: Ditto.
* libm/machine/spu/headers/acos.h: Add domain check.
* libm/machine/spu/headers/acosf.h: Ditto.
* libm/machine/spu/headers/acosh.h: Ditto.
* libm/machine/spu/headers/acoshf.h: Ditto.
* libm/machine/spu/headers/asin.h: Ditto.
* libm/machine/spu/headers/asinf.h: Ditto.
* libm/machine/spu/headers/atanh.h: Ditto.
* libm/machine/spu/headers/atanhf.h: Ditto.
* libm/machine/spu/headers/fmod.h: Ditto.
* libm/machine/spu/headers/fmodf.h: Ditto.
* libm/machine/spu/headers/log.h: Ditto.
* libm/machine/spu/headers/log10.h: Ditto.
* libm/machine/spu/headers/log1p.h: Ditto.
* libm/machine/spu/headers/log1pf.h: Ditto.
* libm/machine/spu/headers/log2.h: Ditto.
* libm/machine/spu/headers/log2f.h: Ditto.
* libm/machine/spu/headers/sqrt.h: Ditto.
* libm/machine/spu/headers/sqrtf.h: Ditto.
* libm/machine/spu/headers/tgamma.h: Ditto.
* libm/machine/spu/headers/tgammaf.h: Ditto.
Diffstat (limited to 'newlib/libm/machine/spu/headers/sqrt.h')
-rw-r--r-- | newlib/libm/machine/spu/headers/sqrt.h | 5 |
1 files changed, 5 insertions, 0 deletions
diff --git a/newlib/libm/machine/spu/headers/sqrt.h b/newlib/libm/machine/spu/headers/sqrt.h index 064c47def..63b034dd9 100644 --- a/newlib/libm/machine/spu/headers/sqrt.h +++ b/newlib/libm/machine/spu/headers/sqrt.h @@ -43,6 +43,7 @@ */ #include <spu_intrinsics.h> #include "headers/vec_literal.h" +#include "headers/dom_chkd_less_than.h" static __inline double _sqrt(double in) { @@ -53,6 +54,7 @@ static __inline double _sqrt(double in) vec_ullong2 mask = VEC_SPLAT_U64(0x7FE0000000000000ULL); vec_double2 x, dx, de, dd, dy, dg, dy2, dhalf; vec_double2 denorm, neg; + vec_double2 vc = { 0.0, 0.0 }; fhalf = VEC_SPLAT_F32(0.5f); dhalf = VEC_SPLAT_F64(0.5); @@ -124,6 +126,9 @@ static __inline double _sqrt(double in) dg = spu_sel(spu_andc(spu_or(dg, neg), denorm), x, nochange); +#ifndef _IEEE_LIBM + dom_chkd_less_than(spu_splats(in), vc); +#endif return (spu_extract(dg, 0)); } #endif /* _SQRT_H_ */ |