blob: b8e1146c69f27edb63a18f6da1d8f5a1d5b00fc9 (
plain) (
blame)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
|
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <time.h>
#include <inttypes.h>
#include <libfixmath/fixmath.h>
#include "hiclock.h"
int main(int argc, char** argv) {
printf("libfixmath test tool\n");
hiclock_init();
uintptr_t iter = (1 << 16);
uintptr_t pass = (1 << 8);
uintptr_t i;
srand(time(NULL));
fix16_t fix_args[iter];
for(i = 0; i < iter; i++)
fix_args[i] = (rand() ^ (rand() << 16));
fix16_t fix_result[iter];
hiclock_t fix_start = hiclock();
for(i = 0; i < pass; i++) {
uintptr_t j;
for(j = 0; j < iter; j++)
fix_result[j] = fix16_atan(fix_args[j]);
}
hiclock_t fix_end = hiclock();
float flt_args[iter];
for(i = 0; i < iter; i++)
flt_args[i] = fix16_to_float(fix_args[i]);
float flt_result[iter];
hiclock_t flt_start = hiclock();
for(i = 0; i < pass; i++) {
uintptr_t j;
for(j = 0; j < iter; j++)
flt_result[j] = atanf(flt_args[j]);
}
hiclock_t flt_end = hiclock();
fix16_t fix_error = 0;
for(i = 0; i < iter; i++)
fix_error += abs(fix16_from_float(flt_result[i]) - fix_result[i]);
hiclock_t flt_duration = (flt_end - flt_start);
hiclock_t fix_duration = (fix_end - fix_start);
printf("Floating Point: %08"PRIuHICLOCK" @ %"PRIu32"Hz\n", flt_duration, HICLOCKS_PER_SEC);
printf("Fixed Point: %08"PRIuHICLOCK" @ %"PRIu32"Hz\n", fix_duration, HICLOCKS_PER_SEC);
printf("Difference: %08"PRIiHICLOCK" (% 3.2f%)\n", (flt_duration - fix_duration), ((fix_duration * 100.0) / flt_duration));
printf("Error: %f%%\n", ((fix16_to_dbl(fix_error) * 100.0) / iter));
return EXIT_SUCCESS;
}
|
