[BACK]Return to sub_n.c CVS log [TXT][DIR] Up to [local] / OpenXM_contrib / gmp / mpn / tests

Diff for /OpenXM_contrib/gmp/mpn/tests/Attic/sub_n.c between version 1.1.1.1 and 1.1.1.2

version 1.1.1.1, 2000/01/10 15:35:26 version 1.1.1.2, 2000/09/09 14:12:46
Line 1 
Line 1 
   /*
   Copyright (C) 1996, 1997, 1998, 1999, 2000 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 2.1 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; see the file COPYING.LIB.  If not, write to
   the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
   MA 02111-1307, USA.
   */
   
 #include <stdio.h>  #include <stdio.h>
 #include "gmp.h"  #include "gmp.h"
 #include "gmp-impl.h"  #include "gmp-impl.h"
   
 #ifndef USG  #if defined (USG) || defined (__SVR4) || defined (_UNICOS) || defined (__hpux)
 #include <sys/time.h>  #include <time.h>
 #include <sys/resource.h>  
   
 unsigned long  int
 cputime ()  cputime ()
 {  {
     struct rusage rus;    if (CLOCKS_PER_SEC < 100000)
       return clock () * 1000 / CLOCKS_PER_SEC;
     getrusage (0, &rus);    return clock () / (CLOCKS_PER_SEC / 1000);
     return rus.ru_utime.tv_sec * 1000 + rus.ru_utime.tv_usec / 1000;  
 }  }
 #else  #else
 #include <time.h>  #include <sys/types.h>
   #include <sys/time.h>
   #include <sys/resource.h>
   
 #ifndef CLOCKS_PER_SEC  int
 #define CLOCKS_PER_SEC 1000000  
 #endif  
   
 #if CLOCKS_PER_SEC >= 10000  
 #define CLOCK_TO_MILLISEC(cl) ((cl) / (CLOCKS_PER_SEC / 1000))  
 #else  
 #define CLOCK_TO_MILLISEC(cl) ((cl) * 1000 / CLOCKS_PER_SEC)  
 #endif  
   
 unsigned long  
 cputime ()  cputime ()
 {  {
   return CLOCK_TO_MILLISEC (clock ());    struct rusage rus;
   
     getrusage (0, &rus);
     return rus.ru_utime.tv_sec * 1000 + rus.ru_utime.tv_usec / 1000;
 }  }
 #endif  #endif
   
Line 40  cputime ()
Line 54  cputime ()
 #if defined (__m88k__)  #if defined (__m88k__)
 #define CLOCK 20 M  #define CLOCK 20 M
 #elif defined (__i386__)  #elif defined (__i386__)
 #define CLOCK (16.666667 M)  #define CLOCK (16666667)
 #elif defined (__m68k__)  #elif defined (__m68k__)
 #define CLOCK (20 M)  #define CLOCK (20 M)
 #elif defined (_IBMR2)  #elif defined (_IBMR2)
Line 61  cputime ()
Line 75  cputime ()
 #endif  #endif
   
 #ifndef OPS  #ifndef OPS
 #define OPS 10000000  #define OPS (CLOCK/5)
 #endif  #endif
 #ifndef SIZE  #ifndef SIZE
 #define SIZE 328  #define SIZE 328
Line 120  main (argc, argv)
Line 134  main (argc, argv)
 {  {
   mp_limb_t s1[SIZE];    mp_limb_t s1[SIZE];
   mp_limb_t s2[SIZE];    mp_limb_t s2[SIZE];
   mp_limb_t dx[SIZE+1];    mp_limb_t dx[SIZE+2];
   mp_limb_t dy[SIZE+1];    mp_limb_t dy[SIZE+2];
   int cyx, cyy;    int cyx, cyy;
   int i;    int i;
   long t0, t;    long t0, t;
Line 130  main (argc, argv)
Line 144  main (argc, argv)
   
   for (test = 0; ; test++)    for (test = 0; ; test++)
     {      {
   #if TIMES == 1 && ! defined (PRINT)
         if (test % (SIZE > 10000 ? 1 : 10000 / SIZE) == 0)
           {
             printf ("\r%d", test);
             fflush (stdout);
           }
   #endif
   
 #ifdef RANDOM  #ifdef RANDOM
       size = (random () % SIZE + 1);        size = random () % SIZE + 1;
 #else  #else
       size = SIZE;        size = SIZE;
 #endif  #endif
   
       mpn_random2 (s1, size);        dx[0] = 0x87654321;
       mpn_random2 (s2, size);        dy[0] = 0x87654321;
         dx[size+1] = 0x12345678;
         dy[size+1] = 0x12345678;
   
       dx[size] = 0x12345678;  #if TIMES != 1
       dy[size] = 0x12345678;        mpn_random (s1, size);
         mpn_random (s2, size);
   
 #ifdef PRINT  #ifndef NOCHECK
       mpn_print (s1, size);  
       mpn_print (s2, size);  
 #endif  
       t0 = cputime();        t0 = cputime();
       for (i = 0; i < TIMES; i++)        for (i = 0; i < TIMES; i++)
         cyx = refmpn_sub_n (dx, s1, s2, size);          refmpn_sub_n (dx+1, s1, s2, size);
       t = cputime() - t0;        t = cputime() - t0;
 #if TIMES != 1        printf ("refmpn_sub_n: %5ldms (%.2f cycles/limb)\n",
       printf ("refmpn_sub_n:   %ldms (%.2f cycles/limb)\n",                t, ((double) t * CLOCK) / (OPS * 1000.0));
               t,  
               ((double) t * CLOCK) / (OPS * 1000.0));  
 #endif  #endif
 #ifdef PRINT  
       printf ("%d ", cyx); mpn_print (dx, size);  
 #endif  
   
       t0 = cputime();        t0 = cputime();
       for (i = 0; i < TIMES; i++)        for (i = 0; i < TIMES; i++)
         cyx = mpn_sub_n (dx, s1, s2, size);          mpn_sub_n (dx+1, s1, s2, size);
       t = cputime() - t0;        t = cputime() - t0;
 #if TIMES != 1        printf ("mpn_sub_n:    %5ldms (%.2f cycles/limb)\n",
       printf ("mpn_sub_n:   %ldms (%.2f cycles/limb)\n",                t, ((double) t * CLOCK) / (OPS * 1000.0));
               t,  
               ((double) t * CLOCK) / (OPS * 1000.0));  
 #endif  #endif
   
   #ifndef NOCHECK
         mpn_random2 (s1, size);
         mpn_random2 (s2, size);
   
 #ifdef PRINT  #ifdef PRINT
       printf ("%d ", cyx); mpn_print (dx, size);        mpn_print (s1, size);
         mpn_print (s2, size);
 #endif  #endif
   
 #ifndef NOCHECK  
       /* Put garbage in the destination.  */        /* Put garbage in the destination.  */
       for (i = 0; i < size; i++)        for (i = 0; i < size; i++)
         {          {
           dx[i] = 0x7654321;            dx[i+1] = 0xdead;
           dy[i] = 0x1234567;            dy[i+1] = 0xbeef;
         }          }
   
       cyx = refmpn_sub_n (dx, s1, s2, size);        cyx = refmpn_sub_n (dx+1, s1, s2, size);
       cyy = mpn_sub_n (dy, s1, s2, size);        cyy = mpn_sub_n (dy+1, s1, s2, size);
       if (cyx != cyy || mpn_cmp (dx, dy, size) != 0  #ifdef PRINT
           || dx[size] != 0x12345678 || dy[size] != 0x12345678)        printf ("%d ", cyx); mpn_print (dx+1, size);
         printf ("%d ", cyy); mpn_print (dy+1, size);
   #endif
         if (cyx != cyy || mpn_cmp (dx, dy, size+2) != 0
             || dx[0] != 0x87654321 || dx[size+1] != 0x12345678)
         {          {
 #ifndef PRINT  #ifndef PRINT
           printf ("%d ", cyx); mpn_print (dx, size);            printf ("%d ", cyx); mpn_print (dx+1, size);
           printf ("%d ", cyy); mpn_print (dy, size);            printf ("%d ", cyy); mpn_print (dy+1, size);
 #endif  #endif
             printf ("TEST NUMBER %d\n", test);
           abort();            abort();
         }          }
 #endif  #endif
Line 201  mpn_print (mp_ptr p, mp_size_t size)
Line 226  mpn_print (mp_ptr p, mp_size_t size)
   
   for (i = size - 1; i >= 0; 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]);        printf ("%0*lX", (int) (2 * sizeof(mp_limb_t)), p[i]);
   #endif
 #ifdef SPACE  #ifdef SPACE
       if (i != 0)        if (i != 0)
         printf (" ");          printf (" ");
Legend:
Removed from v.1.1.1.1  
changed lines
  Added in v.1.1.1.2
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>