diff options
author | das <das@FreeBSD.org> | 2005-04-05 23:27:47 +0000 |
---|---|---|
committer | das <das@FreeBSD.org> | 2005-04-05 23:27:47 +0000 |
commit | bb5d9cb76836662658cc4da11687e6c90002189c (patch) | |
tree | 8ba01e3a2ce69dab2d7306b4c6c44d9f79e44225 /lib/msun | |
parent | 0cf2f689219da669cfcdf2a1664e8a16eaf01e90 (diff) | |
download | FreeBSD-src-bb5d9cb76836662658cc4da11687e6c90002189c.zip FreeBSD-src-bb5d9cb76836662658cc4da11687e6c90002189c.tar.gz |
Fix a long-standing bug in k_rem_pio2(), which led to large errors when
tanf() was called with big arguments close to multiples of pi/2.
Reported by: ucbtest via bde
Diffstat (limited to 'lib/msun')
-rw-r--r-- | lib/msun/src/k_rem_pio2f.c | 2 |
1 files changed, 1 insertions, 1 deletions
diff --git a/lib/msun/src/k_rem_pio2f.c b/lib/msun/src/k_rem_pio2f.c index 655d2a8..e1e9e35 100644 --- a/lib/msun/src/k_rem_pio2f.c +++ b/lib/msun/src/k_rem_pio2f.c @@ -88,7 +88,7 @@ recompute: iq[jz-1] -= i<<(8-q0); ih = iq[jz-1]>>(7-q0); } - else if(q0==0) ih = iq[jz-1]>>8; + else if(q0==0) ih = iq[jz-1]>>7; else if(z>=(float)0.5) ih=2; if(ih>0) { /* q > 0.5 */ |