OLD | NEW |
1 /* | 1 /* |
2 * Copyright (C) 2011 Google Inc. All rights reserved. | 2 * Copyright (C) 2011 Google Inc. All rights reserved. |
3 * | 3 * |
4 * Redistribution and use in source and binary forms, with or without | 4 * Redistribution and use in source and binary forms, with or without |
5 * modification, are permitted provided that the following conditions | 5 * modification, are permitted provided that the following conditions |
6 * are met: | 6 * are met: |
7 * | 7 * |
8 * 1. Redistributions of source code must retain the above copyright | 8 * 1. Redistributions of source code must retain the above copyright |
9 * notice, this list of conditions and the following disclaimer. | 9 * notice, this list of conditions and the following disclaimer. |
10 * 2. Redistributions in binary form must reproduce the above copyright | 10 * 2. Redistributions in binary form must reproduce the above copyright |
(...skipping 34 matching lines...) Expand 10 before | Expand all | Expand 10 after Loading... |
45 namespace blink { | 45 namespace blink { |
46 | 46 |
47 #if ENABLE(ASSERT) | 47 #if ENABLE(ASSERT) |
48 const int kMaxFFTPow2Size = 24; | 48 const int kMaxFFTPow2Size = 24; |
49 #endif | 49 #endif |
50 | 50 |
51 // Normal constructor: allocates for a given fftSize. | 51 // Normal constructor: allocates for a given fftSize. |
52 FFTFrame::FFTFrame(unsigned fftSize) | 52 FFTFrame::FFTFrame(unsigned fftSize) |
53 : m_FFTSize(fftSize) | 53 : m_FFTSize(fftSize) |
54 , m_log2FFTSize(static_cast<unsigned>(log2(fftSize))) | 54 , m_log2FFTSize(static_cast<unsigned>(log2(fftSize))) |
| 55 , m_realData(fftSize / 2) |
| 56 , m_imagData(fftSize / 2) |
55 , m_forwardContext(0) | 57 , m_forwardContext(0) |
56 , m_inverseContext(0) | 58 , m_inverseContext(0) |
57 , m_complexData(fftSize) | 59 , m_complexData(fftSize) |
58 , m_realData(fftSize / 2) | |
59 , m_imagData(fftSize / 2) | |
60 { | 60 { |
61 // We only allow power of two. | 61 // We only allow power of two. |
62 ASSERT(1UL << m_log2FFTSize == m_FFTSize); | 62 ASSERT(1UL << m_log2FFTSize == m_FFTSize); |
63 | 63 |
64 m_forwardContext = contextForSize(fftSize, DFT_R2C); | 64 m_forwardContext = contextForSize(fftSize, DFT_R2C); |
65 m_inverseContext = contextForSize(fftSize, IDFT_C2R); | 65 m_inverseContext = contextForSize(fftSize, IDFT_C2R); |
66 } | 66 } |
67 | 67 |
68 // Creates a blank/empty frame (interpolate() must later be called). | 68 // Creates a blank/empty frame (interpolate() must later be called). |
69 FFTFrame::FFTFrame() | 69 FFTFrame::FFTFrame() |
70 : m_FFTSize(0) | 70 : m_FFTSize(0) |
71 , m_log2FFTSize(0) | 71 , m_log2FFTSize(0) |
72 , m_forwardContext(0) | 72 , m_forwardContext(0) |
73 , m_inverseContext(0) | 73 , m_inverseContext(0) |
74 { | 74 { |
75 } | 75 } |
76 | 76 |
77 // Copy constructor. | 77 // Copy constructor. |
78 FFTFrame::FFTFrame(const FFTFrame& frame) | 78 FFTFrame::FFTFrame(const FFTFrame& frame) |
79 : m_FFTSize(frame.m_FFTSize) | 79 : m_FFTSize(frame.m_FFTSize) |
80 , m_log2FFTSize(frame.m_log2FFTSize) | 80 , m_log2FFTSize(frame.m_log2FFTSize) |
| 81 , m_realData(frame.m_FFTSize / 2) |
| 82 , m_imagData(frame.m_FFTSize / 2) |
81 , m_forwardContext(0) | 83 , m_forwardContext(0) |
82 , m_inverseContext(0) | 84 , m_inverseContext(0) |
83 , m_complexData(frame.m_FFTSize) | 85 , m_complexData(frame.m_FFTSize) |
84 , m_realData(frame.m_FFTSize / 2) | |
85 , m_imagData(frame.m_FFTSize / 2) | |
86 { | 86 { |
87 m_forwardContext = contextForSize(m_FFTSize, DFT_R2C); | 87 m_forwardContext = contextForSize(m_FFTSize, DFT_R2C); |
88 m_inverseContext = contextForSize(m_FFTSize, IDFT_C2R); | 88 m_inverseContext = contextForSize(m_FFTSize, IDFT_C2R); |
89 | 89 |
90 // Copy/setup frame data. | 90 // Copy/setup frame data. |
91 unsigned nbytes = sizeof(float) * (m_FFTSize / 2); | 91 unsigned nbytes = sizeof(float) * (m_FFTSize / 2); |
92 memcpy(realData(), frame.realData(), nbytes); | 92 memcpy(realData(), frame.realData(), nbytes); |
93 memcpy(imagData(), frame.imagData(), nbytes); | 93 memcpy(imagData(), frame.imagData(), nbytes); |
94 } | 94 } |
95 | 95 |
(...skipping 42 matching lines...) Expand 10 before | Expand all | Expand 10 after Loading... |
138 // Compute inverse transform. | 138 // Compute inverse transform. |
139 av_rdft_calc(m_inverseContext, interleavedData); | 139 av_rdft_calc(m_inverseContext, interleavedData); |
140 | 140 |
141 // Scale so that a forward then inverse FFT yields exactly the original data
. For some reason | 141 // Scale so that a forward then inverse FFT yields exactly the original data
. For some reason |
142 // av_rdft_calc above returns values that are half of what I expect. Hence m
ake the scale factor | 142 // av_rdft_calc above returns values that are half of what I expect. Hence m
ake the scale factor |
143 // twice as large to compensate for that. | 143 // twice as large to compensate for that. |
144 const float scale = 2.0 / m_FFTSize; | 144 const float scale = 2.0 / m_FFTSize; |
145 VectorMath::vsmul(interleavedData, 1, &scale, data, 1, m_FFTSize); | 145 VectorMath::vsmul(interleavedData, 1, &scale, data, 1, m_FFTSize); |
146 } | 146 } |
147 | 147 |
148 float* FFTFrame::realData() const | |
149 { | |
150 return const_cast<float*>(m_realData.data()); | |
151 } | |
152 | |
153 float* FFTFrame::imagData() const | |
154 { | |
155 return const_cast<float*>(m_imagData.data()); | |
156 } | |
157 | |
158 float* FFTFrame::getUpToDateComplexData() | 148 float* FFTFrame::getUpToDateComplexData() |
159 { | 149 { |
160 // FIXME: if we can't completely get rid of this method, SSE | 150 // FIXME: if we can't completely get rid of this method, SSE |
161 // optimization could be considered if it shows up hot on profiles. | 151 // optimization could be considered if it shows up hot on profiles. |
162 int len = m_FFTSize / 2; | 152 int len = m_FFTSize / 2; |
163 const float* real = m_realData.data(); | 153 const float* real = m_realData.data(); |
164 const float* imag = m_imagData.data(); | 154 const float* imag = m_imagData.data(); |
165 float* c = m_complexData.data(); | 155 float* c = m_complexData.data(); |
166 for (int i = 0; i < len; ++i) { | 156 for (int i = 0; i < len; ++i) { |
167 int baseComplexIndex = 2 * i; | 157 int baseComplexIndex = 2 * i; |
(...skipping 14 matching lines...) Expand all Loading... |
182 | 172 |
183 RDFTContext* context = av_rdft_init(pow2size, (RDFTransformType)trans); | 173 RDFTContext* context = av_rdft_init(pow2size, (RDFTransformType)trans); |
184 return context; | 174 return context; |
185 } | 175 } |
186 | 176 |
187 } // namespace blink | 177 } // namespace blink |
188 | 178 |
189 #endif // USE(WEBAUDIO_FFMPEG) | 179 #endif // USE(WEBAUDIO_FFMPEG) |
190 | 180 |
191 #endif // ENABLE(WEB_AUDIO) | 181 #endif // ENABLE(WEB_AUDIO) |
OLD | NEW |