1 #include "MFCCF16.h"
2 #include <stdio.h>
3 #include "Error.h"
4 
5 #include "mfccdata_f16.h"
6 
7 #define SNR_THRESHOLD 45
8 
9 /*
10 
11 Reference patterns are generated with
12 a double precision computation.
13 
14 */
15 #define REL_ERROR (2.0e-2)
16 #define ABS_ERROR (2.0e-2)
17 
18 
test_mfcc_f16()19     void MFCCF16::test_mfcc_f16()
20     {
21         const float16_t *inp1=input1.ptr();
22         float16_t *tmpinp=tmpin.ptr();
23         float16_t *outp=output.ptr();
24         float16_t *tmpp=tmp.ptr();
25 
26 
27         memcpy((void*)tmpinp,(void*)inp1,sizeof(float16_t)*this->fftLen);
28         arm_mfcc_f16(&mfcc,tmpinp,outp,tmpp);
29 
30         ASSERT_EMPTY_TAIL(output);
31 
32         ASSERT_SNR(output,ref,(float16_t)SNR_THRESHOLD);
33 
34         ASSERT_CLOSE_ERROR(output,ref,ABS_ERROR,REL_ERROR);
35 
36     }
37 
38 
setUp(Testing::testID_t id,std::vector<Testing::param_t> & params,Client::PatternMgr * mgr)39     void MFCCF16::setUp(Testing::testID_t id,std::vector<Testing::param_t>& params,Client::PatternMgr *mgr)
40     {
41 
42        (void)params;
43 
44        Testing::nbSamples_t nb=MAX_NB_SAMPLES;
45 
46 
47        switch(id)
48        {
49         case MFCCF16::TEST_MFCC_F16_1:
50         {
51             nb = 256;
52             this->fftLen = nb;
53             ref.reload(MFCCF16::REF_MFCC_NOISE_256_F16_ID,mgr,nb);
54             input1.reload(MFCCF16::INPUTS_MFCC_NOISE_256_F16_ID,mgr,nb);
55             arm_mfcc_init_f16(&mfcc,
56                     nb,20,13,mfcc_dct_coefs_config1_f16,
57                     mfcc_filter_pos_config3_f16,mfcc_filter_len_config3_f16,
58                     mfcc_filter_coefs_config3_f16,
59                     mfcc_window_coefs_config3_f16);
60             tmp.create(2*nb,MFCCF16::TMP_MFCC_F16_ID,mgr);
61             tmpin.create(nb,MFCCF16::TMPIN_MFCC_F16_ID,mgr);
62           }
63           break;
64 
65         case MFCCF16::TEST_MFCC_F16_2:
66           {
67             nb = 512;
68             this->fftLen = nb;
69             ref.reload(MFCCF16::REF_MFCC_NOISE_512_F16_ID,mgr,nb);
70             input1.reload(MFCCF16::INPUTS_MFCC_NOISE_512_F16_ID,mgr,nb);
71             arm_mfcc_init_f16(&mfcc,
72                       nb,20,13,mfcc_dct_coefs_config1_f16,
73                       mfcc_filter_pos_config2_f16,mfcc_filter_len_config2_f16,
74                       mfcc_filter_coefs_config2_f16,
75                       mfcc_window_coefs_config2_f16);
76             tmp.create(2*nb,MFCCF16::TMP_MFCC_F16_ID,mgr);
77             tmpin.create(nb,MFCCF16::TMPIN_MFCC_F16_ID,mgr);
78           }
79           break;
80         case MFCCF16::TEST_MFCC_F16_3:
81           {
82             nb = 1024;
83             this->fftLen = nb;
84             ref.reload(MFCCF16::REF_MFCC_NOISE_1024_F16_ID,mgr,nb);
85             input1.reload(MFCCF16::INPUTS_MFCC_NOISE_1024_F16_ID,mgr,nb);
86             arm_mfcc_init_f16(&mfcc,
87                       nb,20,13,mfcc_dct_coefs_config1_f16,
88                       mfcc_filter_pos_config1_f16,mfcc_filter_len_config1_f16,
89                       mfcc_filter_coefs_config1_f16,
90                       mfcc_window_coefs_config1_f16);
91             tmp.create(2*nb,MFCCF16::TMP_MFCC_F16_ID,mgr);
92             tmpin.create(nb,MFCCF16::TMPIN_MFCC_F16_ID,mgr);
93 
94           }
95           break;
96 
97         case MFCCF16::TEST_MFCC_F16_4:
98         {
99             nb = 256;
100             this->fftLen = nb;
101             ref.reload(MFCCF16::REF_MFCC_SINE_256_F16_ID,mgr,nb);
102             input1.reload(MFCCF16::INPUTS_MFCC_SINE_256_F16_ID,mgr,nb);
103             arm_mfcc_init_f16(&mfcc,
104                     nb,20,13,mfcc_dct_coefs_config1_f16,
105                     mfcc_filter_pos_config3_f16,mfcc_filter_len_config3_f16,
106                     mfcc_filter_coefs_config3_f16,
107                     mfcc_window_coefs_config3_f16);
108             tmp.create(2*nb,MFCCF16::TMP_MFCC_F16_ID,mgr);
109             tmpin.create(nb,MFCCF16::TMPIN_MFCC_F16_ID,mgr);
110           }
111           break;
112 
113         case MFCCF16::TEST_MFCC_F16_5:
114           {
115             nb = 512;
116             this->fftLen = nb;
117             ref.reload(MFCCF16::REF_MFCC_SINE_512_F16_ID,mgr,nb);
118             input1.reload(MFCCF16::INPUTS_MFCC_SINE_512_F16_ID,mgr,nb);
119             arm_mfcc_init_f16(&mfcc,
120                       nb,20,13,mfcc_dct_coefs_config1_f16,
121                       mfcc_filter_pos_config2_f16,mfcc_filter_len_config2_f16,
122                       mfcc_filter_coefs_config2_f16,
123                       mfcc_window_coefs_config2_f16);
124             tmp.create(2*nb,MFCCF16::TMP_MFCC_F16_ID,mgr);
125             tmpin.create(nb,MFCCF16::TMPIN_MFCC_F16_ID,mgr);
126           }
127           break;
128         case MFCCF16::TEST_MFCC_F16_6:
129           {
130             nb = 1024;
131             this->fftLen = nb;
132             ref.reload(MFCCF16::REF_MFCC_SINE_1024_F16_ID,mgr,nb);
133             input1.reload(MFCCF16::INPUTS_MFCC_SINE_1024_F16_ID,mgr,nb);
134             arm_mfcc_init_f16(&mfcc,
135                       nb,20,13,mfcc_dct_coefs_config1_f16,
136                       mfcc_filter_pos_config1_f16,mfcc_filter_len_config1_f16,
137                       mfcc_filter_coefs_config1_f16,
138                       mfcc_window_coefs_config1_f16);
139             tmp.create(2*nb,MFCCF16::TMP_MFCC_F16_ID,mgr);
140             tmpin.create(nb,MFCCF16::TMPIN_MFCC_F16_ID,mgr);
141 
142           }
143           break;
144 
145        }
146 
147 
148        output.create(ref.nbSamples(),MFCCF16::OUTPUT_MFCC_F16_ID,mgr);
149 
150     }
151 
tearDown(Testing::testID_t id,Client::PatternMgr * mgr)152     void MFCCF16::tearDown(Testing::testID_t id,Client::PatternMgr *mgr)
153     {
154         (void)id;
155         (void)mgr;
156         //output.dump(mgr);
157     }
158