| OLD | NEW |
| 1 /* | 1 /* |
| 2 * Copyright (C) 2010 Google Inc. All rights reserved. | 2 * Copyright (C) 2010 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 12 matching lines...) Expand all Loading... |
| 23 * LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND | 23 * LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND |
| 24 * ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT | 24 * ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT |
| 25 * (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF | 25 * (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF |
| 26 * THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. | 26 * THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. |
| 27 */ | 27 */ |
| 28 | 28 |
| 29 #include "platform/audio/FFTFrame.h" | 29 #include "platform/audio/FFTFrame.h" |
| 30 #include "platform/audio/VectorMath.h" | 30 #include "platform/audio/VectorMath.h" |
| 31 #include "platform/Logging.h" | 31 #include "platform/Logging.h" |
| 32 #include "wtf/MathExtras.h" | 32 #include "wtf/MathExtras.h" |
| 33 #include "wtf/OwnPtr.h" | 33 #include "wtf/PtrUtil.h" |
| 34 #include <complex> | 34 #include <complex> |
| 35 #include <memory> |
| 35 | 36 |
| 36 #ifndef NDEBUG | 37 #ifndef NDEBUG |
| 37 #include <stdio.h> | 38 #include <stdio.h> |
| 38 #endif | 39 #endif |
| 39 | 40 |
| 40 namespace blink { | 41 namespace blink { |
| 41 | 42 |
| 42 void FFTFrame::doPaddedFFT(const float* data, size_t dataSize) | 43 void FFTFrame::doPaddedFFT(const float* data, size_t dataSize) |
| 43 { | 44 { |
| 44 // Zero-pad the impulse response | 45 // Zero-pad the impulse response |
| 45 AudioFloatArray paddedResponse(fftSize()); // zero-initialized | 46 AudioFloatArray paddedResponse(fftSize()); // zero-initialized |
| 46 paddedResponse.copyToRange(data, 0, dataSize); | 47 paddedResponse.copyToRange(data, 0, dataSize); |
| 47 | 48 |
| 48 // Get the frequency-domain version of padded response | 49 // Get the frequency-domain version of padded response |
| 49 doFFT(paddedResponse.data()); | 50 doFFT(paddedResponse.data()); |
| 50 } | 51 } |
| 51 | 52 |
| 52 PassOwnPtr<FFTFrame> FFTFrame::createInterpolatedFrame(const FFTFrame& frame1, c
onst FFTFrame& frame2, double x) | 53 std::unique_ptr<FFTFrame> FFTFrame::createInterpolatedFrame(const FFTFrame& fram
e1, const FFTFrame& frame2, double x) |
| 53 { | 54 { |
| 54 OwnPtr<FFTFrame> newFrame = adoptPtr(new FFTFrame(frame1.fftSize())); | 55 std::unique_ptr<FFTFrame> newFrame = wrapUnique(new FFTFrame(frame1.fftSize(
))); |
| 55 | 56 |
| 56 newFrame->interpolateFrequencyComponents(frame1, frame2, x); | 57 newFrame->interpolateFrequencyComponents(frame1, frame2, x); |
| 57 | 58 |
| 58 // In the time-domain, the 2nd half of the response must be zero, to avoid c
ircular convolution aliasing... | 59 // In the time-domain, the 2nd half of the response must be zero, to avoid c
ircular convolution aliasing... |
| 59 int fftSize = newFrame->fftSize(); | 60 int fftSize = newFrame->fftSize(); |
| 60 AudioFloatArray buffer(fftSize); | 61 AudioFloatArray buffer(fftSize); |
| 61 newFrame->doInverseFFT(buffer.data()); | 62 newFrame->doInverseFFT(buffer.data()); |
| 62 buffer.zeroRange(fftSize / 2, fftSize); | 63 buffer.zeroRange(fftSize / 2, fftSize); |
| 63 | 64 |
| 64 // Put back into frequency domain. | 65 // Put back into frequency domain. |
| (...skipping 195 matching lines...) Expand 10 before | Expand all | Expand 10 after Loading... |
| 260 | 261 |
| 261 VectorMath::zvmul(realP1, imagP1, realP2, imagP2, realP1, imagP1, halfSize); | 262 VectorMath::zvmul(realP1, imagP1, realP2, imagP2, realP1, imagP1, halfSize); |
| 262 | 263 |
| 263 // Multiply the packed DC/nyquist component | 264 // Multiply the packed DC/nyquist component |
| 264 realP1[0] = real0 * realP2[0]; | 265 realP1[0] = real0 * realP2[0]; |
| 265 imagP1[0] = imag0 * imagP2[0]; | 266 imagP1[0] = imag0 * imagP2[0]; |
| 266 } | 267 } |
| 267 | 268 |
| 268 } // namespace blink | 269 } // namespace blink |
| 269 | 270 |
| OLD | NEW |