FreeCalypso > hg > gsm-codec-lib

changeset 335:03198f6b0427

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
libtwamr: integrate d8_31pf.c
author Mychaela Falconia <falcon@freecalypso.org>
date Thu, 18 Apr 2024 22:35:05 +0000
parents dfaa8f322a8d
children 7f99b8ed30e5
files libtwamr/Makefile libtwamr/d8_31pf.c libtwamr/d8_31pf.h libtwamr/namespace.h
diffstat 4 files changed, 271 insertions(+), 3 deletions(-) [+]
line wrap: on
line diff
--- a/libtwamr/Makefile	Thu Apr 18 22:30:03 2024 +0000
+++ b/libtwamr/Makefile	Thu Apr 18 22:35:05 2024 +0000
@@ -3,9 +3,9 @@
 OBJS=	a_refl.o agc.o autocorr.o az_lsp.o b_cn_cod.o basicop2.o bgnscd.o \
 	bitno.o bits2prm.o c1035pf.o c2_11pf.o c2_9pf.o c3_14pf.o c4_17pf.o \
 	c8_31pf.o c_g_aver.o calc_cor.o calc_en.o cbsearch.o convolve.o cor_h.o\
-	d1035pf.o d2_11pf.o d2_9pf.o d3_14pf.o d4_17pf.o gmed_n.o graytab.o \
-	inv_sqrt.o log2.o oper_32b.o prmno.o s10_8pf.o set_sign.o sqrt_l.o \
-	tls_flags.o window.o
+	d1035pf.o d2_11pf.o d2_9pf.o d3_14pf.o d4_17pf.o d8_31pf.o gmed_n.o \
+	graytab.o inv_sqrt.o log2.o oper_32b.o prmno.o s10_8pf.o set_sign.o \
+	sqrt_l.o tls_flags.o window.o
 LIB=	libtwamr.a
 
 INSTALL_PREFIX=	/usr/local
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/libtwamr/d8_31pf.c	Thu Apr 18 22:35:05 2024 +0000
@@ -0,0 +1,226 @@
+/*
+********************************************************************************
+*
+*      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             : d8_31pf.c
+*      Purpose          : Builds the innovative codevector
+*
+********************************************************************************
+*/
+ 
+/*
+********************************************************************************
+*                         MODULE INCLUDE FILE AND VERSION ID
+********************************************************************************
+*/
+#include "namespace.h"
+#include "d8_31pf.h"
+ 
+/*
+********************************************************************************
+*                         INCLUDE FILES
+********************************************************************************
+*/
+#include "typedef.h"
+#include "basic_op.h"
+#include "no_count.h"
+#include "cnst.h"
+
+/*
+********************************************************************************
+*                         LOCAL VARIABLES AND TABLES
+********************************************************************************
+*/
+#define NB_PULSE  8           /* number of pulses  */
+
+/* define values/representation for output codevector and sign */
+#define POS_CODE  8191 
+#define NEG_CODE  8191 
+
+static void decompress10 (
+   Word16 MSBs,        /* i : MSB part of the index                 */
+   Word16 LSBs,        /* i : LSB part of the index                 */
+   Word16 index1,      /* i : index for first pos in pos_index[]    */ 
+   Word16 index2,      /* i : index for second pos in pos_index[]   */ 
+   Word16 index3,      /* i : index for third pos in pos_index[]    */ 
+   Word16 pos_indx[])  /* o : position of 3 pulses (decompressed)   */
+{
+   Word16 ia, ib, ic;
+
+   /*
+     pos_indx[index1] = ((MSBs-25*(MSBs/25))%5)*2 + (LSBs-4*(LSBs/4))%2;
+     pos_indx[index2] = ((MSBs-25*(MSBs/25))/5)*2 + (LSBs-4*(LSBs/4))/2;
+     pos_indx[index3] = (MSBs/25)*2 + LSBs/4;
+     */
+
+   test ();
+   if (sub(MSBs, 124) > 0)
+   {
+      MSBs = 124;                                              move16 (); 
+   }
+   
+   ia = mult(MSBs, 1311);
+   ia = sub(MSBs, extract_l(L_shr(L_mult(ia, 25), 1)));    
+   ib = shl(sub(ia, extract_l(L_shr(L_mult(mult(ia, 6554), 5), 1))), 1);
+   
+   ic = shl(shr(LSBs, 2), 2);
+   ic = sub(LSBs, ic);
+   pos_indx[index1] = add(ib, (ic & 1));                        logic16 ();
+   
+   ib = shl(mult(ia, 6554), 1);
+   pos_indx[index2] = add(ib, shr(ic, 1));
+   
+   pos_indx[index3] = add(shl(mult(MSBs, 1311), 1), shr(LSBs, 2));    
+
+   return;
+}    
+
+/*************************************************************************
+ *
+ *  FUNCTION:  decompress_code()
+ *
+ *  PURPOSE: decompression of the linear codewords to 4+three indeces  
+ *           one bit from each pulse is made robust to errors by 
+ *           minimizing the phase shift of a bit error.
+ *           4 signs (one for each track) 
+ *           i0,i4,i1 => one index (7+3) bits, 3   LSBs more robust
+ *           i2,i6,i5 => one index (7+3) bits, 3   LSBs more robust
+ *           i3,i7    => one index (5+2) bits, 2-3 LSbs more robust
+ *
+ *************************************************************************/
+static void decompress_code (
+    Word16 indx[],      /* i : position and sign of 8 pulses (compressed) */
+    Word16 sign_indx[], /* o : signs of 4 pulses (signs only)             */
+    Word16 pos_indx[]   /* o : position index of 8 pulses (position only) */
+)
+{
+    Word16 i, ia, ib, MSBs, LSBs, MSBs0_24;
+
+    for (i = 0; i < NB_TRACK_MR102; i++)
+    {
+       sign_indx[i] = indx[i];                                  move16 (); 
+    }
+    
+    /*
+      First index: 10x10x10 -> 2x5x2x5x2x5-> 125x2x2x2 -> 7+1x3 bits 
+      MSBs = indx[NB_TRACK]/8;
+      LSBs = indx[NB_TRACK]%8;
+      */
+    MSBs = shr(indx[NB_TRACK_MR102], 3);
+    LSBs = indx[NB_TRACK_MR102] & 7;                            logic16 ();
+    decompress10 (MSBs, LSBs, 0, 4, 1, pos_indx);               
+    
+    /*
+      Second index: 10x10x10 -> 2x5x2x5x2x5-> 125x2x2x2 -> 7+1x3 bits       
+      MSBs = indx[NB_TRACK+1]/8;
+      LSBs = indx[NB_TRACK+1]%8;
+      */
+    MSBs = shr(indx[NB_TRACK_MR102+1], 3);
+    LSBs = indx[NB_TRACK_MR102+1] & 7;                          logic16 ();
+    decompress10 (MSBs, LSBs, 2, 6, 5, pos_indx);               
+    
+    /*
+      Third index: 10x10 -> 2x5x2x5-> 25x2x2 -> 5+1x2 bits    
+      MSBs = indx[NB_TRACK+2]/4;
+      LSBs = indx[NB_TRACK+2]%4;
+      MSBs0_24 = (MSBs*25+12)/32;
+      if ((MSBs0_24/5)%2==1)
+         pos_indx[3] = (4-(MSBs0_24%5))*2 + LSBs%2;
+      else
+         pos_indx[3] = (MSBs0_24%5)*2 + LSBs%2;
+      pos_indx[7] = (MSBs0_24/5)*2 + LSBs/2;
+      */
+    MSBs = shr(indx[NB_TRACK_MR102+2], 2);
+    LSBs = indx[NB_TRACK_MR102+2] & 3;                          logic16 ();
+
+    MSBs0_24 = shr(add(extract_l(L_shr(L_mult(MSBs, 25), 1)), 12), 5);
+    
+    ia = mult(MSBs0_24, 6554) & 1;
+    ib = sub(MSBs0_24, extract_l(L_shr(L_mult(mult(MSBs0_24, 6554), 5), 1)));
+
+    test ();
+    if (sub(ia, 1) == 0)
+    {
+       ib = sub(4, ib);
+    }
+    pos_indx[3] = add(shl(ib, 1), (LSBs & 1));               logic16 ();           
+    
+    ia = shl(mult(MSBs0_24, 6554), 1);
+    pos_indx[7] = add(ia, shr(LSBs, 1));
+}
+
+/*
+********************************************************************************
+*                         PUBLIC PROGRAM CODE
+********************************************************************************
+*/
+/*************************************************************************
+ *
+ *  FUNCTION:   dec_8i40_31bits()
+ *
+ *  PURPOSE:  Builds the innovative codevector from the received
+ *            index of algebraic codebook.
+ *
+ *   See  c8_31pf.c  for more details about the algebraic codebook structure.
+ *
+ *************************************************************************/
+
+void dec_8i40_31bits (
+    Word16 index[],    /* i : index of 8 pulses (sign+position)         */
+    Word16 cod[]       /* o : algebraic (fixed) codebook excitation     */
+)
+{
+    Word16 i, j, pos1, pos2, sign;
+    Word16 linear_signs[NB_TRACK_MR102];
+    Word16 linear_codewords[NB_PULSE];
+    
+    for (i = 0; i < L_CODE; i++)
+    {
+        cod[i] = 0;                                    move16 (); 
+    }
+    
+    decompress_code (index, linear_signs, linear_codewords);
+    
+    /* decode the positions and signs of pulses and build the codeword */
+
+    for (j = 0; j < NB_TRACK_MR102; j++)
+    {
+       /* compute index i */
+       
+       i = linear_codewords[j];
+       i = extract_l (L_shr (L_mult (i, 4), 1));
+       pos1 = add (i, j);   /* position of pulse "j" */
+       
+       test (); 
+       if (linear_signs[j] == 0)
+       {
+          sign = POS_CODE;                             move16 (); /* +1.0 */
+       }
+       else
+       {
+          sign = -NEG_CODE;                            move16 (); /* -1.0 */
+       }
+       
+       cod[pos1] = sign;                               move16 (); 
+       
+       /* compute index i */
+       
+       i = linear_codewords[add (j, 4)];        
+       i = extract_l (L_shr (L_mult (i, 4), 1));
+       pos2 = add (i, j);      /* position of pulse "j+4" */
+       
+       test (); 
+       if (sub (pos2, pos1) < 0)
+       {
+          sign = negate (sign);
+       }
+       cod[pos2] = add (cod[pos2], sign);              move16 (); 
+    }
+    
+    return;
+}
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/libtwamr/d8_31pf.h	Thu Apr 18 22:35:05 2024 +0000
@@ -0,0 +1,41 @@
+/*
+********************************************************************************
+*
+*      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             : d8_31pf.h
+*      Purpose          : Builds the innovative codevector
+*
+********************************************************************************
+*/
+#ifndef d8_31pf_h
+#define d8_31pf_h "$Id $"
+ 
+/*
+********************************************************************************
+*                         INCLUDE FILES
+********************************************************************************
+*/
+#include "typedef.h"
+ 
+/*
+********************************************************************************
+*                         DEFINITION OF DATA TYPES
+********************************************************************************
+*/
+ 
+/*
+********************************************************************************
+*                         DECLARATION OF PROTOTYPES
+********************************************************************************
+*/ 
+void dec_8i40_31bits (
+    Word16 index[],    /* i : index of 8 pulses (sign+position)         */
+    Word16 cod[]       /* o : algebraic (fixed) codebook excitation     */
+);
+ 
+#endif
--- a/libtwamr/namespace.h	Thu Apr 18 22:30:03 2024 +0000
+++ b/libtwamr/namespace.h	Thu Apr 18 22:35:05 2024 +0000
@@ -112,6 +112,7 @@
 #define	decode_2i40_11bits	AMR__decode_2i40_11bits
 #define	decode_3i40_14bits	AMR__decode_3i40_14bits
 #define	decode_4i40_17bits	AMR__decode_4i40_17bits
+#define	dec_8i40_31bits		AMR__dec_8i40_31bits
 #define	dec_10i40_35bits	AMR__dec_10i40_35bits
 
 #endif	/* include guard */