From 1671fbe1ca0825a77cb9137a08cf4ad8a48d0f95 Mon Sep 17 00:00:00 2001
From: Jeff Johnston <jjohnstn@redhat.com>
Date: Fri, 28 Sep 2007 18:58:10 +0000
Subject: 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.
---
 newlib/libm/machine/spu/headers/sqrt.h | 5 +++++
 1 file changed, 5 insertions(+)

(limited to 'newlib/libm/machine/spu/headers/sqrt.h')

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_ */
-- 
cgit v1.2.3