113  *   a <- (a*pa+b*pb)/(2^31)
121  * Factors pa, pb, qa and qb must be at most 2^31 in absolute value.
127 	int64_t pa, int64_t pb, int64_t qa, int64_t qb)  in co_reduce()  argument
143 		 *   |pb| <= 2^31  in co_reduce()
155 		za = wa * (uint64_t)pa + wb * (uint64_t)pb + (uint64_t)cca;  in co_reduce()
189  *   a <- (a*pa+b*pb)/(2^31) mod m
194  * Factors pa, pb, qa and qb must be at most 2^31 in absolute value.
200 	int64_t pa, int64_t pb, int64_t qa, int64_t qb,  in co_reduce_mod()  argument
209 	fa = ((a[0] * (uint32_t)pa + b[0] * (uint32_t)pb) * m0i) & 0x7FFFFFFF;  in co_reduce_mod()
222 		za = wa * (uint64_t)pa + wb * (uint64_t)pb  in co_reduce_mod()
297 	 * parameters pa, pb, qa and qb such that the new values for a  in br_i31_moddiv()
300 	 *   a' = (a*pa + b*pb) / (2^31)  in br_i31_moddiv()
307 	 * pa with -pa, and pb with -pb. The total length of a and b is  in br_i31_moddiv()
340 		int64_t pa, pb, qa, qb;  in br_i31_moddiv()  local
388 		 *   a' = a*pa + b*pb  in br_i31_moddiv()
395 		pb = 0;  in br_i31_moddiv()
445 			pb -= qb & -(int64_t)cAB;  in br_i31_moddiv()
449 			qb -= pb & -(int64_t)cBA;  in br_i31_moddiv()
456 			pb += pb & ((int64_t)cA - 1);  in br_i31_moddiv()
467 		r = co_reduce(a, b, len, pa, pb, qa, qb);  in br_i31_moddiv()
469 		pb -= pb * ((r & 1) << 1);  in br_i31_moddiv()
472 		co_reduce_mod(u, v, len, pa, pb, qa, qb, m + 1, m0i);  in br_i31_moddiv()