| /* |
| Copyright 1996, 1997, 1998, 1999, 2000, 2001, 2002, 2004, 2006, 2007, 2008 |
| 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 <stdlib.h> |
| #include <string.h> |
| #include <stdio.h> |
| #include "gmp.h" |
| #include "gmp-impl.h" |
| #include "longlong.h" |
| #include "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 <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, 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 RANDOM |
| size = random () % SIZE + 1; |
| #else |
| size = SIZE; |
| #endif |
| |
| rp[-1] = 0x87654321; |
| rp[size] = 0x12345678; |
| |
| #ifdef FIXED_XLIMB |
| xlimb = FIXED_XLIMB; |
| #else |
| mpn_random2 (&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*BITS_PER_MP_LIMB*BITS_PER_MP_LIMB/1e9); |
| #endif |
| |
| #ifndef NOCHECK |
| mpn_random2 (s1, size); |
| #ifdef ZERO |
| memset (rp, 0, size * sizeof *rp); |
| #else |
| mpn_random2 (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] >> (BITS_PER_MP_LIMB/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 (""); |
| } |