diff options
Diffstat (limited to 'gmp/tests/devel/copy.c')
| -rw-r--r-- | gmp/tests/devel/copy.c | 226 |
1 files changed, 226 insertions, 0 deletions
diff --git a/gmp/tests/devel/copy.c b/gmp/tests/devel/copy.c new file mode 100644 index 0000000000..1f951f5389 --- /dev/null +++ b/gmp/tests/devel/copy.c @@ -0,0 +1,226 @@ +/* +Copyright 1999-2001, 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.h" +#include "gmp-impl.h" +#include "tests.h" + +#ifdef OPERATION_copyi +#define func MPN_COPY_INCR +#define reffunc refmpn_copyi +#define funcname "MPN_COPY_INCR" +#endif + +#ifdef OPERATION_copyd +#define func MPN_COPY_DECR +#define reffunc refmpn_copyd +#define funcname "MPN_COPY_DECR" +#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 496 +#endif +#ifndef TIMES +#define TIMES OPS/(SIZE+1) +#endif + +int +main (int argc, char **argv) +{ + mp_ptr s1, dx, dy; + int i; + long t0, t; + unsigned int test; + mp_size_t size; + unsigned int ntests; + + s1 = malloc (SIZE * sizeof (mp_limb_t)); + dx = malloc ((SIZE + 2) * sizeof (mp_limb_t)); + dy = malloc ((SIZE + 2) * sizeof (mp_limb_t)); + + 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 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, size); + + t0 = cputime(); + for (i = 0; i < TIMES; i++) + func (dx+1, s1, size); + t = cputime() - t0; + printf (funcname ": %5ldms (%.3f cycles/limb)\n", + t, ((double) t * CLOCK) / (TIMES * size * 1000.0)); +#endif + +#ifndef NOCHECK + mpn_random2 (s1, size); + +#ifdef PRINT + mpn_print (s1, size); +#endif + + /* Put garbage in the destination. */ + for (i = 0; i < size; i++) + { + dx[i+1] = 0xdead; + dy[i+1] = 0xbeef; + } + + reffunc (dx+1, s1, size); + func (dy+1, s1, size); + +#ifdef PRINT + mpn_print (dx+1, size); + mpn_print (dy+1, size); +#endif + + if (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"); + } +#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 (""); +} |