1b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar#include "timing.h"
2b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar#include <stdio.h>
3b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar
4b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar#define INPUT_TYPE int64_t
5b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar#define INPUT_SIZE 512
6b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar#define FUNCTION_NAME __floatdidf
7b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar
8b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar#ifndef LIBNAME
9b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar#define LIBNAME UNKNOWN
10b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar#endif
11b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar
12b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar#define LIBSTRING		LIBSTRINGX(LIBNAME)
13b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar#define LIBSTRINGX(a)	LIBSTRINGXX(a)
14b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar#define LIBSTRINGXX(a)	#a
15b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar
16b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbardouble FUNCTION_NAME(INPUT_TYPE x);
17b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar
18b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbarint main(int argc, char *argv[]) {
19b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar	INPUT_TYPE input[INPUT_SIZE];
20b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar	int i, j;
21b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar
22b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar	srand(42);
23b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar
24b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar	// Initialize the input array with data of various sizes.
25b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar	for (i=0; i<INPUT_SIZE; ++i)
26b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar		input[i] = (((uint64_t)rand() << 32) | (uint64_t)rand()) >> (rand() & 63);
27b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar
28b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar	double bestTime = __builtin_inf();
29b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar	void *dummyp;
30b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar	for (j=0; j<1024; ++j) {
31b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar
32b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar		uint64_t startTime = mach_absolute_time();
33b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar		for (i=0; i<INPUT_SIZE; ++i)
34b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar			FUNCTION_NAME(input[i]);
35b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar		uint64_t endTime = mach_absolute_time();
36b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar
37b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar		double thisTime = intervalInCycles(startTime, endTime);
38b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar		bestTime = __builtin_fmin(thisTime, bestTime);
39b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar
40b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar		// Move the stack alignment between trials to eliminate (mostly) aliasing effects
41b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar		dummyp = alloca(1);
42b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar	}
43b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar
44b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar	printf("%16s: %f cycles.\n", LIBSTRING, bestTime / (double) INPUT_SIZE);
45b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar
46b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar	return 0;
47b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar}
48