|
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 \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\*/ |