#include "config.h"
#if ENABLE(WEB_AUDIO)
#include "AudioBus.h"
#include "DenormalDisabler.h"
#include "SincResampler.h"
#include "VectorMath.h"
#include <algorithm>
#include <assert.h>
#include <math.h>
namespace WebCore {
constexpr unsigned MaxBusChannels = 32;
RefPtr<AudioBus> AudioBus::create(unsigned numberOfChannels, size_t length, bool allocate)
{
ASSERT(numberOfChannels <= MaxBusChannels);
if (numberOfChannels > MaxBusChannels)
return nullptr;
return adoptRef(*new AudioBus(numberOfChannels, length, allocate));
}
AudioBus::AudioBus(unsigned numberOfChannels, size_t length, bool allocate)
: m_length(length)
{
m_channels.reserveInitialCapacity(numberOfChannels);
for (unsigned i = 0; i < numberOfChannels; ++i) {
auto channel = allocate ? makeUnique<AudioChannel>(length) : makeUnique<AudioChannel>(nullptr, length);
m_channels.uncheckedAppend(WTFMove(channel));
}
m_layout = LayoutCanonical; }
void AudioBus::setChannelMemory(unsigned channelIndex, float* storage, size_t length)
{
if (channelIndex < m_channels.size()) {
channel(channelIndex)->set(storage, length);
m_length = length; }
}
void AudioBus::resizeSmaller(size_t newLength)
{
ASSERT(newLength <= m_length);
if (newLength <= m_length)
m_length = newLength;
for (unsigned i = 0; i < m_channels.size(); ++i)
m_channels[i]->resizeSmaller(newLength);
}
void AudioBus::zero()
{
for (unsigned i = 0; i < m_channels.size(); ++i)
m_channels[i]->zero();
}
AudioChannel* AudioBus::channelByType(unsigned channelType)
{
if (m_layout != LayoutCanonical)
return 0;
switch (numberOfChannels()) {
case 1: if (channelType == ChannelMono || channelType == ChannelLeft)
return channel(0);
return 0;
case 2: switch (channelType) {
case ChannelLeft: return channel(0);
case ChannelRight: return channel(1);
default: return 0;
}
case 4: switch (channelType) {
case ChannelLeft: return channel(0);
case ChannelRight: return channel(1);
case ChannelSurroundLeft: return channel(2);
case ChannelSurroundRight: return channel(3);
default: return 0;
}
case 5: switch (channelType) {
case ChannelLeft: return channel(0);
case ChannelRight: return channel(1);
case ChannelCenter: return channel(2);
case ChannelSurroundLeft: return channel(3);
case ChannelSurroundRight: return channel(4);
default: return 0;
}
case 6: switch (channelType) {
case ChannelLeft: return channel(0);
case ChannelRight: return channel(1);
case ChannelCenter: return channel(2);
case ChannelLFE: return channel(3);
case ChannelSurroundLeft: return channel(4);
case ChannelSurroundRight: return channel(5);
default: return 0;
}
}
ASSERT_NOT_REACHED();
return 0;
}
const AudioChannel* AudioBus::channelByType(unsigned type) const
{
return const_cast<AudioBus*>(this)->channelByType(type);
}
bool AudioBus::topologyMatches(const AudioBus& bus) const
{
if (numberOfChannels() != bus.numberOfChannels())
return false;
if (length() > bus.length())
return false;
return true;
}
RefPtr<AudioBus> AudioBus::createBufferFromRange(const AudioBus* sourceBuffer, unsigned startFrame, unsigned endFrame)
{
size_t numberOfSourceFrames = sourceBuffer->length();
unsigned numberOfChannels = sourceBuffer->numberOfChannels();
bool isRangeSafe = startFrame < endFrame && endFrame <= numberOfSourceFrames;
ASSERT(isRangeSafe);
if (!isRangeSafe)
return nullptr;
size_t rangeLength = endFrame - startFrame;
RefPtr<AudioBus> audioBus = create(numberOfChannels, rangeLength);
audioBus->setSampleRate(sourceBuffer->sampleRate());
for (unsigned i = 0; i < numberOfChannels; ++i)
audioBus->channel(i)->copyFromRange(sourceBuffer->channel(i), startFrame, endFrame);
return audioBus;
}
float AudioBus::maxAbsValue() const
{
float max = 0.0f;
for (unsigned i = 0; i < numberOfChannels(); ++i) {
const AudioChannel* channel = this->channel(i);
max = std::max(max, channel->maxAbsValue());
}
return max;
}
void AudioBus::normalize()
{
float max = maxAbsValue();
if (max)
scale(1.0f / max);
}
void AudioBus::scale(float scale)
{
for (unsigned i = 0; i < numberOfChannels(); ++i)
channel(i)->scale(scale);
}
void AudioBus::copyFromRange(const AudioBus& sourceBus, unsigned startFrame, unsigned endFrame)
{
if (!topologyMatches(sourceBus)) {
ASSERT_NOT_REACHED();
zero();
return;
}
size_t numberOfSourceFrames = sourceBus.length();
bool isRangeSafe = startFrame < endFrame && endFrame <= numberOfSourceFrames;
ASSERT(isRangeSafe);
if (!isRangeSafe) {
zero();
return;
}
unsigned numberOfChannels = this->numberOfChannels();
ASSERT(numberOfChannels <= MaxBusChannels);
if (numberOfChannels > MaxBusChannels) {
zero();
return;
}
for (unsigned i = 0; i < numberOfChannels; ++i)
channel(i)->copyFromRange(sourceBus.channel(i), startFrame, endFrame);
}
void AudioBus::copyFrom(const AudioBus& sourceBus, ChannelInterpretation channelInterpretation)
{
if (&sourceBus == this)
return;
zero();
sumFrom(sourceBus, channelInterpretation);
}
void AudioBus::sumFrom(const AudioBus& sourceBus, ChannelInterpretation channelInterpretation)
{
if (&sourceBus == this)
return;
unsigned numberOfSourceChannels = sourceBus.numberOfChannels();
unsigned numberOfDestinationChannels = numberOfChannels();
if (numberOfDestinationChannels == numberOfSourceChannels) {
for (unsigned i = 0; i < numberOfSourceChannels; ++i)
channel(i)->sumFrom(sourceBus.channel(i));
} else {
switch (channelInterpretation) {
case ChannelInterpretation::Speakers:
if (numberOfSourceChannels < numberOfDestinationChannels)
speakersSumFromByUpMixing(sourceBus);
else
speakersSumFromByDownMixing(sourceBus);
break;
case ChannelInterpretation::Discrete:
discreteSumFrom(sourceBus);
break;
default:
ASSERT_NOT_REACHED();
}
}
}
void AudioBus::speakersSumFromByUpMixing(const AudioBus& sourceBus)
{
unsigned numberOfSourceChannels = sourceBus.numberOfChannels();
unsigned numberOfDestinationChannels = numberOfChannels();
if ((numberOfSourceChannels == 1 && numberOfDestinationChannels == 2)
|| (numberOfSourceChannels == 1 && numberOfDestinationChannels == 4)) {
auto* sourceL = sourceBus.channelByType(ChannelLeft);
channelByType(ChannelLeft)->sumFrom(sourceL);
channelByType(ChannelRight)->sumFrom(sourceL);
} else if (numberOfSourceChannels == 1 && numberOfDestinationChannels == 6) {
channelByType(ChannelCenter)->sumFrom(sourceBus.channelByType(ChannelLeft));
} else if ((numberOfSourceChannels == 2 && numberOfDestinationChannels == 4)
|| (numberOfSourceChannels == 2 && numberOfDestinationChannels == 6)) {
channelByType(ChannelLeft)->sumFrom(sourceBus.channelByType(ChannelLeft));
channelByType(ChannelRight)->sumFrom(sourceBus.channelByType(ChannelRight));
} else if (numberOfSourceChannels == 4 && numberOfDestinationChannels == 6) {
channelByType(ChannelLeft)->sumFrom(sourceBus.channelByType(ChannelLeft));
channelByType(ChannelRight)->sumFrom(sourceBus.channelByType(ChannelRight));
channelByType(ChannelSurroundLeft)->sumFrom(sourceBus.channelByType(ChannelSurroundLeft));
channelByType(ChannelSurroundRight)->sumFrom(sourceBus.channelByType(ChannelSurroundRight));
} else {
discreteSumFrom(sourceBus);
}
}
void AudioBus::speakersSumFromByDownMixing(const AudioBus& sourceBus)
{
unsigned numberOfSourceChannels = sourceBus.numberOfChannels();
unsigned numberOfDestinationChannels = numberOfChannels();
if (numberOfSourceChannels == 2 && numberOfDestinationChannels == 1) {
AudioBus& sourceBusSafe = const_cast<AudioBus&>(sourceBus);
const float* sourceL = sourceBusSafe.channelByType(ChannelLeft)->data();
const float* sourceR = sourceBusSafe.channelByType(ChannelRight)->data();
float* destination = channelByType(ChannelLeft)->mutableData();
VectorMath::multiplyByScalarThenAddToOutput(sourceL, 0.5, destination, length());
VectorMath::multiplyByScalarThenAddToOutput(sourceR, 0.5, destination, length());
} else if (numberOfSourceChannels == 4 && numberOfDestinationChannels == 1) {
auto* sourceL = sourceBus.channelByType(ChannelLeft)->data();
auto* sourceR = sourceBus.channelByType(ChannelRight)->data();
auto* sourceSL = sourceBus.channelByType(ChannelSurroundLeft)->data();
auto* sourceSR = sourceBus.channelByType(ChannelSurroundRight)->data();
auto* destination = channelByType(ChannelLeft)->mutableData();
VectorMath::multiplyByScalarThenAddToOutput(sourceL, 0.25, destination, length());
VectorMath::multiplyByScalarThenAddToOutput(sourceR, 0.25, destination, length());
VectorMath::multiplyByScalarThenAddToOutput(sourceSL, 0.25, destination, length());
VectorMath::multiplyByScalarThenAddToOutput(sourceSR, 0.25, destination, length());
} else if (numberOfSourceChannels == 6 && numberOfDestinationChannels == 1) {
auto* sourceL = sourceBus.channelByType(ChannelLeft)->data();
auto* sourceR = sourceBus.channelByType(ChannelRight)->data();
auto* sourceC = sourceBus.channelByType(ChannelCenter)->data();
auto* sourceSL = sourceBus.channelByType(ChannelSurroundLeft)->data();
auto* sourceSR = sourceBus.channelByType(ChannelSurroundRight)->data();
auto* destination = channelByType(ChannelLeft)->mutableData();
float scaleSqrtHalf = sqrtf(0.5);
VectorMath::multiplyByScalarThenAddToOutput(sourceL, scaleSqrtHalf, destination, length());
VectorMath::multiplyByScalarThenAddToOutput(sourceR, scaleSqrtHalf, destination, length());
VectorMath::add(sourceC, destination, destination, length());
VectorMath::multiplyByScalarThenAddToOutput(sourceSL, 0.5, destination, length());
VectorMath::multiplyByScalarThenAddToOutput(sourceSR, 0.5, destination, length());
} else if (numberOfSourceChannels == 4 && numberOfDestinationChannels == 2) {
auto* sourceL = sourceBus.channelByType(ChannelLeft)->data();
auto* sourceR = sourceBus.channelByType(ChannelRight)->data();
auto* sourceSL = sourceBus.channelByType(ChannelSurroundLeft)->data();
auto* sourceSR = sourceBus.channelByType(ChannelSurroundRight)->data();
auto* destinationL = channelByType(ChannelLeft)->mutableData();
auto* destinationR = channelByType(ChannelRight)->mutableData();
VectorMath::multiplyByScalarThenAddToOutput(sourceL, 0.5, destinationL, length());
VectorMath::multiplyByScalarThenAddToOutput(sourceSL, 0.5, destinationL, length());
VectorMath::multiplyByScalarThenAddToOutput(sourceR, 0.5, destinationR, length());
VectorMath::multiplyByScalarThenAddToOutput(sourceSR, 0.5, destinationR, length());
} else if (numberOfSourceChannels == 6 && numberOfDestinationChannels == 2) {
auto* sourceL = sourceBus.channelByType(ChannelLeft)->data();
auto* sourceR = sourceBus.channelByType(ChannelRight)->data();
auto* sourceC = sourceBus.channelByType(ChannelCenter)->data();
auto* sourceSL = sourceBus.channelByType(ChannelSurroundLeft)->data();
auto* sourceSR = sourceBus.channelByType(ChannelSurroundRight)->data();
float* destinationL = channelByType(ChannelLeft)->mutableData();
float* destinationR = channelByType(ChannelRight)->mutableData();
float scaleSqrtHalf = sqrtf(0.5);
VectorMath::add(sourceL, destinationL, destinationL, length());
VectorMath::multiplyByScalarThenAddToOutput(sourceC, scaleSqrtHalf, destinationL, length());
VectorMath::multiplyByScalarThenAddToOutput(sourceSL, scaleSqrtHalf, destinationL, length());
VectorMath::add(sourceR, destinationR, destinationR, length());
VectorMath::multiplyByScalarThenAddToOutput(sourceC, scaleSqrtHalf, destinationR, length());
VectorMath::multiplyByScalarThenAddToOutput(sourceSR, scaleSqrtHalf, destinationR, length());
} else if (numberOfSourceChannels == 6 && numberOfDestinationChannels == 4) {
auto* sourceL = sourceBus.channelByType(ChannelLeft)->data();
auto* sourceR = sourceBus.channelByType(ChannelRight)->data();
auto* sourceC = sourceBus.channelByType(ChannelCenter)->data();
auto* destinationL = channelByType(ChannelLeft)->mutableData();
auto* destinationR = channelByType(ChannelRight)->mutableData();
auto scaleSqrtHalf = sqrtf(0.5);
VectorMath::add(sourceL, destinationL, destinationL, length());
VectorMath::multiplyByScalarThenAddToOutput(sourceC, scaleSqrtHalf, destinationL, length());
VectorMath::add(sourceR, destinationR, destinationR, length());
VectorMath::multiplyByScalarThenAddToOutput(sourceC, scaleSqrtHalf, destinationR, length());
channel(2)->sumFrom(sourceBus.channel(4));
channel(3)->sumFrom(sourceBus.channel(5));
} else {
discreteSumFrom(sourceBus);
}
}
void AudioBus::discreteSumFrom(const AudioBus& sourceBus)
{
unsigned numberOfSourceChannels = sourceBus.numberOfChannels();
unsigned numberOfDestinationChannels = numberOfChannels();
if (numberOfDestinationChannels < numberOfSourceChannels) {
for (unsigned i = 0; i < numberOfDestinationChannels; ++i)
channel(i)->sumFrom(sourceBus.channel(i));
} else if (numberOfDestinationChannels > numberOfSourceChannels) {
for (unsigned i = 0; i < numberOfSourceChannels; ++i)
channel(i)->sumFrom(sourceBus.channel(i));
}
}
void AudioBus::copyWithGainFrom(const AudioBus& sourceBus, float gain)
{
if (!topologyMatches(sourceBus)) {
ASSERT_NOT_REACHED();
zero();
return;
}
if (sourceBus.isSilent()) {
zero();
return;
}
unsigned numberOfChannels = this->numberOfChannels();
ASSERT(numberOfChannels <= MaxBusChannels);
if (numberOfChannels > MaxBusChannels)
return;
if (this == &sourceBus && gain == 1)
return;
AudioBus& sourceBusSafe = const_cast<AudioBus&>(sourceBus);
const float* sources[MaxBusChannels];
float* destinations[MaxBusChannels];
for (unsigned i = 0; i < numberOfChannels; ++i) {
sources[i] = sourceBusSafe.channel(i)->data();
destinations[i] = channel(i)->mutableData();
}
unsigned framesToProcess = length();
if (gain == 1) {
for (unsigned channelIndex = 0; channelIndex < numberOfChannels; ++channelIndex)
memcpy(destinations[channelIndex], sources[channelIndex], framesToProcess * sizeof(*destinations[channelIndex]));
} else if (!gain) {
for (unsigned channelIndex = 0; channelIndex < numberOfChannels; ++channelIndex)
memset(destinations[channelIndex], 0, framesToProcess * sizeof(*destinations[channelIndex]));
} else {
for (unsigned channelIndex = 0; channelIndex < numberOfChannels; ++channelIndex)
VectorMath::multiplyByScalar(sources[channelIndex], gain, destinations[channelIndex], framesToProcess);
}
}
void AudioBus::copyWithSampleAccurateGainValuesFrom(const AudioBus &sourceBus, float* gainValues, unsigned numberOfGainValues)
{
if (sourceBus.numberOfChannels() != 1 && !topologyMatches(sourceBus)) {
ASSERT_NOT_REACHED();
return;
}
if (!gainValues || numberOfGainValues > sourceBus.length()) {
ASSERT_NOT_REACHED();
return;
}
if (sourceBus.length() == numberOfGainValues && sourceBus.length() == length() && sourceBus.isSilent()) {
zero();
return;
}
const float* source = sourceBus.channel(0)->data();
for (unsigned channelIndex = 0; channelIndex < numberOfChannels(); ++channelIndex) {
if (sourceBus.numberOfChannels() == numberOfChannels())
source = sourceBus.channel(channelIndex)->data();
float* destination = channel(channelIndex)->mutableData();
VectorMath::multiply(source, gainValues, destination, numberOfGainValues);
}
}
RefPtr<AudioBus> AudioBus::createBySampleRateConverting(const AudioBus* sourceBus, bool mixToMono, double newSampleRate)
{
ASSERT(sourceBus && sourceBus->sampleRate());
if (!sourceBus || !sourceBus->sampleRate())
return nullptr;
double sourceSampleRate = sourceBus->sampleRate();
double destinationSampleRate = newSampleRate;
double sampleRateRatio = sourceSampleRate / destinationSampleRate;
unsigned numberOfSourceChannels = sourceBus->numberOfChannels();
if (numberOfSourceChannels == 1)
mixToMono = false;
if (sourceSampleRate == destinationSampleRate) {
if (mixToMono)
return AudioBus::createByMixingToMono(sourceBus);
return AudioBus::createBufferFromRange(sourceBus, 0, sourceBus->length());
}
if (sourceBus->isSilent()) {
RefPtr<AudioBus> silentBus = create(numberOfSourceChannels, sourceBus->length() / sampleRateRatio);
silentBus->setSampleRate(newSampleRate);
return silentBus;
}
const AudioBus* resamplerSourceBus;
RefPtr<AudioBus> mixedMonoBus;
if (mixToMono) {
mixedMonoBus = AudioBus::createByMixingToMono(sourceBus);
resamplerSourceBus = mixedMonoBus.get();
} else {
resamplerSourceBus = sourceBus;
}
int sourceLength = resamplerSourceBus->length();
int destinationLength = sourceLength / sampleRateRatio;
unsigned numberOfDestinationChannels = resamplerSourceBus->numberOfChannels();
RefPtr<AudioBus> destinationBus = create(numberOfDestinationChannels, destinationLength);
for (unsigned i = 0; i < numberOfDestinationChannels; ++i) {
const float* source = resamplerSourceBus->channel(i)->data();
float* destination = destinationBus->channel(i)->mutableData();
SincResampler::processBuffer(source, destination, sourceLength, sampleRateRatio);
}
destinationBus->clearSilentFlag();
destinationBus->setSampleRate(newSampleRate);
return destinationBus;
}
RefPtr<AudioBus> AudioBus::createByMixingToMono(const AudioBus* sourceBus)
{
if (sourceBus->isSilent())
return create(1, sourceBus->length());
switch (sourceBus->numberOfChannels()) {
case 1:
return AudioBus::createBufferFromRange(sourceBus, 0, sourceBus->length());
case 2:
{
unsigned n = sourceBus->length();
RefPtr<AudioBus> destinationBus = create(1, n);
const float* sourceL = sourceBus->channel(0)->data();
const float* sourceR = sourceBus->channel(1)->data();
float* destination = destinationBus->channel(0)->mutableData();
VectorMath::addVectorsThenMultiplyByScalar(sourceL, sourceR, 0.5, destination, n);
destinationBus->clearSilentFlag();
destinationBus->setSampleRate(sourceBus->sampleRate());
return destinationBus;
}
}
ASSERT_NOT_REACHED();
return nullptr;
}
bool AudioBus::isSilent() const
{
for (size_t i = 0; i < m_channels.size(); ++i) {
if (!m_channels[i]->isSilent())
return false;
}
return true;
}
void AudioBus::clearSilentFlag()
{
for (size_t i = 0; i < m_channels.size(); ++i)
m_channels[i]->clearSilentFlag();
}
}
#endif // ENABLE(WEB_AUDIO)