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/fmod.h | 14 ++++++++++++++
 1 file changed, 14 insertions(+)

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

diff --git a/newlib/libm/machine/spu/headers/fmod.h b/newlib/libm/machine/spu/headers/fmod.h
index 1506fb8dc..0236da187 100644
--- a/newlib/libm/machine/spu/headers/fmod.h
+++ b/newlib/libm/machine/spu/headers/fmod.h
@@ -34,6 +34,7 @@
 #define _FMOD_H_	1
 
 #include <spu_intrinsics.h>
+#include <errno.h>
 #include "headers/vec_literal.h"
 
 /* This implementation returns zero if y is a denorm or zero.
@@ -57,6 +58,10 @@ static __inline double _fmod(double x, double y)
   vec_uint4 sign_mask = (vec_uint4)(VEC_SPLAT_U64(0x8000000000000000ULL));
   vec_uint4 implied_1 = (vec_uint4)(VEC_SPLAT_U64(0x0010000000000000ULL));
   vec_uint4 mant_mask = (vec_uint4)(VEC_SPLAT_U64(0x000FFFFFFFFFFFFFULL));
+  vec_ullong2 domain;
+  vec_int4 verrno;
+  vec_double2 vc = { 0.0, 0.0 };
+  vec_int4 fail = { EDOM, EDOM, EDOM, EDOM };
 
   vx = (vec_uint4)spu_promote(x, 0);
   vy = (vec_uint4)spu_promote(y, 0);
@@ -147,6 +152,15 @@ static __inline double _fmod(double x, double y)
   result = spu_sel(spu_andc(result, spu_rlmask(result0, -1)), vx,
                    resultx);
 
+#ifndef _IEEE_LIBM
+  /*
+   * If y is zero, set errno to EDOM
+   */
+  domain = spu_cmpeq(vc, (vec_double2) vy);
+  verrno = spu_splats(errno);
+  errno = spu_extract(spu_sel(verrno, fail, (vector unsigned int) domain), 0);
+#endif
+
   return (spu_extract((vec_double2)result, 0));
 }
 #endif /* _FMOD_H_ */
-- 
cgit v1.2.3