From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: (qmail 90699 invoked by alias); 3 Mar 2020 15:45:49 -0000 Mailing-List: contact newlib-cvs-help@sourceware.org; run by ezmlm Precedence: bulk List-Id: List-Subscribe: List-Archive: List-Post: List-Help: , Sender: newlib-cvs-owner@sourceware.org Received: (qmail 90564 invoked by uid 9078); 3 Mar 2020 15:45:43 -0000 Date: Tue, 03 Mar 2020 15:45:00 -0000 Message-ID: <20200303154543.90563.qmail@sourceware.org> Content-Type: text/plain; charset="us-ascii" MIME-Version: 1.0 Content-Transfer-Encoding: 7bit From: Corinna Vinschen To: newlib-cvs@sourceware.org Subject: [newlib-cygwin] Fix error in float trig. function range reduction X-Act-Checkin: newlib-cygwin X-Git-Author: Fabian Schriever X-Git-Refname: refs/heads/master X-Git-Oldrev: ecf27dd2e0ed1dff4dc919a7c805e951913d953f X-Git-Newrev: d4bcecb3e91f73091195adbfdcb0d4287cdf1753 X-SW-Source: 2020-q1/txt/msg00026.txt https://sourceware.org/git/gitweb.cgi?p=newlib-cygwin.git;h=d4bcecb3e91f73091195adbfdcb0d4287cdf1753 commit d4bcecb3e91f73091195adbfdcb0d4287cdf1753 Author: Fabian Schriever Date: Tue Mar 3 14:49:06 2020 +0100 Fix error in float trig. function range reduction The single-precision trigonometric functions show rather high errors in specific ranges starting at about 30000 radians. For example the sinf procedure produces an error of 7626.55 ULP with the input 5.195880078125e+04 (0x474AF6CD) (compared with MPFR in 128bit precision). For the test we used 100k values evenly spaced in the range of [30k, 70k]. The issues are periodic at higher ranges. This error was introduced when the double precision range reduction was first converted to float. The shift by 8 bits always returns 0 as iq is never higher than 255. The fix reduces the error of the example above to 0.45 ULP, highest error within the test set fell to 1.31 ULP, which is not perfect, but still a significant improvement. Testing other previously erroneous ranges no longer show particularly large accuracy errors. Diff: --- newlib/libm/math/kf_rem_pio2.c | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/newlib/libm/math/kf_rem_pio2.c b/newlib/libm/math/kf_rem_pio2.c index 1573ca9..3ca7cc2 100644 --- a/newlib/libm/math/kf_rem_pio2.c +++ b/newlib/libm/math/kf_rem_pio2.c @@ -101,7 +101,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 */