45#define TEST_LENGTH_SAMPLES 480
46#define MAX_BLOCKSIZE 320
47#define MAX_NUMTAPS 100
49#define TEST_COUNT 5000
52static ne10_float32_t * guarded_in_c = NULL;
53static ne10_float32_t * guarded_in_neon = NULL;
54static ne10_float32_t * in_c = NULL;
55static ne10_float32_t * in_neon = NULL;
57static ne10_float32_t * guarded_out_c = NULL;
58static ne10_float32_t * guarded_out_neon = NULL;
59static ne10_float32_t * out_c = NULL;
60static ne10_float32_t * out_neon = NULL;
62static ne10_float32_t * guarded_fir_state_c = NULL;
63static ne10_float32_t * guarded_fir_state_neon = NULL;
64static ne10_float32_t * fir_state_c = NULL;
65static ne10_float32_t * fir_state_neon = NULL;
67#if defined (SMOKE_TEST)||(REGRESSION_TEST)
68static ne10_float32_t snr = 0.0f;
70#ifdef PERFORMANCE_TEST
71static ne10_int64_t time_c = 0;
72static ne10_int64_t time_neon = 0;
73static ne10_float32_t time_speedup = 0.0f;
74static ne10_float32_t time_savings = 0.0f;
81static ne10_float32_t testCoeffs32_f32[32] =
83 0.068186, 0.064344, -0.162450, 0.057015, 0.029743, 0.010066, 0.047792, 0.021273,
84 -0.096447, -0.211652, -0.086613, 0.057501, -0.187605, -0.167199, -0.026983, -0.025464,
85 -0.061495, 0.110914, -0.081973, -0.055231, -0.074430, -0.196536, 0.016845, -0.096493,
86 0.039625, -0.110273, -0.042966, -0.043804, 0.087350, -0.085191, 0.009420, 0.086440
93static ne10_float32_t testCoeffs8_f32[8] =
95 0.039625, -0.110273, -0.042966, -0.043804, 0.087350, -0.085191, 0.009420, 0.086440
107static ne10_float32_t testCoeffs27_f32[27] =
109 0.010066, 0.047792, 0.021273, -0.096447, -0.211652, -0.086613, 0.057501, -0.187605,
110 -0.167199, -0.026983, -0.025464, -0.061495, 0.110914, -0.081973, -0.055231, -0.074430,
111 -0.196536, 0.016845, -0.096493, 0.039625, -0.110273, -0.042966, -0.043804, 0.087350,
112 -0.085191, 0.009420, 0.086440
115static ne10_float32_t testCoeffs6_f32[6] =
117 -0.042966, -0.043804, 0.087350, -0.085191, 0.009420, 0.086440
125static ne10_float32_t testInput_f32[80] =
127 -0.432565, -1.665584, 0.125332, 0.287676, -1.146471, 1.190915, 1.189164, -0.037633,
128 0.327292, 0.174639, -0.186709, 0.725791, -0.588317, 2.183186, -0.136396, 0.113931,
129 1.066768, 0.059281, -0.095648, -0.832349, 0.294411, -1.336182, 0.714325, 1.623562,
130 -0.691776, 0.857997, 1.254001, -1.593730, -1.440964, 0.571148, -0.399886, 0.689997,
131 0.815622, 0.711908, 1.290250, 0.668601, 1.190838, -1.202457, -0.019790, -0.156717,
132 -1.604086, 0.257304, -1.056473, 1.415141, -0.805090, 0.528743, 0.219321, -0.921902,
133 -2.170674, -0.059188, -1.010634, 0.614463, 0.507741, 1.692430, 0.591283, -0.643595,
134 0.380337, -1.009116, -0.019511, -0.048221, 0.000043, -0.317859, 1.095004, -1.873990,
135 0.428183, 0.895638, 0.730957, 0.577857, 0.040314, 0.677089, 0.568900, -0.255645,
136 -0.377469, -0.295887, -1.475135, -0.234004, 0.118445, 0.314809, 1.443508, -0.350975
144 ne10_uint32_t blockSize;
145 ne10_uint32_t numTaps;
147 ne10_uint32_t numFrames;
148 ne10_float32_t *coeffsF32;
149 ne10_float32_t *inputF32;
153#if defined (SMOKE_TEST)||(REGRESSION_TEST)
154static test_config CONFIG[] = {{0, 1, 1, 10, &testCoeffs6_f32[0], &testInput_f32[0]},
155 {8, 6, 6, 10, &testCoeffs6_f32[0], &testInput_f32[0]},
156 {8, 8, 2, 10, &testCoeffs8_f32[0], &testInput_f32[0]},
157 {8, 27, 4, 10, &testCoeffs27_f32[0], &testInput_f32[0]},
158 {8, 32, 4, 10, &testCoeffs32_f32[0], &testInput_f32[0]},
159 {80, 6, 6, 1, &testCoeffs6_f32[0], &testInput_f32[0]},
160 {80, 8, 2, 1, &testCoeffs8_f32[0], &testInput_f32[0]},
161 {80, 27, 4, 1, &testCoeffs27_f32[0], &testInput_f32[0]},
162 {80, 32, 4, 1, &testCoeffs32_f32[0], &testInput_f32[0]}
164#define NUM_TESTS (sizeof(CONFIG) / sizeof(CONFIG[0]) )
166#ifdef PERFORMANCE_TEST
169 {8, 27, 3, 10, &testCoeffs27_f32[0], &testInput_f32[0]},
170 {8, 32, 4, 10, &testCoeffs32_f32[0], &testInput_f32[0]},
171 {80, 27, 3, 1, &testCoeffs27_f32[0], &testInput_f32[0]},
172 {80, 32, 4, 1, &testCoeffs32_f32[0], &testInput_f32[0]}
174#define NUM_PERF_TESTS (sizeof(CONFIG_PERF) / sizeof(CONFIG_PERF[0]) )
178void test_fir_interpolate_case0()
182 ne10_uint16_t loop = 0;
183 ne10_uint16_t block = 0;
185 ne10_uint16_t length = 0;
188 ne10_result_t status_c = NE10_OK;
189 ne10_result_t status_neon = NE10_OK;
191 fprintf (stdout,
"----------%30s start\n", __FUNCTION__);
194 NE10_SRC_ALLOC (in_c, guarded_in_c, TEST_LENGTH_SAMPLES);
195 NE10_SRC_ALLOC (in_neon, guarded_in_neon, TEST_LENGTH_SAMPLES);
198 NE10_DST_ALLOC (out_c, guarded_out_c, TEST_LENGTH_SAMPLES);
199 NE10_DST_ALLOC (out_neon, guarded_out_neon, TEST_LENGTH_SAMPLES);
202 NE10_DST_ALLOC (fir_state_c, guarded_fir_state_c, MAX_NUMTAPS + MAX_BLOCKSIZE);
203 NE10_DST_ALLOC (fir_state_neon, guarded_fir_state_neon, MAX_NUMTAPS + MAX_BLOCKSIZE);
205#ifdef ENABLE_NE10_FIR_INTERPOLATE_FLOAT_NEON
206#if defined (SMOKE_TEST)||(REGRESSION_TEST)
207 ne10_uint16_t pos = 0;
208 for (loop = 0; loop < NUM_TESTS; loop++)
210 config = &CONFIG[loop];
211 length = config->numFrames * config->blockSize * config->D;
214 status_c = ne10_fir_interpolate_init_float (&SC, config->D, config->numTaps, config->coeffsF32, fir_state_c, config->blockSize);
215 status_neon = ne10_fir_interpolate_init_float (&SN, config->D, config->numTaps, config->coeffsF32, fir_state_neon, config->blockSize);
217 if ( ( (status_c == NE10_ERR) || (status_neon == NE10_ERR)))
219 if (config->numTaps == 27)
221 fprintf (stdout,
"length of input data is wrong!\n");
226 fprintf (stdout,
"initialization error\n");
230 for (i = 0; i < 80; i++)
232 in_c[i] = testInput_f32[i];
233 in_neon[i] = testInput_f32[i];
236 GUARD_ARRAY (out_c, TEST_LENGTH_SAMPLES);
237 GUARD_ARRAY (out_neon, TEST_LENGTH_SAMPLES);
239 for (block = 0; block < config->numFrames; block++)
241 ne10_fir_interpolate_float_c (&SC, in_c + (block * config->blockSize), out_c + (block * config->blockSize * config->D), config->blockSize);
243 for (block = 0; block < config->numFrames; block++)
245 ne10_fir_interpolate_float_neon (&SN, in_neon + (block * config->blockSize), out_neon + (block * config->blockSize * config->D), config->blockSize);
248 CHECK_ARRAY_GUARD (out_c, TEST_LENGTH_SAMPLES);
249 CHECK_ARRAY_GUARD (out_neon, TEST_LENGTH_SAMPLES);
252 snr = CAL_SNR_FLOAT32 (out_c, out_neon, length);
253 assert_false ( (snr < SNR_THRESHOLD));
256#if defined (DEBUG_TRACE)
257 printf (
"--------------------config %d\n", loop);
258 printf (
"snr %f\n", snr);
260 for (pos = 0; pos < length; pos++)
262#if defined (DEBUG_TRACE)
263 printf (
"pos %d \n", pos);
264 printf (
"c %f (0x%04X) neon %f (0x%04X)\n", out_c[pos], * (ne10_uint32_t*) &out_c[pos], out_neon[pos], * (ne10_uint32_t*) &out_neon[pos]);
266 assert_float_vec_equal (&out_c[pos], &out_neon[pos], ERROR_MARGIN_SMALL, 1);
273#ifdef PERFORMANCE_TEST
275 fprintf (stdout,
"%25s%20s%20s%20s%20s\n",
"FIR Length&Taps",
"C Time in ms",
"NEON Time in ms",
"Time Savings",
"Performance Ratio");
276 for (loop = 0; loop < NUM_PERF_TESTS; loop++)
278 config = &CONFIG_PERF[loop];
279 length = config->numFrames * config->blockSize * config->D;
282 status_c = ne10_fir_interpolate_init_float (&SC, config->D, config->numTaps, config->coeffsF32, fir_state_c, config->blockSize);
283 status_neon = ne10_fir_interpolate_init_float (&SN, config->D, config->numTaps, config->coeffsF32, fir_state_neon, config->blockSize);
285 if ( ( (status_c == NE10_ERR) || (status_neon == NE10_ERR)))
287 if (config->numTaps == 27)
289 fprintf (stdout,
"length of input data is wrong!\n");
294 fprintf (stdout,
"initialization error\n");
299 for (i = 0; i < 80; i++)
301 in_c[i] = testInput_f32[i];
302 in_neon[i] = testInput_f32[i];
309 for (k = 0; k < TEST_COUNT; k++)
311 for (block = 0; block < config->numFrames; block++)
313 ne10_fir_interpolate_float_c (&SC, in_c + (block * config->blockSize), out_c + (block * config->blockSize * config->D), config->blockSize);
318#ifdef ENABLE_NE10_FIR_INTERPOLATE_FLOAT_NEON
323 for (k = 0; k < TEST_COUNT; k++)
325 for (block = 0; block < config->numFrames; block++)
327 ne10_fir_interpolate_float_neon (&SN, in_neon + (block * config->blockSize), out_neon + (block * config->blockSize * config->D), config->blockSize);
334 time_speedup = (ne10_float32_t) time_c / time_neon;
335 time_savings = ( ( (ne10_float32_t) (time_c - time_neon)) / time_c) * 100;
336 ne10_log (__FUNCTION__,
"%20d,%4d%20lld%20lld%19.2f%%%18.2f:1\n", config->numTaps, time_c, time_neon, time_savings, time_speedup);
341 free (guarded_in_neon);
342 free (guarded_out_c);
343 free (guarded_out_neon);
344 free (guarded_fir_state_c);
345 free (guarded_fir_state_neon);
346 fprintf (stdout,
"----------%30s end\n", __FUNCTION__);
349void test_fir_interpolate()
351 test_fir_interpolate_case0();
354static void my_test_setup (
void)
356 ne10_log_buffer_ptr = ne10_log_buffer;
359void test_fixture_fir_interpolate (
void)
361 test_fixture_start();
363 fixture_setup (my_test_setup);
365 run_test (test_fir_interpolate);
void ne10_fir_interpolate_float_c(const ne10_fir_interpolate_instance_f32_t *S, ne10_float32_t *pSrc, ne10_float32_t *pDst, ne10_uint32_t blockSize)
Processing function for the floating-point FIR interpolator.
Instance structure for the floating-point FIR Interpolation.