aboutsummaryrefslogtreecommitdiff
path: root/math
AgeCommit message (Collapse)AuthorFilesLines
2013-04-11libm-test.inc: Fix tests where cos(PI/2) != 0.Carlos O'Donell1-2/+50
The value of PI is never exactly PI in any floating point representation, and the value of PI/2 is never PI/2. It is wrong to expect cos(M_PI_2l) to return 0, instead it will return an answer that is non-zero because M_PI_2l doesn't round to exactly PI/2 in the type used. That is to say that the correct answer is to do the following: * Take PI or PI/2. * Round to the floating point representation. * Take the rounded value and compute an infinite precision cos or sin. * Use the rounded result of the infinite precision cos or sin as the answer to the test. I used printf to do the type rounding, and Wolfram's Alpha to do the infinite precision cos calculations. The following changes bring x86-64 and x86 to 1/2 ulp for two tests. It shows that the x86 cos implementation is quite good, and that our test are flawed. Unfortunately given that the rounding errors are type dependent we need to fix this for each type. No regressions on x86-64 or x86. --- 2013-04-11 Carlos O'Donell <carlos@redhat.com> * math/libm-test.inc (cos_test): Fix PI/2 test. (sincos_test): Likewise. * sysdeps/x86_64/fpu/libm-test-ulps: Regenerate. * sysdeps/i386/fpu/libm-test-ulps: Regenerate.
2013-04-11Don't use run-via-rtld-prefix for anything other than testsSiddhesh Poyarekar1-1/+1
run-via-rtld-prefix checks whether the program to be run is a static test and skips if it is. This is fine, except that it assumes that the program to be run is the second $^, which is true only for tests. This change creates an rtld-prefix, which is simply the dynamic linker prefix with the necessary arguments and uses that in the non-test targets.
2013-04-06README.libm-test: Use testrun.sh to run libm test.Carlos O'Donell1-2/+1
Document the use of the convenience testrun.sh script for running the libm test. --- 2013-04-06 Carlos O'Donell <carlos@redhat.com> * math/README.libm-test (How can I generate "libm-test-ulps"?): Use testrun.sh to run libm tests.
2013-04-06New Makefile target `regen-ulps'.Carlos O'Donell3-19/+78
The wiki "Regeneration" page has this to say about update ULPs. "The libm-test-ulps files are semiautomatically updated. To update an ulps baseline, run each of the failing tests (test-float, test-double, etc.) with -u; this will generate a file called ULPs; concatenate each of those files with the existing libm-test-ulps file, after removing any entries for particularly huge numbers of ulps that you do not want to mark as expected. Then run gen-libm-test.pl -n -u FILE where FILE is the concatenated file produced in the previous step. This generates a file called NewUlps which is the new sorted version of libm-test-ulps." The same information is listed in math/README.libm-test, and is a lot of manual work that you often want to run over-and-over again while working on a particular test. The `regen-ulps' convenience target does this automatically for developers. We strictly assume the source tree is readonly and add a new --output-dir option to libm-test.inc to allow for writing out ULPs to $(objpfx). When run the new target does the following: * Starts with the baseline ULPs file. * Runs each of the libm math tests with -u. * Adds new changes seen with -u to the baseline. * Sorts and prepares the test output with gen-libm-test.pl. * Leaves math/NewUlps in your build tree to copy to your source tree, cleanup, and checkin. The math test documentation in math/README.libm-test is updated document the new Makefile target. --- 2013-04-06 Carlos O'Donell <carlos@redhat.com> * Makefile.in (regen-ulps): New target. * math/Makefile [ifneq (no,$(PERL)]: Declare regen-ulps with .PHONY. [ifneq (no,$(PERL)] (run-regen-ulps): New variable. [ifneq (no,$(PERL)] (regen-ulps): New target. [ifeq (no,$(PERL)] (regen-ulps): New target. * math/libm-test.inc (ulps_file_name): Define. (output_dir): New variable. (options): Add "output-dir" option. (parse_opt): Handle 'o' case. (main): If output_dir is non-NULL use it as a prefix otherwise use "". * math/README.libm-test: Update `How can I generate "libm-test-ulps"?'
2013-04-05[BZ #14686, #15336] Fix standard compliance. Don't use hard-coded qNaN values.Thomas Schwinge1-0/+25
2013-04-03Fix missing underflow from cexp (bug 14478).Joseph Myers4-0/+39
2013-04-02Fix cacosh inaccuracy and spurious exceptions (bug 15327).Joseph Myers4-59/+639
2013-04-02New <math.h> macro named issignaling to check for a signaling NaN (sNaN).Thomas Schwinge8-15/+177
It is based on draft TS 18661 and currently enabled as a GNU extension.
2013-03-30Fix casinh inaccuracy for imaginary part < 1.0, real part small (bug 10357).Joseph Myers4-0/+1133
2013-03-27Fix casinh inaccuracy near i, imaginary part > 1 (bug 15307).Joseph Myers4-0/+491
2013-03-22BZ#13889: expl (709.75) wrongly overflows for ldbl-128ibmAdhemerval Zanella1-0/+3
The patch increase the high value to check if expl overflows. Current high mark value is not really correct, the algorithm accepts high values. It also adds a correct wrapper function to check for overflow and underflow.
2013-03-21libm-test: Properly wrap blocks consisting of several statements.Thomas Schwinge1-8/+14
2013-03-21On 32-bit x86, disable certain tests involving sNaN values.Thomas Schwinge1-10/+22
Follow-up to commit 495ded2c8c1eb8c0ac4b54add2dd397852e19cba.
2013-03-21Fix Bessel function spurious overflows for ldbl-128 / ldbl-128ibm (bug 15285).Joseph Myers1-0/+4
2013-03-21Fix casinh inaccuracy for argument with imaginary part 1 (bug 15287).Joseph Myers4-0/+332
2013-03-19Fix types of constants in k_casinh*.c.Joseph Myers2-4/+4
2013-03-16Fix y1l spurious overflows for ldbl-96 (bug 15283).Joseph Myers1-0/+20
2013-03-16Fix spurious underflow exceptions for Bessel functions for ldbl-128 / ↵Joseph Myers1-8/+4
ldbl-128ibm (bug 14155).
2013-03-16Move "-sNaN" value into a separate variable.Thomas Schwinge1-6/+9
2013-03-15 * math/libm-test.inc (ldexp_test): Add missing START/END markers.Thomas Schwinge1-0/+4
2013-03-15Promote a math test for sNaN handling to the top-level.Thomas Schwinge2-1/+337
2013-03-15Use GCC's builtins for generating NaNs.Thomas Schwinge2-18/+19
2013-03-15Add one test, removed a duplicated one, add a few comments about missing tests.Thomas Schwinge1-1/+15
2013-03-15Better distinguish between NaN/qNaN/sNaN.Thomas Schwinge5-600/+600
2013-03-15Enable tests for x86_64, too, which currently are enabled only for x86.Thomas Schwinge2-2/+2
2013-03-15Fix copy'n'pastos.Thomas Schwinge1-2/+2
2013-03-14Fix spurious underflow exceptions for Bessel functions for double (bug 14155).Joseph Myers1-4/+2
2013-02-19Fix some testsuite build warning fixes in libmSiddhesh Poyarekar3-19/+22
2013-02-15Remove bounded-pointers build system support.Joseph Myers1-3/+0
2013-02-14Remove bp-sym.h and BP_SYM uses from C code.Joseph Myers1-3/+2
2013-01-31Fix casinh spurious underflows away from [-i,i] (bug 15062).Joseph Myers4-0/+147
2013-01-17Fix cacos real-part inaccuracy for result real part near 0 (bug 15023).Joseph Myers11-128/+366
2013-01-11Add 64-bit VIS3 optimized GMP routines for sparc.David S. Miller1-1/+2
* math/Makefile: Recognize gmp-sysdep_routines. * sysdeps/sparc/sparc64/multiarch/Makefile: Add VIS3 optimized GMP routines to sysdeps. * sysdeps/sparc/sparc64/multiarch/add_n-vis3.S: New file. * sysdeps/sparc/sparc64/multiarch/add_n.S: New file. * sysdeps/sparc/sparc64/multiarch/addmul_1-vis3.S: New file. * sysdeps/sparc/sparc64/multiarch/addmul_1.S: New file. * sysdeps/sparc/sparc64/multiarch/mul_1-vis3.S: New file. * sysdeps/sparc/sparc64/multiarch/mul_1.S: New file. * sysdeps/sparc/sparc64/multiarch/sub_n-vis3.S: New file. * sysdeps/sparc/sparc64/multiarch/sub_n.S: New file. * sysdeps/sparc/sparc64/multiarch/submul_1-vis3.S: New file. * sysdeps/sparc/sparc64/multiarch/submul_1.S: New file.
2013-01-10Use __extension__ with long long in installed headers.Joseph Myers1-0/+2
2013-01-10Revert "CFLAGS-test-tgmath2.c"Andreas Schwab1-1/+1
This reverts commit 3167bfc62bb38e3057f7694566b5bb157236a5a8.
2013-01-10CFLAGS-test-tgmath2.cAndreas Schwab1-1/+1
2013-01-07Fix casinh, casin overflow (bug 14996).Joseph Myers4-20/+85
2013-01-04Fix casinh, casin inaccuracy from cancellation (bug 14994).Joseph Myers4-18/+88
2013-01-02Update copyright notices with scripts/update-copyrights.Joseph Myers192-195/+192
2012-12-30math: use existing nonnull attribute defineMike Frysinger1-2/+1
Signed-off-by: Mike Frysinger <vapier@gentoo.org>
2012-11-28Fix powl inaccuracy for x86_64 and x86 (bug 13881).Joseph Myers1-4/+31
2012-11-22Disable some fma tests for ldbl-128ibm.Joseph Myers1-0/+6
2012-11-22Fix ldbl-128ibm atanl spurious underflows (bug 14871).Joseph Myers1-0/+8
2012-11-22Fix ldbl-128ibm hypotl internal underflows (bug 14869).Joseph Myers1-0/+3
2012-11-22Fix ldbl-128ibm hypotl inaccuracy for arguments with large ratio (bug 14868).Joseph Myers1-0/+2
2012-11-22Mark some clog / clog10 underflows as expected for IBM long double.Joseph Myers1-4/+4
2012-11-19Fix spurious underflows in ldbl-128 atan implementation.David S. Miller1-0/+5
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. Miller1-0/+24
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 Shawcroft1-3/+9
2012-11-07Fix spurious underflows from pow with results close to 1 (bug 14811).Joseph Myers1-0/+29
generated by cgit v1.1 at 2024-11-28 22:28:02 +0000