cannam@127
|
1 /*
|
cannam@127
|
2 * Copyright (c) 2001 Matteo Frigo
|
cannam@127
|
3 * Copyright (c) 2001 Massachusetts Institute of Technology
|
cannam@127
|
4 *
|
cannam@127
|
5 * This program is free software; you can redistribute it and/or modify
|
cannam@127
|
6 * it under the terms of the GNU General Public License as published by
|
cannam@127
|
7 * the Free Software Foundation; either version 2 of the License, or
|
cannam@127
|
8 * (at your option) any later version.
|
cannam@127
|
9 *
|
cannam@127
|
10 * This program is distributed in the hope that it will be useful,
|
cannam@127
|
11 * but WITHOUT ANY WARRANTY; without even the implied warranty of
|
cannam@127
|
12 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
cannam@127
|
13 * GNU General Public License for more details.
|
cannam@127
|
14 *
|
cannam@127
|
15 * You should have received a copy of the GNU General Public License
|
cannam@127
|
16 * along with this program; if not, write to the Free Software
|
cannam@127
|
17 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
|
cannam@127
|
18 *
|
cannam@127
|
19 */
|
cannam@127
|
20
|
cannam@127
|
21
|
cannam@127
|
22 #include "bench.h"
|
cannam@127
|
23 #include "my-getopt.h"
|
cannam@127
|
24 #include <stdio.h>
|
cannam@127
|
25 #include <stdlib.h>
|
cannam@127
|
26
|
cannam@127
|
27 int verbose;
|
cannam@127
|
28
|
cannam@127
|
29 static const struct my_option options[] =
|
cannam@127
|
30 {
|
cannam@127
|
31 {"accuracy", REQARG, 'a'},
|
cannam@127
|
32 {"accuracy-rounds", REQARG, 405},
|
cannam@127
|
33 {"impulse-accuracy-rounds", REQARG, 406},
|
cannam@127
|
34 {"can-do", REQARG, 'd'},
|
cannam@127
|
35 {"help", NOARG, 'h'},
|
cannam@127
|
36 {"info", REQARG, 'i'},
|
cannam@127
|
37 {"info-all", NOARG, 'I'},
|
cannam@127
|
38 {"print-precision", NOARG, 402},
|
cannam@127
|
39 {"print-time-min", NOARG, 400},
|
cannam@127
|
40 {"random-seed", REQARG, 404},
|
cannam@127
|
41 {"report-benchmark", NOARG, 320},
|
cannam@127
|
42 {"report-mflops", NOARG, 300},
|
cannam@127
|
43 {"report-time", NOARG, 310},
|
cannam@127
|
44 {"report-verbose", NOARG, 330},
|
cannam@127
|
45 {"speed", REQARG, 's'},
|
cannam@127
|
46 {"setup-speed", REQARG, 'S'},
|
cannam@127
|
47 {"time-min", REQARG, 't'},
|
cannam@127
|
48 {"time-repeat", REQARG, 'r'},
|
cannam@127
|
49 {"user-option", REQARG, 'o'},
|
cannam@127
|
50 {"verbose", OPTARG, 'v'},
|
cannam@127
|
51 {"verify", REQARG, 'y'},
|
cannam@127
|
52 {"verify-rounds", REQARG, 401},
|
cannam@127
|
53 {"verify-tolerance", REQARG, 403},
|
cannam@127
|
54 {0, NOARG, 0}
|
cannam@127
|
55 };
|
cannam@127
|
56
|
cannam@127
|
57 int bench_main(int argc, char *argv[])
|
cannam@127
|
58 {
|
cannam@127
|
59 double tmin = 0.0;
|
cannam@127
|
60 double tol;
|
cannam@127
|
61 int repeat = 0;
|
cannam@127
|
62 int rounds = 10;
|
cannam@127
|
63 int iarounds = 0;
|
cannam@127
|
64 int arounds = 1; /* this is too low for precise results */
|
cannam@127
|
65 int c;
|
cannam@127
|
66
|
cannam@127
|
67 report = report_verbose; /* default */
|
cannam@127
|
68 verbose = 0;
|
cannam@127
|
69
|
cannam@127
|
70 tol = SINGLE_PRECISION ? 1.0e-3 : (QUAD_PRECISION ? 1e-29 : 1.0e-10);
|
cannam@127
|
71
|
cannam@127
|
72 main_init(&argc, &argv);
|
cannam@127
|
73
|
cannam@127
|
74 bench_srand(1);
|
cannam@127
|
75
|
cannam@127
|
76 while ((c = my_getopt (argc, argv, options)) != -1) {
|
cannam@127
|
77 switch (c) {
|
cannam@127
|
78 case 't' :
|
cannam@127
|
79 tmin = strtod(my_optarg, 0);
|
cannam@127
|
80 break;
|
cannam@127
|
81 case 'r':
|
cannam@127
|
82 repeat = atoi(my_optarg);
|
cannam@127
|
83 break;
|
cannam@127
|
84 case 's':
|
cannam@127
|
85 timer_init(tmin, repeat);
|
cannam@127
|
86 speed(my_optarg, 0);
|
cannam@127
|
87 break;
|
cannam@127
|
88 case 'S':
|
cannam@127
|
89 timer_init(tmin, repeat);
|
cannam@127
|
90 speed(my_optarg, 1);
|
cannam@127
|
91 break;
|
cannam@127
|
92 case 'd':
|
cannam@127
|
93 report_can_do(my_optarg);
|
cannam@127
|
94 break;
|
cannam@127
|
95 case 'o':
|
cannam@127
|
96 useropt(my_optarg);
|
cannam@127
|
97 break;
|
cannam@127
|
98 case 'v':
|
cannam@127
|
99 if (verbose >= 0) { /* verbose < 0 disables output */
|
cannam@127
|
100 if (my_optarg)
|
cannam@127
|
101 verbose = atoi(my_optarg);
|
cannam@127
|
102 else
|
cannam@127
|
103 ++verbose;
|
cannam@127
|
104 }
|
cannam@127
|
105 break;
|
cannam@127
|
106 case 'y':
|
cannam@127
|
107 verify(my_optarg, rounds, tol);
|
cannam@127
|
108 break;
|
cannam@127
|
109 case 'a':
|
cannam@127
|
110 accuracy(my_optarg, arounds, iarounds);
|
cannam@127
|
111 break;
|
cannam@127
|
112 case 'i':
|
cannam@127
|
113 report_info(my_optarg);
|
cannam@127
|
114 break;
|
cannam@127
|
115 case 'I':
|
cannam@127
|
116 report_info_all();
|
cannam@127
|
117 break;
|
cannam@127
|
118 case 'h':
|
cannam@127
|
119 if (verbose >= 0) my_usage(argv[0], options);
|
cannam@127
|
120 break;
|
cannam@127
|
121
|
cannam@127
|
122 case 300: /* --report-mflops */
|
cannam@127
|
123 report = report_mflops;
|
cannam@127
|
124 break;
|
cannam@127
|
125
|
cannam@127
|
126 case 310: /* --report-time */
|
cannam@127
|
127 report = report_time;
|
cannam@127
|
128 break;
|
cannam@127
|
129
|
cannam@127
|
130 case 320: /* --report-benchmark */
|
cannam@127
|
131 report = report_benchmark;
|
cannam@127
|
132 break;
|
cannam@127
|
133
|
cannam@127
|
134 case 330: /* --report-verbose */
|
cannam@127
|
135 report = report_verbose;
|
cannam@127
|
136 break;
|
cannam@127
|
137
|
cannam@127
|
138 case 400: /* --print-time-min */
|
cannam@127
|
139 timer_init(tmin, repeat);
|
cannam@127
|
140 ovtpvt("%g\n", time_min);
|
cannam@127
|
141 break;
|
cannam@127
|
142
|
cannam@127
|
143 case 401: /* --verify-rounds */
|
cannam@127
|
144 rounds = atoi(my_optarg);
|
cannam@127
|
145 break;
|
cannam@127
|
146
|
cannam@127
|
147 case 402: /* --print-precision */
|
cannam@127
|
148 if (SINGLE_PRECISION)
|
cannam@127
|
149 ovtpvt("single\n");
|
cannam@127
|
150 else if (QUAD_PRECISION)
|
cannam@127
|
151 ovtpvt("quad\n");
|
cannam@127
|
152 else if (LDOUBLE_PRECISION)
|
cannam@127
|
153 ovtpvt("long-double\n");
|
cannam@127
|
154 else if (DOUBLE_PRECISION)
|
cannam@127
|
155 ovtpvt("double\n");
|
cannam@127
|
156 else
|
cannam@127
|
157 ovtpvt("unknown %d\n", sizeof(bench_real));
|
cannam@127
|
158 break;
|
cannam@127
|
159
|
cannam@127
|
160 case 403: /* --verify-tolerance */
|
cannam@127
|
161 tol = strtod(my_optarg, 0);
|
cannam@127
|
162 break;
|
cannam@127
|
163
|
cannam@127
|
164 case 404: /* --random-seed */
|
cannam@127
|
165 bench_srand(atoi(my_optarg));
|
cannam@127
|
166 break;
|
cannam@127
|
167
|
cannam@127
|
168 case 405: /* --accuracy-rounds */
|
cannam@127
|
169 arounds = atoi(my_optarg);
|
cannam@127
|
170 break;
|
cannam@127
|
171
|
cannam@127
|
172 case 406: /* --impulse-accuracy-rounds */
|
cannam@127
|
173 iarounds = atoi(my_optarg);
|
cannam@127
|
174 break;
|
cannam@127
|
175
|
cannam@127
|
176 case '?':
|
cannam@127
|
177 /* my_getopt() already printed an error message. */
|
cannam@127
|
178 cleanup();
|
cannam@127
|
179 return 1;
|
cannam@127
|
180
|
cannam@127
|
181 default:
|
cannam@127
|
182 abort ();
|
cannam@127
|
183 }
|
cannam@127
|
184 }
|
cannam@127
|
185
|
cannam@127
|
186 /* assume that any remaining arguments are problems to be
|
cannam@127
|
187 benchmarked */
|
cannam@127
|
188 while (my_optind < argc) {
|
cannam@127
|
189 timer_init(tmin, repeat);
|
cannam@127
|
190 speed(argv[my_optind++], 0);
|
cannam@127
|
191 }
|
cannam@127
|
192
|
cannam@127
|
193 cleanup();
|
cannam@127
|
194 return 0;
|
cannam@127
|
195 }
|