Gitiles
Code Review Sign In
review.mlplatform.org / ml / ethos-u / ml-embedded-evaluation-kit / f32a86a6969508d7a156decbed0bfc9466ad92fa / . / tests / common / PlatformMathTests.cpp
blob: 923630080f7e8a082781f19f5f79d828a455435a [file] [log] [blame]
liabar01dee53bc2021-10-29 15:59:04 +0100[diff] [blame]1/*
Richard Burtonf32a86a2022-11-15 11:46:11 +0000[diff] [blame^]2 * SPDX-FileCopyrightText: Copyright 2021 - 2022 Arm Limited and/or its affiliates <open-source-office@arm.com>
liabar01dee53bc2021-10-29 15:59:04 +0100[diff] [blame]3 * SPDX-License-Identifier: Apache-2.0
4 *
5 * Licensed under the Apache License, Version 2.0 (the "License");
6 * you may not use this file except in compliance with the License.
7 * You may obtain a copy of the License at
8 *
9 * http://www.apache.org/licenses/LICENSE-2.0
10 *
11 * Unless required by applicable law or agreed to in writing, software
12 * distributed under the License is distributed on an "AS IS" BASIS,
13 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
14 * See the License for the specific language governing permissions and
15 * limitations under the License.
16 */
17#include "PlatformMath.hpp"
18#include <catch.hpp>
19#include <limits>
20#include <numeric>
21
22TEST_CASE("Test CosineF32")
23{
24 /*Test Constants: */
25 std::vector<double> inputA{
26 0.0, 0.01, 0.02, 0.03, 0.04, 0.05, 0.06, 0.07, 0.08, 0.09, 0.1,
27 0.11, 0.12, 0.13, 0.14, 0.15, 0.16, 0.17, 0.18, 0.19, 0.2, 0.21,
28 0.22, 0.23, 0.24, 0.25, 0.26, 0.27, 0.28, 0.29, 0.3, 0.31, 0.32,
29 0.33, 0.34, 0.35, 0.36, 0.37, 0.38, 0.39, 0.4, 0.41, 0.42, 0.43,
30 0.44, 0.45, 0.46, 0.47, 0.48, 0.49, 0.5, 0.51, 0.52, 0.53, 0.54,
31 0.55, 0.56, 0.57, 0.58, 0.59, 0.6, 0.61, 0.62, 0.63, 0.64, 0.65,
32 0.66, 0.67, 0.68, 0.69, 0.7, 0.71, 0.72, 0.73, 0.74, 0.75, 0.76,
33 0.77, 0.78, 0.79, 0.8, 0.81, 0.82, 0.83, 0.84, 0.85, 0.86, 0.87,
34 0.88, 0.89, 0.9, 0.91, 0.92, 0.93, 0.94, 0.95, 0.96, 0.97, 0.98,
35 0.99, 1.0
36 };
37 std::vector<double> expectedResult{
38 1.0, 0.9995065603657316, 0.9980267284282716, 0.99556196460308,
39 0.9921147013144779, 0.9876883405951378, 0.9822872507286887,
40 0.9759167619387474, 0.9685831611286311, 0.9602936856769431,
41 0.9510565162951535, 0.9408807689542255, 0.9297764858882515,
42 0.9177546256839811, 0.9048270524660195, 0.8910065241883679,
43 0.8763066800438636, 0.8607420270039436, 0.8443279255020151,
44 0.8270805742745618, 0.8090169943749475, 0.7901550123756904,
45 0.7705132427757893, 0.7501110696304596, 0.7289686274214116,
46 0.7071067811865476, 0.6845471059286886, 0.6613118653236518,
47 0.6374239897486896, 0.6129070536529766, 0.5877852522924731,
48 0.5620833778521306, 0.5358267949789965, 0.5090414157503712,
49 0.48175367410171516, 0.4539904997395468, 0.42577929156507266,
50 0.39714789063478056, 0.3681245526846781, 0.3387379202452915,
51 0.30901699437494745, 0.2789911060392295, 0.24868988716485496,
52 0.2181432413965427, 0.18738131458572474, 0.15643446504023092,
53 0.12533323356430426, 0.0941083133185145, 0.06279051952931353,
54 0.031410759078128396, 6.123233995736766e-17, -0.03141075907812828,
55 -0.0627905195293134, -0.09410831331851438, -0.12533323356430437,
56 -0.15643446504023104, -0.18738131458572482, -0.21814324139654234,
57 -0.24868988716485463, -0.27899110603922916, -0.30901699437494734,
58 -0.33873792024529137, -0.368124552684678, -0.39714789063478045,
59 -0.4257792915650727, -0.4539904997395467, -0.48175367410171543,
60 -0.5090414157503713, -0.5358267949789969, -0.5620833778521304,
61 -0.587785252292473, -0.6129070536529763, -0.6374239897486897,
62 -0.6613118653236517, -0.6845471059286887, -0.7071067811865475,
63 -0.7289686274214113, -0.7501110696304596, -0.7705132427757891,
64 -0.7901550123756904, -0.8090169943749473, -0.8270805742745619,
65 -0.8443279255020149, -0.8607420270039435, -0.8763066800438634,
66 -0.8910065241883678, -0.9048270524660194, -0.9177546256839811,
67 -0.9297764858882513, -0.9408807689542255, -0.9510565162951535,
68 -0.9602936856769431, -0.9685831611286311, -0.9759167619387474,
69 -0.9822872507286886, -0.9876883405951377, -0.9921147013144778,
70 -0.99556196460308, -0.9980267284282716, -0.9995065603657316, -1.0
71 };
72
73 float tolerance = 10e-7;
74 for (size_t i = 0; i < inputA.size(); i++) {
75 CHECK (expectedResult[i] ==
76 Approx(arm::app::math::MathUtils::CosineF32(M_PI*inputA[i])).margin(tolerance));
77 }
78}
79
80TEST_CASE("Test SineF32")
81{
82/*Test Constants: */
83std::vector<double> inputA{
84 0.0, 0.01, 0.02, 0.03, 0.04, 0.05, 0.06, 0.07, 0.08, 0.09, 0.1,
85 0.11, 0.12, 0.13, 0.14, 0.15, 0.16, 0.17, 0.18, 0.19, 0.2, 0.21,
86 0.22, 0.23, 0.24, 0.25, 0.26, 0.27, 0.28, 0.29, 0.3, 0.31, 0.32,
87 0.33, 0.34, 0.35, 0.36, 0.37, 0.38, 0.39, 0.4, 0.41, 0.42, 0.43,
88 0.44, 0.45, 0.46, 0.47, 0.48, 0.49, 0.5, 0.51, 0.52, 0.53, 0.54,
89 0.55, 0.56, 0.57, 0.58, 0.59, 0.6, 0.61, 0.62, 0.63, 0.64, 0.65,
90 0.66, 0.67, 0.68, 0.69, 0.7, 0.71, 0.72, 0.73, 0.74, 0.75, 0.76,
91 0.77, 0.78, 0.79, 0.8, 0.81, 0.82, 0.83, 0.84, 0.85, 0.86, 0.87,
92 0.88, 0.89, 0.9, 0.91, 0.92, 0.93, 0.94, 0.95, 0.96, 0.97, 0.98,
93 0.99, 1.0
94};
95std::vector<double> expectedResult{
96
97 0.0, 0.03141075907812829, 0.06279051952931337, 0.09410831331851431,
98 0.12533323356430426, 0.15643446504023087, 0.1873813145857246,
99 0.21814324139654256, 0.2486898871648548, 0.2789911060392293,
100 0.3090169943749474, 0.33873792024529137, 0.3681245526846779,
101 0.3971478906347806, 0.4257792915650727, 0.45399049973954675,
102 0.4817536741017153, 0.5090414157503713, 0.5358267949789967,
103 0.5620833778521306, 0.5877852522924731, 0.6129070536529764,
104 0.6374239897486896, 0.6613118653236518, 0.6845471059286886,
105 0.7071067811865475, 0.7289686274214116, 0.7501110696304596,
106 0.7705132427757893, 0.7901550123756903, 0.8090169943749475,
107 0.8270805742745618, 0.8443279255020151, 0.8607420270039436,
108 0.8763066800438637, 0.8910065241883678, 0.9048270524660196,
109 0.9177546256839811, 0.9297764858882513, 0.9408807689542255,
110 0.9510565162951535, 0.960293685676943, 0.9685831611286311,
111 0.9759167619387473, 0.9822872507286886, 0.9876883405951378,
112 0.9921147013144779, 0.99556196460308, 0.9980267284282716,
113 0.9995065603657316, 1.0, 0.9995065603657316,
114 0.9980267284282716, 0.99556196460308, 0.9921147013144778,
115 0.9876883405951377, 0.9822872507286886, 0.9759167619387474,
116 0.9685831611286312, 0.9602936856769431, 0.9510565162951536,
117 0.9408807689542255, 0.9297764858882513, 0.9177546256839813,
118 0.9048270524660195, 0.8910065241883679, 0.8763066800438635,
119 0.8607420270039436, 0.844327925502015, 0.827080574274562,
120 0.8090169943749475, 0.7901550123756905, 0.7705132427757893,
121 0.7501110696304597, 0.7289686274214114, 0.7071067811865476,
122 0.6845471059286888, 0.6613118653236518, 0.6374239897486899,
123 0.6129070536529764, 0.5877852522924732, 0.5620833778521305,
124 0.535826794978997, 0.5090414157503714, 0.4817536741017156,
125 0.45399049973954686, 0.4257792915650729, 0.3971478906347806,
126 0.36812455268467814, 0.3387379202452913, 0.3090169943749475,
127 0.2789911060392291, 0.24868988716485482, 0.21814324139654231,
128 0.18738131458572502, 0.15643446504023098, 0.12533323356430454,
129 0.09410831331851435, 0.06279051952931358, 0.031410759078128236,
130 1.2246467991473532e-16
131 };
132 float tolerance = 10e-4;
133 for (size_t i = 0; i < inputA.size(); i++) {
134 CHECK (expectedResult[i] ==
135 Approx(arm::app::math::MathUtils::SineF32(M_PI*inputA[i])).margin(tolerance));
136 }
137}
138
139TEST_CASE("Test SqrtF32")
140{
141 /*Test Constants: */
142 std::vector<float> inputA{0,1,2,9,M_PI};
Kshitij Sisodiab59ba682021-11-23 17:19:52 +0000[diff] [blame]143 size_t len = inputA.size();
liabar01dee53bc2021-10-29 15:59:04 +0100[diff] [blame]144 std::vector<float> expectedResult{0, 1, 1.414213562, 3, 1.772453851 };
145
146 for (size_t i=0; i < len; i++){
147 CHECK (expectedResult[i] == Approx(arm::app::math::MathUtils::SqrtF32(inputA[i])));
148 }
149}
150
151TEST_CASE("Test MeanF32")
152{
Richard Burtonc2911442022-04-22 09:08:21 +0100[diff] [blame]153 /* Test Constants: */
liabar01dee53bc2021-10-29 15:59:04 +0100[diff] [blame]154 std::vector<float> input
155 {0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 1.000};
156
157 /* Manually calculated mean of above vector */
158 float expectedResult = 0.100;
159 CHECK (expectedResult == Approx(arm::app::math::MathUtils::MeanF32(input.data(), input.size())));
Richard Burtonc2911442022-04-22 09:08:21 +0100[diff] [blame]160
161 /* Mean of 0 */
162 std::vector<float> input2{1, 2, -1, -2};
163 float expectedResult2 = 0.0f;
164 CHECK (expectedResult2 == Approx(arm::app::math::MathUtils::MeanF32(input2.data(), input2.size())));
165
166 /* All 0s */
167 std::vector<float> input3 = std::vector<float>(9, 0);
168 float expectedResult3 = 0.0f;
169 CHECK (expectedResult3 == Approx(arm::app::math::MathUtils::MeanF32(input3.data(), input3.size())));
170
171 /* All 1s */
172 std::vector<float> input4 = std::vector<float>(9, 1);
173 float expectedResult4 = 1.0f;
174 CHECK (expectedResult4 == Approx(arm::app::math::MathUtils::MeanF32(input4.data(), input4.size())));
liabar01dee53bc2021-10-29 15:59:04 +0100[diff] [blame]175}
176
177TEST_CASE("Test StdDevF32")
178{
179 /*Test Constants: */
180 /* Normally distributed sample data generated by numpy normal library */
181 std::vector<float> input
182 {1.62434536, -0.61175641, -0.52817175, -1.07296862, 0.86540763, -2.3015387,
183 1.74481176, -0.7612069, 0.3190391, -0.24937038, 1.46210794, -2.06014071,
184 -0.3224172, -0.38405435, 1.13376944, -1.09989127, -0.17242821, -0.87785842,
185 0.04221375, 0.58281521, -1.10061918, 1.14472371, 0.90159072, 0.50249434,
186 0.90085595, -0.68372786, -0.12289023, -0.93576943, -0.26788808, 0.53035547,
187 -0.69166075,-0.39675353, -0.6871727, -0.84520564, -0.67124613, -0.0126646,
188 -1.11731035, 0.2344157, 1.65980218, 0.74204416, -0.19183555, -0.88762896,
189 -0.74715829, 1.6924546, 0.05080775, -0.63699565, 0.19091548, 2.10025514,
190 0.12015895, 0.61720311
191 };
192 uint32_t inputLen = input.size();
193
194 /*Calculate mean using std library to avoid dependency on MathUtils::MeanF32 */
195 float mean = (std::accumulate(input.begin(), input.end(), 0.0f))/float(inputLen);
196 float output = arm::app::math::MathUtils::StdDevF32(input.data(), inputLen, mean);
197
198 /*Manually calculated standard deviation of above vector*/
199 float expectedResult = 0.969589282958136;
200
201 CHECK (expectedResult == Approx(output));
Richard Burtonc2911442022-04-22 09:08:21 +0100[diff] [blame]202
203 /* All 0s should have 0 std dev. */
204 std::vector<float> input2 = std::vector<float>(4, 0);
205 float expectedResult2 = 0.0f;
206 CHECK (expectedResult2 == Approx(arm::app::math::MathUtils::StdDevF32(input2.data(), input2.size(), 0.0f)));
207
208 /* All 1s should have 0 std dev. */
209 std::vector<float> input3 = std::vector<float>(4, 1);
210 float expectedResult3 = 0.0f;
211 CHECK (expectedResult3 == Approx(arm::app::math::MathUtils::StdDevF32(input3.data(), input3.size(), 1.0f)));
212
213 /* Manually calclualted std value */
214 std::vector<float> input4 {1, 2, 3, 4, 5, 6, 7, 8, 9, 0};
215 float mean2 = (std::accumulate(input4.begin(), input4.end(), 0.0f))/float(input4.size());
216 float expectedResult4 = 2.872281323;
217 CHECK (expectedResult4 == Approx(arm::app::math::MathUtils::StdDevF32(input4.data(), input4.size(), mean2)));
liabar01dee53bc2021-10-29 15:59:04 +0100[diff] [blame]218}
219
220TEST_CASE("Test FFT32")
221{
222 /*Test Constants: */
223 std::vector<float> input_zeros(512, 0);
224 std::vector<float> input_ones(512, 1);
225 /* Random numbers generated using numpy rand with range [0:1] */
226 std::vector<float> input_random{
227 0.42333686,0.6547418,0.8933691,0.91466254,0.5992143,0.99474055,0.97750413,0.97160685,
228 0.72718734,0.8699537,0.643911,0.09764466,0.0050113136,0.46823388,0.13709934,0.44892532,
229 0.59728205,0.04055081,0.5579888,0.18445836,0.66469765,0.82715863,0.91934484,0.7844356,
230 0.23489648,0.021708783,0.67819905,0.75761676,0.48374954,0.14006922,0.87082034,0.7694296,
231 0.80479276,0.8241704,0.95917296,0.5758142,0.16839339,0.34290153,0.5846108,0.6878044,
232 0.1067114,0.5198196,0.4356897,0.68049103,0.12480807,0.3538696,0.06067087,0.056964435,
233 0.5382167,0.07761527,0.6673144,0.9045368,0.11050189,0.03530183,0.07864744,0.98752064,
234 0.014321936,0.101833574,0.43293256,0.87102246,0.52411795,0.90232223,0.49560344,0.6803092,
235 0.2908511,0.14653015,0.99105513,0.7057098,0.09623502,0.039713606,0.88669086,0.56018597,
236 0.90632766,0.99241334,0.18748309,0.38991618,0.6359827,0.05665585,0.732304,0.2703365,
237 0.19014524,0.5017947,0.78862536,0.81253093,0.35050204,0.2832596,0.65221876,0.59856164,
238 0.42758793,0.78865635,0.30943435,0.93780816,0.62568265,0.35397422,0.84209913,0.48590583,
239 0.34837773,0.5811646,0.42924216,0.26692122,0.030709852,0.84459823,0.09085059,0.29297647,
240 0.48539516,0.33488297,0.7877257,0.8728821,0.28454545,0.7109578,0.86097074,0.8536262,
241 0.4978063,0.5760398,0.77506036,0.7716988,0.27041402,0.52340513,0.2055419,0.8728235,
242 0.13492358,0.79122984,0.52998376,0.33897072,0.6426309,0.8766521,0.89287037,0.74047667,
243 0.42341164,0.67437655,0.4682156,0.67123246,0.54287183,0.3580476,0.94756556,0.2699457,
244 0.6131569,0.75043845,0.8115012,0.49610943,0.7108478,0.90941435,0.02233071,0.37346774,
245 0.33732748,0.46691266,0.35784695,0.39391598,0.8556212,0.884142,0.11730601,0.550112,
246 0.31513855,0.69654715,0.58585805,0.4493127,0.78515726,0.8176612,0.9846698,0.32842383,
247 0.41843212,0.48470423,0.6757128,0.95876855,0.5989163,0.13587572,0.72886884,0.88291156,
248 0.34402263,0.66211045,0.86188424,0.21498202,0.26397392,0.67372984,0.91386956,0.7339788,
249 0.91308993,0.1953016,0.1539217,0.214701,0.58234113,0.8019992,0.63969976,0.041050985,
250 0.7293308,0.26341477,0.54768014,0.97596467,0.12385198,0.44149798,0.5519762,0.1697347,
251 0.577215,0.8213594,0.47874716,0.64515114,0.61467725,0.18463866,0.23890929,0.51052976,
252 0.16807361,0.53142565,0.2414274,0.41690814,0.98815554,0.6245643,0.9477003,0.24780034,
253 0.82469565,0.8614785,0.9565832,0.062440686,0.9710724,0.039196696,0.11030199,0.35234734,
254 0.02065066,0.12832293,0.7328055,0.48924434,0.17247158,0.5769348,0.44146806,0.53575355,
255 0.17258933,0.6980237,0.86494404,0.50573164,0.5033998,0.71199447,0.41353586,0.26767612,
256 0.3789118,0.046621118,0.58491063,0.22861995,0.03134273,0.53280216,0.23382367,0.07748905,
257 0.96875405,0.6613716,0.64087844,0.8377165,0.051519375,0.68997836,0.3776376,0.43362603,
258 0.5358754,0.51419014,0.12823892,0.26574057,0.508808,0.15734084,0.78327274,0.5045347,
259 0.5445746,0.89297736,0.8531272,0.91270804,0.87429863,0.3965137,0.13544834,0.74269205,
260 0.80592203,0.045050766,0.13362087,0.17090783,0.02873757,0.99339336,0.6394376,0.48203012,
261 0.70598215,0.37082237,0.39792424,0.89938444,0.312602,0.48755112,0.18220617,0.17303479,
262 0.31954846,0.78080165,0.1755106,0.68262285,0.84665287,0.8520143,0.8459509,0.39417005,
263 0.30087698,0.81362164,0.61927587,0.32739028,0.9023775,0.27578092,0.6830477,0.15842387,
264 0.8473049,0.43057114,0.2019703,0.20560141,0.6237757,0.60283095,0.27645138,0.26605442,
265 0.27985683,0.41353813,0.85139906,0.71711886,0.5444832,0.73613757,0.7397004,0.7406752,
266 0.41016674,0.31896713,0.4541723,0.2795807,0.47941738,0.00504193,0.89091027,0.8097144,
267 0.63033766,0.37252298,0.9132861,0.5102532,0.04104481,0.30368647,0.21573475,0.99520445,
268 0.5047808,0.6868845,0.99881023,0.30377692,0.2554386,0.47201005,0.11120686,0.10077732,
269 0.1853349,0.49159425,0.3938629,0.8989509,0.9887155,0.698771,0.695701,0.78368753,0.52537227,
270 0.19451462,0.3659248,0.1968508,0.7751828,0.33103722,0.40406147,0.37832898,0.68663514,
271 0.32225925,0.41771907,0.034218453,0.42808908,0.20685343,0.1861495,0.045986768,0.8532299,
272 0.17200677,0.44670314,0.56831235,0.5388232,0.5430553,0.69175136,0.6462231,0.42827028,
273 0.10050113,0.30627027,0.9967943,0.6684778,0.5928422,0.63392985,0.99123496,0.79301435,
274 0.7936309,0.42839453,0.39781123,0.22329247,0.0122212395,0.2807108,0.19812097,0.5576105,
275 0.115653396,0.3732018,0.7622857,0.19847734,0.5310287,0.7298145,0.5518292,0.9117333,
276 0.13215758,0.33716795,0.42372775,0.6779287,0.35799992,0.097887225,0.20171605,0.9948177,
277 0.1829232,0.80349857,0.9807098,0.22959666,0.67322475,0.63094735,0.93454355,0.15962408,
278 0.04335433,0.47104993,0.36784375,0.45258796,0.93415564,0.1655446,0.7195017,0.76236975,
279 0.3846913,0.01330617,0.84716374,0.1227003,0.65102947,0.6632434,0.3728453,0.4222391,
280 0.6942989,0.16014872,0.10798196,0.94033676,0.026525471,0.8379024,0.5484514,0.13500613,
281 0.22919805,0.7001831,0.6573261,0.38086265,0.8725666,0.35077834,0.28415123,0.42283052,
282 0.668379,0.9769895,0.37621376,0.646407,0.11188069,0.17129017,0.7441628,0.25617477,
283 0.7751679,0.8565412,0.67631435,0.45213568,0.61896557,0.3387995,0.51607716,0.60779697,
284 0.16428445,0.5080923,0.13012086,0.61184275,0.7690249,0.9578811,0.67365676,0.16241212,
285 0.97157824,0.5595742,0.75936574,0.6043881,0.2149638,0.4925318,0.58727825,0.97953695,
286 0.01605968,0.2819307,0.6448378,0.4265335,0.661541,0.3976571,0.40607136,0.46425515,
287 0.2055872,0.2716193,0.4132582,0.8372537,0.37787434,0.082228854,0.7985557,0.9718134,
288 0.35222608,0.4853643,0.2569464,0.14783978,0.4889042,0.62900156,0.19994198,0.4618481,
289 0.21673755,0.51749533,0.1260157,0.83759904,0.36438805,0.6704668,0.22010763,0.2359318,
290 0.53004104,0.9723652,0.91218954,0.9153926,0.48207277,0.34850466,0.8939421
291 };
292 std::vector<std::vector<float>> input_vectors{input_zeros,
293 input_ones,
294 input_random};
295
296 std::vector<float> output_zeros(512, 0);
297 std::vector<float> output_ones(512, 0);
298 std::vector<float> output_random(512, 0);
299 std::vector<std::vector<float>> output_vectors{output_zeros,
300 output_ones,
301 output_random};
302
303 std::vector<float> expected_result_zeros(512, 0);
304 std::vector<float> expected_result_ones(512, 0);
305 expected_result_ones[0] = 512.0;
306
307 /* Values are stored as [real0, realN/2, real1, im1, real2, im2, ...] */
308 std::vector<float> expected_result_random{
309 2.59510161e+02, 2.59796867, 2.55982143, -5.91349888,
310 -1.80049237, 1.09902763, 4.00943240, 7.76684892,
311 4.32617219, -4.33636417, 2.98128463, -0.83763449,
312 2.92973078, -3.75655459, 2.27203161, -4.61106145,
313 5.55562176, -4.71880166, 2.13693416, -2.20496619,
314 2.31174036, -3.52991041, -6.16870680e-01, 2.43455407,
315 4.76317833, 8.66518565, 1.72350562, -1.33641312,
316 5.82836675, 0.89396187, 8.15031483, 4.34599034,
317 5.99780199, 2.94900065, -2.34462045e-02, -5.03789597,
318 1.21907020e+01, -6.47012928, 1.24434715, 0.23621713,
319 9.20921279e-01, -9.20510398, -7.73267254e-01, -3.72141078,
320 6.28883709, 4.00634065, 4.46491682, 0.74625307,
321 -5.87506158e-01, -4.22833058, 3.15189786, -1.82518672,
322 -6.93782260, 2.4170692, 3.23045185, 3.33383799,
323 5.10059531e-02, -3.4233929, -2.91651323, -0.0258584,
324 5.84499843, -9.51454903, -1.49214047e+01, 5.52200123,
325 4.52179590, -7.08268703, 5.16775420e-01, -2.90878759,
326 5.04314682, -1.16928599, -1.07329243e+01, -0.2719951,
327 -3.95269565, -2.32475678, -4.11031641, -2.20538835,
328 -5.89005095e-01, -5.65483456, -1.08927018e+01, 5.74801823,
329 -5.72520347, -3.94970165, -5.18407515e-01, 1.23622633,
330 9.56297959, 0.24424306, -3.74306351, 9.63476301,
331 -4.74493837, -0.35443496, 6.54760504, 1.16188913,
332 -1.33416950e+01, 7.19088609, -2.56045800, -0.49557866,
333 2.93460322, 6.91076746, -2.84779221e-01, 6.59958391,
334 1.70963995, -7.67293252, 4.18500790, -0.14627552,
335 -1.24855113, -1.43322867, 3.60644904e-01, -4.11521374,
336 -1.00628421e+01, 2.87531563, -9.34809732, -0.58251846,
337 3.61799848, -5.10288284, -5.96239076, 1.99792128,
338 -7.83229243e-02, 1.81741166, 2.32709681, 0.68487206,
339 -3.08398468, 1.5177629 , 1.41015388, 4.51146401,
340 1.90769911, 0.56093423, 3.58389141, -0.14974575,
341 -2.20163907, -1.62177814, -1.91904127, 1.94645907,
342 -1.37722930e-01, -4.30291678, -2.61843435, 0.12691241,
343 3.28959117, 4.7309582 , -2.93995652, 0.39835926,
344 2.89711768, 1.42284586, -5.13129145, -7.26477374,
345 3.74616158, -2.59659457, 3.85748750, -2.93737277,
346 3.17748694, -4.45041455, -2.68466437, 1.37377726,
347 1.60008368, 1.63787578, -1.95661401, -6.34937202,
348 -2.62744282, -5.20892662, 8.90553959e-01, -6.37113573,
349 2.35885332,-10.04547561, 3.29866159e-01, 1.89217741,
350 8.82516491e-01, -3.53298728, -2.22525608, -5.64794388,
351 -5.19226843, 2.5971315 , 4.49346648, -0.20428409,
352 -6.14851885,-11.90893875, 3.75899776, -1.86910056,
353 2.78518535, -2.6359501 , -2.13423317, 4.86509946,
354 -2.37625499, 7.42404308, 6.71175474, 6.06191618,
355 2.59014379, 4.76329698, 9.19140042, 9.69149015,
356 3.33307819, -4.03094924, 2.12988453, -0.15820258,
357 -1.03422801e+01, -3.04462388, 3.59852152, -2.00887343,
358 -3.69998656, 0.90050102, 6.79959099e-01, -1.88604949,
359 1.24235316, 0.41309537, 6.13876866, -7.1040085 ,
360 6.17728674, 1.91667103, -1.32895472, -0.17674661,
361 -6.94720428, 3.10502593, -2.33990738, -1.27840434,
362 3.21442520, 2.14102714, 2.37498837, 3.8158066 ,
363 -2.24107675, -5.52527559, -2.95697930, -0.50367608,
364 3.01687661, 7.08195792, 6.78604790, -3.94154162,
365 2.24402195, 4.60132638, 3.42211139, 4.17689039,
366 -1.17277194, 2.15404472, -2.37481930, 1.42611867,
367 -4.63033506e-01, 3.21563035, 1.38662123, 3.98598717,
368 -3.75283402, -2.47600433, -1.97290542, 2.83361487,
369 -8.45662834e-01, 5.57411581, -9.72981483e-01,-11.394208,
370 1.88220611, -0.80225125, -4.34295854e-01, -8.2954126,
371 3.81795409, -3.17146 , -4.61994107, -1.59820505,
372 -5.98834455, -4.93129451, -5.13862996e-01, -0.15649305,
373 -5.59094391, 6.25244435, -6.59974456,13.17193115,
374 4.48609092, 1.64741879, 7.40985006, 0.44896188,
375 3.81058449, -0.76425931, -5.47938416, 4.01447941,
376 -3.21535548, -1.45542238, 7.22740830e-01, -0.23983128,
377 -4.32373034, 0.1337671 , -5.89365226, -3.18756318,
378 7.90979161, 5.27570134, -3.43094553, -6.00826981,
379 1.17932561, -3.50027177, 1.81306385e-01, 1.1062498 ,
380 7.23650536e-01, -1.55500613, -3.88047911, -2.43746762,
381 -6.81565579, 2.16343352, 2.46366137, 2.38704469,
382 -2.55106395, 6.5091449 , -2.06510578,11.11320924,
383 2.06649835, -1.05026064, 1.63564303, -0.04638729,
384 1.45053876, 0.43730146, 1.25027939, 0.79932743,
385 2.81088838, 6.95136058, -4.41417255, 2.89610628,
386 1.15426258, -2.60704937, -2.77744882, 4.12872365,
387 -2.98288336, 6.75607352, 2.36553382, -2.10540332,
388 -7.30042988, -5.44897893, 3.44048454, 4.29726231,
389 2.18199500, -0.80126759, -4.04051175, 4.57584864,
390 9.56312116e-01, -4.45183318, 3.42348929, -9.84138181,
391 8.69604433, -6.6481311 , 4.68232735e-01, 1.41031176,
392 -1.24085700,10.61672181, 3.56591473e-01,10.51631421,
393 -9.97435470e-01, 2.72157537, 8.63583929, -2.19404252,
394 -1.53605811, -4.41068581, 2.05371873, 1.25665769,
395 1.65289503, 4.52520582, -5.35062642e-01, 0.82084677,
396 -1.10079476e+01, -5.09361474, 9.63129107, 3.90056638,
397 -4.19779738, 0.06565745, 2.42526917, 2.5854233 ,
398 -3.66709357, 3.80502971, -1.01489353e-01, -6.85423228,
399 -1.39361494e+01, -0.43904617, 6.01800968, -1.30751495,
400 -4.75122234, 2.74740671, 5.54971138, 9.43409003,
401 -9.94733058e-01, -3.0096825 , -8.60263376, 0.36653762,
402 3.53318614, -2.69194556, -8.95145740, -4.71570923,
403 5.15417709, 2.68645385, -2.78042293, 8.21739385,
404 5.90225003e-01, 2.13319153, 1.72158888, 0.18114627,
405 3.92269446, 3.3525857 , -3.40313825, 4.39280934,
406 -1.70368966, 1.29121245, -3.11326453, -1.85941318,
407 4.57078881, 0.72531039, 6.00445664, 4.9588524 ,
408 -3.32944491, -0.02080722, -7.42374632, -3.23290026,
409 -8.16145790e-01, 3.55935439, -6.19206533e-01, 2.42859073,
410 2.21486456, 3.76402487, 3.90930695, -3.61610186,
411 -8.12712547e-01,-14.63377988, 1.14460823, -3.14089899,
412 3.18097435, 1.21957751, 2.85181833, 0.89990235,
413 -4.32147361, -5.54219361, -1.12253677, 2.96141081,
414 -4.42577070, 3.17282306, 4.91746710, 1.16977744,
415 -4.55148089, -2.82520179, -1.71684103, 1.91487668,
416 7.70726836e-01, 0.78534837, -5.91048566, -4.8288477 ,
417 -1.35560162, -3.60938315, 1.15812301, 2.44299541,
418 1.36115190, -5.40950935, 7.08292127, 0.27720591,
419 -1.60210828e-01, 2.75862348, -1.57403782, 9.97207524,
420 -2.08957576, 8.70299964, -5.33004663, 4.1547783 ,
421 3.51580675, -5.10788085, 4.37938353, -3.73449894,
422 1.44673271, 0.51941469, 8.52232446e-01, -1.1134965 ,
423 -1.43972745, -1.62952127, -2.50759973, 1.19012213,
424 5.72772282e-01, -2.71833059, -6.84718990, 4.2621535 ,
425 1.58954734, -0.53827818, -1.44624396e-01, 7.63866979,
426 4.10423977e-01, -2.4785678 , -5.02681867, -2.03469811,
427 9.59505727e-01, 2.68589705, 3.20889444,-10.76452533,
428 -3.84771551, 2.49189796, -3.19895938, -3.49948794,
429 -2.67238970, 5.11386526, -3.85957031, -1.40741978,
430 1.76663166e-01,-11.7111276 , 6.39997364e-01, -1.30321198,
431 3.20767633, 1.65750671, -1.16187257e+01, -4.36634782,
432 -3.18675281, -4.89279155, -4.08760307, 2.19269283,
433 -1.58924870, 0.17948212, 4.81376107, 2.01871001,
434 -3.24211095e-01, -0.2790092 , 1.12603878, -3.61503491,
435 -2.86982317, -3.03634532, 8.07713910, 2.21302089,
436 2.91496011, -2.58564072, -2.45868054, -4.69447438,
437 -2.45868054, 4.69447438,
438 2.91496011, 2.58564072, 8.07713910, -2.21302089,
439 -2.86982317, 3.03634532, 1.12603878, 3.61503491,
440 -3.24211095e-01, 0.2790092 , 4.81376107, -2.01871001,
441 -1.58924870, -0.17948212, -4.08760307, -2.19269283,
442 -3.18675281, 4.89279155, -1.16187257e+01, 4.36634782,
443 3.20767633, -1.65750671, 6.39997364e-01, 1.30321198,
444 1.76663166e-01,11.7111276 , -3.85957031, 1.40741978,
445 -2.67238970, -5.11386526, -3.19895938, 3.49948794,
446 -3.84771551, -2.49189796, 3.20889444,10.76452533,
447 9.59505727e-01, -2.68589705, -5.02681867, 2.03469811,
448 4.10423977e-01, 2.4785678 , -1.44624396e-01, -7.63866979,
449 1.58954734, 0.53827818, -6.84718990, -4.2621535 ,
450 5.72772282e-01, 2.71833059, -2.50759973, -1.19012213,
451 -1.43972745, 1.62952127, 8.52232446e-01, 1.1134965 ,
452 1.44673271, -0.51941469, 4.37938353, 3.73449894,
453 3.51580675, 5.10788085, -5.33004663, -4.1547783 ,
454 -2.08957576, -8.70299964, -1.57403782, -9.97207524,
455 -1.60210828e-01, -2.75862348, 7.08292127, -0.27720591,
456 1.36115190, 5.40950935, 1.15812301, -2.44299541,
457 -1.35560162, 3.60938315, -5.91048566, 4.8288477 ,
458 7.70726836e-01, -0.78534837, -1.71684103, -1.91487668,
459 -4.55148089, 2.82520179, 4.91746710, -1.16977744,
460 -4.42577070, -3.17282306, -1.12253677, -2.96141081,
461 -4.32147361, 5.54219361, 2.85181833, -0.89990235,
462 3.18097435, -1.21957751, 1.14460823, 3.14089899,
463 -8.12712547e-01,14.63377988, 3.90930695, 3.61610186,
464 2.21486456, -3.76402487, -6.19206533e-01, -2.42859073,
465 -8.16145790e-01, -3.55935439, -7.42374632, 3.23290026,
466 -3.32944491, 0.02080722, 6.00445664, -4.9588524 ,
467 4.57078881, -0.72531039, -3.11326453, 1.85941318,
468 -1.70368966, -1.29121245, -3.40313825, -4.39280934,
469 3.92269446, -3.3525857 , 1.72158888, -0.18114627,
470 5.90225003e-01, -2.13319153, -2.78042293, -8.21739385,
471 5.15417709, -2.68645385, -8.95145740, 4.71570923,
472 3.53318614, 2.69194556, -8.60263376, -0.36653762,
473 -9.94733058e-01, 3.0096825 , 5.54971138, -9.43409003,
474 -4.75122234, -2.74740671, 6.01800968, 1.30751495,
475 -1.39361494e+01, 0.43904617, -1.01489353e-01, 6.85423228,
476 -3.66709357, -3.80502971, 2.42526917, -2.5854233 ,
477 -4.19779738, -0.06565745, 9.63129107, -3.90056638,
478 -1.10079476e+01, 5.09361474, -5.35062642e-01, -0.82084677,
479 1.65289503, -4.52520582, 2.05371873, -1.25665769,
480 -1.53605811, 4.41068581, 8.63583929, 2.19404252,
481 -9.97435470e-01, -2.72157537, 3.56591473e-01,-10.51631421,
482 -1.24085700,-10.61672181, 4.68232735e-01, -1.41031176,
483 8.69604433, 6.6481311 , 3.42348929, 9.84138181,
484 9.56312116e-01, 4.45183318, -4.04051175, -4.57584864,
485 2.18199500, 0.80126759, 3.44048454, -4.29726231,
486 -7.30042988, 5.44897893, 2.36553382, 2.10540332,
487 -2.98288336, -6.75607352, -2.77744882, -4.12872365,
488 1.15426258, 2.60704937, -4.41417255, -2.89610628,
489 2.81088838, -6.95136058, 1.25027939, -0.79932743,
490 1.45053876, -0.43730146, 1.63564303, 0.04638729,
491 2.06649835, 1.05026064, -2.06510578,-11.11320924,
492 -2.55106395, -6.5091449 , 2.46366137, -2.38704469,
493 -6.81565579, -2.16343352, -3.88047911, 2.43746762,
494 7.23650536e-01, 1.55500613, 1.81306385e-01, -1.1062498 ,
495 1.17932561, 3.50027177, -3.43094553, 6.00826981,
496 7.90979161, -5.27570134, -5.89365226, 3.18756318,
497 -4.32373034, -0.1337671 , 7.22740830e-01, 0.23983128,
498 -3.21535548, 1.45542238, -5.47938416, -4.01447941,
499 3.81058449, 0.76425931, 7.40985006, -0.44896188,
500 4.48609092, -1.64741879, -6.59974456,-13.17193115,
501 -5.59094391, -6.25244435, -5.13862996e-01, 0.15649305,
502 -5.98834455, 4.93129451, -4.61994107, 1.59820505,
503 3.81795409, 3.17146 , -4.34295854e-01, 8.2954126 ,
504 1.88220611, 0.80225125, -9.72981483e-01,11.394208 ,
505 -8.45662834e-01, -5.57411581, -1.97290542, -2.83361487,
506 -3.75283402, 2.47600433, 1.38662123, -3.98598717,
507 -4.63033506e-01, -3.21563035, -2.37481930, -1.42611867,
508 -1.17277194, -2.15404472, 3.42211139, -4.17689039,
509 2.24402195, -4.60132638, 6.78604790, 3.94154162,
510 3.01687661, -7.08195792, -2.95697930, 0.50367608,
511 -2.24107675, 5.52527559, 2.37498837, -3.8158066 ,
512 3.21442520, -2.14102714, -2.33990738, 1.27840434,
513 -6.94720428, -3.10502593, -1.32895472, 0.17674661,
514 6.17728674, -1.91667103, 6.13876866, 7.1040085 ,
515 1.24235316, -0.41309537, 6.79959099e-01, 1.88604949,
516 -3.69998656, -0.90050102, 3.59852152, 2.00887343,
517 -1.03422801e+01, 3.04462388, 2.12988453, 0.15820258,
518 3.33307819, 4.03094924, 9.19140042, -9.69149015,
519 2.59014379, -4.76329698, 6.71175474, -6.06191618,
520 -2.37625499, -7.42404308, -2.13423317, -4.86509946,
521 2.78518535, 2.6359501 , 3.75899776, 1.86910056,
522 -6.14851885,11.90893875, 4.49346648, 0.20428409,
523 -5.19226843, -2.5971315 , -2.22525608, 5.64794388,
524 8.82516491e-01, 3.53298728, 3.29866159e-01, -1.89217741,
525 2.35885332,10.04547561, 8.90553959e-01, 6.37113573,
526 -2.62744282, 5.20892662, -1.95661401, 6.34937202,
527 1.60008368, -1.63787578, -2.68466437, -1.37377726,
528 3.17748694, 4.45041455, 3.85748750, 2.93737277,
529 3.74616158, 2.59659457, -5.13129145, 7.26477374,
530 2.89711768, -1.42284586, -2.93995652, -0.39835926,
531 3.28959117, -4.7309582 , -2.61843435, -0.12691241,
532 -1.37722930e-01, 4.30291678, -1.91904127, -1.94645907,
533 -2.20163907, 1.62177814, 3.58389141, 0.14974575,
534 1.90769911, -0.56093423, 1.41015388, -4.51146401,
535 -3.08398468, -1.5177629 , 2.32709681, -0.68487206,
536 -7.83229243e-02, -1.81741166, -5.96239076, -1.99792128,
537 3.61799848, 5.10288284, -9.34809732, 0.58251846,
538 -1.00628421e+01, -2.87531563, 3.60644904e-01, 4.11521374,
539 -1.24855113, 1.43322867, 4.18500790, 0.14627552,
540 1.70963995, 7.67293252, -2.84779221e-01, -6.59958391,
541 2.93460322, -6.91076746, -2.56045800, 0.49557866,
542 -1.33416950e+01, -7.19088609, 6.54760504, -1.16188913,
543 -4.74493837, 0.35443496, -3.74306351, -9.63476301,
544 9.56297959, -0.24424306, -5.18407515e-01, -1.23622633,
545 -5.72520347, 3.94970165, -1.08927018e+01, -5.74801823,
546 -5.89005095e-01, 5.65483456, -4.11031641, 2.20538835,
547 -3.95269565, 2.32475678, -1.07329243e+01, 0.2719951 ,
548 5.04314682, 1.16928599, 5.16775420e-01, 2.90878759,
549 4.52179590, 7.08268703, -1.49214047e+01, -5.52200123,
550 5.84499843, 9.51454903, -2.91651323, 0.0258584 ,
551 5.10059531e-02, 3.4233929 , 3.23045185, -3.33383799,
552 -6.93782260, -2.4170692 , 3.15189786, 1.82518672,
553 -5.87506158e-01, 4.22833058, 4.46491682, -0.74625307,
554 6.28883709, -4.00634065, -7.73267254e-01, 3.72141078,
555 9.20921279e-01, 9.20510398, 1.24434715, -0.23621713,
556 1.21907020e+01, 6.47012928, -2.34462045e-02, 5.03789597,
557 5.99780199, -2.94900065, 8.15031483, -4.34599034,
558 5.82836675, -0.89396187, 1.72350562, 1.33641312,
559 4.76317833, -8.66518565, -6.16870680e-01, -2.43455407,
560 2.31174036, 3.52991041, 2.13693416, 2.20496619,
561 5.55562176, 4.71880166, 2.27203161, 4.61106145,
562 2.92973078, 3.75655459, 2.98128463, 0.83763449,
563 4.32617219, 4.33636417, 4.00943240, -7.76684892,
564 -1.80049237, -1.09902763, 2.55982143, 5.91349888
565 };
566
567 std::vector<std::vector<float>> expected_results_vectors{expected_result_zeros,
568 expected_result_ones,
569 expected_result_random};
570 arm::app::math::FftInstance fftInstance;
571 /* Iterate over each of the input vectors, calculate FFT and compare with corresponding expected_results vectors */
Kshitij Sisodiab59ba682021-11-23 17:19:52 +0000[diff] [blame]572 for (size_t j = 0; j < input_vectors.size(); j++) {
liabar01dee53bc2021-10-29 15:59:04 +0100[diff] [blame]573 uint16_t fftLen = input_vectors[j].size();
574 arm::app::math::MathUtils::FftInitF32(fftLen, fftInstance);
575 arm::app::math::MathUtils::FftF32(input_vectors[j], output_vectors[j], fftInstance);
576
577 float tolerance = 10e-4;
Kshitij Sisodiab59ba682021-11-23 17:19:52 +0000[diff] [blame]578 for (size_t i = 0; i < fftLen/2; i++) {
liabar01dee53bc2021-10-29 15:59:04 +0100[diff] [blame]579 CHECK (output_vectors[j][i] ==
580 Approx(expected_results_vectors[j][i]).margin(tolerance));
581
582 }
583 }
584}
585
586TEST_CASE("Test VecLogarithmF32")
587{
588 /*Test Constants: */
589 std::vector<float> input =
590 { 0.1e-10, 0.5, 1, M_PI, M_E };
591 std::vector<float> expectedResult =
592 {-25.328436, -0.693147181, 0, 1.144729886, 1};
593 std::vector<float> output(input.size());
594
595 arm::app::math::MathUtils::VecLogarithmF32(input,output);
596
597 for (size_t i = 0; i < input.size(); i++)
598 CHECK (expectedResult[i] == Approx(output[i]));
599}
600
601TEST_CASE("Test DotProductF32")
602{
603 /*Test Constants: */
604 std::vector<float> inputA
605 {1,1,1,0,0,0};
606 std::vector<float> inputB
607 {0,0,0,1,1,1};
608 uint32_t len = inputA.size();
609
610 float expectedResult = 0;
611 float dot_prod = arm::app::math::MathUtils::DotProductF32(inputA.data(), inputB.data(), len);
612 CHECK(dot_prod == expectedResult);
613}
614
615TEST_CASE("Test ComplexMagnitudeSquaredF32")
616{
617 /*Test Constants: */
618 std::vector<float> input
619 {0.0, 0.0, 0.5, 0.5,1,1};
Kshitij Sisodiab59ba682021-11-23 17:19:52 +0000[diff] [blame]620 size_t inputLen = input.size();
liabar01dee53bc2021-10-29 15:59:04 +0100[diff] [blame]621
622 std::vector<float> expectedResult
623 {0.0, 0.5, 2,};
Kshitij Sisodiab59ba682021-11-23 17:19:52 +0000[diff] [blame]624 size_t outputLen = inputLen/2;
liabar01dee53bc2021-10-29 15:59:04 +0100[diff] [blame]625 std::vector<float>output(outputLen);
626
627 /* Pass pointers to input/output vectors as this function over-writes the first half
628 * of the input vector with output results */
629 arm::app::math::MathUtils::ComplexMagnitudeSquaredF32(input.data(), inputLen, output.data(), outputLen);
630
Kshitij Sisodiab178b282022-01-04 13:37:53 +0000[diff] [blame]631 for (size_t i = 0; i < outputLen; i++) {
liabar01dee53bc2021-10-29 15:59:04 +0100[diff] [blame]632 CHECK (expectedResult[i] == Approx(output[i]));
Kshitij Sisodiab178b282022-01-04 13:37:53 +0000[diff] [blame]633 }
634}
635
636/**
637 * @brief Simple function to test the Softmax function
638 *
639 * @param input Input vector
640 * @param goldenOutput Expected output vector
641 */
642static void TestSoftmaxF32(
643 const std::vector<float>& input,
644 const std::vector<float>& goldenOutput)
645{
646 std::vector<float> output = input; /* Function modifies the vector in-place */
647 arm::app::math::MathUtils::SoftmaxF32(output);
648
649 for (size_t i = 0; i < goldenOutput.size(); ++i) {
650 CHECK(goldenOutput[i] == Approx(output[i]));
651 }
652
653 REQUIRE(output.size() == goldenOutput.size());
654}
655
656TEST_CASE("Test SoftmaxF32")
657{
658 SECTION("Simple series") {
659 const std::vector<float> input {
660 0.0, 1.0, 2.0, 3.0, 4.0,
661 5.0, 6.0, 7.0, 8.0, 9.0
662 };
663
664 const std::vector<float> expectedOutput {
665 7.80134161e-05, 2.12062451e-04,
666 5.76445508e-04, 1.56694135e-03,
667 4.25938820e-03, 1.15782175e-02,
668 3.14728583e-02, 8.55520989e-02,
669 2.32554716e-01, 6.32149258e-01
670 };
671
672 TestSoftmaxF32(input, expectedOutput);
673 }
674
675 SECTION("Random series") {
676 const std::vector<float> input {
677 0.8810943246170809, 0.5877587675947015,
678 0.6841546454788743, 0.4155920960071594,
679 0.9799415323651671, 0.5066432973545711,
680 0.3846024252355448, 0.4568689569632123,
681 0.3284413744557605, 0.49152323726213554
682 };
683
684 const std::vector<float> expectedOutput {
685 0.13329595, 0.09940837,
686 0.10946799, 0.08368583,
687 0.14714509, 0.09166319,
688 0.08113220, 0.08721240,
689 0.07670132, 0.09028766
690 };
691
692 TestSoftmaxF32(input, expectedOutput);
693 }
694
695 SECTION("Series with large STD") {
696 const std::vector<float> input {
697 0.001, 1000.000
698 };
699
700 const std::vector<float> expectedOutput {
701 0.000, 1.000
702 };
703
704 TestSoftmaxF32(input, expectedOutput);
705 }
liabar01dee53bc2021-10-29 15:59:04 +0100[diff] [blame]706}