Lines Matching refs:mn
40 mp_size_t mn, fn, size, max_mn;
67 mn = gmp_urandomm_ui (rands, MAX_K) % (fn / 4 + 10);
68 } while (mn == 0);
69 ASSERT (mn <= max_mn);
70 mpn_random2 (mp, mn);
71 ASSERT (mp [mn - 1] != 0);
73 if (fn >= mn)
75 mpn_tdiv_qr (qp, fk, 0, fk, fn, mp, mn);
76 mpn_tdiv_qr (qp, fks1, 0, fks1, fn, mp, mn);
80 MPN_ZERO (fk + fn, mn - fn);
81 MPN_ZERO (fks1 + fn, mn - fn);
84 mpn_random2 (fkm - 1, 2*mn+1+2);
86 fk_after = fkm [2 * mn + 1];
88 mpn_random2 (fks1m - 1, 2*mn+1+2);
90 fk1_after = fks1m [2 * mn + 1];
93 signflip = mpn_fib2m (fkm, fks1m, qp, 1, mp, mn);
94 if (fkm [-1] != fk_before || fkm [2 * mn + 1] != fk_after
95 || fks1m [-1] != fk1_before || fks1m [2 * mn + 1] != fk1_after)
97 printf ("REDZONE violation in test %d, k = %lu, mn = %u\n",
98 test, k, (unsigned) mn);
104 if (fkm[2 * mn + 1] != fk_after)
106 printf ("after fkm:"); mpn_dump (fkm + 2 * mn + 1, 1);
114 if (fks1m[2 * mn + 1] != fk1_after)
116 printf ("after fks1m:"); mpn_dump (fks1m + 2 * mn + 1, 1);
122 if (mpn_cmp (fkm, fk, mn) != 0)
124 if (mpn_sub_n (fk, mp, fk, mn) || mpn_cmp (fkm, fk, mn) != 0)
126 printf ("ERROR(k) in test %d, k = %lu, mn = %u\n",
127 test, k, (unsigned) mn);
128 mpn_dump (fk, mn);
129 mpn_dump (fkm, mn);
130 mpn_dump (mp, mn);
136 if (mpn_cmp (fks1m, fks1, mn) != 0)
138 if (mpn_sub_n (fks1, mp, fks1, mn) || mpn_cmp (fks1m, fks1, mn) != 0)
140 printf ("ERROR(k-1) in test %d, k = %lu, mn = %u\n",
141 test, k, (unsigned) mn);
142 mpn_dump (fks1, mn);
143 mpn_dump (fks1m, mn);
144 mpn_dump (mp, mn);
150 if (signflip != 0 && ! mpn_zero_p (fks1m, mn) && ! mpn_zero_p (fkm, mn))
154 if (! mpn_lshift (fks1m, fks1m, mn, 1) &&
155 mpn_cmp (mp, fks1m, mn) == 0)
157 if (! mpn_lshift (fkm, fkm, mn, 1) &&
158 mpn_cmp (mp, fkm, mn) == 0)
161 printf ("ERROR(sign) in test %d, k = %lu, mn = %u\n",
162 test, k, (unsigned) mn);
189 mp_size_t mn, kn;
193 mn = gmp_urandomm_ui (rands, MAX_MN - 1) + 1;
194 mpn_random2 (mp, mn);
202 mpn_random2 (fka - 1, 2*mn+1+2);
204 fka_after = fka [2 * mn + 1];
206 mpn_random2 (fks1a - 1, 2*mn+1+2);
208 fk1a_after = fks1a [2 * mn + 1];
210 signflip = mpn_fib2m (fka, fks1a, kp, kn, mp, mn);
211 if (fka [-1] != fka_before || fka [2 * mn + 1] != fka_after
212 || fks1a [-1] != fk1a_before || fks1a [2 * mn + 1] != fk1a_after)
220 if (fka[2 * mn + 1] != fka_after)
222 printf ("after fka:"); mpn_dump (fka + 2 * mn + 1, 1);
230 if (fks1a[2 * mn + 1] != fk1a_after)
232 printf ("after fks1a:"); mpn_dump (fks1a + 2 * mn + 1, 1);
238 if (signflip && ! mpn_zero_p (fks1a, mn))
239 mpn_sub_n (fks1a, mp, fks1a, mn);
240 if (mpn_sub_n (fka, fka, fks1a, mn))
241 ASSERT_CARRY (mpn_add_n (fka, fka, mp, mn));
247 mpn_random2 (fkb - 1, 2*mn+1+2);
249 fkb_after = fkb [2 * mn + 1];
251 mpn_random2 (fks1b - 1, 2*mn+1+2);
253 fk1b_after = fks1b [2 * mn + 1];
255 signflip = mpn_fib2m (fkb, fks1b, kp, kn, mp, mn);
256 if (fkb [-1] != fkb_before || fkb [2 * mn + 1] != fkb_after
257 || fks1b [-1] != fk1b_before || fks1b [2 * mn + 1] != fk1b_after)
265 if (fkb[2 * mn + 1] != fkb_after)
267 printf ("after fkb:"); mpn_dump (fkb + 2 * mn + 1, 1);
275 if (fks1b[2 * mn + 1] != fk1b_after)
277 printf ("after fks1b:"); mpn_dump (fks1b + 2 * mn + 1, 1);
283 if (mpn_cmp (fks1a, fkb, mn) != 0)
285 if (mpn_sub_n (fkb, mp, fkb, mn) || mpn_cmp (fks1a, fkb, mn) != 0)
288 mpn_dump (fks1a, mn);
289 mpn_dump (fkb, mn);
290 mpn_dump (mp, mn);
296 if (mpn_cmp (fka, fks1b, mn) != 0)
298 if (mpn_sub_n (fks1b, mp, fks1b, mn) || mpn_cmp (fka, fks1b, mn) != 0)
301 mpn_dump (fka, mn);
302 mpn_dump (fks1b, mn);
303 mpn_dump (mp, mn);
309 if (signflip != 0 && ! mpn_zero_p (fks1b, mn) && ! mpn_zero_p (fkb, mn))
313 if (! mpn_lshift (fks1b, fks1b, mn, 1) &&
314 mpn_cmp (mp, fks1b, mn) == 0)
316 if (! mpn_lshift (fkb, fkb, mn, 1) &&
317 mpn_cmp (mp, fkb, mn) == 0)
Indexes created Sat Mar 28 00:23:22 UTC 2026