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/fmod.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/fmod.h')
| -rw-r--r-- | newlib/libm/machine/spu/headers/fmod.h | 14 |
1 files changed, 14 insertions, 0 deletions
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_ */ |