зеркало из https://github.com/mozilla/gecko-dev.git
837 строки
25 KiB
Diff
837 строки
25 KiB
Diff
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
|
|
@@ -6,18 +6,18 @@
|
|
*
|
|
* Developed at SunSoft, a Sun Microsystems, Inc. business.
|
|
* Permission to use, copy, modify, and distribute this
|
|
* software is freely granted, provided that this notice
|
|
* is preserved.
|
|
* ====================================================
|
|
*/
|
|
|
|
-#include <sys/cdefs.h>
|
|
-__FBSDID("$FreeBSD$");
|
|
+//#include <sys/cdefs.h>
|
|
+//__FBSDID("$FreeBSD$");
|
|
|
|
/* __ieee754_acos(x)
|
|
* Method :
|
|
* acos(x) = pi/2 - asin(x)
|
|
* acos(-x) = pi/2 + asin(x)
|
|
* 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
|
|
@@ -7,18 +7,18 @@
|
|
* Developed at SunSoft, a Sun Microsystems, Inc. business.
|
|
* Permission to use, copy, modify, and distribute this
|
|
* software is freely granted, provided that this notice
|
|
* is preserved.
|
|
* ====================================================
|
|
*
|
|
*/
|
|
|
|
-#include <sys/cdefs.h>
|
|
-__FBSDID("$FreeBSD$");
|
|
+//#include <sys/cdefs.h>
|
|
+//__FBSDID("$FreeBSD$");
|
|
|
|
/* __ieee754_acosh(x)
|
|
* Method :
|
|
* Based on
|
|
* acosh(x) = log [ x + sqrt(x*x-1) ]
|
|
* 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
|
|
@@ -6,18 +6,18 @@
|
|
*
|
|
* Developed at SunSoft, a Sun Microsystems, Inc. business.
|
|
* Permission to use, copy, modify, and distribute this
|
|
* software is freely granted, provided that this notice
|
|
* is preserved.
|
|
* ====================================================
|
|
*/
|
|
|
|
-#include <sys/cdefs.h>
|
|
-__FBSDID("$FreeBSD$");
|
|
+//#include <sys/cdefs.h>
|
|
+//__FBSDID("$FreeBSD$");
|
|
|
|
/* __ieee754_asin(x)
|
|
* Method :
|
|
* Since asin(x) = x + x^3/6 + x^5*3/40 + x^7*15/336 + ...
|
|
* we approximate asin(x) on [0,0.5] by
|
|
* 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
|
|
@@ -7,18 +7,18 @@
|
|
* Developed at SunSoft, a Sun Microsystems, Inc. business.
|
|
* Permission to use, copy, modify, and distribute this
|
|
* software is freely granted, provided that this notice
|
|
* is preserved.
|
|
* ====================================================
|
|
*
|
|
*/
|
|
|
|
-#include <sys/cdefs.h>
|
|
-__FBSDID("$FreeBSD$");
|
|
+//#include <sys/cdefs.h>
|
|
+//__FBSDID("$FreeBSD$");
|
|
|
|
/* __ieee754_atan2(y,x)
|
|
* Method :
|
|
* 1. Reduce y to positive by atan2(y,x)=-atan2(-y,x).
|
|
* 2. Reduce x to positive by (if x and y are unexceptional):
|
|
* 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
|
|
@@ -7,18 +7,18 @@
|
|
* Developed at SunSoft, a Sun Microsystems, Inc. business.
|
|
* Permission to use, copy, modify, and distribute this
|
|
* software is freely granted, provided that this notice
|
|
* is preserved.
|
|
* ====================================================
|
|
*
|
|
*/
|
|
|
|
-#include <sys/cdefs.h>
|
|
-__FBSDID("$FreeBSD$");
|
|
+//#include <sys/cdefs.h>
|
|
+//__FBSDID("$FreeBSD$");
|
|
|
|
/* __ieee754_atanh(x)
|
|
* Method :
|
|
* 1.Reduced x to positive by atanh(-x) = -atanh(x)
|
|
* 2.For x>=0.5
|
|
* 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
|
|
@@ -6,18 +6,18 @@
|
|
*
|
|
* Developed at SunSoft, a Sun Microsystems, Inc. business.
|
|
* Permission to use, copy, modify, and distribute this
|
|
* software is freely granted, provided that this notice
|
|
* is preserved.
|
|
* ====================================================
|
|
*/
|
|
|
|
-#include <sys/cdefs.h>
|
|
-__FBSDID("$FreeBSD$");
|
|
+//#include <sys/cdefs.h>
|
|
+//__FBSDID("$FreeBSD$");
|
|
|
|
/* __ieee754_cosh(x)
|
|
* Method :
|
|
* mathematically cosh(x) if defined to be (exp(x)+exp(-x))/2
|
|
* 1. Replace x by |x| (cosh(x) = cosh(-x)).
|
|
* 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
|
|
@@ -5,18 +5,18 @@
|
|
* Copyright (C) 2004 by Sun Microsystems, Inc. All rights reserved.
|
|
*
|
|
* Permission to use, copy, modify, and distribute this
|
|
* software is freely granted, provided that this notice
|
|
* is preserved.
|
|
* ====================================================
|
|
*/
|
|
|
|
-#include <sys/cdefs.h>
|
|
-__FBSDID("$FreeBSD$");
|
|
+//#include <sys/cdefs.h>
|
|
+//__FBSDID("$FreeBSD$");
|
|
|
|
/* __ieee754_exp(x)
|
|
* Returns the exponential of x.
|
|
*
|
|
* Method
|
|
* 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
|
|
@@ -6,18 +6,18 @@
|
|
*
|
|
* Developed at SunSoft, a Sun Microsystems, Inc. business.
|
|
* Permission to use, copy, modify, and distribute this
|
|
* software is freely granted, provided that this notice
|
|
* is preserved.
|
|
* ====================================================
|
|
*/
|
|
|
|
-#include <sys/cdefs.h>
|
|
-__FBSDID("$FreeBSD$");
|
|
+//#include <sys/cdefs.h>
|
|
+//__FBSDID("$FreeBSD$");
|
|
|
|
/* __ieee754_hypot(x,y)
|
|
*
|
|
* Method :
|
|
* If (assume round-to-nearest) z=x*x+y*y
|
|
* 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
|
|
@@ -6,18 +6,18 @@
|
|
*
|
|
* Developed at SunSoft, a Sun Microsystems, Inc. business.
|
|
* Permission to use, copy, modify, and distribute this
|
|
* software is freely granted, provided that this notice
|
|
* is preserved.
|
|
* ====================================================
|
|
*/
|
|
|
|
-#include <sys/cdefs.h>
|
|
-__FBSDID("$FreeBSD$");
|
|
+//#include <sys/cdefs.h>
|
|
+//__FBSDID("$FreeBSD$");
|
|
|
|
/* __ieee754_log(x)
|
|
* Return the logrithm of x
|
|
*
|
|
* Method :
|
|
* 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
|
|
@@ -6,32 +6,32 @@
|
|
*
|
|
* Developed at SunSoft, a Sun Microsystems, Inc. business.
|
|
* Permission to use, copy, modify, and distribute this
|
|
* software is freely granted, provided that this notice
|
|
* is preserved.
|
|
* ====================================================
|
|
*/
|
|
|
|
-#include <sys/cdefs.h>
|
|
-__FBSDID("$FreeBSD$");
|
|
+//#include <sys/cdefs.h>
|
|
+//__FBSDID("$FreeBSD$");
|
|
|
|
/*
|
|
* Return the base 10 logarithm of x. See e_log.c and k_log.h for most
|
|
* comments.
|
|
*
|
|
* log10(x) = (f - 0.5*f*f + k_log1p(f)) / ln10 + k * log10(2)
|
|
* in not-quite-routine extra precision.
|
|
*/
|
|
|
|
#include <float.h>
|
|
|
|
#include "math_private.h"
|
|
#include "k_log.h"
|
|
-
|
|
+
|
|
static const double
|
|
two54 = 1.80143985094819840000e+16, /* 0x43500000, 0x00000000 */
|
|
ivln10hi = 4.34294481878168880939e-01, /* 0x3fdbcb7b, 0x15200000 */
|
|
ivln10lo = 2.50829467116452752298e-11, /* 0x3dbb9438, 0xca9aadd5 */
|
|
log10_2hi = 3.01029995663611771306e-01, /* 0x3FD34413, 0x509F6000 */
|
|
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
|
|
@@ -6,18 +6,18 @@
|
|
*
|
|
* Developed at SunSoft, a Sun Microsystems, Inc. business.
|
|
* Permission to use, copy, modify, and distribute this
|
|
* software is freely granted, provided that this notice
|
|
* is preserved.
|
|
* ====================================================
|
|
*/
|
|
|
|
-#include <sys/cdefs.h>
|
|
-__FBSDID("$FreeBSD$");
|
|
+//#include <sys/cdefs.h>
|
|
+//__FBSDID("$FreeBSD$");
|
|
|
|
/*
|
|
* Return the base 2 logarithm of x. See e_log.c and k_log.h for most
|
|
* comments.
|
|
*
|
|
* 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
|
|
@@ -4,18 +4,18 @@
|
|
* Copyright (C) 2004 by Sun Microsystems, Inc. All rights reserved.
|
|
*
|
|
* Permission to use, copy, modify, and distribute this
|
|
* software is freely granted, provided that this notice
|
|
* is preserved.
|
|
* ====================================================
|
|
*/
|
|
|
|
-#include <sys/cdefs.h>
|
|
-__FBSDID("$FreeBSD$");
|
|
+//#include <sys/cdefs.h>
|
|
+//__FBSDID("$FreeBSD$");
|
|
|
|
/* __ieee754_pow(x,y) return x**y
|
|
*
|
|
* n
|
|
* Method: Let x = 2 * (1+f)
|
|
* 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_sinh.cpp b/modules/fdlibm/src/e_sinh.cpp
|
|
--- a/modules/fdlibm/src/e_sinh.cpp
|
|
+++ b/modules/fdlibm/src/e_sinh.cpp
|
|
@@ -6,18 +6,18 @@
|
|
*
|
|
* Developed at SunSoft, a Sun Microsystems, Inc. business.
|
|
* Permission to use, copy, modify, and distribute this
|
|
* software is freely granted, provided that this notice
|
|
* is preserved.
|
|
* ====================================================
|
|
*/
|
|
|
|
-#include <sys/cdefs.h>
|
|
-__FBSDID("$FreeBSD$");
|
|
+//#include <sys/cdefs.h>
|
|
+//__FBSDID("$FreeBSD$");
|
|
|
|
/* __ieee754_sinh(x)
|
|
* Method :
|
|
* mathematically sinh(x) if defined to be (exp(x)-exp(-x))/2
|
|
* 1. Replace x by |x| (sinh(-x) = -sinh(x)).
|
|
* 2.
|
|
* E + E/(E+1)
|
|
* 0 <= x <= 22 : sinh(x) := --------------, E=expm1(x)
|
|
diff --git a/modules/fdlibm/src/e_sqrt.cpp b/modules/fdlibm/src/e_sqrt.cpp
|
|
--- a/modules/fdlibm/src/e_sqrt.cpp
|
|
+++ b/modules/fdlibm/src/e_sqrt.cpp
|
|
@@ -6,18 +6,18 @@
|
|
*
|
|
* Developed at SunSoft, a Sun Microsystems, Inc. business.
|
|
* Permission to use, copy, modify, and distribute this
|
|
* software is freely granted, provided that this notice
|
|
* is preserved.
|
|
* ====================================================
|
|
*/
|
|
|
|
-#include <sys/cdefs.h>
|
|
-__FBSDID("$FreeBSD$");
|
|
+//#include <sys/cdefs.h>
|
|
+//__FBSDID("$FreeBSD$");
|
|
|
|
/* __ieee754_sqrt(x)
|
|
* Return correctly rounded sqrt.
|
|
* ------------------------------------------
|
|
* | Use the hardware sqrt if you have one |
|
|
* ------------------------------------------
|
|
* Method:
|
|
* Bit by bit method using integer arithmetic. (Slow, but portable)
|
|
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
|
|
@@ -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)
|
|
* HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
|
|
* 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
|
|
* SUCH DAMAGE.
|
|
*/
|
|
|
|
-#include <sys/cdefs.h>
|
|
-__FBSDID("$FreeBSD$");
|
|
+//#include <sys/cdefs.h>
|
|
+//__FBSDID("$FreeBSD$");
|
|
|
|
#include <complex.h>
|
|
|
|
-#include "math_private.h"
|
|
+ #include "math_private.h"
|
|
|
|
static const uint32_t k = 1799; /* constant for reduction */
|
|
static const double kln2 = 1246.97177782734161156; /* k * ln2 */
|
|
|
|
/*
|
|
* 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
|
|
@@ -6,18 +6,18 @@
|
|
*
|
|
* Developed at SunSoft, a Sun Microsystems, Inc. business.
|
|
* Permission to use, copy, modify, and distribute this
|
|
* software is freely granted, provided that this notice
|
|
* is preserved.
|
|
* ====================================================
|
|
*/
|
|
|
|
-#include <sys/cdefs.h>
|
|
-__FBSDID("$FreeBSD$");
|
|
+//#include <sys/cdefs.h>
|
|
+//__FBSDID("$FreeBSD$");
|
|
|
|
/*
|
|
* k_log1p(f):
|
|
* Return log(1+f) - f for 1+f in ~[sqrt(2)/2, sqrt(2)].
|
|
*
|
|
* 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/s_asinh.cpp b/modules/fdlibm/src/s_asinh.cpp
|
|
--- a/modules/fdlibm/src/s_asinh.cpp
|
|
+++ b/modules/fdlibm/src/s_asinh.cpp
|
|
@@ -5,18 +5,18 @@
|
|
*
|
|
* Developed at SunPro, a Sun Microsystems, Inc. business.
|
|
* Permission to use, copy, modify, and distribute this
|
|
* software is freely granted, provided that this notice
|
|
* is preserved.
|
|
* ====================================================
|
|
*/
|
|
|
|
-#include <sys/cdefs.h>
|
|
-__FBSDID("$FreeBSD$");
|
|
+//#include <sys/cdefs.h>
|
|
+//__FBSDID("$FreeBSD$");
|
|
|
|
/* asinh(x)
|
|
* Method :
|
|
* Based on
|
|
* asinh(x) = sign(x) * log [ |x| + sqrt(x*x+1) ]
|
|
* 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
|
|
@@ -5,18 +5,18 @@
|
|
*
|
|
* Developed at SunPro, a Sun Microsystems, Inc. business.
|
|
* Permission to use, copy, modify, and distribute this
|
|
* software is freely granted, provided that this notice
|
|
* is preserved.
|
|
* ====================================================
|
|
*/
|
|
|
|
-#include <sys/cdefs.h>
|
|
-__FBSDID("$FreeBSD$");
|
|
+//#include <sys/cdefs.h>
|
|
+//__FBSDID("$FreeBSD$");
|
|
|
|
/* atan(x)
|
|
* Method
|
|
* 1. Reduce x to positive by atan(x) = -atan(-x).
|
|
* 2. According to the integer k=4t+0.25 chopped, t=x, the argument
|
|
* 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
|
|
@@ -7,18 +7,18 @@
|
|
* Permission to use, copy, modify, and distribute this
|
|
* software is freely granted, provided that this notice
|
|
* is preserved.
|
|
* ====================================================
|
|
*
|
|
* Optimized by Bruce D. Evans.
|
|
*/
|
|
|
|
-#include <sys/cdefs.h>
|
|
-__FBSDID("$FreeBSD$");
|
|
+//#include <sys/cdefs.h>
|
|
+//__FBSDID("$FreeBSD$");
|
|
|
|
#include "math_private.h"
|
|
|
|
/* cbrt(x)
|
|
* Return cube root of x
|
|
*/
|
|
static const u_int32_t
|
|
B1 = 715094163, /* B1 = (1023-1023/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
|
|
@@ -5,18 +5,18 @@
|
|
*
|
|
* Developed at SunPro, a Sun Microsystems, Inc. business.
|
|
* Permission to use, copy, modify, and distribute this
|
|
* software is freely granted, provided that this notice
|
|
* is preserved.
|
|
* ====================================================
|
|
*/
|
|
|
|
-#include <sys/cdefs.h>
|
|
-__FBSDID("$FreeBSD$");
|
|
+//#include <sys/cdefs.h>
|
|
+//__FBSDID("$FreeBSD$");
|
|
|
|
/*
|
|
* ceil(x)
|
|
* Return x rounded toward -inf to integral value
|
|
* Method:
|
|
* 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
|
|
@@ -8,18 +8,18 @@
|
|
*
|
|
* Developed at SunPro, a Sun Microsystems, Inc. business.
|
|
* Permission to use, copy, modify, and distribute this
|
|
* software is freely granted, provided that this notice
|
|
* is preserved.
|
|
* ====================================================
|
|
*/
|
|
|
|
-#include <sys/cdefs.h>
|
|
-__FBSDID("$FreeBSD$");
|
|
+//#include <sys/cdefs.h>
|
|
+//__FBSDID("$FreeBSD$");
|
|
|
|
#include "math_private.h"
|
|
|
|
static const float huge = 1.0e30;
|
|
|
|
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
|
|
@@ -5,18 +5,18 @@
|
|
*
|
|
* Developed at SunPro, a Sun Microsystems, Inc. business.
|
|
* Permission to use, copy, modify, and distribute this
|
|
* software is freely granted, provided that this notice
|
|
* is preserved.
|
|
* ====================================================
|
|
*/
|
|
|
|
-#include <sys/cdefs.h>
|
|
-__FBSDID("$FreeBSD$");
|
|
+//#include <sys/cdefs.h>
|
|
+//__FBSDID("$FreeBSD$");
|
|
|
|
/*
|
|
* copysign(double x, double y)
|
|
* copysign(x,y) returns a value with the magnitude of x and
|
|
* with the sign bit of y.
|
|
*/
|
|
|
|
#include "math_private.h"
|
|
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
|
|
@@ -5,18 +5,18 @@
|
|
*
|
|
* Developed at SunPro, a Sun Microsystems, Inc. business.
|
|
* Permission to use, copy, modify, and distribute this
|
|
* software is freely granted, provided that this notice
|
|
* is preserved.
|
|
* ====================================================
|
|
*/
|
|
|
|
-#include <sys/cdefs.h>
|
|
-__FBSDID("$FreeBSD$");
|
|
+//#include <sys/cdefs.h>
|
|
+//__FBSDID("$FreeBSD$");
|
|
|
|
/* expm1(x)
|
|
* Returns exp(x)-1, the exponential of x minus 1.
|
|
*
|
|
* Method
|
|
* 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
|
|
@@ -5,18 +5,18 @@
|
|
*
|
|
* Developed at SunPro, a Sun Microsystems, Inc. business.
|
|
* Permission to use, copy, modify, and distribute this
|
|
* software is freely granted, provided that this notice
|
|
* is preserved.
|
|
* ====================================================
|
|
*/
|
|
|
|
-#include <sys/cdefs.h>
|
|
-__FBSDID("$FreeBSD$");
|
|
+//#include <sys/cdefs.h>
|
|
+//__FBSDID("$FreeBSD$");
|
|
|
|
/*
|
|
* fabs(x) returns the absolute value of x.
|
|
*/
|
|
|
|
#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
|
|
@@ -5,18 +5,18 @@
|
|
*
|
|
* Developed at SunPro, a Sun Microsystems, Inc. business.
|
|
* Permission to use, copy, modify, and distribute this
|
|
* software is freely granted, provided that this notice
|
|
* is preserved.
|
|
* ====================================================
|
|
*/
|
|
|
|
-#include <sys/cdefs.h>
|
|
-__FBSDID("$FreeBSD$");
|
|
+//#include <sys/cdefs.h>
|
|
+//__FBSDID("$FreeBSD$");
|
|
|
|
/*
|
|
* floor(x)
|
|
* Return x rounded toward -inf to integral value
|
|
* Method:
|
|
* 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
|
|
@@ -8,18 +8,18 @@
|
|
*
|
|
* Developed at SunPro, a Sun Microsystems, Inc. business.
|
|
* Permission to use, copy, modify, and distribute this
|
|
* software is freely granted, provided that this notice
|
|
* is preserved.
|
|
* ====================================================
|
|
*/
|
|
|
|
-#include <sys/cdefs.h>
|
|
-__FBSDID("$FreeBSD$");
|
|
+//#include <sys/cdefs.h>
|
|
+//__FBSDID("$FreeBSD$");
|
|
|
|
/*
|
|
* floorf(x)
|
|
* Return x rounded toward -inf to integral value
|
|
* Method:
|
|
* 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
|
|
@@ -5,18 +5,18 @@
|
|
*
|
|
* Developed at SunPro, a Sun Microsystems, Inc. business.
|
|
* Permission to use, copy, modify, and distribute this
|
|
* software is freely granted, provided that this notice
|
|
* is preserved.
|
|
* ====================================================
|
|
*/
|
|
|
|
-#include <sys/cdefs.h>
|
|
-__FBSDID("$FreeBSD$");
|
|
+//#include <sys/cdefs.h>
|
|
+//__FBSDID("$FreeBSD$");
|
|
|
|
/* double log1p(double x)
|
|
*
|
|
* Method :
|
|
* 1. Argument Reduction: find k and f such that
|
|
* 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
|
|
@@ -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)
|
|
* HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
|
|
* 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
|
|
* SUCH DAMAGE.
|
|
*/
|
|
|
|
-#include <sys/cdefs.h>
|
|
-__FBSDID("$FreeBSD$");
|
|
+//#include <sys/cdefs.h>
|
|
+//__FBSDID("$FreeBSD$");
|
|
|
|
#include <fenv.h>
|
|
#include "math_private.h"
|
|
|
|
/*
|
|
* We save and restore the floating-point environment to avoid raising
|
|
* an inexact exception. We can get away with using fesetenv()
|
|
* instead of feclearexcept()/feupdateenv() to restore the environment
|
|
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
|
|
@@ -5,18 +5,18 @@
|
|
*
|
|
* Developed at SunPro, a Sun Microsystems, Inc. business.
|
|
* Permission to use, copy, modify, and distribute this
|
|
* software is freely granted, provided that this notice
|
|
* is preserved.
|
|
* ====================================================
|
|
*/
|
|
|
|
-#include <sys/cdefs.h>
|
|
-__FBSDID("$FreeBSD$");
|
|
+//#include <sys/cdefs.h>
|
|
+//__FBSDID("$FreeBSD$");
|
|
|
|
/*
|
|
* rint(x)
|
|
* Return x rounded to integral value according to the prevailing
|
|
* rounding mode.
|
|
* 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
|
|
@@ -8,18 +8,18 @@
|
|
*
|
|
* Developed at SunPro, a Sun Microsystems, Inc. business.
|
|
* Permission to use, copy, modify, and distribute this
|
|
* software is freely granted, provided that this notice
|
|
* is preserved.
|
|
* ====================================================
|
|
*/
|
|
|
|
-#include <sys/cdefs.h>
|
|
-__FBSDID("$FreeBSD$");
|
|
+//#include <sys/cdefs.h>
|
|
+//__FBSDID("$FreeBSD$");
|
|
|
|
#include <float.h>
|
|
#include <stdint.h>
|
|
|
|
#include "math_private.h"
|
|
|
|
static const float
|
|
TWO23[2]={
|
|
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
|
|
@@ -5,18 +5,18 @@
|
|
*
|
|
* Developed at SunPro, a Sun Microsystems, Inc. business.
|
|
* Permission to use, copy, modify, and distribute this
|
|
* software is freely granted, provided that this notice
|
|
* is preserved.
|
|
* ====================================================
|
|
*/
|
|
|
|
-#include <sys/cdefs.h>
|
|
-__FBSDID("$FreeBSD$");
|
|
+//#include <sys/cdefs.h>
|
|
+//__FBSDID("$FreeBSD$");
|
|
|
|
/*
|
|
* scalbn (double x, int n)
|
|
* scalbn(x,n) returns x* 2**n computed by exponent
|
|
* manipulation rather than by actually performing an
|
|
* exponentiation or a multiplication.
|
|
*/
|
|
|
|
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
|
|
@@ -5,18 +5,18 @@
|
|
*
|
|
* Developed at SunPro, a Sun Microsystems, Inc. business.
|
|
* Permission to use, copy, modify, and distribute this
|
|
* software is freely granted, provided that this notice
|
|
* is preserved.
|
|
* ====================================================
|
|
*/
|
|
|
|
-#include <sys/cdefs.h>
|
|
-__FBSDID("$FreeBSD$");
|
|
+//#include <sys/cdefs.h>
|
|
+//__FBSDID("$FreeBSD$");
|
|
|
|
/* Tanh(x)
|
|
* Return the Hyperbolic Tangent of x
|
|
*
|
|
* Method :
|
|
* 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
|
|
@@ -5,18 +5,18 @@
|
|
*
|
|
* Developed at SunPro, a Sun Microsystems, Inc. business.
|
|
* Permission to use, copy, modify, and distribute this
|
|
* software is freely granted, provided that this notice
|
|
* is preserved.
|
|
* ====================================================
|
|
*/
|
|
|
|
-#include <sys/cdefs.h>
|
|
-__FBSDID("$FreeBSD$");
|
|
+//#include <sys/cdefs.h>
|
|
+//__FBSDID("$FreeBSD$");
|
|
|
|
/*
|
|
* trunc(x)
|
|
* Return x rounded toward 0 to integral value
|
|
* Method:
|
|
* 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
|
|
@@ -5,18 +5,18 @@
|
|
*
|
|
* Developed at SunPro, a Sun Microsystems, Inc. business.
|
|
* Permission to use, copy, modify, and distribute this
|
|
* software is freely granted, provided that this notice
|
|
* is preserved.
|
|
* ====================================================
|
|
*/
|
|
|
|
-#include <sys/cdefs.h>
|
|
-__FBSDID("$FreeBSD$");
|
|
+//#include <sys/cdefs.h>
|
|
+//__FBSDID("$FreeBSD$");
|
|
|
|
/*
|
|
* truncf(x)
|
|
* Return x rounded toward 0 to integral value
|
|
* Method:
|
|
* Bit twiddling.
|
|
* Exception:
|
|
* Inexact flag raised if x not equal to truncf(x).
|