Annotation of OpenXM_contrib/gmp/mpn/tests/rshift.c, Revision 1.1.1.2
1.1.1.2 ! maekawa 1: /*
! 2: Copyright (C) 1996, 1998, 1999, 2000 Free Software Foundation, Inc.
! 3:
! 4: This file is part of the GNU MP Library.
! 5:
! 6: The GNU MP Library is free software; you can redistribute it and/or modify
! 7: it under the terms of the GNU Lesser General Public License as published by
! 8: the Free Software Foundation; either version 2.1 of the License, or (at your
! 9: option) any later version.
! 10:
! 11: The GNU MP Library is distributed in the hope that it will be useful, but
! 12: WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
! 13: or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public
! 14: License for more details.
! 15:
! 16: You should have received a copy of the GNU Lesser General Public License
! 17: along with the GNU MP Library; see the file COPYING.LIB. If not, write to
! 18: the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
! 19: MA 02111-1307, USA.
! 20: */
! 21:
1.1 maekawa 22: #include <stdio.h>
23: #include "gmp.h"
24: #include "gmp-impl.h"
25:
1.1.1.2 ! maekawa 26: #if defined (USG) || defined (__SVR4) || defined (_UNICOS) || defined (__hpux)
! 27: #include <time.h>
1.1 maekawa 28:
1.1.1.2 ! maekawa 29: int
1.1 maekawa 30: cputime ()
31: {
1.1.1.2 ! maekawa 32: if (CLOCKS_PER_SEC < 100000)
! 33: return clock () * 1000 / CLOCKS_PER_SEC;
! 34: return clock () / (CLOCKS_PER_SEC / 1000);
1.1 maekawa 35: }
36: #else
1.1.1.2 ! maekawa 37: #include <sys/types.h>
! 38: #include <sys/time.h>
! 39: #include <sys/resource.h>
1.1 maekawa 40:
1.1.1.2 ! maekawa 41: int
1.1 maekawa 42: cputime ()
43: {
1.1.1.2 ! maekawa 44: struct rusage rus;
! 45:
! 46: getrusage (0, &rus);
! 47: return rus.ru_utime.tv_sec * 1000 + rus.ru_utime.tv_usec / 1000;
1.1 maekawa 48: }
49: #endif
50:
51: #define M * 1000000
52:
53: #ifndef CLOCK
54: #if defined (__m88k__)
55: #define CLOCK 20 M
56: #elif defined (__i386__)
1.1.1.2 ! maekawa 57: #define CLOCK (16666667)
1.1 maekawa 58: #elif defined (__m68k__)
59: #define CLOCK (20 M)
60: #elif defined (_IBMR2)
61: #define CLOCK (25 M)
62: #elif defined (__sparc__)
63: #define CLOCK (20 M)
64: #elif defined (__sun__)
65: #define CLOCK (20 M)
66: #elif defined (__mips)
67: #define CLOCK (40 M)
68: #elif defined (__hppa__)
69: #define CLOCK (50 M)
70: #elif defined (__alpha)
71: #define CLOCK (133 M)
72: #else
73: #error "Don't know CLOCK of your machine"
74: #endif
75: #endif
76:
77: #ifndef OPS
1.1.1.2 ! maekawa 78: #define OPS (CLOCK/5)
1.1 maekawa 79: #endif
80: #ifndef SIZE
81: #define SIZE 496
82: #endif
83: #ifndef TIMES
84: #define TIMES OPS/SIZE
85: #else
86: #undef OPS
87: #define OPS (SIZE*TIMES)
88: #endif
89:
90: mp_limb_t
91: refmpn_rshift (wp, up, usize, cnt)
92: register mp_ptr wp;
93: register mp_srcptr up;
94: mp_size_t usize;
95: register unsigned int cnt;
96: {
97: register mp_limb_t high_limb, low_limb;
98: register unsigned sh_1, sh_2;
99: register mp_size_t i;
100: mp_limb_t retval;
101:
102: #ifdef DEBUG
103: if (usize == 0 || cnt == 0)
104: abort ();
105: #endif
106:
107: sh_1 = cnt;
108: #if 0
109: if (sh_1 == 0)
110: {
111: if (wp != up)
112: {
113: /* Copy from low end to high end, to allow specified input/output
114: overlapping. */
115: for (i = 0; i < usize; i++)
116: wp[i] = up[i];
117: }
118: return 0;
119: }
120: #endif
121:
122: wp -= 1;
123: sh_2 = BITS_PER_MP_LIMB - sh_1;
124: high_limb = up[0];
125: retval = high_limb << sh_2;
126: low_limb = high_limb;
127:
128: for (i = 1; i < usize; i++)
129: {
130: high_limb = up[i];
131: wp[i] = (low_limb >> sh_1) | (high_limb << sh_2);
132: low_limb = high_limb;
133: }
134: low_limb >>= sh_1;
135: wp[i] = low_limb;
136:
137: return retval;
138: }
139:
140: #ifndef CNT
141: #define CNT 4
142: #endif
143:
144: main (argc, argv)
145: int argc;
146: char **argv;
147: {
148: mp_limb_t s1[SIZE];
149: mp_limb_t dx[SIZE+2];
150: mp_limb_t dy[SIZE+2];
151: mp_limb_t cyx, cyy;
152: int i;
153: long t0, t;
154: int test;
155: int cnt = CNT;
156: mp_size_t size;
157:
158: for (test = 0; ; test++)
159: {
1.1.1.2 ! maekawa 160: #if TIMES == 1 && ! defined (PRINT)
! 161: if (test % (SIZE > 10000 ? 1 : 10000 / SIZE) == 0)
! 162: {
! 163: printf ("\r%d", test);
! 164: fflush (stdout);
! 165: }
! 166: #endif
! 167:
! 168: #if TIMES == 1
! 169: cnt = random () % (BITS_PER_MP_LIMB - 1) + 1;
! 170: #endif
! 171:
1.1 maekawa 172: #ifdef RANDOM
1.1.1.2 ! maekawa 173: size = random () % SIZE + 1;
1.1 maekawa 174: #else
175: size = SIZE;
176: #endif
177: mpn_random2 (s1, size);
178:
179: dx[size+1] = 0x12345678;
180: dy[size+1] = 0x12345678;
181: dx[0] = 0x87654321;
182: dy[0] = 0x87654321;
183:
1.1.1.2 ! maekawa 184: #if TIMES != 1
! 185: mpn_random (s1, size);
! 186:
! 187: #ifndef NOCHECK
1.1 maekawa 188: t0 = cputime();
189: for (i = 0; i < TIMES; i++)
1.1.1.2 ! maekawa 190: refmpn_rshift (dx+1, s1, size, cnt);
1.1 maekawa 191: t = cputime() - t0;
192: printf ("refmpn_rshift: %5ldms (%.2f cycles/limb)\n",
1.1.1.2 ! maekawa 193: t, ((double) t * CLOCK) / (OPS * 1000.0));
1.1 maekawa 194: #endif
195:
196: t0 = cputime();
197: for (i = 0; i < TIMES; i++)
1.1.1.2 ! maekawa 198: mpn_rshift (dx+1, s1, size, cnt);
1.1 maekawa 199: t = cputime() - t0;
1.1.1.2 ! maekawa 200: printf ("mpn_rshift: %5ldms (%.2f cycles/limb)\n",
! 201: t, ((double) t * CLOCK) / (OPS * 1000.0));
1.1 maekawa 202: #endif
1.1.1.2 ! maekawa 203:
! 204: #ifndef NOCHECK
! 205: mpn_random2 (s1, size);
! 206:
1.1 maekawa 207: #ifdef PRINT
1.1.1.2 ! maekawa 208: printf ("%-*d ", (int) (2 * sizeof(mp_limb_t)), cnt); mpn_print (s1, size);
1.1 maekawa 209: #endif
210:
211: /* Put garbage in the destination. */
1.1.1.2 ! maekawa 212: for (i = 0; i < size; i++)
1.1 maekawa 213: {
1.1.1.2 ! maekawa 214: dx[i+1] = 0xdead;
! 215: dy[i+1] = 0xbeef;
1.1 maekawa 216: }
217:
218: cyx = refmpn_rshift (dx+1, s1, size, cnt);
219: cyy = mpn_rshift (dy+1, s1, size, cnt);
1.1.1.2 ! maekawa 220: #ifdef PRINT
! 221: printf ("%*lX ", (int) (2 * sizeof(mp_limb_t)), cyx);
! 222: mpn_print (dx+1, size);
! 223: printf ("%*lX ", (int) (2 * sizeof(mp_limb_t)), cyy);
! 224: mpn_print (dy+1, size);
! 225: #endif
1.1 maekawa 226: if (cyx != cyy || mpn_cmp (dx, dy, size+2) != 0
1.1.1.2 ! maekawa 227: || dx[0] != 0x87654321 || dx[size+1] != 0x12345678)
1.1 maekawa 228: {
229: #ifndef PRINT
230: printf ("%*lX ", (int) (2 * sizeof(mp_limb_t)), cyx);
231: mpn_print (dx+1, size);
232: printf ("%*lX ", (int) (2 * sizeof(mp_limb_t)), cyy);
233: mpn_print (dy+1, size);
234: #endif
1.1.1.2 ! maekawa 235: printf ("TEST NUMBER %d\n", test);
! 236: if (dy[size+1] != 0x12345678)
! 237: printf ("clobbered at high end\n");
! 238: if (dy[0] != 0x87654321)
! 239: printf ("clobbered at low end\n");
1.1 maekawa 240: abort();
241: }
242: #endif
243: }
244: }
245:
246: mpn_print (mp_ptr p, mp_size_t size)
247: {
248: mp_size_t i;
249:
250: for (i = size - 1; i >= 0; i--)
251: {
1.1.1.2 ! maekawa 252: #ifdef _LONG_LONG_LIMB
! 253: printf ("%0*lX%0*lX", (int) (sizeof(mp_limb_t)),
! 254: (unsigned long) (p[i] >> (BITS_PER_MP_LIMB/2)),
! 255: (int) (sizeof(mp_limb_t)), (unsigned long) (p[i]));
! 256: #else
1.1 maekawa 257: printf ("%0*lX", (int) (2 * sizeof(mp_limb_t)), p[i]);
1.1.1.2 ! maekawa 258: #endif
1.1 maekawa 259: #ifdef SPACE
260: if (i != 0)
261: printf (" ");
262: #endif
263: }
264: puts ("");
265: }
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>
![[BACK]](/icons/cvsweb/back.gif){kind=icon}
Return to rshift.c CVS log ![[TXT]](/icons/cvsweb/text.gif){kind=icon}
![[DIR]](/icons/cvsweb/dir.gif){kind=icon}
Up to [local] / OpenXM_contrib / gmp / mpn / tests