# HG changeset patch # User Mychaela Falconia # Date 1715022164 0 # Node ID 9b699f30e6f3d1f72bf73b084ce1926b40aecdbd # Parent 38ee824804628633f8379cabb945d5fe538dc8ca libtwamr: integrate pre_proc.c diff -r 38ee82480462 -r 9b699f30e6f3 libtwamr/Makefile --- a/libtwamr/Makefile Mon May 06 18:57:05 2024 +0000 +++ b/libtwamr/Makefile Mon May 06 19:02:44 2024 +0000 @@ -10,10 +10,10 @@ gc_pred.o gmed_n.o graytab.o hp_max.o int_lpc.o int_lsf.o inter_36.o \ inv_sqrt.o lag_wind.o levinson.o log2.o lpc.o lsfwt.o lsp.o lsp_avg.o \ lsp_az.o lsp_lsf.o lsp_tab.o mac_32.o oper_32b.o ph_disp.o pitch_fr.o \ - post_pro.o pow2.o pre_big.o prmno.o q_gain_c.o q_gain_p.o q_plsf.o \ - q_plsf3_tab.o q_plsf5_tab.o q_plsf_3.o q_plsf_5.o qgain475.o qgain795.o\ - qua_gain.o qua_gain_tab.o reorder.o residu.o s10_8pf.o set_sign.o \ - sqrt_l.o syn_filt.o tls_flags.o weight_a.o window.o + post_pro.o pow2.o pre_big.o pre_proc.o prmno.o q_gain_c.o q_gain_p.o \ + q_plsf.o q_plsf3_tab.o q_plsf5_tab.o q_plsf_3.o q_plsf_5.o qgain475.o \ + qgain795.o qua_gain.o qua_gain_tab.o reorder.o residu.o s10_8pf.o \ + set_sign.o sqrt_l.o syn_filt.o tls_flags.o weight_a.o window.o HDRS= namespace.h LIB= libtwamr.a diff -r 38ee82480462 -r 9b699f30e6f3 libtwamr/namespace.list --- a/libtwamr/namespace.list Mon May 06 18:57:05 2024 +0000 +++ b/libtwamr/namespace.list Mon May 06 19:02:44 2024 +0000 @@ -28,6 +28,7 @@ Lag_window Levinson Levinson_reset Lsf_lsp Lsp_lsf Reorder_lsf Lsf_wt Lsp_Az Pitch_fr Pitch_fr_reset +Pre_Process Pre_Process_reset Post_Process Post_Process_reset Q_plsf_reset Q_plsf_3 Q_plsf_5 Qua_gain Residu Syn_filt Weight_Ai diff -r 38ee82480462 -r 9b699f30e6f3 libtwamr/pre_proc.c --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/libtwamr/pre_proc.c Mon May 06 19:02:44 2024 +0000 @@ -0,0 +1,116 @@ +/* +******************************************************************************** +* +* GSM AMR-NB speech codec R98 Version 7.6.0 December 12, 2001 +* R99 Version 3.3.0 +* REL-4 Version 4.1.0 +* +******************************************************************************** +* +* File : pre_proc.c +* Purpose : Preprocessing of input speech. +* +******************************************************************************** +*/ + +/* +******************************************************************************** +* MODULE INCLUDE FILE AND VERSION ID +******************************************************************************** +*/ +#include "namespace.h" +#include "pre_proc.h" + +/* +******************************************************************************** +* INCLUDE FILES +******************************************************************************** +*/ +#include "typedef.h" +#include "basic_op.h" +#include "oper_32b.h" +#include "no_count.h" + +/* +******************************************************************************** +* LOCAL VARIABLES AND TABLES +******************************************************************************** +*/ +/* filter coefficients (fc = 80 Hz, coeff. b[] is divided by 2) */ +static const Word16 b[3] = {1899, -3798, 1899}; +static const Word16 a[3] = {4096, 7807, -3733}; + +/* +******************************************************************************** +* PUBLIC PROGRAM CODE +******************************************************************************** +*/ + +/************************************************************************* +* +* Function: Pre_Process_reset +* Purpose: Initializes state memory to zero +* +************************************************************************** +*/ +void Pre_Process_reset (Pre_ProcessState *state) +{ + state->y2_hi = 0; + state->y2_lo = 0; + state->y1_hi = 0; + state->y1_lo = 0; + state->x0 = 0; + state->x1 = 0; +} + +/************************************************************************* + * + * FUNCTION: Pre_Process() + * + * PURPOSE: Preprocessing of input speech. + * + * DESCRIPTION: + * - 2nd order high pass filtering with cut off frequency at 80 Hz. + * - Divide input by two. + * + * + * Algorithm: + * + * y[i] = b[0]*x[i]/2 + b[1]*x[i-1]/2 + b[2]*x[i-2]/2 + * + a[1]*y[i-1] + a[2]*y[i-2]; + * + * + * Input is divided by two in the filtering process. + * + *************************************************************************/ +int Pre_Process ( + Pre_ProcessState *st, + Word16 signal[], /* input/output signal */ + Word16 lg) /* lenght of signal */ +{ + Word16 i, x2; + Word32 L_tmp; + + for (i = 0; i < lg; i++) + { + x2 = st->x1; move16 (); + st->x1 = st->x0; move16 (); + st->x0 = signal[i]; move16 (); + + /* y[i] = b[0]*x[i]/2 + b[1]*x[i-1]/2 + b140[2]*x[i-2]/2 */ + /* + a[1]*y[i-1] + a[2] * y[i-2]; */ + + L_tmp = Mpy_32_16 (st->y1_hi, st->y1_lo, a[1]); + L_tmp = L_add (L_tmp, Mpy_32_16 (st->y2_hi, st->y2_lo, a[2])); + L_tmp = L_mac (L_tmp, st->x0, b[0]); + L_tmp = L_mac (L_tmp, st->x1, b[1]); + L_tmp = L_mac (L_tmp, x2, b[2]); + L_tmp = L_shl (L_tmp, 3); + signal[i] = round (L_tmp); move16 (); + + st->y2_hi = st->y1_hi; move16 (); + st->y2_lo = st->y1_lo; move16 (); + L_Extract (L_tmp, &st->y1_hi, &st->y1_lo); + } + return 0; +} diff -r 38ee82480462 -r 9b699f30e6f3 libtwamr/pre_proc.h --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/libtwamr/pre_proc.h Mon May 06 19:02:44 2024 +0000 @@ -0,0 +1,62 @@ +/* +******************************************************************************** +* +* GSM AMR-NB speech codec R98 Version 7.6.0 December 12, 2001 +* R99 Version 3.3.0 +* REL-4 Version 4.1.0 +* +******************************************************************************** +* +* File : pre_proc.h +* Purpose : Preprocessing of input speech. +* +******************************************************************************** +*/ +#ifndef pre_proc_h +#define pre_proc_h "$Id $" + +/* +******************************************************************************** +* INCLUDE FILES +******************************************************************************** +*/ +#include "typedef.h" + +/* +******************************************************************************** +* LOCAL VARIABLES AND TABLES +******************************************************************************** +*/ + +/* +******************************************************************************** +* DEFINITION OF DATA TYPES +******************************************************************************** +*/ +typedef struct { + Word16 y2_hi; + Word16 y2_lo; + Word16 y1_hi; + Word16 y1_lo; + Word16 x0; + Word16 x1; +} Pre_ProcessState; + +/* +******************************************************************************** +* DECLARATION OF PROTOTYPES +******************************************************************************** +*/ + +void Pre_Process_reset (Pre_ProcessState *st); +/* reset of pre processing state (i.e. set state memory to zero) + returns 0 on success + */ + +int Pre_Process ( + Pre_ProcessState *st, + Word16 signal[], /* Input/output signal */ + Word16 lg /* Lenght of signal */ +); + +#endif