| /bionic/libm/upstream-freebsd/lib/msun/src/ |
| D | s_erff.c | 27 one = 1, variable
89 return (float)(1-i)+one/x; /* erff(+-inf)=+-1 */ in erff()
100 s = one+z*(qq1+z*(qq2+z*qq3)); in erff()
105 s = fabsf(x)-one; in erff()
107 Q = one+s*(qa1+s*(qa2+s*qa3)); in erff()
111 if(hx>=0) return one-tiny; else return tiny-one; in erff()
114 s = one/(x*x); in erff()
117 S=one+s*(sa1+s*(sa2+s*sa3)); in erff()
120 S=one+s*(sb1+s*(sb2+s*sb3)); in erff()
124 if(hx>=0) return one-r/x; else return r/x-one; in erff()
[all …]
|
| D | s_erf.c | 119 one = 1, variable
205 return (double)(1-i)+one/x; /* erf(+-inf)=+-1 */ in erf()
216 s = one+z*(qq1+z*(qq2+z*(qq3+z*(qq4+z*qq5)))); in erf()
221 s = fabs(x)-one; in erf()
223 Q = one+s*(qa1+s*(qa2+s*(qa3+s*(qa4+s*(qa5+s*qa6))))); in erf()
227 if(hx>=0) return one-tiny; else return tiny-one; in erf()
230 s = one/(x*x); in erf()
233 S=one+s*(sa1+s*(sa2+s*(sa3+s*(sa4+s*(sa5+s*(sa6+s*(sa7+ in erf()
237 S=one+s*(sb1+s*(sb2+s*(sb3+s*(sb4+s*(sb5+s*(sb6+s*sb7)))))); in erf()
242 if(hx>=0) return one-r/x; else return r/x-one; in erf()
[all …]
|
| D | s_tanhf.c | 23 static const float one=1.0, two=2.0, huge = 1.0e30; variable
36 if (jx>=0) return one/x+one; /* tanh(+-inf)=+-1 */ in tanhf()
37 else return one/x-one; /* tanh(NaN) = NaN */ in tanhf()
43 if(huge+x>one) return x; /* tanh(tiny) = tiny with inexact */ in tanhf()
47 z = one - two/(t+two); in tanhf()
54 z = one - tiny; /* raise inexact flag */ in tanhf()
|
| D | s_tanh.c | 46 static const double one = 1.0, two = 2.0, huge = 1.0e300; variable
59 if (jx>=0) return one/x+one; /* tanh(+-inf)=+-1 */ in tanh()
60 else return one/x-one; /* tanh(NaN) = NaN */ in tanh()
66 if(huge+x>one) return x; /* tanh(tiny) = tiny with inexact */ in tanh()
70 z = one - two/(t+two); in tanh()
77 z = one - tiny; /* raise inexact flag */ in tanh()
|
| D | e_lgammaf_r.c | 28 one = 1, variable
113 case 4: y = __kernel_sindf(pi*(one-y)); break; in sin_pif()
140 return one/vzero; in __ieee754_lgammaf_r()
148 return one/vzero; in __ieee754_lgammaf_r()
150 if(t==zero) return one/vzero; /* -integer */ in __ieee754_lgammaf_r()
162 if(ix>=0x3f3b4a20) {y = one-x; i= 0;} in __ieee754_lgammaf_r()
163 else if(ix>=0x3e6d3308) {y= x-(tc-one); i=1;} in __ieee754_lgammaf_r()
169 else {y=x-one;i=2;} in __ieee754_lgammaf_r()
183 p2 = one+y*(v1+y*(v2+y*v3)); in __ieee754_lgammaf_r()
192 q = one+y*(r1+y*(r2+y*r3)); in __ieee754_lgammaf_r()
[all …]
|
| D | e_j0f.c | 32 one = 1.0, variable
55 if(ix>=0x7f800000) return one/(x*x); in __ieee754_j0f()
78 if(huge+x>one) { /* raise inexact if x != 0 */ in __ieee754_j0f()
79 if(ix<0x39800000) return one; /* |x|<2**-12 */ in __ieee754_j0f()
80 else return one - x*x/4; in __ieee754_j0f()
85 s = one+z*(S01+z*(S02+z*(S03+z*S04))); in __ieee754_j0f()
87 return one + z*((r/s)-qrtr); in __ieee754_j0f()
90 return((one+u)*(one-u)+z*(r/s)); in __ieee754_j0f()
116 if(ix==0) return -one/vzero; in __ieee754_y0f()
154 v = one+z*(v01+z*(v02+z*(v03+z*v04))); in __ieee754_y0f()
[all …]
|
| D | e_j0.c | 70 one = 1.0, variable
93 if(ix>=0x7ff00000) return one/(x*x); in __ieee754_j0()
116 if(huge+x>one) { /* raise inexact if x != 0 */ in __ieee754_j0()
117 if(ix<0x3e400000) return one; /* |x|<2**-27 */ in __ieee754_j0()
118 else return one - x*x/4; in __ieee754_j0()
123 s = one+z*(S01+z*(S02+z*(S03+z*S04))); in __ieee754_j0()
125 return one + z*((r/s)-qrtr); in __ieee754_j0()
128 return((one+u)*(one-u)+z*(r/s)); in __ieee754_j0()
160 if((ix|lx)==0) return -one/vzero; in __ieee754_y0()
199 v = one+z*(v01+z*(v02+z*(v03+z*v04))); in __ieee754_y0()
[all …]
|
| D | e_lgamma_r.c | 93 one = 1.00000000000000000000e+00, /* 0x3FF00000, 0x00000000 */ variable
192 case 4: y = __kernel_sin(pi*(one-y),zero,0); break; in sin_pi()
219 return one/vzero; in __ieee754_lgamma_r()
227 return one/vzero; in __ieee754_lgamma_r()
229 if(t==zero) return one/vzero; /* -integer */ in __ieee754_lgamma_r()
241 if(ix>=0x3FE76944) {y = one-x; i= 0;} in __ieee754_lgamma_r()
242 else if(ix>=0x3FCDA661) {y= x-(tc-one); i=1;} in __ieee754_lgamma_r()
248 else {y=x-one;i=2;} in __ieee754_lgamma_r()
267 p2 = one+y*(v1+y*(v2+y*(v3+y*(v4+y*v5)))); in __ieee754_lgamma_r()
276 q = one+y*(r1+y*(r2+y*(r3+y*(r4+y*(r5+y*r6))))); in __ieee754_lgamma_r()
[all …]
|
| D | s_asinhf.c | 23 one = 1.0000000000e+00, /* 0x3F800000 */ variable
36 if(huge+x>one) return x; /* return x inexact except 0 */ in asinhf()
42 w = __ieee754_logf((float)2.0*t+one/(__ieee754_sqrtf(x*x+one)+t)); in asinhf()
45 w =log1pf(fabsf(x)+t/(one+__ieee754_sqrtf(one+t))); in asinhf()
|
| D | s_asinh.c | 33 one = 1.00000000000000000000e+00, /* 0x3FF00000, 0x00000000 */ variable
46 if(huge+x>one) return x; /* return x inexact except 0 */ in asinh()
52 w = __ieee754_log(2.0*t+one/(__ieee754_sqrt(x*x+one)+t)); in asinh()
55 w =log1p(fabs(x)+t/(one+__ieee754_sqrt(one+t))); in asinh()
|
| D | s_expm1f.c | 25 one = 1.0, variable
63 return tiny-one; /* return -1 */ in expm1f()
92 r1 = one+hxs*(Q1+hxs*Q2); in expm1f()
103 else return one+(float)2.0*(x-e); in expm1f()
106 y = one-(e-x); in expm1f()
109 return y-one; in expm1f()
111 t = one; in expm1f()
119 y += one; in expm1f()
|
| D | e_j1.c | 70 one = 1.0, variable
94 if(ix>=0x7ff00000) return one/x; in __ieee754_j1()
118 if(huge+x>one) return 0.5*x;/* inexact if x!=0 necessary */ in __ieee754_j1()
122 s = one+z*(s01+z*(s02+z*(s03+z*(s04+z*s05)))); in __ieee754_j1()
157 if((ix|lx)==0) return -one/vzero; in __ieee754_y1()
192 v = one+z*(V0[0]+z*(V0[1]+z*(V0[2]+z*(V0[3]+z*V0[4])))); in __ieee754_y1()
193 return(x*(u/v) + tpi*(__ieee754_j1(x)*__ieee754_log(x)-one/x)); in __ieee754_y1()
282 z = one/(x*x); in pone()
284 s = one+z*(q[0]+z*(q[1]+z*(q[2]+z*(q[3]+z*q[4])))); in pone()
285 return one+ r/s; in pone()
[all …]
|
| D | e_j1f.c | 32 one = 1.0, variable
56 if(ix>=0x7f800000) return one/x; in __ieee754_j1f()
80 if(huge+x>one) return (float)0.5*x;/* inexact if x!=0 necessary */ in __ieee754_j1f()
84 s = one+z*(s01+z*(s02+z*(s03+z*(s04+z*s05)))); in __ieee754_j1f()
113 if(ix==0) return -one/vzero; in __ieee754_y1f()
147 v = one+z*(V0[0]+z*(V0[1]+z*(V0[2]+z*(V0[3]+z*V0[4])))); in __ieee754_y1f()
148 return(x*(u/v) + tpi*(__ieee754_j1f(x)*__ieee754_logf(x)-one/x)); in __ieee754_y1f()
237 z = one/(x*x); in ponef()
239 s = one+z*(q[0]+z*(q[1]+z*(q[2]+z*(q[3]+z*q[4])))); in ponef()
240 return one+ r/s; in ponef()
[all …]
|
| D | s_asinhl.c | 54 one = 1.00000000000000000000e+00, /* 0x3FF00000, 0x00000000 */ variable
79 if (huge + x > one) RETURNI(x); /* return x inexact except 0 */ in asinhl()
85 w = logl(2.0*t+one/(sqrtl(x*x+one)+t)); in asinhl()
88 w =log1pl(fabsl(x)+t/(one+sqrtl(one+t))); in asinhl()
|
| D | s_expm1.c | 117 one = 1.0, variable
157 return tiny-one; /* return -1 */ in expm1()
186 r1 = one+hxs*(Q1+hxs*(Q2+hxs*(Q3+hxs*(Q4+hxs*Q5)))); in expm1()
197 else return one+2.0*(x-e); in expm1()
200 y = one-(e-x); in expm1()
203 return y-one; in expm1()
205 t = one; in expm1()
213 y += one; in expm1()
|
| D | s_atanf.c | 45 one = 1.0, variable
63 if(huge+x>one) return x; /* raise inexact */ in atanf()
70 id = 0; x = ((float)2.0*x-one)/((float)2.0+x); in atanf()
72 id = 1; x = (x-one)/(x+one); in atanf()
76 id = 2; x = (x-(float)1.5)/(one+(float)1.5*x); in atanf()
|
| D | s_atanl.c | 29 one = 1.0, variable
56 if(huge+x>one) return x; /* raise inexact */ in atanl()
63 id = 0; x = (2.0*x-one)/(2.0+x); in atanl()
65 id = 1; x = (x-one)/(x+one); in atanl()
69 id = 2; x = (x-1.5)/(one+1.5*x); in atanl()
|
| D | s_atan.c | 70 one = 1.0, variable
91 if(huge+x>one) return x; /* raise inexact */ in atan()
98 id = 0; x = (2.0*x-one)/(2.0+x); in atan()
100 id = 1; x = (x-one)/(x+one); in atan()
104 id = 2; x = (x-1.5)/(one+1.5*x); in atan()
|
| D | e_acosf.c | 23 one = 1.0000000000e+00, /* 0x3F800000 */ variable
52 q = one+z*qS1; in __ieee754_acosf()
56 z = (one+x)*(float)0.5; in __ieee754_acosf()
58 q = one+z*qS1; in __ieee754_acosf()
65 z = (one-x)*(float)0.5; in __ieee754_acosf()
72 q = one+z*qS1; in __ieee754_acosf()
|
| D | s_tanhl.c | 42 static const double one = 1.0; variable
126 if (jx>=0) return one/x+one; /* tanh(+-inf)=+-1 */ in tanhl()
127 else return one/x-one; /* tanh(NaN) = NaN */ in tanhl()
166 z = one - one/(lo+0.5+hi); in tanhl()
169 z = one - tiny; /* raise inexact flag */ in tanhl()
|
| D | e_asinf.c | 23 one = 1.0000000000e+00, /* 0x3F800000 */ variable
48 if(huge+x>one) return x;/* return x with inexact if x!=0*/ in __ieee754_asinf()
52 q = one+t*qS1; in __ieee754_asinf()
57 w = one-fabsf(x); in __ieee754_asinf()
60 q = one+t*qS1; in __ieee754_asinf()
|
| D | e_sqrtf.c | 23 static const float one = 1.0, tiny=1.0e-30; variable
76 z = one-tiny; /* trigger inexact flag */ in __ieee754_sqrtf()
77 if (z>=one) { in __ieee754_sqrtf()
78 z = one+tiny; in __ieee754_sqrtf()
79 if (z>one) in __ieee754_sqrtf()
|
| D | e_acos.c | 47 one= 1.00000000000000000000e+00, /* 0x3FF00000, 0x00000000 */ variable
84 q = one+z*(qS1+z*(qS2+z*(qS3+z*qS4))); in __ieee754_acos()
88 z = (one+x)*0.5; in __ieee754_acos()
90 q = one+z*(qS1+z*(qS2+z*(qS3+z*qS4))); in __ieee754_acos()
96 z = (one-x)*0.5; in __ieee754_acos()
102 q = one+z*(qS1+z*(qS2+z*(qS3+z*qS4))); in __ieee754_acos()
|
| /bionic/libm/upstream-freebsd/lib/msun/ld128/ |
| D | s_erfl.c | 32 one = 1, variable
205 return (1-i)+one/x; /* erfl(+-inf)=+-1 */ in erfl()
218 s = one+z*(qq1+z*(qq2+z*(qq3+z*(qq4+z*(qq5+z*(qq6+z*(qq7+ in erfl()
224 s = ax-one; in erfl()
227 Q = one+s*(qa1+s*(qa2+s*(qa3+s*(qa4+s*(qa5+s*(qa6+s*(qa7+ in erfl()
232 if(x>=0) return (one-tiny); else return (tiny-one); in erfl()
234 s = one/(ax*ax); in erfl()
239 S=one+s*(sa1+s*(sa2+s*(sa3+s*(sa4+s*(sa5+s*(sa6+s*(sa7+ in erfl()
246 S=one+s*(sb1+s*(sb2+s*(sb3+s*(sb4+s*(sb5+s*(sb6+s*(sb7+ in erfl()
252 if(x>=0) return (one-r/ax); else return (r/ax-one); in erfl()
[all …]
|
| /bionic/libm/upstream-freebsd/lib/msun/bsdsrc/ |
| D | b_tgamma.c | 128 static const double zero = 0., one = 1.0, tiny = 1e-300; variable
147 u.a = one - tiny; /* raise inexact */
148 return (one/x);
164 z = one/(x*x);
169 u.a -= one;
192 y = x - one;
193 ym1 = y - one;
200 yy.a = r.a - one;
204 for (ym1 = y-one; ym1 > LEFT + x0; y = ym1--, yy.a--) {
232 t = (one-x0); t += x;
[all …]
|