| /bionic/libm/upstream-freebsd/lib/msun/src/ |
| D | e_lgammaf_r.c | 94 float y,z; in sin_pif() local
101 y = -x; /* x is assume negative */ in sin_pif()
107 z = floorf(y); in sin_pif()
108 if(z!=y) { /* inexact anyway */ in sin_pif()
109 y *= (float)0.5; in sin_pif()
110 y = (float)2.0*(y - floorf(y)); /* y = |x| mod 2.0 */ in sin_pif()
111 n = (int) (y*(float)4.0); in sin_pif()
114 y = zero; n = 0; /* y must be even */ in sin_pif()
116 if(ix<0x4b000000) z = y+two23; /* exact */ in sin_pif()
119 y = n; in sin_pif()
[all …]
|
| D | e_lgamma_r.c | 161 double y,z; in sin_pi() local
168 y = -x; /* x is assume negative */ in sin_pi()
174 z = floor(y); in sin_pi()
175 if(z!=y) { /* inexact anyway */ in sin_pi()
176 y *= 0.5; in sin_pi()
177 y = 2.0*(y - floor(y)); /* y = |x| mod 2.0 */ in sin_pi()
178 n = (int) (y*4.0); in sin_pi()
181 y = zero; n = 0; /* y must be even */ in sin_pi()
183 if(ix<0x43300000) z = y+two52; /* exact */ in sin_pi()
186 y = n; in sin_pi()
[all …]
|
| D | s_ccoshf.c | 44 float x, y, h; in ccoshf() local
48 y = cimagf(z); in ccoshf()
51 GET_FLOAT_WORD(hy, y); in ccoshf()
58 return (cpackf(coshf(x), x * y)); in ccoshf()
60 return (cpackf(coshf(x) * cosf(y), sinhf(x) * sinf(y))); in ccoshf()
66 return (cpackf(h * cosf(y), copysignf(h, x) * sinf(y))); in ccoshf()
69 z = __ldexp_cexpf(cpackf(fabsf(x), y), -1); in ccoshf()
74 return (cpackf(h * h * cosf(y), h * sinf(y))); in ccoshf()
79 return (cpackf(y - y, copysignf(0, x * (y - y)))); in ccoshf()
83 return (cpackf(x * x, copysignf(0, x) * y)); in ccoshf()
[all …]
|
| D | s_csinh.c | 50 double x, y, h; in csinh() local
54 y = cimag(z); in csinh()
57 EXTRACT_WORDS(hy, ly, y); in csinh()
65 return (cpack(sinh(x), y)); in csinh()
67 return (cpack(sinh(x) * cos(y), cosh(x) * sin(y))); in csinh()
73 return (cpack(copysign(h, x) * cos(y), h * sin(y))); in csinh()
76 z = __ldexp_cexp(cpack(fabs(x), y), -1); in csinh()
81 return (cpack(h * cos(y), h * h * sin(y))); in csinh()
95 return (cpack(copysign(0, x * (y - y)), y - y)); in csinh()
104 return (cpack(x, y)); in csinh()
[all …]
|
| D | s_ccosh.c | 50 double x, y, h; in ccosh() local
54 y = cimag(z); in ccosh()
57 EXTRACT_WORDS(hy, ly, y); in ccosh()
65 return (cpack(cosh(x), x * y)); in ccosh()
67 return (cpack(cosh(x) * cos(y), sinh(x) * sin(y))); in ccosh()
73 return (cpack(h * cos(y), copysign(h, x) * sin(y))); in ccosh()
76 z = __ldexp_cexp(cpack(fabs(x), y), -1); in ccosh()
81 return (cpack(h * h * cos(y), h * sin(y))); in ccosh()
95 return (cpack(y - y, copysign(0, x * (y - y)))); in ccosh()
105 return (cpack(x * x, copysign(0, x) * y)); in ccosh()
[all …]
|
| D | s_csinhf.c | 44 float x, y, h; in csinhf() local
48 y = cimagf(z); in csinhf()
51 GET_FLOAT_WORD(hy, y); in csinhf()
58 return (cpackf(sinhf(x), y)); in csinhf()
60 return (cpackf(sinhf(x) * cosf(y), coshf(x) * sinf(y))); in csinhf()
66 return (cpackf(copysignf(h, x) * cosf(y), h * sinf(y))); in csinhf()
69 z = __ldexp_cexpf(cpackf(fabsf(x), y), -1); in csinhf()
74 return (cpackf(h * cosf(y), h * h * sinf(y))); in csinhf()
79 return (cpackf(copysignf(0, x * (y - y)), y - y)); in csinhf()
83 return (cpackf(x, y)); in csinhf()
[all …]
|
| D | e_rem_pio2.c | 55 __ieee754_rem_pio2(double x, double *y) in __ieee754_rem_pio2() argument
66 {y[0] = x; y[1] = 0; return 0;} in __ieee754_rem_pio2()
74 y[0] = z - pio2_1t; in __ieee754_rem_pio2()
75 y[1] = (z-y[0])-pio2_1t; in __ieee754_rem_pio2()
79 y[0] = z + pio2_1t; in __ieee754_rem_pio2()
80 y[1] = (z-y[0])+pio2_1t; in __ieee754_rem_pio2()
86 y[0] = z - 2*pio2_1t; in __ieee754_rem_pio2()
87 y[1] = (z-y[0])-2*pio2_1t; in __ieee754_rem_pio2()
91 y[0] = z + 2*pio2_1t; in __ieee754_rem_pio2()
92 y[1] = (z-y[0])+2*pio2_1t; in __ieee754_rem_pio2()
[all …]
|
| D | s_cexpf.c | 42 float x, y, exp_x; in cexpf() local
46 y = cimagf(z); in cexpf()
48 GET_FLOAT_WORD(hy, y); in cexpf()
53 return (cpackf(expf(x), y)); in cexpf()
57 return (cpackf(cosf(y), sinf(y))); in cexpf()
62 return (cpackf(y - y, y - y)); in cexpf()
68 return (cpackf(x, y - y)); in cexpf()
87 return (cpackf(exp_x * cosf(y), exp_x * sinf(y))); in cexpf()
|
| D | s_cexp.c | 42 double x, y, exp_x; in cexp() local
46 y = cimag(z); in cexp()
48 EXTRACT_WORDS(hy, ly, y); in cexp()
53 return (cpack(exp(x), y)); in cexp()
57 return (cpack(cos(y), sin(y))); in cexp()
62 return (cpack(y - y, y - y)); in cexp()
68 return (cpack(x, y - y)); in cexp()
87 return (cpack(exp_x * cos(y), exp_x * sin(y))); in cexp()
|
| D | s_ctanhf.c | 42 float x, y; in ctanhf() local
47 y = cimagf(z); in ctanhf()
54 return (cpackf(x, (y == 0 ? y : x * y))); in ctanhf()
57 copysignf(0, isinf(y) ? y : sinf(y) * cosf(y)))); in ctanhf()
60 if (!isfinite(y)) in ctanhf()
61 return (cpackf(y - y, y - y)); in ctanhf()
66 4 * sinf(y) * cosf(y) * exp_mx * exp_mx)); in ctanhf()
69 t = tanf(y); in ctanhf()
|
| D | s_ctanh.c | 77 double x, y; in ctanh() local
82 y = cimag(z); in ctanh()
105 return (cpack(x, (y == 0 ? y : x * y))); in ctanh()
107 return (cpack(x, copysign(0, isinf(y) ? y : sin(y) * cos(y)))); in ctanh()
114 if (!isfinite(y)) in ctanh()
115 return (cpack(y - y, y - y)); in ctanh()
125 4 * sin(y) * cos(y) * exp_mx * exp_mx)); in ctanh()
129 t = tan(y); in ctanh()
|
| D | s_cos.c | 57 double y[2],z=0.0; in cos() local
76 n = __ieee754_rem_pio2(x,y); in cos()
78 case 0: return __kernel_cos(y[0],y[1]); in cos()
79 case 1: return -__kernel_sin(y[0],y[1],1); in cos()
80 case 2: return -__kernel_cos(y[0],y[1]); in cos()
82 return __kernel_sin(y[0],y[1],1); in cos()
|
| D | s_sin.c | 57 double y[2],z=0.0; in sin() local
76 n = __ieee754_rem_pio2(x,y); in sin()
78 case 0: return __kernel_sin(y[0],y[1],1); in sin()
79 case 1: return __kernel_cos(y[0],y[1]); in sin()
80 case 2: return -__kernel_sin(y[0],y[1],1); in sin()
82 return -__kernel_cos(y[0],y[1]); in sin()
|
| D | s_expm1f.c | 44 float y,hi,lo,c,t,e,hxs,hfx,r1,twopk; in expm1f() local
106 y = one-(e-x); in expm1f()
107 if (k == 128) y = y*2.0F*0x1p127F; in expm1f()
108 else y = y*twopk; in expm1f()
109 return y-one; in expm1f()
114 y = t-(e-x); in expm1f()
115 y = y*twopk; in expm1f()
118 y = x-(e+t); in expm1f()
119 y += one; in expm1f()
120 y = y*twopk; in expm1f()
[all …]
|
| D | s_fdim.c | 34 fn(type x, type y) \
39 if (isnan(y)) \
40 return (y); \
41 return (x > y ? x - y : 0.0); \
|
| D | s_remquof.c | 30 remquof(float x, float y, int *quo) in remquof() argument
36 GET_FLOAT_WORD(hy,y); in remquof()
44 return (x*y)/(x*y); in remquof()
107 y = fabsf(y); in remquof()
108 if (y < 0x1p-125f) { in remquof()
109 if (x+x>y || (x+x==y && (q & 1))) { in remquof()
111 x-=y; in remquof()
113 } else if (x>0.5f*y || (x==0.5f*y && (q & 1))) { in remquof()
115 x-=y; in remquof()
|
| D | s_expm1.c | 135 double y,hi,lo,c,t,e,hxs,hfx,r1,twopk; in expm1() local
200 y = one-(e-x); in expm1()
201 if (k == 1024) y = y*2.0*0x1p1023; in expm1()
202 else y = y*twopk; in expm1()
203 return y-one; in expm1()
208 y = t-(e-x); in expm1()
209 y = y*twopk; in expm1()
212 y = x-(e+t); in expm1()
213 y += one; in expm1()
214 y = y*twopk; in expm1()
[all …]
|
| D | s_nextafterf.c | 23 nextafterf(float x, float y) in nextafterf() argument
29 GET_FLOAT_WORD(hy,y); in nextafterf()
35 return x+y; in nextafterf()
36 if(x==y) return y; /* x=y, return y */ in nextafterf()
60 SET_FLOAT_WORD(y,hx); in nextafterf()
61 return y; in nextafterf()
|
| D | s_fmaxl.c | 35 fmaxl(long double x, long double y) in fmaxl() argument
41 u[1].e = y; in fmaxl()
46 return (y); in fmaxl()
52 return (u[0].bits.sign ? y : x); in fmaxl()
54 return (x > y ? x : y); in fmaxl()
|
| D | s_fminl.c | 35 fminl(long double x, long double y) in fminl() argument
41 u[1].e = y; in fminl()
46 return (y); in fminl()
52 return (u[1].bits.sign ? y : x); in fminl()
54 return (x < y ? x : y); in fminl()
|
| D | s_fminf.c | 35 fminf(float x, float y) in fminf() argument
40 u[1].f = y; in fminf()
44 return (y); in fminf()
52 return (x < y ? x : y); in fminf()
|
| /bionic/libc/include/sys/ |
| D | param.h | 39 #define howmany(x, y) (((x)+((y)-1))/(y)) argument
41 #define roundup(x, y) ((((x)+((y)-1))/(y))*(y)) argument
|
| /bionic/libm/upstream-freebsd/lib/msun/bsdsrc/ |
| D | b_tgamma.c | 188 double y, ym1, t; local
190 y = x - one;
191 ym1 = y - one;
192 if (y <= 1.0 + (LEFT + x0)) {
193 yy = ratfun_gam(y - x0, 0);
196 r.a = y;
199 y = ym1;
200 yy.b = r.b = y - yy.a;
202 for (ym1 = y-one; ym1 > LEFT + x0; y = ym1--, yy.a--) {
204 r.b = r.a*yy.b + y*r.b;
[all …]
|
| /bionic/libc/upstream-openbsd/lib/libc/gdtoa/ |
| D | sum.c | 42 ULong carry, *xc, *xa, *xb, *xe, y; local
61 y = (*xa & 0xffff) + (*xb & 0xffff) + carry;
62 carry = (y & 0x10000) >> 16;
65 Storeinc(xc, z, y);
70 y = (*xa & 0xffff) + carry;
71 carry = (y & 0x10000) >> 16;
74 Storeinc(xc, z, y);
78 y = *xa++ + *xb++ + carry;
79 carry = (y & 0x10000) >> 16;
80 *xc++ = y & 0xffff;
[all …]
|
| D | dmisc.c | 119 ULLong borrow, carry, y, ys; local
121 ULong borrow, carry, y, ys; local
150 y = *bx - (ys & 0xffffffffUL) - borrow;
151 borrow = y >> 32 & 1UL;
152 *bx++ = y & 0xffffffffUL;
159 y = (*bx & 0xffff) - (ys & 0xffff) - borrow;
160 borrow = (y & 0x10000) >> 16;
163 Storeinc(bx, z, y);
167 y = *bx - (ys & 0xffff) - borrow;
168 borrow = (y & 0x10000) >> 16;
[all …]
|