1 | /*****************************************************************************
|
---|
2 |
|
---|
3 | FFTRealFixLen.h
|
---|
4 | Copyright (c) 2005 Laurent de Soras
|
---|
5 |
|
---|
6 | --- Legal stuff ---
|
---|
7 |
|
---|
8 | This library is free software; you can redistribute it and/or
|
---|
9 | modify it under the terms of the GNU Lesser General Public
|
---|
10 | License as published by the Free Software Foundation; either
|
---|
11 | version 2.1 of the License, or (at your option) any later version.
|
---|
12 |
|
---|
13 | This library is distributed in the hope that it will be useful,
|
---|
14 | but WITHOUT ANY WARRANTY; without even the implied warranty of
|
---|
15 | MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
|
---|
16 | Lesser General Public License for more details.
|
---|
17 |
|
---|
18 | You should have received a copy of the GNU Lesser General Public
|
---|
19 | License along with this library; if not, write to the Free Software
|
---|
20 | Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
|
---|
21 |
|
---|
22 | *Tab=3***********************************************************************/
|
---|
23 |
|
---|
24 |
|
---|
25 |
|
---|
26 | #if ! defined (FFTRealFixLen_HEADER_INCLUDED)
|
---|
27 | #define FFTRealFixLen_HEADER_INCLUDED
|
---|
28 |
|
---|
29 | #if defined (_MSC_VER)
|
---|
30 | #pragma once
|
---|
31 | #pragma warning (4 : 4250) // "Inherits via dominance."
|
---|
32 | #endif
|
---|
33 |
|
---|
34 |
|
---|
35 |
|
---|
36 | /*\\\ INCLUDE FILES \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\*/
|
---|
37 |
|
---|
38 | #include "Array.h"
|
---|
39 | #include "DynArray.h"
|
---|
40 | #include "FFTRealFixLenParam.h"
|
---|
41 | #include "OscSinCos.h"
|
---|
42 |
|
---|
43 |
|
---|
44 |
|
---|
45 | template <int LL2>
|
---|
46 | class FFTRealFixLen
|
---|
47 | {
|
---|
48 | typedef int CompileTimeCheck1 [(LL2 >= 0) ? 1 : -1];
|
---|
49 | typedef int CompileTimeCheck2 [(LL2 <= 30) ? 1 : -1];
|
---|
50 |
|
---|
51 | /*\\\ PUBLIC \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\*/
|
---|
52 |
|
---|
53 | public:
|
---|
54 |
|
---|
55 | typedef FFTRealFixLenParam::DataType DataType;
|
---|
56 | typedef OscSinCos <DataType> OscType;
|
---|
57 |
|
---|
58 | enum { FFT_LEN_L2 = LL2 };
|
---|
59 | enum { FFT_LEN = 1 << FFT_LEN_L2 };
|
---|
60 |
|
---|
61 | FFTRealFixLen ();
|
---|
62 |
|
---|
63 | inline long get_length () const;
|
---|
64 | void do_fft (DataType f [], const DataType x []);
|
---|
65 | void do_ifft (const DataType f [], DataType x []);
|
---|
66 | void rescale (DataType x []) const;
|
---|
67 |
|
---|
68 |
|
---|
69 |
|
---|
70 | /*\\\ PROTECTED \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\*/
|
---|
71 |
|
---|
72 | protected:
|
---|
73 |
|
---|
74 |
|
---|
75 |
|
---|
76 | /*\\\ PRIVATE \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\*/
|
---|
77 |
|
---|
78 | private:
|
---|
79 |
|
---|
80 | enum { TRIGO_BD_LIMIT = FFTRealFixLenParam::TRIGO_BD_LIMIT };
|
---|
81 |
|
---|
82 | enum { BR_ARR_SIZE_L2 = ((FFT_LEN_L2 - 3) < 0) ? 0 : (FFT_LEN_L2 - 2) };
|
---|
83 | enum { BR_ARR_SIZE = 1 << BR_ARR_SIZE_L2 };
|
---|
84 |
|
---|
85 | enum { TRIGO_BD = ((FFT_LEN_L2 - TRIGO_BD_LIMIT) < 0)
|
---|
86 | ? (int)FFT_LEN_L2
|
---|
87 | : (int)TRIGO_BD_LIMIT };
|
---|
88 | enum { TRIGO_TABLE_ARR_SIZE_L2 = (LL2 < 4) ? 0 : (TRIGO_BD - 2) };
|
---|
89 | enum { TRIGO_TABLE_ARR_SIZE = 1 << TRIGO_TABLE_ARR_SIZE_L2 };
|
---|
90 |
|
---|
91 | enum { NBR_TRIGO_OSC = FFT_LEN_L2 - TRIGO_BD };
|
---|
92 | enum { TRIGO_OSC_ARR_SIZE = (NBR_TRIGO_OSC > 0) ? NBR_TRIGO_OSC : 1 };
|
---|
93 |
|
---|
94 | void build_br_lut ();
|
---|
95 | void build_trigo_lut ();
|
---|
96 | void build_trigo_osc ();
|
---|
97 |
|
---|
98 | DynArray <DataType>
|
---|
99 | _buffer;
|
---|
100 | DynArray <long>
|
---|
101 | _br_data;
|
---|
102 | DynArray <DataType>
|
---|
103 | _trigo_data;
|
---|
104 | Array <OscType, TRIGO_OSC_ARR_SIZE>
|
---|
105 | _trigo_osc;
|
---|
106 |
|
---|
107 |
|
---|
108 |
|
---|
109 | /*\\\ FORBIDDEN MEMBER FUNCTIONS \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\*/
|
---|
110 |
|
---|
111 | private:
|
---|
112 |
|
---|
113 | FFTRealFixLen (const FFTRealFixLen &other);
|
---|
114 | FFTRealFixLen& operator = (const FFTRealFixLen &other);
|
---|
115 | bool operator == (const FFTRealFixLen &other);
|
---|
116 | bool operator != (const FFTRealFixLen &other);
|
---|
117 |
|
---|
118 | }; // class FFTRealFixLen
|
---|
119 |
|
---|
120 |
|
---|
121 |
|
---|
122 | #include "FFTRealFixLen.hpp"
|
---|
123 |
|
---|
124 |
|
---|
125 |
|
---|
126 | #endif // FFTRealFixLen_HEADER_INCLUDED
|
---|
127 |
|
---|
128 |
|
---|
129 |
|
---|
130 | /*\\\ EOF \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\*/
|
---|