summaryrefslogtreecommitdiff
path: root/tests/random.c
diff options
context:
space:
mode:
authorthevenyp <thevenyp@211d60ee-9f03-0410-a15a-8952a2c7a4e4>2009-03-17 17:23:41 +0000
committerthevenyp <thevenyp@211d60ee-9f03-0410-a15a-8952a2c7a4e4>2009-03-17 17:23:41 +0000
commit960f02dc3261c73558e179af2d0be839a75aa885 (patch)
treec4a9e61ccc5c6e4577f7c9a619813ddea86b8f76 /tests/random.c
parent07f353af30b60f3adf1aa55418a1fa185b937b07 (diff)
downloadmpc-960f02dc3261c73558e179af2d0be839a75aa885.tar.gz
Modify #include chain so as to support DLL creation on Cygwin
git-svn-id: svn://scm.gforge.inria.fr/svn/mpc/trunk@457 211d60ee-9f03-0410-a15a-8952a2c7a4e4
Diffstat (limited to 'tests/random.c')
-rw-r--r--tests/random.c331
1 files changed, 164 insertions, 167 deletions
diff --git a/tests/random.c b/tests/random.c
index df0a550..cd16c00 100644
--- a/tests/random.c
+++ b/tests/random.c
@@ -1,167 +1,164 @@
-/* Handle seed for random numbers.
-
-Copyright (C) 2008 Philippe Th\'eveny.
-
-This file is part of the MPC Library.
-
-The MPC 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 2.1 of the License, or (at your
-option) any later version.
-
-The MPC 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 MPC Library; see the file COPYING.LIB. If not, write to
-the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
-MA 02111-1307, USA. */
-
-/* Put test_start at the beginning of your test function and
- test_end at the end.
- These are an adaptation of those of MPFR. */
-
-#include <stdio.h>
-#include <stdlib.h>
-
-#include <gmp.h>
-#include <mpfr.h>
-#include "mpc.h"
-#include "mpc-impl.h"
-
-#include "config.h"
-
-#ifdef TIME_WITH_SYS_TIME
-# include <sys/time.h>
-# include <time.h>
-#else
-# ifdef HAVE_SYS_TIME_H
-# include <sys/time.h>
-# else
-# include <time.h>
-# endif
-#endif
-
-gmp_randstate_t rands;
-char rands_initialized;
-
-void
-test_start (void)
-{
- char *environment_seed;
- unsigned long seed;
-
- if (rands_initialized)
- {
- fprintf (stderr,
- "Put test_start at the beginning of your test function.\n");
- exit (1);
- }
-
- gmp_randinit_default (rands);
- rands_initialized = 1;
-
- environment_seed = getenv ("GMP_CHECK_RANDOMIZE");
- if (environment_seed == NULL)
- gmp_randseed_ui (rands, 0xfac11e);
- else
- {
- seed = atoi (environment_seed);
- if (seed == 0 || seed == 1)
- {
-#if HAVE_GETTIMEOFDAY
- struct timeval tv;
- gettimeofday (&tv, NULL);
- seed = tv.tv_sec + tv.tv_usec;
-#else
- time_t tv;
- time (&tv);
- seed = tv;
-#endif
- gmp_randseed_ui (rands, seed);
- printf ("Seed GMP_CHECK_RANDOMIZE=%lu "
- "(include this in bug reports)\n", seed);
- }
- else
- {
- printf ("Re-seeding with GMP_CHECK_RANDOMIZE=%lu\n", seed);
- gmp_randseed_ui (rands, seed);
- }
- }
-}
-
-void
-test_end (void)
-{
- if (rands_initialized)
- {
- rands_initialized = 0;
- gmp_randclear (rands);
- }
-}
-
-/* Set z to a non zero value random value with absolute values of Re(z) and
- Im(z) either zero (but not both in the same time) or otherwise greater than
- or equal to 2^{emin-1} and less than 2^emax.
- Each part is negative with probability equal to NEGATIVE_PROBABILITY / 256.
- The result has one zero part (but never the two of them) with probability
- equal to ZERO_PROBABILITY / 256.
-*/
-void
-test_default_random (mpc_ptr z, mp_exp_t emin, mp_exp_t emax,
- unsigned int negative_probability,
- unsigned int zero_probability)
-{
- const unsigned long range = (long)emax - (long)emin + 1;
- unsigned long r;
-
- if (!rands_initialized)
- {
- fprintf (stderr,
- "Put test_start at the beginning of your test function.\n");
- exit (1);
- }
-
- do
- {
- mpc_urandom (z, rands);
- } while (mpfr_zero_p (MPC_RE (z)) || mpfr_zero_p (MPC_IM (z)));
-
- if (zero_probability > 256)
- zero_probability = 256;
- r = gmp_urandomb_ui (rands, 19);
- if ((r & 0x1FF) < zero_probability
- || ((r >> 9) & 0x1FF) < zero_probability)
- {
- int zero_re_p = (r & 0x1FF) < zero_probability;
- int zero_im_p = ((r >> 9) & 0x1FF) < zero_probability;
-
- if (zero_re_p && zero_im_p)
- {
- /* we just want one zero part. */
- zero_re_p = (r >> 18) & 1;
- zero_im_p = !zero_re_p;
- }
- if (zero_re_p)
- mpfr_set_ui (MPC_RE (z), 0, GMP_RNDN);
- if (zero_im_p)
- mpfr_set_ui (MPC_IM (z), 0, GMP_RNDN);
- }
- if (!mpfr_zero_p (MPC_RE (z)))
- mpfr_set_exp (MPC_RE (z),
- (mp_exp_t) gmp_urandomm_ui (rands, range) + emin);
-
- if (!mpfr_zero_p (MPC_IM (z)))
- mpfr_set_exp (MPC_IM (z),
- (mp_exp_t) gmp_urandomm_ui (rands, range) + emin);
-
- if (negative_probability > 256)
- negative_probability = 256;
- r = gmp_urandomb_ui (rands, 16);
- if ((r & 0xFF) < negative_probability)
- mpfr_neg (MPC_RE (z), MPC_RE (z), GMP_RNDN);
- if (((r>>8) & 0xFF) < negative_probability)
- mpfr_neg (MPC_IM (z), MPC_IM (z), GMP_RNDN);
-}
+/* Handle seed for random numbers.
+
+Copyright (C) 2008, 2009 Philippe Th\'eveny.
+
+This file is part of the MPC Library.
+
+The MPC 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 2.1 of the License, or (at your
+option) any later version.
+
+The MPC 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 MPC Library; see the file COPYING.LIB. If not, write to
+the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
+MA 02111-1307, USA. */
+
+/* Put test_start at the beginning of your test function and
+ test_end at the end.
+ These are an adaptation of those of MPFR. */
+
+#include <stdio.h>
+#include <stdlib.h>
+
+#include "mpc-tests.h"
+
+#include "config.h"
+
+#ifdef TIME_WITH_SYS_TIME
+# include <sys/time.h>
+# include <time.h>
+#else
+# ifdef HAVE_SYS_TIME_H
+# include <sys/time.h>
+# else
+# include <time.h>
+# endif
+#endif
+
+gmp_randstate_t rands;
+char rands_initialized;
+
+void
+test_start (void)
+{
+ char *environment_seed;
+ unsigned long seed;
+
+ if (rands_initialized)
+ {
+ fprintf (stderr,
+ "Put test_start at the beginning of your test function.\n");
+ exit (1);
+ }
+
+ gmp_randinit_default (rands);
+ rands_initialized = 1;
+
+ environment_seed = getenv ("GMP_CHECK_RANDOMIZE");
+ if (environment_seed == NULL)
+ gmp_randseed_ui (rands, 0xfac11e);
+ else
+ {
+ seed = atoi (environment_seed);
+ if (seed == 0 || seed == 1)
+ {
+#if HAVE_GETTIMEOFDAY
+ struct timeval tv;
+ gettimeofday (&tv, NULL);
+ seed = tv.tv_sec + tv.tv_usec;
+#else
+ time_t tv;
+ time (&tv);
+ seed = tv;
+#endif
+ gmp_randseed_ui (rands, seed);
+ printf ("Seed GMP_CHECK_RANDOMIZE=%lu "
+ "(include this in bug reports)\n", seed);
+ }
+ else
+ {
+ printf ("Re-seeding with GMP_CHECK_RANDOMIZE=%lu\n", seed);
+ gmp_randseed_ui (rands, seed);
+ }
+ }
+}
+
+void
+test_end (void)
+{
+ if (rands_initialized)
+ {
+ rands_initialized = 0;
+ gmp_randclear (rands);
+ }
+}
+
+/* Set z to a non zero value random value with absolute values of Re(z) and
+ Im(z) either zero (but not both in the same time) or otherwise greater than
+ or equal to 2^{emin-1} and less than 2^emax.
+ Each part is negative with probability equal to NEGATIVE_PROBABILITY / 256.
+ The result has one zero part (but never the two of them) with probability
+ equal to ZERO_PROBABILITY / 256.
+*/
+void
+test_default_random (mpc_ptr z, mp_exp_t emin, mp_exp_t emax,
+ unsigned int negative_probability,
+ unsigned int zero_probability)
+{
+ const unsigned long range = (long)emax - (long)emin + 1;
+ unsigned long r;
+
+ if (!rands_initialized)
+ {
+ fprintf (stderr,
+ "Put test_start at the beginning of your test function.\n");
+ exit (1);
+ }
+
+ do
+ {
+ mpc_urandom (z, rands);
+ } while (mpfr_zero_p (MPC_RE (z)) || mpfr_zero_p (MPC_IM (z)));
+
+ if (zero_probability > 256)
+ zero_probability = 256;
+ r = gmp_urandomb_ui (rands, 19);
+ if ((r & 0x1FF) < zero_probability
+ || ((r >> 9) & 0x1FF) < zero_probability)
+ {
+ int zero_re_p = (r & 0x1FF) < zero_probability;
+ int zero_im_p = ((r >> 9) & 0x1FF) < zero_probability;
+
+ if (zero_re_p && zero_im_p)
+ {
+ /* we just want one zero part. */
+ zero_re_p = (r >> 18) & 1;
+ zero_im_p = !zero_re_p;
+ }
+ if (zero_re_p)
+ mpfr_set_ui (MPC_RE (z), 0, GMP_RNDN);
+ if (zero_im_p)
+ mpfr_set_ui (MPC_IM (z), 0, GMP_RNDN);
+ }
+ if (!mpfr_zero_p (MPC_RE (z)))
+ mpfr_set_exp (MPC_RE (z),
+ (mp_exp_t) gmp_urandomm_ui (rands, range) + emin);
+
+ if (!mpfr_zero_p (MPC_IM (z)))
+ mpfr_set_exp (MPC_IM (z),
+ (mp_exp_t) gmp_urandomm_ui (rands, range) + emin);
+
+ if (negative_probability > 256)
+ negative_probability = 256;
+ r = gmp_urandomb_ui (rands, 16);
+ if ((r & 0xFF) < negative_probability)
+ mpfr_neg (MPC_RE (z), MPC_RE (z), GMP_RNDN);
+ if (((r>>8) & 0xFF) < negative_probability)
+ mpfr_neg (MPC_IM (z), MPC_IM (z), GMP_RNDN);
+}