| 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 19 matching lines...) Expand all Loading... |
| 30 | 30 |
| 31 #if ENABLE(WEB_AUDIO) | 31 #if ENABLE(WEB_AUDIO) |
| 32 | 32 |
| 33 #include "core/platform/audio/Biquad.h" | 33 #include "core/platform/audio/Biquad.h" |
| 34 | 34 |
| 35 #include <stdio.h> | 35 #include <stdio.h> |
| 36 #include <algorithm> | 36 #include <algorithm> |
| 37 #include "core/platform/audio/DenormalDisabler.h" | 37 #include "core/platform/audio/DenormalDisabler.h" |
| 38 #include "wtf/MathExtras.h" | 38 #include "wtf/MathExtras.h" |
| 39 | 39 |
| 40 #if OS(DARWIN) | 40 #if OS(MACOSX) |
| 41 #include <Accelerate/Accelerate.h> | 41 #include <Accelerate/Accelerate.h> |
| 42 #endif | 42 #endif |
| 43 | 43 |
| 44 namespace WebCore { | 44 namespace WebCore { |
| 45 | 45 |
| 46 const int kBufferSize = 1024; | 46 const int kBufferSize = 1024; |
| 47 | 47 |
| 48 Biquad::Biquad() | 48 Biquad::Biquad() |
| 49 { | 49 { |
| 50 #if OS(DARWIN) | 50 #if OS(MACOSX) |
| 51 // Allocate two samples more for filter history | 51 // Allocate two samples more for filter history |
| 52 m_inputBuffer.allocate(kBufferSize + 2); | 52 m_inputBuffer.allocate(kBufferSize + 2); |
| 53 m_outputBuffer.allocate(kBufferSize + 2); | 53 m_outputBuffer.allocate(kBufferSize + 2); |
| 54 #endif | 54 #endif |
| 55 | 55 |
| 56 #if USE(WEBAUDIO_IPP) | 56 #if USE(WEBAUDIO_IPP) |
| 57 int bufferSize; | 57 int bufferSize; |
| 58 ippsIIRGetStateSize64f_BiQuad_32f(1, &bufferSize); | 58 ippsIIRGetStateSize64f_BiQuad_32f(1, &bufferSize); |
| 59 m_ippInternalBuffer = ippsMalloc_8u(bufferSize); | 59 m_ippInternalBuffer = ippsMalloc_8u(bufferSize); |
| 60 #endif // USE(WEBAUDIO_IPP) | 60 #endif // USE(WEBAUDIO_IPP) |
| 61 | 61 |
| 62 // Initialize as pass-thru (straight-wire, no filter effect) | 62 // Initialize as pass-thru (straight-wire, no filter effect) |
| 63 setNormalizedCoefficients(1, 0, 0, 1, 0, 0); | 63 setNormalizedCoefficients(1, 0, 0, 1, 0, 0); |
| 64 | 64 |
| 65 reset(); // clear filter memory | 65 reset(); // clear filter memory |
| 66 } | 66 } |
| 67 | 67 |
| 68 Biquad::~Biquad() | 68 Biquad::~Biquad() |
| 69 { | 69 { |
| 70 #if USE(WEBAUDIO_IPP) | 70 #if USE(WEBAUDIO_IPP) |
| 71 ippsFree(m_ippInternalBuffer); | 71 ippsFree(m_ippInternalBuffer); |
| 72 #endif // USE(WEBAUDIO_IPP) | 72 #endif // USE(WEBAUDIO_IPP) |
| 73 } | 73 } |
| 74 | 74 |
| 75 void Biquad::process(const float* sourceP, float* destP, size_t framesToProcess) | 75 void Biquad::process(const float* sourceP, float* destP, size_t framesToProcess) |
| 76 { | 76 { |
| 77 #if OS(DARWIN) | 77 #if OS(MACOSX) |
| 78 // Use vecLib if available | 78 // Use vecLib if available |
| 79 processFast(sourceP, destP, framesToProcess); | 79 processFast(sourceP, destP, framesToProcess); |
| 80 | 80 |
| 81 #elif USE(WEBAUDIO_IPP) | 81 #elif USE(WEBAUDIO_IPP) |
| 82 ippsIIR64f_32f(sourceP, destP, static_cast<int>(framesToProcess), m_biquadSt
ate); | 82 ippsIIR64f_32f(sourceP, destP, static_cast<int>(framesToProcess), m_biquadSt
ate); |
| 83 #else // USE(WEBAUDIO_IPP) | 83 #else // USE(WEBAUDIO_IPP) |
| 84 | 84 |
| 85 int n = framesToProcess; | 85 int n = framesToProcess; |
| 86 | 86 |
| 87 // Create local copies of member variables | 87 // Create local copies of member variables |
| (...skipping 30 matching lines...) Expand all Loading... |
| 118 m_y2 = DenormalDisabler::flushDenormalFloatToZero(y2); | 118 m_y2 = DenormalDisabler::flushDenormalFloatToZero(y2); |
| 119 | 119 |
| 120 m_b0 = b0; | 120 m_b0 = b0; |
| 121 m_b1 = b1; | 121 m_b1 = b1; |
| 122 m_b2 = b2; | 122 m_b2 = b2; |
| 123 m_a1 = a1; | 123 m_a1 = a1; |
| 124 m_a2 = a2; | 124 m_a2 = a2; |
| 125 #endif | 125 #endif |
| 126 } | 126 } |
| 127 | 127 |
| 128 #if OS(DARWIN) | 128 #if OS(MACOSX) |
| 129 | 129 |
| 130 // Here we have optimized version using Accelerate.framework | 130 // Here we have optimized version using Accelerate.framework |
| 131 | 131 |
| 132 void Biquad::processFast(const float* sourceP, float* destP, size_t framesToProc
ess) | 132 void Biquad::processFast(const float* sourceP, float* destP, size_t framesToProc
ess) |
| 133 { | 133 { |
| 134 double filterCoefficients[5]; | 134 double filterCoefficients[5]; |
| 135 filterCoefficients[0] = m_b0; | 135 filterCoefficients[0] = m_b0; |
| 136 filterCoefficients[1] = m_b1; | 136 filterCoefficients[1] = m_b1; |
| 137 filterCoefficients[2] = m_b2; | 137 filterCoefficients[2] = m_b2; |
| 138 filterCoefficients[3] = m_a1; | 138 filterCoefficients[3] = m_a1; |
| (...skipping 33 matching lines...) Expand 10 before | Expand all | Expand 10 after Loading... |
| 172 | 172 |
| 173 // Save history. Note that sourceP and destP reference m_inputBuffer and m_
outputBuffer respectively. | 173 // Save history. Note that sourceP and destP reference m_inputBuffer and m_
outputBuffer respectively. |
| 174 // These buffers are allocated (in the constructor) with space for two extra
samples so it's OK to access | 174 // These buffers are allocated (in the constructor) with space for two extra
samples so it's OK to access |
| 175 // array values two beyond framesToProcess. | 175 // array values two beyond framesToProcess. |
| 176 sourceP[0] = sourceP[framesToProcess - 2 + 2]; | 176 sourceP[0] = sourceP[framesToProcess - 2 + 2]; |
| 177 sourceP[1] = sourceP[framesToProcess - 1 + 2]; | 177 sourceP[1] = sourceP[framesToProcess - 1 + 2]; |
| 178 destP[0] = destP[framesToProcess - 2 + 2]; | 178 destP[0] = destP[framesToProcess - 2 + 2]; |
| 179 destP[1] = destP[framesToProcess - 1 + 2]; | 179 destP[1] = destP[framesToProcess - 1 + 2]; |
| 180 } | 180 } |
| 181 | 181 |
| 182 #endif // OS(DARWIN) | 182 #endif // OS(MACOSX) |
| 183 | 183 |
| 184 | 184 |
| 185 void Biquad::reset() | 185 void Biquad::reset() |
| 186 { | 186 { |
| 187 #if OS(DARWIN) | 187 #if OS(MACOSX) |
| 188 // Two extra samples for filter history | 188 // Two extra samples for filter history |
| 189 double* inputP = m_inputBuffer.data(); | 189 double* inputP = m_inputBuffer.data(); |
| 190 inputP[0] = 0; | 190 inputP[0] = 0; |
| 191 inputP[1] = 0; | 191 inputP[1] = 0; |
| 192 | 192 |
| 193 double* outputP = m_outputBuffer.data(); | 193 double* outputP = m_outputBuffer.data(); |
| 194 outputP[0] = 0; | 194 outputP[0] = 0; |
| 195 outputP[1] = 0; | 195 outputP[1] = 0; |
| 196 | 196 |
| 197 #elif USE(WEBAUDIO_IPP) | 197 #elif USE(WEBAUDIO_IPP) |
| (...skipping 380 matching lines...) Expand 10 before | Expand all | Expand 10 after Loading... |
| 578 Complex denominator = Complex(1, 0) + (a1 + a2 * z) * z; | 578 Complex denominator = Complex(1, 0) + (a1 + a2 * z) * z; |
| 579 Complex response = numerator / denominator; | 579 Complex response = numerator / denominator; |
| 580 magResponse[k] = static_cast<float>(abs(response)); | 580 magResponse[k] = static_cast<float>(abs(response)); |
| 581 phaseResponse[k] = static_cast<float>(atan2(imag(response), real(respons
e))); | 581 phaseResponse[k] = static_cast<float>(atan2(imag(response), real(respons
e))); |
| 582 } | 582 } |
| 583 } | 583 } |
| 584 | 584 |
| 585 } // namespace WebCore | 585 } // namespace WebCore |
| 586 | 586 |
| 587 #endif // ENABLE(WEB_AUDIO) | 587 #endif // ENABLE(WEB_AUDIO) |
| OLD | NEW |