3 * Copyright 2010 Free Software Foundation, Inc.
5 * This file is part of GNU Radio
7 * GNU Radio is free software; you can redistribute it and/or modify
8 * it under the terms of the GNU General Public License as published by
9 * the Free Software Foundation; either version 3, or (at your option)
12 * GNU Radio is distributed in the hope that it will be useful,
13 * but WITHOUT ANY WARRANTY; without even the implied warranty of
14 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
15 * GNU General Public License for more details.
17 * You should have received a copy of the GNU General Public License
18 * along with GNU Radio; see the file COPYING. If not, write to
19 * the Free Software Foundation, Inc., 51 Franklin Street,
20 * Boston, MA 02110-1301, USA.
28 #include <qa_gri_fir_filter_with_buffer_ccc.h>
29 #include <gri_fir_filter_with_buffer_ccc.h>
33 #include <cppunit/TestAssert.h>
38 typedef gr_complex i_type;
39 typedef gr_complex o_type;
40 typedef gr_complex tap_type;
41 typedef gr_complex acc_type;
45 #define ERR_DELTA (1e-5)
47 #define NELEM(x) (sizeof (x) / sizeof (x[0]))
52 return 2.0 * ((float) random () / RANDOM_MAX - 0.5); // uniformly (-1, 1)
56 random_complex (gr_complex *buf, unsigned n)
58 for (unsigned i = 0; i < n; i++){
59 float re = rint (uniform () * 32767);
60 float im = rint (uniform () * 32767);
61 buf[i] = gr_complex (re, im);
66 ref_dotprod (const i_type input[], const tap_type taps[], int ntaps)
69 for (int i = 0; i < ntaps; i++) {
70 sum += input[i] * taps[i];
77 qa_gri_fir_filter_with_buffer_ccc::t1 ()
83 qa_gri_fir_filter_with_buffer_ccc::t2 ()
89 qa_gri_fir_filter_with_buffer_ccc::t3 ()
95 // Test for ntaps in [0,9], and input lengths in [0,17].
96 // This ensures that we are building the shifted taps correctly,
97 // and exercises all corner cases on input alignment and length.
100 qa_gri_fir_filter_with_buffer_ccc::test_decimate(unsigned int decimate)
102 const int MAX_TAPS = 9;
103 const int OUTPUT_LEN = 17;
104 const int INPUT_LEN = MAX_TAPS + OUTPUT_LEN;
106 // Mem aligned buffer not really necessary, but why not?
107 i_type *input = (i_type *)malloc16Align(INPUT_LEN * sizeof(i_type));
108 i_type *dline = (i_type*)malloc16Align(INPUT_LEN * sizeof(i_type));
109 o_type expected_output[OUTPUT_LEN];
110 o_type actual_output[OUTPUT_LEN];
111 tap_type taps[MAX_TAPS];
113 srandom (0); // we want reproducibility
114 memset(dline, 0, INPUT_LEN*sizeof(i_type));
116 for (int n = 0; n <= MAX_TAPS; n++){
117 for (int ol = 0; ol <= OUTPUT_LEN; ol++){
119 // cerr << "@@@ n:ol " << n << ":" << ol << endl;
121 // build random test case
122 random_complex (input, INPUT_LEN);
123 random_complex (taps, MAX_TAPS);
125 // compute expected output values
126 memset(dline, 0, INPUT_LEN*sizeof(i_type));
127 for (int o = 0; o < (int)(ol/decimate); o++){
128 // use an actual delay line for this test
129 for(int dd = 0; dd < (int)decimate; dd++) {
130 for(int oo = INPUT_LEN-1; oo > 0; oo--)
131 dline[oo] = dline[oo-1];
132 dline[0] = input[decimate*o+dd];
134 expected_output[o] = ref_dotprod (dline, taps, n);
138 vector<tap_type> f1_taps(&taps[0], &taps[n]);
139 gri_fir_filter_with_buffer_ccc *f1 = new gri_fir_filter_with_buffer_ccc(f1_taps);
141 // zero the output, then do the filtering
142 memset (actual_output, 0, sizeof (actual_output));
143 f1->filterNdec (actual_output, input, ol/decimate, decimate);
147 // we use a sloppy error margin because on the x86 architecture,
148 // our reference implementation is using 80 bit floating point
149 // arithmetic, while the SSE version is using 32 bit float point
152 for (int o = 0; o < (int)(ol/decimate); o++){
153 CPPUNIT_ASSERT_COMPLEXES_EQUAL(expected_output[o], actual_output[o],
154 abs (expected_output[o]) * ERR_DELTA);