tests/mpn/t-hgcd.c

86 	  fprintf (stderr, "Bad return code from hgcd\n");
90 	  fprintf (stderr, "hgcd:     %d\n", (int) res);
184   struct hgcd_matrix hgcd;
209   mpn_hgcd_matrix_init (&hgcd, asize, hgcd_init_tp);
241 		     &hgcd, hgcd_tp);
246       fprintf (stderr, "Different return value from hgcd and hgcd_ref\n");
255       if (!hgcd_ref_equal (&hgcd, &ref)
279 hgcd_ref_init (struct hgcd_ref *hgcd)
286 	mpz_init (hgcd->m[i][j]);
291 hgcd_ref_clear (struct hgcd_ref *hgcd)
298 	mpz_clear (hgcd->m[i][j]);
317 hgcd_ref (struct hgcd_ref *hgcd, mpz_t a, mpz_t b)
336       mpz_set_ui (hgcd->m[0][0], 1); mpz_set_ui (hgcd->m[0][1], 0);
337       mpz_set_ui (hgcd->m[1][0], 1); mpz_set_ui (hgcd->m[1][1], 1);
345       mpz_set_ui (hgcd->m[0][0], 1); mpz_set_ui (hgcd->m[0][1], 1);
346       mpz_set_ui (hgcd->m[1][0], 0); mpz_set_ui (hgcd->m[1][1], 1);
362 	  mpz_addmul (hgcd->m[0][1], q, hgcd->m[0][0]);
363 	  mpz_addmul (hgcd->m[1][1], q, hgcd->m[1][0]);
369 	  mpz_addmul (hgcd->m[0][0], q, hgcd->m[0][1]);
370 	  mpz_addmul (hgcd->m[1][0], q, hgcd->m[1][1]);
392 hgcd_ref_equal (const struct hgcd_matrix *hgcd, const struct hgcd_ref *ref)
401 	if (!mpz_mpn_equal (ref->m[i][j], hgcd->p[i][j], hgcd->n))