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
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
|
#include <fpm/fixed.hpp>
#include <fpm/math.hpp>
#include <fixmath.h>
#include <iostream>
#include <fstream>
#include <cmath>
#include <string>
static constexpr double PI = 3.1415926535897932384626433832795;
using fpm::fixed_8_24;
using fpm::fixed_16_16;
using fpm::fixed_24_8;
using fixed_20_12 = fpm::fixed<std::int32_t, std::int64_t, 12>;
using std::sin;
using std::cos;
using std::tan;
using std::asin;
using std::acos;
using std::atan;
using std::atan2;
using std::sqrt;
static Fix16 sin(Fix16 x) { return x.sin(); }
static Fix16 cos(Fix16 x) { return x.cos(); }
static Fix16 tan(Fix16 x) { return x.tan(); }
static Fix16 asin(Fix16 x) { return x.asin(); }
static Fix16 acos(Fix16 x) { return x.acos(); }
static Fix16 atan(Fix16 x) { return x.atan(); }
static Fix16 atan2(Fix16 x, Fix16 y) { return x.atan2(y); }
static Fix16 sqrt(Fix16 x) { return x.sqrt(); }
static Fix16 exp(Fix16 x) { return fix16_exp(x); }
static Fix16 log(Fix16 x) { return fix16_log(x); }
static Fix16 log2(Fix16 x) { return fix16_log2(x); }
class csv_output
{
public:
explicit csv_output(const std::string& filename)
: m_stream(filename)
{
m_stream.setf(std::ios::fixed);
m_stream.precision(12);
m_stream << "x,real,Q24.8,Q20.12,Q16.16,Q8.24,fix16\n";
}
void write_row(double x, double y_real, double y_q24_8, double y_q20_12, double y_q16_16, double y_q8_24, double y_fix16)
{
m_stream << x << "," << y_real << "," << y_q24_8 << "," << y_q20_12 << "," << y_q16_16 << "," << y_q8_24 << "," << y_fix16 << "\n";
}
void write_row(double x, double y_real, double y_q24_8, double y_q20_12, double y_q16_16, double y_q8_24)
{
m_stream << x << "," << y_real << "," << y_q24_8 << "," << y_q20_12 << "," << y_q16_16 << "," << y_q8_24 << ",-\n";
}
private:
std::ofstream m_stream;
};
template <typename Callable, typename... Args>
static void check_all(csv_output& output, double value, Callable&& callable, Args&& ...args)
{
output.write_row(value,
callable(std::forward<Args>(args)...),
static_cast<double>(callable(fixed_24_8(std::forward<Args>(args))...)),
static_cast<double>(callable(fixed_20_12(std::forward<Args>(args))...)),
static_cast<double>(callable(fixed_16_16(std::forward<Args>(args))...)),
static_cast<double>(callable(fixed_8_24(std::forward<Args>(args))...)),
static_cast<double>(callable(Fix16(std::forward<Args>(args))...)));
}
template <typename Callable, typename... Args>
static void check_fpm(csv_output& output, double value, Callable&& callable, Args&& ...args)
{
output.write_row(value,
callable(std::forward<Args>(args)...),
static_cast<double>(callable(fixed_24_8(std::forward<Args>(args))...)),
static_cast<double>(callable(fixed_20_12(std::forward<Args>(args))...)),
static_cast<double>(callable(fixed_16_16(std::forward<Args>(args))...)),
static_cast<double>(callable(fixed_8_24(std::forward<Args>(args))...)));
}
int main()
{
csv_output out_sin("sin.csv");
csv_output out_cos("cos.csv");
csv_output out_tan("tan.csv");
csv_output out_atan2("atan2.csv");
for (int angle = -179; angle <= 180; ++angle)
{
const double val = angle * PI / 180.0;
check_all(out_sin, val, [](auto x) { return sin(x); }, val);
check_all(out_cos, val, [](auto x) { return cos(x); }, val);
if ((angle + 90) % 180 != 0)
{
check_all(out_tan, val, [](auto x) { return tan(x); }, val);
}
const auto y = std::sin(val);
const auto x = std::cos(val);
check_all(out_atan2, val, [](auto y, auto x) { return atan2(y, x); }, y, x);
}
csv_output out_asin("asin.csv");
csv_output out_acos("acos.csv");
for (int value = -100; value <= 100; ++value)
{
const double val = value / 100.0;
check_all(out_asin, val, [](auto x) { return asin(x); }, val);
check_all(out_acos, val, [](auto x) { return acos(x); }, val);
}
csv_output out_atan("atan.csv");
for (int value = -5000; value <= 5000; value += 5)
{
const double val = value / 1000.0;
check_all(out_atan, val, [](auto x) { return atan(x); }, val);
}
csv_output out_sqrt("sqrt.csv");
for (int i = 0; i < 1000; ++i)
{
const auto val = i / 10.0;
check_all(out_sqrt, val, [](auto x) { return sqrt(x); }, val);
}
csv_output out_cbrt("cbrt.csv");
for (int i = -1000; i < 1000; ++i)
{
const auto val = i / 10.0;
check_fpm(out_cbrt, val, [](auto x) { return cbrt(x); }, val);
}
csv_output out_exp("exp.csv");
csv_output out_exp2("exp2.csv");
csv_output out_pow("pow.csv");
for (int i = -40; i <= 40; i++)
{
const auto val = i / 10.0;
check_all(out_exp, val, [](auto x) { return exp(x); }, val);
check_fpm(out_exp2, val, [](auto x) { return exp2(x); }, val);
check_fpm(out_pow, val, [](auto x) { return pow(decltype(x){3.36}, x); }, val);
}
csv_output out_log("log.csv");
csv_output out_log2("log2.csv");
csv_output out_log10("log10.csv");
for (int i = 1; i < 1000; i++)
{
const auto val = i / 10.0;
check_all(out_log, val, [](auto x) { return log(x); }, val);
check_all(out_log2, val, [](auto x) { return log2(x); }, val);
check_fpm(out_log10, val, [](auto x) { return log10(x); }, val);
}
}
|
