furnace/extern/libsndfile-modified/src/GSM610/code.c
tildearrow 061991fe60 desubmodulize libsndfile - PLEASE READ
it appears a one-character typo in the cmake_minimum_required line prevents it from compiling under CMake 4.0.
in order to fix that, I had to take this thing out of submodules...

it is recommended to do this after you pull;

git submodule deinit extern/libsndfile
2025-04-02 15:09:53 -05:00

87 lines
2.3 KiB
C

/*
* Copyright 1992 by Jutta Degener and Carsten Bormann, Technische
* Universitaet Berlin. See the accompanying file "COPYRIGHT" for
* details. THERE IS ABSOLUTELY NO WARRANTY FOR THIS SOFTWARE.
*/
#include <stdlib.h>
#include <string.h>
#include "gsm610_priv.h"
/*
* 4.2 FIXED POINT IMPLEMENTATION OF THE RPE-LTP CODER
*/
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)) ;
}