[Gmp-commit] /var/hg/gmp: Provide forgotten file.
mercurial at gmplib.org
mercurial at gmplib.org
Thu Nov 8 18:23:01 UTC 2018
details: /var/hg/gmp/rev/395fb75310ab
changeset: 17683:395fb75310ab
user: Torbjorn Granlund <tg at gmplib.org>
date: Thu Nov 08 19:22:54 2018 +0100
description:
Provide forgotten file.
* tests/devel/cnd_aors_n.c: New file.
diffstat:
tests/devel/cnd_aors_n.c | 257 +++++++++++++++++++++++++++++++++++++++++++++++
1 files changed, 257 insertions(+), 0 deletions(-)
diffs (261 lines):
diff -r 3230453a6131 -r 395fb75310ab tests/devel/cnd_aors_n.c
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/tests/devel/cnd_aors_n.c Thu Nov 08 19:22:54 2018 +0100
@@ -0,0 +1,257 @@
+/*
+Copyright 1996-2004, 2009, 2011 Free Software Foundation, Inc.
+
+This file is part of the GNU MP Library test suite.
+
+The GNU MP Library test suite is free software; you can redistribute it
+and/or modify it under the terms of the GNU 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 test suite 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 General
+Public License for more details.
+
+You should have received a copy of the GNU General Public License along with
+the GNU MP Library test suite. If not, see https://www.gnu.org/licenses/. */
+
+#include <stdlib.h>
+#include <string.h>
+#include <stdio.h>
+#include "gmp-impl.h"
+#include "tests/tests.h"
+
+#ifdef OPERATION_cnd_add_n
+#define func __gmpn_cnd_add_n
+#define reffunc refmpn_cnd_add_n
+#define funcname "mpn_cnd_add_n"
+#endif
+
+#ifdef OPERATION_cnd_sub_n
+#define func __gmpn_cnd_sub_n
+#define reffunc refmpn_cnd_sub_n
+#define funcname "mpn_cnd_sub_n"
+#endif
+
+#if defined (USG) || defined (__SVR4) || defined (_UNICOS) || defined (__hpux)
+#include <time.h>
+
+int
+cputime ()
+{
+ if (CLOCKS_PER_SEC < 100000)
+ return clock () * 1000 / CLOCKS_PER_SEC;
+ return clock () / (CLOCKS_PER_SEC / 1000);
+}
+#else
+#include <sys/types.h>
+#include <sys/time.h>
+#include <sys/resource.h>
+
+int
+cputime ()
+{
+ struct rusage rus;
+
+ getrusage (0, &rus);
+ return rus.ru_utime.tv_sec * 1000 + rus.ru_utime.tv_usec / 1000;
+}
+#endif
+
+static void print_posneg (mp_limb_t);
+static void mpn_print (mp_ptr, mp_size_t);
+
+#define LXW ((int) (2 * sizeof (mp_limb_t)))
+#define M * 1000000
+
+#ifndef CLOCK
+#error "Don't know CLOCK of your machine"
+#endif
+
+#ifndef OPS
+#define OPS (CLOCK/5)
+#endif
+#ifndef SIZE
+#define SIZE 328
+#endif
+#ifndef TIMES
+#define TIMES OPS/(SIZE+1)
+#endif
+
+int
+main (int argc, char **argv)
+{
+ mp_ptr s1, s2, dx, dy;
+ mp_limb_t cyx, cyy;
+ int i;
+#if TIMES != 1
+ long t0, t;
+#endif
+ unsigned int test;
+ mp_size_t size;
+ unsigned int ntests;
+
+ s1 = malloc ((SIZE + 2) * sizeof (mp_limb_t));
+ s2 = malloc ((SIZE + 2) * sizeof (mp_limb_t));
+ dx = malloc ((SIZE + 2) * sizeof (mp_limb_t));
+ dy = malloc ((SIZE + 2) * sizeof (mp_limb_t));
+
+ s1 += 1;
+ s2 += 1;
+
+ ntests = ~(unsigned) 0;
+ if (argc == 2)
+ ntests = strtol (argv[1], 0, 0);
+
+ for (test = 1; test <= ntests; test++)
+ {
+#if TIMES == 1 && ! defined (PRINT)
+ if (test % (SIZE > 100000 ? 1 : 100000 / SIZE) == 0)
+ {
+ printf ("\r%u", test);
+ fflush (stdout);
+ }
+#endif
+
+#ifdef PLAIN_RANDOM
+#define MPN_RANDOM mpn_random
+#else
+#define MPN_RANDOM mpn_random2
+#endif
+
+#ifdef RANDOM
+ size = random () % SIZE + 1;
+#else
+ size = SIZE;
+#endif
+
+ dx[0] = 0x87654321;
+ dy[0] = 0x87654321;
+ dx[size+1] = 0x12345678;
+ dy[size+1] = 0x12345678;
+
+#if TIMES != 1
+ mpn_random (s1 - 1, size + 2);
+ mpn_random (s2 - 1, size + 2);
+
+ t0 = cputime();
+ for (i = 0; i < TIMES; i++)
+ func (i & 1, dx+1, s1, s2, size);
+ t = cputime() - t0;
+ printf (funcname ": %5ldms (%.3f cycles/limb)\n",
+ t, ((double) t * CLOCK) / (TIMES * size * 1000.0));
+#endif
+
+#ifndef NOCHECK
+#ifndef ZEROup
+ MPN_RANDOM (s1 - 1, size + 2);
+#else
+ MPN_ZERO (s1, size);
+#endif
+#ifndef ZEROvp
+ MPN_RANDOM (s2 - 1, size + 2);
+#else
+ MPN_ZERO (s2, size);
+#endif
+
+#ifdef PRINT
+ mpn_print (s1, size);
+ mpn_print (s2, size);
+#endif
+
+ /* Put garbage in the destination. */
+ for (i = 0; i < size; i++)
+ {
+ dx[i+1] = 0xdead;
+ dy[i+1] = 0xbeef;
+ }
+
+ int cond = random() & 1;
+
+ cyx = reffunc (cond, dx+1, s1, s2, size);
+ cyy = func (cond, dy+1, s1, s2, size);
+
+#ifdef PRINT
+ mpn_print (&cyx, 1);
+ mpn_print (dx+1, size);
+ mpn_print (&cyy, 1);
+ mpn_print (dy+1, size);
+#endif
+
+ if (cyx != cyy || mpn_cmp (dx, dy, size+2) != 0
+ || dx[0] != 0x87654321 || dx[size+1] != 0x12345678)
+ {
+ mp_size_t s, e;
+ for (s = 0;; s++)
+ if ((unsigned long long) (dx+1)[s] != (unsigned long long) (dy+1)[s])
+ break;
+ for (e = size - 1;; e--)
+ if ((unsigned long long) (dx+1)[e] != (unsigned long long) (dy+1)[e])
+ break;
+#ifndef PRINT
+ for (i = s; i <= e; i++)
+ {
+ printf ("%6d: ", i);
+ printf ("%0*llX ", LXW, (unsigned long long) (dx+1)[i]);
+ printf ("%0*llX ", LXW, (unsigned long long) (dy+1)[i]);
+ print_posneg ((dy+1)[i] - (dx+1)[i]);
+ printf ("\n");
+ }
+ printf ("%6s: ", "retval");
+ printf ("%0*llX ", LXW, (unsigned long long) cyx);
+ printf ("%0*llX ", LXW, (unsigned long long) cyy);
+ print_posneg (cyx - cyy);
+#endif
+ printf ("\n");
+ if (dy[0] != 0x87654321)
+ printf ("clobbered at low end\n");
+ if (dy[size+1] != 0x12345678)
+ printf ("clobbered at high end\n");
+ printf ("TEST NUMBER %u\n", test);
+ abort();
+ }
+#endif
+ }
+ exit (0);
+}
+
+static void
+print_posneg (mp_limb_t d)
+{
+ char buf[LXW + 2];
+ if (d == 0)
+ printf (" %*X", LXW, 0);
+ else if (-d < d)
+ {
+ sprintf (buf, "%llX", (unsigned long long) -d);
+ printf ("%*s-%s", LXW - (int) strlen (buf), "", buf);
+ }
+ else
+ {
+ sprintf (buf, "%llX", (unsigned long long) d);
+ printf ("%*s+%s", LXW - (int) strlen (buf), "", buf);
+ }
+}
+
+static void
+mpn_print (mp_ptr p, mp_size_t size)
+{
+ mp_size_t i;
+
+ for (i = size - 1; i >= 0; i--)
+ {
+#ifdef _LONG_LONG_LIMB
+ printf ("%0*lX%0*lX", (int) (sizeof(mp_limb_t)),
+ (unsigned long) (p[i] >> (GMP_LIMB_BITS/2)),
+ (int) (sizeof(mp_limb_t)), (unsigned long) (p[i]));
+#else
+ printf ("%0*lX", (int) (2 * sizeof(mp_limb_t)), p[i]);
+#endif
+#ifdef SPACE
+ if (i != 0)
+ printf (" ");
+#endif
+ }
+ puts ("");
+}
More information about the gmp-commit
mailing list