Home
last modified time | relevance | path

Searched refs:y (Results 1 – 25 of 140) sorted by relevance

123456

/bionic/libm/upstream-freebsd/lib/msun/ld128/
De_lgammal_r.c190 long double y,z; in sin_pil() local
194 y = -x; in sin_pil()
196 vz = y+0x1.p112; in sin_pil()
198 if (z == y) in sin_pil()
201 vz = y+0x1.p110; in sin_pil()
204 if (z > y) { in sin_pil()
209 y = y - z + n * 0.25; in sin_pil()
212 case 0: y = __kernel_sinl(pi*y,zero,0); break; in sin_pil()
214 case 2: y = __kernel_cosl(pi*(0.5-y),zero); break; in sin_pil()
216 case 4: y = __kernel_sinl(pi*(one-y),zero,0); break; in sin_pil()
[all …]
/bionic/libm/upstream-freebsd/lib/msun/src/
De_lgammaf_r.c88 float y,z; in sin_pif() local
91 y = -x; in sin_pif()
93 vz = y+0x1p23F; /* depend on 0 <= y < 0x1p23 */ in sin_pif()
95 if (z == y) in sin_pif()
98 vz = y+0x1p21F; in sin_pif()
101 if (z > y) { in sin_pif()
106 y = y - z + n * 0.25F; /* y mod 2 */ in sin_pif()
109 case 0: y = __kernel_sindf(pi*y); break; in sin_pif()
111 case 2: y = __kernel_cosdf(pi*((float)0.5-y)); break; in sin_pif()
113 case 4: y = __kernel_sindf(pi*(one-y)); break; in sin_pif()
[all …]
De_lgamma_r.c167 double y,z; in sin_pi() local
170 y = -x; in sin_pi()
172 vz = y+0x1p52; /* depend on 0 <= y < 0x1p52 */ in sin_pi()
174 if (z == y) in sin_pi()
177 vz = y+0x1p50; in sin_pi()
180 if (z > y) { in sin_pi()
185 y = y - z + n * 0.25; /* y mod 2 */ in sin_pi()
188 case 0: y = __kernel_sin(pi*y,zero,0); break; in sin_pi()
190 case 2: y = __kernel_cos(pi*(0.5-y),zero); break; in sin_pi()
192 case 4: y = __kernel_sin(pi*(one-y),zero,0); break; in sin_pi()
[all …]
Ds_ccoshf.c44 float x, y, h; in ccoshf() local
48 y = cimagf(z); in ccoshf()
51 GET_FLOAT_WORD(hy, y); in ccoshf()
58 return (CMPLXF(coshf(x), x * y)); in ccoshf()
60 return (CMPLXF(coshf(x) * cosf(y), sinhf(x) * sinf(y))); in ccoshf()
66 return (CMPLXF(h * cosf(y), copysignf(h, x) * sinf(y))); in ccoshf()
69 z = __ldexp_cexpf(CMPLXF(fabsf(x), y), -1); in ccoshf()
74 return (CMPLXF(h * h * cosf(y), h * sinf(y))); in ccoshf()
79 return (CMPLXF(y - y, x * copysignf(0, y))); in ccoshf()
82 return (CMPLXF(x * x, copysignf(0, x) * y)); in ccoshf()
[all …]
Ds_csinhf.c44 float x, y, h; in csinhf() local
48 y = cimagf(z); in csinhf()
51 GET_FLOAT_WORD(hy, y); in csinhf()
58 return (CMPLXF(sinhf(x), y)); in csinhf()
60 return (CMPLXF(sinhf(x) * cosf(y), coshf(x) * sinf(y))); in csinhf()
66 return (CMPLXF(copysignf(h, x) * cosf(y), h * sinf(y))); in csinhf()
69 z = __ldexp_cexpf(CMPLXF(fabsf(x), y), -1); in csinhf()
74 return (CMPLXF(h * cosf(y), h * h * sinf(y))); in csinhf()
79 return (CMPLXF(x, y - y)); in csinhf()
82 return (CMPLXF(x + x, y)); in csinhf()
[all …]
Ds_ccosh.c52 double x, y, h; in ccosh() local
56 y = cimag(z); in ccosh()
59 EXTRACT_WORDS(hy, ly, y); in ccosh()
67 return (CMPLX(cosh(x), x * y)); in ccosh()
69 return (CMPLX(cosh(x) * cos(y), sinh(x) * sin(y))); in ccosh()
75 return (CMPLX(h * cos(y), copysign(h, x) * sin(y))); in ccosh()
78 z = __ldexp_cexp(CMPLX(fabs(x), y), -1); in ccosh()
83 return (CMPLX(h * h * cos(y), h * sin(y))); in ccosh()
98 return (CMPLX(y - y, x * copysign(0, y))); in ccosh()
108 return (CMPLX(x * x, copysign(0, x) * y)); in ccosh()
[all …]
Ds_csinh.c52 double x, y, h; in csinh() local
56 y = cimag(z); in csinh()
59 EXTRACT_WORDS(hy, ly, y); in csinh()
67 return (CMPLX(sinh(x), y)); in csinh()
69 return (CMPLX(sinh(x) * cos(y), cosh(x) * sin(y))); in csinh()
75 return (CMPLX(copysign(h, x) * cos(y), h * sin(y))); in csinh()
78 z = __ldexp_cexp(CMPLX(fabs(x), y), -1); in csinh()
83 return (CMPLX(h * cos(y), h * h * sin(y))); in csinh()
97 return (CMPLX(x, y - y)); in csinh()
105 return (CMPLX(x + x, y)); in csinh()
[all …]
De_rem_pio2.c55 __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 …]
Dcatrig.c153 do_hard_work(double x, double y, double *rx, int *B_is_usable, double *B, in do_hard_work() argument
159 R = hypot(x, y + 1); /* |z+I| */ in do_hard_work()
160 S = hypot(x, y - 1); /* |z-I| */ in do_hard_work()
177 if (y == 1 && x < DBL_EPSILON * DBL_EPSILON / 128) { in do_hard_work()
183 } else if (x >= DBL_EPSILON * fabs(y - 1)) { in do_hard_work()
188 Am1 = f(x, 1 + y, R) + f(x, 1 - y, S); in do_hard_work()
190 } else if (y < 1) { in do_hard_work()
195 *rx = x / sqrt((1 - y) * (1 + y)); in do_hard_work()
200 *rx = log1p((y - 1) + sqrt((y - 1) * (y + 1))); in do_hard_work()
206 *new_y = y; in do_hard_work()
[all …]
Dcatrigf.c88 do_hard_work(float x, float y, float *rx, int *B_is_usable, float *B, in do_hard_work() argument
94 R = hypotf(x, y + 1); in do_hard_work()
95 S = hypotf(x, y - 1); in do_hard_work()
102 if (y == 1 && x < FLT_EPSILON * FLT_EPSILON / 128) { in do_hard_work()
104 } else if (x >= FLT_EPSILON * fabsf(y - 1)) { in do_hard_work()
105 Am1 = f(x, 1 + y, R) + f(x, 1 - y, S); in do_hard_work()
107 } else if (y < 1) { in do_hard_work()
108 *rx = x / sqrtf((1 - y) * (1 + y)); in do_hard_work()
110 *rx = log1pf((y - 1) + sqrtf((y - 1) * (y + 1))); in do_hard_work()
116 *new_y = y; in do_hard_work()
[all …]
Ds_cexpf.c42 float x, y, exp_x; in cexpf() local
46 y = cimagf(z); in cexpf()
48 GET_FLOAT_WORD(hy, y); in cexpf()
53 return (CMPLXF(expf(x), y)); in cexpf()
57 return (CMPLXF(cosf(y), sinf(y))); in cexpf()
62 return (CMPLXF(y - y, y - y)); in cexpf()
68 return (CMPLXF(x, y - y)); in cexpf()
87 return (CMPLXF(exp_x * cosf(y), exp_x * sinf(y))); in cexpf()
Ds_ctanhf.c42 float x, y; in ctanhf() local
47 y = cimagf(z); in ctanhf()
54 return (CMPLXF((x + 0) * (y + 0), in ctanhf()
55 y == 0 ? y : (x + 0) * (y + 0))); in ctanhf()
58 copysignf(0, isinf(y) ? y : sinf(y) * cosf(y)))); in ctanhf()
61 if (!isfinite(y)) in ctanhf()
62 return (CMPLXF(y - y, y - y)); in ctanhf()
67 4 * sinf(y) * cosf(y) * exp_mx * exp_mx)); in ctanhf()
70 t = tanf(y); in ctanhf()
Ds_cexp.c42 double x, y, exp_x; in cexp() local
46 y = cimag(z); in cexp()
48 EXTRACT_WORDS(hy, ly, y); in cexp()
53 return (CMPLX(exp(x), y)); in cexp()
57 return (CMPLX(cos(y), sin(y))); in cexp()
62 return (CMPLX(y - y, y - y)); in cexp()
68 return (CMPLX(x, y - y)); in cexp()
87 return (CMPLX(exp_x * cos(y), exp_x * sin(y))); in cexp()
Ds_ctanh.c77 double x, y; in ctanh() local
82 y = cimag(z); in ctanh()
105 return (CMPLX((x + 0) * (y + 0), in ctanh()
106 y == 0 ? y : (x + 0) * (y + 0))); in ctanh()
108 return (CMPLX(x, copysign(0, isinf(y) ? y : sin(y) * cos(y)))); in ctanh()
115 if (!isfinite(y)) in ctanh()
116 return (CMPLX(y - y, y - y)); in ctanh()
126 4 * sin(y) * cos(y) * exp_mx * exp_mx)); in ctanh()
130 t = tan(y); in ctanh()
Ds_cos.c57 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()
Ds_sin.c57 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()
Ds_expm1f.c44 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 …]
Ds_remquof.c30 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()
Ds_expm1.c135 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 …]
Ds_fdim.c34 fn(type x, type y) \
39 if (isnan(y)) \
40 return (y); \
41 return (x > y ? x - y : 0.0); \
Ds_nextafterf.c23 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()
/bionic/libc/include/sys/
Dparam.h43 #define howmany(x, y) (((x)+((y)-1))/(y)) argument
45 #define roundup(x, y) ((((x)+((y)-1))/(y))*(y)) argument
/bionic/libm/upstream-freebsd/lib/msun/bsdsrc/
Db_tgamma.c188 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/
Dsum.c42 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 …]
Ddmisc.c119 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 …]

123456