arch/generic/mul_n.c

54 impn_mul_n_basecase (mp_ptr prodp, mp_srcptr up, mp_srcptr vp, mp_size_t size)
56 impn_mul_n_basecase (prodp, up, vp, size)
58      mp_srcptr up;
73 	MPN_COPY (prodp, up, size);
79     cy_limb = mpn_mul_1 (prodp, up, size, v_limb);
93 	    cy_limb = mpn_add_n (prodp, prodp, up, size);
96 	cy_limb = mpn_addmul_1 (prodp, up, size, v_limb);
106 	     mp_srcptr up, mp_srcptr vp, mp_size_t size, mp_ptr tspace)
108 impn_mul_n (prodp, up, vp, size, tspace)
110      mp_srcptr up;
131       MPN_MUL_N_RECURSE (prodp, up, vp, esize, tspace);
132       cy_limb = mpn_addmul_1 (prodp + esize, up, esize, vp[esize]);
134       cy_limb = mpn_addmul_1 (prodp + esize, vp, size, up[esize]);
163       MPN_MUL_N_RECURSE (prodp + size, up + hsize, vp + hsize, hsize, tspace);
167       if (mpn_cmp (up + hsize, up, hsize) >= 0)
169 	  mpn_sub_n (prodp, up + hsize, up, hsize);
174 	  mpn_sub_n (prodp, up, up + hsize, hsize);
207       MPN_MUL_N_RECURSE (tspace, up, vp, hsize, tspace + size);
224 impn_sqr_n_basecase (mp_ptr prodp, mp_srcptr up, mp_size_t size)
226 impn_sqr_n_basecase (prodp, up, size)
228      mp_srcptr up;
238   v_limb = up[0];
242 	MPN_COPY (prodp, up, size);
248     cy_limb = mpn_mul_1 (prodp, up, size, v_limb);
257       v_limb = up[i];
262 	    cy_limb = mpn_add_n (prodp, prodp, up, size);
265 	cy_limb = mpn_addmul_1 (prodp, up, size, v_limb);
275 	     mp_srcptr up, mp_size_t size, mp_ptr tspace)
277 impn_sqr_n (prodp, up, size, tspace)
279      mp_srcptr up;
299       MPN_SQR_N_RECURSE (prodp, up, esize, tspace);
300       cy_limb = mpn_addmul_1 (prodp + esize, up, esize, up[esize]);
302       cy_limb = mpn_addmul_1 (prodp + esize, up, size, up[esize]);
315       MPN_SQR_N_RECURSE (prodp + size, up + hsize, hsize, tspace);
319       if (mpn_cmp (up + hsize, up, hsize) >= 0)
321 	  mpn_sub_n (prodp, up + hsize, up, hsize);
325 	  mpn_sub_n (prodp, up, up + hsize, hsize);
345       MPN_SQR_N_RECURSE (tspace, up, hsize, tspace + size);
363 mpn_mul_n (mp_ptr prodp, mp_srcptr up, mp_srcptr vp, mp_size_t size)
365 mpn_mul_n (prodp, up, vp, size)
367      mp_srcptr up;
374   if (up == vp)
378 	  impn_sqr_n_basecase (prodp, up, size);
384 	  impn_sqr_n (prodp, up, size, tspace);
391 	  impn_mul_n_basecase (prodp, up, vp, size);
397 	  impn_mul_n (prodp, up, vp, size, tspace);