![[BACK]](/icons/cvsweb/back.gif){kind=icon}
Return to lshift.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|
Return to lshift.c CVS log | Up to [local] / OpenXM_contrib / gmp / mpn / tests |
1.1 maekawa 1: #include <stdio.h>
2: #include "gmp.h"
3: #include "gmp-impl.h"
4: #include "longlong.h"
5:
6: #ifndef USG
7: #include <sys/time.h>
8: #include <sys/resource.h>
9:
10: unsigned long
11: cputime ()
12: {
13: struct rusage rus;
14:
15: getrusage (0, &rus);
16: return rus.ru_utime.tv_sec * 1000 + rus.ru_utime.tv_usec / 1000;
17: }
18: #else
19: #include <time.h>
20:
21: #ifndef CLOCKS_PER_SEC
22: #define CLOCKS_PER_SEC 1000000
23: #endif
24:
25: #if CLOCKS_PER_SEC >= 10000
26: #define CLOCK_TO_MILLISEC(cl) ((cl) / (CLOCKS_PER_SEC / 1000))
27: #else
28: #define CLOCK_TO_MILLISEC(cl) ((cl) * 1000 / CLOCKS_PER_SEC)
29: #endif
30:
31: unsigned long
32: cputime ()
33: {
34: return CLOCK_TO_MILLISEC (clock ());
35: }
36: #endif
37:
38: #define M * 1000000
39:
40: #ifndef CLOCK
41: #if defined (__m88k__)
42: #define CLOCK 20 M
43: #elif defined (__i386__)
44: #define CLOCK (16.666667 M)
45: #elif defined (__m68k__)
46: #define CLOCK (20 M)
47: #elif defined (_IBMR2)
48: #define CLOCK (25 M)
49: #elif defined (__sparc__)
50: #define CLOCK (20 M)
51: #elif defined (__sun__)
52: #define CLOCK (20 M)
53: #elif defined (__mips)
54: #define CLOCK (40 M)
55: #elif defined (__hppa__)
56: #define CLOCK (50 M)
57: #elif defined (__alpha)
58: #define CLOCK (133 M)
59: #else
60: #error "Don't know CLOCK of your machine"
61: #endif
62: #endif
63:
64: #ifndef OPS
65: #define OPS 10000000
66: #endif
67: #ifndef SIZE
68: #define SIZE 496
69: #endif
70: #ifndef TIMES
71: #define TIMES OPS/SIZE
72: #else
73: #undef OPS
74: #define OPS (SIZE*TIMES)
75: #endif
76:
77: mp_limb_t
78: refmpn_lshift (wp, up, usize, cnt)
79: register mp_ptr wp;
80: register mp_srcptr up;
81: mp_size_t usize;
82: register unsigned int cnt;
83: {
84: register mp_limb_t high_limb, low_limb;
85: register unsigned sh_1, sh_2;
86: register mp_size_t i;
87: mp_limb_t retval;
88:
89: #ifdef DEBUG
90: if (usize == 0 || cnt == 0)
91: abort ();
92: #endif
93:
94: sh_1 = cnt;
95: #if 0
96: if (sh_1 == 0)
97: {
98: if (wp != up)
99: {
100: /* Copy from high end to low end, to allow specified input/output
101: overlapping. */
102: for (i = usize - 1; i >= 0; i--)
103: wp[i] = up[i];
104: }
105: return 0;
106: }
107: #endif
108:
109: wp += 1;
110: sh_2 = BITS_PER_MP_LIMB - sh_1;
111: i = usize - 1;
112: low_limb = up[i];
113: retval = low_limb >> sh_2;
114: high_limb = low_limb;
115: while (--i >= 0)
116: {
117: low_limb = up[i];
118: wp[i] = (high_limb << sh_1) | (low_limb >> sh_2);
119: high_limb = low_limb;
120: }
121: wp[i] = high_limb << sh_1;
122:
123: return retval;
124: }
125:
126: #ifndef CNT
127: #define CNT 4
128: #endif
129:
130: main (argc, argv)
131: int argc;
132: char **argv;
133: {
134: mp_limb_t s1[SIZE];
135: mp_limb_t dx[SIZE+2];
136: mp_limb_t dy[SIZE+2];
137: mp_limb_t cyx, cyy;
138: int i;
139: long t0, t;
140: int test;
141: int cnt = CNT;
142: mp_size_t size;
143:
144: for (test = 0; ; test++)
145: {
146: #ifdef RANDOM
147: size = (random () % SIZE + 1);
148: #else
149: size = SIZE;
150: #endif
151: mpn_random2 (s1, size);
152:
153: dx[size+1] = 0x12345678;
154: dy[size+1] = 0x12345678;
155: dx[0] = 0x87654321;
156: dy[0] = 0x87654321;
157:
158: #ifdef PRINT
159: mpn_print (s1, size);
160: #endif
161: t0 = cputime();
162: for (i = 0; i < TIMES; i++)
163: cyx = refmpn_lshift (dx+1, s1, size, cnt);
164: t = cputime() - t0;
165: #if TIMES != 1
166: printf ("refmpn_lshift: %5ldms (%.2f cycles/limb)\n",
167: t,
168: ((double) t * CLOCK) / (OPS * 1000.0));
169: #endif
170: #ifdef PRINT
171: printf ("%*lX ", (int) (2 * sizeof(mp_limb_t)), cyx); mpn_print (dx+1, size);
172: #endif
173:
174: t0 = cputime();
175: for (i = 0; i < TIMES; i++)
176: cyy = mpn_lshift (dx+1, s1, size, cnt);
177: t = cputime() - t0;
178: #if TIMES != 1
179: printf ("mpn_lshift: %5ldms (%.2f cycles/limb)\n",
180: t,
181: ((double) t * CLOCK) / (OPS * 1000.0));
182: #endif
183: #ifdef PRINT
184: printf ("%*lX ", (int) (2 * sizeof(mp_limb_t)), cyy); mpn_print (dx+1, size);
185: #endif
186:
187: #ifndef NOCHECK
188: /* Put garbage in the destination. */
189: for (i = 1; i <= size; i++)
190: {
191: dx[i] = 0x7654321;
192: dy[i] = 0x1234567;
193: }
194:
195: cyx = refmpn_lshift (dx+1, s1, size, cnt);
196: cyy = mpn_lshift (dy+1, s1, size, cnt);
197:
198: if (cyx != cyy || mpn_cmp (dx, dy, size+2) != 0
199: || dx[size+1] != 0x12345678 || dx[0] != 0x87654321)
200: {
201: #ifndef PRINT
202: printf ("%*lX ", (int) (2 * sizeof(mp_limb_t)), cyx);
203: mpn_print (dx+1, size);
204: printf ("%*lX ", (int) (2 * sizeof(mp_limb_t)), cyy);
205: mpn_print (dy+1, size);
206: #endif
207: abort();
208: }
209: #endif
210: }
211: }
212:
213: mpn_print (mp_ptr p, mp_size_t size)
214: {
215: mp_size_t i;
216:
217: for (i = size - 1; i >= 0; i--)
218: {
219: printf ("%0*lX", (int) (2 * sizeof(mp_limb_t)), p[i]);
220: #ifdef SPACE
221: if (i != 0)
222: printf (" ");
223: #endif
224: }
225: puts ("");
226: }