/netbsd/external/lgpl3/gmp/dist/mpz/ |
H A D | fib_ui.c | 59 mp_ptr fp, xp, yp; in mpz_fib_ui() local 83 mpn_trace ("yp", yp, size)); in mpz_fib_ui() 92 yp[size] = 0; in mpz_fib_ui() 93 ASSERT_NOCARRY (mpn_add_n_sub_n (xp, yp, xp, yp, size+1)); in mpz_fib_ui() 95 ASSERT (yp[size] <= 1); in mpz_fib_ui() 104 c2 -= mpn_sub_n (yp, fp, yp, size); in mpz_fib_ui() 105 yp[size] = c2; in mpz_fib_ui() 137 c = mpn_addlsh1_n (yp, xp, yp, size); in mpz_fib_ui() 139 c = mpn_lshift (yp, yp, size, 1); in mpz_fib_ui() 140 c += mpn_add_n (yp, yp, xp, size); in mpz_fib_ui() [all …]
|
H A D | lucnum_ui.c | 96 mp_ptr yp; in mpz_lucnum_ui() local 101 yp = TMP_ALLOC_LIMBS (yalloc); in mpz_lucnum_ui() 104 xsize = mpn_fib2_ui (xp, yp, n/2); in mpz_lucnum_ui() 108 ysize -= (yp[ysize-1] == 0); in mpz_lucnum_ui() 109 ASSERT (yp[ysize-1] != 0); in mpz_lucnum_ui() 113 c = mpn_addlsh1_n (xp, yp, xp, xsize); in mpz_lucnum_ui() 116 c += mpn_add_n (xp, xp, yp, xsize); in mpz_lucnum_ui() 124 c = mpn_mul (lp, xp, xsize, yp, ysize); in mpz_lucnum_ui()
|
H A D | powm.c | 197 mp_ptr r2, xp, yp, odd_inv_2exp; in mpz_powm() local 255 yp = tp; in mpz_powm() 257 mpn_mul (yp, xp, ncnt, mp, nodd); in mpz_powm() 259 mpn_mul (yp, mp, nodd, xp, ncnt); in mpz_powm() 261 mpn_add (rp, yp, n, rp, nodd); in mpz_powm()
|
/netbsd/external/lgpl3/gmp/dist/tests/mpn/ |
H A D | t-instrument.c | 137 refmpn_zero (yp, size); in check() 141 mpn_add_n (wp, xp, yp, size); in check() 158 mpn_and_n (wp, xp, yp, size); in check() 164 mpn_andn_n (wp, xp, yp, size); in check() 225 mpn_ior_n (wp, xp, yp, size); in check() 231 mpn_iorn_n (wp, xp, yp, size); in check() 285 mpn_nand_n (wp, xp, yp, size); in check() 291 mpn_nior_n (wp, xp, yp, size); in check() 342 mpn_sub_n (wp, xp, yp, size); in check() 383 mpn_xor_n (wp, xp, yp, size); in check() [all …]
|
H A D | t-fat.c | 80 mp_limb_t wp[2], xp[2], yp[2], r; in check_functions() local 87 yp[0] = 456; in check_functions() 88 mpn_add_n (wp, xp, yp, (mp_size_t) 1); in check_functions() 203 yp[0] = 7; in check_functions() 204 mpn_mul_basecase (wp, xp, (mp_size_t) 1, yp, (mp_size_t) 1); in check_functions() 213 yp[0] = 7; in check_functions() 214 mpn_mullo_basecase (wp, xp, yp, (mp_size_t) 1); in check_functions() 264 yp[0] = 666; in check_functions() 265 mpn_sub_n (wp, xp, yp, (mp_size_t) 1); in check_functions()
|
/netbsd/external/lgpl3/mpfr/dist/src/ |
H A D | grandom.c | 59 mpz_init (yp); in mpfr_grandom() 73 mpz_set_ui (yp, 0); in mpfr_grandom() 83 mpz_mul (b, yp, yp); in mpfr_grandom() 110 mpz_mul_2exp (yp, yp, 32); in mpfr_grandom() 114 mpz_add (yp, yp, y); in mpfr_grandom() 118 mpz_mul (b, yp, yp); in mpfr_grandom() 128 mpz_mul (b, yp, yp); in mpfr_grandom() 169 mpz_mul_2exp (yp, yp, 32); in mpfr_grandom() 173 mpz_add (yp, yp, y); in mpfr_grandom() 176 mpz_mul (b, yp, yp); in mpfr_grandom() [all …]
|
H A D | set_uj.c | 43 mp_limb_t yp[uintmaxpml]; in mpfr_set_uj_2exp() local 63 MPFR_TMP_INIT1 (yp, y, uintmax_bit_size); in mpfr_set_uj_2exp() 71 yp[0] = limb << cnt; in mpfr_set_uj_2exp() 80 yp[i] = j; /* Only the low bits are copied */ in mpfr_set_uj_2exp() 86 limb = yp[--k]; in mpfr_set_uj_2exp() 90 len = numberof (yp) - k; in mpfr_set_uj_2exp() 96 mpn_lshift (yp+len, yp, k, cnt); /* Normalize the high limb */ in mpfr_set_uj_2exp() 98 mpn_copyd (yp+len, yp, k); /* Must use copyd */ in mpfr_set_uj_2exp() 102 yp[0] = MPFR_LIMB_ZERO; in mpfr_set_uj_2exp() 104 MPN_ZERO (yp, len); /* Zero the last limbs */ in mpfr_set_uj_2exp() [all …]
|
H A D | round_raw_generic.c | 64 mp_limb_t *yp, in mpfr_round_raw_generic() argument 120 mpn_copyd (yp + (nw - xsize), xp, xsize); in mpfr_round_raw_generic() 121 MPN_ZERO(yp, nw - xsize); in mpfr_round_raw_generic() 170 yp[0] &= himask; in mpfr_round_raw_generic() 199 yp[0] &= himask; in mpfr_round_raw_generic() 217 mpn_copyi (yp, xp + xsize - nw, nw); in mpfr_round_raw_generic() 218 yp[0] &= himask; in mpfr_round_raw_generic() 236 yp[0] &= himask; in mpfr_round_raw_generic() 252 yp[0] &= himask; in mpfr_round_raw_generic() 269 mpn_copyi (yp, xp + xsize - nw, nw); in mpfr_round_raw_generic() [all …]
|
H A D | mul_ui.c | 77 mp_limb_t *yp; in mpfr_mul_ui() local 82 yp = MPFR_MANT (y); in mpfr_mul_ui() 87 yp = MPFR_TMP_LIMBS_ALLOC (xn + 1); in mpfr_mul_ui() 90 yp[xn] = mpn_mul_1 (yp, MPFR_MANT (x), xn, u); in mpfr_mul_ui() 96 MPFR_ASSERTD (yp[xn] != 0); in mpfr_mul_ui() 97 if (MPFR_LIKELY (MPFR_LIMB_MSB (yp[xn]) == 0)) in mpfr_mul_ui() 99 count_leading_zeros (cnt, yp[xn]); in mpfr_mul_ui() 100 mpn_lshift (yp, yp, xn + 1, cnt); in mpfr_mul_ui() 109 MPFR_RNDRAW (inexact, y, yp, (mpfr_prec_t) (xn + 1) * GMP_NUMB_BITS, in mpfr_mul_ui()
|
H A D | div_ui.c | 42 mp_limb_t *xp, *yp, *tmp, c, d; in mpfr_div_ui() local 107 yp = MPFR_MANT (y); in mpfr_div_ui() 158 MPN_COPY(yp, tmp, yn); in mpfr_div_ui() 216 mpn_lshift (yp, tmp + 1, yn, shlz); in mpfr_div_ui() 217 yp[0] |= tmp[0] >> (GMP_NUMB_BITS - shlz); in mpfr_div_ui() 228 rb = yp[0] & (MPFR_LIMB_ONE << (sh - 1)); in mpfr_div_ui() 236 d = yp[0] & MPFR_LIMB_MASK (sh); in mpfr_div_ui() 237 yp[0] ^= d; /* clear the lowest sh bits */ in mpfr_div_ui() 280 if (yp[0] & (MPFR_LIMB_ONE << sh)) in mpfr_div_ui() 293 MPFR_UNLIKELY (mpn_add_1 (yp, yp, yn, MPFR_LIMB_ONE << sh))) in mpfr_div_ui() [all …]
|
H A D | odd_p.c | 34 mp_limb_t *yp; in mpfr_odd_p() local 63 yp = MPFR_MANT(y); in mpfr_odd_p() 65 if (expo % GMP_NUMB_BITS == 0 ? (yp[yn] & 1) == 0 in mpfr_odd_p() 66 : MPFR_LIMB_LSHIFT(yp[yn], (expo % GMP_NUMB_BITS) - 1) != MPFR_LIMB_HIGHBIT) in mpfr_odd_p() 69 if (yp[yn] != 0) in mpfr_odd_p()
|
H A D | get_float128.c | 59 mp_limb_t *yp; in mpfr_get_float128() local 77 yp = MPFR_MANT (y); in mpfr_get_float128() 85 r += yp[i]; in mpfr_get_float128()
|
/netbsd/external/lgpl3/gmp/dist/mpn/generic/ |
H A D | mullo_n.c | 142 ASSERT (! MPN_OVERLAP_P (rp, n, yp, n)); in mpn_dc_mullo_n() 168 mpn_mul_n (tp, xp, yp, n2); in mpn_dc_mullo_n() 173 mpn_mul_basecase (tp + n, xp + n2, n1, yp, n1); in mpn_dc_mullo_n() 175 mpn_mullo_basecase (tp + n, xp + n2, yp, n1); in mpn_dc_mullo_n() 184 mpn_mullo_basecase (tp + n, xp, yp + n2, n1); in mpn_dc_mullo_n() 207 ASSERT (! MPN_OVERLAP_P (rp, n, yp, n)); in mpn_mullo_n() 213 mpn_mul_basecase (tp, xp, n, yp, n); in mpn_mullo_n() 218 mpn_mullo_basecase (rp, xp, yp, n); in mpn_mullo_n() 228 mpn_dc_mullo_n (rp, xp, yp, n, tp); in mpn_mullo_n() 235 mpn_fft_mul (tp, xp, n, yp, n); in mpn_mullo_n() [all …]
|
H A D | redc_n.c | 52 mp_ptr xp, yp, scratch; in mpn_redc_n() local 67 yp = scratch + n; in mpn_redc_n() 68 mpn_mulmod_bnm1 (yp, rn, xp, n, mp, n, scratch + n + rn); in mpn_redc_n() 72 cy = mpn_sub_n (yp + rn, yp, up, 2*n - rn); /* undo wrap around */ in mpn_redc_n() 73 MPN_DECR_U (yp + 2*n - rn, rn, cy); in mpn_redc_n() 75 cy = mpn_sub_n (rp, up + n, yp + n, n); in mpn_redc_n()
|
H A D | add_err1_n.c | 55 mp_ptr ep, mp_srcptr yp, in mpn_add_err1_n() argument 63 ASSERT (! MPN_OVERLAP_P (rp, n, yp, n)); in mpn_add_err1_n() 66 ASSERT (! MPN_OVERLAP_P (ep, 2, yp, n)); in mpn_add_err1_n() 69 yp += n - 1; in mpn_add_err1_n() 74 yl = *yp--; in mpn_add_err1_n()
|
H A D | sub_err1_n.c | 55 mp_ptr ep, mp_srcptr yp, in mpn_sub_err1_n() argument 63 ASSERT (! MPN_OVERLAP_P (rp, n, yp, n)); in mpn_sub_err1_n() 66 ASSERT (! MPN_OVERLAP_P (ep, 2, yp, n)); in mpn_sub_err1_n() 69 yp += n - 1; in mpn_sub_err1_n() 74 yl = *yp--; in mpn_sub_err1_n()
|
H A D | bsqrtinv.c | 55 mpn_bsqrtinv (mp_ptr rp, mp_srcptr yp, mp_bitcnt_t bnb, mp_ptr tp) in mpn_bsqrtinv() argument 70 if ((yp[0] & 3) != 1) in mpn_bsqrtinv() 75 if ((yp[0] & 7) != 1) in mpn_bsqrtinv() 92 mpn_mullo_n (rp, yp, tp2, bn); in mpn_bsqrtinv()
|
/netbsd/external/lgpl3/gmp/dist/tune/ |
H A D | speed-ext.c | 89 mean_calls (mp_ptr wp, mp_srcptr xp, mp_srcptr yp, mp_size_t size) in mean_calls() argument 95 c = mpn_add_n (wp, xp, yp, size); in mean_calls() 107 mean_open (mp_ptr wp, mp_srcptr xp, mp_srcptr yp, mp_size_t size) in mean_open() argument 115 y = yp[0]; in mean_open() 126 y = yp[i]; in mean_open() 148 mean_open2 (mp_ptr wp, mp_srcptr xp, mp_srcptr yp, mp_size_t size) in mean_open2() argument 156 y = yp[0]; in mean_open2() 165 ynext = yp[i+1]; in mean_open2() 202 speed_operand_src (s, s->yp, s->size); \ 209 mean_fun (wp, s->xp, s->yp, s->size); \
|
H A D | speed.h | 753 mp_ptr xp, yp; \ 764 yp = s->yp; \ 818 yp = s->yp; \ 878 yp = s->yp; \ 1031 mp_limb_t yp[N]; \ 2719 yp = s->yp; \ 2731 if (yp != s->yp) MPN_COPY (yp, s->yp, s->size); \ 3273 mp_limb_t yp[2]; \ 3375 mp_limb_t yp[2]; \ 3393 if (yp[1] == 0) \ [all …]
|
/netbsd/external/bsd/tmux/dist/ |
H A D | cmd-display-menu.c | 104 const char *xp, *yp; in cmd_display_menu_get_position() local 252 yp = args_get(args, 'y'); in cmd_display_menu_get_position() 253 if (yp == NULL || strcmp(yp, "C") == 0) in cmd_display_menu_get_position() 254 yp = "#{popup_centre_y}"; in cmd_display_menu_get_position() 255 else if (strcmp(yp, "P") == 0) in cmd_display_menu_get_position() 256 yp = "#{popup_pane_bottom}"; in cmd_display_menu_get_position() 257 else if (strcmp(yp, "M") == 0) in cmd_display_menu_get_position() 258 yp = "#{popup_mouse_top}"; in cmd_display_menu_get_position() 259 else if (strcmp(yp, "S") == 0) in cmd_display_menu_get_position() 261 else if (strcmp(yp, "W") == 0) in cmd_display_menu_get_position() [all …]
|
/netbsd/external/lgpl3/gmp/dist/tests/ |
H A D | refmpz.c | 49 mp_ptr xp, yp; in refmpz_hamdist() local 64 yp = refmpn_malloc_limbs (tsize); in refmpz_hamdist() 65 refmpn_zero (yp, tsize); in refmpz_hamdist() 66 refmpn_copy (yp, PTR(y), ysize); in refmpz_hamdist() 72 refmpn_neg (yp, yp, tsize); in refmpz_hamdist() 74 ret = refmpn_hamdist (xp, yp, tsize); in refmpz_hamdist() 77 free (yp); in refmpz_hamdist()
|
/netbsd/external/lgpl3/gmp/dist/mpn/x86/p6/ |
H A D | mul_basecase.asm | 78 movl (%eax), %eax C yp[0] 159 mull %ecx C xp[0] * yp[0] 167 mull %ecx C xp[1] * yp[0] 173 movl 4(%ecx), %ecx C yp[1] 179 mull %ecx C xp[1] * yp[1] 187 mull %ecx C xp[0] * yp[1] 293 C eax yp 301 leal 4(%eax,%edx,4), %ebp C yp end 328 movl (%edx,%ebp,4), %ebp C yp limb - multiplier 510 movl (%edx,%ebp,4), %ebp C yp next multiplier [all …]
|
/netbsd/external/lgpl3/gmp/dist/mpn/x86/k7/ |
H A D | mul_basecase.asm | 85 movl (%eax), %eax C yp low limb 107 movl %eax, %ecx C yp low limb 165 mull %ecx C xp[1] * yp[0] 171 movl 4(%ecx), %ecx C yp[1] 177 mull %ecx C xp[1] * yp[1] 186 mull %ecx C xp[0] * yp[1] 207 C The first limb of yp is processed with a simple mpn_mul_1 style loop 298 C eax yp 306 leal 4(%eax,%edx,4), %ebp C yp end 318 movl (%ebp,%edx,4), %ebp C yp second lowest limb - multiplier [all …]
|
/netbsd/external/lgpl3/gmp/dist/mpn/x86/k6/ |
H A D | mul_basecase.asm | 53 C mp_srcptr yp, mp_size_t ysize); 55 C Calculate xp,xsize multiplied by yp,ysize, storing the result in 165 mull %ecx C xp[0] * yp[0] 223 C eax yp[0] 400 C added to give the location of the next limb of yp, which is the multiplier 405 C point is one limb upwards for each limb of yp processed). 480 movl 4(%ebx), %ebp C multiplier (yp second limb) 482 leal 4(%ebx,%edx,4), %ebx C yp adjust for ysize indexing
|
/netbsd/external/lgpl3/mpfr/dist/tune/ |
H A D | speed.c | 97 s->yp[size-1] |= MPFR_LIMB_HIGHBIT; \ 98 MPFR_TMP_INIT1 (s->yp, y, s->size); \ 104 speed_operand_src (s, s->yp, size); \ 154 s.yp = s.xp + size; in domeasure() 155 mpn_random (s.yp, size); in domeasure()
|