annotate src/fftw-3.3.5/libbench2/bench-main.c @ 42:2cd0e3b3e1fd

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