1*a58d3d2aSXin Li /***********************************************************************
2*a58d3d2aSXin Li Copyright (c) 2006-2011, Skype Limited. All rights reserved.
3*a58d3d2aSXin Li Redistribution and use in source and binary forms, with or without
4*a58d3d2aSXin Li modification, are permitted provided that the following conditions
5*a58d3d2aSXin Li are met:
6*a58d3d2aSXin Li - Redistributions of source code must retain the above copyright notice,
7*a58d3d2aSXin Li this list of conditions and the following disclaimer.
8*a58d3d2aSXin Li - Redistributions in binary form must reproduce the above copyright
9*a58d3d2aSXin Li notice, this list of conditions and the following disclaimer in the
10*a58d3d2aSXin Li documentation and/or other materials provided with the distribution.
11*a58d3d2aSXin Li - Neither the name of Internet Society, IETF or IETF Trust, nor the
12*a58d3d2aSXin Li names of specific contributors, may be used to endorse or promote
13*a58d3d2aSXin Li products derived from this software without specific prior written
14*a58d3d2aSXin Li permission.
15*a58d3d2aSXin Li THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
16*a58d3d2aSXin Li AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
17*a58d3d2aSXin Li IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
18*a58d3d2aSXin Li ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
19*a58d3d2aSXin Li LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
20*a58d3d2aSXin Li CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
21*a58d3d2aSXin Li SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
22*a58d3d2aSXin Li INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
23*a58d3d2aSXin Li CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
24*a58d3d2aSXin Li ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
25*a58d3d2aSXin Li POSSIBILITY OF SUCH DAMAGE.
26*a58d3d2aSXin Li ***********************************************************************/
27*a58d3d2aSXin Li
28*a58d3d2aSXin Li #ifdef HAVE_CONFIG_H
29*a58d3d2aSXin Li #include "config.h"
30*a58d3d2aSXin Li #endif
31*a58d3d2aSXin Li
32*a58d3d2aSXin Li #include "main.h"
33*a58d3d2aSXin Li
34*a58d3d2aSXin Li /* Encode side-information parameters to payload */
silk_encode_indices(silk_encoder_state * psEncC,ec_enc * psRangeEnc,opus_int FrameIndex,opus_int encode_LBRR,opus_int condCoding)35*a58d3d2aSXin Li void silk_encode_indices(
36*a58d3d2aSXin Li silk_encoder_state *psEncC, /* I/O Encoder state */
37*a58d3d2aSXin Li ec_enc *psRangeEnc, /* I/O Compressor data structure */
38*a58d3d2aSXin Li opus_int FrameIndex, /* I Frame number */
39*a58d3d2aSXin Li opus_int encode_LBRR, /* I Flag indicating LBRR data is being encoded */
40*a58d3d2aSXin Li opus_int condCoding /* I The type of conditional coding to use */
41*a58d3d2aSXin Li )
42*a58d3d2aSXin Li {
43*a58d3d2aSXin Li opus_int i, k, typeOffset;
44*a58d3d2aSXin Li opus_int encode_absolute_lagIndex, delta_lagIndex;
45*a58d3d2aSXin Li opus_int16 ec_ix[ MAX_LPC_ORDER ];
46*a58d3d2aSXin Li opus_uint8 pred_Q8[ MAX_LPC_ORDER ];
47*a58d3d2aSXin Li const SideInfoIndices *psIndices;
48*a58d3d2aSXin Li
49*a58d3d2aSXin Li if( encode_LBRR ) {
50*a58d3d2aSXin Li psIndices = &psEncC->indices_LBRR[ FrameIndex ];
51*a58d3d2aSXin Li } else {
52*a58d3d2aSXin Li psIndices = &psEncC->indices;
53*a58d3d2aSXin Li }
54*a58d3d2aSXin Li
55*a58d3d2aSXin Li /*******************************************/
56*a58d3d2aSXin Li /* Encode signal type and quantizer offset */
57*a58d3d2aSXin Li /*******************************************/
58*a58d3d2aSXin Li typeOffset = 2 * psIndices->signalType + psIndices->quantOffsetType;
59*a58d3d2aSXin Li celt_assert( typeOffset >= 0 && typeOffset < 6 );
60*a58d3d2aSXin Li celt_assert( encode_LBRR == 0 || typeOffset >= 2 );
61*a58d3d2aSXin Li if( encode_LBRR || typeOffset >= 2 ) {
62*a58d3d2aSXin Li ec_enc_icdf( psRangeEnc, typeOffset - 2, silk_type_offset_VAD_iCDF, 8 );
63*a58d3d2aSXin Li } else {
64*a58d3d2aSXin Li ec_enc_icdf( psRangeEnc, typeOffset, silk_type_offset_no_VAD_iCDF, 8 );
65*a58d3d2aSXin Li }
66*a58d3d2aSXin Li
67*a58d3d2aSXin Li /****************/
68*a58d3d2aSXin Li /* Encode gains */
69*a58d3d2aSXin Li /****************/
70*a58d3d2aSXin Li /* first subframe */
71*a58d3d2aSXin Li if( condCoding == CODE_CONDITIONALLY ) {
72*a58d3d2aSXin Li /* conditional coding */
73*a58d3d2aSXin Li silk_assert( psIndices->GainsIndices[ 0 ] >= 0 && psIndices->GainsIndices[ 0 ] < MAX_DELTA_GAIN_QUANT - MIN_DELTA_GAIN_QUANT + 1 );
74*a58d3d2aSXin Li ec_enc_icdf( psRangeEnc, psIndices->GainsIndices[ 0 ], silk_delta_gain_iCDF, 8 );
75*a58d3d2aSXin Li } else {
76*a58d3d2aSXin Li /* independent coding, in two stages: MSB bits followed by 3 LSBs */
77*a58d3d2aSXin Li silk_assert( psIndices->GainsIndices[ 0 ] >= 0 && psIndices->GainsIndices[ 0 ] < N_LEVELS_QGAIN );
78*a58d3d2aSXin Li ec_enc_icdf( psRangeEnc, silk_RSHIFT( psIndices->GainsIndices[ 0 ], 3 ), silk_gain_iCDF[ psIndices->signalType ], 8 );
79*a58d3d2aSXin Li ec_enc_icdf( psRangeEnc, psIndices->GainsIndices[ 0 ] & 7, silk_uniform8_iCDF, 8 );
80*a58d3d2aSXin Li }
81*a58d3d2aSXin Li
82*a58d3d2aSXin Li /* remaining subframes */
83*a58d3d2aSXin Li for( i = 1; i < psEncC->nb_subfr; i++ ) {
84*a58d3d2aSXin Li silk_assert( psIndices->GainsIndices[ i ] >= 0 && psIndices->GainsIndices[ i ] < MAX_DELTA_GAIN_QUANT - MIN_DELTA_GAIN_QUANT + 1 );
85*a58d3d2aSXin Li ec_enc_icdf( psRangeEnc, psIndices->GainsIndices[ i ], silk_delta_gain_iCDF, 8 );
86*a58d3d2aSXin Li }
87*a58d3d2aSXin Li
88*a58d3d2aSXin Li /****************/
89*a58d3d2aSXin Li /* Encode NLSFs */
90*a58d3d2aSXin Li /****************/
91*a58d3d2aSXin Li ec_enc_icdf( psRangeEnc, psIndices->NLSFIndices[ 0 ], &psEncC->psNLSF_CB->CB1_iCDF[ ( psIndices->signalType >> 1 ) * psEncC->psNLSF_CB->nVectors ], 8 );
92*a58d3d2aSXin Li silk_NLSF_unpack( ec_ix, pred_Q8, psEncC->psNLSF_CB, psIndices->NLSFIndices[ 0 ] );
93*a58d3d2aSXin Li celt_assert( psEncC->psNLSF_CB->order == psEncC->predictLPCOrder );
94*a58d3d2aSXin Li for( i = 0; i < psEncC->psNLSF_CB->order; i++ ) {
95*a58d3d2aSXin Li if( psIndices->NLSFIndices[ i+1 ] >= NLSF_QUANT_MAX_AMPLITUDE ) {
96*a58d3d2aSXin Li ec_enc_icdf( psRangeEnc, 2 * NLSF_QUANT_MAX_AMPLITUDE, &psEncC->psNLSF_CB->ec_iCDF[ ec_ix[ i ] ], 8 );
97*a58d3d2aSXin Li ec_enc_icdf( psRangeEnc, psIndices->NLSFIndices[ i+1 ] - NLSF_QUANT_MAX_AMPLITUDE, silk_NLSF_EXT_iCDF, 8 );
98*a58d3d2aSXin Li } else if( psIndices->NLSFIndices[ i+1 ] <= -NLSF_QUANT_MAX_AMPLITUDE ) {
99*a58d3d2aSXin Li ec_enc_icdf( psRangeEnc, 0, &psEncC->psNLSF_CB->ec_iCDF[ ec_ix[ i ] ], 8 );
100*a58d3d2aSXin Li ec_enc_icdf( psRangeEnc, -psIndices->NLSFIndices[ i+1 ] - NLSF_QUANT_MAX_AMPLITUDE, silk_NLSF_EXT_iCDF, 8 );
101*a58d3d2aSXin Li } else {
102*a58d3d2aSXin Li ec_enc_icdf( psRangeEnc, psIndices->NLSFIndices[ i+1 ] + NLSF_QUANT_MAX_AMPLITUDE, &psEncC->psNLSF_CB->ec_iCDF[ ec_ix[ i ] ], 8 );
103*a58d3d2aSXin Li }
104*a58d3d2aSXin Li }
105*a58d3d2aSXin Li
106*a58d3d2aSXin Li /* Encode NLSF interpolation factor */
107*a58d3d2aSXin Li if( psEncC->nb_subfr == MAX_NB_SUBFR ) {
108*a58d3d2aSXin Li silk_assert( psIndices->NLSFInterpCoef_Q2 >= 0 && psIndices->NLSFInterpCoef_Q2 < 5 );
109*a58d3d2aSXin Li ec_enc_icdf( psRangeEnc, psIndices->NLSFInterpCoef_Q2, silk_NLSF_interpolation_factor_iCDF, 8 );
110*a58d3d2aSXin Li }
111*a58d3d2aSXin Li
112*a58d3d2aSXin Li if( psIndices->signalType == TYPE_VOICED )
113*a58d3d2aSXin Li {
114*a58d3d2aSXin Li /*********************/
115*a58d3d2aSXin Li /* Encode pitch lags */
116*a58d3d2aSXin Li /*********************/
117*a58d3d2aSXin Li /* lag index */
118*a58d3d2aSXin Li encode_absolute_lagIndex = 1;
119*a58d3d2aSXin Li if( condCoding == CODE_CONDITIONALLY && psEncC->ec_prevSignalType == TYPE_VOICED ) {
120*a58d3d2aSXin Li /* Delta Encoding */
121*a58d3d2aSXin Li delta_lagIndex = psIndices->lagIndex - psEncC->ec_prevLagIndex;
122*a58d3d2aSXin Li if( delta_lagIndex < -8 || delta_lagIndex > 11 ) {
123*a58d3d2aSXin Li delta_lagIndex = 0;
124*a58d3d2aSXin Li } else {
125*a58d3d2aSXin Li delta_lagIndex = delta_lagIndex + 9;
126*a58d3d2aSXin Li encode_absolute_lagIndex = 0; /* Only use delta */
127*a58d3d2aSXin Li }
128*a58d3d2aSXin Li silk_assert( delta_lagIndex >= 0 && delta_lagIndex < 21 );
129*a58d3d2aSXin Li ec_enc_icdf( psRangeEnc, delta_lagIndex, silk_pitch_delta_iCDF, 8 );
130*a58d3d2aSXin Li }
131*a58d3d2aSXin Li if( encode_absolute_lagIndex ) {
132*a58d3d2aSXin Li /* Absolute encoding */
133*a58d3d2aSXin Li opus_int32 pitch_high_bits, pitch_low_bits;
134*a58d3d2aSXin Li pitch_high_bits = silk_DIV32_16( psIndices->lagIndex, silk_RSHIFT( psEncC->fs_kHz, 1 ) );
135*a58d3d2aSXin Li pitch_low_bits = psIndices->lagIndex - silk_SMULBB( pitch_high_bits, silk_RSHIFT( psEncC->fs_kHz, 1 ) );
136*a58d3d2aSXin Li silk_assert( pitch_low_bits < psEncC->fs_kHz / 2 );
137*a58d3d2aSXin Li silk_assert( pitch_high_bits < 32 );
138*a58d3d2aSXin Li ec_enc_icdf( psRangeEnc, pitch_high_bits, silk_pitch_lag_iCDF, 8 );
139*a58d3d2aSXin Li ec_enc_icdf( psRangeEnc, pitch_low_bits, psEncC->pitch_lag_low_bits_iCDF, 8 );
140*a58d3d2aSXin Li }
141*a58d3d2aSXin Li psEncC->ec_prevLagIndex = psIndices->lagIndex;
142*a58d3d2aSXin Li
143*a58d3d2aSXin Li /* Countour index */
144*a58d3d2aSXin Li silk_assert( psIndices->contourIndex >= 0 );
145*a58d3d2aSXin Li silk_assert( ( psIndices->contourIndex < 34 && psEncC->fs_kHz > 8 && psEncC->nb_subfr == 4 ) ||
146*a58d3d2aSXin Li ( psIndices->contourIndex < 11 && psEncC->fs_kHz == 8 && psEncC->nb_subfr == 4 ) ||
147*a58d3d2aSXin Li ( psIndices->contourIndex < 12 && psEncC->fs_kHz > 8 && psEncC->nb_subfr == 2 ) ||
148*a58d3d2aSXin Li ( psIndices->contourIndex < 3 && psEncC->fs_kHz == 8 && psEncC->nb_subfr == 2 ) );
149*a58d3d2aSXin Li ec_enc_icdf( psRangeEnc, psIndices->contourIndex, psEncC->pitch_contour_iCDF, 8 );
150*a58d3d2aSXin Li
151*a58d3d2aSXin Li /********************/
152*a58d3d2aSXin Li /* Encode LTP gains */
153*a58d3d2aSXin Li /********************/
154*a58d3d2aSXin Li /* PERIndex value */
155*a58d3d2aSXin Li silk_assert( psIndices->PERIndex >= 0 && psIndices->PERIndex < 3 );
156*a58d3d2aSXin Li ec_enc_icdf( psRangeEnc, psIndices->PERIndex, silk_LTP_per_index_iCDF, 8 );
157*a58d3d2aSXin Li
158*a58d3d2aSXin Li /* Codebook Indices */
159*a58d3d2aSXin Li for( k = 0; k < psEncC->nb_subfr; k++ ) {
160*a58d3d2aSXin Li silk_assert( psIndices->LTPIndex[ k ] >= 0 && psIndices->LTPIndex[ k ] < ( 8 << psIndices->PERIndex ) );
161*a58d3d2aSXin Li ec_enc_icdf( psRangeEnc, psIndices->LTPIndex[ k ], silk_LTP_gain_iCDF_ptrs[ psIndices->PERIndex ], 8 );
162*a58d3d2aSXin Li }
163*a58d3d2aSXin Li
164*a58d3d2aSXin Li /**********************/
165*a58d3d2aSXin Li /* Encode LTP scaling */
166*a58d3d2aSXin Li /**********************/
167*a58d3d2aSXin Li if( condCoding == CODE_INDEPENDENTLY ) {
168*a58d3d2aSXin Li silk_assert( psIndices->LTP_scaleIndex >= 0 && psIndices->LTP_scaleIndex < 3 );
169*a58d3d2aSXin Li ec_enc_icdf( psRangeEnc, psIndices->LTP_scaleIndex, silk_LTPscale_iCDF, 8 );
170*a58d3d2aSXin Li }
171*a58d3d2aSXin Li silk_assert( !condCoding || psIndices->LTP_scaleIndex == 0 );
172*a58d3d2aSXin Li }
173*a58d3d2aSXin Li
174*a58d3d2aSXin Li psEncC->ec_prevSignalType = psIndices->signalType;
175*a58d3d2aSXin Li
176*a58d3d2aSXin Li /***************/
177*a58d3d2aSXin Li /* Encode seed */
178*a58d3d2aSXin Li /***************/
179*a58d3d2aSXin Li silk_assert( psIndices->Seed >= 0 && psIndices->Seed < 4 );
180*a58d3d2aSXin Li ec_enc_icdf( psRangeEnc, psIndices->Seed, silk_uniform4_iCDF, 8 );
181*a58d3d2aSXin Li }
182