d/freebsd-skq
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Fixed the unexpectedly large maximum error after the previous commit.

Browse Source 
It was because I forgot to translate the part of the double precision
algorithm that chops t so that t*t is exact.  Now the maximum error
is the same as for double precision (almost exactly 2.0/3 ulps).
This commit is contained in:
bde 2005-12-11 17:58:14 +00:00
parent 3c8b4e6094
commit 94455e43e1
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
lib/msun/src/s_cbrtf.c
View File

@ -67,9 +67,9 @@ cbrtf(float x)
/* chop t to 12 bits and make it larger than cbrt(x) */
GET_FLOAT_WORD(high,t);
SET_FLOAT_WORD(t,high+0x00001000);
SET_FLOAT_WORD(t,(high&0xfffff000)+0x00001000);
/* one step Newton iteration to 24 bits with error less than 0.984 ulps */
/* one step Newton iteration to 24 bits with error less than 0.667 ulps */
s=t*t; /* t*t is exact */
r=x/s;
w=t+t;