Supporters update.
[dcpomatic.git] / test / audio_filter_test.cc
1 /*
2     Copyright (C) 2014-2021 Carl Hetherington <cth@carlh.net>
3
4     This file is part of DCP-o-matic.
5
6     DCP-o-matic is free software; you can redistribute it and/or modify
7     it under the terms of the GNU General Public License as published by
8     the Free Software Foundation; either version 2 of the License, or
9     (at your option) any later version.
10
11     DCP-o-matic is distributed in the hope that it will be useful,
12     but WITHOUT ANY WARRANTY; without even the implied warranty of
13     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
14     GNU General Public License for more details.
15
16     You should have received a copy of the GNU General Public License
17     along with DCP-o-matic.  If not, see <http://www.gnu.org/licenses/>.
18
19 */
20
21
22 /** @file  test/audio_filter_test.cc
23  *  @brief Test AudioFilter, LowPassAudioFilter, HighPassAudioFilter classes.
24  *  @ingroup selfcontained
25  */
26
27
28 #include <boost/test/unit_test.hpp>
29 #include "lib/audio_filter.h"
30 #include "lib/audio_buffers.h"
31
32
33 using std::make_shared;
34
35
36 static void
37 audio_filter_impulse_test_one (AudioFilter& f, int block_size, int num_blocks)
38 {
39         int c = 0;
40
41         for (int i = 0; i < num_blocks; ++i) {
42
43                 auto in = make_shared<AudioBuffers>(1, block_size);
44                 for (int j = 0; j < block_size; ++j) {
45                         in->data()[0][j] = c + j;
46                 }
47
48                 auto out = f.run (in);
49
50                 for (int j = 0; j < out->frames(); ++j) {
51                         BOOST_CHECK_EQUAL (out->data()[0][j], c + j);
52                 }
53
54                 c += block_size;
55         }
56 }
57
58
59 /** Create a filter with an impulse as a kernel and check that it
60  *  passes data through unaltered.
61  */
62 BOOST_AUTO_TEST_CASE (audio_filter_impulse_kernel_test)
63 {
64         AudioFilter f (0.02);
65
66         f._ir.resize(f._M + 1);
67         f._ir[0] = 1;
68         for (int i = 1; i <= f._M; ++i) {
69                 f._ir[i] = 0;
70         }
71
72         audio_filter_impulse_test_one (f, 32, 1);
73         audio_filter_impulse_test_one (f, 256, 1);
74         audio_filter_impulse_test_one (f, 2048, 1);
75 }
76
77
78 /** Create filters and pass them impulses as input and check that
79  *  the filter kernels comes back.
80  */
81 BOOST_AUTO_TEST_CASE (audio_filter_impulse_input_test)
82 {
83         LowPassAudioFilter lpf (0.02, 0.3);
84
85         auto in = make_shared<AudioBuffers>(1, 1751);
86         in->make_silent ();
87         in->data(0)[0] = 1;
88
89         auto out = lpf.run (in);
90         for (int j = 0; j < out->frames(); ++j) {
91                 if (j <= lpf._M) {
92                         BOOST_CHECK_EQUAL (out->data(0)[j], lpf._ir[j]);
93                 } else {
94                         BOOST_CHECK_EQUAL (out->data(0)[j], 0);
95                 }
96         }
97
98         HighPassAudioFilter hpf (0.02, 0.3);
99
100         in = make_shared<AudioBuffers>(1, 9133);
101         in->make_silent ();
102         in->data(0)[0] = 1;
103
104         out = hpf.run (in);
105         for (int j = 0; j < out->frames(); ++j) {
106                 if (j <= hpf._M) {
107                         BOOST_CHECK_EQUAL (out->data(0)[j], hpf._ir[j]);
108                 } else {
109                         BOOST_CHECK_EQUAL (out->data(0)[j], 0);
110                 }
111         }
112 }