Question regarding libm

[Manuel Novoa III]

On Fri, Sep 30, 2005 at 10:01:42PM +0200, Jan-Benedict Glaw wrote:
> On Thu, 2005-09-29 22:02:06 +0200, Jan-Benedict Glaw <jbglaw at lug-owl.de> wrote:
> > On Thu, 2005-09-29 13:39:21 -0600, Manuel Novoa III <mjn3 at codepoet.org> wrote:
> 
> > > Currently, _fpmaxtostr() assumes a base 2 representation.  But there's
> > > a compile-time test for that using FLT_RADIX.  Otherwise, it makes no
> > > assumptions about internal representation of float.  It _does_ assume
> > > some IEEE float behavior though.  If you could isolate the problem down
> > > to a particular area of the code, I'd be happy to take a look.
> 
> Actually, some assumptions about the float representation are done in
> some places. I just figured out that floor() does quite some fiddling
> with the bits. From my view, the xxx_WORD() macros (in math_private.h)
> are quite IEEE centric. Maybe I'll add another version using the VAX
> representation.
> 
> I guess that's why all my sine values are so off, too...

We didn't write the libm code.  I started a rewrite a while back but
have been too busy to work on it for some time.

Manuel