Fix formatting, this is hard to explain, so I'll show one example.
- float ynf(int n, float x) /* wrapper ynf */ +float +ynf(int n, float x) /* wrapper ynf */ This is because the __STDC__ stuff was indented. Reviewed by: md5
This commit is contained in:
parent
027d0ad454
commit
59b19ff14a
@ -27,7 +27,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
#include "math.h"
|
||||
#include "math_private.h"
|
||||
|
||||
float __ieee754_gammaf_r(float x, int *signgamp)
|
||||
float
|
||||
__ieee754_gammaf_r(float x, int *signgamp)
|
||||
{
|
||||
return __ieee754_lgammaf_r(x,signgamp);
|
||||
}
|
||||
|
@ -49,7 +49,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
#include "math.h"
|
||||
#include "math_private.h"
|
||||
|
||||
double __ieee754_hypot(double x, double y)
|
||||
double
|
||||
__ieee754_hypot(double x, double y)
|
||||
{
|
||||
double a=x,b=y,t1,t2,y1,y2,w;
|
||||
int32_t j,k,ha,hb;
|
||||
|
@ -20,7 +20,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
#include "math.h"
|
||||
#include "math_private.h"
|
||||
|
||||
float __ieee754_hypotf(float x, float y)
|
||||
float
|
||||
__ieee754_hypotf(float x, float y)
|
||||
{
|
||||
float a=x,b=y,t1,t2,y1,y2,w;
|
||||
int32_t j,k,ha,hb;
|
||||
|
@ -81,7 +81,8 @@ S04 = 1.16614003333790000205e-09; /* 0x3E1408BC, 0xF4745D8F */
|
||||
|
||||
static const double zero = 0.0;
|
||||
|
||||
double __ieee754_j0(double x)
|
||||
double
|
||||
__ieee754_j0(double x)
|
||||
{
|
||||
double z, s,c,ss,cc,r,u,v;
|
||||
int32_t hx,ix;
|
||||
@ -141,7 +142,8 @@ v02 = 7.60068627350353253702e-05, /* 0x3F13ECBB, 0xF578C6C1 */
|
||||
v03 = 2.59150851840457805467e-07, /* 0x3E91642D, 0x7FF202FD */
|
||||
v04 = 4.41110311332675467403e-10; /* 0x3DFE5018, 0x3BD6D9EF */
|
||||
|
||||
double __ieee754_y0(double x)
|
||||
double
|
||||
__ieee754_y0(double x)
|
||||
{
|
||||
double z, s,c,ss,cc,u,v;
|
||||
int32_t hx,ix,lx;
|
||||
|
@ -39,7 +39,8 @@ S04 = 1.1661400734e-09; /* 0x30a045e8 */
|
||||
|
||||
static const float zero = 0.0;
|
||||
|
||||
float __ieee754_j0f(float x)
|
||||
float
|
||||
__ieee754_j0f(float x)
|
||||
{
|
||||
float z, s,c,ss,cc,r,u,v;
|
||||
int32_t hx,ix;
|
||||
@ -99,7 +100,8 @@ v02 = 7.6006865129e-05, /* 0x389f65e0 */
|
||||
v03 = 2.5915085189e-07, /* 0x348b216c */
|
||||
v04 = 4.4111031494e-10; /* 0x2ff280c2 */
|
||||
|
||||
float __ieee754_y0f(float x)
|
||||
float
|
||||
__ieee754_y0f(float x)
|
||||
{
|
||||
float z, s,c,ss,cc,u,v;
|
||||
int32_t hx,ix;
|
||||
|
@ -82,7 +82,8 @@ s05 = 1.23542274426137913908e-11; /* 0x3DAB2ACF, 0xCFB97ED8 */
|
||||
|
||||
static const double zero = 0.0;
|
||||
|
||||
double __ieee754_j1(double x)
|
||||
double
|
||||
__ieee754_j1(double x)
|
||||
{
|
||||
double z, s,c,ss,cc,r,u,v,y;
|
||||
int32_t hx,ix;
|
||||
@ -138,7 +139,8 @@ static const double V0[5] = {
|
||||
1.66559246207992079114e-11, /* 0x3DB25039, 0xDACA772A */
|
||||
};
|
||||
|
||||
double __ieee754_y1(double x)
|
||||
double
|
||||
__ieee754_y1(double x)
|
||||
{
|
||||
double z, s,c,ss,cc,u,v;
|
||||
int32_t hx,ix,lx;
|
||||
|
@ -40,7 +40,8 @@ s05 = 1.2354227016e-11; /* 0x2d59567e */
|
||||
|
||||
static const float zero = 0.0;
|
||||
|
||||
float __ieee754_j1f(float x)
|
||||
float
|
||||
__ieee754_j1f(float x)
|
||||
{
|
||||
float z, s,c,ss,cc,r,u,v,y;
|
||||
int32_t hx,ix;
|
||||
@ -96,7 +97,8 @@ static const float V0[5] = {
|
||||
1.6655924903e-11, /* 0x2d9281cf */
|
||||
};
|
||||
|
||||
float __ieee754_y1f(float x)
|
||||
float
|
||||
__ieee754_y1f(float x)
|
||||
{
|
||||
float z, s,c,ss,cc,u,v;
|
||||
int32_t hx,ix;
|
||||
|
@ -50,7 +50,8 @@ one = 1.00000000000000000000e+00; /* 0x3FF00000, 0x00000000 */
|
||||
|
||||
static const double zero = 0.00000000000000000000e+00;
|
||||
|
||||
double __ieee754_jn(int n, double x)
|
||||
double
|
||||
__ieee754_jn(int n, double x)
|
||||
{
|
||||
int32_t i,hx,ix,lx, sgn;
|
||||
double a, b, temp, di;
|
||||
@ -205,7 +206,8 @@ static const double zero = 0.00000000000000000000e+00;
|
||||
if(sgn==1) return -b; else return b;
|
||||
}
|
||||
|
||||
double __ieee754_yn(int n, double x)
|
||||
double
|
||||
__ieee754_yn(int n, double x)
|
||||
{
|
||||
int32_t i,hx,ix,lx;
|
||||
int32_t sign;
|
||||
|
@ -27,7 +27,8 @@ one = 1.0000000000e+00; /* 0x3F800000 */
|
||||
|
||||
static const float zero = 0.0000000000e+00;
|
||||
|
||||
float __ieee754_jnf(int n, float x)
|
||||
float
|
||||
__ieee754_jnf(int n, float x)
|
||||
{
|
||||
int32_t i,hx,ix, sgn;
|
||||
float a, b, temp, di;
|
||||
@ -159,7 +160,8 @@ static const float zero = 0.0000000000e+00;
|
||||
if(sgn==1) return -b; else return b;
|
||||
}
|
||||
|
||||
float __ieee754_ynf(int n, float x)
|
||||
float
|
||||
__ieee754_ynf(int n, float x)
|
||||
{
|
||||
int32_t i,hx,ix,ib;
|
||||
int32_t sign;
|
||||
|
@ -25,7 +25,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
|
||||
extern int signgam;
|
||||
|
||||
double __ieee754_lgamma(double x)
|
||||
double
|
||||
__ieee754_lgamma(double x)
|
||||
{
|
||||
return __ieee754_lgamma_r(x,&signgam);
|
||||
}
|
||||
|
@ -199,7 +199,8 @@ static const double zero= 0.00000000000000000000e+00;
|
||||
}
|
||||
|
||||
|
||||
double __ieee754_lgamma_r(double x, int *signgamp)
|
||||
double
|
||||
__ieee754_lgamma_r(double x, int *signgamp)
|
||||
{
|
||||
double t,y,z,nadj,p,p1,p2,p3,q,r,w;
|
||||
int i,hx,lx,ix;
|
||||
|
@ -28,7 +28,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
|
||||
extern int signgam;
|
||||
|
||||
float __ieee754_lgammaf(float x)
|
||||
float
|
||||
__ieee754_lgammaf(float x)
|
||||
{
|
||||
return __ieee754_lgammaf_r(x,&signgam);
|
||||
}
|
||||
|
@ -135,7 +135,8 @@ static const float zero= 0.0000000000e+00;
|
||||
}
|
||||
|
||||
|
||||
float __ieee754_lgammaf_r(float x, int *signgamp)
|
||||
float
|
||||
__ieee754_lgammaf_r(float x, int *signgamp)
|
||||
{
|
||||
float t,y,z,nadj,p,p1,p2,p3,q,r,w;
|
||||
int i,hx,ix;
|
||||
|
@ -82,7 +82,8 @@ Lg7 = 1.479819860511658591e-01; /* 3FC2F112 DF3E5244 */
|
||||
|
||||
static const double zero = 0.0;
|
||||
|
||||
double __generic___ieee754_log(double x)
|
||||
double
|
||||
__generic___ieee754_log(double x)
|
||||
{
|
||||
double hfsq,f,s,z,R,w,t1,t2,dk;
|
||||
int32_t k,hx,i,j;
|
||||
|
@ -58,7 +58,8 @@ log10_2lo = 3.69423907715893078616e-13; /* 0x3D59FEF3, 0x11F12B36 */
|
||||
|
||||
static const double zero = 0.0;
|
||||
|
||||
double __generic___ieee754_log10(double x)
|
||||
double
|
||||
__generic___ieee754_log10(double x)
|
||||
{
|
||||
double y,z;
|
||||
int32_t i,k,hx;
|
||||
|
@ -28,7 +28,8 @@ log10_2lo = 7.9034151668e-07; /* 0x355427db */
|
||||
|
||||
static const float zero = 0.0;
|
||||
|
||||
float __ieee754_log10f(float x)
|
||||
float
|
||||
__ieee754_log10f(float x)
|
||||
{
|
||||
float y,z;
|
||||
int32_t i,k,hx;
|
||||
|
@ -34,7 +34,8 @@ Lg7 = 1.4798198640e-01; /* 3E178897 */
|
||||
|
||||
static const float zero = 0.0;
|
||||
|
||||
float __ieee754_logf(float x)
|
||||
float
|
||||
__ieee754_logf(float x)
|
||||
{
|
||||
float hfsq,f,s,z,R,w,t1,t2,dk;
|
||||
int32_t k,ix,i,j;
|
||||
|
@ -95,7 +95,8 @@ 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 __ieee754_pow(double x, double y)
|
||||
double
|
||||
__ieee754_pow(double x, double y)
|
||||
{
|
||||
double z,ax,z_h,z_l,p_h,p_l;
|
||||
double y1,t1,t2,r,s,t,u,v,w;
|
||||
|
@ -53,7 +53,8 @@ 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 __ieee754_powf(float x, float y)
|
||||
float
|
||||
__ieee754_powf(float x, float y)
|
||||
{
|
||||
float z,ax,z_h,z_l,p_h,p_l;
|
||||
float y1,t1,t2,r,s,t,u,v,w;
|
||||
|
@ -29,7 +29,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
static const double zero = 0.0;
|
||||
|
||||
|
||||
double __generic___ieee754_remainder(double x, double p)
|
||||
double
|
||||
__generic___ieee754_remainder(double x, double p)
|
||||
{
|
||||
int32_t hx,hp;
|
||||
u_int32_t sx,lx,lp;
|
||||
|
@ -23,7 +23,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
static const float zero = 0.0;
|
||||
|
||||
|
||||
float __ieee754_remainderf(float x, float p)
|
||||
float
|
||||
__ieee754_remainderf(float x, float p)
|
||||
{
|
||||
int32_t hx,hp;
|
||||
u_int32_t sx;
|
||||
|
@ -24,9 +24,11 @@ static char rcsid[] = "$FreeBSD$";
|
||||
#include "math_private.h"
|
||||
|
||||
#ifdef _SCALB_INT
|
||||
double __generic___ieee754_scalb(double x, int fn)
|
||||
double
|
||||
__generic___ieee754_scalb(double x, int fn)
|
||||
#else
|
||||
double __generic___ieee754_scalb(double x, double fn)
|
||||
double
|
||||
__generic___ieee754_scalb(double x, double fn)
|
||||
#endif
|
||||
{
|
||||
#ifdef _SCALB_INT
|
||||
|
@ -21,9 +21,11 @@ static char rcsid[] = "$FreeBSD$";
|
||||
#include "math_private.h"
|
||||
|
||||
#ifdef _SCALB_INT
|
||||
float __ieee754_scalbf(float x, int fn)
|
||||
float
|
||||
__ieee754_scalbf(float x, int fn)
|
||||
#else
|
||||
float __ieee754_scalbf(float x, float fn)
|
||||
float
|
||||
__ieee754_scalbf(float x, float fn)
|
||||
#endif
|
||||
{
|
||||
#ifdef _SCALB_INT
|
||||
|
@ -37,7 +37,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
|
||||
static const double one = 1.0, shuge = 1.0e307;
|
||||
|
||||
double __ieee754_sinh(double x)
|
||||
double
|
||||
__ieee754_sinh(double x)
|
||||
{
|
||||
double t,w,h;
|
||||
int32_t ix,jx;
|
||||
|
@ -22,7 +22,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
|
||||
static const float one = 1.0, shuge = 1.0e37;
|
||||
|
||||
float __ieee754_sinhf(float x)
|
||||
float
|
||||
__ieee754_sinhf(float x)
|
||||
{
|
||||
float t,w,h;
|
||||
int32_t ix,jx;
|
||||
|
@ -89,7 +89,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
|
||||
static const double one = 1.0, tiny=1.0e-300;
|
||||
|
||||
double __generic___ieee754_sqrt(double x)
|
||||
double
|
||||
__generic___ieee754_sqrt(double x)
|
||||
{
|
||||
double z;
|
||||
int32_t sign = (int)0x80000000;
|
||||
|
@ -22,7 +22,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
|
||||
static const float one = 1.0, tiny=1.0e-30;
|
||||
|
||||
float __ieee754_sqrtf(float x)
|
||||
float
|
||||
__ieee754_sqrtf(float x)
|
||||
{
|
||||
float z;
|
||||
int32_t sign = (int)0x80000000;
|
||||
|
@ -61,7 +61,8 @@ C4 = -2.75573143513906633035e-07, /* 0xBE927E4F, 0x809C52AD */
|
||||
C5 = 2.08757232129817482790e-09, /* 0x3E21EE9E, 0xBDB4B1C4 */
|
||||
C6 = -1.13596475577881948265e-11; /* 0xBDA8FAE9, 0xBE8838D4 */
|
||||
|
||||
double __kernel_cos(double x, double y)
|
||||
double
|
||||
__kernel_cos(double x, double y)
|
||||
{
|
||||
double a,hz,z,r,qx;
|
||||
int32_t ix;
|
||||
|
@ -29,7 +29,8 @@ C4 = -2.7557314297e-07, /* 0xb493f27c */
|
||||
C5 = 2.0875723372e-09, /* 0x310f74f6 */
|
||||
C6 = -1.1359647598e-11; /* 0xad47d74e */
|
||||
|
||||
float __kernel_cosf(float x, float y)
|
||||
float
|
||||
__kernel_cosf(float x, float y)
|
||||
{
|
||||
float a,hz,z,r,qx;
|
||||
int32_t ix;
|
||||
|
@ -54,7 +54,8 @@ S4 = 2.75573137070700676789e-06, /* 0x3EC71DE3, 0x57B1FE7D */
|
||||
S5 = -2.50507602534068634195e-08, /* 0xBE5AE5E6, 0x8A2B9CEB */
|
||||
S6 = 1.58969099521155010221e-10; /* 0x3DE5D93A, 0x5ACFD57C */
|
||||
|
||||
double __kernel_sin(double x, double y, int iy)
|
||||
double
|
||||
__kernel_sin(double x, double y, int iy)
|
||||
{
|
||||
double z,r,v;
|
||||
int32_t ix;
|
||||
|
@ -29,7 +29,8 @@ S4 = 2.7557314297e-06, /* 0x3638ef1b */
|
||||
S5 = -2.5050759689e-08, /* 0xb2d72f34 */
|
||||
S6 = 1.5896910177e-10; /* 0x2f2ec9d3 */
|
||||
|
||||
float __kernel_sinf(float x, float y, int iy)
|
||||
float
|
||||
__kernel_sinf(float x, float y, int iy)
|
||||
{
|
||||
float z,r,v;
|
||||
int32_t ix;
|
||||
|
@ -77,7 +77,8 @@ static const double zero = 0.0; /* used as const */
|
||||
*/
|
||||
|
||||
|
||||
double __kernel_standard(double x, double y, int type)
|
||||
double
|
||||
__kernel_standard(double x, double y, int type)
|
||||
{
|
||||
struct exception exc;
|
||||
#ifndef HUGE_VAL /* this is the only routine that uses HUGE_VAL */
|
||||
|
@ -70,7 +70,8 @@ T[] = {
|
||||
2.59073051863633712884e-05, /* 0x3EFB2A70, 0x74BF7AD4 */
|
||||
};
|
||||
|
||||
double __kernel_tan(double x, double y, int iy)
|
||||
double
|
||||
__kernel_tan(double x, double y, int iy)
|
||||
{
|
||||
double z,r,v,w,s;
|
||||
int32_t ix,hx;
|
||||
|
@ -39,7 +39,8 @@ T[] = {
|
||||
2.5907305826e-05, /* 0x37d95384 */
|
||||
};
|
||||
|
||||
float __kernel_tanf(float x, float y, int iy)
|
||||
float
|
||||
__kernel_tanf(float x, float y, int iy)
|
||||
{
|
||||
float z,r,v,w,s;
|
||||
int32_t ix,hx;
|
||||
|
@ -33,7 +33,8 @@ one = 1.00000000000000000000e+00, /* 0x3FF00000, 0x00000000 */
|
||||
ln2 = 6.93147180559945286227e-01, /* 0x3FE62E42, 0xFEFA39EF */
|
||||
huge= 1.00000000000000000000e+300;
|
||||
|
||||
double asinh(double x)
|
||||
double
|
||||
asinh(double x)
|
||||
{
|
||||
double t,w;
|
||||
int32_t hx,ix;
|
||||
|
@ -25,7 +25,8 @@ one = 1.0000000000e+00, /* 0x3F800000 */
|
||||
ln2 = 6.9314718246e-01, /* 0x3f317218 */
|
||||
huge= 1.0000000000e+30;
|
||||
|
||||
float asinhf(float x)
|
||||
float
|
||||
asinhf(float x)
|
||||
{
|
||||
float t,w;
|
||||
int32_t hx,ix;
|
||||
|
@ -69,7 +69,8 @@ static const double aT[] = {
|
||||
one = 1.0,
|
||||
huge = 1.0e300;
|
||||
|
||||
double __generic_atan(double x)
|
||||
double
|
||||
__generic_atan(double x)
|
||||
{
|
||||
double w,s1,s2,z;
|
||||
int32_t ix,hx,id;
|
||||
|
@ -52,7 +52,8 @@ static const float aT[] = {
|
||||
one = 1.0,
|
||||
huge = 1.0e30;
|
||||
|
||||
float atanf(float x)
|
||||
float
|
||||
atanf(float x)
|
||||
{
|
||||
float w,s1,s2,z;
|
||||
int32_t ix,hx,id;
|
||||
|
@ -31,7 +31,8 @@ E = 1.41428571428571436819e+00, /* 99/70 = 0x3FF6A0EA, 0x0EA0EA0F */
|
||||
F = 1.60714285714285720630e+00, /* 45/28 = 0x3FF9B6DB, 0x6DB6DB6E */
|
||||
G = 3.57142857142857150787e-01; /* 5/14 = 0x3FD6DB6D, 0xB6DB6DB7 */
|
||||
|
||||
double cbrt(double x)
|
||||
double
|
||||
cbrt(double x)
|
||||
{
|
||||
int32_t hx;
|
||||
double r,s,t=0.0,w;
|
||||
|
@ -34,7 +34,8 @@ E = 1.4142856598e+00, /* 99/70 = 0x3fb50750 */
|
||||
F = 1.6071428061e+00, /* 45/28 = 0x3fcdb6db */
|
||||
G = 3.5714286566e-01; /* 5/14 = 0x3eb6db6e */
|
||||
|
||||
float cbrtf(float x)
|
||||
float
|
||||
cbrtf(float x)
|
||||
{
|
||||
float r,s,t;
|
||||
int32_t hx;
|
||||
|
@ -28,7 +28,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
|
||||
static const double huge = 1.0e300;
|
||||
|
||||
double __generic_ceil(double x)
|
||||
double
|
||||
__generic_ceil(double x)
|
||||
{
|
||||
int32_t i0,i1,j0;
|
||||
u_int32_t i,j;
|
||||
|
@ -22,7 +22,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
|
||||
static const float huge = 1.0e30;
|
||||
|
||||
float ceilf(float x)
|
||||
float
|
||||
ceilf(float x)
|
||||
{
|
||||
int32_t i0,j0;
|
||||
u_int32_t i;
|
||||
|
@ -23,7 +23,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
#include "math.h"
|
||||
#include "math_private.h"
|
||||
|
||||
double __generic_copysign(double x, double y)
|
||||
double
|
||||
__generic_copysign(double x, double y)
|
||||
{
|
||||
u_int32_t hx,hy;
|
||||
GET_HIGH_WORD(hx,x);
|
||||
|
@ -26,7 +26,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
#include "math.h"
|
||||
#include "math_private.h"
|
||||
|
||||
float copysignf(float x, float y)
|
||||
float
|
||||
copysignf(float x, float y)
|
||||
{
|
||||
u_int32_t ix,iy;
|
||||
GET_FLOAT_WORD(ix,x);
|
||||
|
@ -48,7 +48,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
#include "math.h"
|
||||
#include "math_private.h"
|
||||
|
||||
double __generic_cos(double x)
|
||||
double
|
||||
__generic_cos(double x)
|
||||
{
|
||||
double y[2],z=0.0;
|
||||
int32_t n, ix;
|
||||
|
@ -22,7 +22,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
|
||||
static const float one=1.0;
|
||||
|
||||
float cosf(float x)
|
||||
float
|
||||
cosf(float x)
|
||||
{
|
||||
float y[2],z=0.0;
|
||||
int32_t n,ix;
|
||||
|
@ -187,7 +187,8 @@ sb5 = 2.55305040643316442583e+03, /* 0x40A3F219, 0xCEDF3BE6 */
|
||||
sb6 = 4.74528541206955367215e+02, /* 0x407DA874, 0xE79FE763 */
|
||||
sb7 = -2.24409524465858183362e+01; /* 0xC03670E2, 0x42712D62 */
|
||||
|
||||
double erf(double x)
|
||||
double
|
||||
erf(double x)
|
||||
{
|
||||
int32_t hx,ix,i;
|
||||
double R,S,P,Q,s,y,z,r;
|
||||
@ -238,7 +239,8 @@ sb7 = -2.24409524465858183362e+01; /* 0xC03670E2, 0x42712D62 */
|
||||
if(hx>=0) return one-r/x; else return r/x-one;
|
||||
}
|
||||
|
||||
double erfc(double x)
|
||||
double
|
||||
erfc(double x)
|
||||
{
|
||||
int32_t hx,ix;
|
||||
double R,S,P,Q,s,y,z,r;
|
||||
|
@ -95,7 +95,8 @@ sb5 = 2.5530502930e+03, /* 0x451f90ce */
|
||||
sb6 = 4.7452853394e+02, /* 0x43ed43a7 */
|
||||
sb7 = -2.2440952301e+01; /* 0xc1b38712 */
|
||||
|
||||
float erff(float x)
|
||||
float
|
||||
erff(float x)
|
||||
{
|
||||
int32_t hx,ix,i;
|
||||
float R,S,P,Q,s,y,z,r;
|
||||
@ -147,7 +148,8 @@ sb7 = -2.2440952301e+01; /* 0xc1b38712 */
|
||||
if(hx>=0) return one-r/x; else return r/x-one;
|
||||
}
|
||||
|
||||
float erfcf(float x)
|
||||
float
|
||||
erfcf(float x)
|
||||
{
|
||||
int32_t hx,ix;
|
||||
float R,S,P,Q,s,y,z,r;
|
||||
|
@ -127,7 +127,8 @@ Q3 = -7.93650757867487942473e-05, /* BF14CE19 9EAADBB7 */
|
||||
Q4 = 4.00821782732936239552e-06, /* 3ED0CFCA 86E65239 */
|
||||
Q5 = -2.01099218183624371326e-07; /* BE8AFDB7 6E09C32D */
|
||||
|
||||
double expm1(double x)
|
||||
double
|
||||
expm1(double x)
|
||||
{
|
||||
double y,hi,lo,c,t,e,hxs,hfx,r1;
|
||||
int32_t k,xsb;
|
||||
|
@ -35,7 +35,8 @@ Q3 = -7.9365076090e-05, /* 0xb8a670cd */
|
||||
Q4 = 4.0082177293e-06, /* 0x36867e54 */
|
||||
Q5 = -2.0109921195e-07; /* 0xb457edbb */
|
||||
|
||||
float expm1f(float x)
|
||||
float
|
||||
expm1f(float x)
|
||||
{
|
||||
float y,hi,lo,c,t,e,hxs,hfx,r1;
|
||||
int32_t k,xsb;
|
||||
|
@ -21,7 +21,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
#include "math.h"
|
||||
#include "math_private.h"
|
||||
|
||||
double fabs(double x)
|
||||
double
|
||||
fabs(double x)
|
||||
{
|
||||
u_int32_t high;
|
||||
GET_HIGH_WORD(high,x);
|
||||
|
@ -24,7 +24,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
#include "math.h"
|
||||
#include "math_private.h"
|
||||
|
||||
float fabsf(float x)
|
||||
float
|
||||
fabsf(float x)
|
||||
{
|
||||
u_int32_t ix;
|
||||
GET_FLOAT_WORD(ix,x);
|
||||
|
@ -28,7 +28,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
|
||||
static const double huge = 1.0e300;
|
||||
|
||||
double __generic_floor(double x)
|
||||
double
|
||||
__generic_floor(double x)
|
||||
{
|
||||
int32_t i0,i1,j0;
|
||||
u_int32_t i,j;
|
||||
|
@ -31,7 +31,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
|
||||
static const float huge = 1.0e30;
|
||||
|
||||
float floorf(float x)
|
||||
float
|
||||
floorf(float x)
|
||||
{
|
||||
int32_t i0,j0;
|
||||
u_int32_t i;
|
||||
|
@ -30,7 +30,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
static const double
|
||||
two54 = 1.80143985094819840000e+16; /* 0x43500000, 0x00000000 */
|
||||
|
||||
double frexp(double x, int *eptr)
|
||||
double
|
||||
frexp(double x, int *eptr)
|
||||
{
|
||||
int32_t hx, ix, lx;
|
||||
EXTRACT_WORDS(hx,lx,x);
|
||||
|
@ -23,7 +23,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
static const float
|
||||
two25 = 3.3554432000e+07; /* 0x4c000000 */
|
||||
|
||||
float frexpf(float x, int *eptr)
|
||||
float
|
||||
frexpf(float x, int *eptr)
|
||||
{
|
||||
int32_t hx,ix;
|
||||
GET_FLOAT_WORD(hx,x);
|
||||
|
@ -18,7 +18,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
#include "math_private.h"
|
||||
#include <errno.h>
|
||||
|
||||
double ldexp(double value, int exp)
|
||||
double
|
||||
ldexp(double value, int exp)
|
||||
{
|
||||
if(!finite(value)||value==0.0) return value;
|
||||
value = scalbn(value,exp);
|
||||
|
@ -21,7 +21,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
#include "math_private.h"
|
||||
#include <errno.h>
|
||||
|
||||
float ldexpf(float value, int exp)
|
||||
float
|
||||
ldexpf(float value, int exp)
|
||||
{
|
||||
if(!finitef(value)||value==(float)0.0) return value;
|
||||
value = scalbnf(value,exp);
|
||||
|
@ -96,7 +96,8 @@ Lp7 = 1.479819860511658591e-01; /* 3FC2F112 DF3E5244 */
|
||||
|
||||
static const double zero = 0.0;
|
||||
|
||||
double log1p(double x)
|
||||
double
|
||||
log1p(double x)
|
||||
{
|
||||
double hfsq,f,c,s,z,R,u;
|
||||
int32_t k,hx,hu,ax;
|
||||
|
@ -34,7 +34,8 @@ Lp7 = 1.4798198640e-01; /* 3E178897 */
|
||||
|
||||
static const float zero = 0.0;
|
||||
|
||||
float log1pf(float x)
|
||||
float
|
||||
log1pf(float x)
|
||||
{
|
||||
float hfsq,f,c,s,z,R,u;
|
||||
int32_t k,hx,hu,ax;
|
||||
|
@ -23,7 +23,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
#include "math.h"
|
||||
#include "math_private.h"
|
||||
|
||||
double __generic_logb(double x)
|
||||
double
|
||||
__generic_logb(double x)
|
||||
{
|
||||
int32_t lx,ix;
|
||||
EXTRACT_WORDS(ix,lx,x);
|
||||
|
@ -20,7 +20,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
#include "math.h"
|
||||
#include "math_private.h"
|
||||
|
||||
float logbf(float x)
|
||||
float
|
||||
logbf(float x)
|
||||
{
|
||||
int32_t ix;
|
||||
GET_FLOAT_WORD(ix,x);
|
||||
|
@ -29,7 +29,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
|
||||
static const double one = 1.0;
|
||||
|
||||
double modf(double x, double *iptr)
|
||||
double
|
||||
modf(double x, double *iptr)
|
||||
{
|
||||
int32_t i0,i1,j0;
|
||||
u_int32_t i;
|
||||
|
@ -22,7 +22,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
|
||||
static const float one = 1.0;
|
||||
|
||||
float modff(float x, float *iptr)
|
||||
float
|
||||
modff(float x, float *iptr)
|
||||
{
|
||||
int32_t i0,j0;
|
||||
u_int32_t i;
|
||||
|
@ -24,7 +24,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
#include "math.h"
|
||||
#include "math_private.h"
|
||||
|
||||
double nextafter(double x, double y)
|
||||
double
|
||||
nextafter(double x, double y)
|
||||
{
|
||||
int32_t hx,hy,ix,iy;
|
||||
u_int32_t lx,ly;
|
||||
|
@ -20,7 +20,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
#include "math.h"
|
||||
#include "math_private.h"
|
||||
|
||||
float nextafterf(float x, float y)
|
||||
float
|
||||
nextafterf(float x, float y)
|
||||
{
|
||||
int32_t hx,hy,ix,iy;
|
||||
|
||||
|
@ -40,7 +40,8 @@ TWO52[2]={
|
||||
-4.50359962737049600000e+15, /* 0xC3300000, 0x00000000 */
|
||||
};
|
||||
|
||||
double __generic_rint(double x)
|
||||
double
|
||||
__generic_rint(double x)
|
||||
{
|
||||
int32_t i0,j0,sx;
|
||||
u_int32_t i,i1;
|
||||
|
@ -32,7 +32,8 @@ TWO23[2]={
|
||||
-8.3886080000e+06, /* 0xcb000000 */
|
||||
};
|
||||
|
||||
float rintf(float x)
|
||||
float
|
||||
rintf(float x)
|
||||
{
|
||||
int32_t i0,j0,sx;
|
||||
u_int32_t i,i1;
|
||||
|
@ -30,7 +30,8 @@ twom54 = 5.55111512312578270212e-17, /* 0x3C900000, 0x00000000 */
|
||||
huge = 1.0e+300,
|
||||
tiny = 1.0e-300;
|
||||
|
||||
double __generic_scalbn (double x, int n)
|
||||
double
|
||||
__generic_scalbn (double x, int n)
|
||||
{
|
||||
int32_t k,hx,lx;
|
||||
EXTRACT_WORDS(hx,lx,x);
|
||||
|
@ -26,7 +26,8 @@ twom25 = 2.9802322388e-08, /* 0x33000000 */
|
||||
huge = 1.0e+30,
|
||||
tiny = 1.0e-30;
|
||||
|
||||
float scalbnf (float x, int n)
|
||||
float
|
||||
scalbnf (float x, int n)
|
||||
{
|
||||
int32_t k,ix;
|
||||
GET_FLOAT_WORD(ix,x);
|
||||
|
@ -23,7 +23,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
#include "math.h"
|
||||
#include "math_private.h"
|
||||
|
||||
double __generic_significand(double x)
|
||||
double
|
||||
__generic_significand(double x)
|
||||
{
|
||||
return __ieee754_scalb(x,(double) -ilogb(x));
|
||||
}
|
||||
|
@ -20,7 +20,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
#include "math.h"
|
||||
#include "math_private.h"
|
||||
|
||||
float significandf(float x)
|
||||
float
|
||||
significandf(float x)
|
||||
{
|
||||
return __ieee754_scalbf(x,(float) -ilogbf(x));
|
||||
}
|
||||
|
@ -48,7 +48,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
#include "math.h"
|
||||
#include "math_private.h"
|
||||
|
||||
double __generic_sin(double x)
|
||||
double
|
||||
__generic_sin(double x)
|
||||
{
|
||||
double y[2],z=0.0;
|
||||
int32_t n, ix;
|
||||
|
@ -20,7 +20,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
#include "math.h"
|
||||
#include "math_private.h"
|
||||
|
||||
float sinf(float x)
|
||||
float
|
||||
sinf(float x)
|
||||
{
|
||||
float y[2],z=0.0;
|
||||
int32_t n, ix;
|
||||
|
@ -47,7 +47,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
#include "math.h"
|
||||
#include "math_private.h"
|
||||
|
||||
double __generic_tan(double x)
|
||||
double
|
||||
__generic_tan(double x)
|
||||
{
|
||||
double y[2],z=0.0;
|
||||
int32_t n, ix;
|
||||
|
@ -20,7 +20,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
#include "math.h"
|
||||
#include "math_private.h"
|
||||
|
||||
float tanf(float x)
|
||||
float
|
||||
tanf(float x)
|
||||
{
|
||||
float y[2],z=0.0;
|
||||
int32_t n, ix;
|
||||
|
@ -43,7 +43,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
|
||||
static const double one=1.0, two=2.0, tiny = 1.0e-300;
|
||||
|
||||
double tanh(double x)
|
||||
double
|
||||
tanh(double x)
|
||||
{
|
||||
double t,z;
|
||||
int32_t jx,ix;
|
||||
|
@ -22,7 +22,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
|
||||
static const float one=1.0, two=2.0, tiny = 1.0e-30;
|
||||
|
||||
float tanhf(float x)
|
||||
float
|
||||
tanhf(float x)
|
||||
{
|
||||
float t,z;
|
||||
int32_t jx,ix;
|
||||
|
@ -22,7 +22,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
#include "math_private.h"
|
||||
|
||||
|
||||
double acos(double x) /* wrapper acos */
|
||||
double
|
||||
acos(double x) /* wrapper acos */
|
||||
{
|
||||
#ifdef _IEEE_LIBM
|
||||
return __ieee754_acos(x);
|
||||
|
@ -25,7 +25,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
#include "math_private.h"
|
||||
|
||||
|
||||
float acosf(float x) /* wrapper acosf */
|
||||
float
|
||||
acosf(float x) /* wrapper acosf */
|
||||
{
|
||||
#ifdef _IEEE_LIBM
|
||||
return __ieee754_acosf(x);
|
||||
|
@ -21,7 +21,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
#include "math.h"
|
||||
#include "math_private.h"
|
||||
|
||||
double acosh(double x) /* wrapper acosh */
|
||||
double
|
||||
acosh(double x) /* wrapper acosh */
|
||||
{
|
||||
#ifdef _IEEE_LIBM
|
||||
return __ieee754_acosh(x);
|
||||
|
@ -25,7 +25,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
#include "math.h"
|
||||
#include "math_private.h"
|
||||
|
||||
float acoshf(float x) /* wrapper acoshf */
|
||||
float
|
||||
acoshf(float x) /* wrapper acoshf */
|
||||
{
|
||||
#ifdef _IEEE_LIBM
|
||||
return __ieee754_acoshf(x);
|
||||
|
@ -23,7 +23,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
#include "math_private.h"
|
||||
|
||||
|
||||
double asin(double x) /* wrapper asin */
|
||||
double
|
||||
asin(double x) /* wrapper asin */
|
||||
{
|
||||
#ifdef _IEEE_LIBM
|
||||
return __ieee754_asin(x);
|
||||
|
@ -26,7 +26,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
#include "math_private.h"
|
||||
|
||||
|
||||
float asinf(float x) /* wrapper asinf */
|
||||
float
|
||||
asinf(float x) /* wrapper asinf */
|
||||
{
|
||||
#ifdef _IEEE_LIBM
|
||||
return __ieee754_asinf(x);
|
||||
|
@ -22,7 +22,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
#include "math_private.h"
|
||||
|
||||
|
||||
double atan2(double y, double x) /* wrapper atan2 */
|
||||
double
|
||||
atan2(double y, double x) /* wrapper atan2 */
|
||||
{
|
||||
#ifdef _IEEE_LIBM
|
||||
return __ieee754_atan2(y,x);
|
||||
|
@ -25,7 +25,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
#include "math_private.h"
|
||||
|
||||
|
||||
float atan2f(float y, float x) /* wrapper atan2f */
|
||||
float
|
||||
atan2f(float y, float x) /* wrapper atan2f */
|
||||
{
|
||||
#ifdef _IEEE_LIBM
|
||||
return __ieee754_atan2f(y,x);
|
||||
|
@ -22,7 +22,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
#include "math_private.h"
|
||||
|
||||
|
||||
double atanh(double x) /* wrapper atanh */
|
||||
double
|
||||
atanh(double x) /* wrapper atanh */
|
||||
{
|
||||
#ifdef _IEEE_LIBM
|
||||
return __ieee754_atanh(x);
|
||||
|
@ -25,7 +25,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
#include "math_private.h"
|
||||
|
||||
|
||||
float atanhf(float x) /* wrapper atanhf */
|
||||
float
|
||||
atanhf(float x) /* wrapper atanhf */
|
||||
{
|
||||
#ifdef _IEEE_LIBM
|
||||
return __ieee754_atanhf(x);
|
||||
|
@ -21,7 +21,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
#include "math.h"
|
||||
#include "math_private.h"
|
||||
|
||||
double cosh(double x) /* wrapper cosh */
|
||||
double
|
||||
cosh(double x) /* wrapper cosh */
|
||||
{
|
||||
#ifdef _IEEE_LIBM
|
||||
return __ieee754_cosh(x);
|
||||
|
@ -24,7 +24,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
#include "math.h"
|
||||
#include "math_private.h"
|
||||
|
||||
float coshf(float x) /* wrapper coshf */
|
||||
float
|
||||
coshf(float x) /* wrapper coshf */
|
||||
{
|
||||
#ifdef _IEEE_LIBM
|
||||
return __ieee754_coshf(x);
|
||||
|
@ -25,7 +25,8 @@ static const double
|
||||
o_threshold= 7.09782712893383973096e+02, /* 0x40862E42, 0xFEFA39EF */
|
||||
u_threshold= -7.45133219101941108420e+02; /* 0xc0874910, 0xD52D3051 */
|
||||
|
||||
double exp(double x) /* wrapper exp */
|
||||
double
|
||||
exp(double x) /* wrapper exp */
|
||||
{
|
||||
#ifdef _IEEE_LIBM
|
||||
return __ieee754_exp(x);
|
||||
|
@ -28,7 +28,8 @@ static const float
|
||||
o_threshold= 8.8721679688e+01, /* 0x42b17180 */
|
||||
u_threshold= -1.0397208405e+02; /* 0xc2cff1b5 */
|
||||
|
||||
float expf(float x) /* wrapper expf */
|
||||
float
|
||||
expf(float x) /* wrapper expf */
|
||||
{
|
||||
#ifdef _IEEE_LIBM
|
||||
return __ieee754_expf(x);
|
||||
|
@ -22,7 +22,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
#include "math_private.h"
|
||||
|
||||
|
||||
double fmod(double x, double y) /* wrapper fmod */
|
||||
double
|
||||
fmod(double x, double y) /* wrapper fmod */
|
||||
{
|
||||
#ifdef _IEEE_LIBM
|
||||
return __ieee754_fmod(x,y);
|
||||
|
@ -25,7 +25,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
#include "math_private.h"
|
||||
|
||||
|
||||
float fmodf(float x, float y) /* wrapper fmodf */
|
||||
float
|
||||
fmodf(float x, float y) /* wrapper fmodf */
|
||||
{
|
||||
#ifdef _IEEE_LIBM
|
||||
return __ieee754_fmodf(x,y);
|
||||
|
@ -25,7 +25,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
|
||||
extern int signgam;
|
||||
|
||||
double gamma(double x)
|
||||
double
|
||||
gamma(double x)
|
||||
{
|
||||
#ifdef _IEEE_LIBM
|
||||
return __ieee754_gamma_r(x,&signgam);
|
||||
|
@ -22,7 +22,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
#include "math_private.h"
|
||||
|
||||
|
||||
double gamma_r(double x, int *signgamp) /* wrapper lgamma_r */
|
||||
double
|
||||
gamma_r(double x, int *signgamp) /* wrapper lgamma_r */
|
||||
{
|
||||
#ifdef _IEEE_LIBM
|
||||
return __ieee754_gamma_r(x,signgamp);
|
||||
|
@ -22,7 +22,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
|
||||
extern int signgam;
|
||||
|
||||
float gammaf(float x)
|
||||
float
|
||||
gammaf(float x)
|
||||
{
|
||||
#ifdef _IEEE_LIBM
|
||||
return __ieee754_gammaf_r(x,&signgam);
|
||||
|
@ -25,7 +25,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
#include "math_private.h"
|
||||
|
||||
|
||||
float gammaf_r(float x, int *signgamp) /* wrapper lgammaf_r */
|
||||
float
|
||||
gammaf_r(float x, int *signgamp) /* wrapper lgammaf_r */
|
||||
{
|
||||
#ifdef _IEEE_LIBM
|
||||
return __ieee754_gammaf_r(x,signgamp);
|
||||
|
@ -22,7 +22,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
#include "math_private.h"
|
||||
|
||||
|
||||
double hypot(double x, double y)/* wrapper hypot */
|
||||
double
|
||||
hypot(double x, double y)/* wrapper hypot */
|
||||
{
|
||||
#ifdef _IEEE_LIBM
|
||||
return __ieee754_hypot(x,y);
|
||||
|
@ -25,7 +25,8 @@ static char rcsid[] = "$FreeBSD$";
|
||||
#include "math_private.h"
|
||||
|
||||
|
||||
float hypotf(float x, float y) /* wrapper hypotf */
|
||||
float
|
||||
hypotf(float x, float y) /* wrapper hypotf */
|
||||
{
|
||||
#ifdef _IEEE_LIBM
|
||||
return __ieee754_hypotf(x,y);
|
||||
|
Some files were not shown because too many files have changed in this diff Show More
Loading…
Reference in New Issue
Block a user