mirror of
https://github.com/mozilla/gecko-dev.git
synced 2024-11-23 12:51:06 +00:00
Bug 1804298: Update the fdlibm patches to be vendored by mach vendor r=arai
For mach vendor, patches are relative to the directory the third party code is vendored into. Differential Revision: https://phabricator.services.mozilla.com/D163971
This commit is contained in:
parent
ead4517cc3
commit
bec1a99cde
@ -1,6 +1,6 @@
|
||||
diff --git a/modules/fdlibm/src/fdlibm.h b/modules/fdlibm/src/fdlibm.h
|
||||
--- a/modules/fdlibm/src/fdlibm.h
|
||||
+++ b/modules/fdlibm/src/fdlibm.h
|
||||
diff --git a/fdlibm.h b/fdlibm.h
|
||||
--- a/fdlibm.h
|
||||
+++ b/fdlibm.h
|
||||
@@ -12,504 +12,51 @@
|
||||
/*
|
||||
* from: @(#)fdlibm.h 5.1 93/09/24
|
||||
|
@ -1,6 +1,6 @@
|
||||
diff --git a/modules/fdlibm/src/fdlibm.h b/modules/fdlibm/src/fdlibm.h
|
||||
--- a/modules/fdlibm/src/fdlibm.h
|
||||
+++ b/modules/fdlibm/src/fdlibm.h
|
||||
diff --git a/fdlibm.h b/fdlibm.h
|
||||
--- a/fdlibm.h
|
||||
+++ b/fdlibm.h
|
||||
@@ -9,18 +9,18 @@
|
||||
* ====================================================
|
||||
*/
|
||||
|
@ -1,6 +1,6 @@
|
||||
diff --git a/modules/fdlibm/src/fdlibm.h b/modules/fdlibm/src/fdlibm.h
|
||||
--- a/modules/fdlibm/src/fdlibm.h
|
||||
+++ b/modules/fdlibm/src/fdlibm.h
|
||||
diff --git a/fdlibm.h b/fdlibm.h
|
||||
--- a/fdlibm.h
|
||||
+++ b/fdlibm.h
|
||||
@@ -12,16 +12,18 @@
|
||||
/*
|
||||
* from: @(#)fdlibm.h 5.1 93/09/24
|
||||
|
@ -1,6 +1,6 @@
|
||||
diff --git a/modules/fdlibm/src/e_acos.cpp b/modules/fdlibm/src/e_acos.cpp
|
||||
--- a/modules/fdlibm/src/e_acos.cpp
|
||||
+++ b/modules/fdlibm/src/e_acos.cpp
|
||||
diff --git a/e_acos.cpp b/e_acos.cpp
|
||||
--- a/e_acos.cpp
|
||||
+++ b/e_acos.cpp
|
||||
@@ -35,17 +35,16 @@ __FBSDID("$FreeBSD$");
|
||||
* if x is NaN, return x itself;
|
||||
* if |x|>1, return NaN with invalid signal.
|
||||
@ -19,9 +19,9 @@ diff --git a/modules/fdlibm/src/e_acos.cpp b/modules/fdlibm/src/e_acos.cpp
|
||||
pio2_hi = 1.57079632679489655800e+00; /* 0x3FF921FB, 0x54442D18 */
|
||||
static volatile double
|
||||
pio2_lo = 6.12323399573676603587e-17; /* 0x3C91A626, 0x33145C07 */
|
||||
diff --git a/modules/fdlibm/src/e_acosh.cpp b/modules/fdlibm/src/e_acosh.cpp
|
||||
--- a/modules/fdlibm/src/e_acosh.cpp
|
||||
+++ b/modules/fdlibm/src/e_acosh.cpp
|
||||
diff --git a/e_acosh.cpp b/e_acosh.cpp
|
||||
--- a/e_acosh.cpp
|
||||
+++ b/e_acosh.cpp
|
||||
@@ -26,17 +26,16 @@ __FBSDID("$FreeBSD$");
|
||||
*
|
||||
* Special cases:
|
||||
@ -40,9 +40,9 @@ diff --git a/modules/fdlibm/src/e_acosh.cpp b/modules/fdlibm/src/e_acosh.cpp
|
||||
|
||||
double
|
||||
__ieee754_acosh(double x)
|
||||
diff --git a/modules/fdlibm/src/e_asin.cpp b/modules/fdlibm/src/e_asin.cpp
|
||||
--- a/modules/fdlibm/src/e_asin.cpp
|
||||
+++ b/modules/fdlibm/src/e_asin.cpp
|
||||
diff --git a/e_asin.cpp b/e_asin.cpp
|
||||
--- a/e_asin.cpp
|
||||
+++ b/e_asin.cpp
|
||||
@@ -41,17 +41,16 @@ __FBSDID("$FreeBSD$");
|
||||
* Special cases:
|
||||
* if x is NaN, return x itself;
|
||||
@ -61,9 +61,9 @@ diff --git a/modules/fdlibm/src/e_asin.cpp b/modules/fdlibm/src/e_asin.cpp
|
||||
pio2_hi = 1.57079632679489655800e+00, /* 0x3FF921FB, 0x54442D18 */
|
||||
pio2_lo = 6.12323399573676603587e-17, /* 0x3C91A626, 0x33145C07 */
|
||||
pio4_hi = 7.85398163397448278999e-01, /* 0x3FE921FB, 0x54442D18 */
|
||||
diff --git a/modules/fdlibm/src/e_atan2.cpp b/modules/fdlibm/src/e_atan2.cpp
|
||||
--- a/modules/fdlibm/src/e_atan2.cpp
|
||||
+++ b/modules/fdlibm/src/e_atan2.cpp
|
||||
diff --git a/e_atan2.cpp b/e_atan2.cpp
|
||||
--- a/e_atan2.cpp
|
||||
+++ b/e_atan2.cpp
|
||||
@@ -39,17 +39,16 @@ __FBSDID("$FreeBSD$");
|
||||
* The hexadecimal values are the intended ones for the following
|
||||
* constants. The decimal values may be used, provided that the
|
||||
@ -82,9 +82,9 @@ diff --git a/modules/fdlibm/src/e_atan2.cpp b/modules/fdlibm/src/e_atan2.cpp
|
||||
zero = 0.0,
|
||||
pi_o_4 = 7.8539816339744827900E-01, /* 0x3FE921FB, 0x54442D18 */
|
||||
pi_o_2 = 1.5707963267948965580E+00, /* 0x3FF921FB, 0x54442D18 */
|
||||
diff --git a/modules/fdlibm/src/e_atanh.cpp b/modules/fdlibm/src/e_atanh.cpp
|
||||
--- a/modules/fdlibm/src/e_atanh.cpp
|
||||
+++ b/modules/fdlibm/src/e_atanh.cpp
|
||||
diff --git a/e_atanh.cpp b/e_atanh.cpp
|
||||
--- a/e_atanh.cpp
|
||||
+++ b/e_atanh.cpp
|
||||
@@ -30,17 +30,16 @@ __FBSDID("$FreeBSD$");
|
||||
* atanh(x) is NaN if |x| > 1 with signal;
|
||||
* atanh(NaN) is that NaN with no signal;
|
||||
@ -103,9 +103,9 @@ diff --git a/modules/fdlibm/src/e_atanh.cpp b/modules/fdlibm/src/e_atanh.cpp
|
||||
double
|
||||
__ieee754_atanh(double x)
|
||||
{
|
||||
diff --git a/modules/fdlibm/src/e_cosh.cpp b/modules/fdlibm/src/e_cosh.cpp
|
||||
--- a/modules/fdlibm/src/e_cosh.cpp
|
||||
+++ b/modules/fdlibm/src/e_cosh.cpp
|
||||
diff --git a/e_cosh.cpp b/e_cosh.cpp
|
||||
--- a/e_cosh.cpp
|
||||
+++ b/e_cosh.cpp
|
||||
@@ -32,17 +32,16 @@ __FBSDID("$FreeBSD$");
|
||||
*
|
||||
* Special cases:
|
||||
@ -124,9 +124,9 @@ diff --git a/modules/fdlibm/src/e_cosh.cpp b/modules/fdlibm/src/e_cosh.cpp
|
||||
__ieee754_cosh(double x)
|
||||
{
|
||||
double t,w;
|
||||
diff --git a/modules/fdlibm/src/e_exp.cpp b/modules/fdlibm/src/e_exp.cpp
|
||||
--- a/modules/fdlibm/src/e_exp.cpp
|
||||
+++ b/modules/fdlibm/src/e_exp.cpp
|
||||
diff --git a/e_exp.cpp b/e_exp.cpp
|
||||
--- a/e_exp.cpp
|
||||
+++ b/e_exp.cpp
|
||||
@@ -73,17 +73,16 @@ __FBSDID("$FreeBSD$");
|
||||
* The hexadecimal values are the intended ones for the following
|
||||
* constants. The decimal values may be used, provided that the
|
||||
@ -145,9 +145,9 @@ diff --git a/modules/fdlibm/src/e_exp.cpp b/modules/fdlibm/src/e_exp.cpp
|
||||
o_threshold= 7.09782712893383973096e+02, /* 0x40862E42, 0xFEFA39EF */
|
||||
u_threshold= -7.45133219101941108420e+02, /* 0xc0874910, 0xD52D3051 */
|
||||
ln2HI[2] ={ 6.93147180369123816490e-01, /* 0x3fe62e42, 0xfee00000 */
|
||||
diff --git a/modules/fdlibm/src/e_hypot.cpp b/modules/fdlibm/src/e_hypot.cpp
|
||||
--- a/modules/fdlibm/src/e_hypot.cpp
|
||||
+++ b/modules/fdlibm/src/e_hypot.cpp
|
||||
diff --git a/e_hypot.cpp b/e_hypot.cpp
|
||||
--- a/e_hypot.cpp
|
||||
+++ b/e_hypot.cpp
|
||||
@@ -43,17 +43,16 @@ __FBSDID("$FreeBSD$");
|
||||
*
|
||||
* Accuracy:
|
||||
@ -166,9 +166,9 @@ diff --git a/modules/fdlibm/src/e_hypot.cpp b/modules/fdlibm/src/e_hypot.cpp
|
||||
double a,b,t1,t2,y1,y2,w;
|
||||
int32_t j,k,ha,hb;
|
||||
|
||||
diff --git a/modules/fdlibm/src/e_log.cpp b/modules/fdlibm/src/e_log.cpp
|
||||
--- a/modules/fdlibm/src/e_log.cpp
|
||||
+++ b/modules/fdlibm/src/e_log.cpp
|
||||
diff --git a/e_log.cpp b/e_log.cpp
|
||||
--- a/e_log.cpp
|
||||
+++ b/e_log.cpp
|
||||
@@ -62,17 +62,16 @@ __FBSDID("$FreeBSD$");
|
||||
* The hexadecimal values are the intended ones for the following
|
||||
* constants. The decimal values may be used, provided that the
|
||||
@ -187,9 +187,9 @@ diff --git a/modules/fdlibm/src/e_log.cpp b/modules/fdlibm/src/e_log.cpp
|
||||
two54 = 1.80143985094819840000e+16, /* 43500000 00000000 */
|
||||
Lg1 = 6.666666666666735130e-01, /* 3FE55555 55555593 */
|
||||
Lg2 = 3.999999999940941908e-01, /* 3FD99999 9997FA04 */
|
||||
diff --git a/modules/fdlibm/src/e_log10.cpp b/modules/fdlibm/src/e_log10.cpp
|
||||
--- a/modules/fdlibm/src/e_log10.cpp
|
||||
+++ b/modules/fdlibm/src/e_log10.cpp
|
||||
diff --git a/e_log10.cpp b/e_log10.cpp
|
||||
--- a/e_log10.cpp
|
||||
+++ b/e_log10.cpp
|
||||
@@ -19,17 +19,16 @@ __FBSDID("$FreeBSD$");
|
||||
* comments.
|
||||
*
|
||||
@ -208,9 +208,9 @@ diff --git a/modules/fdlibm/src/e_log10.cpp b/modules/fdlibm/src/e_log10.cpp
|
||||
ivln10hi = 4.34294481878168880939e-01, /* 0x3fdbcb7b, 0x15200000 */
|
||||
ivln10lo = 2.50829467116452752298e-11, /* 0x3dbb9438, 0xca9aadd5 */
|
||||
log10_2hi = 3.01029995663611771306e-01, /* 0x3FD34413, 0x509F6000 */
|
||||
diff --git a/modules/fdlibm/src/e_log2.cpp b/modules/fdlibm/src/e_log2.cpp
|
||||
--- a/modules/fdlibm/src/e_log2.cpp
|
||||
+++ b/modules/fdlibm/src/e_log2.cpp
|
||||
diff --git a/e_log2.cpp b/e_log2.cpp
|
||||
--- a/e_log2.cpp
|
||||
+++ b/e_log2.cpp
|
||||
@@ -21,17 +21,16 @@ __FBSDID("$FreeBSD$");
|
||||
* This reduces x to {k, 1+f} exactly as in e_log.c, then calls the kernel,
|
||||
* then does the combining and scaling steps
|
||||
@ -229,9 +229,9 @@ diff --git a/modules/fdlibm/src/e_log2.cpp b/modules/fdlibm/src/e_log2.cpp
|
||||
ivln2hi = 1.44269504072144627571e+00, /* 0x3ff71547, 0x65200000 */
|
||||
ivln2lo = 1.67517131648865118353e-10; /* 0x3de705fc, 0x2eefa200 */
|
||||
|
||||
diff --git a/modules/fdlibm/src/e_pow.cpp b/modules/fdlibm/src/e_pow.cpp
|
||||
--- a/modules/fdlibm/src/e_pow.cpp
|
||||
+++ b/modules/fdlibm/src/e_pow.cpp
|
||||
diff --git a/e_pow.cpp b/e_pow.cpp
|
||||
--- a/e_pow.cpp
|
||||
+++ b/e_pow.cpp
|
||||
@@ -53,17 +53,16 @@ __FBSDID("$FreeBSD$");
|
||||
* Constants :
|
||||
* The hexadecimal values are the intended ones for the following
|
||||
@ -250,9 +250,9 @@ diff --git a/modules/fdlibm/src/e_pow.cpp b/modules/fdlibm/src/e_pow.cpp
|
||||
dp_l[] = { 0.0, 1.35003920212974897128e-08,}, /* 0x3E4CFDEB, 0x43CFD006 */
|
||||
zero = 0.0,
|
||||
half = 0.5,
|
||||
diff --git a/modules/fdlibm/src/e_rem_pio2.cpp b/modules/fdlibm/src/e_rem_pio2.cpp
|
||||
--- a/modules/fdlibm/src/e_rem_pio2.cpp
|
||||
+++ b/modules/fdlibm/src/e_rem_pio2.cpp
|
||||
diff --git a/e_rem_pio2.cpp b/e_rem_pio2.cpp
|
||||
--- a/e_rem_pio2.cpp
|
||||
+++ b/e_rem_pio2.cpp
|
||||
@@ -19,17 +19,16 @@ __FBSDID("$FreeBSD$");
|
||||
/* __ieee754_rem_pio2(x,y)
|
||||
*
|
||||
@ -271,9 +271,9 @@ diff --git a/modules/fdlibm/src/e_rem_pio2.cpp b/modules/fdlibm/src/e_rem_pio2.c
|
||||
* pio2_1t: pi/2 - pio2_1
|
||||
* pio2_2: second 33 bit of pi/2
|
||||
* pio2_2t: pi/2 - (pio2_1+pio2_2)
|
||||
diff --git a/modules/fdlibm/src/e_sinh.cpp b/modules/fdlibm/src/e_sinh.cpp
|
||||
--- a/modules/fdlibm/src/e_sinh.cpp
|
||||
+++ b/modules/fdlibm/src/e_sinh.cpp
|
||||
diff --git a/e_sinh.cpp b/e_sinh.cpp
|
||||
--- a/e_sinh.cpp
|
||||
+++ b/e_sinh.cpp
|
||||
@@ -29,17 +29,16 @@ __FBSDID("$FreeBSD$");
|
||||
*
|
||||
* Special cases:
|
||||
@ -292,9 +292,9 @@ diff --git a/modules/fdlibm/src/e_sinh.cpp b/modules/fdlibm/src/e_sinh.cpp
|
||||
__ieee754_sinh(double x)
|
||||
{
|
||||
double t,h;
|
||||
diff --git a/modules/fdlibm/src/k_cos.cpp b/modules/fdlibm/src/k_cos.cpp
|
||||
--- a/modules/fdlibm/src/k_cos.cpp
|
||||
+++ b/modules/fdlibm/src/k_cos.cpp
|
||||
diff --git a/k_cos.cpp b/k_cos.cpp
|
||||
--- a/k_cos.cpp
|
||||
+++ b/k_cos.cpp
|
||||
@@ -48,17 +48,16 @@ __FBSDID("$FreeBSD$");
|
||||
* and tmp having the same precision as x. If they have extra
|
||||
* precision due to compiler bugs, then the extra precision is
|
||||
@ -313,9 +313,9 @@ diff --git a/modules/fdlibm/src/k_cos.cpp b/modules/fdlibm/src/k_cos.cpp
|
||||
C2 = -1.38888888888741095749e-03, /* 0xBF56C16C, 0x16C15177 */
|
||||
C3 = 2.48015872894767294178e-05, /* 0x3EFA01A0, 0x19CB1590 */
|
||||
C4 = -2.75573143513906633035e-07, /* 0xBE927E4F, 0x809C52AD */
|
||||
diff --git a/modules/fdlibm/src/k_exp.cpp b/modules/fdlibm/src/k_exp.cpp
|
||||
--- a/modules/fdlibm/src/k_exp.cpp
|
||||
+++ b/modules/fdlibm/src/k_exp.cpp
|
||||
diff --git a/k_exp.cpp b/k_exp.cpp
|
||||
--- a/k_exp.cpp
|
||||
+++ b/k_exp.cpp
|
||||
@@ -26,17 +26,16 @@
|
||||
* SUCH DAMAGE.
|
||||
*/
|
||||
@ -334,9 +334,9 @@ diff --git a/modules/fdlibm/src/k_exp.cpp b/modules/fdlibm/src/k_exp.cpp
|
||||
/*
|
||||
* Compute exp(x), scaled to avoid spurious overflow. An exponent is
|
||||
* returned separately in 'expt'.
|
||||
diff --git a/modules/fdlibm/src/k_rem_pio2.cpp b/modules/fdlibm/src/k_rem_pio2.cpp
|
||||
--- a/modules/fdlibm/src/k_rem_pio2.cpp
|
||||
+++ b/modules/fdlibm/src/k_rem_pio2.cpp
|
||||
diff --git a/k_rem_pio2.cpp b/k_rem_pio2.cpp
|
||||
--- a/k_rem_pio2.cpp
|
||||
+++ b/k_rem_pio2.cpp
|
||||
@@ -126,17 +126,16 @@ __FBSDID("$FreeBSD$");
|
||||
* The hexadecimal values are the intended ones for the following
|
||||
* constants. The decimal values may be used, provided that the
|
||||
@ -355,9 +355,9 @@ diff --git a/modules/fdlibm/src/k_rem_pio2.cpp b/modules/fdlibm/src/k_rem_pio2.c
|
||||
* Table of constants for 2/pi, 396 Hex digits (476 decimal) of 2/pi
|
||||
*
|
||||
* integer array, contains the (24*i)-th to (24*i+23)-th
|
||||
diff --git a/modules/fdlibm/src/k_sin.cpp b/modules/fdlibm/src/k_sin.cpp
|
||||
--- a/modules/fdlibm/src/k_sin.cpp
|
||||
+++ b/modules/fdlibm/src/k_sin.cpp
|
||||
diff --git a/k_sin.cpp b/k_sin.cpp
|
||||
--- a/k_sin.cpp
|
||||
+++ b/k_sin.cpp
|
||||
@@ -39,17 +39,16 @@ __FBSDID("$FreeBSD$");
|
||||
* ~ sin(x) + (1-x*x/2)*y
|
||||
* For better accuracy, let
|
||||
@ -376,9 +376,9 @@ diff --git a/modules/fdlibm/src/k_sin.cpp b/modules/fdlibm/src/k_sin.cpp
|
||||
S2 = 8.33333333332248946124e-03, /* 0x3F811111, 0x1110F8A6 */
|
||||
S3 = -1.98412698298579493134e-04, /* 0xBF2A01A0, 0x19C161D5 */
|
||||
S4 = 2.75573137070700676789e-06, /* 0x3EC71DE3, 0x57B1FE7D */
|
||||
diff --git a/modules/fdlibm/src/k_tan.cpp b/modules/fdlibm/src/k_tan.cpp
|
||||
--- a/modules/fdlibm/src/k_tan.cpp
|
||||
+++ b/modules/fdlibm/src/k_tan.cpp
|
||||
diff --git a/k_tan.cpp b/k_tan.cpp
|
||||
--- a/k_tan.cpp
|
||||
+++ b/k_tan.cpp
|
||||
@@ -44,17 +44,16 @@ __FBSDID("$FreeBSD$");
|
||||
* 3 2
|
||||
* tan(x+y) = x + (T1*x + (x *(r+y)+y))
|
||||
@ -397,9 +397,9 @@ diff --git a/modules/fdlibm/src/k_tan.cpp b/modules/fdlibm/src/k_tan.cpp
|
||||
2.18694882948595424599e-02, /* 3F9664F4, 8406D637 */
|
||||
8.86323982359930005737e-03, /* 3F8226E3, E96E8493 */
|
||||
3.59207910759131235356e-03, /* 3F6D6D22, C9560328 */
|
||||
diff --git a/modules/fdlibm/src/math_private.h b/modules/fdlibm/src/math_private.h
|
||||
--- a/modules/fdlibm/src/math_private.h
|
||||
+++ b/modules/fdlibm/src/math_private.h
|
||||
diff --git a/math_private.h b/math_private.h
|
||||
--- a/math_private.h
|
||||
+++ b/math_private.h
|
||||
@@ -15,16 +15,18 @@
|
||||
*/
|
||||
|
||||
@ -419,9 +419,9 @@ diff --git a/modules/fdlibm/src/math_private.h b/modules/fdlibm/src/math_private
|
||||
* to dig two 32 bit words out of the 64 bit IEEE floating point
|
||||
* value. That is non-ANSI, and, moreover, the gcc instruction
|
||||
* scheduler gets it wrong. We instead use the following macros.
|
||||
diff --git a/modules/fdlibm/src/s_asinh.cpp b/modules/fdlibm/src/s_asinh.cpp
|
||||
--- a/modules/fdlibm/src/s_asinh.cpp
|
||||
+++ b/modules/fdlibm/src/s_asinh.cpp
|
||||
diff --git a/s_asinh.cpp b/s_asinh.cpp
|
||||
--- a/s_asinh.cpp
|
||||
+++ b/s_asinh.cpp
|
||||
@@ -21,17 +21,16 @@ __FBSDID("$FreeBSD$");
|
||||
* asinh(x) := x if 1+x*x=1,
|
||||
* := sign(x)*(log(x)+ln2)) for large |x|, else
|
||||
@ -440,9 +440,9 @@ diff --git a/modules/fdlibm/src/s_asinh.cpp b/modules/fdlibm/src/s_asinh.cpp
|
||||
huge= 1.00000000000000000000e+300;
|
||||
|
||||
double
|
||||
diff --git a/modules/fdlibm/src/s_atan.cpp b/modules/fdlibm/src/s_atan.cpp
|
||||
--- a/modules/fdlibm/src/s_atan.cpp
|
||||
+++ b/modules/fdlibm/src/s_atan.cpp
|
||||
diff --git a/s_atan.cpp b/s_atan.cpp
|
||||
--- a/s_atan.cpp
|
||||
+++ b/s_atan.cpp
|
||||
@@ -30,17 +30,16 @@ __FBSDID("$FreeBSD$");
|
||||
* The hexadecimal values are the intended ones for the following
|
||||
* constants. The decimal values may be used, provided that the
|
||||
@ -461,9 +461,9 @@ diff --git a/modules/fdlibm/src/s_atan.cpp b/modules/fdlibm/src/s_atan.cpp
|
||||
9.82793723247329054082e-01, /* atan(1.5)hi 0x3FEF730B, 0xD281F69B */
|
||||
1.57079632679489655800e+00, /* atan(inf)hi 0x3FF921FB, 0x54442D18 */
|
||||
};
|
||||
diff --git a/modules/fdlibm/src/s_cbrt.cpp b/modules/fdlibm/src/s_cbrt.cpp
|
||||
--- a/modules/fdlibm/src/s_cbrt.cpp
|
||||
+++ b/modules/fdlibm/src/s_cbrt.cpp
|
||||
diff --git a/s_cbrt.cpp b/s_cbrt.cpp
|
||||
--- a/s_cbrt.cpp
|
||||
+++ b/s_cbrt.cpp
|
||||
@@ -11,17 +11,16 @@
|
||||
*
|
||||
* Optimized by Bruce D. Evans.
|
||||
@ -482,9 +482,9 @@ diff --git a/modules/fdlibm/src/s_cbrt.cpp b/modules/fdlibm/src/s_cbrt.cpp
|
||||
static const u_int32_t
|
||||
B1 = 715094163, /* B1 = (1023-1023/3-0.03306235651)*2**20 */
|
||||
B2 = 696219795; /* B2 = (1023-1023/3-54/3-0.03306235651)*2**20 */
|
||||
diff --git a/modules/fdlibm/src/s_ceil.cpp b/modules/fdlibm/src/s_ceil.cpp
|
||||
--- a/modules/fdlibm/src/s_ceil.cpp
|
||||
+++ b/modules/fdlibm/src/s_ceil.cpp
|
||||
diff --git a/s_ceil.cpp b/s_ceil.cpp
|
||||
--- a/s_ceil.cpp
|
||||
+++ b/s_ceil.cpp
|
||||
@@ -19,17 +19,16 @@ __FBSDID("$FreeBSD$");
|
||||
* Method:
|
||||
* Bit twiddling.
|
||||
@ -503,9 +503,9 @@ diff --git a/modules/fdlibm/src/s_ceil.cpp b/modules/fdlibm/src/s_ceil.cpp
|
||||
ceil(double x)
|
||||
{
|
||||
int32_t i0,i1,j0;
|
||||
diff --git a/modules/fdlibm/src/s_ceilf.cpp b/modules/fdlibm/src/s_ceilf.cpp
|
||||
--- a/modules/fdlibm/src/s_ceilf.cpp
|
||||
+++ b/modules/fdlibm/src/s_ceilf.cpp
|
||||
diff --git a/s_ceilf.cpp b/s_ceilf.cpp
|
||||
--- a/s_ceilf.cpp
|
||||
+++ b/s_ceilf.cpp
|
||||
@@ -11,17 +11,16 @@
|
||||
* software is freely granted, provided that this notice
|
||||
* is preserved.
|
||||
@ -524,9 +524,9 @@ diff --git a/modules/fdlibm/src/s_ceilf.cpp b/modules/fdlibm/src/s_ceilf.cpp
|
||||
ceilf(float x)
|
||||
{
|
||||
int32_t i0,j0;
|
||||
diff --git a/modules/fdlibm/src/s_copysign.cpp b/modules/fdlibm/src/s_copysign.cpp
|
||||
--- a/modules/fdlibm/src/s_copysign.cpp
|
||||
+++ b/modules/fdlibm/src/s_copysign.cpp
|
||||
diff --git a/s_copysign.cpp b/s_copysign.cpp
|
||||
--- a/s_copysign.cpp
|
||||
+++ b/s_copysign.cpp
|
||||
@@ -14,17 +14,16 @@
|
||||
__FBSDID("$FreeBSD$");
|
||||
|
||||
@ -545,9 +545,9 @@ diff --git a/modules/fdlibm/src/s_copysign.cpp b/modules/fdlibm/src/s_copysign.c
|
||||
u_int32_t hx,hy;
|
||||
GET_HIGH_WORD(hx,x);
|
||||
GET_HIGH_WORD(hy,y);
|
||||
diff --git a/modules/fdlibm/src/s_cos.cpp b/modules/fdlibm/src/s_cos.cpp
|
||||
--- a/modules/fdlibm/src/s_cos.cpp
|
||||
+++ b/modules/fdlibm/src/s_cos.cpp
|
||||
diff --git a/s_cos.cpp b/s_cos.cpp
|
||||
--- a/s_cos.cpp
|
||||
+++ b/s_cos.cpp
|
||||
@@ -41,17 +41,16 @@ __FBSDID("$FreeBSD$");
|
||||
* trig(NaN) is that NaN;
|
||||
*
|
||||
@ -566,9 +566,9 @@ diff --git a/modules/fdlibm/src/s_cos.cpp b/modules/fdlibm/src/s_cos.cpp
|
||||
cos(double x)
|
||||
{
|
||||
double y[2],z=0.0;
|
||||
diff --git a/modules/fdlibm/src/s_expm1.cpp b/modules/fdlibm/src/s_expm1.cpp
|
||||
--- a/modules/fdlibm/src/s_expm1.cpp
|
||||
+++ b/modules/fdlibm/src/s_expm1.cpp
|
||||
diff --git a/s_expm1.cpp b/s_expm1.cpp
|
||||
--- a/s_expm1.cpp
|
||||
+++ b/s_expm1.cpp
|
||||
@@ -105,17 +105,16 @@ __FBSDID("$FreeBSD$");
|
||||
* The hexadecimal values are the intended ones for the following
|
||||
* constants. The decimal values may be used, provided that the
|
||||
@ -587,9 +587,9 @@ diff --git a/modules/fdlibm/src/s_expm1.cpp b/modules/fdlibm/src/s_expm1.cpp
|
||||
o_threshold = 7.09782712893383973096e+02,/* 0x40862E42, 0xFEFA39EF */
|
||||
ln2_hi = 6.93147180369123816490e-01,/* 0x3fe62e42, 0xfee00000 */
|
||||
ln2_lo = 1.90821492927058770002e-10,/* 0x3dea39ef, 0x35793c76 */
|
||||
diff --git a/modules/fdlibm/src/s_fabs.cpp b/modules/fdlibm/src/s_fabs.cpp
|
||||
--- a/modules/fdlibm/src/s_fabs.cpp
|
||||
+++ b/modules/fdlibm/src/s_fabs.cpp
|
||||
diff --git a/s_fabs.cpp b/s_fabs.cpp
|
||||
--- a/s_fabs.cpp
|
||||
+++ b/s_fabs.cpp
|
||||
@@ -12,17 +12,16 @@
|
||||
|
||||
#include <sys/cdefs.h>
|
||||
@ -608,9 +608,9 @@ diff --git a/modules/fdlibm/src/s_fabs.cpp b/modules/fdlibm/src/s_fabs.cpp
|
||||
u_int32_t high;
|
||||
GET_HIGH_WORD(high,x);
|
||||
SET_HIGH_WORD(x,high&0x7fffffff);
|
||||
diff --git a/modules/fdlibm/src/s_floor.cpp b/modules/fdlibm/src/s_floor.cpp
|
||||
--- a/modules/fdlibm/src/s_floor.cpp
|
||||
+++ b/modules/fdlibm/src/s_floor.cpp
|
||||
diff --git a/s_floor.cpp b/s_floor.cpp
|
||||
--- a/s_floor.cpp
|
||||
+++ b/s_floor.cpp
|
||||
@@ -19,17 +19,16 @@ __FBSDID("$FreeBSD$");
|
||||
* Method:
|
||||
* Bit twiddling.
|
||||
@ -629,9 +629,9 @@ diff --git a/modules/fdlibm/src/s_floor.cpp b/modules/fdlibm/src/s_floor.cpp
|
||||
floor(double x)
|
||||
{
|
||||
int32_t i0,i1,j0;
|
||||
diff --git a/modules/fdlibm/src/s_floorf.cpp b/modules/fdlibm/src/s_floorf.cpp
|
||||
--- a/modules/fdlibm/src/s_floorf.cpp
|
||||
+++ b/modules/fdlibm/src/s_floorf.cpp
|
||||
diff --git a/s_floorf.cpp b/s_floorf.cpp
|
||||
--- a/s_floorf.cpp
|
||||
+++ b/s_floorf.cpp
|
||||
@@ -20,17 +20,16 @@ __FBSDID("$FreeBSD$");
|
||||
* floorf(x)
|
||||
* Return x rounded toward -inf to integral value
|
||||
@ -650,9 +650,9 @@ diff --git a/modules/fdlibm/src/s_floorf.cpp b/modules/fdlibm/src/s_floorf.cpp
|
||||
floorf(float x)
|
||||
{
|
||||
int32_t i0,j0;
|
||||
diff --git a/modules/fdlibm/src/s_log1p.cpp b/modules/fdlibm/src/s_log1p.cpp
|
||||
--- a/modules/fdlibm/src/s_log1p.cpp
|
||||
+++ b/modules/fdlibm/src/s_log1p.cpp
|
||||
diff --git a/s_log1p.cpp b/s_log1p.cpp
|
||||
--- a/s_log1p.cpp
|
||||
+++ b/s_log1p.cpp
|
||||
@@ -75,17 +75,16 @@ __FBSDID("$FreeBSD$");
|
||||
* if(u==1.0) return x ; else
|
||||
* return log(u)*(x/(u-1.0));
|
||||
@ -671,9 +671,9 @@ diff --git a/modules/fdlibm/src/s_log1p.cpp b/modules/fdlibm/src/s_log1p.cpp
|
||||
two54 = 1.80143985094819840000e+16, /* 43500000 00000000 */
|
||||
Lp1 = 6.666666666666735130e-01, /* 3FE55555 55555593 */
|
||||
Lp2 = 3.999999999940941908e-01, /* 3FD99999 9997FA04 */
|
||||
diff --git a/modules/fdlibm/src/s_nearbyint.cpp b/modules/fdlibm/src/s_nearbyint.cpp
|
||||
--- a/modules/fdlibm/src/s_nearbyint.cpp
|
||||
+++ b/modules/fdlibm/src/s_nearbyint.cpp
|
||||
diff --git a/s_nearbyint.cpp b/s_nearbyint.cpp
|
||||
--- a/s_nearbyint.cpp
|
||||
+++ b/s_nearbyint.cpp
|
||||
@@ -25,17 +25,18 @@
|
||||
* OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
|
||||
* SUCH DAMAGE.
|
||||
@ -694,9 +694,9 @@ diff --git a/modules/fdlibm/src/s_nearbyint.cpp b/modules/fdlibm/src/s_nearbyint
|
||||
* because the only exception defined for rint() is overflow, and
|
||||
* rounding can't overflow as long as emax >= p.
|
||||
*
|
||||
diff --git a/modules/fdlibm/src/s_rint.cpp b/modules/fdlibm/src/s_rint.cpp
|
||||
--- a/modules/fdlibm/src/s_rint.cpp
|
||||
+++ b/modules/fdlibm/src/s_rint.cpp
|
||||
diff --git a/s_rint.cpp b/s_rint.cpp
|
||||
--- a/s_rint.cpp
|
||||
+++ b/s_rint.cpp
|
||||
@@ -20,17 +20,16 @@ __FBSDID("$FreeBSD$");
|
||||
* Method:
|
||||
* Using floating addition.
|
||||
@ -715,9 +715,9 @@ diff --git a/modules/fdlibm/src/s_rint.cpp b/modules/fdlibm/src/s_rint.cpp
|
||||
-4.50359962737049600000e+15, /* 0xC3300000, 0x00000000 */
|
||||
};
|
||||
|
||||
diff --git a/modules/fdlibm/src/s_rintf.cpp b/modules/fdlibm/src/s_rintf.cpp
|
||||
--- a/modules/fdlibm/src/s_rintf.cpp
|
||||
+++ b/modules/fdlibm/src/s_rintf.cpp
|
||||
diff --git a/s_rintf.cpp b/s_rintf.cpp
|
||||
--- a/s_rintf.cpp
|
||||
+++ b/s_rintf.cpp
|
||||
@@ -14,17 +14,16 @@
|
||||
*/
|
||||
|
||||
@ -736,9 +736,9 @@ diff --git a/modules/fdlibm/src/s_rintf.cpp b/modules/fdlibm/src/s_rintf.cpp
|
||||
-8.3886080000e+06, /* 0xcb000000 */
|
||||
};
|
||||
|
||||
diff --git a/modules/fdlibm/src/s_scalbn.cpp b/modules/fdlibm/src/s_scalbn.cpp
|
||||
--- a/modules/fdlibm/src/s_scalbn.cpp
|
||||
+++ b/modules/fdlibm/src/s_scalbn.cpp
|
||||
diff --git a/s_scalbn.cpp b/s_scalbn.cpp
|
||||
--- a/s_scalbn.cpp
|
||||
+++ b/s_scalbn.cpp
|
||||
@@ -2,19 +2,20 @@
|
||||
* Copyright (c) 2005-2020 Rich Felker, et al.
|
||||
*
|
||||
@ -761,9 +761,9 @@ diff --git a/modules/fdlibm/src/s_scalbn.cpp b/modules/fdlibm/src/s_scalbn.cpp
|
||||
if (n > 1023) {
|
||||
y *= 0x1p1023;
|
||||
n -= 1023;
|
||||
diff --git a/modules/fdlibm/src/s_sin.cpp b/modules/fdlibm/src/s_sin.cpp
|
||||
--- a/modules/fdlibm/src/s_sin.cpp
|
||||
+++ b/modules/fdlibm/src/s_sin.cpp
|
||||
diff --git a/s_sin.cpp b/s_sin.cpp
|
||||
--- a/s_sin.cpp
|
||||
+++ b/s_sin.cpp
|
||||
@@ -41,17 +41,16 @@ __FBSDID("$FreeBSD$");
|
||||
* trig(NaN) is that NaN;
|
||||
*
|
||||
@ -782,9 +782,9 @@ diff --git a/modules/fdlibm/src/s_sin.cpp b/modules/fdlibm/src/s_sin.cpp
|
||||
sin(double x)
|
||||
{
|
||||
double y[2],z=0.0;
|
||||
diff --git a/modules/fdlibm/src/s_tan.cpp b/modules/fdlibm/src/s_tan.cpp
|
||||
--- a/modules/fdlibm/src/s_tan.cpp
|
||||
+++ b/modules/fdlibm/src/s_tan.cpp
|
||||
diff --git a/s_tan.cpp b/s_tan.cpp
|
||||
--- a/s_tan.cpp
|
||||
+++ b/s_tan.cpp
|
||||
@@ -40,17 +40,16 @@ __FBSDID("$FreeBSD$");
|
||||
* trig(NaN) is that NaN;
|
||||
*
|
||||
@ -803,9 +803,9 @@ diff --git a/modules/fdlibm/src/s_tan.cpp b/modules/fdlibm/src/s_tan.cpp
|
||||
tan(double x)
|
||||
{
|
||||
double y[2],z=0.0;
|
||||
diff --git a/modules/fdlibm/src/s_tanh.cpp b/modules/fdlibm/src/s_tanh.cpp
|
||||
--- a/modules/fdlibm/src/s_tanh.cpp
|
||||
+++ b/modules/fdlibm/src/s_tanh.cpp
|
||||
diff --git a/s_tanh.cpp b/s_tanh.cpp
|
||||
--- a/s_tanh.cpp
|
||||
+++ b/s_tanh.cpp
|
||||
@@ -34,17 +34,16 @@ __FBSDID("$FreeBSD$");
|
||||
*
|
||||
* Special cases:
|
||||
@ -824,9 +824,9 @@ diff --git a/modules/fdlibm/src/s_tanh.cpp b/modules/fdlibm/src/s_tanh.cpp
|
||||
double
|
||||
tanh(double x)
|
||||
{
|
||||
diff --git a/modules/fdlibm/src/s_trunc.cpp b/modules/fdlibm/src/s_trunc.cpp
|
||||
--- a/modules/fdlibm/src/s_trunc.cpp
|
||||
+++ b/modules/fdlibm/src/s_trunc.cpp
|
||||
diff --git a/s_trunc.cpp b/s_trunc.cpp
|
||||
--- a/s_trunc.cpp
|
||||
+++ b/s_trunc.cpp
|
||||
@@ -19,17 +19,16 @@ __FBSDID("$FreeBSD$");
|
||||
* Method:
|
||||
* Bit twiddling.
|
||||
@ -845,9 +845,9 @@ diff --git a/modules/fdlibm/src/s_trunc.cpp b/modules/fdlibm/src/s_trunc.cpp
|
||||
trunc(double x)
|
||||
{
|
||||
int32_t i0,i1,j0;
|
||||
diff --git a/modules/fdlibm/src/s_truncf.cpp b/modules/fdlibm/src/s_truncf.cpp
|
||||
--- a/modules/fdlibm/src/s_truncf.cpp
|
||||
+++ b/modules/fdlibm/src/s_truncf.cpp
|
||||
diff --git a/s_truncf.cpp b/s_truncf.cpp
|
||||
--- a/s_truncf.cpp
|
||||
+++ b/s_truncf.cpp
|
||||
@@ -17,17 +17,16 @@ __FBSDID("$FreeBSD$");
|
||||
* truncf(x)
|
||||
* Return x rounded toward 0 to integral value
|
||||
|
@ -1,6 +1,6 @@
|
||||
diff --git a/modules/fdlibm/src/math_private.h b/modules/fdlibm/src/math_private.h
|
||||
--- a/modules/fdlibm/src/math_private.h
|
||||
+++ b/modules/fdlibm/src/math_private.h
|
||||
diff --git a/math_private.h b/math_private.h
|
||||
--- a/math_private.h
|
||||
+++ b/math_private.h
|
||||
@@ -12,16 +12,17 @@
|
||||
/*
|
||||
* from: @(#)fdlibm.h 5.1 93/09/24
|
||||
|
@ -1,6 +1,6 @@
|
||||
diff --git a/modules/fdlibm/src/math_private.h b/modules/fdlibm/src/math_private.h
|
||||
--- a/modules/fdlibm/src/math_private.h
|
||||
+++ b/modules/fdlibm/src/math_private.h
|
||||
diff --git a/math_private.h b/math_private.h
|
||||
--- a/math_private.h
|
||||
+++ b/math_private.h
|
||||
@@ -14,52 +14,38 @@
|
||||
* $FreeBSD$
|
||||
*/
|
||||
|
@ -1,6 +1,6 @@
|
||||
diff --git a/modules/fdlibm/src/math_private.h b/modules/fdlibm/src/math_private.h
|
||||
--- a/modules/fdlibm/src/math_private.h
|
||||
+++ b/modules/fdlibm/src/math_private.h
|
||||
diff --git a/math_private.h b/math_private.h
|
||||
--- a/math_private.h
|
||||
+++ b/math_private.h
|
||||
@@ -867,16 +867,53 @@ irintl(long double x)
|
||||
#define __ieee754_j1f j1f
|
||||
#define __ieee754_y0f y0f
|
||||
|
@ -1,6 +1,6 @@
|
||||
diff --git a/modules/fdlibm/src/e_acos.cpp b/modules/fdlibm/src/e_acos.cpp
|
||||
--- a/modules/fdlibm/src/e_acos.cpp
|
||||
+++ b/modules/fdlibm/src/e_acos.cpp
|
||||
diff --git a/e_acos.cpp b/e_acos.cpp
|
||||
--- a/e_acos.cpp
|
||||
+++ b/e_acos.cpp
|
||||
@@ -99,12 +99,8 @@ __ieee754_acos(double x)
|
||||
c = (z-df*df)/(s+df);
|
||||
p = z*(pS0+z*(pS1+z*(pS2+z*(pS3+z*(pS4+z*pS5)))));
|
||||
@ -14,9 +14,9 @@ diff --git a/modules/fdlibm/src/e_acos.cpp b/modules/fdlibm/src/e_acos.cpp
|
||||
-#if LDBL_MANT_DIG == 53
|
||||
-__weak_reference(acos, acosl);
|
||||
-#endif
|
||||
diff --git a/modules/fdlibm/src/e_acosh.cpp b/modules/fdlibm/src/e_acosh.cpp
|
||||
--- a/modules/fdlibm/src/e_acosh.cpp
|
||||
+++ b/modules/fdlibm/src/e_acosh.cpp
|
||||
diff --git a/e_acosh.cpp b/e_acosh.cpp
|
||||
--- a/e_acosh.cpp
|
||||
+++ b/e_acosh.cpp
|
||||
@@ -56,12 +56,8 @@ __ieee754_acosh(double x)
|
||||
} else if (hx > 0x40000000) { /* 2**28 > x > 2 */
|
||||
t=x*x;
|
||||
@ -30,9 +30,9 @@ diff --git a/modules/fdlibm/src/e_acosh.cpp b/modules/fdlibm/src/e_acosh.cpp
|
||||
-#if LDBL_MANT_DIG == 53
|
||||
-__weak_reference(acosh, acoshl);
|
||||
-#endif
|
||||
diff --git a/modules/fdlibm/src/e_asin.cpp b/modules/fdlibm/src/e_asin.cpp
|
||||
--- a/modules/fdlibm/src/e_asin.cpp
|
||||
+++ b/modules/fdlibm/src/e_asin.cpp
|
||||
diff --git a/e_asin.cpp b/e_asin.cpp
|
||||
--- a/e_asin.cpp
|
||||
+++ b/e_asin.cpp
|
||||
@@ -105,12 +105,8 @@ __ieee754_asin(double x)
|
||||
c = (t-w*w)/(s+w);
|
||||
r = p/q;
|
||||
@ -46,9 +46,9 @@ diff --git a/modules/fdlibm/src/e_asin.cpp b/modules/fdlibm/src/e_asin.cpp
|
||||
-#if LDBL_MANT_DIG == 53
|
||||
-__weak_reference(asin, asinl);
|
||||
-#endif
|
||||
diff --git a/modules/fdlibm/src/e_atan2.cpp b/modules/fdlibm/src/e_atan2.cpp
|
||||
--- a/modules/fdlibm/src/e_atan2.cpp
|
||||
+++ b/modules/fdlibm/src/e_atan2.cpp
|
||||
diff --git a/e_atan2.cpp b/e_atan2.cpp
|
||||
--- a/e_atan2.cpp
|
||||
+++ b/e_atan2.cpp
|
||||
@@ -117,12 +117,8 @@ __ieee754_atan2(double y, double x)
|
||||
switch (m) {
|
||||
case 0: return z ; /* atan(+,+) */
|
||||
@ -62,9 +62,9 @@ diff --git a/modules/fdlibm/src/e_atan2.cpp b/modules/fdlibm/src/e_atan2.cpp
|
||||
-#if LDBL_MANT_DIG == 53
|
||||
-__weak_reference(atan2, atan2l);
|
||||
-#endif
|
||||
diff --git a/modules/fdlibm/src/e_atanh.cpp b/modules/fdlibm/src/e_atanh.cpp
|
||||
--- a/modules/fdlibm/src/e_atanh.cpp
|
||||
+++ b/modules/fdlibm/src/e_atanh.cpp
|
||||
diff --git a/e_atanh.cpp b/e_atanh.cpp
|
||||
--- a/e_atanh.cpp
|
||||
+++ b/e_atanh.cpp
|
||||
@@ -56,12 +56,8 @@ __ieee754_atanh(double x)
|
||||
SET_HIGH_WORD(x,ix);
|
||||
if(ix<0x3fe00000) { /* x < 0.5 */
|
||||
@ -78,9 +78,9 @@ diff --git a/modules/fdlibm/src/e_atanh.cpp b/modules/fdlibm/src/e_atanh.cpp
|
||||
-#if LDBL_MANT_DIG == 53
|
||||
-__weak_reference(atanh, atanhl);
|
||||
-#endif
|
||||
diff --git a/modules/fdlibm/src/e_cosh.cpp b/modules/fdlibm/src/e_cosh.cpp
|
||||
--- a/modules/fdlibm/src/e_cosh.cpp
|
||||
+++ b/modules/fdlibm/src/e_cosh.cpp
|
||||
diff --git a/e_cosh.cpp b/e_cosh.cpp
|
||||
--- a/e_cosh.cpp
|
||||
+++ b/e_cosh.cpp
|
||||
@@ -73,12 +73,8 @@ __ieee754_cosh(double x)
|
||||
|
||||
/* |x| in [log(maxdouble), overflowthresold] */
|
||||
@ -94,9 +94,9 @@ diff --git a/modules/fdlibm/src/e_cosh.cpp b/modules/fdlibm/src/e_cosh.cpp
|
||||
-#if (LDBL_MANT_DIG == 53)
|
||||
-__weak_reference(cosh, coshl);
|
||||
-#endif
|
||||
diff --git a/modules/fdlibm/src/e_exp.cpp b/modules/fdlibm/src/e_exp.cpp
|
||||
--- a/modules/fdlibm/src/e_exp.cpp
|
||||
+++ b/modules/fdlibm/src/e_exp.cpp
|
||||
diff --git a/e_exp.cpp b/e_exp.cpp
|
||||
--- a/e_exp.cpp
|
||||
+++ b/e_exp.cpp
|
||||
@@ -152,12 +152,8 @@ __ieee754_exp(double x) /* default IEEE double exp */
|
||||
else y = one-((lo-(x*c)/(2.0-c))-hi);
|
||||
if(k >= -1021) {
|
||||
@ -110,9 +110,9 @@ diff --git a/modules/fdlibm/src/e_exp.cpp b/modules/fdlibm/src/e_exp.cpp
|
||||
-#if (LDBL_MANT_DIG == 53)
|
||||
-__weak_reference(exp, expl);
|
||||
-#endif
|
||||
diff --git a/modules/fdlibm/src/e_hypot.cpp b/modules/fdlibm/src/e_hypot.cpp
|
||||
--- a/modules/fdlibm/src/e_hypot.cpp
|
||||
+++ b/modules/fdlibm/src/e_hypot.cpp
|
||||
diff --git a/e_hypot.cpp b/e_hypot.cpp
|
||||
--- a/e_hypot.cpp
|
||||
+++ b/e_hypot.cpp
|
||||
@@ -117,12 +117,8 @@ __ieee754_hypot(double x, double y)
|
||||
w = sqrt(t1*y1-(w*(-w)-(t1*y2+t2*b)));
|
||||
}
|
||||
@ -126,9 +126,9 @@ diff --git a/modules/fdlibm/src/e_hypot.cpp b/modules/fdlibm/src/e_hypot.cpp
|
||||
-#if LDBL_MANT_DIG == 53
|
||||
-__weak_reference(hypot, hypotl);
|
||||
-#endif
|
||||
diff --git a/modules/fdlibm/src/e_log.cpp b/modules/fdlibm/src/e_log.cpp
|
||||
--- a/modules/fdlibm/src/e_log.cpp
|
||||
+++ b/modules/fdlibm/src/e_log.cpp
|
||||
diff --git a/e_log.cpp b/e_log.cpp
|
||||
--- a/e_log.cpp
|
||||
+++ b/e_log.cpp
|
||||
@@ -135,12 +135,8 @@ __ieee754_log(double x)
|
||||
hfsq=0.5*f*f;
|
||||
if(k==0) return f-(hfsq-s*(hfsq+R)); else
|
||||
@ -142,9 +142,9 @@ diff --git a/modules/fdlibm/src/e_log.cpp b/modules/fdlibm/src/e_log.cpp
|
||||
-#if (LDBL_MANT_DIG == 53)
|
||||
-__weak_reference(log, logl);
|
||||
-#endif
|
||||
diff --git a/modules/fdlibm/src/e_log10.cpp b/modules/fdlibm/src/e_log10.cpp
|
||||
--- a/modules/fdlibm/src/e_log10.cpp
|
||||
+++ b/modules/fdlibm/src/e_log10.cpp
|
||||
diff --git a/e_log10.cpp b/e_log10.cpp
|
||||
--- a/e_log10.cpp
|
||||
+++ b/e_log10.cpp
|
||||
@@ -82,12 +82,8 @@ __ieee754_log10(double x)
|
||||
* with some parallelism and it reduces the error for many args.
|
||||
*/
|
||||
@ -158,9 +158,9 @@ diff --git a/modules/fdlibm/src/e_log10.cpp b/modules/fdlibm/src/e_log10.cpp
|
||||
-#if (LDBL_MANT_DIG == 53)
|
||||
-__weak_reference(log10, log10l);
|
||||
-#endif
|
||||
diff --git a/modules/fdlibm/src/e_log2.cpp b/modules/fdlibm/src/e_log2.cpp
|
||||
--- a/modules/fdlibm/src/e_log2.cpp
|
||||
+++ b/modules/fdlibm/src/e_log2.cpp
|
||||
diff --git a/e_log2.cpp b/e_log2.cpp
|
||||
--- a/e_log2.cpp
|
||||
+++ b/e_log2.cpp
|
||||
@@ -105,12 +105,8 @@ __ieee754_log2(double x)
|
||||
|
||||
/* spadd(val_hi, val_lo, y), except for not using double_t: */
|
||||
@ -174,9 +174,9 @@ diff --git a/modules/fdlibm/src/e_log2.cpp b/modules/fdlibm/src/e_log2.cpp
|
||||
-#if (LDBL_MANT_DIG == 53)
|
||||
-__weak_reference(log2, log2l);
|
||||
-#endif
|
||||
diff --git a/modules/fdlibm/src/e_pow.cpp b/modules/fdlibm/src/e_pow.cpp
|
||||
--- a/modules/fdlibm/src/e_pow.cpp
|
||||
+++ b/modules/fdlibm/src/e_pow.cpp
|
||||
diff --git a/e_pow.cpp b/e_pow.cpp
|
||||
--- a/e_pow.cpp
|
||||
+++ b/e_pow.cpp
|
||||
@@ -302,12 +302,8 @@ __ieee754_pow(double x, double y)
|
||||
r = (z*t1)/(t1-two)-(w+z*w);
|
||||
z = one-(r-z);
|
||||
@ -190,9 +190,9 @@ diff --git a/modules/fdlibm/src/e_pow.cpp b/modules/fdlibm/src/e_pow.cpp
|
||||
-#if (LDBL_MANT_DIG == 53)
|
||||
-__weak_reference(pow, powl);
|
||||
-#endif
|
||||
diff --git a/modules/fdlibm/src/e_sinh.cpp b/modules/fdlibm/src/e_sinh.cpp
|
||||
--- a/modules/fdlibm/src/e_sinh.cpp
|
||||
+++ b/modules/fdlibm/src/e_sinh.cpp
|
||||
diff --git a/e_sinh.cpp b/e_sinh.cpp
|
||||
--- a/e_sinh.cpp
|
||||
+++ b/e_sinh.cpp
|
||||
@@ -67,12 +67,8 @@ __ieee754_sinh(double x)
|
||||
|
||||
/* |x| in [log(maxdouble), overflowthresold] */
|
||||
@ -206,9 +206,9 @@ diff --git a/modules/fdlibm/src/e_sinh.cpp b/modules/fdlibm/src/e_sinh.cpp
|
||||
-#if (LDBL_MANT_DIG == 53)
|
||||
-__weak_reference(sinh, sinhl);
|
||||
-#endif
|
||||
diff --git a/modules/fdlibm/src/s_asinh.cpp b/modules/fdlibm/src/s_asinh.cpp
|
||||
--- a/modules/fdlibm/src/s_asinh.cpp
|
||||
+++ b/modules/fdlibm/src/s_asinh.cpp
|
||||
diff --git a/s_asinh.cpp b/s_asinh.cpp
|
||||
--- a/s_asinh.cpp
|
||||
+++ b/s_asinh.cpp
|
||||
@@ -50,12 +50,8 @@ asinh(double x)
|
||||
t = fabs(x);
|
||||
w = __ieee754_log(2.0*t+one/(__ieee754_sqrt(x*x+one)+t));
|
||||
@ -222,9 +222,9 @@ diff --git a/modules/fdlibm/src/s_asinh.cpp b/modules/fdlibm/src/s_asinh.cpp
|
||||
-#if LDBL_MANT_DIG == 53
|
||||
-__weak_reference(asinh, asinhl);
|
||||
-#endif
|
||||
diff --git a/modules/fdlibm/src/s_atan.cpp b/modules/fdlibm/src/s_atan.cpp
|
||||
--- a/modules/fdlibm/src/s_atan.cpp
|
||||
+++ b/modules/fdlibm/src/s_atan.cpp
|
||||
diff --git a/s_atan.cpp b/s_atan.cpp
|
||||
--- a/s_atan.cpp
|
||||
+++ b/s_atan.cpp
|
||||
@@ -112,12 +112,8 @@ atan(double x)
|
||||
s1 = z*(aT[0]+w*(aT[2]+w*(aT[4]+w*(aT[6]+w*(aT[8]+w*aT[10])))));
|
||||
s2 = w*(aT[1]+w*(aT[3]+w*(aT[5]+w*(aT[7]+w*aT[9]))));
|
||||
@ -238,9 +238,9 @@ diff --git a/modules/fdlibm/src/s_atan.cpp b/modules/fdlibm/src/s_atan.cpp
|
||||
-#if LDBL_MANT_DIG == 53
|
||||
-__weak_reference(atan, atanl);
|
||||
-#endif
|
||||
diff --git a/modules/fdlibm/src/s_cbrt.cpp b/modules/fdlibm/src/s_cbrt.cpp
|
||||
--- a/modules/fdlibm/src/s_cbrt.cpp
|
||||
+++ b/modules/fdlibm/src/s_cbrt.cpp
|
||||
diff --git a/s_cbrt.cpp b/s_cbrt.cpp
|
||||
--- a/s_cbrt.cpp
|
||||
+++ b/s_cbrt.cpp
|
||||
@@ -106,12 +106,8 @@ cbrt(double x)
|
||||
s=t*t; /* t*t is exact */
|
||||
r=x/s; /* error <= 0.5 ulps; |r| < |t| */
|
||||
@ -254,9 +254,9 @@ diff --git a/modules/fdlibm/src/s_cbrt.cpp b/modules/fdlibm/src/s_cbrt.cpp
|
||||
-#if (LDBL_MANT_DIG == 53)
|
||||
-__weak_reference(cbrt, cbrtl);
|
||||
-#endif
|
||||
diff --git a/modules/fdlibm/src/s_ceil.cpp b/modules/fdlibm/src/s_ceil.cpp
|
||||
--- a/modules/fdlibm/src/s_ceil.cpp
|
||||
+++ b/modules/fdlibm/src/s_ceil.cpp
|
||||
diff --git a/s_ceil.cpp b/s_ceil.cpp
|
||||
--- a/s_ceil.cpp
|
||||
+++ b/s_ceil.cpp
|
||||
@@ -65,12 +65,8 @@ ceil(double x)
|
||||
}
|
||||
}
|
||||
@ -270,9 +270,9 @@ diff --git a/modules/fdlibm/src/s_ceil.cpp b/modules/fdlibm/src/s_ceil.cpp
|
||||
-#if LDBL_MANT_DIG == 53
|
||||
-__weak_reference(ceil, ceill);
|
||||
-#endif
|
||||
diff --git a/modules/fdlibm/src/s_cos.cpp b/modules/fdlibm/src/s_cos.cpp
|
||||
--- a/modules/fdlibm/src/s_cos.cpp
|
||||
+++ b/modules/fdlibm/src/s_cos.cpp
|
||||
diff --git a/s_cos.cpp b/s_cos.cpp
|
||||
--- a/s_cos.cpp
|
||||
+++ b/s_cos.cpp
|
||||
@@ -77,12 +77,8 @@ cos(double x)
|
||||
case 0: return __kernel_cos(y[0],y[1]);
|
||||
case 1: return -__kernel_sin(y[0],y[1],1);
|
||||
@ -286,9 +286,9 @@ diff --git a/modules/fdlibm/src/s_cos.cpp b/modules/fdlibm/src/s_cos.cpp
|
||||
-#if (LDBL_MANT_DIG == 53)
|
||||
-__weak_reference(cos, cosl);
|
||||
-#endif
|
||||
diff --git a/modules/fdlibm/src/s_expm1.cpp b/modules/fdlibm/src/s_expm1.cpp
|
||||
--- a/modules/fdlibm/src/s_expm1.cpp
|
||||
+++ b/modules/fdlibm/src/s_expm1.cpp
|
||||
diff --git a/s_expm1.cpp b/s_expm1.cpp
|
||||
--- a/s_expm1.cpp
|
||||
+++ b/s_expm1.cpp
|
||||
@@ -210,12 +210,8 @@ expm1(double x)
|
||||
SET_HIGH_WORD(t,((0x3ff-k)<<20)); /* 2^-k */
|
||||
y = x-(e+t);
|
||||
@ -302,9 +302,9 @@ diff --git a/modules/fdlibm/src/s_expm1.cpp b/modules/fdlibm/src/s_expm1.cpp
|
||||
-#if (LDBL_MANT_DIG == 53)
|
||||
-__weak_reference(expm1, expm1l);
|
||||
-#endif
|
||||
diff --git a/modules/fdlibm/src/s_floor.cpp b/modules/fdlibm/src/s_floor.cpp
|
||||
--- a/modules/fdlibm/src/s_floor.cpp
|
||||
+++ b/modules/fdlibm/src/s_floor.cpp
|
||||
diff --git a/s_floor.cpp b/s_floor.cpp
|
||||
--- a/s_floor.cpp
|
||||
+++ b/s_floor.cpp
|
||||
@@ -66,12 +66,8 @@ floor(double x)
|
||||
}
|
||||
}
|
||||
@ -318,9 +318,9 @@ diff --git a/modules/fdlibm/src/s_floor.cpp b/modules/fdlibm/src/s_floor.cpp
|
||||
-#if LDBL_MANT_DIG == 53
|
||||
-__weak_reference(floor, floorl);
|
||||
-#endif
|
||||
diff --git a/modules/fdlibm/src/s_log1p.cpp b/modules/fdlibm/src/s_log1p.cpp
|
||||
--- a/modules/fdlibm/src/s_log1p.cpp
|
||||
+++ b/modules/fdlibm/src/s_log1p.cpp
|
||||
diff --git a/s_log1p.cpp b/s_log1p.cpp
|
||||
--- a/s_log1p.cpp
|
||||
+++ b/s_log1p.cpp
|
||||
@@ -168,12 +168,8 @@ log1p(double x)
|
||||
return k*ln2_hi-((R-(k*ln2_lo+c))-f);
|
||||
}
|
||||
@ -334,9 +334,9 @@ diff --git a/modules/fdlibm/src/s_log1p.cpp b/modules/fdlibm/src/s_log1p.cpp
|
||||
-#if (LDBL_MANT_DIG == 53)
|
||||
-__weak_reference(log1p, log1pl);
|
||||
-#endif
|
||||
diff --git a/modules/fdlibm/src/s_rint.cpp b/modules/fdlibm/src/s_rint.cpp
|
||||
--- a/modules/fdlibm/src/s_rint.cpp
|
||||
+++ b/modules/fdlibm/src/s_rint.cpp
|
||||
diff --git a/s_rint.cpp b/s_rint.cpp
|
||||
--- a/s_rint.cpp
|
||||
+++ b/s_rint.cpp
|
||||
@@ -80,12 +80,8 @@ rint(double x)
|
||||
if((i1&i)==0) return x; /* x is integral */
|
||||
i>>=1;
|
||||
@ -350,9 +350,9 @@ diff --git a/modules/fdlibm/src/s_rint.cpp b/modules/fdlibm/src/s_rint.cpp
|
||||
-#if (LDBL_MANT_DIG == 53)
|
||||
-__weak_reference(rint, rintl);
|
||||
-#endif
|
||||
diff --git a/modules/fdlibm/src/s_scalbn.cpp b/modules/fdlibm/src/s_scalbn.cpp
|
||||
--- a/modules/fdlibm/src/s_scalbn.cpp
|
||||
+++ b/modules/fdlibm/src/s_scalbn.cpp
|
||||
diff --git a/s_scalbn.cpp b/s_scalbn.cpp
|
||||
--- a/s_scalbn.cpp
|
||||
+++ b/s_scalbn.cpp
|
||||
@@ -36,13 +36,8 @@ double scalbn(double x, int n)
|
||||
if (n < -1022)
|
||||
n = -1022;
|
||||
@ -367,9 +367,9 @@ diff --git a/modules/fdlibm/src/s_scalbn.cpp b/modules/fdlibm/src/s_scalbn.cpp
|
||||
-__weak_reference(scalbn, ldexpl);
|
||||
-__weak_reference(scalbn, scalbnl);
|
||||
-#endif
|
||||
diff --git a/modules/fdlibm/src/s_sin.cpp b/modules/fdlibm/src/s_sin.cpp
|
||||
--- a/modules/fdlibm/src/s_sin.cpp
|
||||
+++ b/modules/fdlibm/src/s_sin.cpp
|
||||
diff --git a/s_sin.cpp b/s_sin.cpp
|
||||
--- a/s_sin.cpp
|
||||
+++ b/s_sin.cpp
|
||||
@@ -77,12 +77,8 @@ sin(double x)
|
||||
case 0: return __kernel_sin(y[0],y[1],1);
|
||||
case 1: return __kernel_cos(y[0],y[1]);
|
||||
@ -383,9 +383,9 @@ diff --git a/modules/fdlibm/src/s_sin.cpp b/modules/fdlibm/src/s_sin.cpp
|
||||
-#if (LDBL_MANT_DIG == 53)
|
||||
-__weak_reference(sin, sinl);
|
||||
-#endif
|
||||
diff --git a/modules/fdlibm/src/s_tan.cpp b/modules/fdlibm/src/s_tan.cpp
|
||||
--- a/modules/fdlibm/src/s_tan.cpp
|
||||
+++ b/modules/fdlibm/src/s_tan.cpp
|
||||
diff --git a/s_tan.cpp b/s_tan.cpp
|
||||
--- a/s_tan.cpp
|
||||
+++ b/s_tan.cpp
|
||||
@@ -71,12 +71,8 @@ tan(double x)
|
||||
|
||||
/* argument reduction needed */
|
||||
@ -399,9 +399,9 @@ diff --git a/modules/fdlibm/src/s_tan.cpp b/modules/fdlibm/src/s_tan.cpp
|
||||
-#if (LDBL_MANT_DIG == 53)
|
||||
-__weak_reference(tan, tanl);
|
||||
-#endif
|
||||
diff --git a/modules/fdlibm/src/s_tanh.cpp b/modules/fdlibm/src/s_tanh.cpp
|
||||
--- a/modules/fdlibm/src/s_tanh.cpp
|
||||
+++ b/modules/fdlibm/src/s_tanh.cpp
|
||||
diff --git a/s_tanh.cpp b/s_tanh.cpp
|
||||
--- a/s_tanh.cpp
|
||||
+++ b/s_tanh.cpp
|
||||
@@ -72,12 +72,8 @@ tanh(double x)
|
||||
z= -t/(t+two);
|
||||
}
|
||||
@ -415,9 +415,9 @@ diff --git a/modules/fdlibm/src/s_tanh.cpp b/modules/fdlibm/src/s_tanh.cpp
|
||||
-#if (LDBL_MANT_DIG == 53)
|
||||
-__weak_reference(tanh, tanhl);
|
||||
-#endif
|
||||
diff --git a/modules/fdlibm/src/s_trunc.cpp b/modules/fdlibm/src/s_trunc.cpp
|
||||
--- a/modules/fdlibm/src/s_trunc.cpp
|
||||
+++ b/modules/fdlibm/src/s_trunc.cpp
|
||||
diff --git a/s_trunc.cpp b/s_trunc.cpp
|
||||
--- a/s_trunc.cpp
|
||||
+++ b/s_trunc.cpp
|
||||
@@ -55,12 +55,8 @@ trunc(double x)
|
||||
i = ((u_int32_t)(0xffffffff))>>(j0-20);
|
||||
if((i1&i)==0) return x; /* x is integral */
|
||||
|
@ -1,6 +1,6 @@
|
||||
diff --git a/modules/fdlibm/src/e_acos.cpp b/modules/fdlibm/src/e_acos.cpp
|
||||
--- a/modules/fdlibm/src/e_acos.cpp
|
||||
+++ b/modules/fdlibm/src/e_acos.cpp
|
||||
diff --git a/e_acos.cpp b/e_acos.cpp
|
||||
--- a/e_acos.cpp
|
||||
+++ b/e_acos.cpp
|
||||
@@ -6,18 +6,18 @@
|
||||
*
|
||||
* Developed at SunSoft, a Sun Microsystems, Inc. business.
|
||||
@ -22,9 +22,9 @@ diff --git a/modules/fdlibm/src/e_acos.cpp b/modules/fdlibm/src/e_acos.cpp
|
||||
* For |x|<=0.5
|
||||
* acos(x) = pi/2 - (x + x*x^2*R(x^2)) (see asin.c)
|
||||
* For x>0.5
|
||||
diff --git a/modules/fdlibm/src/e_acosh.cpp b/modules/fdlibm/src/e_acosh.cpp
|
||||
--- a/modules/fdlibm/src/e_acosh.cpp
|
||||
+++ b/modules/fdlibm/src/e_acosh.cpp
|
||||
diff --git a/e_acosh.cpp b/e_acosh.cpp
|
||||
--- a/e_acosh.cpp
|
||||
+++ b/e_acosh.cpp
|
||||
@@ -7,18 +7,18 @@
|
||||
* Developed at SunSoft, a Sun Microsystems, Inc. business.
|
||||
* Permission to use, copy, modify, and distribute this
|
||||
@ -46,9 +46,9 @@ diff --git a/modules/fdlibm/src/e_acosh.cpp b/modules/fdlibm/src/e_acosh.cpp
|
||||
* we have
|
||||
* acosh(x) := log(x)+ln2, if x is large; else
|
||||
* acosh(x) := log(2x-1/(sqrt(x*x-1)+x)) if x>2; else
|
||||
diff --git a/modules/fdlibm/src/e_asin.cpp b/modules/fdlibm/src/e_asin.cpp
|
||||
--- a/modules/fdlibm/src/e_asin.cpp
|
||||
+++ b/modules/fdlibm/src/e_asin.cpp
|
||||
diff --git a/e_asin.cpp b/e_asin.cpp
|
||||
--- a/e_asin.cpp
|
||||
+++ b/e_asin.cpp
|
||||
@@ -6,18 +6,18 @@
|
||||
*
|
||||
* Developed at SunSoft, a Sun Microsystems, Inc. business.
|
||||
@ -70,9 +70,9 @@ diff --git a/modules/fdlibm/src/e_asin.cpp b/modules/fdlibm/src/e_asin.cpp
|
||||
* asin(x) = x + x*x^2*R(x^2)
|
||||
* where
|
||||
* R(x^2) is a rational approximation of (asin(x)-x)/x^3
|
||||
diff --git a/modules/fdlibm/src/e_atan2.cpp b/modules/fdlibm/src/e_atan2.cpp
|
||||
--- a/modules/fdlibm/src/e_atan2.cpp
|
||||
+++ b/modules/fdlibm/src/e_atan2.cpp
|
||||
diff --git a/e_atan2.cpp b/e_atan2.cpp
|
||||
--- a/e_atan2.cpp
|
||||
+++ b/e_atan2.cpp
|
||||
@@ -7,18 +7,18 @@
|
||||
* Developed at SunSoft, a Sun Microsystems, Inc. business.
|
||||
* Permission to use, copy, modify, and distribute this
|
||||
@ -94,9 +94,9 @@ diff --git a/modules/fdlibm/src/e_atan2.cpp b/modules/fdlibm/src/e_atan2.cpp
|
||||
* ARG (x+iy) = arctan(y/x) ... if x > 0,
|
||||
* ARG (x+iy) = pi - arctan[y/(-x)] ... if x < 0,
|
||||
*
|
||||
diff --git a/modules/fdlibm/src/e_atanh.cpp b/modules/fdlibm/src/e_atanh.cpp
|
||||
--- a/modules/fdlibm/src/e_atanh.cpp
|
||||
+++ b/modules/fdlibm/src/e_atanh.cpp
|
||||
diff --git a/e_atanh.cpp b/e_atanh.cpp
|
||||
--- a/e_atanh.cpp
|
||||
+++ b/e_atanh.cpp
|
||||
@@ -7,18 +7,18 @@
|
||||
* Developed at SunSoft, a Sun Microsystems, Inc. business.
|
||||
* Permission to use, copy, modify, and distribute this
|
||||
@ -118,9 +118,9 @@ diff --git a/modules/fdlibm/src/e_atanh.cpp b/modules/fdlibm/src/e_atanh.cpp
|
||||
* 1 2x x
|
||||
* atanh(x) = --- * log(1 + -------) = 0.5 * log1p(2 * --------)
|
||||
* 2 1 - x 1 - x
|
||||
diff --git a/modules/fdlibm/src/e_cosh.cpp b/modules/fdlibm/src/e_cosh.cpp
|
||||
--- a/modules/fdlibm/src/e_cosh.cpp
|
||||
+++ b/modules/fdlibm/src/e_cosh.cpp
|
||||
diff --git a/e_cosh.cpp b/e_cosh.cpp
|
||||
--- a/e_cosh.cpp
|
||||
+++ b/e_cosh.cpp
|
||||
@@ -6,18 +6,18 @@
|
||||
*
|
||||
* Developed at SunSoft, a Sun Microsystems, Inc. business.
|
||||
@ -142,9 +142,9 @@ diff --git a/modules/fdlibm/src/e_cosh.cpp b/modules/fdlibm/src/e_cosh.cpp
|
||||
* 2.
|
||||
* [ exp(x) - 1 ]^2
|
||||
* 0 <= x <= ln2/2 : cosh(x) := 1 + -------------------
|
||||
diff --git a/modules/fdlibm/src/e_exp.cpp b/modules/fdlibm/src/e_exp.cpp
|
||||
--- a/modules/fdlibm/src/e_exp.cpp
|
||||
+++ b/modules/fdlibm/src/e_exp.cpp
|
||||
diff --git a/e_exp.cpp b/e_exp.cpp
|
||||
--- a/e_exp.cpp
|
||||
+++ b/e_exp.cpp
|
||||
@@ -5,18 +5,18 @@
|
||||
* Copyright (C) 2004 by Sun Microsystems, Inc. All rights reserved.
|
||||
*
|
||||
@ -166,9 +166,9 @@ diff --git a/modules/fdlibm/src/e_exp.cpp b/modules/fdlibm/src/e_exp.cpp
|
||||
* 1. Argument reduction:
|
||||
* Reduce x to an r so that |r| <= 0.5*ln2 ~ 0.34658.
|
||||
* Given x, find r and integer k such that
|
||||
diff --git a/modules/fdlibm/src/e_hypot.cpp b/modules/fdlibm/src/e_hypot.cpp
|
||||
--- a/modules/fdlibm/src/e_hypot.cpp
|
||||
+++ b/modules/fdlibm/src/e_hypot.cpp
|
||||
diff --git a/e_hypot.cpp b/e_hypot.cpp
|
||||
--- a/e_hypot.cpp
|
||||
+++ b/e_hypot.cpp
|
||||
@@ -6,18 +6,18 @@
|
||||
*
|
||||
* Developed at SunSoft, a Sun Microsystems, Inc. business.
|
||||
@ -190,9 +190,9 @@ diff --git a/modules/fdlibm/src/e_hypot.cpp b/modules/fdlibm/src/e_hypot.cpp
|
||||
* has error less than sqrt(2)/2 ulp, than
|
||||
* sqrt(z) has error less than 1 ulp (exercise).
|
||||
*
|
||||
diff --git a/modules/fdlibm/src/e_log.cpp b/modules/fdlibm/src/e_log.cpp
|
||||
--- a/modules/fdlibm/src/e_log.cpp
|
||||
+++ b/modules/fdlibm/src/e_log.cpp
|
||||
diff --git a/e_log.cpp b/e_log.cpp
|
||||
--- a/e_log.cpp
|
||||
+++ b/e_log.cpp
|
||||
@@ -6,18 +6,18 @@
|
||||
*
|
||||
* Developed at SunSoft, a Sun Microsystems, Inc. business.
|
||||
@ -214,9 +214,9 @@ diff --git a/modules/fdlibm/src/e_log.cpp b/modules/fdlibm/src/e_log.cpp
|
||||
* 1. Argument Reduction: find k and f such that
|
||||
* x = 2^k * (1+f),
|
||||
* where sqrt(2)/2 < 1+f < sqrt(2) .
|
||||
diff --git a/modules/fdlibm/src/e_log10.cpp b/modules/fdlibm/src/e_log10.cpp
|
||||
--- a/modules/fdlibm/src/e_log10.cpp
|
||||
+++ b/modules/fdlibm/src/e_log10.cpp
|
||||
diff --git a/e_log10.cpp b/e_log10.cpp
|
||||
--- a/e_log10.cpp
|
||||
+++ b/e_log10.cpp
|
||||
@@ -6,32 +6,32 @@
|
||||
*
|
||||
* Developed at SunSoft, a Sun Microsystems, Inc. business.
|
||||
@ -253,9 +253,9 @@ diff --git a/modules/fdlibm/src/e_log10.cpp b/modules/fdlibm/src/e_log10.cpp
|
||||
log10_2lo = 3.69423907715893078616e-13; /* 0x3D59FEF3, 0x11F12B36 */
|
||||
|
||||
static const double zero = 0.0;
|
||||
diff --git a/modules/fdlibm/src/e_log2.cpp b/modules/fdlibm/src/e_log2.cpp
|
||||
--- a/modules/fdlibm/src/e_log2.cpp
|
||||
+++ b/modules/fdlibm/src/e_log2.cpp
|
||||
diff --git a/e_log2.cpp b/e_log2.cpp
|
||||
--- a/e_log2.cpp
|
||||
+++ b/e_log2.cpp
|
||||
@@ -6,18 +6,18 @@
|
||||
*
|
||||
* Developed at SunSoft, a Sun Microsystems, Inc. business.
|
||||
@ -277,9 +277,9 @@ diff --git a/modules/fdlibm/src/e_log2.cpp b/modules/fdlibm/src/e_log2.cpp
|
||||
* This reduces x to {k, 1+f} exactly as in e_log.c, then calls the kernel,
|
||||
* then does the combining and scaling steps
|
||||
* log2(x) = (f - 0.5*f*f + k_log1p(f)) / ln2 + k
|
||||
diff --git a/modules/fdlibm/src/e_pow.cpp b/modules/fdlibm/src/e_pow.cpp
|
||||
--- a/modules/fdlibm/src/e_pow.cpp
|
||||
+++ b/modules/fdlibm/src/e_pow.cpp
|
||||
diff --git a/e_pow.cpp b/e_pow.cpp
|
||||
--- a/e_pow.cpp
|
||||
+++ b/e_pow.cpp
|
||||
@@ -4,18 +4,18 @@
|
||||
* Copyright (C) 2004 by Sun Microsystems, Inc. All rights reserved.
|
||||
*
|
||||
@ -301,9 +301,9 @@ diff --git a/modules/fdlibm/src/e_pow.cpp b/modules/fdlibm/src/e_pow.cpp
|
||||
* 1. Compute and return log2(x) in two pieces:
|
||||
* log2(x) = w1 + w2,
|
||||
* where w1 has 53-24 = 29 bit trailing zeros.
|
||||
diff --git a/modules/fdlibm/src/e_rem_pio2.cpp b/modules/fdlibm/src/e_rem_pio2.cpp
|
||||
--- a/modules/fdlibm/src/e_rem_pio2.cpp
|
||||
+++ b/modules/fdlibm/src/e_rem_pio2.cpp
|
||||
diff --git a/e_rem_pio2.cpp b/e_rem_pio2.cpp
|
||||
--- a/e_rem_pio2.cpp
|
||||
+++ b/e_rem_pio2.cpp
|
||||
@@ -8,18 +8,18 @@
|
||||
* Permission to use, copy, modify, and distribute this
|
||||
* software is freely granted, provided that this notice
|
||||
@ -325,9 +325,9 @@ diff --git a/modules/fdlibm/src/e_rem_pio2.cpp b/modules/fdlibm/src/e_rem_pio2.c
|
||||
*/
|
||||
|
||||
#include <float.h>
|
||||
diff --git a/modules/fdlibm/src/e_sinh.cpp b/modules/fdlibm/src/e_sinh.cpp
|
||||
--- a/modules/fdlibm/src/e_sinh.cpp
|
||||
+++ b/modules/fdlibm/src/e_sinh.cpp
|
||||
diff --git a/e_sinh.cpp b/e_sinh.cpp
|
||||
--- a/e_sinh.cpp
|
||||
+++ b/e_sinh.cpp
|
||||
@@ -6,18 +6,18 @@
|
||||
*
|
||||
* Developed at SunSoft, a Sun Microsystems, Inc. business.
|
||||
@ -349,9 +349,9 @@ diff --git a/modules/fdlibm/src/e_sinh.cpp b/modules/fdlibm/src/e_sinh.cpp
|
||||
* 2.
|
||||
* E + E/(E+1)
|
||||
* 0 <= x <= 22 : sinh(x) := --------------, E=expm1(x)
|
||||
diff --git a/modules/fdlibm/src/k_cos.cpp b/modules/fdlibm/src/k_cos.cpp
|
||||
--- a/modules/fdlibm/src/k_cos.cpp
|
||||
+++ b/modules/fdlibm/src/k_cos.cpp
|
||||
diff --git a/k_cos.cpp b/k_cos.cpp
|
||||
--- a/k_cos.cpp
|
||||
+++ b/k_cos.cpp
|
||||
@@ -6,18 +6,18 @@
|
||||
*
|
||||
* Developed at SunSoft, a Sun Microsystems, Inc. business.
|
||||
@ -373,9 +373,9 @@ diff --git a/modules/fdlibm/src/k_cos.cpp b/modules/fdlibm/src/k_cos.cpp
|
||||
* Input y is the tail of x.
|
||||
*
|
||||
* Algorithm
|
||||
diff --git a/modules/fdlibm/src/k_exp.cpp b/modules/fdlibm/src/k_exp.cpp
|
||||
--- a/modules/fdlibm/src/k_exp.cpp
|
||||
+++ b/modules/fdlibm/src/k_exp.cpp
|
||||
diff --git a/k_exp.cpp b/k_exp.cpp
|
||||
--- a/k_exp.cpp
|
||||
+++ b/k_exp.cpp
|
||||
@@ -21,22 +21,22 @@
|
||||
* DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
|
||||
* OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
|
||||
@ -402,9 +402,9 @@ diff --git a/modules/fdlibm/src/k_exp.cpp b/modules/fdlibm/src/k_exp.cpp
|
||||
* Compute exp(x), scaled to avoid spurious overflow. An exponent is
|
||||
* returned separately in 'expt'.
|
||||
*
|
||||
diff --git a/modules/fdlibm/src/k_log.h b/modules/fdlibm/src/k_log.h
|
||||
--- a/modules/fdlibm/src/k_log.h
|
||||
+++ b/modules/fdlibm/src/k_log.h
|
||||
diff --git a/k_log.h b/k_log.h
|
||||
--- a/k_log.h
|
||||
+++ b/k_log.h
|
||||
@@ -6,18 +6,18 @@
|
||||
*
|
||||
* Developed at SunSoft, a Sun Microsystems, Inc. business.
|
||||
@ -426,9 +426,9 @@ diff --git a/modules/fdlibm/src/k_log.h b/modules/fdlibm/src/k_log.h
|
||||
* The following describes the overall strategy for computing
|
||||
* logarithms in base e. The argument reduction and adding the final
|
||||
* term of the polynomial are done by the caller for increased accuracy
|
||||
diff --git a/modules/fdlibm/src/k_rem_pio2.cpp b/modules/fdlibm/src/k_rem_pio2.cpp
|
||||
--- a/modules/fdlibm/src/k_rem_pio2.cpp
|
||||
+++ b/modules/fdlibm/src/k_rem_pio2.cpp
|
||||
diff --git a/k_rem_pio2.cpp b/k_rem_pio2.cpp
|
||||
--- a/k_rem_pio2.cpp
|
||||
+++ b/k_rem_pio2.cpp
|
||||
@@ -6,18 +6,18 @@
|
||||
*
|
||||
* Developed at SunSoft, a Sun Microsystems, Inc. business.
|
||||
@ -450,9 +450,9 @@ diff --git a/modules/fdlibm/src/k_rem_pio2.cpp b/modules/fdlibm/src/k_rem_pio2.c
|
||||
* __kernel_rem_pio2 return the last three digits of N with
|
||||
* y = x - N*pi/2
|
||||
* so that |y| < pi/2.
|
||||
diff --git a/modules/fdlibm/src/k_sin.cpp b/modules/fdlibm/src/k_sin.cpp
|
||||
--- a/modules/fdlibm/src/k_sin.cpp
|
||||
+++ b/modules/fdlibm/src/k_sin.cpp
|
||||
diff --git a/k_sin.cpp b/k_sin.cpp
|
||||
--- a/k_sin.cpp
|
||||
+++ b/k_sin.cpp
|
||||
@@ -6,18 +6,18 @@
|
||||
*
|
||||
* Developed at SunSoft, a Sun Microsystems, Inc. business.
|
||||
@ -474,9 +474,9 @@ diff --git a/modules/fdlibm/src/k_sin.cpp b/modules/fdlibm/src/k_sin.cpp
|
||||
* Input iy indicates whether y is 0. (if iy=0, y assume to be 0).
|
||||
*
|
||||
* Algorithm
|
||||
diff --git a/modules/fdlibm/src/k_tan.cpp b/modules/fdlibm/src/k_tan.cpp
|
||||
--- a/modules/fdlibm/src/k_tan.cpp
|
||||
+++ b/modules/fdlibm/src/k_tan.cpp
|
||||
diff --git a/k_tan.cpp b/k_tan.cpp
|
||||
--- a/k_tan.cpp
|
||||
+++ b/k_tan.cpp
|
||||
@@ -6,18 +6,18 @@
|
||||
*
|
||||
* Permission to use, copy, modify, and distribute this
|
||||
@ -498,9 +498,9 @@ diff --git a/modules/fdlibm/src/k_tan.cpp b/modules/fdlibm/src/k_tan.cpp
|
||||
* Input k indicates whether tan (if k = 1) or -1/tan (if k = -1) is returned.
|
||||
*
|
||||
* Algorithm
|
||||
diff --git a/modules/fdlibm/src/s_asinh.cpp b/modules/fdlibm/src/s_asinh.cpp
|
||||
--- a/modules/fdlibm/src/s_asinh.cpp
|
||||
+++ b/modules/fdlibm/src/s_asinh.cpp
|
||||
diff --git a/s_asinh.cpp b/s_asinh.cpp
|
||||
--- a/s_asinh.cpp
|
||||
+++ b/s_asinh.cpp
|
||||
@@ -5,18 +5,18 @@
|
||||
*
|
||||
* Developed at SunPro, a Sun Microsystems, Inc. business.
|
||||
@ -522,9 +522,9 @@ diff --git a/modules/fdlibm/src/s_asinh.cpp b/modules/fdlibm/src/s_asinh.cpp
|
||||
* we have
|
||||
* asinh(x) := x if 1+x*x=1,
|
||||
* := sign(x)*(log(x)+ln2)) for large |x|, else
|
||||
diff --git a/modules/fdlibm/src/s_atan.cpp b/modules/fdlibm/src/s_atan.cpp
|
||||
--- a/modules/fdlibm/src/s_atan.cpp
|
||||
+++ b/modules/fdlibm/src/s_atan.cpp
|
||||
diff --git a/s_atan.cpp b/s_atan.cpp
|
||||
--- a/s_atan.cpp
|
||||
+++ b/s_atan.cpp
|
||||
@@ -5,18 +5,18 @@
|
||||
*
|
||||
* Developed at SunPro, a Sun Microsystems, Inc. business.
|
||||
@ -546,9 +546,9 @@ diff --git a/modules/fdlibm/src/s_atan.cpp b/modules/fdlibm/src/s_atan.cpp
|
||||
* is further reduced to one of the following intervals and the
|
||||
* arctangent of t is evaluated by the corresponding formula:
|
||||
*
|
||||
diff --git a/modules/fdlibm/src/s_cbrt.cpp b/modules/fdlibm/src/s_cbrt.cpp
|
||||
--- a/modules/fdlibm/src/s_cbrt.cpp
|
||||
+++ b/modules/fdlibm/src/s_cbrt.cpp
|
||||
diff --git a/s_cbrt.cpp b/s_cbrt.cpp
|
||||
--- a/s_cbrt.cpp
|
||||
+++ b/s_cbrt.cpp
|
||||
@@ -7,18 +7,18 @@
|
||||
* Permission to use, copy, modify, and distribute this
|
||||
* software is freely granted, provided that this notice
|
||||
@ -570,9 +570,9 @@ diff --git a/modules/fdlibm/src/s_cbrt.cpp b/modules/fdlibm/src/s_cbrt.cpp
|
||||
* Return cube root of x
|
||||
*/
|
||||
static const u_int32_t
|
||||
diff --git a/modules/fdlibm/src/s_ceil.cpp b/modules/fdlibm/src/s_ceil.cpp
|
||||
--- a/modules/fdlibm/src/s_ceil.cpp
|
||||
+++ b/modules/fdlibm/src/s_ceil.cpp
|
||||
diff --git a/s_ceil.cpp b/s_ceil.cpp
|
||||
--- a/s_ceil.cpp
|
||||
+++ b/s_ceil.cpp
|
||||
@@ -5,18 +5,18 @@
|
||||
*
|
||||
* Developed at SunPro, a Sun Microsystems, Inc. business.
|
||||
@ -594,9 +594,9 @@ diff --git a/modules/fdlibm/src/s_ceil.cpp b/modules/fdlibm/src/s_ceil.cpp
|
||||
* Bit twiddling.
|
||||
* Exception:
|
||||
* Inexact flag raised if x not equal to ceil(x).
|
||||
diff --git a/modules/fdlibm/src/s_ceilf.cpp b/modules/fdlibm/src/s_ceilf.cpp
|
||||
--- a/modules/fdlibm/src/s_ceilf.cpp
|
||||
+++ b/modules/fdlibm/src/s_ceilf.cpp
|
||||
diff --git a/s_ceilf.cpp b/s_ceilf.cpp
|
||||
--- a/s_ceilf.cpp
|
||||
+++ b/s_ceilf.cpp
|
||||
@@ -8,18 +8,18 @@
|
||||
*
|
||||
* Developed at SunPro, a Sun Microsystems, Inc. business.
|
||||
@ -618,9 +618,9 @@ diff --git a/modules/fdlibm/src/s_ceilf.cpp b/modules/fdlibm/src/s_ceilf.cpp
|
||||
float
|
||||
ceilf(float x)
|
||||
{
|
||||
diff --git a/modules/fdlibm/src/s_copysign.cpp b/modules/fdlibm/src/s_copysign.cpp
|
||||
--- a/modules/fdlibm/src/s_copysign.cpp
|
||||
+++ b/modules/fdlibm/src/s_copysign.cpp
|
||||
diff --git a/s_copysign.cpp b/s_copysign.cpp
|
||||
--- a/s_copysign.cpp
|
||||
+++ b/s_copysign.cpp
|
||||
@@ -5,18 +5,18 @@
|
||||
*
|
||||
* Developed at SunPro, a Sun Microsystems, Inc. business.
|
||||
@ -642,9 +642,9 @@ diff --git a/modules/fdlibm/src/s_copysign.cpp b/modules/fdlibm/src/s_copysign.c
|
||||
*/
|
||||
|
||||
#include "math_private.h"
|
||||
diff --git a/modules/fdlibm/src/s_cos.cpp b/modules/fdlibm/src/s_cos.cpp
|
||||
--- a/modules/fdlibm/src/s_cos.cpp
|
||||
+++ b/modules/fdlibm/src/s_cos.cpp
|
||||
diff --git a/s_cos.cpp b/s_cos.cpp
|
||||
--- a/s_cos.cpp
|
||||
+++ b/s_cos.cpp
|
||||
@@ -5,18 +5,18 @@
|
||||
*
|
||||
* Developed at SunPro, a Sun Microsystems, Inc. business.
|
||||
@ -666,9 +666,9 @@ diff --git a/modules/fdlibm/src/s_cos.cpp b/modules/fdlibm/src/s_cos.cpp
|
||||
* __kernel_sin ... sine function on [-pi/4,pi/4]
|
||||
* __kernel_cos ... cosine function on [-pi/4,pi/4]
|
||||
* __ieee754_rem_pio2 ... argument reduction routine
|
||||
diff --git a/modules/fdlibm/src/s_expm1.cpp b/modules/fdlibm/src/s_expm1.cpp
|
||||
--- a/modules/fdlibm/src/s_expm1.cpp
|
||||
+++ b/modules/fdlibm/src/s_expm1.cpp
|
||||
diff --git a/s_expm1.cpp b/s_expm1.cpp
|
||||
--- a/s_expm1.cpp
|
||||
+++ b/s_expm1.cpp
|
||||
@@ -5,18 +5,18 @@
|
||||
*
|
||||
* Developed at SunPro, a Sun Microsystems, Inc. business.
|
||||
@ -690,9 +690,9 @@ diff --git a/modules/fdlibm/src/s_expm1.cpp b/modules/fdlibm/src/s_expm1.cpp
|
||||
* 1. Argument reduction:
|
||||
* Given x, find r and integer k such that
|
||||
*
|
||||
diff --git a/modules/fdlibm/src/s_fabs.cpp b/modules/fdlibm/src/s_fabs.cpp
|
||||
--- a/modules/fdlibm/src/s_fabs.cpp
|
||||
+++ b/modules/fdlibm/src/s_fabs.cpp
|
||||
diff --git a/s_fabs.cpp b/s_fabs.cpp
|
||||
--- a/s_fabs.cpp
|
||||
+++ b/s_fabs.cpp
|
||||
@@ -5,18 +5,18 @@
|
||||
*
|
||||
* Developed at SunPro, a Sun Microsystems, Inc. business.
|
||||
@ -714,9 +714,9 @@ diff --git a/modules/fdlibm/src/s_fabs.cpp b/modules/fdlibm/src/s_fabs.cpp
|
||||
#include "math_private.h"
|
||||
|
||||
double
|
||||
diff --git a/modules/fdlibm/src/s_floor.cpp b/modules/fdlibm/src/s_floor.cpp
|
||||
--- a/modules/fdlibm/src/s_floor.cpp
|
||||
+++ b/modules/fdlibm/src/s_floor.cpp
|
||||
diff --git a/s_floor.cpp b/s_floor.cpp
|
||||
--- a/s_floor.cpp
|
||||
+++ b/s_floor.cpp
|
||||
@@ -5,18 +5,18 @@
|
||||
*
|
||||
* Developed at SunPro, a Sun Microsystems, Inc. business.
|
||||
@ -738,9 +738,9 @@ diff --git a/modules/fdlibm/src/s_floor.cpp b/modules/fdlibm/src/s_floor.cpp
|
||||
* Bit twiddling.
|
||||
* Exception:
|
||||
* Inexact flag raised if x not equal to floor(x).
|
||||
diff --git a/modules/fdlibm/src/s_floorf.cpp b/modules/fdlibm/src/s_floorf.cpp
|
||||
--- a/modules/fdlibm/src/s_floorf.cpp
|
||||
+++ b/modules/fdlibm/src/s_floorf.cpp
|
||||
diff --git a/s_floorf.cpp b/s_floorf.cpp
|
||||
--- a/s_floorf.cpp
|
||||
+++ b/s_floorf.cpp
|
||||
@@ -8,18 +8,18 @@
|
||||
*
|
||||
* Developed at SunPro, a Sun Microsystems, Inc. business.
|
||||
@ -762,9 +762,9 @@ diff --git a/modules/fdlibm/src/s_floorf.cpp b/modules/fdlibm/src/s_floorf.cpp
|
||||
* Bit twiddling.
|
||||
* Exception:
|
||||
* Inexact flag raised if x not equal to floorf(x).
|
||||
diff --git a/modules/fdlibm/src/s_log1p.cpp b/modules/fdlibm/src/s_log1p.cpp
|
||||
--- a/modules/fdlibm/src/s_log1p.cpp
|
||||
+++ b/modules/fdlibm/src/s_log1p.cpp
|
||||
diff --git a/s_log1p.cpp b/s_log1p.cpp
|
||||
--- a/s_log1p.cpp
|
||||
+++ b/s_log1p.cpp
|
||||
@@ -5,18 +5,18 @@
|
||||
*
|
||||
* Developed at SunPro, a Sun Microsystems, Inc. business.
|
||||
@ -786,9 +786,9 @@ diff --git a/modules/fdlibm/src/s_log1p.cpp b/modules/fdlibm/src/s_log1p.cpp
|
||||
* 1+x = 2^k * (1+f),
|
||||
* where sqrt(2)/2 < 1+f < sqrt(2) .
|
||||
*
|
||||
diff --git a/modules/fdlibm/src/s_nearbyint.cpp b/modules/fdlibm/src/s_nearbyint.cpp
|
||||
--- a/modules/fdlibm/src/s_nearbyint.cpp
|
||||
+++ b/modules/fdlibm/src/s_nearbyint.cpp
|
||||
diff --git a/s_nearbyint.cpp b/s_nearbyint.cpp
|
||||
--- a/s_nearbyint.cpp
|
||||
+++ b/s_nearbyint.cpp
|
||||
@@ -21,18 +21,18 @@
|
||||
* DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
|
||||
* OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
|
||||
@ -810,9 +810,9 @@ diff --git a/modules/fdlibm/src/s_nearbyint.cpp b/modules/fdlibm/src/s_nearbyint
|
||||
/*
|
||||
* We save and restore the floating-point environment to avoid raising
|
||||
* an inexact exception. We can get away with using fesetenv()
|
||||
diff --git a/modules/fdlibm/src/s_rint.cpp b/modules/fdlibm/src/s_rint.cpp
|
||||
--- a/modules/fdlibm/src/s_rint.cpp
|
||||
+++ b/modules/fdlibm/src/s_rint.cpp
|
||||
diff --git a/s_rint.cpp b/s_rint.cpp
|
||||
--- a/s_rint.cpp
|
||||
+++ b/s_rint.cpp
|
||||
@@ -5,18 +5,18 @@
|
||||
*
|
||||
* Developed at SunPro, a Sun Microsystems, Inc. business.
|
||||
@ -834,9 +834,9 @@ diff --git a/modules/fdlibm/src/s_rint.cpp b/modules/fdlibm/src/s_rint.cpp
|
||||
* Method:
|
||||
* Using floating addition.
|
||||
* Exception:
|
||||
diff --git a/modules/fdlibm/src/s_rintf.cpp b/modules/fdlibm/src/s_rintf.cpp
|
||||
--- a/modules/fdlibm/src/s_rintf.cpp
|
||||
+++ b/modules/fdlibm/src/s_rintf.cpp
|
||||
diff --git a/s_rintf.cpp b/s_rintf.cpp
|
||||
--- a/s_rintf.cpp
|
||||
+++ b/s_rintf.cpp
|
||||
@@ -8,18 +8,18 @@
|
||||
*
|
||||
* Developed at SunPro, a Sun Microsystems, Inc. business.
|
||||
@ -858,9 +858,9 @@ diff --git a/modules/fdlibm/src/s_rintf.cpp b/modules/fdlibm/src/s_rintf.cpp
|
||||
|
||||
static const float
|
||||
TWO23[2]={
|
||||
diff --git a/modules/fdlibm/src/s_sin.cpp b/modules/fdlibm/src/s_sin.cpp
|
||||
--- a/modules/fdlibm/src/s_sin.cpp
|
||||
+++ b/modules/fdlibm/src/s_sin.cpp
|
||||
diff --git a/s_sin.cpp b/s_sin.cpp
|
||||
--- a/s_sin.cpp
|
||||
+++ b/s_sin.cpp
|
||||
@@ -5,18 +5,18 @@
|
||||
*
|
||||
* Developed at SunPro, a Sun Microsystems, Inc. business.
|
||||
@ -882,9 +882,9 @@ diff --git a/modules/fdlibm/src/s_sin.cpp b/modules/fdlibm/src/s_sin.cpp
|
||||
* __kernel_sin ... sine function on [-pi/4,pi/4]
|
||||
* __kernel_cos ... cose function on [-pi/4,pi/4]
|
||||
* __ieee754_rem_pio2 ... argument reduction routine
|
||||
diff --git a/modules/fdlibm/src/s_tan.cpp b/modules/fdlibm/src/s_tan.cpp
|
||||
--- a/modules/fdlibm/src/s_tan.cpp
|
||||
+++ b/modules/fdlibm/src/s_tan.cpp
|
||||
diff --git a/s_tan.cpp b/s_tan.cpp
|
||||
--- a/s_tan.cpp
|
||||
+++ b/s_tan.cpp
|
||||
@@ -5,18 +5,18 @@
|
||||
*
|
||||
* Developed at SunPro, a Sun Microsystems, Inc. business.
|
||||
@ -906,9 +906,9 @@ diff --git a/modules/fdlibm/src/s_tan.cpp b/modules/fdlibm/src/s_tan.cpp
|
||||
* __kernel_tan ... tangent function on [-pi/4,pi/4]
|
||||
* __ieee754_rem_pio2 ... argument reduction routine
|
||||
*
|
||||
diff --git a/modules/fdlibm/src/s_tanh.cpp b/modules/fdlibm/src/s_tanh.cpp
|
||||
--- a/modules/fdlibm/src/s_tanh.cpp
|
||||
+++ b/modules/fdlibm/src/s_tanh.cpp
|
||||
diff --git a/s_tanh.cpp b/s_tanh.cpp
|
||||
--- a/s_tanh.cpp
|
||||
+++ b/s_tanh.cpp
|
||||
@@ -5,18 +5,18 @@
|
||||
*
|
||||
* Developed at SunPro, a Sun Microsystems, Inc. business.
|
||||
@ -930,9 +930,9 @@ diff --git a/modules/fdlibm/src/s_tanh.cpp b/modules/fdlibm/src/s_tanh.cpp
|
||||
* x -x
|
||||
* e - e
|
||||
* 0. tanh(x) is defined to be -----------
|
||||
diff --git a/modules/fdlibm/src/s_trunc.cpp b/modules/fdlibm/src/s_trunc.cpp
|
||||
--- a/modules/fdlibm/src/s_trunc.cpp
|
||||
+++ b/modules/fdlibm/src/s_trunc.cpp
|
||||
diff --git a/s_trunc.cpp b/s_trunc.cpp
|
||||
--- a/s_trunc.cpp
|
||||
+++ b/s_trunc.cpp
|
||||
@@ -5,18 +5,18 @@
|
||||
*
|
||||
* Developed at SunPro, a Sun Microsystems, Inc. business.
|
||||
@ -954,9 +954,9 @@ diff --git a/modules/fdlibm/src/s_trunc.cpp b/modules/fdlibm/src/s_trunc.cpp
|
||||
* Bit twiddling.
|
||||
* Exception:
|
||||
* Inexact flag raised if x not equal to trunc(x).
|
||||
diff --git a/modules/fdlibm/src/s_truncf.cpp b/modules/fdlibm/src/s_truncf.cpp
|
||||
--- a/modules/fdlibm/src/s_truncf.cpp
|
||||
+++ b/modules/fdlibm/src/s_truncf.cpp
|
||||
diff --git a/s_truncf.cpp b/s_truncf.cpp
|
||||
--- a/s_truncf.cpp
|
||||
+++ b/s_truncf.cpp
|
||||
@@ -5,18 +5,18 @@
|
||||
*
|
||||
* Developed at SunPro, a Sun Microsystems, Inc. business.
|
||||
|
@ -1,6 +1,6 @@
|
||||
diff --git a/modules/fdlibm/src/k_exp.cpp b/modules/fdlibm/src/k_exp.cpp
|
||||
--- a/modules/fdlibm/src/k_exp.cpp
|
||||
+++ b/modules/fdlibm/src/k_exp.cpp
|
||||
diff --git a/k_exp.cpp b/k_exp.cpp
|
||||
--- a/k_exp.cpp
|
||||
+++ b/k_exp.cpp
|
||||
@@ -24,18 +24,16 @@
|
||||
* LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
|
||||
* OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
|
||||
|
@ -1,6 +1,6 @@
|
||||
diff --git a/modules/fdlibm/src/math_private.h b/modules/fdlibm/src/math_private.h
|
||||
--- a/modules/fdlibm/src/math_private.h
|
||||
+++ b/modules/fdlibm/src/math_private.h
|
||||
diff --git a/math_private.h b/math_private.h
|
||||
--- a/math_private.h
|
||||
+++ b/math_private.h
|
||||
@@ -12,16 +12,17 @@
|
||||
/*
|
||||
* from: @(#)fdlibm.h 5.1 93/09/24
|
||||
|
@ -1,6 +1,6 @@
|
||||
diff --git a/modules/fdlibm/src/math_private.h b/modules/fdlibm/src/math_private.h
|
||||
--- a/modules/fdlibm/src/math_private.h
|
||||
+++ b/modules/fdlibm/src/math_private.h
|
||||
diff --git a/math_private.h b/math_private.h
|
||||
--- a/math_private.h
|
||||
+++ b/math_private.h
|
||||
@@ -33,16 +33,23 @@
|
||||
* to dig two 32 bit words out of the 64 bit IEEE floating point
|
||||
* value. That is non-ANSI, and, moreover, the gcc instruction
|
||||
|
@ -1,6 +1,6 @@
|
||||
diff --git a/modules/fdlibm/src/math_private.h b/modules/fdlibm/src/math_private.h
|
||||
--- a/modules/fdlibm/src/math_private.h
|
||||
+++ b/modules/fdlibm/src/math_private.h
|
||||
diff --git a/math_private.h b/math_private.h
|
||||
--- a/math_private.h
|
||||
+++ b/math_private.h
|
||||
@@ -323,16 +323,27 @@ do { \
|
||||
if (sizeof(type) >= sizeof(long double)) \
|
||||
(lval) = (rval); \
|
||||
|
@ -1,6 +1,6 @@
|
||||
diff --git a/modules/fdlibm/src/e_exp.cpp b/modules/fdlibm/src/e_exp.cpp
|
||||
--- a/modules/fdlibm/src/e_exp.cpp
|
||||
+++ b/modules/fdlibm/src/e_exp.cpp
|
||||
diff --git a/e_exp.cpp b/e_exp.cpp
|
||||
--- a/e_exp.cpp
|
||||
+++ b/e_exp.cpp
|
||||
@@ -146,14 +146,17 @@ __ieee754_exp(double x) /* default IEEE double exp */
|
||||
if(k >= -1021)
|
||||
INSERT_WORDS(twopk,((u_int32_t)(0x3ff+k))<<20, 0);
|
||||
@ -20,9 +20,9 @@ diff --git a/modules/fdlibm/src/e_exp.cpp b/modules/fdlibm/src/e_exp.cpp
|
||||
return y*twopk*twom1000;
|
||||
}
|
||||
}
|
||||
diff --git a/modules/fdlibm/src/s_expm1.cpp b/modules/fdlibm/src/s_expm1.cpp
|
||||
--- a/modules/fdlibm/src/s_expm1.cpp
|
||||
+++ b/modules/fdlibm/src/s_expm1.cpp
|
||||
diff --git a/s_expm1.cpp b/s_expm1.cpp
|
||||
--- a/s_expm1.cpp
|
||||
+++ b/s_expm1.cpp
|
||||
@@ -192,17 +192,20 @@ expm1(double x)
|
||||
e -= hxs;
|
||||
if(k== -1) return 0.5*(x-e)-0.5;
|
||||
|
@ -1,6 +1,6 @@
|
||||
diff --git a/modules/fdlibm/src/s_nearbyint.cpp b/modules/fdlibm/src/s_nearbyint.cpp
|
||||
--- a/modules/fdlibm/src/s_nearbyint.cpp
|
||||
+++ b/modules/fdlibm/src/s_nearbyint.cpp
|
||||
diff --git a/s_nearbyint.cpp b/s_nearbyint.cpp
|
||||
--- a/s_nearbyint.cpp
|
||||
+++ b/s_nearbyint.cpp
|
||||
@@ -54,9 +54,8 @@ fn(type x) \
|
||||
fegetenv(&env); \
|
||||
ret = rint(x); \
|
||||
|
@ -1,6 +1,6 @@
|
||||
diff --git a/modules/fdlibm/src/math_private.h b/modules/fdlibm/src/math_private.h
|
||||
--- a/modules/fdlibm/src/math_private.h
|
||||
+++ b/modules/fdlibm/src/math_private.h
|
||||
diff --git a/math_private.h b/math_private.h
|
||||
--- a/math_private.h
|
||||
+++ b/math_private.h
|
||||
@@ -309,17 +309,17 @@ do { \
|
||||
/* The above works on non-i386 too, but we use this to check v. */
|
||||
#define LD80C(m, ex, v) { .e = (v), }
|
||||
|
@ -1,6 +1,6 @@
|
||||
diff --git a/modules/fdlibm/src/e_exp.cpp b/modules/fdlibm/src/e_exp.cpp
|
||||
--- a/modules/fdlibm/src/e_exp.cpp
|
||||
+++ b/modules/fdlibm/src/e_exp.cpp
|
||||
diff --git a/e_exp.cpp b/e_exp.cpp
|
||||
--- a/e_exp.cpp
|
||||
+++ b/e_exp.cpp
|
||||
@@ -91,16 +91,18 @@ ln2LO[2] ={ 1.90821492927058770002e-10, /* 0x3dea39ef, 0x35793c76 */
|
||||
-1.90821492927058770002e-10,},/* 0xbdea39ef, 0x35793c76 */
|
||||
invln2 = 1.44269504088896338700e+00, /* 0x3ff71547, 0x652b82fe */
|
||||
|
@ -1,6 +1,6 @@
|
||||
diff --git a/modules/fdlibm/src/e_acos.cpp b/modules/fdlibm/src/e_acos.cpp
|
||||
--- a/modules/fdlibm/src/e_acos.cpp
|
||||
+++ b/modules/fdlibm/src/e_acos.cpp
|
||||
diff --git a/e_acos.cpp b/e_acos.cpp
|
||||
--- a/e_acos.cpp
|
||||
+++ b/e_acos.cpp
|
||||
@@ -33,16 +33,17 @@
|
||||
*
|
||||
* Special cases:
|
||||
@ -45,9 +45,9 @@ diff --git a/modules/fdlibm/src/e_acos.cpp b/modules/fdlibm/src/e_acos.cpp
|
||||
r = p/q;
|
||||
w = r*s+c;
|
||||
return 2.0*(df+w);
|
||||
diff --git a/modules/fdlibm/src/e_acosh.cpp b/modules/fdlibm/src/e_acosh.cpp
|
||||
--- a/modules/fdlibm/src/e_acosh.cpp
|
||||
+++ b/modules/fdlibm/src/e_acosh.cpp
|
||||
diff --git a/e_acosh.cpp b/e_acosh.cpp
|
||||
--- a/e_acosh.cpp
|
||||
+++ b/e_acosh.cpp
|
||||
@@ -24,16 +24,17 @@
|
||||
* acosh(x) := log(2x-1/(sqrt(x*x-1)+x)) if x>2; else
|
||||
* acosh(x) := log1p(t+sqrt(2.0*t+t*t)); where t=x-1.
|
||||
@ -83,9 +83,9 @@ diff --git a/modules/fdlibm/src/e_acosh.cpp b/modules/fdlibm/src/e_acosh.cpp
|
||||
+ return log1p(t+std::sqrt(2.0*t+t*t));
|
||||
}
|
||||
}
|
||||
diff --git a/modules/fdlibm/src/e_asin.cpp b/modules/fdlibm/src/e_asin.cpp
|
||||
--- a/modules/fdlibm/src/e_asin.cpp
|
||||
+++ b/modules/fdlibm/src/e_asin.cpp
|
||||
diff --git a/e_asin.cpp b/e_asin.cpp
|
||||
--- a/e_asin.cpp
|
||||
+++ b/e_asin.cpp
|
||||
@@ -39,16 +39,17 @@
|
||||
* = pio4_hi+(pio4-2f)-(2s*z*R(z)-(pio2_lo+2c))
|
||||
*
|
||||
@ -123,9 +123,9 @@ diff --git a/modules/fdlibm/src/e_asin.cpp b/modules/fdlibm/src/e_asin.cpp
|
||||
SET_LOW_WORD(w,0);
|
||||
c = (t-w*w)/(s+w);
|
||||
r = p/q;
|
||||
diff --git a/modules/fdlibm/src/e_hypot.cpp b/modules/fdlibm/src/e_hypot.cpp
|
||||
--- a/modules/fdlibm/src/e_hypot.cpp
|
||||
+++ b/modules/fdlibm/src/e_hypot.cpp
|
||||
diff --git a/e_hypot.cpp b/e_hypot.cpp
|
||||
--- a/e_hypot.cpp
|
||||
+++ b/e_hypot.cpp
|
||||
@@ -41,16 +41,17 @@
|
||||
* hypot(x,y) is INF if x or y is +INF or -INF; else
|
||||
* hypot(x,y) is NAN if x or y is NAN.
|
||||
@ -172,9 +172,9 @@ diff --git a/modules/fdlibm/src/e_hypot.cpp b/modules/fdlibm/src/e_hypot.cpp
|
||||
return t1*w;
|
||||
} else return w;
|
||||
}
|
||||
diff --git a/modules/fdlibm/src/e_pow.cpp b/modules/fdlibm/src/e_pow.cpp
|
||||
--- a/modules/fdlibm/src/e_pow.cpp
|
||||
+++ b/modules/fdlibm/src/e_pow.cpp
|
||||
diff --git a/e_pow.cpp b/e_pow.cpp
|
||||
--- a/e_pow.cpp
|
||||
+++ b/e_pow.cpp
|
||||
@@ -52,16 +52,17 @@
|
||||
*
|
||||
* Constants :
|
||||
@ -212,9 +212,9 @@ diff --git a/modules/fdlibm/src/e_pow.cpp b/modules/fdlibm/src/e_pow.cpp
|
||||
if(lx==0) {
|
||||
if(ix==0x7ff00000||ix==0||ix==0x3ff00000){
|
||||
z = ax; /*x is +-0,+-inf,+-1*/
|
||||
diff --git a/modules/fdlibm/src/s_asinh.cpp b/modules/fdlibm/src/s_asinh.cpp
|
||||
--- a/modules/fdlibm/src/s_asinh.cpp
|
||||
+++ b/modules/fdlibm/src/s_asinh.cpp
|
||||
diff --git a/s_asinh.cpp b/s_asinh.cpp
|
||||
--- a/s_asinh.cpp
|
||||
+++ b/s_asinh.cpp
|
||||
@@ -19,16 +19,17 @@
|
||||
* asinh(x) = sign(x) * log [ |x| + sqrt(x*x+1) ]
|
||||
* we have
|
||||
|
@ -1,6 +1,6 @@
|
||||
diff --git a/modules/fdlibm/src/math_private.h b/modules/fdlibm/src/math_private.h
|
||||
--- a/modules/fdlibm/src/math_private.h
|
||||
+++ b/modules/fdlibm/src/math_private.h
|
||||
diff --git a/math_private.h b/math_private.h
|
||||
--- a/math_private.h
|
||||
+++ b/math_private.h
|
||||
@@ -612,46 +612,16 @@ rnint(__double_t x)
|
||||
* magic number would need to be variable. Assuming that the
|
||||
* rounding precision is always the default is too fragile. This
|
||||
|
@ -1,6 +1,6 @@
|
||||
diff --git a/modules/fdlibm/src/math_private.h b/modules/fdlibm/src/math_private.h
|
||||
--- a/modules/fdlibm/src/math_private.h
|
||||
+++ b/modules/fdlibm/src/math_private.h
|
||||
diff --git a/math_private.h b/math_private.h
|
||||
--- a/math_private.h
|
||||
+++ b/math_private.h
|
||||
@@ -20,16 +20,24 @@
|
||||
#include <cfloat>
|
||||
#include <stdint.h>
|
||||
|
@ -1,6 +1,6 @@
|
||||
diff --git a/modules/fdlibm/src/e_rem_pio2.cpp b/modules/fdlibm/src/e_rem_pio2.cpp
|
||||
--- a/modules/fdlibm/src/e_rem_pio2.cpp
|
||||
+++ b/modules/fdlibm/src/e_rem_pio2.cpp
|
||||
diff --git a/e_rem_pio2.cpp b/e_rem_pio2.cpp
|
||||
--- a/e_rem_pio2.cpp
|
||||
+++ b/e_rem_pio2.cpp
|
||||
@@ -43,17 +43,17 @@ invpio2 = 6.36619772367581382433e-01, /* 0x3FE45F30, 0x6DC9C883 */
|
||||
pio2_1 = 1.57079632673412561417e+00, /* 0x3FF921FB, 0x54400000 */
|
||||
pio2_1t = 6.07710050650619224932e-11, /* 0x3DD0B461, 0x1A626331 */
|
||||
@ -20,9 +20,9 @@ diff --git a/modules/fdlibm/src/e_rem_pio2.cpp b/modules/fdlibm/src/e_rem_pio2.c
|
||||
double tx[3],ty[2];
|
||||
int32_t e0,i,j,nx,n,ix,hx;
|
||||
u_int32_t low;
|
||||
diff --git a/modules/fdlibm/src/k_rem_pio2.cpp b/modules/fdlibm/src/k_rem_pio2.cpp
|
||||
--- a/modules/fdlibm/src/k_rem_pio2.cpp
|
||||
+++ b/modules/fdlibm/src/k_rem_pio2.cpp
|
||||
diff --git a/k_rem_pio2.cpp b/k_rem_pio2.cpp
|
||||
--- a/k_rem_pio2.cpp
|
||||
+++ b/k_rem_pio2.cpp
|
||||
@@ -305,17 +305,18 @@ __kernel_rem_pio2(double *x, double *y, int e0, int nx, int prec)
|
||||
q0 = e0-24*(jv+1);
|
||||
|
||||
@ -43,9 +43,9 @@ diff --git a/modules/fdlibm/src/k_rem_pio2.cpp b/modules/fdlibm/src/k_rem_pio2.c
|
||||
for(i=0,j=jz,z=q[jz];j>0;i++,j--) {
|
||||
fw = (double)((int32_t)(twon24* z));
|
||||
iq[i] = (int32_t)(z-two24*fw);
|
||||
diff --git a/modules/fdlibm/src/s_cos.cpp b/modules/fdlibm/src/s_cos.cpp
|
||||
--- a/modules/fdlibm/src/s_cos.cpp
|
||||
+++ b/modules/fdlibm/src/s_cos.cpp
|
||||
diff --git a/s_cos.cpp b/s_cos.cpp
|
||||
--- a/s_cos.cpp
|
||||
+++ b/s_cos.cpp
|
||||
@@ -43,17 +43,17 @@
|
||||
* Accuracy:
|
||||
* TRIG(x) returns trig(x) nearly rounded
|
||||
@ -65,9 +65,9 @@ diff --git a/modules/fdlibm/src/s_cos.cpp b/modules/fdlibm/src/s_cos.cpp
|
||||
int32_t n, ix;
|
||||
|
||||
/* High word of x. */
|
||||
diff --git a/modules/fdlibm/src/s_sin.cpp b/modules/fdlibm/src/s_sin.cpp
|
||||
--- a/modules/fdlibm/src/s_sin.cpp
|
||||
+++ b/modules/fdlibm/src/s_sin.cpp
|
||||
diff --git a/s_sin.cpp b/s_sin.cpp
|
||||
--- a/s_sin.cpp
|
||||
+++ b/s_sin.cpp
|
||||
@@ -43,17 +43,17 @@
|
||||
* Accuracy:
|
||||
* TRIG(x) returns trig(x) nearly rounded
|
||||
@ -87,9 +87,9 @@ diff --git a/modules/fdlibm/src/s_sin.cpp b/modules/fdlibm/src/s_sin.cpp
|
||||
int32_t n, ix;
|
||||
|
||||
/* High word of x. */
|
||||
diff --git a/modules/fdlibm/src/s_tan.cpp b/modules/fdlibm/src/s_tan.cpp
|
||||
--- a/modules/fdlibm/src/s_tan.cpp
|
||||
+++ b/modules/fdlibm/src/s_tan.cpp
|
||||
diff --git a/s_tan.cpp b/s_tan.cpp
|
||||
--- a/s_tan.cpp
|
||||
+++ b/s_tan.cpp
|
||||
@@ -42,17 +42,17 @@
|
||||
* Accuracy:
|
||||
* TRIG(x) returns trig(x) nearly rounded
|
||||
|
@ -1,6 +1,6 @@
|
||||
diff --git a/modules/fdlibm/src/math_private.h b/modules/fdlibm/src/math_private.h
|
||||
--- a/modules/fdlibm/src/math_private.h
|
||||
+++ b/modules/fdlibm/src/math_private.h
|
||||
diff --git a/math_private.h b/math_private.h
|
||||
--- a/math_private.h
|
||||
+++ b/math_private.h
|
||||
@@ -27,16 +27,17 @@
|
||||
|
||||
/*
|
||||
|
Loading…
Reference in New Issue
Block a user