tests/mpz/reuse.c

212   mpz_t in1, in2, in3;
229   mpz_init (in1);
266       MAKE_RANDOM_OP (in1, size_range, 0);
273 	     mpz_sizeinbase(in1,2),
281 	mpz_neg (in1, in1);
294 	  (dss[i].fptr) (ref1, in1, in2);
297 	  mpz_set (res1, in1);
301 	    FAIL (dss, i, in1, in2, NULL);
304 	  INVOKE_RSS (dss[i], res1, in1, res1);
307 	    FAIL (dss, i, in1, in2, NULL);
315 	  (ddss_div[i].fptr) (ref1, ref2, in1, in2);
319 	  mpz_set (res1, in1);
325 	    FAIL (ddss_div, i, in1, in2, NULL);
328 	  mpz_set (res2, in1);
333 	    FAIL (ddss_div, i, in1, in2, NULL);
337 	  INVOKE_RRSS (ddss_div[i], res1, res2, in1, res1);
341 	    FAIL (ddss_div, i, in1, in2, NULL);
345 	  INVOKE_RRSS (ddss_div[i], res1, res2, in1, res2);
349 	    FAIL (ddss_div, i, in1, in2, NULL);
354 	  if (ds[i].nonneg && mpz_sgn (in1) < 0)
357 	  (ds[i].fptr) (ref1, in1);
360 	  mpz_set (res1, in1);
364 	    FAIL (ds, i, in1, in2, NULL);
374 	  (dsi[i].fptr) (ref1, in1, in2i);
377 	  mpz_set (res1, in1);
381 	    FAIL (dsi, i, in1, in2, NULL);
388 	      r1 = (dsi_div[i].fptr) (ref1, in1, in2i);
391 	      mpz_set (res1, in1);
395 		FAIL (dsi_div, i, in1, in2, NULL);
400 	      r1 = (ddsi_div[i].fptr) (ref1, ref2, in1, in2i);
403 	      mpz_set (res1, in1);
408 		FAIL (ddsi_div, i, in1, in2, NULL);
411 	      mpz_set (res2, in1);
415 		FAIL (ddsi_div, i, in1, in2, NULL);
419       if (mpz_sgn (in1) >= 0)
421 	  mpz_sqrtrem (ref1, ref2, in1);
425 	  mpz_set (res1, in1);
430 	    FAIL2 (mpz_sqrtrem, in1, NULL, NULL);
432 	  mpz_set (res2, in1);
437 	    FAIL2 (mpz_sqrtrem, in1, NULL, NULL);
439 	  mpz_set (res1, in1);
443 	    FAIL2 (mpz_sqrtrem, in1, NULL, NULL);
446       if (mpz_sgn (in1) >= 0)
448 	  mpz_root (ref1, in1, in2i % 0x100 + 1);
451 	  mpz_set (res1, in1);
455 	    FAIL2 (mpz_root, in1, in2, NULL);
458       if (mpz_sgn (in1) >= 0)
460 	  mpz_rootrem (ref1, ref2, in1, in2i % 0x100 + 1);
464 	  mpz_set (res1, in1);
469 	    FAIL2 (mpz_rootrem, in1, in2, NULL);
471 	  mpz_set (res2, in1);
476 	    FAIL2 (mpz_rootrem, in1, in2, NULL);
481 	  mpz_gcdext (ref1, ref2, ref3, in1, in2);
503 	  mpz_set (res1, in1);
509 	  mpz_set (res2, in1);
515 	  mpz_set (res3, in1);
521 	  GCDEXT_CHECK3 (in1, res1);
526 	  GCDEXT_CHECK3 (in1, res2);
531 	  GCDEXT_CHECK3 (in1, res3);
533 	  mpz_set (res1, in1);
538 	  mpz_set (res1, in1);
544 	  mpz_set (res2, in1);
549 	  mpz_set (res2, in1);
555 	  mpz_set (res3, in1);
560 	  mpz_set (res3, in1);
563 	  mpz_set (res1, in1);
568 	  mpz_set (res2, in1);
573 	  GCDEXT_CHECK2 (in1, res1);
577 	  GCDEXT_CHECK2 (in1, res2);
579 	  /* Identical inputs, gcd(in1, in1). Then the result should be
580 	     gcd = abs(in1), s = 0, t = sgn(in1). */
581 	  mpz_abs (ref1, in1);
583 	  mpz_set_si (ref3, mpz_sgn (in1));
602 	  mpz_set (res1, in1);
608 	  mpz_set (res2, in1);
614 	  mpz_set (res3, in1);
617 	  mpz_set (res1, in1);
623 	  mpz_set (res2, in1);
631 	  && (mpz_sgn (in2) >= 0 || mpz_invert (t, in1, in3)))
633 	  mpz_powm (ref1, in1, in2, in3);
636 	  mpz_set (res1, in1);
640 	    FAIL2 (mpz_powm, in1, in2, in3);
643 	  mpz_powm (res1, in1, res1, in3);
646 	    FAIL2 (mpz_powm, in1, in2, in3);
649 	  mpz_powm (res1, in1, in2, res1);
652 	    FAIL2 (mpz_powm, in1, in2, in3);
658 	  mpz_powm_ui (ref1, in1, in2i, in3);
661 	  mpz_set (res1, in1);
665 	    FAIL2 (mpz_powm_ui, in1, in2, in3);
668 	  mpz_powm_ui (res1, in1, in2i, res1);
671 	    FAIL2 (mpz_powm_ui, in1, in2, in3);
675 	r1 = mpz_gcd_ui (ref1, in1, in2i);
678 	mpz_set (res1, in1);
682 	  FAIL2 (mpz_gcd_ui, in1, in2, NULL);
689 	  refretval = mpz_remove (ref1, in1, in2);
692 	  mpz_set (res1, in1);
696 	    FAIL2 (mpz_remove, in1, in2, NULL);
699 	  retval = mpz_remove (res1, in1, res1);
702 	    FAIL2 (mpz_remove, in1, in2, NULL);
708 	  mpz_mul (t, in1, in2);
728 	  mpz_mul (t, in1, in2);
750   mpz_clear (in1);
769 dump (const char *name, mpz_t in1, mpz_t in2, mpz_t in3)
772   mpz_out_str (stdout, -16, in1);