1 /* Test that routines allow reusing a source variable as destination. */
13 typedef void (*dss_func) (mpz_ptr, mpz_srcptr, mpz_srcptr);
15 typedef void (*dss_func) ();
18 dss_func dss_funcs[] =
20 mpz_add, mpz_and, mpz_cdiv_q, mpz_cdiv_r, mpz_fdiv_q, mpz_fdiv_r,
21 mpz_gcd, mpz_ior, mpz_mul, mpz_sub, mpz_tdiv_q, mpz_tdiv_r
24 char *dss_func_names[] =
26 "mpz_add", "mpz_and", "mpz_cdiv_q", "mpz_cdiv_r", "mpz_fdiv_q", "mpz_fdiv_r",
27 "mpz_gcd", "mpz_ior", "mpz_mul", "mpz_sub", "mpz_tdiv_q", "mpz_tdiv_r"
30 char dss_func_division[] = {0, 0, 1, 1, 1, 1, 0, 0, 0, 0, 1, 1};
33 mpz_divexact /* requires special operands */
41 int pass, reps = 10000;
43 mpz_t res1, res2, res3;
46 reps = atoi (argv[1]);
55 for (pass = 1; pass <= reps; pass++)
57 mpz_random (in1, urandom () % SIZE - SIZE/2);
58 mpz_random (in2, urandom () % SIZE - SIZE/2);
60 for (i = 0; i < sizeof (dss_funcs) / sizeof (dss_func); i++)
62 if (dss_func_division[i] && mpz_cmp_ui (in2, 0) == 0)
65 (dss_funcs[i]) (res1, in1, in2);
68 (dss_funcs[i]) (out1, out1, in2);
72 (dss_funcs[i]) (out1, in1, out1);
75 if (mpz_cmp (res1, res2) != 0)
76 dump_abort (dss_func_names[i], in1, in2);
77 if (mpz_cmp (res1, res3) != 0)
78 dump_abort (dss_func_names[i], in1, in2);
85 dump_abort (name, in1, in2)
89 printf ("failure in %s (", name);
90 mpz_out_str (stdout, -16, in1);
92 mpz_out_str (stdout, -16, in2);
98 void mpz_add_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
99 void mpz_div_2exp _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
100 void mpz_mod_2exp _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
101 void mpz_mul_2exp _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
102 void mpz_mul_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
103 void mpz_pow_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
104 void mpz_sub_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
105 void mpz_tdiv_q_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
106 void mpz_tdiv_r_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
108 void mpz_abs _PROTO ((mpz_ptr, mpz_srcptr));
109 void mpz_com _PROTO ((mpz_ptr, mpz_srcptr));
110 void mpz_sqrt _PROTO ((mpz_ptr, mpz_srcptr));
111 void mpz_neg _PROTO ((mpz_ptr, mpz_srcptr));
113 void mpz_tdiv_qr_ui _PROTO ((mpz_ptr, mpz_ptr, mpz_srcptr, unsigned long int));
115 void mpz_powm_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int, mpz_srcptr));
117 void mpz_gcdext _PROTO ((mpz_ptr, mpz_ptr, mpz_ptr, mpz_srcptr, mpz_srcptr));
119 void mpz_cdiv_qr _PROTO ((mpz_ptr, mpz_ptr, mpz_srcptr, mpz_srcptr));
120 void mpz_fdiv_qr _PROTO ((mpz_ptr, mpz_ptr, mpz_srcptr, mpz_srcptr));
121 void mpz_tdiv_qr _PROTO ((mpz_ptr, mpz_ptr, mpz_srcptr, mpz_srcptr));
123 void mpz_powm _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr, mpz_srcptr));
125 void mpz_sqrtrem _PROTO ((mpz_ptr, mpz_ptr, mpz_srcptr));
127 unsigned long int mpz_cdiv_qr_ui _PROTO ((mpz_ptr, mpz_ptr, mpz_srcptr, unsigned long int));
128 unsigned long int mpz_fdiv_qr_ui _PROTO ((mpz_ptr, mpz_ptr, mpz_srcptr, unsigned long int));
130 unsigned long int mpz_cdiv_q_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
131 unsigned long int mpz_cdiv_r_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
132 unsigned long int mpz_fdiv_q_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
133 unsigned long int mpz_fdiv_r_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
134 unsigned long int mpz_gcd_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));