Cicm_Fft::Cicm_Fft(long aWindowSize)
{
m_window_size = Tools::clip_power_of_two(aWindowSize * 2);
m_array_size = m_window_size / 2;
m_order = Tools::log2(m_window_size);
m_scale = 1. / (double)(m_window_size * 2.);
#ifdef CICM_VDSP
vDSP_ctoz((Cicm_Complex *)aRealVector, 2, &m_input_complexes, 1, m_array_size);
vDSP_fft_zrip(m_fft_setup, &m_input_complexes, 1, m_log2_size, FFT_FORWARD);
#endif
#ifdef CICM_IPPS
Cicm_fft_get_size(m_order, &m_spec_size, &m_init_size, &m_buff_size);
m_fft_buff = Cicm_buffer_malloc(m_buff_size);
m_fft_init = Cicm_buffer_malloc(m_init_size);
m_fft_spec = Cicm_buffer_malloc(m_spec_size);
Cicm_fft_init_handle(&m_fft_handle, m_order, m_fft_spec, m_fft_init);
#endif
#ifdef CICM_FFTW_GSL
m_real_vector = (Cicm_Signal *)Cicm_signal_malloc(m_window_size);
m_input_complexes = (Cicm_Packed *)Cicm_packed_malloc(m_window_size);
m_output_complexes = (Cicm_Packed *)Cicm_packed_malloc(m_window_size);
m_handle_forward = Cicm_fft_init_handle_forward(m_window_size, m_real_vector, m_input_complexes);
m_handle_inverse = Cicm_fft_init_handle_inverse(m_window_size, m_input_complexes, m_real_vector);
#endif
}
void ZeroLatencyConvolver::setImpulseResponse(float* anImpulseResponse, long aSize)
{
int responseSize = 0;
m_number_of_convolutions = 0;
if(aSize > 0)
{
Cicm_Signal* impulseResponse;
if(aSize < 128) // FIR //
{
responseSize = 128;
m_number_of_convolutions = 1;
Cicm_signal_malloc(impulseResponse, responseSize);
for(int i = 0; i < responseSize; i++)
{
if(i < aSize)
impulseResponse[i] = anImpulseResponse[i];
else
impulseResponse[i] = 0.;
}
m_fir->setImpulseResponse(impulseResponse);
}
else if(aSize < 384) // FFT 128 //
{
responseSize = 384;
m_number_of_convolutions = 3;
Cicm_signal_malloc(impulseResponse, responseSize - 128);
for(int i = 128; i < responseSize - 128; i++)
{
if(i < aSize)
impulseResponse[i-128] = anImpulseResponse[i];
else
impulseResponse[i-128] = 0.;
}
m_fir->setImpulseResponse(anImpulseResponse);
m_fft_128->loadImpulseResponse(impulseResponse);
}
else if(aSize < 896) // FFT 256 //
{
responseSize = 896;
m_number_of_convolutions = 5;
Cicm_signal_malloc(impulseResponse, responseSize - 384);
for(int i = 384; i < responseSize - 384; i++)
{
if(i < aSize)
impulseResponse[i-384] = anImpulseResponse[i];
else
impulseResponse[i-384] = 0.;
}
m_fir->setImpulseResponse(anImpulseResponse);
m_fft_128->loadImpulseResponse(anImpulseResponse+128);
m_fft_256->loadImpulseResponse(impulseResponse);
}
else if(aSize < 1920) // FFT 512 //
{
responseSize = 1920;
m_number_of_convolutions = 7;
Cicm_signal_malloc(impulseResponse, responseSize - 896);
for(int i = 896; i < responseSize - 896; i++)
{
if(i < aSize)
impulseResponse[i-896] = anImpulseResponse[i];
else
impulseResponse[i-896] = 0.;
}
m_fir->setImpulseResponse(anImpulseResponse);
m_fft_128->loadImpulseResponse(anImpulseResponse+128);
m_fft_256->loadImpulseResponse(anImpulseResponse+384);
m_fft_512->loadImpulseResponse(impulseResponse);
}
else if(aSize < 3968) // FFT 1024 //
{
responseSize = 3968;
m_number_of_convolutions = 9;
Cicm_signal_malloc(impulseResponse, responseSize - 1920);
for(int i = 1920; i < responseSize - 1920; i++)
{
if(i < aSize)
impulseResponse[i-1920] = anImpulseResponse[i];
else
impulseResponse[i-1920] = 0.;
}
m_fir->setImpulseResponse(anImpulseResponse);
m_fft_128->loadImpulseResponse(anImpulseResponse+128);
m_fft_256->loadImpulseResponse(anImpulseResponse+284);
m_fft_512->loadImpulseResponse(anImpulseResponse+896);
m_fft_1024->loadImpulseResponse(impulseResponse);
}
else if(aSize < 8064) // FFT 2048 //
{
responseSize = 8064;
m_number_of_convolutions = 11;
Cicm_signal_malloc(impulseResponse, responseSize - 3968);
for(int i = 3968; i < responseSize - 3968; i++)
{
if(i < aSize)
impulseResponse[i-3968] = anImpulseResponse[i];
else
impulseResponse[i-3968] = 0.;
}
m_fir->setImpulseResponse(anImpulseResponse);
m_fft_128->loadImpulseResponse(anImpulseResponse+128);
m_fft_256->loadImpulseResponse(anImpulseResponse+284);
m_fft_512->loadImpulseResponse(anImpulseResponse+896);
m_fft_1024->loadImpulseResponse(anImpulseResponse+1920);
m_fft_2048->loadImpulseResponse(impulseResponse);
}
else if(aSize < 16256) // FFT 4096 //
{
responseSize = 16256;
m_number_of_convolutions = 13;
Cicm_signal_malloc(impulseResponse, responseSize - 8064);
for(int i = 8064; i < responseSize - 8064; i++)
{
if(i < aSize)
impulseResponse[i-8064] = anImpulseResponse[i];
else
impulseResponse[i-8064] = 0.;
}
m_fir->setImpulseResponse(anImpulseResponse);
m_fft_128->loadImpulseResponse(anImpulseResponse+128);
m_fft_256->loadImpulseResponse(anImpulseResponse+284);
m_fft_512->loadImpulseResponse(anImpulseResponse+896);
m_fft_1024->loadImpulseResponse(anImpulseResponse+1920);
m_fft_2048->loadImpulseResponse(anImpulseResponse+3968);
m_fft_4096->loadImpulseResponse(impulseResponse);
}
else if(aSize < 32640) // FFT 8192 //
{
responseSize = 32640;
m_number_of_convolutions = 15;
Cicm_signal_malloc(impulseResponse, responseSize - 16256);
for(int i = 16256; i < responseSize - 16256; i++)
{
if(i < aSize)
impulseResponse[i-16256] = anImpulseResponse[i];
else
impulseResponse[i-16256] = 0.;
}
m_fir->setImpulseResponse(anImpulseResponse);
m_fft_128->loadImpulseResponse(anImpulseResponse+128);
m_fft_256->loadImpulseResponse(anImpulseResponse+284);
m_fft_512->loadImpulseResponse(anImpulseResponse+896);
m_fft_1024->loadImpulseResponse(anImpulseResponse+1920);
m_fft_2048->loadImpulseResponse(anImpulseResponse+3968);
m_fft_4096->loadImpulseResponse(anImpulseResponse+8064);
m_fft_8192->loadImpulseResponse(impulseResponse);
}
else
{
responseSize = 32640;
m_number_of_convolutions = 15;
while(responseSize < aSize)
{
responseSize += 16384;
m_number_of_convolutions++;
}
if(m_fft_16384->getNumberOfInstances() != m_number_of_convolutions - 15)
{
delete m_fft_16384;
m_fft_16384 = new FftConvolution(16384, m_number_of_convolutions - 15);
}
Cicm_signal_malloc(impulseResponse, responseSize - 32640);
for(int i = 32640; i < responseSize - 32640; i++)
{
if(i < aSize)
impulseResponse[i-32640] = anImpulseResponse[i];
else
impulseResponse[i-32640] = 0.;
}
m_fir->setImpulseResponse(anImpulseResponse);
m_fft_128->loadImpulseResponse(anImpulseResponse+128);
m_fft_256->loadImpulseResponse(anImpulseResponse+284);
m_fft_512->loadImpulseResponse(anImpulseResponse+896);
m_fft_1024->loadImpulseResponse(anImpulseResponse+1920);
m_fft_2048->loadImpulseResponse(anImpulseResponse+3968);
m_fft_4096->loadImpulseResponse(anImpulseResponse+8064);
m_fft_8192->loadImpulseResponse(anImpulseResponse+16256);
m_fft_16384->loadImpulseResponse(impulseResponse);
}
Cicm_free(impulseResponse);
}
}
