aboutsummaryrefslogtreecommitdiff
path: root/math
AgeCommit message (Collapse)Author
2012-11-19Fix spurious underflows in ldbl-128 atan implementation.David S. Miller
With help from Joseph Myers. * sysdeps/ieee754/ldbl-128/s_atanl.c (__atanl): Handle tiny and very large arguments properly. * math/libm-test.inc (atan_test): New tests. (atan2_test): New tests. * sysdeps/sparc/fpu/libm-test-ulps: Update. * sysdeps/x86_64/fpu/libm-test-ulps: Update.
2012-11-18Correct tinyness handling in long-double and float y0/y1.David S. Miller
With help from Joseph Myers. * sysdeps/ieee754/flt-32/e_j0f.c (__ieee754_y0f): Adjust tinyness cutoff to 2**-13. * sysdeps/ieee754/flt-32/e_j1f.c (__ieee754_y1f): Adjust tinyness cutoff to 2**-25. * sysdeps/ieee754/ldbl-128/e_j0l.c (U0): New constant. ( __ieee754_y0l): Avoid arithmetic underflow when 'x' is very small. * sysdeps/ieee754/ldbl-128/e_j1l.c (__ieee754_y1l): Likewise. * math/libm-test.inc (y0_test): New tests. (y1_test): New tests. * sysdeps/i386/fpu/libm-test-ulps: Update. * sysdeps/x86_64/fpu/libm-test-ulps: Update. * sysdeps/sparc/fpu/libm-test-ulps: Update.
2012-11-15Fix buffer overflow in libm-test.inc check_complex()Marcus Shawcroft
2012-11-07Fix spurious underflows from pow with results close to 1 (bug 14811).Joseph Myers
2012-11-06Fix fma underflows with small x * y (bug 14793).Joseph Myers
2012-11-04Fix fma overflow results outside round-to-nearest mode (bug 14797).Joseph Myers
2012-11-04[BZ #157] Remove include/stub-tag.h for good.Thomas Schwinge
2012-11-03Make fma use of Dekker and Knuth algorithms use round-to-nearest (bug 14796).Joseph Myers
2012-11-02Make math.h FP_* macros usable in #if (bug 3439).Joseph Myers
2012-11-01Fix fma (a, b, c) for small a * b (bugs 14784, 14785).Joseph Myers
2012-10-31Fix ldbl-128ibm atan2l for x near 1.Joseph Myers
2012-10-31Fix fma underflow exceptions in after-rounding edge cases.Joseph Myers
2012-10-30Fix fma missing underflows and bad results for some subnormal results (bugs ↵Joseph Myers
14152, 14783).
2012-10-18Comment on slow compilation of math/test-tgmath2.cMarkus Trippelsdorf
The math/test-tgmath2.c testcase uses very long lines (>3 million), so it may take a while to compile it. See: http://llvm.org/bugs/show_bug.cgi?id=14106 and http://gcc.gnu.org/bugzilla/show_bug.cgi?id=54402
2012-10-17Don't call feclearexcept in gamma function tests.Joseph Myers
2012-10-09Make libm-test.inc test comparison macros.Joseph Myers
2012-10-08Improve libm test coverage of classification macros / functions.Joseph Myers
2012-10-05Improve libm-test.inc comment.Joseph Myers
2012-10-01Fix sign of inexact zero return from fma (bug 14645).Joseph Myers
2012-09-29Fix sign of exact zero return from fma (bug 14638).Joseph Myers
2012-09-25Fix inaccuracy of clog, clog10 near |z| = 1 (bug 13629).Joseph Myers
2012-09-25Fix wrong ussage of sincos for subnormal argumentsLiubov Dmitrieva
2012-09-03Add optimized sinf and cosf routines for x86 and x86-64Liubov Dmitrieva
* sysdeps/i386/i686/fpu/multiarch/Makefile (sysdep_routines): Add s_sinf-sse2, s_conf-sse2. * sysdeps/i386/i686/fpu/multiarch/s_sinf-sse2.S: New file. * sysdeps/i386/i686/fpu/multiarch/s_cosf-sse2.S: New file. * sysdeps/i386/i686/fpu/multiarch/s_sinf.c: New file. * sysdeps/i386/i686/fpu/multiarch/s_cosf.c: New file. * sysdeps/ieee754/flt-32/s_sinf.c (SINF, SINF_FUNC): Add macros for using routine as __sinf_ia32. Use macro for function declaration and weak_alias. * sysdeps/ieee754/flt-32/s_cosf.c (COSF, COSF_FUNC): Add macros for using routine as __cosf_ia32. Use macro for function declaration and weak_alias. * sysdeps/i386/i686/fpu/multiarch/e_expf-sse2.S: Fix Copyright. * sysdeps/i386/i686/fpu/multiarch/e_expf.c: Fix Copyright. * sysdeps/x86_64/fpu/s_sinf.S: New file. * sysdeps/x86_64/fpu/s_cosf.S: New file. * sysdeps/x86_64/fpu/libm-test-ulps: Update. * math/libm-test.inc (cos_test): Add more test cases. (sin_test): Likewise. (sincos_test): Likewise.
2012-07-31Improve clog, clog10 handling of values with real or imaginary part slightly ↵Joseph Myers
above 1 (bug 13629).
2012-07-26Improve clog, clog10 handling of values with real or imaginary part 1 (bug ↵Joseph Myers
13629).
2012-07-25Set up errno properly for yn.Marek Polacek
2012-07-22Add missing includes.Thomas Schwinge
Follow-up to commits 38de94a5efbc3de186d6f287d666f74e5b4c8247, 76da7265320010c7a273ed99f53938c0f32d5fad.
2012-07-12Fix ynl return value with LDBL_MIN.Marek Polacek
2012-07-09Fix clog, clog10 spurious underflow exceptions (bug 14337).Joseph Myers
2012-07-06Fix tanf underflow close to pi/4 (bug 14154).Joseph Myers
2012-07-06Fix expm1 spurious underflow exceptions (bug 6778).Joseph Myers
2012-07-05Fix csqrt underflow (bugs 14157, 14331).Joseph Myers
2012-07-04Fix ctan, ctanh of subnormals in round-upwards mode (bug 14328).Joseph Myers
2012-07-03Fix float range reduction problems (bug 14283).Joseph Myers
2012-07-03Add cosh tests for bug 14273.Joseph Myers
2012-06-15Fix pow (0.0, -Inf) (bug 14241).Joseph Myers
2012-06-14Allow more spurious underflow exceptions from csqrt and Bessel functions.Joseph Myers
2012-06-01Fix fmod for subnormals (bug 14048).Joseph Myers
2012-05-31Use more decimal places for long double constants in math.h.Joseph Myers
2012-05-26Fix fma test for underflow.Marek Polacek
2012-05-25Fix acosf underflow (bug 14153).Joseph Myers
2012-05-24Add tests for underflow exception bugs 10846, 14036.Joseph Myers
2012-05-24Don't include exceptions in libm-test-ulps test names.Joseph Myers
2012-05-24Test underflow exceptions in libm-test.inc.Joseph Myers
2012-05-19Fix ccos, csin, ccosh, csinh overflows (bug 14123).Joseph Myers
2012-05-17math: support platforms with limited FP rounding or exception supportChris Metcalf
For some tests, just claim that fetestexcept() always returns true, so the rest of the test can be compiled. For libm-test, provide known bogus values for unsupported rounding modes, so fesetround() will return failure. Elsewhere, just add some #ifdefs to avoid code that uses particular FP exceptions if the exceptions aren't supported.
2012-05-10Add a missing #include.Roland McGrath
2012-05-10Fix for logb/logbf/logbl (bugs 13954/13955/13956)Adhemerval Zanella
POSIX 2008 states that if the input for 'logb[f|l]' is a subnormal number it should be treated as if it were normalized. This means the implementation should calculate the log2 of the mantissa and add it to the subnormal exponent (-126 for float and -1022 for double and IBM long double). This patch takes care of that.
2012-05-07Fix x86/x86_64 expm1l inaccuracy and exceptions (bugs 13885, 13923).Joseph Myers
2012-05-06Fix exp10 inaccuracy and exceptions (bugs 13884, 13914).Joseph Myers