# HG changeset patch # User Mychaela Falconia # Date 1713479705 0 # Node ID 03198f6b042721021b49cc11ac72933a6c676f82 # Parent dfaa8f322a8d540cc1b2575b0130ced063f7bb89 libtwamr: integrate d8_31pf.c diff -r dfaa8f322a8d -r 03198f6b0427 libtwamr/Makefile --- 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 diff -r dfaa8f322a8d -r 03198f6b0427 libtwamr/d8_31pf.c --- /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; +} diff -r dfaa8f322a8d -r 03198f6b0427 libtwamr/d8_31pf.h --- /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 diff -r dfaa8f322a8d -r 03198f6b0427 libtwamr/namespace.h --- 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 */