c@225
|
1 /* -*- c-basic-offset: 4 indent-tabs-mode: nil -*- vi:set ts=8 sts=4 sw=4: */
|
c@225
|
2 /*
|
c@225
|
3 QM DSP Library
|
c@225
|
4
|
c@225
|
5 Centre for Digital Music, Queen Mary, University of London.
|
c@309
|
6 This file 2005-2006 Christian Landone.
|
c@309
|
7
|
c@309
|
8 This program is free software; you can redistribute it and/or
|
c@309
|
9 modify it under the terms of the GNU General Public License as
|
c@309
|
10 published by the Free Software Foundation; either version 2 of the
|
c@309
|
11 License, or (at your option) any later version. See the file
|
c@309
|
12 COPYING included with this distribution for more information.
|
c@225
|
13 */
|
c@225
|
14
|
c@225
|
15 #ifndef DECIMATOR_H
|
c@225
|
16 #define DECIMATOR_H
|
c@225
|
17
|
c@375
|
18 /**
|
c@375
|
19 * Decimator carries out a fast downsample by a power-of-two
|
c@375
|
20 * factor. Only a limited number of factors are supported, from two to
|
c@375
|
21 * whatever getHighestSupportedFactor() returns. This is much faster
|
c@375
|
22 * than Resampler but has a worse signal-noise ratio.
|
c@375
|
23 */
|
c@225
|
24 class Decimator
|
c@225
|
25 {
|
c@225
|
26 public:
|
c@247
|
27 void process( const double* src, double* dst );
|
c@280
|
28 void process( const float* src, float* dst );
|
c@225
|
29
|
c@279
|
30 /**
|
c@279
|
31 * Construct a Decimator to operate on input blocks of length
|
c@279
|
32 * inLength, with decimation factor decFactor. inLength should be
|
c@279
|
33 * a multiple of decFactor. Output blocks will be of length
|
c@279
|
34 * inLength / decFactor.
|
c@279
|
35 *
|
c@279
|
36 * decFactor must be a power of two. The highest supported factor
|
c@279
|
37 * is obtained through getHighestSupportedFactor(); for higher
|
c@279
|
38 * factors, you will need to chain more than one decimator.
|
c@279
|
39 */
|
c@225
|
40 Decimator( unsigned int inLength, unsigned int decFactor );
|
c@225
|
41 virtual ~Decimator();
|
c@225
|
42
|
c@247
|
43 int getFactor() const { return m_decFactor; }
|
c@247
|
44 static int getHighestSupportedFactor() { return 8; }
|
c@247
|
45
|
c@225
|
46 private:
|
c@225
|
47 void resetFilter();
|
c@225
|
48 void deInitialise();
|
c@225
|
49 void initialise( unsigned int inLength, unsigned int decFactor );
|
c@280
|
50 void doAntiAlias( const double* src, double* dst, unsigned int length );
|
c@280
|
51 void doAntiAlias( const float* src, double* dst, unsigned int length );
|
c@225
|
52
|
c@225
|
53 unsigned int m_inputLength;
|
c@225
|
54 unsigned int m_outputLength;
|
c@225
|
55 unsigned int m_decFactor;
|
c@225
|
56
|
c@225
|
57 double Input;
|
c@225
|
58 double Output ;
|
c@225
|
59
|
c@225
|
60 double o1,o2,o3,o4,o5,o6,o7;
|
c@225
|
61
|
c@225
|
62 double a[ 9 ];
|
c@225
|
63 double b[ 9 ];
|
c@225
|
64
|
c@225
|
65 double* decBuffer;
|
c@225
|
66 };
|
c@225
|
67
|
c@225
|
68 #endif //
|