summaryrefslogtreecommitdiffstats
path: root/newlib/libm/machine/spu/headers/sqrt.h
diff options
context:
space:
mode:
authorJeff Johnston <jjohnstn@redhat.com>2007-09-28 18:58:10 +0000
committerJeff Johnston <jjohnstn@redhat.com>2007-09-28 18:58:10 +0000
commit1671fbe1ca0825a77cb9137a08cf4ad8a48d0f95 (patch)
tree79e8a71cd7363e82effb317acda5cd10c04050c7 /newlib/libm/machine/spu/headers/sqrt.h
parente30a7b84aa1527df63f20fe0cb3605afab652d94 (diff)
downloadcygnal-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.h5
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_ */