Date: Thu, 9 May 2013 21:14:25 +0200 From: Szabolcs Nagy <nsz@...t70.net> To: musl@...ts.openwall.com Subject: Re: Using float_t and double_t in math functions * Rich Felker <dalias@...ifal.cx> [2013-05-09 10:57:14 -0400]: > > Do you want to do it, or do you want me to? I don't mind but you're > more familiar with the code and probably better aware of where it's > okay to change. (BTW, it's probably not safe to change arg-reduction > code, right?) > i will do it, but it will take some time there are some simple cases like sinf where it should clearly work (polynomial evals, most temporaries) but in other cases the excess precision can actually hurt (correctly rounded operations, double-double representation, over/underflow flags should be raised,..) and there are some non-trivial cases eg. acosh(x) for large x is ln(2*x), or to avoid overflow: log(x) + 0.693147180559945309417232121458176568; which could be rewritten as log(x) + (double_t)0.693147180559945309417232121458176568L; so gcc recognizes the const and uses fldln2 on i386 (which probably saves 1-2 bytes) (the double rounding is not an issue for special consts) so i think this will need to be done with extensive testing and inspection of the generated code
Powered by blists - more mailing lists
Confused about mailing lists and their use? Read about mailing lists on Wikipedia and check out these guidelines on proper formatting of your messages.