1b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar#include "timing.h" 2b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar#include <stdio.h> 3b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar 4b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar#define INPUT_TYPE uint64_t 5b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar#define INPUT_SIZE 512 6b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar#define FUNCTION_NAME __floatundisf 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 Dunbarfloat 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