mpn/generic/toom_interpolate_16pts.c

40 #define DO_mpn_sublsh_n(dst,src,n,s,ws) mpn_sublsh_n(dst,dst,src,n,s)
43 DO_mpn_sublsh_n(mp_ptr dst, mp_srcptr src, mp_size_t n, unsigned int s, mp_ptr ws)
46   return mpn_submul_1(dst,src,n,CNST_LIMB(1) <<(s));
50   return    __cy + mpn_sub_n(dst,dst,ws,n);
56 #define DO_mpn_addlsh_n(dst,src,n,s,ws) mpn_addlsh_n(dst,dst,src,n,s)
59 DO_mpn_addlsh_n(mp_ptr dst, mp_srcptr src, mp_size_t n, unsigned int s, mp_ptr ws)
62   return mpn_addmul_1(dst,src,n,CNST_LIMB(1) <<(s));
66   return    __cy + mpn_add_n(dst,dst,ws,n);
72 #define DO_mpn_subrsh(dst,nd,src,ns,s,ws) mpn_subrsh(dst,nd,src,ns,s)
75 #define DO_mpn_subrsh(dst,nd,src,ns,s,ws)				\
78   MPN_DECR_U (dst, nd, src[0] >> s);					\
79   __cy = DO_mpn_sublsh_n (dst, src + 1, ns - 1, GMP_NUMB_BITS - s, ws);	\
80   MPN_DECR_U (dst + ns - 1, nd - ns + 1, __cy);				\
141 #define mpn_divexact_by255(dst,src,size) \
142   (255 & 1 * mpn_bdiv_dbm1 (dst, src, size, __GMP_CAST (mp_limb_t, GMP_NUMB_MASK / 255)))
145 #define mpn_divexact_by255(dst,src,size) mpn_pi1_bdiv_q_1(dst,src,size,CNST_LIMB(255),BINVERT_255,0)
147 #define mpn_divexact_by255(dst,src,size) mpn_divexact_1(dst,src,size,CNST_LIMB(255))
154 #define mpn_divexact_by255x4(dst,src,size) mpn_pi1_bdiv_q_1(dst,src,size,CNST_LIMB(255),BINVERT_255,2)
156 #define mpn_divexact_by255x4(dst,src,size) mpn_divexact_1(dst,src,size,CNST_LIMB(255)<<2)
162 #define mpn_divexact_by9x16(dst,src,size) mpn_pi1_bdiv_q_1(dst,src,size,CNST_LIMB(9),BINVERT_9,4)
164 #define mpn_divexact_by9x16(dst,src,size) mpn_divexact_1(dst,src,size,CNST_LIMB(9)<<4)
170 #define mpn_divexact_by42525x16(dst,src,size) mpn_pi1_bdiv_q_1(dst,src,size,CNST_LIMB(42525),BINVERT_42525,4)
172 #define mpn_divexact_by42525x16(dst,src,size) mpn_divexact_1(dst,src,size,CNST_LIMB(42525)<<4)
178 #define mpn_divexact_by2835x64(dst,src,size) mpn_pi1_bdiv_q_1(dst,src,size,CNST_LIMB(2835),BINVERT_2835,6)
180 #define mpn_divexact_by2835x64(dst,src,size) mpn_divexact_1(dst,src,size,CNST_LIMB(2835)<<6)
190 #define mpn_divexact_by255x182712915(dst,src,size)				\
192     mpn_pi1_bdiv_q_1(dst,src,size,CNST_LIMB(182712915),BINVERT_182712915,0);	\
193     mpn_divexact_by255(dst,dst,size);						\
196 #define mpn_divexact_by255x182712915(dst,src,size)	\
198     mpn_divexact_1(dst,src,size,CNST_LIMB(182712915));	\
199     mpn_divexact_by255(dst,dst,size);			\
204 #define mpn_divexact_by255x182712915(dst,src,size) \
205   mpn_pi1_bdiv_q_1(dst,src,size,255*CNST_LIMB(182712915),BINVERT_255x182712915,0)
207 #define mpn_divexact_by255x182712915(dst,src,size) mpn_divexact_1(dst,src,size,255*CNST_LIMB(182712915))
218 #define mpn_divexact_by255x188513325(dst,src,size)			\
220     mpn_pi1_bdiv_q_1(dst,src,size,CNST_LIMB(188513325),BINVERT_188513325,0);	\
221     mpn_divexact_by255(dst,dst,size);					\
224 #define mpn_divexact_by255x188513325(dst,src,size)	\
226     mpn_divexact_1(dst,src,size,CNST_LIMB(188513325));	\
227     mpn_divexact_by255(dst,dst,size);			\
232 #define mpn_divexact_by255x188513325(dst,src,size) \
233   mpn_pi1_bdiv_q_1(dst,src,size,255*CNST_LIMB(188513325),BINVERT_255x188513325,0)
235 #define mpn_divexact_by255x188513325(dst,src,size) mpn_divexact_1(dst,src,size,255*CNST_LIMB(188513325))