This is the mail archive of the mailing list for the glibc project.

Index Nav: [Date Index] [Subject Index] [Author Index] [Thread Index]
Message Nav: [Date Prev] [Date Next] [Thread Prev] [Thread Next]

Re: [] libc/2269: triginometric argument reductionerror in libm

On 20 May 2001, Ulrich Drepper wrote:

> Geoff Keating <> writes:
> > Someone needs to fix libc/sysdeps/i386/fpu/s_sin.S so that it does
> > proper range reduction, I think.  It works fine on non-x86 chips.
> It's not worth doing this.  What is the result of such a large
> argument good for?  It has absolutely no meaning.

Please go back and read my original posting.

(1) The sin and other trig functions have exact values for all normalized
IEEE floats.  The IEEE standard requires that the functions return the
closest IEEE value to this correct value.  The current library is
returning values that are grossly wrong.

(2) There are certainly scientific codes that depend upon getting this
right.  It's why (as you can see from what I sent out before) every major
vendor gets it right and has been getting it right for almost twenty

(3) If you suffer from the common misconception that IEEE  floats are
somehow "approximate" please go back and read the standard again.

Bruce Allen

Index Nav: [Date Index] [Subject Index] [Author Index] [Thread Index]
Message Nav: [Date Prev] [Date Next] [Thread Prev] [Thread Next]