Put external headers before internal ones. While there, replace a lot of occurences of "openlibm.h" with <openlibm.h>. It should be thought of as an external header, as it is installed along with the library.
@@ -32,7 +32,7 @@
#ifdef _WIN32
#define __fenv_static
#endif
-#include "openlibm_fenv.h"
+#include <openlibm_fenv.h>
#ifdef __GNUC_GNU_INLINE__
#error "This file must be compiled with C99 'inline' semantics"
@@ -23,8 +23,8 @@
*/
#include <float.h>
+#include <openlibm.h>
-#include "openlibm.h"
#include "math_private.h"
#include "fpmath.h"
@@ -18,7 +18,8 @@
* ld128 version of k_tan.c. See ../src/k_tan.c for most comments.
+
/*
@@ -28,10 +28,10 @@
//__FBSDID("$FreeBSD: src/lib/msun/ld128/s_exp2l.c,v 1.3 2008/02/13 10:44:44 bde Exp $");
#include <stdint.h>
#define TBLBITS 7
@@ -26,7 +26,7 @@
* $FreeBSD: src/lib/msun/ld128/s_nanl.c,v 1.3 2008/03/02 20:16:55 das Exp $
@@ -31,7 +31,8 @@
* only coshl(0)=1 is exact for finite x.
static const long double one = 1.0, half=0.5, huge = 1.0e4900L;
@@ -23,10 +23,9 @@
#define BIAS (LDBL_MAX_EXP - 1)
* ld80 version of k_tan.c. See ../src/k_tan.c for most comments.
@@ -32,8 +32,8 @@
#include "amd64/bsd_ieeefp.h"
@@ -26,9 +26,8 @@
* $FreeBSD: src/lib/msun/ld80/s_nanl.c,v 1.2 2007/12/18 23:46:31 das Exp $
DLLEXPORT long double
@@ -24,7 +24,7 @@
//#include <machine/ieee.h>
@@ -39,8 +39,8 @@
static const double
@@ -16,7 +16,8 @@
#include "cdefs-compat.h"
//__FBSDID("$FreeBSD: src/lib/msun/src/e_acosf.c,v 1.11 2008/08/03 17:39:54 das Exp $");
static const float
@@ -29,7 +29,8 @@
* acosh(NaN) is NaN without signal.
//__FBSDID("$FreeBSD: src/lib/msun/src/e_acoshf.c,v 1.8 2008/02/22 02:30:34 das Exp $");
@@ -21,9 +21,9 @@
#include "invtrig.h"
static const long double
@@ -45,8 +45,8 @@
//__FBSDID("$FreeBSD: src/lib/msun/src/e_asinf.c,v 1.13 2008/08/08 00:21:27 das Exp $");
@@ -43,8 +43,8 @@
static volatile double
//__FBSDID("$FreeBSD: src/lib/msun/src/e_atan2f.c,v 1.12 2008/08/03 17:39:54 das Exp $");
static volatile float
@@ -22,9 +22,9 @@
static volatile long double
@@ -33,7 +33,8 @@
*
static const double one = 1.0, huge = 1e300;
//__FBSDID("$FreeBSD: src/lib/msun/src/e_atanhf.c,v 1.7 2008/02/22 02:30:34 das Exp $");
static const float one = 1.0, huge = 1e30;
@@ -35,7 +35,8 @@
* only cosh(0)=1 is exact for finite x.
static const double one = 1.0, half=0.5, huge = 1.0e300;
//__FBSDID("$FreeBSD: src/lib/msun/src/e_coshf.c,v 1.9 2011/10/21 06:28:47 das Exp $");
static const float one = 1.0, half=0.5, huge = 1.0e30;
@@ -77,8 +77,8 @@
@@ -17,8 +17,8 @@
//__FBSDID("$FreeBSD: src/lib/msun/src/e_expf.c,v 1.16 2011/10/21 06:26:38 das Exp $");
@@ -20,7 +20,8 @@
* Method: shift and subtract
static const double one = 1.0, Zero[] = {0.0, -0.0,};
@@ -22,7 +22,8 @@
static const float one = 1.0, Zero[] = {0.0, -0.0,};
@@ -14,10 +14,11 @@
//__FBSDID("$FreeBSD: src/lib/msun/src/e_fmodl.c,v 1.2 2008/07/31 20:09:47 das Exp $");
@@ -21,7 +21,8 @@
* Method: call __ieee754_gamma_r
extern int signgam;
* Method: See __ieee754_lgamma_r
DLLEXPORT double
* Method: call __ieee754_gammaf_r
@@ -23,7 +23,8 @@
* Method: See __ieee754_lgammaf_r
DLLEXPORT float
@@ -47,8 +47,8 @@
//__FBSDID("$FreeBSD: src/lib/msun/src/e_hypotf.c,v 1.14 2011/10/15 07:00:28 das Exp $");
@@ -16,9 +16,9 @@
/* long double version of hypot(). See e_hypot.c for most comments. */
#define GET_LDBL_MAN(h, l, v) do { \
@@ -61,7 +61,8 @@
* 3. Special cases: y0(0)=-inf, y0(x<0)=NaN, y0(inf)=0.
static double pzero(double), qzero(double);
//__FBSDID("$FreeBSD: src/lib/msun/src/e_j0f.c,v 1.8 2008/02/22 02:30:35 das Exp $");
static float pzerof(float), qzerof(float);
* by method mentioned above.
static double pone(double), qone(double);
//__FBSDID("$FreeBSD: src/lib/msun/src/e_j1f.c,v 1.8 2008/02/22 02:30:35 das Exp $");
static float ponef(float), qonef(float);
@@ -40,7 +40,8 @@
//__FBSDID("$FreeBSD: src/lib/msun/src/e_jnf.c,v 1.11 2010/11/13 10:54:10 uqs Exp $");
* Method: call __ieee754_lgamma_r
@@ -83,7 +83,8 @@
* Method: call __ieee754_lgammaf_r
//__FBSDID("$FreeBSD: src/lib/msun/src/e_lgammaf_r.c,v 1.12 2011/10/15 07:00:28 das Exp $");
@@ -1,5 +1,7 @@
@@ -65,7 +65,8 @@
* to produce the hexadecimal values shown.
* in not-quite-routine extra precision.
#include "k_log.h"
* Float version of e_log10.c. See the latter for most comments.
#include "k_logf.h"
@@ -24,7 +24,8 @@
* Float version of e_log2.c. See the latter for most comments.
//__FBSDID("$FreeBSD: src/lib/msun/src/e_logf.c,v 1.11 2008/03/29 16:37:59 das Exp $");
@@ -57,7 +57,8 @@
//__FBSDID("$FreeBSD: src/lib/msun/src/e_powf.c,v 1.16 2011/10/21 06:26:07 das Exp $");
@@ -25,8 +25,8 @@
@@ -24,8 +24,8 @@
static const double zero = 0.0;
//__FBSDID("$FreeBSD: src/lib/msun/src/e_remainderf.c,v 1.8 2008/02/12 17:11:36 bde Exp $");
static const float zero = 0.0;
* should use scalbn() instead.
#ifdef _SCALB_INT
//__FBSDID("$FreeBSD: src/lib/msun/src/e_scalbf.c,v 1.13 2008/02/22 02:30:35 das Exp $");
@@ -32,7 +32,8 @@
* only sinh(0)=0 is exact for finite x.
static const double one = 1.0, shuge = 1.0e307;
//__FBSDID("$FreeBSD: src/lib/msun/src/e_sinhf.c,v 1.10 2011/10/21 06:28:47 das Exp $");
static const float one = 1.0, shuge = 1.0e37;
@@ -85,8 +85,8 @@
static const double one = 1.0, tiny=1.0e-300;
@@ -13,7 +13,8 @@
* ====================================================
static const float one = 1.0, tiny=1.0e-30;
//__FBSDID("$FreeBSD: src/lib/msun/src/e_sqrtl.c,v 1.1 2008/03/02 01:47:58 das Exp $");
/* Return (x + ulp) for normal positive x. Assumes no overflow. */
@@ -53,7 +53,8 @@
* any extra precision in w.
@@ -19,7 +19,8 @@
//__FBSDID("$FreeBSD: src/lib/msun/src/k_cosf.c,v 1.18 2009/06/03 08:16:34 ed Exp $");
/* |cos(x) - c(x)| < 2**-34.1 (~[-5.37e-11, 5.295e-11]). */
@@ -130,8 +130,8 @@
static const int init_jk[] = {3,4,4,6}; /* initial value for jk */
@@ -44,7 +44,8 @@
* sin(x) = x + (S1*x + (x *(r-y/2)+y))
//__FBSDID("$FreeBSD: src/lib/msun/src/k_sinf.c,v 1.16 2009/06/03 08:16:34 ed Exp $");
/* |sin(x)/x - s(x)| < 2**-37.5 (~[-4.89e-12, 4.824e-12]). */
@@ -49,8 +49,10 @@
* = 1 - 2*(tan(y) - (tan(y)^2)/(1+tan(y)))
static const double xxx[] = {
3.33333333333334091986e-01, /* 3FD55555, 55555563 */
1.33333333333201242699e-01, /* 3FC11111, 1110FE7A */
//__FBSDID("$FreeBSD: src/lib/msun/src/k_tanf.c,v 1.23 2009/06/03 08:16:34 ed Exp $");
/* |tan(x)/x - t(x)| < 2**-25.5 (~[-2e-08, 2e-08]). */
@@ -60,7 +60,7 @@
* := sign(x)*log1p(|x| + x^2/(1 + sqrt(1+x^2)))
//__FBSDID("$FreeBSD: src/lib/msun/src/s_asinhf.c,v 1.9 2008/02/22 02:30:35 das Exp $");
@@ -34,8 +34,8 @@
static const double atanhi[] = {
//__FBSDID("$FreeBSD: src/lib/msun/src/s_atanf.c,v 1.10 2008/08/01 01:24:25 das Exp $");
static const float atanhi[] = {
@@ -20,9 +20,9 @@
@@ -15,7 +15,8 @@
//__FBSDID("$FreeBSD: src/lib/msun/src/s_cbrt.c,v 1.17 2011/03/12 16:50:39 kargl Exp $");
/* cbrt(x)
@@ -17,7 +17,8 @@
//__FBSDID("$FreeBSD: src/lib/msun/src/s_cbrtf.c,v 1.18 2008/02/22 02:30:35 das Exp $");
/* cbrtf(x)
@@ -18,11 +18,11 @@
//__FBSDID("$FreeBSD: src/lib/msun/src/s_cbrtl.c,v 1.1 2011/03/12 19:37:35 kargl Exp $");
// VBS
//#include <ieeefp.h>
-#include "fpmath.h"
+#include "fpmath.h"
#if defined(_WIN32) && defined(__i386__)
#include "i387/bsd_ieeefp.h"
static const double huge = 1.0e300;
//__FBSDID("$FreeBSD: src/lib/msun/src/s_ceilf.c,v 1.8 2008/02/22 02:30:35 das Exp $");
static const float huge = 1.0e30;
* with the sign bit of y.
//#define INLINE_REM_PIO2
//#include "e_rem_pio2.c"
@@ -18,8 +18,8 @@
//__FBSDID("$FreeBSD: src/lib/msun/src/s_cosf.c,v 1.18 2008/02/25 22:19:17 bde Exp $");
//#define INLINE_KERNEL_COSDF
//#define INLINE_KERNEL_SINDF
//#define INLINE_REM_PIO2F
@@ -33,8 +33,8 @@
#if LDBL_MANT_DIG == 64
#include "../ld80/e_rem_pio2l.h"
@@ -107,8 +107,8 @@
* erfc/erf(NaN) is NaN
//__FBSDID("$FreeBSD: src/lib/msun/src/s_erff.c,v 1.8 2008/02/22 02:30:35 das Exp $");
@@ -28,8 +28,8 @@
//__FBSDID("$FreeBSD: src/lib/msun/src/s_exp2.c,v 1.7 2008/02/22 02:27:34 das Exp $");
#define TBLBITS 8
//__FBSDID("$FreeBSD: src/lib/msun/src/s_exp2f.c,v 1.9 2008/02/22 02:27:34 das Exp $");
#define TBLBITS 4
@@ -109,8 +109,8 @@
//__FBSDID("$FreeBSD: src/lib/msun/src/s_expm1f.c,v 1.12 2011/10/21 06:26:38 das Exp $");
@@ -14,7 +14,8 @@
* fabs(x) returns the absolute value of x.
* fabsf(x) returns the absolute value of x.
* no branching!
DLLEXPORT int