Home
last modified time | relevance | path

Searched refs:y (Results 1 – 25 of 145) 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 …]
Dcatrigl.c109 do_hard_work(long double x, long double y, long double *rx, int *B_is_usable, in do_hard_work() argument
115 R = hypotl(x, y + 1); in do_hard_work()
116 S = hypotl(x, y - 1); in do_hard_work()
123 if (y == 1 && x < LDBL_EPSILON * LDBL_EPSILON / 128) { in do_hard_work()
125 } else if (x >= LDBL_EPSILON * fabsl(y - 1)) { in do_hard_work()
126 Am1 = f(x, 1 + y, R) + f(x, 1 - y, S); in do_hard_work()
128 } else if (y < 1) { in do_hard_work()
129 *rx = x / sqrtl((1 - y) * (1 + y)); in do_hard_work()
131 *rx = log1pl((y - 1) + sqrtl((y - 1) * (y + 1))); in do_hard_work()
137 *new_y = y; in do_hard_work()
[all …]
Dcatrig.c155 do_hard_work(double x, double y, double *rx, int *B_is_usable, double *B, in do_hard_work() argument
161 R = hypot(x, y + 1); /* |z+I| */ in do_hard_work()
162 S = hypot(x, y - 1); /* |z-I| */ in do_hard_work()
179 if (y == 1 && x < DBL_EPSILON * DBL_EPSILON / 128) { in do_hard_work()
185 } else if (x >= DBL_EPSILON * fabs(y - 1)) { in do_hard_work()
190 Am1 = f(x, 1 + y, R) + f(x, 1 - y, S); in do_hard_work()
192 } else if (y < 1) { in do_hard_work()
197 *rx = x / sqrt((1 - y) * (1 + y)); in do_hard_work()
202 *rx = log1p((y - 1) + sqrt((y - 1) * (y + 1))); in do_hard_work()
208 *new_y = y; in do_hard_work()
[all …]
Dcatrigf.c90 do_hard_work(float x, float y, float *rx, int *B_is_usable, float *B, in do_hard_work() argument
96 R = hypotf(x, y + 1); in do_hard_work()
97 S = hypotf(x, y - 1); in do_hard_work()
104 if (y == 1 && x < FLT_EPSILON * FLT_EPSILON / 128) { in do_hard_work()
106 } else if (x >= FLT_EPSILON * fabsf(y - 1)) { in do_hard_work()
107 Am1 = f(x, 1 + y, R) + f(x, 1 - y, S); in do_hard_work()
109 } else if (y < 1) { in do_hard_work()
110 *rx = x / sqrtf((1 - y) * (1 + y)); in do_hard_work()
112 *rx = log1pf((y - 1) + sqrtf((y - 1) * (y + 1))); in do_hard_work()
118 *new_y = y; in do_hard_work()
[all …]
Ds_ccoshf.c46 float x, y, h; in ccoshf() local
50 y = cimagf(z); in ccoshf()
53 GET_FLOAT_WORD(hy, y); in ccoshf()
60 return (CMPLXF(coshf(x), x * y)); in ccoshf()
62 return (CMPLXF(coshf(x) * cosf(y), sinhf(x) * sinf(y))); in ccoshf()
68 return (CMPLXF(h * cosf(y), copysignf(h, x) * sinf(y))); in ccoshf()
71 z = __ldexp_cexpf(CMPLXF(fabsf(x), y), -1); in ccoshf()
76 return (CMPLXF(h * h * cosf(y), h * sinf(y))); in ccoshf()
81 return (CMPLXF(y - y, x * copysignf(0, y))); in ccoshf()
84 return (CMPLXF(x * x, copysignf(0, x) * y)); in ccoshf()
[all …]
Ds_ccosh.c54 double x, y, h; in ccosh() local
58 y = cimag(z); in ccosh()
61 EXTRACT_WORDS(hy, ly, y); in ccosh()
69 return (CMPLX(cosh(x), x * y)); in ccosh()
71 return (CMPLX(cosh(x) * cos(y), sinh(x) * sin(y))); in ccosh()
77 return (CMPLX(h * cos(y), copysign(h, x) * sin(y))); in ccosh()
80 z = __ldexp_cexp(CMPLX(fabs(x), y), -1); in ccosh()
85 return (CMPLX(h * h * cos(y), h * sin(y))); in ccosh()
100 return (CMPLX(y - y, x * copysign(0, y))); in ccosh()
110 return (CMPLX(x * x, copysign(0, x) * y)); in ccosh()
[all …]
Ds_csinhf.c46 float x, y, h; in csinhf() local
50 y = cimagf(z); in csinhf()
53 GET_FLOAT_WORD(hy, y); in csinhf()
60 return (CMPLXF(sinhf(x), y)); in csinhf()
62 return (CMPLXF(sinhf(x) * cosf(y), coshf(x) * sinf(y))); in csinhf()
68 return (CMPLXF(copysignf(h, x) * cosf(y), h * sinf(y))); in csinhf()
71 z = __ldexp_cexpf(CMPLXF(fabsf(x), y), -1); in csinhf()
76 return (CMPLXF(h * cosf(y), h * h * sinf(y))); in csinhf()
81 return (CMPLXF(x, y - y)); in csinhf()
84 return (CMPLXF(x + x, y)); in csinhf()
[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 …]
Ds_csinh.c54 double x, y, h; in csinh() local
58 y = cimag(z); in csinh()
61 EXTRACT_WORDS(hy, ly, y); in csinh()
69 return (CMPLX(sinh(x), y)); in csinh()
71 return (CMPLX(sinh(x) * cos(y), cosh(x) * sin(y))); in csinh()
77 return (CMPLX(copysign(h, x) * cos(y), h * sin(y))); in csinh()
80 z = __ldexp_cexp(CMPLX(fabs(x), y), -1); in csinh()
85 return (CMPLX(h * cos(y), h * h * sin(y))); in csinh()
99 return (CMPLX(x, y - y)); in csinh()
107 return (CMPLX(x + x, y)); in csinh()
[all …]
Ds_ctanhf.c44 float x, y; in ctanhf() local
49 y = cimagf(z); in ctanhf()
56 return (CMPLXF(nan_mix(x, y), in ctanhf()
57 y == 0 ? y : nan_mix(x, y))); in ctanhf()
60 copysignf(0, isinf(y) ? y : sinf(y) * cosf(y)))); in ctanhf()
63 if (!isfinite(y)) in ctanhf()
64 return (CMPLXF(y - y, y - y)); in ctanhf()
69 4 * sinf(y) * cosf(y) * exp_mx * exp_mx)); in ctanhf()
72 t = tanf(y); in ctanhf()
Ds_cexpf.c44 float x, y, exp_x; in cexpf() local
48 y = cimagf(z); in cexpf()
50 GET_FLOAT_WORD(hy, y); in cexpf()
55 return (CMPLXF(expf(x), y)); in cexpf()
59 return (CMPLXF(cosf(y), sinf(y))); in cexpf()
64 return (CMPLXF(y - y, y - y)); in cexpf()
70 return (CMPLXF(x, y - y)); in cexpf()
89 return (CMPLXF(exp_x * cosf(y), exp_x * sinf(y))); in cexpf()
Ds_ctanh.c79 double x, y; in ctanh() local
84 y = cimag(z); in ctanh()
107 return (CMPLX(nan_mix(x, y), in ctanh()
108 y == 0 ? y : nan_mix(x, y))); in ctanh()
110 return (CMPLX(x, copysign(0, isinf(y) ? y : sin(y) * cos(y)))); in ctanh()
117 if (!isfinite(y)) in ctanh()
118 return (CMPLX(y - y, y - y)); in ctanh()
128 4 * sin(y) * cos(y) * exp_mx * exp_mx)); in ctanh()
132 t = tan(y); in ctanh()
Ds_cexp.c44 double x, y, exp_x; in cexp() local
48 y = cimag(z); in cexp()
50 EXTRACT_WORDS(hy, ly, y); in cexp()
55 return (CMPLX(exp(x), y)); in cexp()
59 return (CMPLX(cos(y), sin(y))); in cexp()
64 return (CMPLX(y - y, y - y)); in cexp()
70 return (CMPLX(x, y - y)); in cexp()
89 return (CMPLX(exp_x * cos(y), exp_x * sin(y))); in cexp()
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_sincos.c29 double y[2]; in sincos() local
57 n = __ieee754_rem_pio2(x, y); in sincos()
61 __kernel_sincos(y[0], y[1], 1, sn, cs); in sincos()
64 __kernel_sincos(y[0], y[1], 1, cs, sn); in sincos()
68 __kernel_sincos(y[0], y[1], 1, sn, cs); in sincos()
73 __kernel_sincos(y[0], y[1], 1, cs, sn); in sincos()
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_sincosl.c54 long double y[2]; in sincosl() local
84 e0 = __ieee754_rem_pio2l(x, y); in sincosl()
88 __kernel_sincosl(y[0], y[1], 1, sn, cs); in sincosl()
91 __kernel_sincosl(y[0], y[1], 1, cs, sn); in sincosl()
95 __kernel_sincosl(y[0], y[1], 1, sn, cs); in sincosl()
100 __kernel_sincosl(y[0], y[1], 1, cs, sn); in sincosl()
/bionic/libm/upstream-freebsd/lib/msun/bsdsrc/
Db_tgamma.c186 double y, ym1, t; local
188 y = x - one;
189 ym1 = y - one;
190 if (y <= 1.0 + (LEFT + x0)) {
191 yy = ratfun_gam(y - x0, 0);
194 r.a = y;
197 y = ym1;
198 yy.b = r.b = y - yy.a;
200 for (ym1 = y-one; ym1 > LEFT + x0; y = ym1--, yy.a--) {
202 r.b = r.a*yy.b + y*r.b;
[all …]
/bionic/libm/
Dbuiltins.cpp55 float fmaf(float x, float y, float z) { return __builtin_fmaf(x, y, z); } in fmaf() argument
56 double fma(double x, double y, double z) { return __builtin_fma(x, y, z); } in fma() argument
58 float fmaxf(float x, float y) { return __builtin_fmaxf(x, y); } in fmaxf() argument
59 double fmax(double x, double y) { return __builtin_fmax(x, y); } in fmax() argument
61 float fminf(float x, float y) { return __builtin_fminf(x, y); } in fminf() argument
62 double fmin(double x, double y) { return __builtin_fmin(x, y); } in fmin() argument
/bionic/libc/upstream-openbsd/lib/libc/gdtoa/
Dmisc.c121 (y) ULong *y; in lo0bits() argument
123 (ULong *y)
127 ULong x = *y;
133 *y = x >> 1;
136 *y = x >> 2;
162 *y = x;
177 ULLong carry, y; local
179 ULong carry, *x, y; local
192 y = *x * (ULLong)m + carry;
193 carry = y >> 32;
[all …]
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.c120 ULLong borrow, carry, y, ys; local
122 ULong borrow, carry, y, ys; local
151 y = *bx - (ys & 0xffffffffUL) - borrow;
152 borrow = y >> 32 & 1UL;
153 *bx++ = y & 0xffffffffUL;
160 y = (*bx & 0xffff) - (ys & 0xffff) - borrow;
161 borrow = (y & 0x10000) >> 16;
164 Storeinc(bx, z, y);
168 y = *bx - (ys & 0xffff) - borrow;
169 borrow = (y & 0x10000) >> 16;
[all …]

123456