Return to reuse.c CVS log

Up to [local] / OpenXM_contrib / gmp / mpz / tests

Annotation of OpenXM_contrib/gmp/mpz/tests/reuse.c, Revision 1.1.1.1

1.1       maekawa     1: /* Test that routines allow reusing a source variable as destination.  */
                      2:
                      3: #include <stdio.h>
                      4: #include "gmp.h"
                      5: #include "gmp-impl.h"
                      6: #include "urandom.h"
                      7:
                      8: #ifndef SIZE
                      9: #define SIZE 16
                     10: #endif
                     11:
                     12: #if __STDC__
                     13: typedef void (*dss_func) (mpz_ptr, mpz_srcptr, mpz_srcptr);
                     14: #else
                     15: typedef void (*dss_func) ();
                     16: #endif
                     17:
                     18: dss_func dss_funcs[] =
                     19: {
                     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
                     22: };
                     23:
                     24: char *dss_func_names[] =
                     25: {
                     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"
                     28: };
                     29:
                     30: char dss_func_division[] = {0, 0, 1, 1, 1, 1, 0, 0, 0, 0, 1, 1};
                     31:
                     32: #if 0
                     33: mpz_divexact /* requires special operands */
                     34: #endif
                     35:
                     36: main (argc, argv)
                     37:      int argc;
                     38:      char **argv;
                     39: {
                     40:   int i;
                     41:   int pass, reps = 10000;
                     42:   mpz_t in1, in2, out1;
                     43:   mpz_t res1, res2, res3;
                     44:
                     45:   if (argc == 2)
                     46:      reps = atoi (argv[1]);
                     47:
                     48:   mpz_init (in1);
                     49:   mpz_init (in2);
                     50:   mpz_init (out1);
                     51:   mpz_init (res1);
                     52:   mpz_init (res2);
                     53:   mpz_init (res3);
                     54:
                     55:   for (pass = 1; pass <= reps; pass++)
                     56:     {
                     57:       mpz_random (in1, urandom () % SIZE - SIZE/2);
                     58:       mpz_random (in2, urandom () % SIZE - SIZE/2);
                     59:
                     60:       for (i = 0; i < sizeof (dss_funcs) / sizeof (dss_func); i++)
                     61:        {
                     62:          if (dss_func_division[i] && mpz_cmp_ui (in2, 0) == 0)
                     63:            continue;
                     64:
                     65:          (dss_funcs[i]) (res1, in1, in2);
                     66:
                     67:          mpz_set (out1, in1);
                     68:          (dss_funcs[i]) (out1, out1, in2);
                     69:          mpz_set (res2, out1);
                     70:
                     71:          mpz_set (out1, in2);
                     72:          (dss_funcs[i]) (out1, in1, out1);
                     73:          mpz_set (res3, out1);
                     74:
                     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);
                     79:        }
                     80:     }
                     81:
                     82:   exit (0);
                     83: }
                     84:
                     85: dump_abort (name, in1, in2)
                     86:      char *name;
                     87:      mpz_t in1, in2;
                     88: {
                     89:   printf ("failure in %s (", name);
                     90:   mpz_out_str (stdout, -16, in1);
                     91:   printf (" ");
                     92:   mpz_out_str (stdout, -16, in2);
                     93:   printf (")\n");
                     94:   abort ();
                     95: }
                     96:
                     97: #if 0
                     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));
                    107:
                    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));
                    112:
                    113: void mpz_tdiv_qr_ui    _PROTO ((mpz_ptr, mpz_ptr, mpz_srcptr, unsigned long int));
                    114:
                    115: void mpz_powm_ui       _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int, mpz_srcptr));
                    116:
                    117: void mpz_gcdext                _PROTO ((mpz_ptr, mpz_ptr, mpz_ptr, mpz_srcptr, mpz_srcptr));
                    118:
                    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));
                    122:
                    123: void mpz_powm          _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr, mpz_srcptr));
                    124:
                    125: void mpz_sqrtrem       _PROTO ((mpz_ptr, mpz_ptr, mpz_srcptr));
                    126:
                    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));
                    129:
                    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));
                    135: #endif