about summary refs log tree commit diff
path: root/sysdeps/ieee754
Commit message (Collapse)AuthorAgeFilesLines
* Simplify calculation of 2^-m in __mpexpSiddhesh Poyarekar2013-01-182-29/+34
|
* Remove unnecessary multiplication with RADIXISiddhesh Poyarekar2013-01-181-3/+3
|
* Fix header commentSiddhesh Poyarekar2013-01-171-2/+0
|
* Consolidate constant defines into mpa.hSiddhesh Poyarekar2013-01-1613-120/+26
|
* Fix the value of TWOSiddhesh Poyarekar2013-01-161-1/+1
|
* Use libc_fe*() in 32-bit nearbyint()David S. Miller2013-01-151-4/+4
| | | | | * sysdeps/ieee754/dbl-64/s_nearbyint.c (__nearbyint): Use libc_feholdexcept and libc_fesetenv.
* Minor tweak to mp multiplicationSiddhesh Poyarekar2013-01-141-13/+15
| | | | | Add a local variable to remove extra copies to/from memory in the Z array.
* Fix code formatting in mpa.cSiddhesh Poyarekar2013-01-141-228/+463
| | | | This includes the overridden mpa.c in power4.
* Revert "Use ieee754/dbl-64/wordsize-64 on powerpc64"Andreas Schwab2013-01-1017-117/+0
| | | | This reverts commit 7a9d2c397195d9a19a95dbcdb59fb5df03f2d940.
* Revert "sysdeps/ieee754/ldbl-128ibm/e_expl.c"Andreas Schwab2013-01-102-64/+16
| | | | This reverts commit 35eb879e3b7849e86f06658bdb888f7858b30759.
* Fix formatting in mpexp.cSiddhesh Poyarekar2013-01-101-54/+95
|
* Make __mpexp_twomm1 an array of doublesSiddhesh Poyarekar2013-01-102-87/+26
| | | | Cleanup to get rid of the `number` muck which is not necessary here.
* sysdeps/ieee754/ldbl-128ibm/e_expl.cAndreas Schwab2013-01-102-16/+64
|
* Use ieee754/dbl-64/wordsize-64 on powerpc64Andreas Schwab2013-01-1017-0/+117
| | | | | | | | | | | | | | | | | | | | | | | | * sysdeps/ieee754/ldbl-opt/wordsize-64/s_ceil.c: New file. * sysdeps/ieee754/ldbl-opt/wordsize-64/s_finite.c: New file. * sysdeps/ieee754/ldbl-opt/wordsize-64/s_floor.c: New file. * sysdeps/ieee754/ldbl-opt/wordsize-64/s_frexp.c: New file. * sysdeps/ieee754/ldbl-opt/wordsize-64/s_isinf.c: New file. * sysdeps/ieee754/ldbl-opt/wordsize-64/s_isnan.c: New file. * sysdeps/ieee754/ldbl-opt/wordsize-64/s_llround.c: New file. * sysdeps/ieee754/ldbl-opt/wordsize-64/s_logb.c: New file. * sysdeps/ieee754/ldbl-opt/wordsize-64/s_lround.c: New file. * sysdeps/ieee754/ldbl-opt/wordsize-64/s_modf.c: New file. * sysdeps/ieee754/ldbl-opt/wordsize-64/s_nearbyint.c: New file. * sysdeps/ieee754/ldbl-opt/wordsize-64/s_remquo.c: New file. * sysdeps/ieee754/ldbl-opt/wordsize-64/s_rint.c: New file. * sysdeps/ieee754/ldbl-opt/wordsize-64/s_round.c: New file. * sysdeps/ieee754/ldbl-opt/wordsize-64/s_scalbln.c: New file. * sysdeps/ieee754/ldbl-opt/wordsize-64/s_scalbn.c: New file. * sysdeps/ieee754/ldbl-opt/wordsize-64/s_trunc.c: New file. * sysdeps/unix/sysv/linux/powerpc/powerpc64/Implies: Add ieee754/ldbl-opt/wordsize-64. * sysdeps/powerpc/powerpc64/Implies: Add ieee754/dbl-64/wordsize-64.
* Update comments in mpa.cSiddhesh Poyarekar2013-01-091-67/+47
| | | | Fixed comment style and clearer wording in some cases.
* Remove unused __cr and __cpymnSiddhesh Poyarekar2013-01-041-33/+0
|
* Fix code formattingSiddhesh Poyarekar2013-01-041-10/+15
|
* Clean up comment for MP_NOSiddhesh Poyarekar2013-01-041-12/+26
|
* Remove some commented codeSiddhesh Poyarekar2013-01-041-13/+0
|
* Remove commented declarationsSiddhesh Poyarekar2013-01-041-6/+0
|
* Remove argument variable name from function declarationSiddhesh Poyarekar2013-01-041-1/+1
|
* Update copyright notices with scripts/update-copyrights.Joseph Myers2013-01-02196-201/+196
|
* Fix values in __mpexp_twomm1Siddhesh Poyarekar2013-01-021-8/+8
|
* Split mantissa calculation loop and add branch predictionSiddhesh Poyarekar2013-01-021-24/+43
|
* Add assert for potential access beyond array bounds in m1npSiddhesh Poyarekar2013-01-021-4/+17
| | | | | | | | | | | The mpexp code has an access into m1np: for (i=n-1; i>0; i--,n--) { if (m1np[i][p]+m2>0) break; } which could break for p >= 18 or i >= 7. Fortunately this code is never called due to the way the exp function is implemented since values having exponent less than -55 return 1.0. Make sure that if it gets called in future, it is trapped.
* Move more constants into static variablesSiddhesh Poyarekar2013-01-024-22/+8
| | | | Code cleanup.
* Add script to update copyright notices and reformat some to facilitate its use.Joseph Myers2013-01-0128-32/+28
|
* Favour normal numbersSiddhesh Poyarekar2012-12-311-4/+5
|
* Demystify the magic number 134217729.0Siddhesh Poyarekar2012-12-293-8/+10
| | | | | The number 134217729.0 gets used in various places in e_pow.c but there is no explanation of what that number is. Add that explanation.
* Remove unnecessary variable mptwoim1Siddhesh Poyarekar2012-12-282-27/+3
| | | | Code cleanup.
* Replace more constants with their valuesSiddhesh Poyarekar2012-12-285-97/+61
| | | | Code cleanup.
* Replace constants with preprocessor definesSiddhesh Poyarekar2012-12-285-120/+28
| | | | libm Code cleanup.
* Remove redundant __mpexp_nnSiddhesh Poyarekar2012-12-272-28/+2
| | | | It's an array that stores integral float values of the offset.
* Move mpone out to a global constSiddhesh Poyarekar2012-12-277-68/+11
| | | | Code cleanup.
* Include stdlib.h in sysdeps/ieee754/ldbl-128ibm/x2y2m1l.c.Joseph Myers2012-12-051-0/+1
|
* Fix ldbl-128ibm "set but not used" warnings.Joseph Myers2012-12-045-5/+9
|
* Fix "conflicting types for built-in function" warnings from nldbl-*.c.Joseph Myers2012-12-041-0/+99
|
* Fix powl inaccuracy for ldbl-128ibm (bug 14914).Joseph Myers2012-12-041-7/+7
|
* Use hex float 64-bit values in ldbl-96 asinl (bug 14803).Joseph Myers2012-11-281-3/+6
|
* Fix ldbl-128ibm atanl spurious underflows (bug 14871).Joseph Myers2012-11-221-0/+16
|
* Fix sign of inexact zero results for ldbl-128ibm fmal.Joseph Myers2012-11-221-1/+7
|
* Fix ldbl-128ibm powl spurious underflows.Joseph Myers2012-11-221-1/+5
|
* Fix ldbl-128ibm hypotl internal underflows (bug 14869).Joseph Myers2012-11-221-1/+1
|
* Fix ldbl-128ibm hypotl inaccuracy for arguments with large ratio (bug 14868).Joseph Myers2012-11-221-1/+1
|
* Fix expm1l spurious underflows for ldbl-128ibm.Joseph Myers2012-11-221-1/+1
|
* Fix set-but-not-used warnings in ldbl-128 nearbyintl, rintl.Joseph Myers2012-11-202-2/+2
|
* Fix spurious underflows in ldbl-128 atan implementation.David S. Miller2012-11-191-0/+17
| | | | | | | | | | 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.
* Correct tinyness handling in long-double and float y0/y1.David S. Miller2012-11-184-2/+7
| | | | | | | | | | | | | | | | | 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.
* Fix BZ #14811 for ldbl-128 too.David S. Miller2012-11-161-1/+5
| | | | | | | [BZ #14811] * sysdeps/ieee754/ldbl-128/e_powl.c (__ieee754_powl): Saturate nonzero exponents with absolute value below 0x1p-128 to +/- 0x1p-128.
* Don't generate underflow for very small values in log1pl.David S. Miller2012-11-161-0/+6
| | | | | * sysdeps/ieee754/ldbl-128/s_log1pl.c (__log1pl): If xm1 is smaller than LDBL_EPSILON/2.0L, just return xm1.