#include "config.h"
#if ENABLE(WEB_AUDIO)
#if USE(WEBAUDIO_IPP)
#include "FFTFrame.h"
#include "VectorMath.h"
#include <wtf/MathExtras.h>
namespace WebCore {
const unsigned maximumFFTPower2Size = 24;
FFTFrame::FFTFrame(unsigned fftSize)
: m_FFTSize(fftSize)
, m_log2FFTSize(static_cast<unsigned>(log2(fftSize)))
, m_complexData(fftSize)
, m_realData(fftSize / 2)
, m_imagData(fftSize / 2)
{
ASSERT(1UL << m_log2FFTSize == m_FFTSize);
ASSERT(m_log2FFTSize <= maximumFFTPower2Size);
ippsDFTInitAlloc_R_32f(&m_DFTSpec, m_FFTSize, IPP_FFT_NODIV_BY_ANY, ippAlgHintFast);
int bufferSize = 0;
ippsDFTGetBufSize_R_32f(m_DFTSpec, &bufferSize);
m_buffer = ippsMalloc_8u(bufferSize);
}
FFTFrame::FFTFrame()
: m_FFTSize(0)
, m_log2FFTSize(0)
{
}
FFTFrame::FFTFrame(const FFTFrame& frame)
: m_FFTSize(frame.m_FFTSize)
, m_log2FFTSize(frame.m_log2FFTSize)
, m_complexData(frame.m_FFTSize)
, m_realData(frame.m_FFTSize / 2)
, m_imagData(frame.m_FFTSize / 2)
{
ippsDFTInitAlloc_R_32f(&m_DFTSpec, m_FFTSize, IPP_FFT_NODIV_BY_ANY, ippAlgHintFast);
int bufferSize = 0;
ippsDFTGetBufSize_R_32f(m_DFTSpec, &bufferSize);
m_buffer = ippsMalloc_8u(bufferSize);
unsigned numberOfBytes = sizeof(float) * m_FFTSize;
memcpy(realData(), frame.realData(), numberOfBytes);
memcpy(imagData(), frame.imagData(), numberOfBytes);
}
void FFTFrame::initialize()
{
}
void FFTFrame::cleanup()
{
}
FFTFrame::~FFTFrame()
{
ippsFree(m_buffer);
ippsDFTFree_R_32f(m_DFTSpec);
}
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();
unsigned halfSize = fftSize() / 2;
float real0 = realP1[0];
float imag0 = imagP1[0];
VectorMath::zvmul(realP1, imagP1, realP2, imagP2, realP1, imagP1, halfSize);
realP1[0] = real0 * realP2[0];
imagP1[0] = imag0 * imagP2[0];
float scale = 0.5f;
VectorMath::vsmul(realP1, 1, &scale, realP1, 1, halfSize);
VectorMath::vsmul(imagP1, 1, &scale, imagP1, 1, halfSize);
}
void FFTFrame::doFFT(const float* data)
{
Ipp32f* complexP = m_complexData.data();
ippsDFTFwd_RToPerm_32f(reinterpret_cast<Ipp32f*>(const_cast<float*>(data)), complexP, m_DFTSpec, m_buffer);
const Ipp32f scale = 2.0f;
ippsMulC_32f_I(scale, complexP, m_FFTSize);
Ipp32f* realP = m_realData.data();
Ipp32f* imagP = m_imagData.data();
ippsCplxToReal_32fc(reinterpret_cast<Ipp32fc*>(complexP), realP, imagP, m_FFTSize >> 1);
}
void FFTFrame::doInverseFFT(float* data)
{
Ipp32f* complexP = getUpToDateComplexData();
ippsDFTInv_PermToR_32f(complexP, reinterpret_cast<Ipp32f*>(data), m_DFTSpec, m_buffer);
const float scale = 1.0 / (2 * m_FFTSize);
ippsMulC_32f_I(scale, reinterpret_cast<Ipp32f*>(data), m_FFTSize);
}
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 >> 1;
Ipp32f* realP = m_realData.data();
Ipp32f* imagP = m_imagData.data();
Ipp32fc* complexP = reinterpret_cast<Ipp32fc*>(m_complexData.data());
ippsRealToCplx_32f(realP, imagP, complexP, len);
return const_cast<float*>(m_complexData.data());
}
}
#endif // USE(WEBAUDIO_IPP)
#endif // ENABLE(WEB_AUDIO)