From f7c5ad70888cbdb35e5da5043d8af08b80a1754f Mon Sep 17 00:00:00 2001 From: Jameson Nash Date: Sat, 3 Aug 2013 05:15:47 -0400 Subject: [PATCH] add DLLEXPORT to files in src and `int isopenlibm() = 1` function" --- src/Make.files | 2 +- src/common.c | 4 ++++ src/e_acos.c | 2 +- src/e_acosf.c | 2 +- src/e_acosh.c | 2 +- src/e_acoshf.c | 2 +- src/e_acosl.c | 2 +- src/e_asin.c | 2 +- src/e_asinf.c | 2 +- src/e_asinl.c | 2 +- src/e_atan2.c | 2 +- src/e_atan2f.c | 2 +- src/e_atan2l.c | 2 +- src/e_atanh.c | 2 +- src/e_atanhf.c | 2 +- src/e_cosh.c | 2 +- src/e_coshf.c | 2 +- src/e_exp.c | 2 +- src/e_expf.c | 2 +- src/e_fmod.c | 2 +- src/e_fmodf.c | 2 +- src/e_fmodl.c | 2 +- src/e_gamma.c | 2 +- src/e_gamma_r.c | 2 +- src/e_gammaf.c | 2 +- src/e_gammaf_r.c | 2 +- src/e_hypot.c | 2 +- src/e_hypotf.c | 2 +- src/e_hypotl.c | 2 +- src/e_j0.c | 4 ++-- src/e_j0f.c | 4 ++-- src/e_j1.c | 4 ++-- src/e_j1f.c | 4 ++-- src/e_jn.c | 4 ++-- src/e_jnf.c | 4 ++-- src/e_lgamma.c | 2 +- src/e_lgamma_r.c | 2 +- src/e_lgammaf.c | 2 +- src/e_lgammaf_r.c | 2 +- src/e_log.c | 2 +- src/e_log10.c | 2 +- src/e_log10f.c | 2 +- src/e_log2.c | 2 +- src/e_log2f.c | 2 +- src/e_logf.c | 2 +- src/e_pow.c | 2 +- src/e_powf.c | 2 +- src/e_rem_pio2.c | 2 +- src/e_rem_pio2f.c | 2 +- src/e_remainder.c | 2 +- src/e_remainderf.c | 2 +- src/e_remainderl.c | 3 ++- src/e_scalb.c | 4 ++-- src/e_scalbf.c | 4 ++-- src/e_sinh.c | 2 +- src/e_sinhf.c | 2 +- src/e_sqrt.c | 2 +- src/e_sqrtf.c | 2 +- src/e_sqrtl.c | 3 ++- src/k_cos.c | 2 +- src/k_cosf.c | 2 +- src/k_exp.c | 4 ++-- src/k_expf.c | 4 ++-- src/k_rem_pio2.c | 2 +- src/k_sin.c | 2 +- src/k_sinf.c | 2 +- src/k_tan.c | 2 +- src/k_tanf.c | 2 +- src/math_private.h | 10 ++++++++++ src/s_asinh.c | 2 +- src/s_asinhf.c | 2 +- src/s_atan.c | 2 +- src/s_atanf.c | 2 +- src/s_atanl.c | 2 +- src/s_carg.c | 3 ++- src/s_cargf.c | 3 ++- src/s_cargl.c | 3 ++- src/s_cbrt.c | 2 +- src/s_cbrtf.c | 2 +- src/s_cbrtl.c | 2 +- src/s_ccosh.c | 4 ++-- src/s_ccoshf.c | 4 ++-- src/s_ceil.c | 2 +- src/s_ceilf.c | 2 +- src/s_ceill.c | 2 +- src/s_cexp.c | 2 +- src/s_cexpf.c | 2 +- src/s_cimag.c | 2 +- src/s_cimagf.c | 2 +- src/s_cimagl.c | 2 +- src/s_conj.c | 2 +- src/s_conjf.c | 2 +- src/s_conjl.c | 2 +- src/s_copysign.c | 2 +- src/s_copysignf.c | 2 +- src/s_copysignl.c | 4 ++-- src/s_cos.c | 2 +- src/s_cosf.c | 2 +- src/s_cosl.c | 2 +- src/s_cpow.c | 3 ++- src/s_cpowf.c | 3 ++- src/s_cpowl.c | 3 ++- src/s_cproj.c | 2 +- src/s_cprojf.c | 2 +- src/s_cprojl.c | 2 +- src/s_creal.c | 3 ++- src/s_crealf.c | 3 ++- src/s_creall.c | 3 ++- src/s_csinh.c | 4 ++-- src/s_csinhf.c | 4 ++-- src/s_csqrt.c | 2 +- src/s_csqrtf.c | 2 +- src/s_csqrtl.c | 2 +- src/s_ctanh.c | 4 ++-- src/s_ctanhf.c | 4 ++-- src/s_erf.c | 4 ++-- src/s_erff.c | 4 ++-- src/s_exp2.c | 2 +- src/s_exp2f.c | 2 +- src/s_expm1.c | 2 +- src/s_expm1f.c | 2 +- src/s_fabs.c | 2 +- src/s_fabsf.c | 2 +- src/s_fabsl.c | 4 ++-- src/s_fdim.c | 4 ++-- src/s_finite.c | 3 ++- src/s_finitef.c | 3 ++- src/s_floor.c | 2 +- src/s_floorf.c | 2 +- src/s_floorl.c | 3 ++- src/s_fma.c | 2 +- src/s_fmaf.c | 2 +- src/s_fmal.c | 4 ++-- src/s_fmax.c | 4 ++-- src/s_fmaxf.c | 4 ++-- src/s_fmaxl.c | 4 ++-- src/s_fmin.c | 4 ++-- src/s_fminf.c | 4 ++-- src/s_fminl.c | 4 ++-- src/s_fpclassify.c | 8 ++++---- src/s_frexp.c | 2 +- src/s_frexpf.c | 2 +- src/s_frexpl.c | 4 ++-- src/s_ilogb.c | 3 ++- src/s_ilogbf.c | 3 ++- src/s_ilogbl.c | 4 ++-- src/s_isfinite.c | 8 ++++---- src/s_isinf.c | 8 ++++---- src/s_isnan.c | 8 ++++---- src/s_isnormal.c | 8 ++++---- src/s_log1p.c | 2 +- src/s_log1pf.c | 2 +- src/s_logb.c | 2 +- src/s_logbf.c | 2 +- src/s_logbl.c | 4 ++-- src/s_lrint.c | 3 ++- src/s_lround.c | 3 ++- src/s_modf.c | 2 +- src/s_modff.c | 2 +- src/s_modfl.c | 3 ++- src/s_nan.c | 6 +++--- src/s_nearbyint.c | 3 ++- src/s_nextafter.c | 2 +- src/s_nextafterf.c | 2 +- src/s_nextafterl.c | 2 +- src/s_nexttoward.c | 2 +- src/s_nexttowardf.c | 2 +- src/s_remquo.c | 2 +- src/s_remquof.c | 2 +- src/s_remquol.c | 2 +- src/s_rint.c | 2 +- src/s_rintf.c | 2 +- src/s_rintl.c | 2 +- src/s_round.c | 3 ++- src/s_roundf.c | 3 ++- src/s_roundl.c | 3 ++- src/s_scalbln.c | 7 ++++--- src/s_scalbn.c | 2 +- src/s_scalbnf.c | 2 +- src/s_scalbnl.c | 4 ++-- src/s_signbit.c | 8 ++++---- src/s_significand.c | 2 +- src/s_significandf.c | 2 +- src/s_sin.c | 2 +- src/s_sincos.c | 4 ++-- src/s_sincosf.c | 6 +++--- src/s_sincosl.c | 12 ++++++------ src/s_sinf.c | 2 +- src/s_sinl.c | 2 +- src/s_tan.c | 2 +- src/s_tanf.c | 2 +- src/s_tanh.c | 2 +- src/s_tanhf.c | 2 +- src/s_tanl.c | 2 +- src/s_tgammaf.c | 3 ++- src/s_trunc.c | 2 +- src/s_truncf.c | 2 +- src/s_truncl.c | 3 ++- src/w_cabs.c | 3 ++- src/w_cabsf.c | 2 +- src/w_cabsl.c | 3 ++- src/w_drem.c | 3 ++- src/w_dremf.c | 2 +- 203 files changed, 306 insertions(+), 263 deletions(-) create mode 100644 src/common.c diff --git a/src/Make.files b/src/Make.files index d363bfe..04abe4d 100644 --- a/src/Make.files +++ b/src/Make.files @@ -1,4 +1,4 @@ -$(CUR_SRCS) = \ +$(CUR_SRCS) = common.c \ e_acos.c e_acosf.c e_acosh.c e_acoshf.c e_asin.c e_asinf.c \ e_atan2.c e_atan2f.c e_atanh.c e_atanhf.c e_cosh.c e_coshf.c e_exp.c \ e_expf.c e_fmod.c e_fmodf.c e_gamma.c e_gamma_r.c e_gammaf.c \ diff --git a/src/common.c b/src/common.c new file mode 100644 index 0000000..ce86102 --- /dev/null +++ b/src/common.c @@ -0,0 +1,4 @@ +#include "math_private.h" +DLLEXPORT int isopenlibm() { + return 1; +} diff --git a/src/e_acos.c b/src/e_acos.c index 810f5e9..3eea520 100644 --- a/src/e_acos.c +++ b/src/e_acos.c @@ -61,7 +61,7 @@ qS2 = 2.02094576023350569471e+00, /* 0x40002AE5, 0x9C598AC8 */ qS3 = -6.88283971605453293030e-01, /* 0xBFE6066C, 0x1B8D0159 */ qS4 = 7.70381505559019352791e-02; /* 0x3FB3B8C5, 0xB12E9282 */ -double +DLLEXPORT double __ieee754_acos(double x) { double z,p,q,r,w,s,c,df; diff --git a/src/e_acosf.c b/src/e_acosf.c index bfce86b..8e32fae 100644 --- a/src/e_acosf.c +++ b/src/e_acosf.c @@ -31,7 +31,7 @@ pS1 = -4.2743422091e-02, pS2 = -8.6563630030e-03, qS1 = -7.0662963390e-01; -float +DLLEXPORT float __ieee754_acosf(float x) { float z,p,q,r,w,s,c,df; diff --git a/src/e_acosh.c b/src/e_acosh.c index 3a05602..6180b24 100644 --- a/src/e_acosh.c +++ b/src/e_acosh.c @@ -36,7 +36,7 @@ static const double one = 1.0, ln2 = 6.93147180559945286227e-01; /* 0x3FE62E42, 0xFEFA39EF */ -double +DLLEXPORT double __ieee754_acosh(double x) { double t; diff --git a/src/e_acoshf.c b/src/e_acoshf.c index 753d7ce..65de3f2 100644 --- a/src/e_acoshf.c +++ b/src/e_acoshf.c @@ -23,7 +23,7 @@ static const float one = 1.0, ln2 = 6.9314718246e-01; /* 0x3f317218 */ -float +DLLEXPORT float __ieee754_acoshf(float x) { float t; diff --git a/src/e_acosl.c b/src/e_acosl.c index 015e37b..fe6270f 100644 --- a/src/e_acosl.c +++ b/src/e_acosl.c @@ -40,7 +40,7 @@ static const long double pi = 3.14159265358979323846264338327950280e+00L; #endif -long double +DLLEXPORT long double acosl(long double x) { union IEEEl2bits u; diff --git a/src/e_asin.c b/src/e_asin.c index 3b01a19..a0e76f5 100644 --- a/src/e_asin.c +++ b/src/e_asin.c @@ -67,7 +67,7 @@ qS2 = 2.02094576023350569471e+00, /* 0x40002AE5, 0x9C598AC8 */ qS3 = -6.88283971605453293030e-01, /* 0xBFE6066C, 0x1B8D0159 */ qS4 = 7.70381505559019352791e-02; /* 0x3FB3B8C5, 0xB12E9282 */ -double +DLLEXPORT double __ieee754_asin(double x) { double t=0.0,w,p,q,c,r,s; diff --git a/src/e_asinf.c b/src/e_asinf.c index 679ce98..4db4550 100644 --- a/src/e_asinf.c +++ b/src/e_asinf.c @@ -31,7 +31,7 @@ qS1 = -7.0662963390e-01; static const double pio2 = 1.570796326794896558e+00; -float +DLLEXPORT float __ieee754_asinf(float x) { double s; diff --git a/src/e_asinl.c b/src/e_asinl.c index 2deef4b..6de8f0a 100644 --- a/src/e_asinl.c +++ b/src/e_asinl.c @@ -30,7 +30,7 @@ static const long double one = 1.00000000000000000000e+00, huge = 1.000e+300; -long double +DLLEXPORT long double asinl(long double x) { union IEEEl2bits u; diff --git a/src/e_atan2.c b/src/e_atan2.c index 6a2a453..d676c7e 100644 --- a/src/e_atan2.c +++ b/src/e_atan2.c @@ -57,7 +57,7 @@ pi = 3.1415926535897931160E+00; /* 0x400921FB, 0x54442D18 */ static volatile double pi_lo = 1.2246467991473531772E-16; /* 0x3CA1A626, 0x33145C07 */ -double +DLLEXPORT double __ieee754_atan2(double y, double x) { double z; diff --git a/src/e_atan2f.c b/src/e_atan2f.c index 749b7ef..2dfc31f 100644 --- a/src/e_atan2f.c +++ b/src/e_atan2f.c @@ -29,7 +29,7 @@ pi = 3.1415927410e+00; /* 0x40490fdb */ static volatile float pi_lo = -8.7422776573e-08; /* 0xb3bbbd2e */ -float +DLLEXPORT float __ieee754_atan2f(float y, float x) { float z; diff --git a/src/e_atan2l.c b/src/e_atan2l.c index 0baac0e..3bad12f 100644 --- a/src/e_atan2l.c +++ b/src/e_atan2l.c @@ -43,7 +43,7 @@ static const long double pi = 3.14159265358979323846264338327950280e+00L; #endif -long double +DLLEXPORT long double atan2l(long double y, long double x) { union IEEEl2bits ux, uy; diff --git a/src/e_atanh.c b/src/e_atanh.c index eb2d1ed..d09ba50 100644 --- a/src/e_atanh.c +++ b/src/e_atanh.c @@ -39,7 +39,7 @@ static const double one = 1.0, huge = 1e300; static const double zero = 0.0; -double +DLLEXPORT double __ieee754_atanh(double x) { double t; diff --git a/src/e_atanhf.c b/src/e_atanhf.c index 4f33cf5..e69bf48 100644 --- a/src/e_atanhf.c +++ b/src/e_atanhf.c @@ -23,7 +23,7 @@ static const float one = 1.0, huge = 1e30; static const float zero = 0.0; -float +DLLEXPORT float __ieee754_atanhf(float x) { float t; diff --git a/src/e_cosh.c b/src/e_cosh.c index dd495b4..a8a5c6c 100644 --- a/src/e_cosh.c +++ b/src/e_cosh.c @@ -40,7 +40,7 @@ static const double one = 1.0, half=0.5, huge = 1.0e300; -double +DLLEXPORT double __ieee754_cosh(double x) { double t,w; diff --git a/src/e_coshf.c b/src/e_coshf.c index 1b34509..f5e8202 100644 --- a/src/e_coshf.c +++ b/src/e_coshf.c @@ -21,7 +21,7 @@ static const float one = 1.0, half=0.5, huge = 1.0e30; -float +DLLEXPORT float __ieee754_coshf(float x) { float t,w; diff --git a/src/e_exp.c b/src/e_exp.c index 1097c44..2ad215e 100644 --- a/src/e_exp.c +++ b/src/e_exp.c @@ -101,7 +101,7 @@ P5 = 4.13813679705723846039e-08; /* 0x3E663769, 0x72BEA4D0 */ static volatile double twom1000= 9.33263618503218878990e-302; /* 2**-1000=0x01700000,0*/ -double +DLLEXPORT double __ieee754_exp(double x) /* default IEEE double exp */ { double y,hi=0.0,lo=0.0,c,t,twopk; diff --git a/src/e_expf.c b/src/e_expf.c index ca1629c..15d622b 100644 --- a/src/e_expf.c +++ b/src/e_expf.c @@ -41,7 +41,7 @@ P2 = -2.7667332906e-3; /* -0xb55215.0p-32 */ static volatile float twom100 = 7.8886090522e-31; /* 2**-100=0x0d800000 */ -float +DLLEXPORT float __ieee754_expf(float x) { float y,hi=0.0,lo=0.0,c,t,twopk; diff --git a/src/e_fmod.c b/src/e_fmod.c index 5633776..83076fb 100644 --- a/src/e_fmod.c +++ b/src/e_fmod.c @@ -25,7 +25,7 @@ static const double one = 1.0, Zero[] = {0.0, -0.0,}; -double +DLLEXPORT double __ieee754_fmod(double x, double y) { int32_t n,hx,hy,hz,ix,iy,sx,i; diff --git a/src/e_fmodf.c b/src/e_fmodf.c index aed7852..14ee260 100644 --- a/src/e_fmodf.c +++ b/src/e_fmodf.c @@ -27,7 +27,7 @@ static const float one = 1.0, Zero[] = {0.0, -0.0,}; -float +DLLEXPORT float __ieee754_fmodf(float x, float y) { int32_t n,hx,hy,hz,ix,iy,sx,i; diff --git a/src/e_fmodl.c b/src/e_fmodl.c index 8f4c5b2..452d155 100644 --- a/src/e_fmodl.c +++ b/src/e_fmodl.c @@ -61,7 +61,7 @@ static const long double one = 1.0, Zero[] = {0.0, -0.0,}; * - The high part of the mantissa fits in an int64_t with enough room * for an explicit integer bit in front of the fractional bits. */ -long double +DLLEXPORT long double fmodl(long double x, long double y) { union IEEEl2bits ux, uy; diff --git a/src/e_gamma.c b/src/e_gamma.c index 26cb25e..f256a15 100644 --- a/src/e_gamma.c +++ b/src/e_gamma.c @@ -26,7 +26,7 @@ extern int signgam; -double +DLLEXPORT double __ieee754_gamma(double x) { return __ieee754_gamma_r(x,&signgam); diff --git a/src/e_gamma_r.c b/src/e_gamma_r.c index b45b69e..3d4e767 100644 --- a/src/e_gamma_r.c +++ b/src/e_gamma_r.c @@ -25,7 +25,7 @@ #include "openlibm.h" #include "math_private.h" -double +DLLEXPORT double __ieee754_gamma_r(double x, int *signgamp) { return __ieee754_lgamma_r(x,signgamp); diff --git a/src/e_gammaf.c b/src/e_gammaf.c index 05f88d8..7350ab3 100644 --- a/src/e_gammaf.c +++ b/src/e_gammaf.c @@ -27,7 +27,7 @@ extern int signgam; -float +DLLEXPORT float __ieee754_gammaf(float x) { return __ieee754_gammaf_r(x,&signgam); diff --git a/src/e_gammaf_r.c b/src/e_gammaf_r.c index f22df40..68fefd2 100644 --- a/src/e_gammaf_r.c +++ b/src/e_gammaf_r.c @@ -26,7 +26,7 @@ #include "openlibm.h" #include "math_private.h" -float +DLLEXPORT float __ieee754_gammaf_r(float x, int *signgamp) { return __ieee754_lgammaf_r(x,signgamp); diff --git a/src/e_hypot.c b/src/e_hypot.c index db6d0cf..db1745f 100644 --- a/src/e_hypot.c +++ b/src/e_hypot.c @@ -51,7 +51,7 @@ #include "openlibm.h" #include "math_private.h" -double +DLLEXPORT double __ieee754_hypot(double x, double y) { double a,b,t1,t2,y1,y2,w; diff --git a/src/e_hypotf.c b/src/e_hypotf.c index 531f1db..7c2885a 100644 --- a/src/e_hypotf.c +++ b/src/e_hypotf.c @@ -19,7 +19,7 @@ #include "openlibm.h" #include "math_private.h" -float +DLLEXPORT float __ieee754_hypotf(float x, float y) { float a,b,t1,t2,y1,y2,w; diff --git a/src/e_hypotl.c b/src/e_hypotl.c index f7a35ad..269233b 100644 --- a/src/e_hypotl.c +++ b/src/e_hypotl.c @@ -45,7 +45,7 @@ typedef u_int64_t man_t; typedef u_int32_t man_t; #endif -long double +DLLEXPORT long double hypotl(long double x, long double y) { long double a=x,b=y,t1,t2,y1,y2,w; diff --git a/src/e_j0.c b/src/e_j0.c index 8a581a7..aa757bb 100644 --- a/src/e_j0.c +++ b/src/e_j0.c @@ -81,7 +81,7 @@ S04 = 1.16614003333790000205e-09; /* 0x3E1408BC, 0xF4745D8F */ static const double zero = 0.0; -double +DLLEXPORT double __ieee754_j0(double x) { double z, s,c,ss,cc,r,u,v; @@ -142,7 +142,7 @@ v02 = 7.60068627350353253702e-05, /* 0x3F13ECBB, 0xF578C6C1 */ v03 = 2.59150851840457805467e-07, /* 0x3E91642D, 0x7FF202FD */ v04 = 4.41110311332675467403e-10; /* 0x3DFE5018, 0x3BD6D9EF */ -double +DLLEXPORT double __ieee754_y0(double x) { double z, s,c,ss,cc,u,v; diff --git a/src/e_j0f.c b/src/e_j0f.c index c3edd2e..61b59e6 100644 --- a/src/e_j0f.c +++ b/src/e_j0f.c @@ -38,7 +38,7 @@ S04 = 1.1661400734e-09; /* 0x30a045e8 */ static const float zero = 0.0; -float +DLLEXPORT float __ieee754_j0f(float x) { float z, s,c,ss,cc,r,u,v; @@ -99,7 +99,7 @@ v02 = 7.6006865129e-05, /* 0x389f65e0 */ v03 = 2.5915085189e-07, /* 0x348b216c */ v04 = 4.4111031494e-10; /* 0x2ff280c2 */ -float +DLLEXPORT float __ieee754_y0f(float x) { float z, s,c,ss,cc,u,v; diff --git a/src/e_j1.c b/src/e_j1.c index 925d498..dfb8f78 100644 --- a/src/e_j1.c +++ b/src/e_j1.c @@ -82,7 +82,7 @@ s05 = 1.23542274426137913908e-11; /* 0x3DAB2ACF, 0xCFB97ED8 */ static const double zero = 0.0; -double +DLLEXPORT double __ieee754_j1(double x) { double z, s,c,ss,cc,r,u,v,y; @@ -139,7 +139,7 @@ static const double V0[5] = { 1.66559246207992079114e-11, /* 0x3DB25039, 0xDACA772A */ }; -double +DLLEXPORT double __ieee754_y1(double x) { double z, s,c,ss,cc,u,v; diff --git a/src/e_j1f.c b/src/e_j1f.c index 47631e9..aeba0f3 100644 --- a/src/e_j1f.c +++ b/src/e_j1f.c @@ -39,7 +39,7 @@ s05 = 1.2354227016e-11; /* 0x2d59567e */ static const float zero = 0.0; -float +DLLEXPORT float __ieee754_j1f(float x) { float z, s,c,ss,cc,r,u,v,y; @@ -96,7 +96,7 @@ static const float V0[5] = { 1.6655924903e-11, /* 0x2d9281cf */ }; -float +DLLEXPORT float __ieee754_y1f(float x) { float z, s,c,ss,cc,u,v; diff --git a/src/e_jn.c b/src/e_jn.c index dd0d2ac..883e45c 100644 --- a/src/e_jn.c +++ b/src/e_jn.c @@ -50,7 +50,7 @@ one = 1.00000000000000000000e+00; /* 0x3FF00000, 0x00000000 */ static const double zero = 0.00000000000000000000e+00; -double +DLLEXPORT double __ieee754_jn(int n, double x) { int32_t i,hx,ix,lx, sgn; @@ -211,7 +211,7 @@ __ieee754_jn(int n, double x) if(sgn==1) return -b; else return b; } -double +DLLEXPORT double __ieee754_yn(int n, double x) { int32_t i,hx,ix,lx; diff --git a/src/e_jnf.c b/src/e_jnf.c index e6e6e98..30853e4 100644 --- a/src/e_jnf.c +++ b/src/e_jnf.c @@ -25,7 +25,7 @@ one = 1.0000000000e+00; /* 0x3F800000 */ static const float zero = 0.0000000000e+00; -float +DLLEXPORT float __ieee754_jnf(int n, float x) { int32_t i,hx,ix, sgn; @@ -163,7 +163,7 @@ __ieee754_jnf(int n, float x) if(sgn==1) return -b; else return b; } -float +DLLEXPORT float __ieee754_ynf(int n, float x) { int32_t i,hx,ix,ib; diff --git a/src/e_lgamma.c b/src/e_lgamma.c index 4a2e880..469b38f 100644 --- a/src/e_lgamma.c +++ b/src/e_lgamma.c @@ -26,7 +26,7 @@ extern int signgam; -double +DLLEXPORT double __ieee754_lgamma(double x) { return __ieee754_lgamma_r(x,&signgam); diff --git a/src/e_lgamma_r.c b/src/e_lgamma_r.c index 1a4550d..ad2c2ed 100644 --- a/src/e_lgamma_r.c +++ b/src/e_lgamma_r.c @@ -201,7 +201,7 @@ static const double zero= 0.00000000000000000000e+00; } -double +DLLEXPORT double __ieee754_lgamma_r(double x, int *signgamp) { double t,y,z,nadj,p,p1,p2,p3,q,r,w; diff --git a/src/e_lgammaf.c b/src/e_lgammaf.c index a3e5c1b..8ac900c 100644 --- a/src/e_lgammaf.c +++ b/src/e_lgammaf.c @@ -27,7 +27,7 @@ extern int signgam; -float +DLLEXPORT float __ieee754_lgammaf(float x) { return __ieee754_lgammaf_r(x,&signgam); diff --git a/src/e_lgammaf_r.c b/src/e_lgammaf_r.c index 70492c6..5d9fc87 100644 --- a/src/e_lgammaf_r.c +++ b/src/e_lgammaf_r.c @@ -134,7 +134,7 @@ static const float zero= 0.0000000000e+00; } -float +DLLEXPORT float __ieee754_lgammaf_r(float x, int *signgamp) { float t,y,z,nadj,p,p1,p2,p3,q,r,w; diff --git a/src/e_log.c b/src/e_log.c index 81c3f2b..ee73d74 100644 --- a/src/e_log.c +++ b/src/e_log.c @@ -82,7 +82,7 @@ Lg7 = 1.479819860511658591e-01; /* 3FC2F112 DF3E5244 */ static const double zero = 0.0; -double +DLLEXPORT double __ieee754_log(double x) { double hfsq,f,s,z,R,w,t1,t2,dk; diff --git a/src/e_log10.c b/src/e_log10.c index f5fda16..a6f9dd6 100644 --- a/src/e_log10.c +++ b/src/e_log10.c @@ -35,7 +35,7 @@ log10_2lo = 3.69423907715893078616e-13; /* 0x3D59FEF3, 0x11F12B36 */ static const double zero = 0.0; -double +DLLEXPORT double __ieee754_log10(double x) { double f,hfsq,hi,lo,r,val_hi,val_lo,w,y,y2; diff --git a/src/e_log10f.c b/src/e_log10f.c index d628260..9ebce9a 100644 --- a/src/e_log10f.c +++ b/src/e_log10f.c @@ -32,7 +32,7 @@ log10_2lo = 7.9034151668e-07; /* 0x355427db */ static const float zero = 0.0; -float +DLLEXPORT float __ieee754_log10f(float x) { float f,hfsq,hi,lo,r,y; diff --git a/src/e_log2.c b/src/e_log2.c index cfc2fd1..0784679 100644 --- a/src/e_log2.c +++ b/src/e_log2.c @@ -35,7 +35,7 @@ ivln2lo = 1.67517131648865118353e-10; /* 0x3de705fc, 0x2eefa200 */ static const double zero = 0.0; -double +DLLEXPORT double __ieee754_log2(double x) { double f,hfsq,hi,lo,r,val_hi,val_lo,w,y; diff --git a/src/e_log2f.c b/src/e_log2f.c index 173efaf..99fa397 100644 --- a/src/e_log2f.c +++ b/src/e_log2f.c @@ -30,7 +30,7 @@ ivln2lo = -1.7605285393e-04; /* 0xb9389ad4 */ static const float zero = 0.0; -float +DLLEXPORT float __ieee754_log2f(float x) { float f,hfsq,hi,lo,r,y; diff --git a/src/e_logf.c b/src/e_logf.c index 8e88abd..c360cb7 100644 --- a/src/e_logf.c +++ b/src/e_logf.c @@ -31,7 +31,7 @@ Lg4 = 0xf89e26.0p-26; /* 0.24279078841 */ static const float zero = 0.0; -float +DLLEXPORT float __ieee754_logf(float x) { float hfsq,f,s,z,R,w,t1,t2,dk; diff --git a/src/e_pow.c b/src/e_pow.c index 863a846..f1a0c32 100644 --- a/src/e_pow.c +++ b/src/e_pow.c @@ -93,7 +93,7 @@ ivln2 = 1.44269504088896338700e+00, /* 0x3FF71547, 0x652B82FE =1/ln2 */ ivln2_h = 1.44269502162933349609e+00, /* 0x3FF71547, 0x60000000 =24b 1/ln2*/ ivln2_l = 1.92596299112661746887e-08; /* 0x3E54AE0B, 0xF85DDF44 =1/ln2 tail*/ -double +DLLEXPORT double __ieee754_pow(double x, double y) { double z,ax,z_h,z_l,p_h,p_l; diff --git a/src/e_powf.c b/src/e_powf.c index 08886a0..c1e9e8b 100644 --- a/src/e_powf.c +++ b/src/e_powf.c @@ -52,7 +52,7 @@ ivln2 = 1.4426950216e+00, /* 0x3fb8aa3b =1/ln2 */ ivln2_h = 1.4426879883e+00, /* 0x3fb8aa00 =16b 1/ln2*/ ivln2_l = 7.0526075433e-06; /* 0x36eca570 =1/ln2 tail*/ -float +DLLEXPORT float __ieee754_powf(float x, float y) { float z,ax,z_h,z_l,p_h,p_l; diff --git a/src/e_rem_pio2.c b/src/e_rem_pio2.c index 02687cc..2c6aefc 100644 --- a/src/e_rem_pio2.c +++ b/src/e_rem_pio2.c @@ -52,7 +52,7 @@ pio2_3t = 8.47842766036889956997e-32; /* 0x397B839A, 0x252049C1 */ extern #endif //__inline int -int +DLLEXPORT int __ieee754_rem_pio2(double x, double *y) { double z,w,t,r,fn; diff --git a/src/e_rem_pio2f.c b/src/e_rem_pio2f.c index 8fe2790..9f2dc44 100644 --- a/src/e_rem_pio2f.c +++ b/src/e_rem_pio2f.c @@ -44,7 +44,7 @@ pio2_1t = 1.58932547735281966916e-08; /* 0x3E5110b4, 0x611A6263 */ extern #endif //__inline int -int +DLLEXPORT int __ieee754_rem_pio2f(float x, double *y) { double w,r,fn; diff --git a/src/e_remainder.c b/src/e_remainder.c index 70356d5..1b1410f 100644 --- a/src/e_remainder.c +++ b/src/e_remainder.c @@ -31,7 +31,7 @@ static const double zero = 0.0; -double +DLLEXPORT double __ieee754_remainder(double x, double p) { int32_t hx,hp; diff --git a/src/e_remainderf.c b/src/e_remainderf.c index b2ca56c..867c40a 100644 --- a/src/e_remainderf.c +++ b/src/e_remainderf.c @@ -22,7 +22,7 @@ static const float zero = 0.0; -float +DLLEXPORT float __ieee754_remainderf(float x, float p) { int32_t hx,hp; diff --git a/src/e_remainderl.c b/src/e_remainderl.c index ddfe9cf..af92ff5 100644 --- a/src/e_remainderl.c +++ b/src/e_remainderl.c @@ -28,8 +28,9 @@ //__FBSDID("$FreeBSD: src/lib/msun/src/e_remainderl.c,v 1.1 2008/03/30 20:47:42 das Exp $"); #include +#include "math_private.h" -long double +DLLEXPORT long double remainderl(long double x, long double y) { int quo; diff --git a/src/e_scalb.c b/src/e_scalb.c index e78baf6..d575100 100644 --- a/src/e_scalb.c +++ b/src/e_scalb.c @@ -24,10 +24,10 @@ #include "math_private.h" #ifdef _SCALB_INT -double +DLLEXPORT double __ieee754_scalb(double x, int fn) #else -double +DLLEXPORT double __ieee754_scalb(double x, double fn) #endif { diff --git a/src/e_scalbf.c b/src/e_scalbf.c index 73ebcb6..d506719 100644 --- a/src/e_scalbf.c +++ b/src/e_scalbf.c @@ -20,10 +20,10 @@ #include "math_private.h" #ifdef _SCALB_INT -float +DLLEXPORT float __ieee754_scalbf(float x, int fn) #else -float +DLLEXPORT float __ieee754_scalbf(float x, float fn) #endif { diff --git a/src/e_sinh.c b/src/e_sinh.c index 2d8b5c2..5fc451e 100644 --- a/src/e_sinh.c +++ b/src/e_sinh.c @@ -37,7 +37,7 @@ static const double one = 1.0, shuge = 1.0e307; -double +DLLEXPORT double __ieee754_sinh(double x) { double t,h; diff --git a/src/e_sinhf.c b/src/e_sinhf.c index 4dbbdb1..9a97834 100644 --- a/src/e_sinhf.c +++ b/src/e_sinhf.c @@ -21,7 +21,7 @@ static const float one = 1.0, shuge = 1.0e37; -float +DLLEXPORT float __ieee754_sinhf(float x) { float t,h; diff --git a/src/e_sqrt.c b/src/e_sqrt.c index ca2294e..df6eec4 100644 --- a/src/e_sqrt.c +++ b/src/e_sqrt.c @@ -91,7 +91,7 @@ static const double one = 1.0, tiny=1.0e-300; -double +DLLEXPORT double __ieee754_sqrt(double x) { double z; diff --git a/src/e_sqrtf.c b/src/e_sqrtf.c index e6afa9f..e842efe 100644 --- a/src/e_sqrtf.c +++ b/src/e_sqrtf.c @@ -18,7 +18,7 @@ static const float one = 1.0, tiny=1.0e-30; -float +DLLEXPORT float __ieee754_sqrtf(float x) { float z; diff --git a/src/e_sqrtl.c b/src/e_sqrtl.c index ff1ced7..7b2ea83 100644 --- a/src/e_sqrtl.c +++ b/src/e_sqrtl.c @@ -32,6 +32,7 @@ #include "fpmath.h" #include "openlibm.h" +#include "math_private.h" /* Return (x + ulp) for normal positive x. Assumes no overflow. */ static inline long double @@ -74,7 +75,7 @@ dec(long double x) * if possible. */ -long double +DLLEXPORT long double sqrtl(long double x) { union IEEEl2bits u; diff --git a/src/k_cos.c b/src/k_cos.c index eabd230..b468594 100644 --- a/src/k_cos.c +++ b/src/k_cos.c @@ -65,7 +65,7 @@ C4 = -2.75573143513906633035e-07, /* 0xBE927E4F, 0x809C52AD */ C5 = 2.08757232129817482790e-09, /* 0x3E21EE9E, 0xBDB4B1C4 */ C6 = -1.13596475577881948265e-11; /* 0xBDA8FAE9, 0xBE8838D4 */ -double +DLLEXPORT double __kernel_cos(double x, double y) { double hz,z,r,w; diff --git a/src/k_cosf.c b/src/k_cosf.c index f4b8fc7..6daeb65 100644 --- a/src/k_cosf.c +++ b/src/k_cosf.c @@ -34,7 +34,7 @@ C3 = 0x199342e0ee5069.0p-68; /* 0.0000243904487962774090654 */ extern #endif //__inline float -float +DLLEXPORT float __kernel_cosdf(double x) { double r, w, z; diff --git a/src/k_exp.c b/src/k_exp.c index b512c78..8735088 100644 --- a/src/k_exp.c +++ b/src/k_exp.c @@ -71,7 +71,7 @@ __frexp_exp(double x, int *expt) * has filtered out very large x, for which overflow would be inevitable. */ -double +DLLEXPORT double __ldexp_exp(double x, int expt) { double exp_x, scale; @@ -83,7 +83,7 @@ __ldexp_exp(double x, int expt) return (exp_x * scale); } -double complex +DLLEXPORT double complex __ldexp_cexp(double complex z, int expt) { double x, y, exp_x, scale1, scale2; diff --git a/src/k_expf.c b/src/k_expf.c index f4e16ff..c94da5c 100644 --- a/src/k_expf.c +++ b/src/k_expf.c @@ -54,7 +54,7 @@ __frexp_expf(float x, int *expt) return (exp_x); } -float +DLLEXPORT float __ldexp_expf(float x, int expt) { float exp_x, scale; @@ -66,7 +66,7 @@ __ldexp_expf(float x, int expt) return (exp_x * scale); } -float complex +DLLEXPORT float complex __ldexp_cexpf(float complex z, int expt) { float x, y, exp_x, scale1, scale2; diff --git a/src/k_rem_pio2.c b/src/k_rem_pio2.c index 6ce50b9..ba6fd3c 100644 --- a/src/k_rem_pio2.c +++ b/src/k_rem_pio2.c @@ -290,7 +290,7 @@ one = 1.0, two24 = 1.67772160000000000000e+07, /* 0x41700000, 0x00000000 */ twon24 = 5.96046447753906250000e-08; /* 0x3E700000, 0x00000000 */ -int +DLLEXPORT int __kernel_rem_pio2(double *x, double *y, int e0, int nx, int prec) { int32_t jz,jx,jv,jp,jk,carry,n,iq[20],i,j,k,m,q0,ih; diff --git a/src/k_sin.c b/src/k_sin.c index 4969844..8b7c1e5 100644 --- a/src/k_sin.c +++ b/src/k_sin.c @@ -56,7 +56,7 @@ S4 = 2.75573137070700676789e-06, /* 0x3EC71DE3, 0x57B1FE7D */ S5 = -2.50507602534068634195e-08, /* 0xBE5AE5E6, 0x8A2B9CEB */ S6 = 1.58969099521155010221e-10; /* 0x3DE5D93A, 0x5ACFD57C */ -double +DLLEXPORT double __kernel_sin(double x, double y, int iy) { double z,r,v,w; diff --git a/src/k_sinf.c b/src/k_sinf.c index bf9039a..374d610 100644 --- a/src/k_sinf.c +++ b/src/k_sinf.c @@ -33,7 +33,7 @@ S4 = 0x16cd878c3b46a7.0p-71; /* 0.0000027183114939898219064 */ extern #endif //__inline float -float +DLLEXPORT float __kernel_sindf(double x) { double r, s, w, z; diff --git a/src/k_tan.c b/src/k_tan.c index d56297f..090304c 100644 --- a/src/k_tan.c +++ b/src/k_tan.c @@ -75,7 +75,7 @@ static const double xxx[] = { #define T xxx /* INDENT ON */ -double +DLLEXPORT double __kernel_tan(double x, double y, int iy) { double z, r, v, w, s; int32_t ix, hx; diff --git a/src/k_tanf.c b/src/k_tanf.c index c4c25d0..3dac110 100644 --- a/src/k_tanf.c +++ b/src/k_tanf.c @@ -36,7 +36,7 @@ T[] = { extern #endif //__inline float -float +DLLEXPORT float __kernel_tandf(double x, int iy) { double z,r,w,s,t,u; diff --git a/src/math_private.h b/src/math_private.h index 5fe8cd6..43781b4 100644 --- a/src/math_private.h +++ b/src/math_private.h @@ -442,4 +442,14 @@ long double __kernel_sinl(long double, long double, int); long double __kernel_cosl(long double, long double); long double __kernel_tanl(long double, long double, int); +#ifdef _WIN32 +# ifdef IMPORT_EXPORTS +# define DLLEXPORT __declspec(dllimport) +# else +# define DLLEXPORT __declspec(dllexport) +# endif +#else +#define DLLEXPORT __attribute__ ((visibility("default"))) +#endif + #endif /* !_MATH_PRIVATE_H_ */ diff --git a/src/s_asinh.c b/src/s_asinh.c index 36f0489..70d1002 100644 --- a/src/s_asinh.c +++ b/src/s_asinh.c @@ -32,7 +32,7 @@ one = 1.00000000000000000000e+00, /* 0x3FF00000, 0x00000000 */ ln2 = 6.93147180559945286227e-01, /* 0x3FE62E42, 0xFEFA39EF */ huge= 1.00000000000000000000e+300; -double +DLLEXPORT double asinh(double x) { double t,w; diff --git a/src/s_asinhf.c b/src/s_asinhf.c index 9a753c5..8a988df 100644 --- a/src/s_asinhf.c +++ b/src/s_asinhf.c @@ -24,7 +24,7 @@ one = 1.0000000000e+00, /* 0x3F800000 */ ln2 = 6.9314718246e-01, /* 0x3f317218 */ huge= 1.0000000000e+30; -float +DLLEXPORT float asinhf(float x) { float t,w; diff --git a/src/s_atan.c b/src/s_atan.c index 02c7a8d..60d8df8 100644 --- a/src/s_atan.c +++ b/src/s_atan.c @@ -70,7 +70,7 @@ static const double aT[] = { one = 1.0, huge = 1.0e300; -double +DLLEXPORT double atan(double x) { double w,s1,s2,z; diff --git a/src/s_atanf.c b/src/s_atanf.c index db78fc5..505c3e8 100644 --- a/src/s_atanf.c +++ b/src/s_atanf.c @@ -45,7 +45,7 @@ static const float one = 1.0, huge = 1.0e30; -float +DLLEXPORT float atanf(float x) { float w,s1,s2,z; diff --git a/src/s_atanl.c b/src/s_atanl.c index e4e8459..ee47f4e 100644 --- a/src/s_atanl.c +++ b/src/s_atanl.c @@ -29,7 +29,7 @@ static const long double one = 1.0, huge = 1.0e300; -long double +DLLEXPORT long double atanl(long double x) { union IEEEl2bits u; diff --git a/src/s_carg.c b/src/s_carg.c index 50d22b1..195feb4 100644 --- a/src/s_carg.c +++ b/src/s_carg.c @@ -29,8 +29,9 @@ #include #include +#include "math_private.h" -double +DLLEXPORT double carg(double complex z) { diff --git a/src/s_cargf.c b/src/s_cargf.c index 9ef6f1c..7592b1f 100644 --- a/src/s_cargf.c +++ b/src/s_cargf.c @@ -29,8 +29,9 @@ #include #include +#include "math_private.h" -float +DLLEXPORT float cargf(float complex z) { diff --git a/src/s_cargl.c b/src/s_cargl.c index 7b831e3..2a725c9 100644 --- a/src/s_cargl.c +++ b/src/s_cargl.c @@ -29,8 +29,9 @@ #include #include +#include "math_private.h" -long double +DLLEXPORT long double cargl(long double complex z) { diff --git a/src/s_cbrt.c b/src/s_cbrt.c index 57f13d3..2fccb92 100644 --- a/src/s_cbrt.c +++ b/src/s_cbrt.c @@ -33,7 +33,7 @@ P2 = 1.621429720105354466140, /* 0x3ff9f160, 0x4a49d6c2 */ P3 = -0.758397934778766047437, /* 0xbfe844cb, 0xbee751d9 */ P4 = 0.145996192886612446982; /* 0x3fc2b000, 0xd4e4edd7 */ -double +DLLEXPORT double cbrt(double x) { int32_t hx; diff --git a/src/s_cbrtf.c b/src/s_cbrtf.c index 9422505..ccb7957 100644 --- a/src/s_cbrtf.c +++ b/src/s_cbrtf.c @@ -27,7 +27,7 @@ static const unsigned B1 = 709958130, /* B1 = (127-127.0/3-0.03306235651)*2**23 */ B2 = 642849266; /* B2 = (127-127.0/3-24/3-0.03306235651)*2**23 */ -float +DLLEXPORT float cbrtf(float x) { double r,T; diff --git a/src/s_cbrtl.c b/src/s_cbrtl.c index e7c2d10..08bfacc 100644 --- a/src/s_cbrtl.c +++ b/src/s_cbrtl.c @@ -33,7 +33,7 @@ static const unsigned B1 = 709958130; /* B1 = (127-127.0/3-0.03306235651)*2**23 */ -long double +DLLEXPORT long double cbrtl(long double x) { union IEEEl2bits u, v; diff --git a/src/s_ccosh.c b/src/s_ccosh.c index e6ed10c..8f55f42 100644 --- a/src/s_ccosh.c +++ b/src/s_ccosh.c @@ -44,7 +44,7 @@ static const double huge = 0x1p1023; -double complex +DLLEXPORT double complex ccosh(double complex z) { double x, y, h; @@ -146,7 +146,7 @@ ccosh(double complex z) return (cpack((x * x) * (y - y), (x + x) * (y - y))); } -double complex +DLLEXPORT double complex ccos(double complex z) { diff --git a/src/s_ccoshf.c b/src/s_ccoshf.c index ca6131b..6fcf836 100644 --- a/src/s_ccoshf.c +++ b/src/s_ccoshf.c @@ -38,7 +38,7 @@ static const float huge = 0x1p127; -float complex +DLLEXPORT float complex ccoshf(float complex z) { float x, y, h; @@ -96,7 +96,7 @@ ccoshf(float complex z) return (cpackf((x * x) * (y - y), (x + x) * (y - y))); } -float complex +DLLEXPORT float complex ccosf(float complex z) { diff --git a/src/s_ceil.c b/src/s_ceil.c index 7646153..7c02043 100644 --- a/src/s_ceil.c +++ b/src/s_ceil.c @@ -29,7 +29,7 @@ static const double huge = 1.0e300; -double +DLLEXPORT double ceil(double x) { int32_t i0,i1,j0; diff --git a/src/s_ceilf.c b/src/s_ceilf.c index 10659ab..8fb20e5 100644 --- a/src/s_ceilf.c +++ b/src/s_ceilf.c @@ -21,7 +21,7 @@ static const float huge = 1.0e30; -float +DLLEXPORT float ceilf(float x) { int32_t i0,j0; diff --git a/src/s_ceill.c b/src/s_ceill.c index 2035e98..6afc2be 100644 --- a/src/s_ceill.c +++ b/src/s_ceill.c @@ -52,7 +52,7 @@ static const long double huge = 1.0e300; -long double +DLLEXPORT long double ceill(long double x) { union IEEEl2bits u = { .e = x }; diff --git a/src/s_cexp.c b/src/s_cexp.c index b1c6b91..38e6653 100644 --- a/src/s_cexp.c +++ b/src/s_cexp.c @@ -36,7 +36,7 @@ static const u_int32_t exp_ovfl = 0x40862e42, /* high bits of MAX_EXP * ln2 ~= 710 */ cexp_ovfl = 0x4096b8e4; /* (MAX_EXP - MIN_DENORM_EXP) * ln2 */ -double complex +DLLEXPORT double complex cexp(double complex z) { double x, y, exp_x; diff --git a/src/s_cexpf.c b/src/s_cexpf.c index b2453c8..90eab17 100644 --- a/src/s_cexpf.c +++ b/src/s_cexpf.c @@ -36,7 +36,7 @@ static const u_int32_t exp_ovfl = 0x42b17218, /* MAX_EXP * ln2 ~= 88.722839355 */ cexp_ovfl = 0x43400074; /* (MAX_EXP - MIN_DENORM_EXP) * ln2 */ -float complex +DLLEXPORT float complex cexpf(float complex z) { float x, y, exp_x; diff --git a/src/s_cimag.c b/src/s_cimag.c index a14afa7..3b69c14 100644 --- a/src/s_cimag.c +++ b/src/s_cimag.c @@ -29,7 +29,7 @@ #include #include "math_private.h" -double +DLLEXPORT double cimag(double complex z) { const double_complex z1 = { .f = z }; diff --git a/src/s_cimagf.c b/src/s_cimagf.c index dabb09a..1453cac 100644 --- a/src/s_cimagf.c +++ b/src/s_cimagf.c @@ -29,7 +29,7 @@ #include #include "math_private.h" -float +DLLEXPORT float cimagf(float complex z) { const float_complex z1 = { .f = z }; diff --git a/src/s_cimagl.c b/src/s_cimagl.c index 89f4bc3..5f8876c 100644 --- a/src/s_cimagl.c +++ b/src/s_cimagl.c @@ -29,7 +29,7 @@ #include #include "math_private.h" -long double +DLLEXPORT long double cimagl(long double complex z) { const long_double_complex z1 = { .f = z }; diff --git a/src/s_conj.c b/src/s_conj.c index 4c1ae3c..b7a63d1 100644 --- a/src/s_conj.c +++ b/src/s_conj.c @@ -30,7 +30,7 @@ #include "math_private.h" -double complex +DLLEXPORT double complex conj(double complex z) { diff --git a/src/s_conjf.c b/src/s_conjf.c index 96b489a..02793c4 100644 --- a/src/s_conjf.c +++ b/src/s_conjf.c @@ -30,7 +30,7 @@ #include "math_private.h" -float complex +DLLEXPORT float complex conjf(float complex z) { diff --git a/src/s_conjl.c b/src/s_conjl.c index 7047ff4..774b8f6 100644 --- a/src/s_conjl.c +++ b/src/s_conjl.c @@ -30,7 +30,7 @@ #include "math_private.h" -long double complex +DLLEXPORT long double complex conjl(long double complex z) { diff --git a/src/s_copysign.c b/src/s_copysign.c index 7708ace..12b777d 100644 --- a/src/s_copysign.c +++ b/src/s_copysign.c @@ -22,7 +22,7 @@ #include "openlibm.h" #include "math_private.h" -double +DLLEXPORT double copysign(double x, double y) { u_int32_t hx,hy; diff --git a/src/s_copysignf.c b/src/s_copysignf.c index 34ace7e..f7e8415 100644 --- a/src/s_copysignf.c +++ b/src/s_copysignf.c @@ -25,7 +25,7 @@ #include "openlibm.h" #include "math_private.h" -float +DLLEXPORT float copysignf(float x, float y) { u_int32_t ix,iy; diff --git a/src/s_copysignl.c b/src/s_copysignl.c index b6f6f28..b648443 100644 --- a/src/s_copysignl.c +++ b/src/s_copysignl.c @@ -27,10 +27,10 @@ */ #include - +#include "math_private.h" #include "fpmath.h" -long double +DLLEXPORT long double copysignl(long double x, long double y) { union IEEEl2bits ux, uy; diff --git a/src/s_cos.c b/src/s_cos.c index 20a2fe7..2d3c4bb 100644 --- a/src/s_cos.c +++ b/src/s_cos.c @@ -51,7 +51,7 @@ #include "math_private.h" //#include "e_rem_pio2.c" -double +DLLEXPORT double cos(double x) { double y[2],z=0.0; diff --git a/src/s_cosf.c b/src/s_cosf.c index 1faf534..1d4e725 100644 --- a/src/s_cosf.c +++ b/src/s_cosf.c @@ -35,7 +35,7 @@ c2pio2 = 2*M_PI_2, /* 0x400921FB, 0x54442D18 */ c3pio2 = 3*M_PI_2, /* 0x4012D97C, 0x7F3321D2 */ c4pio2 = 4*M_PI_2; /* 0x401921FB, 0x54442D18 */ -float +DLLEXPORT float cosf(float x) { double y; diff --git a/src/s_cosl.c b/src/s_cosl.c index 00b9c87..d9147ae 100644 --- a/src/s_cosl.c +++ b/src/s_cosl.c @@ -44,7 +44,7 @@ #error "Unsupported long double format" #endif -long double +DLLEXPORT long double cosl(long double x) { union IEEEl2bits z; diff --git a/src/s_cpow.c b/src/s_cpow.c index 044fe90..9c618b4 100644 --- a/src/s_cpow.c +++ b/src/s_cpow.c @@ -47,8 +47,9 @@ #include #include #include +#include "math_private.h" -double complex +DLLEXPORT double complex cpow(double complex a, double complex z) { double complex w; diff --git a/src/s_cpowf.c b/src/s_cpowf.c index 994fa06..a01398b 100644 --- a/src/s_cpowf.c +++ b/src/s_cpowf.c @@ -46,8 +46,9 @@ #include #include +#include "math_private.h" -float complex +DLLEXPORT float complex cpowf(float complex a, float complex z) { float complex w; diff --git a/src/s_cpowl.c b/src/s_cpowl.c index 69d14d2..3103763 100644 --- a/src/s_cpowl.c +++ b/src/s_cpowl.c @@ -47,8 +47,9 @@ #include #include +#include "math_private.h" -long double complex +DLLEXPORT long double complex cpowl(long double complex a, long double complex z) { long double complex w; diff --git a/src/s_cproj.c b/src/s_cproj.c index 16679d8..67f0bd0 100644 --- a/src/s_cproj.c +++ b/src/s_cproj.c @@ -32,7 +32,7 @@ #include "math_private.h" -double complex +DLLEXPORT double complex cproj(double complex z) { diff --git a/src/s_cprojf.c b/src/s_cprojf.c index 410ebe8..4750b6b 100644 --- a/src/s_cprojf.c +++ b/src/s_cprojf.c @@ -32,7 +32,7 @@ #include "math_private.h" -float complex +DLLEXPORT float complex cprojf(float complex z) { diff --git a/src/s_cprojl.c b/src/s_cprojl.c index 3edfed1..a048107 100644 --- a/src/s_cprojl.c +++ b/src/s_cprojl.c @@ -32,7 +32,7 @@ #include "math_private.h" -long double complex +DLLEXPORT long double complex cprojl(long double complex z) { diff --git a/src/s_creal.c b/src/s_creal.c index ad14cdf..faf5999 100644 --- a/src/s_creal.c +++ b/src/s_creal.c @@ -27,8 +27,9 @@ */ #include +#include "math_private.h" -double +DLLEXPORT double creal(double complex z) { return z; diff --git a/src/s_crealf.c b/src/s_crealf.c index a5c1562..f17d595 100644 --- a/src/s_crealf.c +++ b/src/s_crealf.c @@ -27,8 +27,9 @@ */ #include +#include "math_private.h" -float +DLLEXPORT float crealf(float complex z) { return z; diff --git a/src/s_creall.c b/src/s_creall.c index 1531d16..2d1ded7 100644 --- a/src/s_creall.c +++ b/src/s_creall.c @@ -27,8 +27,9 @@ */ #include +#include "math_private.h" -long double +DLLEXPORT long double creall(long double complex z) { return z; diff --git a/src/s_csinh.c b/src/s_csinh.c index 69a6415..8402c49 100644 --- a/src/s_csinh.c +++ b/src/s_csinh.c @@ -44,7 +44,7 @@ static const double huge = 0x1p1023; -double complex +DLLEXPORT double complex csinh(double complex z) { double x, y, h; @@ -147,7 +147,7 @@ csinh(double complex z) return (cpack((x * x) * (y - y), (x + x) * (y - y))); } -double complex +DLLEXPORT double complex csin(double complex z) { diff --git a/src/s_csinhf.c b/src/s_csinhf.c index 583ff18..78ed180 100644 --- a/src/s_csinhf.c +++ b/src/s_csinhf.c @@ -38,7 +38,7 @@ static const float huge = 0x1p127; -float complex +DLLEXPORT float complex csinhf(float complex z) { float x, y, h; @@ -96,7 +96,7 @@ csinhf(float complex z) return (cpackf((x * x) * (y - y), (x + x) * (y - y))); } -float complex +DLLEXPORT float complex csinf(float complex z) { diff --git a/src/s_csqrt.c b/src/s_csqrt.c index 8e20850..efd0701 100644 --- a/src/s_csqrt.c +++ b/src/s_csqrt.c @@ -47,7 +47,7 @@ /* We risk spurious overflow for components >= DBL_MAX / (1 + sqrt(2)). */ #define THRESH 0x1.a827999fcef32p+1022 -double complex +DLLEXPORT double complex csqrt(double complex z) { double complex result; diff --git a/src/s_csqrtf.c b/src/s_csqrtf.c index f94c997..7a79c88 100644 --- a/src/s_csqrtf.c +++ b/src/s_csqrtf.c @@ -43,7 +43,7 @@ #pragma STDC CX_LIMITED_RANGE ON #endif -float complex +DLLEXPORT float complex csqrtf(float complex z) { float a = crealf(z), b = cimagf(z); diff --git a/src/s_csqrtl.c b/src/s_csqrtl.c index 8983c2a..c542a13 100644 --- a/src/s_csqrtl.c +++ b/src/s_csqrtl.c @@ -46,7 +46,7 @@ /* We risk spurious overflow for components >= LDBL_MAX / (1 + sqrt(2)). */ #define THRESH (LDBL_MAX / 2.414213562373095048801688724209698L) -long double complex +DLLEXPORT long double complex csqrtl(long double complex z) { long double complex result; diff --git a/src/s_ctanh.c b/src/s_ctanh.c index 073d661..d26d128 100644 --- a/src/s_ctanh.c +++ b/src/s_ctanh.c @@ -71,7 +71,7 @@ #include "math_private.h" -double complex +DLLEXPORT double complex ctanh(double complex z) { double x, y; @@ -134,7 +134,7 @@ ctanh(double complex z) return (cpack((beta * rho * s) / denom, t / denom)); } -double complex +DLLEXPORT double complex ctan(double complex z) { diff --git a/src/s_ctanhf.c b/src/s_ctanhf.c index a2a6c49..c123b54 100644 --- a/src/s_ctanhf.c +++ b/src/s_ctanhf.c @@ -36,7 +36,7 @@ #include "math_private.h" -float complex +DLLEXPORT float complex ctanhf(float complex z) { float x, y; @@ -74,7 +74,7 @@ ctanhf(float complex z) return (cpackf((beta * rho * s) / denom, t / denom)); } -float complex +DLLEXPORT float complex ctanf(float complex z) { diff --git a/src/s_erf.c b/src/s_erf.c index 6f4fa23..d6504b5 100644 --- a/src/s_erf.c +++ b/src/s_erf.c @@ -186,7 +186,7 @@ sb5 = 2.55305040643316442583e+03, /* 0x40A3F219, 0xCEDF3BE6 */ sb6 = 4.74528541206955367215e+02, /* 0x407DA874, 0xE79FE763 */ sb7 = -2.24409524465858183362e+01; /* 0xC03670E2, 0x42712D62 */ -double +DLLEXPORT double erf(double x) { int32_t hx,ix,i; @@ -238,7 +238,7 @@ erf(double x) if(hx>=0) return one-r/x; else return r/x-one; } -double +DLLEXPORT double erfc(double x) { int32_t hx,ix; diff --git a/src/s_erff.c b/src/s_erff.c index 3351c47..29a8d80 100644 --- a/src/s_erff.c +++ b/src/s_erff.c @@ -94,7 +94,7 @@ sb5 = 2.5530502930e+03, /* 0x451f90ce */ sb6 = 4.7452853394e+02, /* 0x43ed43a7 */ sb7 = -2.2440952301e+01; /* 0xc1b38712 */ -float +DLLEXPORT float erff(float x) { int32_t hx,ix,i; @@ -147,7 +147,7 @@ erff(float x) if(hx>=0) return one-r/x; else return r/x-one; } -float +DLLEXPORT float erfcf(float x) { int32_t hx,ix; diff --git a/src/s_exp2.c b/src/s_exp2.c index f76993c..598268b 100644 --- a/src/s_exp2.c +++ b/src/s_exp2.c @@ -337,7 +337,7 @@ static const double tbl[TBLSIZE * 2] = { * Gal, S. and Bachelis, B. An Accurate Elementary Mathematical Library * for the IEEE Floating Point Standard. TOMS 17(1), 26-46 (1991). */ -double +DLLEXPORT double exp2(double x) { double r, t, twopk, twopkp1000, z; diff --git a/src/s_exp2f.c b/src/s_exp2f.c index 4623889..9409bd4 100644 --- a/src/s_exp2f.c +++ b/src/s_exp2f.c @@ -90,7 +90,7 @@ static const double exp2ft[TBLSIZE] = { * Tang, P. Table-driven Implementation of the Exponential Function * in IEEE Floating-Point Arithmetic. TOMS 15(2), 144-157 (1989). */ -float +DLLEXPORT float exp2f(float x) { double tv, twopk, u, z; diff --git a/src/s_expm1.c b/src/s_expm1.c index 8f002c8..ff99652 100644 --- a/src/s_expm1.c +++ b/src/s_expm1.c @@ -128,7 +128,7 @@ Q3 = -7.93650757867487942473e-05, /* BF14CE19 9EAADBB7 */ Q4 = 4.00821782732936239552e-06, /* 3ED0CFCA 86E65239 */ Q5 = -2.01099218183624371326e-07; /* BE8AFDB7 6E09C32D */ -double +DLLEXPORT double expm1(double x) { double y,hi,lo,c,t,e,hxs,hfx,r1,twopk; diff --git a/src/s_expm1f.c b/src/s_expm1f.c index 1ea85a8..9d29ac4 100644 --- a/src/s_expm1f.c +++ b/src/s_expm1f.c @@ -37,7 +37,7 @@ invln2 = 1.4426950216e+00,/* 0x3fb8aa3b */ Q1 = -3.3333212137e-2, /* -0x888868.0p-28 */ Q2 = 1.5807170421e-3; /* 0xcf3010.0p-33 */ -float +DLLEXPORT float expm1f(float x) { float y,hi,lo,c,t,e,hxs,hfx,r1,twopk; diff --git a/src/s_fabs.c b/src/s_fabs.c index 064a6e2..6c53c80 100644 --- a/src/s_fabs.c +++ b/src/s_fabs.c @@ -17,7 +17,7 @@ #include "openlibm.h" #include "math_private.h" -double +DLLEXPORT double fabs(double x) { u_int32_t high; diff --git a/src/s_fabsf.c b/src/s_fabsf.c index 6b9f4bc..f5b4525 100644 --- a/src/s_fabsf.c +++ b/src/s_fabsf.c @@ -23,7 +23,7 @@ #include "openlibm.h" #include "math_private.h" -float +DLLEXPORT float fabsf(float x) { u_int32_t ix; diff --git a/src/s_fabsl.c b/src/s_fabsl.c index c4d1b44..68e922b 100644 --- a/src/s_fabsl.c +++ b/src/s_fabsl.c @@ -29,10 +29,10 @@ */ #include - +#include "math_private.h" #include "fpmath.h" -long double +DLLEXPORT long double fabsl(long double x) { union IEEEl2bits u; diff --git a/src/s_fdim.c b/src/s_fdim.c index 43f4985..12dee6b 100644 --- a/src/s_fdim.c +++ b/src/s_fdim.c @@ -26,11 +26,11 @@ #include "cdefs-compat.h" //__FBSDID("$FreeBSD: src/lib/msun/src/s_fdim.c,v 1.1 2004/06/30 07:04:01 das Exp $"); - #include +#include "math_private.h" #define DECL(type, fn) \ -type \ +DLLEXPORT type \ fn(type x, type y) \ { \ \ diff --git a/src/s_finite.c b/src/s_finite.c index c5962f7..46f84e7 100644 --- a/src/s_finite.c +++ b/src/s_finite.c @@ -21,7 +21,8 @@ #include "openlibm.h" #include "math_private.h" - int finite(double x) +DLLEXPORT int +finite(double x) { int32_t hx; GET_HIGH_WORD(hx,x); diff --git a/src/s_finitef.c b/src/s_finitef.c index 779f24b..6d189a5 100644 --- a/src/s_finitef.c +++ b/src/s_finitef.c @@ -24,7 +24,8 @@ #include "openlibm.h" #include "math_private.h" - int finitef(float x) +DLLEXPORT int +finitef(float x) { int32_t ix; GET_FLOAT_WORD(ix,x); diff --git a/src/s_floor.c b/src/s_floor.c index 39e08c9..95891ef 100644 --- a/src/s_floor.c +++ b/src/s_floor.c @@ -29,7 +29,7 @@ static const double huge = 1.0e300; -double +DLLEXPORT double floor(double x) { int32_t i0,i1,j0; diff --git a/src/s_floorf.c b/src/s_floorf.c index 7676bd2..678bb11 100644 --- a/src/s_floorf.c +++ b/src/s_floorf.c @@ -30,7 +30,7 @@ static const float huge = 1.0e30; -float +DLLEXPORT float floorf(float x) { int32_t i0,j0; diff --git a/src/s_floorl.c b/src/s_floorl.c index 2d4ad46..df578cc 100644 --- a/src/s_floorl.c +++ b/src/s_floorl.c @@ -25,6 +25,7 @@ #include #include +#include "math_private.h" #include #include "fpmath.h" @@ -51,7 +52,7 @@ static const long double huge = 1.0e300; -long double +DLLEXPORT long double floorl(long double x) { union IEEEl2bits u = { .e = x }; diff --git a/src/s_fma.c b/src/s_fma.c index 757eea4..1c93f55 100644 --- a/src/s_fma.c +++ b/src/s_fma.c @@ -174,7 +174,7 @@ dd_mul(double a, double b) * Hardware instructions should be used on architectures that support it, * since this implementation will likely be several times slower. */ -double +DLLEXPORT double fma(double x, double y, double z) { double xs, ys, zs, adj; diff --git a/src/s_fmaf.c b/src/s_fmaf.c index 37b1618..da0e019 100644 --- a/src/s_fmaf.c +++ b/src/s_fmaf.c @@ -39,7 +39,7 @@ * direct double-precision arithmetic suffices, except where double * rounding occurs. */ -float +DLLEXPORT float fmaf(float x, float y, float z) { double xy, result; diff --git a/src/s_fmal.c b/src/s_fmal.c index c484448..1d6b6b8 100644 --- a/src/s_fmal.c +++ b/src/s_fmal.c @@ -30,7 +30,7 @@ #include #include #include - +#include "math_private.h" #include "fpmath.h" /* @@ -162,7 +162,7 @@ dd_mul(long double a, long double b) * Dekker, T. A Floating-Point Technique for Extending the * Available Precision. Numer. Math. 18, 224-242 (1971). */ -long double +DLLEXPORT long double fmal(long double x, long double y, long double z) { long double xs, ys, zs, adj; diff --git a/src/s_fmax.c b/src/s_fmax.c index a1c2132..5df93d8 100644 --- a/src/s_fmax.c +++ b/src/s_fmax.c @@ -28,10 +28,10 @@ //__FBSDID("$FreeBSD: src/lib/msun/src/s_fmax.c,v 1.1 2004/06/30 07:04:01 das Exp $"); #include - +#include "math_private.h" #include "fpmath.h" -double +DLLEXPORT double fmax(double x, double y) { union IEEEd2bits u[2]; diff --git a/src/s_fmaxf.c b/src/s_fmaxf.c index abe6d27..e1744c0 100644 --- a/src/s_fmaxf.c +++ b/src/s_fmaxf.c @@ -28,10 +28,10 @@ //__FBSDID("$FreeBSD: src/lib/msun/src/s_fmaxf.c,v 1.1 2004/06/30 07:04:01 das Exp $"); #include - +#include "math_private.h" #include "fpmath.h" -float +DLLEXPORT float fmaxf(float x, float y) { union IEEEf2bits u[2]; diff --git a/src/s_fmaxl.c b/src/s_fmaxl.c index 76d4809..b77fbca 100644 --- a/src/s_fmaxl.c +++ b/src/s_fmaxl.c @@ -28,10 +28,10 @@ //__FBSDID("$FreeBSD: src/lib/msun/src/s_fmaxl.c,v 1.1 2004/06/30 07:04:01 das Exp $"); #include - +#include "math_private.h" #include "fpmath.h" -long double +DLLEXPORT long double fmaxl(long double x, long double y) { union IEEEl2bits u[2]; diff --git a/src/s_fmin.c b/src/s_fmin.c index bea5722..20474e6 100644 --- a/src/s_fmin.c +++ b/src/s_fmin.c @@ -28,10 +28,10 @@ //__FBSDID("$FreeBSD: src/lib/msun/src/s_fmin.c,v 1.1 2004/06/30 07:04:01 das Exp $"); #include - +#include "math_private.h" #include "fpmath.h" -double +DLLEXPORT double fmin(double x, double y) { union IEEEd2bits u[2]; diff --git a/src/s_fminf.c b/src/s_fminf.c index bcbb97c..f537186 100644 --- a/src/s_fminf.c +++ b/src/s_fminf.c @@ -28,10 +28,10 @@ //__FBSDID("$FreeBSD: src/lib/msun/src/s_fminf.c,v 1.1 2004/06/30 07:04:01 das Exp $"); #include - +#include "math_private.h" #include "fpmath.h" -float +DLLEXPORT float fminf(float x, float y) { union IEEEf2bits u[2]; diff --git a/src/s_fminl.c b/src/s_fminl.c index 94f8b52..db30e18 100644 --- a/src/s_fminl.c +++ b/src/s_fminl.c @@ -28,10 +28,10 @@ //__FBSDID("$FreeBSD: src/lib/msun/src/s_fminl.c,v 1.1 2004/06/30 07:04:01 das Exp $"); #include - +#include "math_private.h" #include "fpmath.h" -long double +DLLEXPORT long double fminl(long double x, long double y) { union IEEEl2bits u[2]; diff --git a/src/s_fpclassify.c b/src/s_fpclassify.c index 1f77597..79efb62 100644 --- a/src/s_fpclassify.c +++ b/src/s_fpclassify.c @@ -26,10 +26,10 @@ */ #include - +#include "math_private.h" #include "fpmath.h" -int +DLLEXPORT int __fpclassifyd(double d) { union IEEEd2bits u; @@ -51,7 +51,7 @@ __fpclassifyd(double d) } -int +DLLEXPORT int __fpclassifyf(float f) { union IEEEf2bits u; @@ -72,7 +72,7 @@ __fpclassifyf(float f) } } -int +DLLEXPORT int __fpclassifyl(long double e) { union IEEEl2bits u; diff --git a/src/s_frexp.c b/src/s_frexp.c index ae8d647..3a7d3bd 100644 --- a/src/s_frexp.c +++ b/src/s_frexp.c @@ -31,7 +31,7 @@ static const double two54 = 1.80143985094819840000e+16; /* 0x43500000, 0x00000000 */ -double +DLLEXPORT double frexp(double x, int *eptr) { int32_t hx, ix, lx; diff --git a/src/s_frexpf.c b/src/s_frexpf.c index b1f51cb..1f31aa1 100644 --- a/src/s_frexpf.c +++ b/src/s_frexpf.c @@ -22,7 +22,7 @@ static const float two25 = 3.3554432000e+07; /* 0x4c000000 */ -float +DLLEXPORT float frexpf(float x, int *eptr) { int32_t hx,ix; diff --git a/src/s_frexpl.c b/src/s_frexpl.c index b93ecce..c8c61f7 100644 --- a/src/s_frexpl.c +++ b/src/s_frexpl.c @@ -28,14 +28,14 @@ #include #include - +#include "math_private.h" #include "fpmath.h" #if LDBL_MAX_EXP != 0x4000 #error "Unsupported long double format" #endif -long double +DLLEXPORT long double frexpl(long double x, int *ex) { union IEEEl2bits u; diff --git a/src/s_ilogb.c b/src/s_ilogb.c index d8bc4bf..14e4265 100644 --- a/src/s_ilogb.c +++ b/src/s_ilogb.c @@ -25,7 +25,8 @@ #include "openlibm.h" #include "math_private.h" - int ilogb(double x) +DLLEXPORT int +ilogb(double x) { int32_t hx,lx,ix; diff --git a/src/s_ilogbf.c b/src/s_ilogbf.c index f7501cd..d83857c 100644 --- a/src/s_ilogbf.c +++ b/src/s_ilogbf.c @@ -21,7 +21,8 @@ #include "openlibm.h" #include "math_private.h" - int ilogbf(float x) +DLLEXPORT int +ilogbf(float x) { int32_t hx,ix; diff --git a/src/s_ilogbl.c b/src/s_ilogbl.c index f75c971..6bbe124 100644 --- a/src/s_ilogbl.c +++ b/src/s_ilogbl.c @@ -16,10 +16,10 @@ #include #include #include - +#include "math_private.h" #include "fpmath.h" -int +DLLEXPORT int ilogbl(long double x) { union IEEEl2bits u; diff --git a/src/s_isfinite.c b/src/s_isfinite.c index cb4ea17..0ae3549 100644 --- a/src/s_isfinite.c +++ b/src/s_isfinite.c @@ -27,10 +27,10 @@ */ #include - +#include "math_private.h" #include "fpmath.h" -int +DLLEXPORT int __isfinite(double d) { union IEEEd2bits u; @@ -39,7 +39,7 @@ __isfinite(double d) return (u.bits.exp != 2047); } -int +DLLEXPORT int __isfinitef(float f) { union IEEEf2bits u; @@ -48,7 +48,7 @@ __isfinitef(float f) return (u.bits.exp != 255); } -int +DLLEXPORT int __isfinitel(long double e) { union IEEEl2bits u; diff --git a/src/s_isinf.c b/src/s_isinf.c index d59d995..08c3cce 100644 --- a/src/s_isinf.c +++ b/src/s_isinf.c @@ -26,12 +26,12 @@ */ #include - +#include "math_private.h" #include "fpmath.h" /* Provided by libc */ #if 1 -int +DLLEXPORT int (isinf) (double d) { union IEEEd2bits u; @@ -41,7 +41,7 @@ int } #endif -int +DLLEXPORT int __isinff(float f) { union IEEEf2bits u; @@ -50,7 +50,7 @@ __isinff(float f) return (u.bits.exp == 255 && u.bits.man == 0); } -int +DLLEXPORT int __isinfl(long double e) { union IEEEl2bits u; diff --git a/src/s_isnan.c b/src/s_isnan.c index 7477e81..5d903e9 100644 --- a/src/s_isnan.c +++ b/src/s_isnan.c @@ -27,12 +27,12 @@ */ #include - +#include "math_private.h" #include "fpmath.h" /* Provided by libc */ #if 1 -int +DLLEXPORT int (isnan) (double d) { union IEEEd2bits u; @@ -42,7 +42,7 @@ int } #endif -int +DLLEXPORT int __isnanf(float f) { union IEEEf2bits u; @@ -51,7 +51,7 @@ __isnanf(float f) return (u.bits.exp == 255 && u.bits.man != 0); } -int +DLLEXPORT int __isnanl(long double e) { union IEEEl2bits u; diff --git a/src/s_isnormal.c b/src/s_isnormal.c index 2b624c3..bcdb86b 100644 --- a/src/s_isnormal.c +++ b/src/s_isnormal.c @@ -27,10 +27,10 @@ */ #include - +#include "math_private.h" #include "fpmath.h" -int +DLLEXPORT int __isnormal(double d) { union IEEEd2bits u; @@ -39,7 +39,7 @@ __isnormal(double d) return (u.bits.exp != 0 && u.bits.exp != 2047); } -int +DLLEXPORT int __isnormalf(float f) { union IEEEf2bits u; @@ -48,7 +48,7 @@ __isnormalf(float f) return (u.bits.exp != 0 && u.bits.exp != 255); } -int +DLLEXPORT int __isnormall(long double e) { union IEEEl2bits u; diff --git a/src/s_log1p.c b/src/s_log1p.c index 5472650..2752cd2 100644 --- a/src/s_log1p.c +++ b/src/s_log1p.c @@ -97,7 +97,7 @@ Lp7 = 1.479819860511658591e-01; /* 3FC2F112 DF3E5244 */ static const double zero = 0.0; -double +DLLEXPORT double log1p(double x) { double hfsq,f,c,s,z,R,u; diff --git a/src/s_log1pf.c b/src/s_log1pf.c index 5f91b84..04120ea 100644 --- a/src/s_log1pf.c +++ b/src/s_log1pf.c @@ -35,7 +35,7 @@ Lp7 = 1.4798198640e-01; /* 3E178897 */ static const float zero = 0.0; -float +DLLEXPORT float log1pf(float x) { float hfsq,f,c,s,z,R,u; diff --git a/src/s_logb.c b/src/s_logb.c index 0407271..4fb376a 100644 --- a/src/s_logb.c +++ b/src/s_logb.c @@ -27,7 +27,7 @@ static const double two54 = 1.80143985094819840000e+16; /* 43500000 00000000 */ -double +DLLEXPORT double logb(double x) { int32_t lx,ix; diff --git a/src/s_logbf.c b/src/s_logbf.c index 586a917..26a36b5 100644 --- a/src/s_logbf.c +++ b/src/s_logbf.c @@ -22,7 +22,7 @@ static const float two25 = 3.355443200e+07; /* 0x4c000000 */ -float +DLLEXPORT float logbf(float x) { int32_t ix; diff --git a/src/s_logbl.c b/src/s_logbl.c index a2c2e89..44a6eab 100644 --- a/src/s_logbl.c +++ b/src/s_logbl.c @@ -13,10 +13,10 @@ #include #include #include - +#include "math_private.h" #include "fpmath.h" -long double +DLLEXPORT long double logbl(long double x) { union IEEEl2bits u; diff --git a/src/s_lrint.c b/src/s_lrint.c index 7630fca..ff7bce7 100644 --- a/src/s_lrint.c +++ b/src/s_lrint.c @@ -27,6 +27,7 @@ #include "cdefs-compat.h" #include #include +#include "math_private.h" #ifndef type //__FBSDID("$FreeBSD: src/lib/msun/src/s_lrint.c,v 1.1 2005/01/11 23:12:55 das Exp $"); @@ -43,7 +44,7 @@ * significant bits than 'type'. Hence, we bend over backwards for the * sake of correctness; an MD implementation could be more efficient. */ -dtype +DLLEXPORT dtype fn(type x) { fenv_t env; diff --git a/src/s_lround.c b/src/s_lround.c index f852ccc..f89d031 100644 --- a/src/s_lround.c +++ b/src/s_lround.c @@ -31,6 +31,7 @@ //VBS end #include #include +#include "math_private.h" #ifndef type //__FBSDID("$FreeBSD: src/lib/msun/src/s_lround.c,v 1.2 2005/04/08 00:52:16 das Exp $"); @@ -55,7 +56,7 @@ static const type dtype_max = DTYPE_MAX + 0.5; #define INRANGE(x) (dtype_max - DTYPE_MAX != 0.5 || \ ((x) > dtype_min && (x) < dtype_max)) -dtype +DLLEXPORT dtype fn(type x) { diff --git a/src/s_modf.c b/src/s_modf.c index 52ca7ba..2615729 100644 --- a/src/s_modf.c +++ b/src/s_modf.c @@ -25,7 +25,7 @@ static const double one = 1.0; -double +DLLEXPORT double modf(double x, double *iptr) { int32_t i0,i1,j0; diff --git a/src/s_modff.c b/src/s_modff.c index 9de60ff..7deee06 100644 --- a/src/s_modff.c +++ b/src/s_modff.c @@ -21,7 +21,7 @@ static const float one = 1.0; -float +DLLEXPORT float modff(float x, float *iptr) { int32_t i0,j0; diff --git a/src/s_modfl.c b/src/s_modfl.c index a5bd53d..d6d9ab6 100644 --- a/src/s_modfl.c +++ b/src/s_modfl.c @@ -38,6 +38,7 @@ #include #include +#include "math_private.h" #include #include "fpmath.h" @@ -54,7 +55,7 @@ static const long double zero[] = { 0.0L, -0.0L }; -long double +DLLEXPORT long double modfl(long double x, long double *iptr) { union IEEEl2bits ux; diff --git a/src/s_nan.c b/src/s_nan.c index c16bade..ef8aad7 100644 --- a/src/s_nan.c +++ b/src/s_nan.c @@ -61,7 +61,7 @@ static __inline int digittoint(int c) { * consider valid, so we might be violating the C standard. But it's * impossible to use nan(3) portably anyway, so this seems good enough. */ -void +DLLEXPORT void _scan_nan(u_int32_t *words, int num_words, const char *s) { int si; /* index into s */ @@ -89,7 +89,7 @@ _scan_nan(u_int32_t *words, int num_words, const char *s) } } -double +DLLEXPORT double nan(const char *s) { union { @@ -106,7 +106,7 @@ nan(const char *s) return (u.d); } -float +DLLEXPORT float nanf(const char *s) { union { diff --git a/src/s_nearbyint.c b/src/s_nearbyint.c index ea1271d..6de1d77 100644 --- a/src/s_nearbyint.c +++ b/src/s_nearbyint.c @@ -29,6 +29,7 @@ #include #include +#include "math_private.h" /* * We save and restore the floating-point environment to avoid raising @@ -38,7 +39,7 @@ * rounding can't overflow as long as emax >= p. */ #define DECL(type, fn, rint) \ -type \ +DLLEXPORT type \ fn(type x) \ { \ type ret; \ diff --git a/src/s_nextafter.c b/src/s_nextafter.c index cabdd32..f3c695e 100644 --- a/src/s_nextafter.c +++ b/src/s_nextafter.c @@ -25,7 +25,7 @@ #include "openlibm.h" #include "math_private.h" -double +DLLEXPORT double nextafter(double x, double y) { volatile double t; diff --git a/src/s_nextafterf.c b/src/s_nextafterf.c index 233a15b..1a070a6 100644 --- a/src/s_nextafterf.c +++ b/src/s_nextafterf.c @@ -19,7 +19,7 @@ #include "openlibm.h" #include "math_private.h" -float +DLLEXPORT float nextafterf(float x, float y) { volatile float t; diff --git a/src/s_nextafterl.c b/src/s_nextafterl.c index 72c5874..2874350 100644 --- a/src/s_nextafterl.c +++ b/src/s_nextafterl.c @@ -30,7 +30,7 @@ #error "Unsupported long double format" #endif -long double +DLLEXPORT long double nextafterl(long double x, long double y) { volatile long double t; diff --git a/src/s_nexttoward.c b/src/s_nexttoward.c index 6112555..eca0370 100644 --- a/src/s_nexttoward.c +++ b/src/s_nexttoward.c @@ -29,7 +29,7 @@ #error "Unsupported long double format" #endif -double +DLLEXPORT double nexttoward(double x, long double y) { union IEEEl2bits uy; diff --git a/src/s_nexttowardf.c b/src/s_nexttowardf.c index bb93378..4c1902f 100644 --- a/src/s_nexttowardf.c +++ b/src/s_nexttowardf.c @@ -20,7 +20,7 @@ #define LDBL_INFNAN_EXP (LDBL_MAX_EXP * 2 - 1) -float +DLLEXPORT float nexttowardf(float x, long double y) { union IEEEl2bits uy; diff --git a/src/s_remquo.c b/src/s_remquo.c index 85f19ab..7b63b76 100644 --- a/src/s_remquo.c +++ b/src/s_remquo.c @@ -28,7 +28,7 @@ static const double Zero[] = {0.0, -0.0,}; * method. In practice, this is far more bits than are needed to use * remquo in reduction algorithms. */ -double +DLLEXPORT double remquo(double x, double y, int *quo) { int32_t n,hx,hy,hz,ix,iy,sx,i; diff --git a/src/s_remquof.c b/src/s_remquof.c index 9d0f79a..579dec0 100644 --- a/src/s_remquof.c +++ b/src/s_remquof.c @@ -26,7 +26,7 @@ static const float Zero[] = {0.0, -0.0,}; * method. In practice, this is far more bits than are needed to use * remquo in reduction algorithms. */ -float +DLLEXPORT float remquof(float x, float y, int *quo) { int32_t n,hx,hy,hz,ix,iy,sx,i; diff --git a/src/s_remquol.c b/src/s_remquol.c index 3909d66..0100cae 100644 --- a/src/s_remquol.c +++ b/src/s_remquol.c @@ -64,7 +64,7 @@ static const long double Zero[] = {0.0L, -0.0L}; * - The high part of the mantissa fits in an int64_t with enough room * for an explicit integer bit in front of the fractional bits. */ -long double +DLLEXPORT long double remquol(long double x, long double y, int *quo) { union IEEEl2bits ux, uy; diff --git a/src/s_rint.c b/src/s_rint.c index 68abf3f..b246e00 100644 --- a/src/s_rint.c +++ b/src/s_rint.c @@ -34,7 +34,7 @@ TWO52[2]={ -4.50359962737049600000e+15, /* 0xC3300000, 0x00000000 */ }; -double +DLLEXPORT double rint(double x) { int32_t i0,j0,sx; diff --git a/src/s_rintf.c b/src/s_rintf.c index d59a92e..91d05c0 100644 --- a/src/s_rintf.c +++ b/src/s_rintf.c @@ -28,7 +28,7 @@ TWO23[2]={ -8.3886080000e+06, /* 0xcb000000 */ }; -float +DLLEXPORT float rintf(float x) { int32_t i0,j0,sx; diff --git a/src/s_rintl.c b/src/s_rintl.c index 5cef540..3dd11c2 100644 --- a/src/s_rintl.c +++ b/src/s_rintl.c @@ -55,7 +55,7 @@ shift[2] = { }; static const float zero[2] = { 0.0, -0.0 }; -long double +DLLEXPORT long double rintl(long double x) { union IEEEl2bits u; diff --git a/src/s_round.c b/src/s_round.c index b86e01c..ec758fc 100644 --- a/src/s_round.c +++ b/src/s_round.c @@ -28,8 +28,9 @@ //__FBSDID("$FreeBSD: src/lib/msun/src/s_round.c,v 1.4 2005/12/02 13:45:06 bde Exp $"); #include +#include "math_private.h" -double +DLLEXPORT double round(double x) { double t; diff --git a/src/s_roundf.c b/src/s_roundf.c index 010ac2b..cbb40f3 100644 --- a/src/s_roundf.c +++ b/src/s_roundf.c @@ -28,8 +28,9 @@ //__FBSDID("$FreeBSD: src/lib/msun/src/s_roundf.c,v 1.4 2005/12/02 13:45:06 bde Exp $"); #include +#include "math_private.h" -float +DLLEXPORT float roundf(float x) { float t; diff --git a/src/s_roundl.c b/src/s_roundl.c index fba814a..24e58aa 100644 --- a/src/s_roundl.c +++ b/src/s_roundl.c @@ -28,8 +28,9 @@ //__FBSDID("$FreeBSD: src/lib/msun/src/s_roundl.c,v 1.2 2005/12/02 13:45:06 bde Exp $"); #include +#include "math_private.h" -long double +DLLEXPORT long double roundl(long double x) { long double t; diff --git a/src/s_scalbln.c b/src/s_scalbln.c index 345f56e..ba0f699 100644 --- a/src/s_scalbln.c +++ b/src/s_scalbln.c @@ -29,8 +29,9 @@ #include #include +#include "math_private.h" -double +DLLEXPORT double scalbln (double x, long n) { int in; @@ -45,7 +46,7 @@ scalbln (double x, long n) return (scalbn(x, in)); } -float +DLLEXPORT float scalblnf (float x, long n) { int in; @@ -60,7 +61,7 @@ scalblnf (float x, long n) return (scalbnf(x, in)); } -long double +DLLEXPORT long double scalblnl (long double x, long n) { int in; diff --git a/src/s_scalbn.c b/src/s_scalbn.c index 234eabe..de2c98f 100644 --- a/src/s_scalbn.c +++ b/src/s_scalbn.c @@ -29,7 +29,7 @@ twom54 = 5.55111512312578270212e-17, /* 0x3C900000, 0x00000000 */ huge = 1.0e+300, tiny = 1.0e-300; -double +DLLEXPORT double scalbn (double x, int n) { int32_t k,hx,lx; diff --git a/src/s_scalbnf.c b/src/s_scalbnf.c index cf16667..e353e24 100644 --- a/src/s_scalbnf.c +++ b/src/s_scalbnf.c @@ -25,7 +25,7 @@ twom25 = 2.9802322388e-08, /* 0x33000000 */ huge = 1.0e+30, tiny = 1.0e-30; -float +DLLEXPORT float scalbnf (float x, int n) { int32_t k,ix; diff --git a/src/s_scalbnl.c b/src/s_scalbnl.c index db1563a..2b7b4fe 100644 --- a/src/s_scalbnl.c +++ b/src/s_scalbnl.c @@ -26,7 +26,7 @@ #include "cdefs-compat.h" #include #include - +#include "math_private.h" #include "fpmath.h" #if LDBL_MAX_EXP != 0x4000 @@ -37,7 +37,7 @@ static const long double huge = 0x1p16000L, tiny = 0x1p-16000L; -long double +DLLEXPORT long double scalbnl (long double x, int n) { union IEEEl2bits u; diff --git a/src/s_signbit.c b/src/s_signbit.c index b07713a..bf3463b 100644 --- a/src/s_signbit.c +++ b/src/s_signbit.c @@ -27,10 +27,10 @@ */ #include - +#include "math_private.h" #include "fpmath.h" -int +DLLEXPORT int __signbit(double d) { union IEEEd2bits u; @@ -39,7 +39,7 @@ __signbit(double d) return (u.bits.sign); } -int +DLLEXPORT int __signbitf(float f) { union IEEEf2bits u; @@ -48,7 +48,7 @@ __signbitf(float f) return (u.bits.sign); } -int +DLLEXPORT int __signbitl(long double e) { union IEEEl2bits u; diff --git a/src/s_significand.c b/src/s_significand.c index 5971716..d6d06ea 100644 --- a/src/s_significand.c +++ b/src/s_significand.c @@ -22,7 +22,7 @@ #include "openlibm.h" #include "math_private.h" -double +DLLEXPORT double significand(double x) { return __ieee754_scalb(x,(double) -ilogb(x)); diff --git a/src/s_significandf.c b/src/s_significandf.c index 1c808b4..781f786 100644 --- a/src/s_significandf.c +++ b/src/s_significandf.c @@ -19,7 +19,7 @@ #include "openlibm.h" #include "math_private.h" -float +DLLEXPORT float significandf(float x) { return __ieee754_scalbf(x,(float) -ilogbf(x)); diff --git a/src/s_sin.c b/src/s_sin.c index 3beca2b..31ed70a 100644 --- a/src/s_sin.c +++ b/src/s_sin.c @@ -51,7 +51,7 @@ #include "math_private.h" //#include "e_rem_pio2.c" -double +DLLEXPORT double sin(double x) { double y[2],z=0.0; diff --git a/src/s_sincos.c b/src/s_sincos.c index f752a71..c10a15d 100644 --- a/src/s_sincos.c +++ b/src/s_sincos.c @@ -58,7 +58,7 @@ C4 = -2.75573143513906633035e-07, /* 0xBE927E4F, 0x809C52AD */ C5 = 2.08757232129817482790e-09, /* 0x3E21EE9E, 0xBDB4B1C4 */ C6 = -1.13596475577881948265e-11; /* 0xBDA8FAE9, 0xBE8838D4 */ -void +DLLEXPORT void __kernel_sincos( double x, double y, int iy, double * k_s, double * k_c ) { /* Inline calculation of sin/cos, as we can save @@ -86,7 +86,7 @@ __kernel_sincos( double x, double y, int iy, double * k_s, double * k_c ) *k_s = x-((z*(half*y-v*r)-y)-v*S1); } -void +DLLEXPORT void sincos(double x, double * s, double * c) { double y[2]; diff --git a/src/s_sincosf.c b/src/s_sincosf.c index aaf9242..26ebca0 100644 --- a/src/s_sincosf.c +++ b/src/s_sincosf.c @@ -40,7 +40,7 @@ C1 = 0x155553e1053a42.0p-57, /* 0.0416666233237390631894 */ C2 = -0x16c087e80f1e27.0p-62, /* -0.00138867637746099294692 */ C3 = 0x199342e0ee5069.0p-68; /* 0.0000243904487962774090654 */ -void +DLLEXPORT void __kernel_sincosdf( double x, float * s, float * c ) { double r, w, z, v; @@ -62,7 +62,7 @@ __kernel_sincosdf( double x, float * s, float * c ) *s = k_s; } -void +DLLEXPORT void sincosf(float x, float * s, float * c) { // Worst approximation of sin and cos NA *s = x; @@ -159,4 +159,4 @@ sincosf(float x, float * s, float * c) { } } -} \ No newline at end of file +} diff --git a/src/s_sincosl.c b/src/s_sincosl.c index c8295d6..2ca1e97 100644 --- a/src/s_sincosl.c +++ b/src/s_sincosl.c @@ -23,9 +23,9 @@ #error "Unsupported long double format" #endif - void - sincosl( long double x, long double * s, long double * c ) - { - *s = cosl( x ); - *c = sinl( x ); - } \ No newline at end of file +DLLEXPORT void +sincosl( long double x, long double * s, long double * c ) +{ + *s = cosl( x ); + *c = sinl( x ); +} diff --git a/src/s_sinf.c b/src/s_sinf.c index 6fed17a..97a2c92 100644 --- a/src/s_sinf.c +++ b/src/s_sinf.c @@ -35,7 +35,7 @@ s2pio2 = 2*M_PI_2, /* 0x400921FB, 0x54442D18 */ s3pio2 = 3*M_PI_2, /* 0x4012D97C, 0x7F3321D2 */ s4pio2 = 4*M_PI_2; /* 0x401921FB, 0x54442D18 */ -float +DLLEXPORT float sinf(float x) { double y; diff --git a/src/s_sinl.c b/src/s_sinl.c index 5ca150a..5c73c33 100644 --- a/src/s_sinl.c +++ b/src/s_sinl.c @@ -39,7 +39,7 @@ #error "Unsupported long double format" #endif -long double +DLLEXPORT long double sinl(long double x) { union IEEEl2bits z; diff --git a/src/s_tan.c b/src/s_tan.c index 3a6fbc4..4392776 100644 --- a/src/s_tan.c +++ b/src/s_tan.c @@ -50,7 +50,7 @@ #include "math_private.h" //#include "e_rem_pio2.c" -double +DLLEXPORT double tan(double x) { double y[2],z=0.0; diff --git a/src/s_tanf.c b/src/s_tanf.c index 7c3e0f1..ccf37a9 100644 --- a/src/s_tanf.c +++ b/src/s_tanf.c @@ -33,7 +33,7 @@ t2pio2 = 2*M_PI_2, /* 0x400921FB, 0x54442D18 */ t3pio2 = 3*M_PI_2, /* 0x4012D97C, 0x7F3321D2 */ t4pio2 = 4*M_PI_2; /* 0x401921FB, 0x54442D18 */ -float +DLLEXPORT float tanf(float x) { double y; diff --git a/src/s_tanh.c b/src/s_tanh.c index cb57d92..a400fd6 100644 --- a/src/s_tanh.c +++ b/src/s_tanh.c @@ -42,7 +42,7 @@ static const double one = 1.0, two = 2.0, tiny = 1.0e-300, huge = 1.0e300; -double +DLLEXPORT double tanh(double x) { double t,z; diff --git a/src/s_tanhf.c b/src/s_tanhf.c index e752612..85b5d18 100644 --- a/src/s_tanhf.c +++ b/src/s_tanhf.c @@ -20,7 +20,7 @@ #include "math_private.h" static const float one=1.0, two=2.0, tiny = 1.0e-30, huge = 1.0e30; -float +DLLEXPORT float tanhf(float x) { float t,z; diff --git a/src/s_tanl.c b/src/s_tanl.c index 0e5b5c2..ab3046a 100644 --- a/src/s_tanl.c +++ b/src/s_tanl.c @@ -45,7 +45,7 @@ #error "Unsupported long double format" #endif -long double +DLLEXPORT long double tanl(long double x) { union IEEEl2bits z; diff --git a/src/s_tgammaf.c b/src/s_tgammaf.c index 8567686..2d37d49 100644 --- a/src/s_tgammaf.c +++ b/src/s_tgammaf.c @@ -28,6 +28,7 @@ //__FBSDID("$FreeBSD: src/lib/msun/src/s_tgammaf.c,v 1.1 2008/02/18 17:27:10 das Exp $"); #include +#include "math_private.h" /* * We simply call tgamma() rather than bloating the math library with @@ -35,7 +36,7 @@ * essentially useless, since the function is superexponential and * floats have very limited range. */ -float +DLLEXPORT float tgammaf(float x) { diff --git a/src/s_trunc.c b/src/s_trunc.c index 3289f5c..f3d56ae 100644 --- a/src/s_trunc.c +++ b/src/s_trunc.c @@ -29,7 +29,7 @@ static const double huge = 1.0e300; -double +DLLEXPORT double trunc(double x) { int32_t i0,i1,j0; diff --git a/src/s_truncf.c b/src/s_truncf.c index 34b881d..470d894 100644 --- a/src/s_truncf.c +++ b/src/s_truncf.c @@ -27,7 +27,7 @@ static const float huge = 1.0e30F; -float +DLLEXPORT float truncf(float x) { int32_t i0,j0; diff --git a/src/s_truncl.c b/src/s_truncl.c index 2c237e5..02f3c61 100644 --- a/src/s_truncl.c +++ b/src/s_truncl.c @@ -25,6 +25,7 @@ #include #include +#include "math_private.h" #include #include "fpmath.h" @@ -38,7 +39,7 @@ static const long double huge = 1.0e300; static const float zero[] = { 0.0, -0.0 }; -long double +DLLEXPORT long double truncl(long double x) { union IEEEl2bits u = { .e = x }; diff --git a/src/w_cabs.c b/src/w_cabs.c index 22f52bb..c992e4f 100644 --- a/src/w_cabs.c +++ b/src/w_cabs.c @@ -11,8 +11,9 @@ #include #include #include +#include "math_private.h" -double +DLLEXPORT double cabs(double complex z) { return hypot(creal(z), cimag(z)); diff --git a/src/w_cabsf.c b/src/w_cabsf.c index 8df96b7..20405a3 100644 --- a/src/w_cabsf.c +++ b/src/w_cabsf.c @@ -14,7 +14,7 @@ static const char rcsid[] = #include #include "math_private.h" -float +DLLEXPORT float cabsf(z) float complex z; { diff --git a/src/w_cabsl.c b/src/w_cabsl.c index c3d8fbe..7ec31b3 100644 --- a/src/w_cabsl.c +++ b/src/w_cabsl.c @@ -12,8 +12,9 @@ #include #include +#include "math_private.h" -long double +DLLEXPORT long double cabsl(long double complex z) { return hypotl(creall(z), cimagl(z)); diff --git a/src/w_drem.c b/src/w_drem.c index 010f143..8a56f63 100644 --- a/src/w_drem.c +++ b/src/w_drem.c @@ -6,8 +6,9 @@ */ #include +#include "math_private.h" -double +DLLEXPORT double drem(x, y) double x, y; { diff --git a/src/w_dremf.c b/src/w_dremf.c index 37d584a..09a76b4 100644 --- a/src/w_dremf.c +++ b/src/w_dremf.c @@ -9,7 +9,7 @@ #include "openlibm.h" #include "math_private.h" -float +DLLEXPORT float dremf(float x, float y) { return remainderf(x, y);