FreeCalypso > hg > gsm-codec-lib

changeset 397:9b699f30e6f3

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
libtwamr: integrate pre_proc.c
author Mychaela Falconia <falcon@freecalypso.org>
date Mon, 06 May 2024 19:02:44 +0000
parents 38ee82480462
children df14b0c17e6d
files libtwamr/Makefile libtwamr/namespace.list libtwamr/pre_proc.c libtwamr/pre_proc.h
diffstat 4 files changed, 183 insertions(+), 4 deletions(-) [+]
line wrap: on
line diff
--- 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
 
--- 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
--- /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;
+}
--- /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