int
main ()
{
int i;
int main_result;
word so[N];
word LARc[M];
main_result = 0;
for (i = 0; i < N; i++)
so[i] = inData[i];
Gsm_LPC_Analysis (so, LARc);
for (i = 0; i < N; i++)
main_result += (so[i] != outData[i]);
for (i = 0; i < M; i++)
main_result += (LARc[i] != outLARc[i]);
printf ("%d\n", main_result);
return main_result;
}
void Gsm_Coder (
struct gsm_state * State,
int16_t * s, /* [0..159] samples IN */
/*
* The RPE-LTD coder works on a frame by frame basis. The length of
* the frame is equal to 160 samples. Some computations are done
* once per frame to produce at the output of the coder the
* LARc [1..8] parameters which are the coded LAR coefficients and
* also to realize the inverse filtering operation for the entire
* frame (160 samples of signal d [0..159]). These parts produce at
* the output of the coder:
*/
int16_t * LARc, /* [0..7] LAR coefficients OUT */
/*
* Procedure 4.2.11 to 4.2.18 are to be executed four times per
* frame. That means once for each sub-segment RPE-LTP analysis of
* 40 samples. These parts produce at the output of the coder:
*/
int16_t *Nc, /* [0..3] LTP lag OUT */
int16_t *bc, /* [0..3] coded LTP gain OUT */
int16_t *Mc, /* [0..3] RPE grid selection OUT */
int16_t *xmaxc, /* [0..3] Coded maximum amplitude OUT */
int16_t *xMc /* [13*4] normalized RPE samples OUT */
)
{
int k ;
int16_t *dp = State->dp0 + 120 ; /* [-120...-1] */
int16_t *dpp = dp ; /* [0...39] */
int16_t so [160] ;
Gsm_Preprocess (State, s, so) ;
Gsm_LPC_Analysis (State, so, LARc) ;
Gsm_Short_Term_Analysis_Filter (State, LARc, so) ;
for (k = 0 ; k <= 3 ; k++, xMc += 13)
{ Gsm_Long_Term_Predictor (State,
so+k*40, /* d [0..39] IN */
dp, /* dp [-120..-1] IN */
State->e + 5, /* e [0..39] OUT */
dpp, /* dpp [0..39] OUT */
Nc++,
bc++) ;
Gsm_RPE_Encoding (/*-S,-*/
State->e + 5, /* e ][0..39][IN/OUT */
xmaxc++, Mc++, xMc) ;
/*
* Gsm_Update_of_reconstructed_short_time_residual_signal
* (dpp, State->e + 5, dp) ;
*/
{
register int i ;
for (i = 0 ; i <= 39 ; i++)
dp [i] = GSM_ADD (State->e [5 + i], dpp [i]) ;
}
dp += 40 ;
dpp += 40 ;
}
memcpy ((char *) State->dp0, (char *) (State->dp0 + 160),
120 * sizeof (*State->dp0)) ;
}