/* Copyright 1996-2002, 2004, 2006-2008 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 #include #include #include "gmp.h" #include "gmp-impl.h" #include "longlong.h" #include "tests/tests.h" #ifdef OPERATION_mul_1 #define func __gmpn_mul_1 #define reffunc refmpn_mul_1 #define funcname "mpn_mul_1" #endif #ifdef OPERATION_addmul_1 #define func __gmpn_addmul_1 #define reffunc refmpn_addmul_1 #define funcname "mpn_addmul_1" #endif #ifdef OPERATION_submul_1 #define func __gmpn_submul_1 #define reffunc refmpn_submul_1 #define funcname "mpn_submul_1" #endif #if defined (USG) || defined (__SVR4) || defined (_UNICOS) || defined (__hpux) #include int cputime () { if (CLOCKS_PER_SEC < 100000) return clock () * 1000 / CLOCKS_PER_SEC; return clock () / (CLOCKS_PER_SEC / 1000); } #else #include #include #include 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, ref, rp; mp_limb_t cy_ref, cy_try; int i; long t0, t; unsigned int test; mp_limb_t xlimb; mp_size_t size; double cyc; unsigned int ntests; s1 = malloc (SIZE * sizeof (mp_limb_t)); ref = malloc (SIZE * sizeof (mp_limb_t)); rp = malloc ((SIZE + 2) * sizeof (mp_limb_t)); rp++; 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 % (1 + 0x80000 / (SIZE + 20)) == 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 rp[-1] = 0x87654321; rp[size] = 0x12345678; #ifdef FIXED_XLIMB xlimb = FIXED_XLIMB; #else MPN_RANDOM (&xlimb, 1); #endif #if TIMES != 1 mpn_random (s1, size); mpn_random (rp, size); MPN_COPY (ref, rp, size); t0 = cputime(); for (i = 0; i < TIMES; i++) func (ref, s1, size, xlimb); t = cputime() - t0; cyc = ((double) t * CLOCK) / (TIMES * size * 1000.0); printf (funcname ": %5ldms (%.3f cycles/limb) [%.2f Gb/s]\n", t, cyc, CLOCK/cyc*GMP_LIMB_BITS*GMP_LIMB_BITS/1e9); #endif #ifndef NOCHECK MPN_RANDOM (s1, size); #ifdef ZERO memset (rp, 0, size * sizeof *rp); #else MPN_RANDOM (rp, size); #endif #if defined (PRINT) || defined (XPRINT) printf ("xlimb="); mpn_print (&xlimb, 1); #endif #ifdef PRINT #ifndef OPERATION_mul_1 printf ("%*s ", (int) (2 * sizeof(mp_limb_t)), ""); mpn_print (rp, size); #endif printf ("%*s ", (int) (2 * sizeof(mp_limb_t)), ""); mpn_print (s1, size); #endif MPN_COPY (ref, rp, size); cy_ref = reffunc (ref, s1, size, xlimb); cy_try = func (rp, s1, size, xlimb); #ifdef PRINT mpn_print (&cy_ref, 1); mpn_print (ref, size); mpn_print (&cy_try, 1); mpn_print (rp, size); #endif if (cy_ref != cy_try || mpn_cmp (ref, rp, size) != 0 || rp[-1] != 0x87654321 || rp[size] != 0x12345678) { printf ("\n ref%*s try%*s diff\n", LXW - 3, "", 2 * LXW - 6, ""); for (i = 0; i < size; i++) { printf ("%6d: ", i); printf ("%0*llX ", LXW, (unsigned long long) ref[i]); printf ("%0*llX ", LXW, (unsigned long long) rp[i]); print_posneg (rp[i] - ref[i]); printf ("\n"); } printf ("retval: "); printf ("%0*llX ", LXW, (unsigned long long) cy_ref); printf ("%0*llX ", LXW, (unsigned long long) cy_try); print_posneg (cy_try - cy_ref); printf ("\n"); if (rp[-1] != 0x87654321) printf ("clobbered at low end\n"); if (rp[size] != 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 (""); }