You cannot select more than 25 topics Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.

148 lines
3.2 KiB
C++

/*
* Project 25 IMBE Encoder/Decoder Fixed-Point implementation
* Developed by Pavel Yazev E-mail: pyazev@gmail.com
* Version 1.0 (c) Copyright 2009
*
* This is free software; you can redistribute it and/or modify it
* under the terms of the GNU General Public License as published by
* the Free Software Foundation; either version 3, or (at your option)
* any later version.
*
* The software is distributed in the hope that it will be useful, but
* WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
* See the GNU General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with this; see the file COPYING. If not, write to the Free
* Software Foundation, Inc., 51 Franklin Street, Boston, MA
* 02110-1301, USA.
*/
#include "typedef.h"
#include "basic_op.h"
#include "imbe.h"
#include "aux_sub.h"
#include "dsp_sub.h"
#include "math_sub.h"
#include "uv_synt.h"
#include "rand_gen.h"
#include "tbls.h"
#include "encode.h"
#include "imbe_vocoder_impl.h"
void imbe_vocoder_impl::uv_synt_init(void)
{
fft_init();
v_zap(uv_mem, 105);
}
void imbe_vocoder_impl::uv_synt(IMBE_PARAM *imbe_param, Word16 *snd)
{
Cmplx16 Uw[FFTLENGTH];
Word16 i, index_a, index_b, index_aux, ha, hb, *v_uv_dsn_ptr, *sa_ptr, sa;
Word32 fund_freq, fund_freq_2, fund_freq_acc_a, fund_freq_acc_b;
sa_ptr = imbe_param->sa;
v_uv_dsn_ptr = imbe_param->v_uv_dsn;
fund_freq = imbe_param->fund_freq;
fund_freq_2 = L_shr(fund_freq, 1);
v_zap((Word16 *)&Uw, 2 * FFTLENGTH);
fund_freq_acc_a = L_sub(fund_freq, fund_freq_2);
fund_freq_acc_b = L_add(fund_freq, fund_freq_2);
for(i = 0; i < imbe_param->num_harms; i++)
{
ha = extract_h(fund_freq_acc_a);
hb = extract_h(fund_freq_acc_b);
index_a = (ha >> 8) + ((ha & 0xFF)?1:0);
index_b = (hb >> 8) + ((hb & 0xFF)?1:0);
index_aux = 256 - index_a;
sa = shl(*sa_ptr, 3);
if(*v_uv_dsn_ptr++)
{
while(index_a < index_b)
{
Uw[index_a].re = Uw[index_a].im = 0;
Uw[index_aux].re = Uw[index_aux].im = 0;
index_a++;
index_aux--;
}
}
else
{
while(index_a < index_b)
{
Uw[index_a].re = mult(sa, rand_gen());
Uw[index_a].im = mult(sa, rand_gen());
//Uw[index_a].re = sa;
//Uw[index_a].im = sa;
Uw[index_aux].re = Uw[index_a].re;
Uw[index_aux].im = negate(Uw[index_a].im);
index_a++;
index_aux--;
}
}
fund_freq_acc_a = L_add(fund_freq_acc_a, fund_freq);
fund_freq_acc_b = L_add(fund_freq_acc_b, fund_freq);
sa_ptr++;
}
/*
j = 128;
for(i = 0; i < 128; i++)
Uw_tmp[j++] = Uw[i];
j = 128;
for(i = 0; i < 128; i++)
Uw_tmp[i] = Uw[j++];
for(i = 0; i < 256; i++)
Uw[i] = Uw_tmp[i];
*/
fft((Word16 *)&Uw, FFTLENGTH, -1);
for(i = 0; i < 105; i++)
snd[i] = uv_mem[i];
index_aux = 73;
for(i = 105; i < FRAME; i++)
snd[i] = shl(Uw[index_aux++].re, 3);
// Weighted Overlap Add Algorithm
index_aux = 24;
index_a = 0;
index_b = 48;
for(i = 56; i < 105; i++)
{
snd[i] = extract_h(L_add(L_mult(snd[i], ws[index_b]), L_mult(shl(Uw[index_aux].re, 3), ws[index_a])));
index_aux++;
index_a++;
index_b--;
}
index_aux = 128;
for(i = 0; i < 105; i++)
uv_mem[i] = shl(Uw[index_aux++].re, 3);
}