Home
last modified time | relevance | path

Searched refs:yp (Results 1 – 25 of 130) sorted by relevance

123456

/netbsd-src/external/lgpl3/gmp/dist/mpz/
H A Dfib_ui.c59 mp_ptr fp, xp, yp; in mpz_fib_ui() local
77 TMP_ALLOC_LIMBS_2 (xp,xalloc, yp,xalloc); in mpz_fib_ui()
78 size = mpn_fib2_ui (xp, yp, n2); in mpz_fib_ui()
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()
96 ysize = size + yp[size]; in mpz_fib_ui()
101 c = c2 + mpn_add_n (xp, fp, yp, size); in mpz_fib_ui()
104 c2 -= mpn_sub_n (yp, fp, yp, size); in mpz_fib_ui()
[all …]
H A Dlucnum_ui.c96 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 Dpowm.c197 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-src/external/lgpl3/gmp/dist/tests/mpn/
H A Dt-instrument.c133 mp_limb_t wp[100], xp[100], yp[100]; in check() local
137 refmpn_zero (yp, size); in check()
141 mpn_add_n (wp, xp, yp, size); in check()
146 mpn_add_nc (wp, xp, yp, size, CNST_LIMB(0)); in check()
152 mpn_addlsh1_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()
169 mpn_addmul_1 (wp, xp, size, yp[0]); in check()
174 mpn_addmul_1c (wp, xp, size, yp[0], CNST_LIMB(0)); in check()
220 notdead += mpn_hamdist (xp, yp, size); in check()
[all …]
H A Dt-fat.c80 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-src/external/lgpl3/mpfr/dist/src/
H A Dgrandom.c36 mpz_t x, y, xp, yp, t, a, b, s; in mpfr_grandom() local
59 mpz_init (yp); in mpfr_grandom()
73 mpz_set_ui (yp, 0); in mpfr_grandom()
81 mpz_urandomb (yp, rstate, tprec); in mpfr_grandom()
83 mpz_mul (b, yp, yp); in mpfr_grandom()
94 mpz_addmul_ui (s, yp, 2); 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()
[all …]
H A Dset_uj.c43 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 Dround_raw_generic.c64 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()
169 mpn_copyi (yp, xp + xsize - nw, nw); in mpfr_round_raw_generic()
170 yp[0] &= himask; in mpfr_round_raw_generic()
195 carry = mpn_add_1 (yp, xp + xsize - nw, nw, 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()
235 mpn_copyi (yp, xp + xsize - nw, nw); in mpfr_round_raw_generic()
[all …]
H A Dmul_ui.c77 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 Ddiv_ui.c42 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()
229 sb |= (yp[0] & MPFR_LIMB_MASK(sh - 1)) | w | c; 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()
[all …]
H A Dodd_p.c34 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 Dget_float128.c59 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-src/external/lgpl3/gmp/dist/mpn/generic/
H A Dmullo_n.c137 mpn_dc_mullo_n (mp_ptr rp, mp_srcptr xp, mp_srcptr yp, mp_size_t n, mp_ptr tp) in mpn_dc_mullo_n() argument
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()
177 mpn_dc_mullo_n (tp + n, xp + n2, yp, n1, tp + n); in mpn_dc_mullo_n()
182 mpn_mul_basecase (tp + n, xp, n1, yp + n2, n1); in mpn_dc_mullo_n()
184 mpn_mullo_basecase (tp + n, xp, yp + n2, n1); in mpn_dc_mullo_n()
186 mpn_dc_mullo_n (tp + n, xp, yp + n2, n1, tp + n); in mpn_dc_mullo_n()
203 mpn_mullo_n (mp_ptr rp, mp_srcptr xp, mp_srcptr yp, mp_size_t n) in mpn_mullo_n() argument
[all …]
H A Dredc_n.c52 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 Dsub_err1_n.c55 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 Dadd_err1_n.c55 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 Dbsqrtinv.c55 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-src/external/lgpl3/gmp/dist/tune/
H A Dspeed-ext.c89 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 Dspeed.h114 mp_ptr yp; /* second argument */ member
753 mp_ptr xp, yp; \
764 yp = s->yp; \
768 else if (s->r == 2) { yp = wp; } \
769 else if (s->r == 3) { xp = wp; yp = wp; } \
770 else if (s->r == 4) { yp = xp; } \
777 if (xp == wp || yp == wp) \
781 speed_operand_src (s, yp, s->size); \
801 mp_ptr xp, yp; \
818 yp = s->yp; \
[all …]
H A Dcommon.c1056 SPEED_ROUTINE_MPN_ADDSUB_N_CALL (mpn_add_n_sub_n (ap, sp, s->xp, s->yp, s->size)); in speed_mpn_add_n_sub_n()
1148 SPEED_ROUTINE_MPN_BINARY_N_CALL (mpn_addlsh_n (wp, xp, yp, s->size, 7)); in speed_mpn_addlsh_n()
1155 SPEED_ROUTINE_MPN_BINARY_N_CALL (mpn_sublsh_n (wp, xp, yp, s->size, 7)); in speed_mpn_sublsh_n()
1183 SPEED_ROUTINE_MPN_BINARY_N_CALL (mpn_rsblsh_n (wp, xp, yp, s->size, 7)); in speed_mpn_rsblsh_n()
1204 SPEED_ROUTINE_MPN_BINARY_N_CALL (mpn_cnd_add_n (1, wp, xp, yp, s->size)); in speed_mpn_cnd_add_n()
1209 SPEED_ROUTINE_MPN_BINARY_N_CALL (mpn_cnd_sub_n (1, wp, xp, yp, s->size)); in speed_mpn_cnd_sub_n()
1217 SPEED_ROUTINE_MPN_BINARY_N_CALL (mpn_and_n (wp, xp, yp, s->size)); in speed_mpn_and_n()
1222 SPEED_ROUTINE_MPN_BINARY_N_CALL (mpn_andn_n (wp, xp, yp, s->size)); in speed_mpn_andn_n()
1227 SPEED_ROUTINE_MPN_BINARY_N_CALL (mpn_nand_n (wp, xp, yp, s->size)); in speed_mpn_nand_n()
1232 SPEED_ROUTINE_MPN_BINARY_N_CALL (mpn_ior_n (wp, xp, yp, s->size)); in speed_mpn_ior_n()
[all …]
/netbsd-src/external/lgpl3/gmp/dist/tests/
H A Drefmpz.c49 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-src/external/lgpl3/gmp/dist/mpn/x86/p6/
H A Dmul_basecase.asm78 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]
210 C The first limb of yp is processed with a simple mpn_mul_1 loop running at
293 C eax yp
301 leal 4(%eax,%edx,4), %ebp C yp end
313 movl (%ebp,%edx,4), %ebp C yp second lowest limb - multiplier
[all …]
/netbsd-src/external/lgpl3/gmp/dist/mpn/x86/k7/
H A Dmul_basecase.asm85 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-src/external/lgpl3/gmp/dist/mpn/x86/k6/
H A Dmul_basecase.asm53 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-src/external/lgpl3/mpfr/dist/tune/
H A Dspeed.c97 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()

123456