1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
|
/*
* Copyright (C) 2018-2020 Mindmaze Holdings SA
*
* This file is part of the rtfilter library
*
* The rtfilter library is free software: you can redistribute it and/or
* modify it under the terms of the version 3 of the GNU Lesser General
* Public License as published by the Free Software Foundation.
*
* This program is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU Lesser General Public License for more details.
*
* You should have received a copy of the GNU Lesser General Public License
* along with this program. If not, see <http://www.gnu.org/licenses/>.
*/
#include <check.h>
#include <math.h>
#include <rtfilter.h>
#include <stdlib.h>
#include "float-comparison.h"
#include "rtfilter.h"
#include "testcases.h"
#include "test-common.h"
#define NCH 4
#define NS 128
#define FS 2000
static
struct rtf_filter* create_lowpass(int nch, double fs, int order, double lowcut,
int data_type)
{
double fc;
struct rtf_filter * f;
ck_assert(nch > 0);
ck_assert(lowcut > 0.);
ck_assert(fs > lowcut);
ck_assert(order > 0);
fc = lowcut / fs;
f = rtf_create_butterworth(nch, data_type, fc, order, 0);
ck_assert(f != NULL);
return f;
}
START_TEST(test_lowpass_coeffs)
{
struct rtf_filter * f;
f = create_lowpass(NCH, FS, 2, 50., RTF_FLOAT);
ck_assert(f != NULL);
struct rtf_coeffs * _coeffs = rtf_get_coeffs(f);
struct real_coeffs * coeffs = &_coeffs->coeffs.real;
/* compare generated coeffs against values computed by matlab */
ck_assert(!_coeffs->is_complex);
/* test generated num */
ck_assert_int_eq(coeffs->num_len, 3);
ck_assert(compare_floats(coeffs->num[0], 0.005542717408) == 0);
ck_assert(compare_floats(coeffs->num[1], 0.011085434817) == 0);
ck_assert(compare_floats(coeffs->num[2], 0.0055427174084) == 0);
/* test generated denum */
ck_assert(coeffs->denum[0] == 1.); /* exact value */
ck_assert_int_eq(coeffs->denum_len, 3);
ck_assert(compare_floats(coeffs->denum[1], -1.778631806374) == 0);
ck_assert(compare_floats(coeffs->denum[2], 0.8008026480675) == 0);
rtf_destroy_filter(f);
}
END_TEST
START_TEST(test_lowpass_smoke)
{
struct rtf_filter * f;
f = create_lowpass(NCH, FS, 2, 50., RTF_FLOAT);
ck_assert(f != NULL);
smoke_test(f, NCH, NS);
rtf_destroy_filter(f);
}
END_TEST
START_TEST(test_lowpass_flat)
{
struct rtf_filter * f;
f = create_lowpass(NCH, FS, 2, 50., RTF_FLOAT);
ck_assert(f != NULL);
flat_test(f, NCH, NS);
rtf_destroy_filter(f);
}
END_TEST
TCase* create_lowpass_tcase(void)
{
TCase * tc;
tc = tcase_create("lowpass");
tcase_add_test(tc, test_lowpass_coeffs);
tcase_add_test(tc, test_lowpass_flat);
tcase_add_test(tc, test_lowpass_smoke);
return tc;
}
|
