| OLD | NEW |
| 1 /* | 1 /* |
| 2 * Copyright (C) 2012 Google Inc. All rights reserved. | 2 * Copyright (C) 2012 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 "config.h" | 29 #include "config.h" |
| 30 | 30 |
| 31 #if ENABLE(WEB_AUDIO) | 31 #if ENABLE(WEB_AUDIO) |
| 32 | 32 |
| 33 #include "modules/webaudio/WaveTable.h" | 33 #include "modules/webaudio/PeriodicWave.h" |
| 34 | 34 |
| 35 #include "core/platform/audio/FFTFrame.h" | 35 #include "core/platform/audio/FFTFrame.h" |
| 36 #include "core/platform/audio/VectorMath.h" | 36 #include "core/platform/audio/VectorMath.h" |
| 37 #include "modules/webaudio/OscillatorNode.h" | 37 #include "modules/webaudio/OscillatorNode.h" |
| 38 #include "wtf/OwnPtr.h" |
| 38 #include <algorithm> | 39 #include <algorithm> |
| 39 #include "wtf/OwnPtr.h" | |
| 40 | 40 |
| 41 const unsigned WaveTableSize = 4096; // This must be a power of two. | 41 const unsigned PeriodicWaveSize = 4096; // This must be a power of two. |
| 42 const unsigned NumberOfRanges = 36; // There should be 3 * log2(WaveTableSize) 1
/3 octave ranges. | 42 const unsigned NumberOfRanges = 36; // There should be 3 * log2(PeriodicWaveSize
) 1/3 octave ranges. |
| 43 const float CentsPerRange = 1200 / 3; // 1/3 Octave. | 43 const float CentsPerRange = 1200 / 3; // 1/3 Octave. |
| 44 | 44 |
| 45 namespace WebCore { | 45 namespace WebCore { |
| 46 | 46 |
| 47 using namespace VectorMath; | 47 using namespace VectorMath; |
| 48 | 48 |
| 49 PassRefPtr<WaveTable> WaveTable::create(float sampleRate, Float32Array* real, Fl
oat32Array* imag) | 49 PassRefPtr<PeriodicWave> PeriodicWave::create(float sampleRate, Float32Array* re
al, Float32Array* imag) |
| 50 { | 50 { |
| 51 bool isGood = real && imag && real->length() == imag->length(); | 51 bool isGood = real && imag && real->length() == imag->length(); |
| 52 ASSERT(isGood); | 52 ASSERT(isGood); |
| 53 if (isGood) { | 53 if (isGood) { |
| 54 RefPtr<WaveTable> waveTable = adoptRef(new WaveTable(sampleRate)); | 54 RefPtr<PeriodicWave> periodicWave = adoptRef(new PeriodicWave(sampleRate
)); |
| 55 size_t numberOfComponents = real->length(); | 55 size_t numberOfComponents = real->length(); |
| 56 waveTable->createBandLimitedTables(real->data(), imag->data(), numberOfC
omponents); | 56 periodicWave->createBandLimitedTables(real->data(), imag->data(), number
OfComponents); |
| 57 return waveTable; | 57 return periodicWave; |
| 58 } | 58 } |
| 59 return 0; | 59 return 0; |
| 60 } | 60 } |
| 61 | 61 |
| 62 PassRefPtr<WaveTable> WaveTable::createSine(float sampleRate) | 62 PassRefPtr<PeriodicWave> PeriodicWave::createSine(float sampleRate) |
| 63 { | 63 { |
| 64 RefPtr<WaveTable> waveTable = adoptRef(new WaveTable(sampleRate)); | 64 RefPtr<PeriodicWave> periodicWave = adoptRef(new PeriodicWave(sampleRate)); |
| 65 waveTable->generateBasicWaveform(OscillatorNode::SINE); | 65 periodicWave->generateBasicWaveform(OscillatorNode::SINE); |
| 66 return waveTable; | 66 return periodicWave; |
| 67 } | 67 } |
| 68 | 68 |
| 69 PassRefPtr<WaveTable> WaveTable::createSquare(float sampleRate) | 69 PassRefPtr<PeriodicWave> PeriodicWave::createSquare(float sampleRate) |
| 70 { | 70 { |
| 71 RefPtr<WaveTable> waveTable = adoptRef(new WaveTable(sampleRate)); | 71 RefPtr<PeriodicWave> periodicWave = adoptRef(new PeriodicWave(sampleRate)); |
| 72 waveTable->generateBasicWaveform(OscillatorNode::SQUARE); | 72 periodicWave->generateBasicWaveform(OscillatorNode::SQUARE); |
| 73 return waveTable; | 73 return periodicWave; |
| 74 } | 74 } |
| 75 | 75 |
| 76 PassRefPtr<WaveTable> WaveTable::createSawtooth(float sampleRate) | 76 PassRefPtr<PeriodicWave> PeriodicWave::createSawtooth(float sampleRate) |
| 77 { | 77 { |
| 78 RefPtr<WaveTable> waveTable = adoptRef(new WaveTable(sampleRate)); | 78 RefPtr<PeriodicWave> periodicWave = adoptRef(new PeriodicWave(sampleRate)); |
| 79 waveTable->generateBasicWaveform(OscillatorNode::SAWTOOTH); | 79 periodicWave->generateBasicWaveform(OscillatorNode::SAWTOOTH); |
| 80 return waveTable; | 80 return periodicWave; |
| 81 } | 81 } |
| 82 | 82 |
| 83 PassRefPtr<WaveTable> WaveTable::createTriangle(float sampleRate) | 83 PassRefPtr<PeriodicWave> PeriodicWave::createTriangle(float sampleRate) |
| 84 { | 84 { |
| 85 RefPtr<WaveTable> waveTable = adoptRef(new WaveTable(sampleRate)); | 85 RefPtr<PeriodicWave> periodicWave = adoptRef(new PeriodicWave(sampleRate)); |
| 86 waveTable->generateBasicWaveform(OscillatorNode::TRIANGLE); | 86 periodicWave->generateBasicWaveform(OscillatorNode::TRIANGLE); |
| 87 return waveTable; | 87 return periodicWave; |
| 88 } | 88 } |
| 89 | 89 |
| 90 WaveTable::WaveTable(float sampleRate) | 90 PeriodicWave::PeriodicWave(float sampleRate) |
| 91 : m_sampleRate(sampleRate) | 91 : m_sampleRate(sampleRate) |
| 92 , m_waveTableSize(WaveTableSize) | 92 , m_periodicWaveSize(PeriodicWaveSize) |
| 93 , m_numberOfRanges(NumberOfRanges) | 93 , m_numberOfRanges(NumberOfRanges) |
| 94 , m_centsPerRange(CentsPerRange) | 94 , m_centsPerRange(CentsPerRange) |
| 95 { | 95 { |
| 96 ScriptWrappable::init(this); | 96 ScriptWrappable::init(this); |
| 97 float nyquist = 0.5 * m_sampleRate; | 97 float nyquist = 0.5 * m_sampleRate; |
| 98 m_lowestFundamentalFrequency = nyquist / maxNumberOfPartials(); | 98 m_lowestFundamentalFrequency = nyquist / maxNumberOfPartials(); |
| 99 m_rateScale = m_waveTableSize / m_sampleRate; | 99 m_rateScale = m_periodicWaveSize / m_sampleRate; |
| 100 } | 100 } |
| 101 | 101 |
| 102 void WaveTable::waveDataForFundamentalFrequency(float fundamentalFrequency, floa
t* &lowerWaveData, float* &higherWaveData, float& tableInterpolationFactor) | 102 void PeriodicWave::waveDataForFundamentalFrequency(float fundamentalFrequency, f
loat* &lowerWaveData, float* &higherWaveData, float& tableInterpolationFactor) |
| 103 { | 103 { |
| 104 // Negative frequencies are allowed, in which case we alias to the positive
frequency. | 104 // Negative frequencies are allowed, in which case we alias to the positive
frequency. |
| 105 fundamentalFrequency = fabsf(fundamentalFrequency); | 105 fundamentalFrequency = fabsf(fundamentalFrequency); |
| 106 | 106 |
| 107 // Calculate the pitch range. | 107 // Calculate the pitch range. |
| 108 float ratio = fundamentalFrequency > 0 ? fundamentalFrequency / m_lowestFund
amentalFrequency : 0.5; | 108 float ratio = fundamentalFrequency > 0 ? fundamentalFrequency / m_lowestFund
amentalFrequency : 0.5; |
| 109 float centsAboveLowestFrequency = log2f(ratio) * 1200; | 109 float centsAboveLowestFrequency = log2f(ratio) * 1200; |
| 110 | 110 |
| 111 // Add one to round-up to the next range just in time to truncate partials b
efore aliasing occurs. | 111 // Add one to round-up to the next range just in time to truncate partials b
efore aliasing occurs. |
| 112 float pitchRange = 1 + centsAboveLowestFrequency / m_centsPerRange; | 112 float pitchRange = 1 + centsAboveLowestFrequency / m_centsPerRange; |
| 113 | 113 |
| 114 pitchRange = std::max(pitchRange, 0.0f); | 114 pitchRange = std::max(pitchRange, 0.0f); |
| 115 pitchRange = std::min(pitchRange, static_cast<float>(m_numberOfRanges - 1)); | 115 pitchRange = std::min(pitchRange, static_cast<float>(m_numberOfRanges - 1)); |
| 116 | 116 |
| 117 // The words "lower" and "higher" refer to the table data having the lower a
nd higher numbers of partials. | 117 // The words "lower" and "higher" refer to the table data having the lower a
nd higher numbers of partials. |
| 118 // It's a little confusing since the range index gets larger the more partia
ls we cull out. | 118 // It's a little confusing since the range index gets larger the more partia
ls we cull out. |
| 119 // So the lower table data will have a larger range index. | 119 // So the lower table data will have a larger range index. |
| 120 unsigned rangeIndex1 = static_cast<unsigned>(pitchRange); | 120 unsigned rangeIndex1 = static_cast<unsigned>(pitchRange); |
| 121 unsigned rangeIndex2 = rangeIndex1 < m_numberOfRanges - 1 ? rangeIndex1 + 1
: rangeIndex1; | 121 unsigned rangeIndex2 = rangeIndex1 < m_numberOfRanges - 1 ? rangeIndex1 + 1
: rangeIndex1; |
| 122 | 122 |
| 123 lowerWaveData = m_bandLimitedTables[rangeIndex2]->data(); | 123 lowerWaveData = m_bandLimitedTables[rangeIndex2]->data(); |
| 124 higherWaveData = m_bandLimitedTables[rangeIndex1]->data(); | 124 higherWaveData = m_bandLimitedTables[rangeIndex1]->data(); |
| 125 | 125 |
| 126 // Ranges from 0 -> 1 to interpolate between lower -> higher. | 126 // Ranges from 0 -> 1 to interpolate between lower -> higher. |
| 127 tableInterpolationFactor = pitchRange - rangeIndex1; | 127 tableInterpolationFactor = pitchRange - rangeIndex1; |
| 128 } | 128 } |
| 129 | 129 |
| 130 unsigned WaveTable::maxNumberOfPartials() const | 130 unsigned PeriodicWave::maxNumberOfPartials() const |
| 131 { | 131 { |
| 132 return m_waveTableSize / 2; | 132 return m_periodicWaveSize / 2; |
| 133 } | 133 } |
| 134 | 134 |
| 135 unsigned WaveTable::numberOfPartialsForRange(unsigned rangeIndex) const | 135 unsigned PeriodicWave::numberOfPartialsForRange(unsigned rangeIndex) const |
| 136 { | 136 { |
| 137 // Number of cents below nyquist where we cull partials. | 137 // Number of cents below nyquist where we cull partials. |
| 138 float centsToCull = rangeIndex * m_centsPerRange; | 138 float centsToCull = rangeIndex * m_centsPerRange; |
| 139 | 139 |
| 140 // A value from 0 -> 1 representing what fraction of the partials to keep. | 140 // A value from 0 -> 1 representing what fraction of the partials to keep. |
| 141 float cullingScale = pow(2, -centsToCull / 1200); | 141 float cullingScale = pow(2, -centsToCull / 1200); |
| 142 | 142 |
| 143 // The very top range will have all the partials culled. | 143 // The very top range will have all the partials culled. |
| 144 unsigned numberOfPartials = cullingScale * maxNumberOfPartials(); | 144 unsigned numberOfPartials = cullingScale * maxNumberOfPartials(); |
| 145 | 145 |
| 146 return numberOfPartials; | 146 return numberOfPartials; |
| 147 } | 147 } |
| 148 | 148 |
| 149 // Convert into time-domain wave tables. | 149 // Convert into time-domain wave buffers. |
| 150 // One table is created for each range for non-aliasing playback at different pl
ayback rates. | 150 // One table is created for each range for non-aliasing playback at different pl
ayback rates. |
| 151 // Thus, higher ranges have more high-frequency partials culled out. | 151 // Thus, higher ranges have more high-frequency partials culled out. |
| 152 void WaveTable::createBandLimitedTables(const float* realData, const float* imag
Data, unsigned numberOfComponents) | 152 void PeriodicWave::createBandLimitedTables(const float* realData, const float* i
magData, unsigned numberOfComponents) |
| 153 { | 153 { |
| 154 float normalizationScale = 1; | 154 float normalizationScale = 1; |
| 155 | 155 |
| 156 unsigned fftSize = m_waveTableSize; | 156 unsigned fftSize = m_periodicWaveSize; |
| 157 unsigned halfSize = fftSize / 2; | 157 unsigned halfSize = fftSize / 2; |
| 158 unsigned i; | 158 unsigned i; |
| 159 | 159 |
| 160 numberOfComponents = std::min(numberOfComponents, halfSize); | 160 numberOfComponents = std::min(numberOfComponents, halfSize); |
| 161 | 161 |
| 162 m_bandLimitedTables.reserveCapacity(m_numberOfRanges); | 162 m_bandLimitedTables.reserveCapacity(m_numberOfRanges); |
| 163 | 163 |
| 164 for (unsigned rangeIndex = 0; rangeIndex < m_numberOfRanges; ++rangeIndex) { | 164 for (unsigned rangeIndex = 0; rangeIndex < m_numberOfRanges; ++rangeIndex) { |
| 165 // This FFTFrame is used to cull partials (represented by frequency bins
). | 165 // This FFTFrame is used to cull partials (represented by frequency bins
). |
| 166 FFTFrame frame(fftSize); | 166 FFTFrame frame(fftSize); |
| 167 float* realP = frame.realData(); | 167 float* realP = frame.realData(); |
| 168 float* imagP = frame.imagData(); | 168 float* imagP = frame.imagData(); |
| 169 | 169 |
| 170 // Copy from loaded frequency data and scale. | 170 // Copy from loaded frequency data and scale. |
| 171 float scale = fftSize; | 171 float scale = fftSize; |
| 172 vsmul(realData, 1, &scale, realP, 1, numberOfComponents); | 172 vsmul(realData, 1, &scale, realP, 1, numberOfComponents); |
| 173 vsmul(imagData, 1, &scale, imagP, 1, numberOfComponents); | 173 vsmul(imagData, 1, &scale, imagP, 1, numberOfComponents); |
| 174 | 174 |
| 175 // If fewer components were provided than 1/2 FFT size, then clear the r
emaining bins. | 175 // If fewer components were provided than 1/2 FFT size, then clear the r
emaining bins. |
| 176 for (i = numberOfComponents; i < halfSize; ++i) { | 176 for (i = numberOfComponents; i < halfSize; ++i) { |
| 177 realP[i] = 0; | 177 realP[i] = 0; |
| 178 imagP[i] = 0; | 178 imagP[i] = 0; |
| 179 } | 179 } |
| 180 | 180 |
| 181 // Generate complex conjugate because of the way the inverse FFT is defi
ned. | 181 // Generate complex conjugate because of the way the inverse FFT is defi
ned. |
| 182 float minusOne = -1; | 182 float minusOne = -1; |
| 183 vsmul(imagP, 1, &minusOne, imagP, 1, halfSize); | 183 vsmul(imagP, 1, &minusOne, imagP, 1, halfSize); |
| 184 | 184 |
| 185 // Find the starting bin where we should start culling. | 185 // Find the starting bin where we should start culling. |
| 186 // We need to clear out the highest frequencies to band-limit the wavefo
rm. | 186 // We need to clear out the highest frequencies to band-limit the wavefo
rm. |
| 187 unsigned numberOfPartials = numberOfPartialsForRange(rangeIndex); | 187 unsigned numberOfPartials = numberOfPartialsForRange(rangeIndex); |
| 188 | 188 |
| 189 // Cull the aliasing partials for this pitch range. | 189 // Cull the aliasing partials for this pitch range. |
| 190 for (i = numberOfPartials + 1; i < halfSize; ++i) { | 190 for (i = numberOfPartials + 1; i < halfSize; ++i) { |
| 191 realP[i] = 0; | 191 realP[i] = 0; |
| 192 imagP[i] = 0; | 192 imagP[i] = 0; |
| 193 } | 193 } |
| 194 // Clear packed-nyquist if necessary. | 194 // Clear packed-nyquist if necessary. |
| 195 if (numberOfPartials < halfSize) | 195 if (numberOfPartials < halfSize) |
| 196 imagP[0] = 0; | 196 imagP[0] = 0; |
| 197 | 197 |
| 198 // Clear any DC-offset. | 198 // Clear any DC-offset. |
| 199 realP[0] = 0; | 199 realP[0] = 0; |
| 200 | 200 |
| 201 // Create the band-limited table. | 201 // Create the band-limited table. |
| 202 OwnPtr<AudioFloatArray> table = adoptPtr(new AudioFloatArray(m_waveTable
Size)); | 202 OwnPtr<AudioFloatArray> table = adoptPtr(new AudioFloatArray(m_periodicW
aveSize)); |
| 203 m_bandLimitedTables.append(table.release()); | 203 m_bandLimitedTables.append(table.release()); |
| 204 | 204 |
| 205 // Apply an inverse FFT to generate the time-domain table data. | 205 // Apply an inverse FFT to generate the time-domain table data. |
| 206 float* data = m_bandLimitedTables[rangeIndex]->data(); | 206 float* data = m_bandLimitedTables[rangeIndex]->data(); |
| 207 frame.doInverseFFT(data); | 207 frame.doInverseFFT(data); |
| 208 | 208 |
| 209 // For the first range (which has the highest power), calculate its peak
value then compute normalization scale. | 209 // For the first range (which has the highest power), calculate its peak
value then compute normalization scale. |
| 210 if (!rangeIndex) { | 210 if (!rangeIndex) { |
| 211 float maxValue; | 211 float maxValue; |
| 212 vmaxmgv(data, 1, &maxValue, m_waveTableSize); | 212 vmaxmgv(data, 1, &maxValue, m_periodicWaveSize); |
| 213 | 213 |
| 214 if (maxValue) | 214 if (maxValue) |
| 215 normalizationScale = 1.0f / maxValue; | 215 normalizationScale = 1.0f / maxValue; |
| 216 } | 216 } |
| 217 | 217 |
| 218 // Apply normalization scale. | 218 // Apply normalization scale. |
| 219 vsmul(data, 1, &normalizationScale, data, 1, m_waveTableSize); | 219 vsmul(data, 1, &normalizationScale, data, 1, m_periodicWaveSize); |
| 220 } | 220 } |
| 221 } | 221 } |
| 222 | 222 |
| 223 void WaveTable::generateBasicWaveform(int shape) | 223 void PeriodicWave::generateBasicWaveform(int shape) |
| 224 { | 224 { |
| 225 unsigned fftSize = waveTableSize(); | 225 unsigned fftSize = periodicWaveSize(); |
| 226 unsigned halfSize = fftSize / 2; | 226 unsigned halfSize = fftSize / 2; |
| 227 | 227 |
| 228 AudioFloatArray real(halfSize); | 228 AudioFloatArray real(halfSize); |
| 229 AudioFloatArray imag(halfSize); | 229 AudioFloatArray imag(halfSize); |
| 230 float* realP = real.data(); | 230 float* realP = real.data(); |
| 231 float* imagP = imag.data(); | 231 float* imagP = imag.data(); |
| 232 | 232 |
| 233 // Clear DC and Nyquist. | 233 // Clear DC and Nyquist. |
| 234 realP[0] = 0; | 234 realP[0] = 0; |
| 235 imagP[0] = 0; | 235 imagP[0] = 0; |
| (...skipping 39 matching lines...) Expand 10 before | Expand all | Expand 10 after Loading... |
| 275 realP[n] = a; | 275 realP[n] = a; |
| 276 imagP[n] = b; | 276 imagP[n] = b; |
| 277 } | 277 } |
| 278 | 278 |
| 279 createBandLimitedTables(realP, imagP, halfSize); | 279 createBandLimitedTables(realP, imagP, halfSize); |
| 280 } | 280 } |
| 281 | 281 |
| 282 } // namespace WebCore | 282 } // namespace WebCore |
| 283 | 283 |
| 284 #endif // ENABLE(WEB_AUDIO) | 284 #endif // ENABLE(WEB_AUDIO) |
| OLD | NEW |