#include "config.h"
#if ENABLE(WEB_AUDIO)
#if !OS(DARWIN) && USE(WEBAUDIO_MKL)
#include "FFTFrame.h"
#include "mkl_vml.h"
#include <wtf/MathExtras.h>
namespace {
DFTI_DESCRIPTOR_HANDLE createDescriptorHandle(int fftSize)
{
DFTI_DESCRIPTOR_HANDLE handle = 0;
MKL_LONG status = DftiCreateDescriptor(&handle, DFTI_SINGLE, DFTI_REAL, 1, fftSize);
ASSERT(DftiErrorClass(status, DFTI_NO_ERROR));
status = DftiSetValue(handle, DFTI_PLACEMENT, DFTI_NOT_INPLACE);
ASSERT(DftiErrorClass(status, DFTI_NO_ERROR));
status = DftiSetValue(handle, DFTI_PACKED_FORMAT, DFTI_PERM_FORMAT);
ASSERT(DftiErrorClass(status, DFTI_NO_ERROR));
status = DftiSetValue(handle, DFTI_FORWARD_SCALE, 2.0);
ASSERT(DftiErrorClass(status, DFTI_NO_ERROR));
double scale = 1.0 / (2.0 * fftSize);
status = DftiSetValue(handle, DFTI_BACKWARD_SCALE, scale);
ASSERT(DftiErrorClass(status, DFTI_NO_ERROR));
status = DftiCommitDescriptor(handle);
ASSERT(DftiErrorClass(status, DFTI_NO_ERROR));
return handle;
}
}
namespace WebCore {
const int kMaxFFTPow2Size = 24;
DFTI_DESCRIPTOR_HANDLE* FFTFrame::descriptorHandles = 0;
FFTFrame::FFTFrame(unsigned fftSize)
: m_FFTSize(fftSize)
, m_log2FFTSize(static_cast<unsigned>(log2(fftSize)))
, m_handle(0)
, m_complexData(fftSize)
, m_realData(fftSize / 2)
, m_imagData(fftSize / 2)
{
ASSERT(1UL << m_log2FFTSize == m_FFTSize);
m_handle = descriptorHandleForSize(fftSize);
}
FFTFrame::FFTFrame()
: m_FFTSize(0)
, m_log2FFTSize(0)
, m_handle(0)
{
}
FFTFrame::FFTFrame(const FFTFrame& frame)
: m_FFTSize(frame.m_FFTSize)
, m_log2FFTSize(frame.m_log2FFTSize)
, m_handle(0)
, m_complexData(frame.m_FFTSize)
, m_realData(frame.m_FFTSize / 2)
, m_imagData(frame.m_FFTSize / 2)
{
m_handle = descriptorHandleForSize(m_FFTSize);
unsigned nbytes = sizeof(float) * (m_FFTSize / 2);
memcpy(realData(), frame.realData(), nbytes);
memcpy(imagData(), frame.imagData(), nbytes);
}
FFTFrame::~FFTFrame()
{
}
void FFTFrame::multiply(const FFTFrame& frame)
{
FFTFrame& frame1 = *this;
FFTFrame& frame2 = const_cast<FFTFrame&>(frame);
float* realP1 = frame1.realData();
float* imagP1 = frame1.imagData();
const float* realP2 = frame2.realData();
const float* imagP2 = frame2.imagData();
float scale = 0.5f;
realP1[0] *= scale * realP2[0];
imagP1[0] *= scale * imagP2[0];
float* interleavedData1 = frame1.getUpToDateComplexData();
float* interleavedData2 = frame2.getUpToDateComplexData();
unsigned halfSize = m_FFTSize / 2;
vcMul(halfSize - 1,
reinterpret_cast<MKL_Complex8*>(interleavedData1) + 1,
reinterpret_cast<MKL_Complex8*>(interleavedData2) + 1,
reinterpret_cast<MKL_Complex8*>(interleavedData1) + 1);
for (unsigned i = 1; i < halfSize; ++i) {
int baseComplexIndex = 2 * i;
realP1[i] = scale * interleavedData1[baseComplexIndex];
imagP1[i] = scale * interleavedData1[baseComplexIndex + 1];
}
}
void FFTFrame::doFFT(float* data)
{
MKL_LONG status = DftiComputeForward(m_handle, data, m_complexData.data());
ASSERT_UNUSED(status, DftiErrorClass(status, DFTI_NO_ERROR));
int len = m_FFTSize / 2;
for (int i = 0; i < len; ++i) {
int baseComplexIndex = 2 * i;
m_realData[i] = m_complexData[baseComplexIndex];
m_imagData[i] = m_complexData[baseComplexIndex + 1];
}
}
void FFTFrame::doInverseFFT(float* data)
{
float* interleavedData = getUpToDateComplexData();
MKL_LONG status = DftiComputeBackward(m_handle, interleavedData, data);
ASSERT_UNUSED(status, DftiErrorClass(status, DFTI_NO_ERROR));
}
void FFTFrame::initialize()
{
}
void FFTFrame::cleanup()
{
if (!descriptorHandles)
return;
for (int i = 0; i < kMaxFFTPow2Size; ++i) {
if (descriptorHandles[i]) {
MKL_LONG status = DftiFreeDescriptor(&descriptorHandles[i]);
ASSERT_UNUSED(status, DftiErrorClass(status, DFTI_NO_ERROR));
}
}
delete[] descriptorHandles;
descriptorHandles = 0;
}
float* FFTFrame::realData() const
{
return const_cast<float*>(m_realData.data());
}
float* FFTFrame::imagData() const
{
return const_cast<float*>(m_imagData.data());
}
float* FFTFrame::getUpToDateComplexData()
{
int len = m_FFTSize / 2;
for (int i = 0; i < len; ++i) {
int baseComplexIndex = 2 * i;
m_complexData[baseComplexIndex] = m_realData[i];
m_complexData[baseComplexIndex + 1] = m_imagData[i];
}
return const_cast<float*>(m_complexData.data());
}
DFTI_DESCRIPTOR_HANDLE FFTFrame::descriptorHandleForSize(unsigned fftSize)
{
if (!descriptorHandles) {
descriptorHandles = new DFTI_DESCRIPTOR_HANDLE[kMaxFFTPow2Size];
for (int i = 0; i < kMaxFFTPow2Size; ++i)
descriptorHandles[i] = 0;
}
ASSERT(fftSize);
int pow2size = static_cast<int>(log2(fftSize));
ASSERT(pow2size < kMaxFFTPow2Size);
if (!descriptorHandles[pow2size])
descriptorHandles[pow2size] = createDescriptorHandle(fftSize);
return descriptorHandles[pow2size];
}
}
#endif // !OS(DARWIN) && USE(WEBAUDIO_MKL)
#endif // ENABLE(WEB_AUDIO)