Home | History | Annotate | Download | only in src
      1 /*
      2  * Copyright 1992 by Jutta Degener and Carsten Bormann, Technische
      3  * Universitaet Berlin.  See the accompanying file "COPYRIGHT" for
      4  * details.  THERE IS ABSOLUTELY NO WARRANTY FOR THIS SOFTWARE.
      5  */
      6 
      7 /* $Header: /tmp_amd/presto/export/kbs/jutta/src/gsm/RCS/decode.c,v 1.1 1992/10/28 00:15:50 jutta Exp $ */
      8 
      9 #include <stdio.h>
     10 
     11 #include	"private.h"
     12 #include	"gsm.h"
     13 #include	"proto.h"
     14 
     15 /*
     16  *  4.3 FIXED POINT IMPLEMENTATION OF THE RPE-LTP DECODER
     17  */
     18 
     19 static void Postprocessing P2((S,s),
     20 	struct gsm_state	* S,
     21 	register word 		* s)
     22 {
     23 	register int		k;
     24 	register word		msr = S->msr;
     25 	register longword	ltmp;	/* for GSM_ADD */
     26 	register word		tmp;
     27 
     28 	for (k = 160; k--; s++) {
     29 		tmp = GSM_MULT_R( msr, 28180 );
     30 		msr = GSM_ADD(*s, tmp);  	   /* Deemphasis 	     */
     31 		*s  = GSM_ADD(msr, msr) & 0xFFF8;  /* Truncation & Upscaling */
     32 	}
     33 	S->msr = msr;
     34 }
     35 
     36 void Gsm_Decoder P8((S,LARcr, Ncr,bcr,Mcr,xmaxcr,xMcr,s),
     37 	struct gsm_state	* S,
     38 
     39 	word		* LARcr,	/* [0..7]		IN	*/
     40 
     41 	word		* Ncr,		/* [0..3] 		IN 	*/
     42 	word		* bcr,		/* [0..3]		IN	*/
     43 	word		* Mcr,		/* [0..3] 		IN 	*/
     44 	word		* xmaxcr,	/* [0..3]		IN 	*/
     45 	word		* xMcr,		/* [0..13*4]		IN	*/
     46 
     47 	word		* s)		/* [0..159]		OUT 	*/
     48 {
     49 	int		j, k;
     50 	word		erp[40], wt[160];
     51 	word		* drp = S->dp0 + 120;
     52 
     53 	for (j=0; j <= 3; j++, xmaxcr++, bcr++, Ncr++, Mcr++, xMcr += 13) {
     54 
     55 		Gsm_RPE_Decoding( S, *xmaxcr, *Mcr, xMcr, erp );
     56 		Gsm_Long_Term_Synthesis_Filtering( S, *Ncr, *bcr, erp, drp );
     57 
     58 		for (k = 0; k <= 39; k++) wt[ j * 40 + k ] =  drp[ k ];
     59 	}
     60 
     61 	Gsm_Short_Term_Synthesis_Filter( S, LARcr, wt, s );
     62 	Postprocessing(S, s);
     63 }
     64