mpn/generic/toom_interpolate_16pts.c

288 			mp_size_t n, mp_size_t spt, int half, mp_ptr wsi)
307     cy = DO_mpn_sublsh_n (r3, r0, spt, 14, wsi);
309     DO_mpn_subrsh(r6, n3p1, r0, spt, 2, wsi);
311     cy = DO_mpn_sublsh_n (r2, r0, spt, 28, wsi);
313     DO_mpn_subrsh(r5, n3p1, r0, spt, 4, wsi);
315     cy = DO_mpn_sublsh_n (r1 + BIT_CORRECTION, r0, spt, 42 - CORRECTION_BITS, wsi);
326     DO_mpn_subrsh(r7, n3p1 + BIT_CORRECTION, r0, spt, 6, wsi);
334   r5[n3] -= DO_mpn_sublsh_n (r5 + n, pp, 2 * n, 28, wsi);
335   DO_mpn_subrsh(r2 + n, 2 * n + 1, pp, 2 * n, 4, wsi);
340   mpn_sub_n (wsi, r5, r2, n3p1); /* can be negative */
342   MP_PTR_SWAP(r5, wsi);
345   r6[n3] -= DO_mpn_sublsh_n (r6 + n, pp, 2 * n, 14, wsi);
346   DO_mpn_subrsh(r3 + n, 2 * n + 1, pp, 2 * n, 2, wsi);
351   ASSERT_NOCARRY(mpn_add_n (wsi, r3, r6, n3p1));
353   MP_PTR_SWAP(r3, wsi);
356   cy = DO_mpn_sublsh_n (r7 + n + BIT_CORRECTION, pp, 2 * n, 42 - CORRECTION_BITS, wsi);
359   cy = DO_mpn_sublsh_n (r1 + n, pp + 1, 2 * n - 1, GMP_NUMB_BITS - 6, wsi);
364   DO_mpn_subrsh(r1 + n, 2 * n + 1, pp, 2 * n, 6, wsi);
370   mpn_sub_n (wsi, r7, r1, n3p1); /* can be negative */
372   MP_PTR_SWAP(r7, wsi);
380   DO_mpn_sublsh_n (r5, r6, n3p1, 2, wsi); /* can be negative */
381   DO_mpn_sublsh_n (r5, r6, n3p1,10, wsi); /* can be negative */
388   DO_mpn_sublsh_n (r7, r6, n3p1, 4, wsi); /* can be negative */
389   DO_mpn_sublsh_n (r7, r6, n3p1,12, wsi); /* can be negative */
390   DO_mpn_sublsh_n (r7, r6, n3p1,20, wsi); /* can be negative */
404   DO_mpn_sublsh_n (r6, r7, n3p1, 12, wsi); /* can be negative */
409   DO_mpn_addlsh_n (r6, r5, n3p1, 8, wsi); /* can give a carry */
410   DO_mpn_sublsh_n (r6, r5, n3p1, 4, wsi); /* can be negative */
417   ASSERT_NOCARRY(DO_mpn_sublsh_n (r3, r4, n3p1, 7, wsi));
419   ASSERT_NOCARRY(DO_mpn_sublsh_n (r2, r4, n3p1, 13, wsi));
423   DO_mpn_sublsh_n (r1, r4, n3p1, 19, wsi);
435   ASSERT_NOCARRY(DO_mpn_addlsh_n (r3, r1, n3p1, 7, wsi));
436   ASSERT_NOCARRY(DO_mpn_sublsh_n (r3, r1, n3p1, 12, wsi));