diff options
| author | Rich Felker <dalias@aerifal.cx> | 2012-04-30 03:26:53 -0400 |
|---|---|---|
| committer | Rich Felker <dalias@aerifal.cx> | 2012-04-30 03:26:53 -0400 |
| commit | f6819755779a084bf2f82cb90175a4d9a018de73 (patch) | |
| tree | 94048c26dd47a1f82507421df1d8815e14675197 /include/math.h | |
| parent | 63374ee233c44891db80e6600d7a5a8c82e4ccca (diff) | |
| download | musl-f6819755779a084bf2f82cb90175a4d9a018de73.zip musl-f6819755779a084bf2f82cb90175a4d9a018de73.tar.gz musl-f6819755779a084bf2f82cb90175a4d9a018de73.tar.bz2 | |
first try at writing an efficient and "correct" exp10
this is a nonstandard function so it's not clear what conditions it
should satisfy. my intent is that it be fast and exact for positive
integral exponents when the result fits in the destination type, and
fast and correctly rounded for small negative integral exponents.
otherwise we aim for at most 1ulp error; it seems to differ from pow
by at most 1ulp and it's often 2-5 times faster than pow.
Diffstat (limited to 'include/math.h')
| -rw-r--r-- | include/math.h | 4 |
1 files changed, 4 insertions, 0 deletions
diff --git a/include/math.h b/include/math.h index bc0159e..65b1e91 100644 --- a/include/math.h +++ b/include/math.h @@ -387,6 +387,10 @@ float y1f(float); long double y1l(long double); float ynf(int, float); long double ynl(int, long double); + +double exp10(double); +float exp10f(float); +long double exp10l(long double); #endif #ifdef __cplusplus |