tests/devel/aors_n.c

   1 /*
   2 Copyright 1996, 1997, 1998, 1999, 2000, 2001, 2002, 2003, 2004, 2009
   3 Free Software Foundation, Inc.
   4
   5 This file is part of the GNU MP Library.
   6
   7 The GNU MP Library is free software; you can redistribute it and/or modify
   8 it under the terms of the GNU Lesser General Public License as published by
   9 the Free Software Foundation; either version 3 of the License, or (at your
  10 option) any later version.
  11
  12 The GNU MP Library is distributed in the hope that it will be useful, but
  13 WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
  14 or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public
  15 License for more details.
  16
  17 You should have received a copy of the GNU Lesser General Public License
  18 along with the GNU MP Library.  If not, see http://www.gnu.org/licenses/.  */
  19
  20 #include <stdlib.h>
  21 #include <stdio.h>
  22 #include "gmp.h"
  23 #include "gmp-impl.h"
  24 #include "tests.h"
  25
  26 #ifdef OPERATION_add_n
  27 #define func __gmpn_add_n
  28 #define reffunc refmpn_add_n
  29 #define funcname "mpn_add_n"
  30 #endif
  31
  32 #ifdef OPERATION_sub_n
  33 #define func __gmpn_sub_n
  34 #define reffunc refmpn_sub_n
  35 #define funcname "mpn_sub_n"
  36 #endif
  37
  38 #ifdef OPERATION_addlsh1_n
  39 #define func __gmpn_addlsh1_n
  40 #define reffunc refmpn_addlsh1_n
  41 #define funcname "mpn_addlsh1_n"
  42 #endif
  43
  44 #ifdef OPERATION_sublsh1_n
  45 #define func __gmpn_sublsh1_n
  46 #define reffunc refmpn_sublsh1_n
  47 #define funcname "mpn_sublsh1_n"
  48 #endif
  49
  50 #ifdef OPERATION_rsh1add_n
  51 #define func __gmpn_rsh1add_n
  52 #define reffunc refmpn_rsh1add_n
  53 #define funcname "mpn_rsh1add_n"
  54 #endif
  55
  56 #ifdef OPERATION_rsh1sub_n
  57 #define func __gmpn_rsh1sub_n
  58 #define reffunc refmpn_rsh1sub_n
  59 #define funcname "mpn_rsh1sub_n"
  60 #endif
  61
  62 #if defined (USG) || defined (__SVR4) || defined (_UNICOS) || defined (__hpux)
  63 #include <time.h>
  64
  65 int
  66 cputime ()
  67 {
  68   if (CLOCKS_PER_SEC < 100000)
  69     return clock () * 1000 / CLOCKS_PER_SEC;
  70   return clock () / (CLOCKS_PER_SEC / 1000);
  71 }
  72 #else
  73 #include <sys/types.h>
  74 #include <sys/time.h>
  75 #include <sys/resource.h>
  76
  77 int
  78 cputime ()
  79 {
  80   struct rusage rus;
  81
  82   getrusage (0, &rus);
  83   return rus.ru_utime.tv_sec * 1000 + rus.ru_utime.tv_usec / 1000;
  84 }
  85 #endif
  86
  87 static void mpn_print (mp_ptr, mp_size_t);
  88
  89 #define M * 1000000
  90
  91 #ifndef CLOCK
  92 #error "Don't know CLOCK of your machine"
  93 #endif
  94
  95 #ifndef OPS
  96 #define OPS (CLOCK/5)
  97 #endif
  98 #ifndef SIZE
  99 #define SIZE 328
 100 #endif
 101 #ifndef TIMES
 102 #define TIMES OPS/(SIZE+1)
 103 #endif
 104
 105 int
 106 main (int argc, char **argv)
 107 {
 108   mp_ptr s1, s2, dx, dy;
 109   mp_limb_t cyx, cyy;
 110   int i;
 111 #if TIMES != 1
 112   long t0, t;
 113 #endif
 114   unsigned int test;
 115   mp_size_t size;
 116   unsigned int ntests;
 117
 118   s1 = malloc (SIZE * sizeof (mp_limb_t));
 119   s2 = malloc (SIZE * sizeof (mp_limb_t));
 120   dx = malloc ((SIZE + 2) * sizeof (mp_limb_t));
 121   dy = malloc ((SIZE + 2) * sizeof (mp_limb_t));
 122
 123   ntests = ~(unsigned) 0;
 124   if (argc == 2)
 125     ntests = strtol (argv[1], 0, 0);
 126
 127   for (test = 1; test <= ntests; test++)
 128     {
 129 #if TIMES == 1 && ! defined (PRINT)
 130       if (test % (SIZE > 100000 ? 1 : 100000 / SIZE) == 0)
 131         {
 132           printf ("\r%u", test);
 133           fflush (stdout);
 134         }
 135 #endif
 136
 137 #ifdef RANDOM
 138       size = random () % SIZE + 1;
 139 #else
 140       size = SIZE;
 141 #endif
 142
 143       dx[0] = 0x87654321;
 144       dy[0] = 0x87654321;
 145       dx[size+1] = 0x12345678;
 146       dy[size+1] = 0x12345678;
 147
 148 #if TIMES != 1
 149       mpn_random (s1, size);
 150       mpn_random (s2, size);
 151
 152       t0 = cputime();
 153       for (i = 0; i < TIMES; i++)
 154         func (dx+1, s1, s2, size);
 155       t = cputime() - t0;
 156       printf (funcname ":    %5ldms (%.3f cycles/limb)\n",
 157               t, ((double) t * CLOCK) / (TIMES * size * 1000.0));
 158 #endif
 159
 160 #ifndef NOCHECK
 161       mpn_random2 (s1, size);
 162       mpn_random2 (s2, size);
 163
 164 #ifdef PRINT
 165       mpn_print (s1, size);
 166       mpn_print (s2, size);
 167 #endif
 168
 169       /* Put garbage in the destination.  */
 170       for (i = 0; i < size; i++)
 171         {
 172           dx[i+1] = 0xdead;
 173           dy[i+1] = 0xbeef;
 174         }
 175
 176       cyx = reffunc (dx+1, s1, s2, size);
 177       cyy = func (dy+1, s1, s2, size);
 178
 179 #ifdef PRINT
 180       mpn_print (&cyx, 1);
 181       mpn_print (dx+1, size);
 182       mpn_print (&cyy, 1);
 183       mpn_print (dy+1, size);
 184 #endif
 185
 186       if (cyx != cyy || mpn_cmp (dx, dy, size+2) != 0
 187           || dx[0] != 0x87654321 || dx[size+1] != 0x12345678)
 188         {
 189 #ifndef PRINT
 190           mpn_print (&cyx, 1);
 191           mpn_print (dx+1, size);
 192           mpn_print (&cyy, 1);
 193           mpn_print (dy+1, size);
 194 #endif
 195           printf ("\n");
 196           if (dy[0] != 0x87654321)
 197             printf ("clobbered at low end\n");
 198           if (dy[size+1] != 0x12345678)
 199             printf ("clobbered at high end\n");
 200           printf ("TEST NUMBER %u\n", test);
 201           abort();
 202         }
 203 #endif
 204     }
 205   exit (0);
 206 }
 207
 208 static void
 209 mpn_print (mp_ptr p, mp_size_t size)
 210 {
 211   mp_size_t i;
 212
 213   for (i = size - 1; i >= 0; i--)
 214     {
 215 #ifdef _LONG_LONG_LIMB
 216       printf ("%0*lX%0*lX", (int) (sizeof(mp_limb_t)),
 217               (unsigned long) (p[i] >> (GMP_LIMB_BITS/2)),
 218               (int) (sizeof(mp_limb_t)), (unsigned long) (p[i]));
 219 #else
 220       printf ("%0*lX", (int) (2 * sizeof(mp_limb_t)), p[i]);
 221 #endif
 222 #ifdef SPACE
 223       if (i != 0)
 224         printf (" ");
 225 #endif
 226     }
 227   puts ("");
 228 }