Imported gcc-4.4.3

[msp430-gcc.git] / gmp / tests / rand / t-urndmm.c

diff --git a/gmp/tests/rand/t-urndmm.c b/gmp/tests/rand/t-urndmm.c
new file mode 100644 (file)
index 0000000..11eeef2
--- /dev/null
+++ b/gmp/tests/rand/t-urndmm.c
@@ -0,0 +1,159 @@
+/* Test mpz_urandomm.
+
+Copyright 2002 Free Software Foundation, Inc.
+
+This file is part of the GNU MP Library.
+
+The GNU MP Library is free software; you can redistribute it and/or modify
+it under the terms of the GNU Lesser General Public License as published by
+the Free Software Foundation; either version 3 of the License, or (at your
+option) any later version.
+
+The GNU MP Library is distributed in the hope that it will be useful, but
+WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
+or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public
+License for more details.
+
+You should have received a copy of the GNU Lesser General Public License
+along with the GNU MP Library.  If not, see http://www.gnu.org/licenses/.  */
+
+#include <stdio.h>
+#include "gmp.h"
+#include "gmp-impl.h"
+#include "tests.h"
+
+#ifndef TRUE
+#define TRUE (1)
+#endif
+#ifndef FALSE
+#define FALSE (0)
+#endif
+
+int
+check_params (void)
+{
+  gmp_randstate_t r1, r2;
+  mpz_t a, b, m;
+  int i;
+  int result;
+
+  result = TRUE;
+
+  mpz_init (a);
+  mpz_init (b);
+  mpz_init (m);
+
+  if (result)
+    {
+      /* Test the consistency between urandomm and urandomb. */
+      gmp_randinit_default (r1);
+      gmp_randinit_default (r2);
+      gmp_randseed_ui (r1, 85L);
+      gmp_randseed_ui (r2, 85L);
+      mpz_set_ui (m, 0L);
+      mpz_setbit (m, 80L);
+      for (i = 0; i < 100; i++)
+       {
+         mpz_urandomm (a, r1, m);
+         mpz_urandomb (b, r2, 80L);
+         if (mpz_cmp (a, b) != 0)
+           {
+             result = FALSE;
+             printf ("mpz_urandomm != mpz_urandomb\n");
+             break;
+           }
+       }
+      gmp_randclear (r1);
+      gmp_randclear (r2);
+    }
+
+  if (result)
+    {
+      /* Test that mpz_urandomm returns the correct result with a
+        broken LC.  */
+      mpz_set_ui (a, 0L);
+      gmp_randinit_lc_2exp (r1, a, 0xffL, 8L);
+      mpz_set_ui (m, 5L);
+      /* Warning: This code hangs in gmp 4.1 and below */
+      for (i = 0; i < 100; i++)
+       {
+         mpz_urandomm (a, r1, m);
+         if (mpz_cmp_ui (a, 2L) != 0)
+           {
+             result = FALSE;
+             gmp_printf ("mpz_urandomm returns %Zd instead of 2\n", a);
+             break;
+           }
+       }
+      gmp_randclear (r1);
+    }
+
+  if (result)
+    {
+      /* Test that the results are always in range for either
+         positive or negative values of m.  */
+      gmp_randinit_default (r1);
+      mpz_set_ui (m, 5L);
+      mpz_set_si (b, -5L);
+      for (i = 0; i < 100; i++)
+       {
+         mpz_urandomm (a, r1, m);
+         if (mpz_cmp_ui (a, 5L) >= 0 || mpz_sgn (a) < 0)
+           {
+             result = FALSE;
+             gmp_printf ("Out-of-range or non-positive value: %Zd\n", a);
+             break;
+           }
+         mpz_urandomm (a, r1, b);
+         if (mpz_cmp_ui (a, 5L) >= 0 || mpz_sgn (a) < 0)
+           {
+             result = FALSE;
+             gmp_printf ("Out-of-range or non-positive value (from negative modulus): %Zd\n", a);
+             break;
+           }
+       }
+      gmp_randclear (r1);
+    }
+
+  if (result)
+    {
+      /* Test that m=1 forces always result=0.  */
+      gmp_randinit_default (r1);
+      mpz_set_ui (m, 1L);
+      for (i = 0; i < 100; i++)
+       {
+         mpz_urandomm (a, r1, m);
+         if (mpz_sgn (a) != 0)
+           {
+             result = FALSE;
+             gmp_printf ("mpz_urandomm fails with m=1 (result=%Zd)\n", a);
+             break;
+           }
+       }
+      gmp_randclear (r1);
+    }
+
+  mpz_clear (a);
+  mpz_clear (b);
+  mpz_clear (m);
+  return result;
+}
+
+int
+main (int argc, char *argv[])
+{
+  int result = TRUE;
+
+  tests_start ();
+
+  if (result)
+    if (!check_params ())
+      result = FALSE;
+
+  tests_end ();
+
+  if (result)
+    return 0; /* pass */
+  else
+    return 1; /* fail */
+}