michaelh/zephyr
1
0
Fork 0
Code Releases Activity

tests: lib: cmsis_dsp: filtering: Add DECIM F32 tests for 1.9.0

Browse source 
This commit adds the filtering DECIM F32 test patterns and
implementations for the CMSIS-DSP 1.9.0.

Signed-off-by: Stephanos Ioannidis <root@stephanos.io>
This commit is contained in:
Stephanos Ioannidis 2021-08-20 21:26:32 +09:00 • committed by Carles Cufí
commit ee02e7ae64
4 changed files with 1703 additions and 0 deletions
Download patch file Download diff file Expand all files Collapse all files

1
tests/lib/cmsis_dsp/filtering/CMakeLists.txt
View file

@ -23,6 +23,7 @@ if(CONFIG_CMSIS_DSP_TEST_FILTERING_DECIM)
target_sources(app PRIVATE
src/decim_q15.c
src/decim_q31.c
src/decim_f32.c
)
endif()

154
tests/lib/cmsis_dsp/filtering/src/decim_f32.c Normal file
View file

@ -0,0 +1,154 @@
/*
* Copyright (c) 2021 Stephanos Ioannidis <root@stephanos.io>
* Copyright (C) 2010-2021 ARM Limited or its affiliates. All rights reserved.
*
* SPDX-License-Identifier: Apache-2.0
*/
#include <ztest.h>
#include <zephyr.h>
#include <stdlib.h>
#include <arm_math.h>
#include "../../common/test_common.h"
#include "decim_f32.pat"
#define SNR_ERROR_THRESH ((float32_t)100)
#define REL_ERROR_THRESH (8.0e-4)
#define STATE_BUF_LEN (16 + 768 - 1)
static void test_arm_fir_decimate_f32(void)
{
uint32_t decim_factor, tap_count;
size_t sample_index, block_size, ref_size;
size_t sample_count = ARRAY_SIZE(in_config_decim) / 4;
size_t length = ARRAY_SIZE(ref_decim);
const uint32_t *config = in_config_decim;
const float32_t *input = (const float32_t *)in_val_decim;
const float32_t *coeff = (const float32_t *)in_coeff_decim;
const float32_t *ref = (const float32_t *)ref_decim;
float32_t *state, *output_buf, *output;
arm_status status;
arm_fir_decimate_instance_f32 inst;
/* Allocate buffers */
state = malloc(STATE_BUF_LEN * sizeof(float32_t));
zassert_not_null(state, ASSERT_MSG_BUFFER_ALLOC_FAILED);
output_buf = malloc(length * sizeof(float32_t));
zassert_not_null(output_buf, ASSERT_MSG_BUFFER_ALLOC_FAILED);
output = output_buf;
/* Enumerate samples */
for (sample_index = 0; sample_index < sample_count; sample_index++) {
/* Resolve sample configurations */
decim_factor = config[0];
tap_count = config[1];
block_size = config[2];
ref_size = config[3];
/* Initialise instance */
status = arm_fir_decimate_init_f32(&inst, tap_count,
decim_factor, coeff, state,
block_size);
zassert_equal(status, ARM_MATH_SUCCESS,
ASSERT_MSG_INCORRECT_COMP_RESULT);
/* Run test function */
arm_fir_decimate_f32(&inst, input, output, block_size);
/* Increment pointers */
input += block_size;
output += ref_size;
coeff += tap_count;
config += 4;
}
/* Validate output */
zassert_true(
test_snr_error_f32(length, output_buf, ref, SNR_ERROR_THRESH),
ASSERT_MSG_SNR_LIMIT_EXCEED);
zassert_true(
test_rel_error_f32(length, output_buf, ref, REL_ERROR_THRESH),
ASSERT_MSG_REL_ERROR_LIMIT_EXCEED);
/* Free buffers */
free(state);
free(output_buf);
}
static void test_arm_fir_interpolate_f32(void)
{
uint32_t filter_length, tap_count;
size_t sample_index, block_size, ref_size;
size_t sample_count = ARRAY_SIZE(in_config_interp) / 4;
size_t length = ARRAY_SIZE(ref_interp);
const uint32_t *config = in_config_interp;
const float32_t *input = (const float32_t *)in_val_interp;
const float32_t *coeff = (const float32_t *)in_coeff_interp;
const float32_t *ref = (const float32_t *)ref_interp;
float32_t *state, *output_buf, *output;
arm_status status;
arm_fir_interpolate_instance_f32 inst;
/* Allocate buffers */
state = malloc(STATE_BUF_LEN * sizeof(float32_t));
zassert_not_null(state, ASSERT_MSG_BUFFER_ALLOC_FAILED);
output_buf = malloc(length * sizeof(float32_t));
zassert_not_null(output_buf, ASSERT_MSG_BUFFER_ALLOC_FAILED);
output = output_buf;
/* Enumerate samples */
for (sample_index = 0; sample_index < sample_count; sample_index++) {
/* Resolve sample configurations */
filter_length = config[0];
tap_count = config[1];
block_size = config[2];
ref_size = config[3];
/* Initialise instance */
status = arm_fir_interpolate_init_f32(&inst, filter_length,
tap_count, coeff,
state, block_size);
zassert_equal(status, ARM_MATH_SUCCESS,
ASSERT_MSG_INCORRECT_COMP_RESULT);
/* Run test function */
arm_fir_interpolate_f32(&inst, input, output, block_size);
/* Increment pointers */
input += block_size;
output += ref_size;
coeff += tap_count;
config += 4;
}
/* Validate output */
zassert_true(
test_snr_error_f32(length, output_buf, ref, SNR_ERROR_THRESH),
ASSERT_MSG_SNR_LIMIT_EXCEED);
zassert_true(
test_rel_error_f32(length, output_buf, ref, REL_ERROR_THRESH),
ASSERT_MSG_REL_ERROR_LIMIT_EXCEED);
/* Free buffers */
free(state);
free(output_buf);
}
void test_filtering_decim_f32(void)
{
ztest_test_suite(filtering_decim_f32,
ztest_unit_test(test_arm_fir_decimate_f32),
ztest_unit_test(test_arm_fir_interpolate_f32)
);
ztest_run_test_suite(filtering_decim_f32);
}

1546
tests/lib/cmsis_dsp/filtering/src/decim_f32.pat generated Normal file
View file

File diff suppressed because it is too large Load diff

2
tests/lib/cmsis_dsp/filtering/src/main.c
View file

@ -15,6 +15,7 @@ extern void test_filtering_biquad_f64(void);
extern void test_filtering_decim_q15(void);
extern void test_filtering_decim_q31(void);
extern void test_filtering_decim_f32(void);
extern void test_filtering_fir_q7(void);
extern void test_filtering_fir_q15(void);
@ -43,6 +44,7 @@ void test_main(void)
#ifdef CONFIG_CMSIS_DSP_TEST_FILTERING_DECIM
test_filtering_decim_q15();
test_filtering_decim_q31();
test_filtering_decim_f32();
#endif
#ifdef CONFIG_CMSIS_DSP_TEST_FILTERING_FIR