| /src/lib/libm/src/ |
| w_j0.c | 19 * wrapper j0(double x), y0(double x)
41 y0(double x) /* wrapper y0 */ function
58 return __kernel_standard(x,x,35); /* y0(x>X_TLOSS) */
|
| /src/crypto/external/bsd/heimdal/dist/lib/hcrypto/libtommath/ |
| bn_mp_karatsuba_mul.c | 51 mp_int x0, x1, y0, y1, t1, x0y0, x1y1; local
68 if (mp_init_size (&y0, B) != MP_OKAY)
71 goto Y0;
82 x0.used = y0.used = B;
97 tmpy = y0.dp;
118 mp_clamp (&y0);
122 if (mp_mul (&x0, &y0, &x0y0) != MP_OKAY)
123 goto X1Y1; /* x0y0 = x0*y0 */
127 /* now calc x1+x0 and y1+y0 */
130 if (s_mp_add (&y1, &y0, &x0) != MP_OKAY
[all...] |
| /src/external/lgpl3/gmp/dist/mpn/generic/ |
| brootinv.c | 83 mp_limb_t kinv, k2, r0, y0; local
104 y0 = yp[0];
106 r0 = y0 ^ (((y0 << 1) ^ (y0 << 2)) & (k2 << 3) & 8); /* 4 bits */
107 r0 = kinv * (k2 * r0 * 2 - y0 * powsquaredlimb(r0, k2 & 0x3f)); /* 8 bits */
108 r0 = kinv * (k2 * r0 * 2 - y0 * powsquaredlimb(r0, k2 & 0x3fff)); /* 16 bits */
114 r0 = kinv * (k2 * r0 * 2 - y0 * powsquaredlimb(r0, k2));
|
| /src/sys/arch/m68k/fpe/ |
| fpu_div.c | 157 uint32_t r0, r1, r2, d0, d1, d2, y0, y1, y2; local
209 FPU_SUBCS(d1, r1, y1); FPU_SUBC(d0, r0, y0)
244 y0 = y->fp_mant[0];
|
| fpu_sqrt.c | 194 uint32_t y0, y1, y2; local
261 * The calculation `t = y + bit' breaks down into `t0 = y0, ...,
272 /* if (x >= (t0 = y0 | bit)) { */ /* always true */
275 y0 = bit << 1;
280 t0 = y0 | bit; /* t = y + bit */
284 y0 |= bit << 1; /* y += bit << 1 */
291 /* calculate q1. note (y0&1)==0. */
292 #define t0 y0
304 y0 |= 1; /* y += bit << 1 */
322 /* calculate q2. note (y1&1)==0; y0 (aka t0) is fixed. *
[all...] |
| /src/sys/arch/powerpc/fpu/ |
| fpu_div.c | 161 u_int r0, r1, r2, r3, d0, d1, d2, d3, y0, y1, y2, y3; local
234 FPU_SUBCS(d1, r1, y1); FPU_SUBC(d0, r0, y0)
270 y0 = y->fp_mant[0];
|
| fpu_sqrt.c | 198 u_int y0, y1, y2, y3; local
281 * The calculation `t = y + bit' breaks down into `t0 = y0, ...,
292 /* if (x >= (t0 = y0 | bit)) { */ /* always true */
295 y0 = bit << 1;
300 t0 = y0 | bit; /* t = y + bit */
304 y0 |= bit << 1; /* y += bit << 1 */
311 /* calculate q1. note (y0&1)==0. */
312 #define t0 y0
324 y0 |= 1; /* y += bit << 1 */
342 /* calculate q2. note (y1&1)==0; y0 (aka t0) is fixed. *
[all...] |
| /src/sys/arch/sparc/fpu/ |
| fpu_div.c | 157 u_int r0, r1, r2, r3, d0, d1, d2, d3, y0, y1, y2, y3; local
205 FPU_SUBCS(d1, r1, y1); FPU_SUBC(d0, r0, y0)
241 y0 = y->fp_mant[0];
|
| fpu_sqrt.c | 194 u_int y0, y1, y2, y3; local
261 * The calculation `t = y + bit' breaks down into `t0 = y0, ...,
272 /* if (x >= (t0 = y0 | bit)) { */ /* always true */
275 y0 = bit << 1;
280 t0 = y0 | bit; /* t = y + bit */
284 y0 |= bit << 1; /* y += bit << 1 */
291 /* calculate q1. note (y0&1)==0. */
292 #define t0 y0
304 y0 |= 1; /* y += bit << 1 */
322 /* calculate q2. note (y1&1)==0; y0 (aka t0) is fixed. *
[all...] |
| /src/external/gpl3/gcc/dist/libgcc/config/arc/ieee-754/ |
| divtab-arc-df.c | 55 long double y0 = a0 + (x-1)*a1; local
57 long double approx = 2. * y0 - y0 * x * y0;
|
| /src/external/gpl3/gcc.old/dist/libgcc/config/arc/ieee-754/ |
| divtab-arc-df.c | 55 long double y0 = a0 + (x-1)*a1; local
57 long double approx = 2. * y0 - y0 * x * y0;
|
| /src/external/lgpl3/mpfr/dist/tests/ |
| tlog10.c | 59 mpfr_t t, x, y0, y1, y2; local
69 /* The precisions of t and y0 should be large enough to avoid
71 mpfr_inits2 (64, t, y0, (mpfr_ptr) 0);
72 mpfr_set_exp_t (y0, e, MPFR_RNDN);
74 mpfr_div (y0, y0, t, MPFR_RNDN);
76 mpfr_mul (y0, y0, t, MPFR_RNDN);
81 mpfr_set (y1, y0, MPFR_RNDN);
100 mpfr_clears (t, y0, (mpfr_ptr) 0)
[all...] |
| tfrac.c | 176 mp_limb_t y0; local
195 y0 = MPFR_MANT(y)[0];
196 while (y0 != 0 && (y0 >> 1) << 1 == y0)
197 y0 >>= 1;
198 if (y0 > 0x2000)
|
| /src/sys/crypto/aes/arch/x86/ |
| aes_sse2_4x32_dec.c | 112 __m128i x, y0, y1, y2, y3, y4, y5, y6; local
115 y0 = x & _mm_set1_epi32(0x000000FF);
122 return y0 | y1 | y2 | y3 | y4 | y5 | y6;
|
| aes_sse2_4x32_enc.c | 51 __m128i x, y0, y1, y2, y3, y4, y5, y6; local
54 y0 = x & _mm_set1_epi32(0x000000FF);
61 return y0 | y1 | y2 | y3 | y4 | y5 | y6;
|
| /src/games/hack/ |
| hack.c | 566 int x, y, i, x0 = 0, y0 = 0, m0 = 0, i0 = 9; local
610 if (corrct == 1 && DIST(x, y, x0, y0) != 1)
615 y0 = y;
647 if (u.dx == y0 - u.uy && u.dy == u.ux - x0)
652 if ((u.dx == u.dy && y0 == u.uy) ||
653 (u.dx != u.dy && y0 != u.uy))
658 if ((x0 - u.ux == y0 - u.uy && !u.dy) ||
659 (x0 - u.ux != y0 - u.uy && u.dy))
667 u.dx = x0 - u.ux, u.dy = y0 - u.uy;
|
| /src/lib/libm/noieee_src/ |
| n_j0.c | 67 /* double j0(double x), y0(double x)
92 * Method -- y0(x):
95 * y0(x) = 2/pi*(j0(x)*(ln(x/2)+Euler) + x^2/4 - ...)
96 * therefore y0(x)-2/pi*j0(x)*ln(x) is an even function.
97 * We use the following function to approximate y0,
98 * y0(x) = U(z)/V(z) + (2/pi)*(j0(x)*ln(x)), z= x^2
104 * y0(tiny) = u0 + (2/pi)*ln(tiny), (choose tiny<2**-27)
106 * y0(x) = sqrt(2/(pi*x))*(p0(x)*cos(x0)+q0(x)*sin(x0))
109 * 3. Special cases: y0(0)=-inf, y0(x<0)=NaN, y0(inf)=0
211 y0(double x) function
[all...] |
| /src/sys/crypto/chacha/ |
| chacha_ref.c | 72 uint32_t y0,y1,y2,y3,y4,y5,y6,y7,y8,y9,y10,y11,y12,y13,y14,y15;
74 x0 = y0 = le32dec(c + 0);
92 CHACHA_QUARTERROUND( y0, y4, y8,y12);
96 CHACHA_QUARTERROUND( y0, y5,y10,y15);
102 le32enc(out + 0, x0 + y0);
132 uint32_t y0,y1,y2,y3,y4,y5,y6,y7,y8,y9,y10,y11,y12,y13,y14,y15;
153 y0 = x0;
170 CHACHA_QUARTERROUND( y0, y4, y8,y12);
174 CHACHA_QUARTERROUND( y0, y5,y10,y15);
179 le32enc(s + 0, x0 + y0);
71 uint32_t y0,y1,y2,y3,y4,y5,y6,y7,y8,y9,y10,y11,y12,y13,y14,y15; local
130 uint32_t y0,y1,y2,y3,y4,y5,y6,y7,y8,y9,y10,y11,y12,y13,y14,y15; local
253 uint32_t y0,y1,y2,y3,y4,y5,y6,y7,y8,y9,y10,y11,y12,y13,y14,y15; local
376 uint32_t y0,y1,y2,y3,y4,y5,y6,y7,y8,y9,y10,y11,y12,y13,y14,y15; local
[all...] |
| /src/crypto/external/apache2/openssl/dist/crypto/ec/ |
| ec2_smpl.c | 357 BIGNUM *x0, *y0, *x1, *y1, *x2, *y2, *s, *t; local
385 y0 = BN_CTX_get(ctx);
398 if (!BN_copy(y0, a->Y))
401 if (!EC_POINT_get_affine_coordinates(group, a, x0, y0, ctx))
417 if (!BN_GF2m_add(s, y0, y1))
430 if (BN_GF2m_cmp(y0, y1) || BN_is_zero(x1)) {
|
| /src/crypto/external/bsd/openssl/dist/crypto/ec/ |
| ec2_smpl.c | 359 BIGNUM *x0, *y0, *x1, *y1, *x2, *y2, *s, *t; local
387 y0 = BN_CTX_get(ctx);
400 if (!BN_copy(y0, a->Y))
403 if (!EC_POINT_get_affine_coordinates(group, a, x0, y0, ctx))
419 if (!BN_GF2m_add(s, y0, y1))
432 if (BN_GF2m_cmp(y0, y1) || BN_is_zero(x1)) {
|
| /src/crypto/external/bsd/openssl.old/dist/crypto/ec/ |
| ec2_smpl.c | 353 BIGNUM *x0, *y0, *x1, *y1, *x2, *y2, *s, *t; local
376 y0 = BN_CTX_get(ctx);
389 if (!BN_copy(y0, a->Y))
392 if (!EC_POINT_get_affine_coordinates(group, a, x0, y0, ctx))
408 if (!BN_GF2m_add(s, y0, y1))
421 if (BN_GF2m_cmp(y0, y1) || BN_is_zero(x1)) {
|
| /src/external/gpl3/gcc/dist/libquadmath/math/ |
| lgammaq_neg.c | 499 __float128 y0 = 1 - x0_hi; local
500 __float128 y0_eps = -x0_hi + (1 - y0) - x0_lo;
504 = log (gamma (Y0 + Y0_EPS) / gamma (Y + Y_EPS)). XDIFF
505 accurately approximates the difference Y0 + Y0_EPS - Y -
514 ny0 = y0 + n_up;
515 ny0_eps = y0 - (ny0 - n_up) + y0_eps;
516 y0 = ny0;
526 = (xdiff * log1pq ((y0 - e_hi - e_lo + y0_eps) / e_hi)
528 /* Compute the sum of (B_2k / 2k(2k-1))(Y0^-(2k-1) - Y^-(2k-1)). */
529 __float128 y0r = 1 / y0, yr = 1 / y
[all...] |
| /src/external/gpl3/gcc.old/dist/gcc/config/sh/ |
| divtab.cc | 87 double y0 = (constant - (int) floor ((x * factor * 64.))) / 16384.; local
88 double y1 = 2 * y0 -y0 * y0 * (x + BIAS / (1.*(1LL<<30)));
107 y2d = - floor ((y1 - y0 / (1<<30-14)) * y2d1) / (1LL<<44);
|
| /src/external/gpl3/gcc.old/dist/libquadmath/math/ |
| lgammaq_neg.c | 499 __float128 y0 = 1 - x0_hi; local
500 __float128 y0_eps = -x0_hi + (1 - y0) - x0_lo;
504 = log (gamma (Y0 + Y0_EPS) / gamma (Y + Y_EPS)). XDIFF
505 accurately approximates the difference Y0 + Y0_EPS - Y -
514 ny0 = y0 + n_up;
515 ny0_eps = y0 - (ny0 - n_up) + y0_eps;
516 y0 = ny0;
526 = (xdiff * log1pq ((y0 - e_hi - e_lo + y0_eps) / e_hi)
528 /* Compute the sum of (B_2k / 2k(2k-1))(Y0^-(2k-1) - Y^-(2k-1)). */
529 __float128 y0r = 1 / y0, yr = 1 / y
[all...] |
| /src/external/lgpl3/mpfr/dist/src/ |
| strtofr.c | 490 mp_limb_t *y0, *y; local
503 y0 = MPFR_TMP_LIMBS_ALLOC (2 * ysize + extra_limbs + 1);
504 y = y0 + ysize; /* y has (ysize + extra_limbs + 1) allocated limbs */
681 z = y0;
741 /* y0 = y * K^ysize */
742 MPN_ZERO (y0, ysize);
799 MPFR_ASSERTD (MPFR_LIMB_MSB (y0[2 * ysize - 1]) != 0);
801 mpn_tdiv_qr (result + ysize, result, (mp_size_t) 0, y0,
|
