1 // SPDX-License-Identifier: BSD-3-Clause
2 //
3 // Copyright(c) 2021 Intel Corporation. All rights reserved.
4 //
5 // Author: Shriram Shastry <malladi.sastry@linux.intel.com>
6
7 #include <stdio.h>
8 #include <stdint.h>
9 #include <stdarg.h>
10 #include <stddef.h>
11 #include <setjmp.h>
12 #include <math.h>
13 #include <cmocka.h>
14 #include <string.h>
15
16 #include <sof/math/log.h>
17 #include <sof/audio/format.h>
18 #include <rtos/string.h>
19 #include <sof/common.h>
20 #include "log2_tables.h"
21 /* 'Error[max] = 0.0000236785999981,THD(-dBc) = -92.5128795787487235' */
22 #define CMP_TOLERANCE 0.0000236785691029
23
24 /* testvector in Q32.0 */
25 static const uint32_t uv[100] = {
26 1ULL, 43383509ULL, 86767017ULL, 130150525ULL, 173534033ULL, 216917541ULL,
27 260301049ULL, 303684557ULL, 347068065ULL, 390451573ULL, 433835081ULL, 477218589ULL,
28 520602097ULL, 563985605ULL, 607369113ULL, 650752621ULL, 694136129ULL, 737519638ULL,
29 780903146ULL, 824286654ULL, 867670162ULL, 911053670ULL, 954437178ULL, 997820686ULL,
30 1041204194ULL, 1084587702ULL, 1127971210ULL, 1171354718ULL, 1214738226ULL, 1258121734ULL,
31 1301505242ULL, 1344888750ULL, 1388272258ULL, 1431655766ULL, 1475039274ULL, 1518422782ULL,
32 1561806290ULL, 1605189798ULL, 1648573306ULL, 1691956814ULL, 1735340322ULL, 1778723830ULL,
33 1822107338ULL, 1865490846ULL, 1908874354ULL, 1952257862ULL, 1995641370ULL, 2039024878ULL,
34 2082408386ULL, 2125791894ULL, 2169175403ULL, 2212558911ULL, 2255942419ULL, 2299325927ULL,
35 2342709435ULL, 2386092943ULL, 2429476451ULL, 2472859959ULL, 2516243467ULL, 2559626975ULL,
36 2603010483ULL, 2646393991ULL, 2689777499ULL, 2733161007ULL, 2776544515ULL, 2819928023ULL,
37 2863311531ULL, 2906695039ULL, 2950078547ULL, 2993462055ULL, 3036845563ULL, 3080229071ULL,
38 3123612579ULL, 3166996087ULL, 3210379595ULL, 3253763103ULL, 3297146611ULL, 3340530119ULL,
39 3383913627ULL, 3427297135ULL, 3470680643ULL, 3514064151ULL, 3557447659ULL, 3600831168ULL,
40 3644214676ULL, 3687598184ULL, 3730981692ULL, 3774365200ULL, 3817748708ULL, 3861132216ULL,
41 3904515724ULL, 3947899232ULL, 3991282740ULL, 4034666248ULL, 4078049756ULL, 4121433264ULL,
42 4164816772ULL, 4208200280ULL, 4251583788ULL, 4294967295ULL};
43
test_math_arithmetic_base2log_fixed(void ** state)44 static void test_math_arithmetic_base2log_fixed(void **state)
45 {
46 (void)state;
47
48 uint32_t u[100];
49 int i;
50
51 memcpy_s((void *)&u[0], sizeof(u), (void *)&uv[0], 100U * sizeof(uint32_t));
52 for (i = 0; i < ARRAY_SIZE(log2_lookup_table); i++) {
53 float y = Q_CONVERT_QTOF(base2_logarithm(u[i]), 0);
54 float diff = fabs(log2_lookup_table[i] - (double)y / (1 << 16));
55
56 if (diff > CMP_TOLERANCE) {
57 printf("%s: diff for %.16f: value = %.16f, log2 = %.16f\n", __func__, diff,
58 (double)u[i], (double)y / (1 << 16));
59 assert_true(diff <= CMP_TOLERANCE);
60 }
61 }
62 }
63
main(void)64 int main(void)
65 {
66 const struct CMUnitTest tests[] = {
67 cmocka_unit_test(test_math_arithmetic_base2log_fixed)
68 };
69
70 cmocka_set_message_output(CM_OUTPUT_TAP);
71
72 return cmocka_run_group_tests(tests, NULL, NULL);
73 }
74
