Merge from vendor branch LIBPCAP:
[dragonfly.git] / contrib / bc / lib / testmul.c
1 /* compute the crossover for recursive and simple multiplication */
2
3 #include <stdio.h>
4 #include <time.h>
5 #include "number.h"
6 #ifndef VARARGS
7 #include <stdarg.h>
8 #else
9 #include <varargs.h>
10 #endif
11
12 /* from number.c ... */
13 extern int mul_base_digits;
14 /* extern int mul_small_digits; */
15 extern bc_num _one_;
16
17 /* global variables */
18 int test_n = 1000;
19 int test_time = 30 * CLOCKS_PER_SEC;  /* 30 seconds */ 
20
21 /* Other things for number.c. */
22 int std_only;
23
24 void
25 out_of_memory()
26 {
27   fprintf (stderr, "Fatal error: Out of memory for malloc.\n");
28   exit (1);
29 }
30
31 /* Runtime error will  print a message and stop the machine. */
32
33 #ifndef VARARGS
34 #ifdef __STDC__
35 void
36 rt_error (char *mesg, ...)
37 #else
38 void
39 rt_error (mesg)
40      char *mesg;
41 #endif
42 #else
43 void
44 rt_error (mesg, va_alist)
45      char *mesg;
46 #endif
47 {
48   va_list args;
49   char error_mesg [255];
50
51 #ifndef VARARGS   
52   va_start (args, mesg);
53 #else
54   va_start (args);
55 #endif
56   vsprintf (error_mesg, mesg, args);
57   va_end (args);
58   
59   fprintf (stderr, "Runtime error: %s\n", error_mesg);
60 }
61
62 /* A runtime warning tells of some action taken by the processor that
63    may change the program execution but was not enough of a problem
64    to stop the execution. */
65
66 #ifndef VARARGS
67 #ifdef __STDC__
68 void
69 rt_warn (char *mesg, ...)
70 #else
71 void
72 rt_warn (mesg)
73      char *mesg;
74 #endif
75 #else
76 void
77 rt_warn (mesg, va_alist)
78      char *mesg;
79 #endif
80 {
81   va_list args;
82   char error_mesg [255];
83
84 #ifndef VARARGS   
85   va_start (args, mesg);
86 #else
87   va_start (args);
88 #endif
89   vsprintf (error_mesg, mesg, args);
90   va_end (args);
91
92   fprintf (stderr, "Runtime warning: %s\n", error_mesg);
93 }
94
95 void
96 out_char (int ch)
97 {
98   putchar (ch);
99 }
100
101 /* Time stuff !!! */
102
103 int
104 timeit ( bc_num a, bc_num b, int *n)
105 {
106   clock_t first;
107   int i, res;
108   bc_num c;
109
110   bc_init_num (&c);
111   first = clock();
112   *n = 0;
113   do {
114     for (i=0; i<test_n; i++)
115       bc_multiply(a,b,&c,0);
116     *n += test_n;
117      res = (int) (clock() - first);
118   } while (res < test_time);
119   return res;
120 }
121
122 int debug = 0;  /* Print debugging messages? */
123
124 int main (int argc, char **argv)
125 {
126   bc_num ten, num, expo, big;
127
128   int min, max, mid;
129
130 #if 0
131   int smallsize;
132 #endif
133
134   int n1, n2;
135   clock_t t1, t2;
136   float permul1, permul2;
137
138   /* args? */
139   if (argc > 1)
140     if (strcmp (argv[1], "-d") == 0)
141       debug = 1;
142
143   bc_init_numbers();
144   bc_init_num (&ten);
145   bc_init_num (&num);
146   bc_init_num (&expo);
147   bc_init_num (&big);
148   bc_int2num (&ten, 10);
149
150   if (debug)
151     fprintf (stderr, "Timings are for %d multiplies\n"
152                      "Minimum time is %d seconds\n", test_n,
153              test_time/CLOCKS_PER_SEC);
154
155   /* Two of the same size */
156   min = 10;
157   max = 500;
158
159   if (debug)
160     fprintf (stderr, "Testing numbers of the same length.\n");
161
162   while (min < max) {
163     mid = (min+max)/2;
164     if (debug) fprintf (stderr,"Checking %d...\n", mid);
165
166     bc_int2num (&expo, mid);
167     bc_raise (ten, expo, &num, 0);
168     bc_sub (num, _one_, &num, 0);
169
170     mul_base_digits = 2*mid+1;
171     t1 = timeit (num, num, &n1);
172     permul1 = (float)t1/(float)n1;
173
174     mul_base_digits = 2*mid-1;
175     t2 = timeit (num, num, &n2);
176     permul2 = (float)t2/(float)n2;
177
178     if (permul1 < permul2)
179       min = mid+1;
180     else
181       max = mid-1;
182
183     if (debug) {
184       fprintf (stderr, "n1 = %d :: n2 = %d\n", n1, n2);
185       fprintf (stderr, "p1 = %f :: p2 = %f\n", permul1, permul2);
186     }
187   }  
188
189   if (debug)
190     fprintf (stderr, "Base digits crossover at %d digits\n", min);
191   printf ("#define MUL_BASE_DIGITS %d\n", 2*min);
192
193
194 #if 0
195   mul_base_digits = min;
196
197   /* Small one times a big one. */
198
199   smallsize = min/2;
200   bc_int2num (&expo, smallsize);
201   bc_raise (ten, expo, &big, 0);
202   bc_sub (num, _one_, &big, 0);
203
204   min = min / 2;
205   max = 500;
206
207   if (debug)
208     fprintf (stderr, "Testing numbers of the different length.\n");
209
210   while (min < max) {
211     mid = (min+max)/2;
212     if (debug) fprintf (stderr, "Checking %d...\n", mid);
213
214     bc_int2num (&expo, mid-smallsize);
215     bc_raise (ten, expo, &num, 0);
216     bc_sub (num, _one_, &num, 0);
217
218     mul_small_digits = mid+1;
219     t1 = timeit (big, num, &n1);
220     permul1 = (float)t1/(float)n1;
221
222     mul_small_digits = mid-1;
223     t2 = timeit (big, num, &n2);
224     permul2 = (float)t2/(float)n2;
225
226     if (permul1 < permul2)
227       min = mid+1;
228     else
229       max = mid-1;
230
231     if (debug) {
232       fprintf (stderr, "n1 = %d :: n2 = %d\n", n1, n2);
233       fprintf (stderr, "p1 = %f :: p2 = %f\n", permul1, permul2);
234     }
235   }  
236   
237   if (debug)
238     fprintf (stderr, "Non equal digits crossover at %d total digits\n", min);
239   printf ("#define MUL_SMALL_DIGITS = %d\n", min);
240
241 #endif
242
243   return 0;
244 }