Mercurial > hg > sv-dependency-builds

annotate src/opus-1.3/silk/fixed/encode_frame_FIX.c @ 79:91c729825bca pa_catalina

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Update build for AUDIO_COMPONENT_FIX
author Chris Cannam
date Wed, 30 Oct 2019 12:40:34 +0000
parents 7aeed7906520
children
rev   line source
Chris@69 1 /***********************************************************************
Chris@69 2 Copyright (c) 2006-2011, Skype Limited. All rights reserved.
Chris@69 3 Redistribution and use in source and binary forms, with or without
Chris@69 4 modification, are permitted provided that the following conditions
Chris@69 5 are met:
Chris@69 6 - Redistributions of source code must retain the above copyright notice,
Chris@69 7 this list of conditions and the following disclaimer.
Chris@69 8 - Redistributions in binary form must reproduce the above copyright
Chris@69 9 notice, this list of conditions and the following disclaimer in the
Chris@69 10 documentation and/or other materials provided with the distribution.
Chris@69 11 - Neither the name of Internet Society, IETF or IETF Trust, nor the
Chris@69 12 names of specific contributors, may be used to endorse or promote
Chris@69 13 products derived from this software without specific prior written
Chris@69 14 permission.
Chris@69 15 THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
Chris@69 16 AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
Chris@69 17 IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
Chris@69 18 ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
Chris@69 19 LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
Chris@69 20 CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
Chris@69 21 SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
Chris@69 22 INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
Chris@69 23 CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
Chris@69 24 ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
Chris@69 25 POSSIBILITY OF SUCH DAMAGE.
Chris@69 26 ***********************************************************************/
Chris@69 27
Chris@69 28 #ifdef HAVE_CONFIG_H
Chris@69 29 #include "config.h"
Chris@69 30 #endif
Chris@69 31
Chris@69 32 #include <stdlib.h>
Chris@69 33 #include "main_FIX.h"
Chris@69 34 #include "stack_alloc.h"
Chris@69 35 #include "tuning_parameters.h"
Chris@69 36
Chris@69 37 /* Low Bitrate Redundancy (LBRR) encoding. Reuse all parameters but encode with lower bitrate */
Chris@69 38 static OPUS_INLINE void silk_LBRR_encode_FIX(
Chris@69 39 silk_encoder_state_FIX *psEnc, /* I/O Pointer to Silk FIX encoder state */
Chris@69 40 silk_encoder_control_FIX *psEncCtrl, /* I/O Pointer to Silk FIX encoder control struct */
Chris@69 41 const opus_int16 x16[], /* I Input signal */
Chris@69 42 opus_int condCoding /* I The type of conditional coding used so far for this frame */
Chris@69 43 );
Chris@69 44
Chris@69 45 void silk_encode_do_VAD_FIX(
Chris@69 46 silk_encoder_state_FIX *psEnc, /* I/O Pointer to Silk FIX encoder state */
Chris@69 47 opus_int activity /* I Decision of Opus voice activity detector */
Chris@69 48 )
Chris@69 49 {
Chris@69 50 const opus_int activity_threshold = SILK_FIX_CONST( SPEECH_ACTIVITY_DTX_THRES, 8 );
Chris@69 51
Chris@69 52 /****************************/
Chris@69 53 /* Voice Activity Detection */
Chris@69 54 /****************************/
Chris@69 55 silk_VAD_GetSA_Q8( &psEnc->sCmn, psEnc->sCmn.inputBuf + 1, psEnc->sCmn.arch );
Chris@69 56 /* If Opus VAD is inactive and Silk VAD is active: lower Silk VAD to just under the threshold */
Chris@69 57 if( activity == VAD_NO_ACTIVITY && psEnc->sCmn.speech_activity_Q8 >= activity_threshold ) {
Chris@69 58 psEnc->sCmn.speech_activity_Q8 = activity_threshold - 1;
Chris@69 59 }
Chris@69 60
Chris@69 61 /**************************************************/
Chris@69 62 /* Convert speech activity into VAD and DTX flags */
Chris@69 63 /**************************************************/
Chris@69 64 if( psEnc->sCmn.speech_activity_Q8 < activity_threshold ) {
Chris@69 65 psEnc->sCmn.indices.signalType = TYPE_NO_VOICE_ACTIVITY;
Chris@69 66 psEnc->sCmn.noSpeechCounter++;
Chris@69 67 if( psEnc->sCmn.noSpeechCounter <= NB_SPEECH_FRAMES_BEFORE_DTX ) {
Chris@69 68 psEnc->sCmn.inDTX = 0;
Chris@69 69 } else if( psEnc->sCmn.noSpeechCounter > MAX_CONSECUTIVE_DTX + NB_SPEECH_FRAMES_BEFORE_DTX ) {
Chris@69 70 psEnc->sCmn.noSpeechCounter = NB_SPEECH_FRAMES_BEFORE_DTX;
Chris@69 71 psEnc->sCmn.inDTX = 0;
Chris@69 72 }
Chris@69 73 psEnc->sCmn.VAD_flags[ psEnc->sCmn.nFramesEncoded ] = 0;
Chris@69 74 } else {
Chris@69 75 psEnc->sCmn.noSpeechCounter = 0;
Chris@69 76 psEnc->sCmn.inDTX = 0;
Chris@69 77 psEnc->sCmn.indices.signalType = TYPE_UNVOICED;
Chris@69 78 psEnc->sCmn.VAD_flags[ psEnc->sCmn.nFramesEncoded ] = 1;
Chris@69 79 }
Chris@69 80 }
Chris@69 81
Chris@69 82 /****************/
Chris@69 83 /* Encode frame */
Chris@69 84 /****************/
Chris@69 85 opus_int silk_encode_frame_FIX(
Chris@69 86 silk_encoder_state_FIX *psEnc, /* I/O Pointer to Silk FIX encoder state */
Chris@69 87 opus_int32 *pnBytesOut, /* O Pointer to number of payload bytes; */
Chris@69 88 ec_enc *psRangeEnc, /* I/O compressor data structure */
Chris@69 89 opus_int condCoding, /* I The type of conditional coding to use */
Chris@69 90 opus_int maxBits, /* I If > 0: maximum number of output bits */
Chris@69 91 opus_int useCBR /* I Flag to force constant-bitrate operation */
Chris@69 92 )
Chris@69 93 {
Chris@69 94 silk_encoder_control_FIX sEncCtrl;
Chris@69 95 opus_int i, iter, maxIter, found_upper, found_lower, ret = 0;
Chris@69 96 opus_int16 *x_frame;
Chris@69 97 ec_enc sRangeEnc_copy, sRangeEnc_copy2;
Chris@69 98 silk_nsq_state sNSQ_copy, sNSQ_copy2;
Chris@69 99 opus_int32 seed_copy, nBits, nBits_lower, nBits_upper, gainMult_lower, gainMult_upper;
Chris@69 100 opus_int32 gainsID, gainsID_lower, gainsID_upper;
Chris@69 101 opus_int16 gainMult_Q8;
Chris@69 102 opus_int16 ec_prevLagIndex_copy;
Chris@69 103 opus_int ec_prevSignalType_copy;
Chris@69 104 opus_int8 LastGainIndex_copy2;
Chris@69 105 opus_int gain_lock[ MAX_NB_SUBFR ] = {0};
Chris@69 106 opus_int16 best_gain_mult[ MAX_NB_SUBFR ];
Chris@69 107 opus_int best_sum[ MAX_NB_SUBFR ];
Chris@69 108 SAVE_STACK;
Chris@69 109
Chris@69 110 /* This is totally unnecessary but many compilers (including gcc) are too dumb to realise it */
Chris@69 111 LastGainIndex_copy2 = nBits_lower = nBits_upper = gainMult_lower = gainMult_upper = 0;
Chris@69 112
Chris@69 113 psEnc->sCmn.indices.Seed = psEnc->sCmn.frameCounter++ & 3;
Chris@69 114
Chris@69 115 /**************************************************************/
Chris@69 116 /* Set up Input Pointers, and insert frame in input buffer */
Chris@69 117 /*************************************************************/
Chris@69 118 /* start of frame to encode */
Chris@69 119 x_frame = psEnc->x_buf + psEnc->sCmn.ltp_mem_length;
Chris@69 120
Chris@69 121 /***************************************/
Chris@69 122 /* Ensure smooth bandwidth transitions */
Chris@69 123 /***************************************/
Chris@69 124 silk_LP_variable_cutoff( &psEnc->sCmn.sLP, psEnc->sCmn.inputBuf + 1, psEnc->sCmn.frame_length );
Chris@69 125
Chris@69 126 /*******************************************/
Chris@69 127 /* Copy new frame to front of input buffer */
Chris@69 128 /*******************************************/
Chris@69 129 silk_memcpy( x_frame + LA_SHAPE_MS * psEnc->sCmn.fs_kHz, psEnc->sCmn.inputBuf + 1, psEnc->sCmn.frame_length * sizeof( opus_int16 ) );
Chris@69 130
Chris@69 131 if( !psEnc->sCmn.prefillFlag ) {
Chris@69 132 VARDECL( opus_int16, res_pitch );
Chris@69 133 VARDECL( opus_uint8, ec_buf_copy );
Chris@69 134 opus_int16 *res_pitch_frame;
Chris@69 135
Chris@69 136 ALLOC( res_pitch,
Chris@69 137 psEnc->sCmn.la_pitch + psEnc->sCmn.frame_length
Chris@69 138 + psEnc->sCmn.ltp_mem_length, opus_int16 );
Chris@69 139 /* start of pitch LPC residual frame */
Chris@69 140 res_pitch_frame = res_pitch + psEnc->sCmn.ltp_mem_length;
Chris@69 141
Chris@69 142 /*****************************************/
Chris@69 143 /* Find pitch lags, initial LPC analysis */
Chris@69 144 /*****************************************/
Chris@69 145 silk_find_pitch_lags_FIX( psEnc, &sEncCtrl, res_pitch, x_frame - psEnc->sCmn.ltp_mem_length, psEnc->sCmn.arch );
Chris@69 146
Chris@69 147 /************************/
Chris@69 148 /* Noise shape analysis */
Chris@69 149 /************************/
Chris@69 150 silk_noise_shape_analysis_FIX( psEnc, &sEncCtrl, res_pitch_frame, x_frame, psEnc->sCmn.arch );
Chris@69 151
Chris@69 152 /***************************************************/
Chris@69 153 /* Find linear prediction coefficients (LPC + LTP) */
Chris@69 154 /***************************************************/
Chris@69 155 silk_find_pred_coefs_FIX( psEnc, &sEncCtrl, res_pitch_frame, x_frame, condCoding );
Chris@69 156
Chris@69 157 /****************************************/
Chris@69 158 /* Process gains */
Chris@69 159 /****************************************/
Chris@69 160 silk_process_gains_FIX( psEnc, &sEncCtrl, condCoding );
Chris@69 161
Chris@69 162 /****************************************/
Chris@69 163 /* Low Bitrate Redundant Encoding */
Chris@69 164 /****************************************/
Chris@69 165 silk_LBRR_encode_FIX( psEnc, &sEncCtrl, x_frame, condCoding );
Chris@69 166
Chris@69 167 /* Loop over quantizer and entropy coding to control bitrate */
Chris@69 168 maxIter = 6;
Chris@69 169 gainMult_Q8 = SILK_FIX_CONST( 1, 8 );
Chris@69 170 found_lower = 0;
Chris@69 171 found_upper = 0;
Chris@69 172 gainsID = silk_gains_ID( psEnc->sCmn.indices.GainsIndices, psEnc->sCmn.nb_subfr );
Chris@69 173 gainsID_lower = -1;
Chris@69 174 gainsID_upper = -1;
Chris@69 175 /* Copy part of the input state */
Chris@69 176 silk_memcpy( &sRangeEnc_copy, psRangeEnc, sizeof( ec_enc ) );
Chris@69 177 silk_memcpy( &sNSQ_copy, &psEnc->sCmn.sNSQ, sizeof( silk_nsq_state ) );
Chris@69 178 seed_copy = psEnc->sCmn.indices.Seed;
Chris@69 179 ec_prevLagIndex_copy = psEnc->sCmn.ec_prevLagIndex;
Chris@69 180 ec_prevSignalType_copy = psEnc->sCmn.ec_prevSignalType;
Chris@69 181 ALLOC( ec_buf_copy, 1275, opus_uint8 );
Chris@69 182 for( iter = 0; ; iter++ ) {
Chris@69 183 if( gainsID == gainsID_lower ) {
Chris@69 184 nBits = nBits_lower;
Chris@69 185 } else if( gainsID == gainsID_upper ) {
Chris@69 186 nBits = nBits_upper;
Chris@69 187 } else {
Chris@69 188 /* Restore part of the input state */
Chris@69 189 if( iter > 0 ) {
Chris@69 190 silk_memcpy( psRangeEnc, &sRangeEnc_copy, sizeof( ec_enc ) );
Chris@69 191 silk_memcpy( &psEnc->sCmn.sNSQ, &sNSQ_copy, sizeof( silk_nsq_state ) );
Chris@69 192 psEnc->sCmn.indices.Seed = seed_copy;
Chris@69 193 psEnc->sCmn.ec_prevLagIndex = ec_prevLagIndex_copy;
Chris@69 194 psEnc->sCmn.ec_prevSignalType = ec_prevSignalType_copy;
Chris@69 195 }
Chris@69 196
Chris@69 197 /*****************************************/
Chris@69 198 /* Noise shaping quantization */
Chris@69 199 /*****************************************/
Chris@69 200 if( psEnc->sCmn.nStatesDelayedDecision > 1 || psEnc->sCmn.warping_Q16 > 0 ) {
Chris@69 201 silk_NSQ_del_dec( &psEnc->sCmn, &psEnc->sCmn.sNSQ, &psEnc->sCmn.indices, x_frame, psEnc->sCmn.pulses,
Chris@69 202 sEncCtrl.PredCoef_Q12[ 0 ], sEncCtrl.LTPCoef_Q14, sEncCtrl.AR_Q13, sEncCtrl.HarmShapeGain_Q14,
Chris@69 203 sEncCtrl.Tilt_Q14, sEncCtrl.LF_shp_Q14, sEncCtrl.Gains_Q16, sEncCtrl.pitchL, sEncCtrl.Lambda_Q10, sEncCtrl.LTP_scale_Q14,
Chris@69 204 psEnc->sCmn.arch );
Chris@69 205 } else {
Chris@69 206 silk_NSQ( &psEnc->sCmn, &psEnc->sCmn.sNSQ, &psEnc->sCmn.indices, x_frame, psEnc->sCmn.pulses,
Chris@69 207 sEncCtrl.PredCoef_Q12[ 0 ], sEncCtrl.LTPCoef_Q14, sEncCtrl.AR_Q13, sEncCtrl.HarmShapeGain_Q14,
Chris@69 208 sEncCtrl.Tilt_Q14, sEncCtrl.LF_shp_Q14, sEncCtrl.Gains_Q16, sEncCtrl.pitchL, sEncCtrl.Lambda_Q10, sEncCtrl.LTP_scale_Q14,
Chris@69 209 psEnc->sCmn.arch);
Chris@69 210 }
Chris@69 211
Chris@69 212 if ( iter == maxIter && !found_lower ) {
Chris@69 213 silk_memcpy( &sRangeEnc_copy2, psRangeEnc, sizeof( ec_enc ) );
Chris@69 214 }
Chris@69 215
Chris@69 216 /****************************************/
Chris@69 217 /* Encode Parameters */
Chris@69 218 /****************************************/
Chris@69 219 silk_encode_indices( &psEnc->sCmn, psRangeEnc, psEnc->sCmn.nFramesEncoded, 0, condCoding );
Chris@69 220
Chris@69 221 /****************************************/
Chris@69 222 /* Encode Excitation Signal */
Chris@69 223 /****************************************/
Chris@69 224 silk_encode_pulses( psRangeEnc, psEnc->sCmn.indices.signalType, psEnc->sCmn.indices.quantOffsetType,
Chris@69 225 psEnc->sCmn.pulses, psEnc->sCmn.frame_length );
Chris@69 226
Chris@69 227 nBits = ec_tell( psRangeEnc );
Chris@69 228
Chris@69 229 /* If we still bust after the last iteration, do some damage control. */
Chris@69 230 if ( iter == maxIter && !found_lower && nBits > maxBits ) {
Chris@69 231 silk_memcpy( psRangeEnc, &sRangeEnc_copy2, sizeof( ec_enc ) );
Chris@69 232
Chris@69 233 /* Keep gains the same as the last frame. */
Chris@69 234 psEnc->sShape.LastGainIndex = sEncCtrl.lastGainIndexPrev;
Chris@69 235 for ( i = 0; i < psEnc->sCmn.nb_subfr; i++ ) {
Chris@69 236 psEnc->sCmn.indices.GainsIndices[ i ] = 4;
Chris@69 237 }
Chris@69 238 if (condCoding != CODE_CONDITIONALLY) {
Chris@69 239 psEnc->sCmn.indices.GainsIndices[ 0 ] = sEncCtrl.lastGainIndexPrev;
Chris@69 240 }
Chris@69 241 psEnc->sCmn.ec_prevLagIndex = ec_prevLagIndex_copy;
Chris@69 242 psEnc->sCmn.ec_prevSignalType = ec_prevSignalType_copy;
Chris@69 243 /* Clear all pulses. */
Chris@69 244 for ( i = 0; i < psEnc->sCmn.frame_length; i++ ) {
Chris@69 245 psEnc->sCmn.pulses[ i ] = 0;
Chris@69 246 }
Chris@69 247
Chris@69 248 silk_encode_indices( &psEnc->sCmn, psRangeEnc, psEnc->sCmn.nFramesEncoded, 0, condCoding );
Chris@69 249
Chris@69 250 silk_encode_pulses( psRangeEnc, psEnc->sCmn.indices.signalType, psEnc->sCmn.indices.quantOffsetType,
Chris@69 251 psEnc->sCmn.pulses, psEnc->sCmn.frame_length );
Chris@69 252
Chris@69 253 nBits = ec_tell( psRangeEnc );
Chris@69 254 }
Chris@69 255
Chris@69 256 if( useCBR == 0 && iter == 0 && nBits <= maxBits ) {
Chris@69 257 break;
Chris@69 258 }
Chris@69 259 }
Chris@69 260
Chris@69 261 if( iter == maxIter ) {
Chris@69 262 if( found_lower && ( gainsID == gainsID_lower || nBits > maxBits ) ) {
Chris@69 263 /* Restore output state from earlier iteration that did meet the bitrate budget */
Chris@69 264 silk_memcpy( psRangeEnc, &sRangeEnc_copy2, sizeof( ec_enc ) );
Chris@69 265 celt_assert( sRangeEnc_copy2.offs <= 1275 );
Chris@69 266 silk_memcpy( psRangeEnc->buf, ec_buf_copy, sRangeEnc_copy2.offs );
Chris@69 267 silk_memcpy( &psEnc->sCmn.sNSQ, &sNSQ_copy2, sizeof( silk_nsq_state ) );
Chris@69 268 psEnc->sShape.LastGainIndex = LastGainIndex_copy2;
Chris@69 269 }
Chris@69 270 break;
Chris@69 271 }
Chris@69 272
Chris@69 273 if( nBits > maxBits ) {
Chris@69 274 if( found_lower == 0 && iter >= 2 ) {
Chris@69 275 /* Adjust the quantizer's rate/distortion tradeoff and discard previous "upper" results */
Chris@69 276 sEncCtrl.Lambda_Q10 = silk_ADD_RSHIFT32( sEncCtrl.Lambda_Q10, sEncCtrl.Lambda_Q10, 1 );
Chris@69 277 found_upper = 0;
Chris@69 278 gainsID_upper = -1;
Chris@69 279 } else {
Chris@69 280 found_upper = 1;
Chris@69 281 nBits_upper = nBits;
Chris@69 282 gainMult_upper = gainMult_Q8;
Chris@69 283 gainsID_upper = gainsID;
Chris@69 284 }
Chris@69 285 } else if( nBits < maxBits - 5 ) {
Chris@69 286 found_lower = 1;
Chris@69 287 nBits_lower = nBits;
Chris@69 288 gainMult_lower = gainMult_Q8;
Chris@69 289 if( gainsID != gainsID_lower ) {
Chris@69 290 gainsID_lower = gainsID;
Chris@69 291 /* Copy part of the output state */
Chris@69 292 silk_memcpy( &sRangeEnc_copy2, psRangeEnc, sizeof( ec_enc ) );
Chris@69 293 celt_assert( psRangeEnc->offs <= 1275 );
Chris@69 294 silk_memcpy( ec_buf_copy, psRangeEnc->buf, psRangeEnc->offs );
Chris@69 295 silk_memcpy( &sNSQ_copy2, &psEnc->sCmn.sNSQ, sizeof( silk_nsq_state ) );
Chris@69 296 LastGainIndex_copy2 = psEnc->sShape.LastGainIndex;
Chris@69 297 }
Chris@69 298 } else {
Chris@69 299 /* Within 5 bits of budget: close enough */
Chris@69 300 break;
Chris@69 301 }
Chris@69 302
Chris@69 303 if ( !found_lower && nBits > maxBits ) {
Chris@69 304 int j;
Chris@69 305 for ( i = 0; i < psEnc->sCmn.nb_subfr; i++ ) {
Chris@69 306 int sum=0;
Chris@69 307 for ( j = i*psEnc->sCmn.subfr_length; j < (i+1)*psEnc->sCmn.subfr_length; j++ ) {
Chris@69 308 sum += abs( psEnc->sCmn.pulses[j] );
Chris@69 309 }
Chris@69 310 if ( iter == 0 || (sum < best_sum[i] && !gain_lock[i]) ) {
Chris@69 311 best_sum[i] = sum;
Chris@69 312 best_gain_mult[i] = gainMult_Q8;
Chris@69 313 } else {
Chris@69 314 gain_lock[i] = 1;
Chris@69 315 }
Chris@69 316 }
Chris@69 317 }
Chris@69 318 if( ( found_lower & found_upper ) == 0 ) {
Chris@69 319 /* Adjust gain according to high-rate rate/distortion curve */
Chris@69 320 if( nBits > maxBits ) {
Chris@69 321 if (gainMult_Q8 < 16384) {
Chris@69 322 gainMult_Q8 *= 2;
Chris@69 323 } else {
Chris@69 324 gainMult_Q8 = 32767;
Chris@69 325 }
Chris@69 326 } else {
Chris@69 327 opus_int32 gain_factor_Q16;
Chris@69 328 gain_factor_Q16 = silk_log2lin( silk_LSHIFT( nBits - maxBits, 7 ) / psEnc->sCmn.frame_length + SILK_FIX_CONST( 16, 7 ) );
Chris@69 329 gainMult_Q8 = silk_SMULWB( gain_factor_Q16, gainMult_Q8 );
Chris@69 330 }
Chris@69 331
Chris@69 332 } else {
Chris@69 333 /* Adjust gain by interpolating */
Chris@69 334 gainMult_Q8 = gainMult_lower + silk_DIV32_16( silk_MUL( gainMult_upper - gainMult_lower, maxBits - nBits_lower ), nBits_upper - nBits_lower );
Chris@69 335 /* New gain multplier must be between 25% and 75% of old range (note that gainMult_upper < gainMult_lower) */
Chris@69 336 if( gainMult_Q8 > silk_ADD_RSHIFT32( gainMult_lower, gainMult_upper - gainMult_lower, 2 ) ) {
Chris@69 337 gainMult_Q8 = silk_ADD_RSHIFT32( gainMult_lower, gainMult_upper - gainMult_lower, 2 );
Chris@69 338 } else
Chris@69 339 if( gainMult_Q8 < silk_SUB_RSHIFT32( gainMult_upper, gainMult_upper - gainMult_lower, 2 ) ) {
Chris@69 340 gainMult_Q8 = silk_SUB_RSHIFT32( gainMult_upper, gainMult_upper - gainMult_lower, 2 );
Chris@69 341 }
Chris@69 342 }
Chris@69 343
Chris@69 344 for( i = 0; i < psEnc->sCmn.nb_subfr; i++ ) {
Chris@69 345 opus_int16 tmp;
Chris@69 346 if ( gain_lock[i] ) {
Chris@69 347 tmp = best_gain_mult[i];
Chris@69 348 } else {
Chris@69 349 tmp = gainMult_Q8;
Chris@69 350 }
Chris@69 351 sEncCtrl.Gains_Q16[ i ] = silk_LSHIFT_SAT32( silk_SMULWB( sEncCtrl.GainsUnq_Q16[ i ], tmp ), 8 );
Chris@69 352 }
Chris@69 353
Chris@69 354 /* Quantize gains */
Chris@69 355 psEnc->sShape.LastGainIndex = sEncCtrl.lastGainIndexPrev;
Chris@69 356 silk_gains_quant( psEnc->sCmn.indices.GainsIndices, sEncCtrl.Gains_Q16,
Chris@69 357 &psEnc->sShape.LastGainIndex, condCoding == CODE_CONDITIONALLY, psEnc->sCmn.nb_subfr );
Chris@69 358
Chris@69 359 /* Unique identifier of gains vector */
Chris@69 360 gainsID = silk_gains_ID( psEnc->sCmn.indices.GainsIndices, psEnc->sCmn.nb_subfr );
Chris@69 361 }
Chris@69 362 }
Chris@69 363
Chris@69 364 /* Update input buffer */
Chris@69 365 silk_memmove( psEnc->x_buf, &psEnc->x_buf[ psEnc->sCmn.frame_length ],
Chris@69 366 ( psEnc->sCmn.ltp_mem_length + LA_SHAPE_MS * psEnc->sCmn.fs_kHz ) * sizeof( opus_int16 ) );
Chris@69 367
Chris@69 368 /* Exit without entropy coding */
Chris@69 369 if( psEnc->sCmn.prefillFlag ) {
Chris@69 370 /* No payload */
Chris@69 371 *pnBytesOut = 0;
Chris@69 372 RESTORE_STACK;
Chris@69 373 return ret;
Chris@69 374 }
Chris@69 375
Chris@69 376 /* Parameters needed for next frame */
Chris@69 377 psEnc->sCmn.prevLag = sEncCtrl.pitchL[ psEnc->sCmn.nb_subfr - 1 ];
Chris@69 378 psEnc->sCmn.prevSignalType = psEnc->sCmn.indices.signalType;
Chris@69 379
Chris@69 380 /****************************************/
Chris@69 381 /* Finalize payload */
Chris@69 382 /****************************************/
Chris@69 383 psEnc->sCmn.first_frame_after_reset = 0;
Chris@69 384 /* Payload size */
Chris@69 385 *pnBytesOut = silk_RSHIFT( ec_tell( psRangeEnc ) + 7, 3 );
Chris@69 386
Chris@69 387 RESTORE_STACK;
Chris@69 388 return ret;
Chris@69 389 }
Chris@69 390
Chris@69 391 /* Low-Bitrate Redundancy (LBRR) encoding. Reuse all parameters but encode excitation at lower bitrate */
Chris@69 392 static OPUS_INLINE void silk_LBRR_encode_FIX(
Chris@69 393 silk_encoder_state_FIX *psEnc, /* I/O Pointer to Silk FIX encoder state */
Chris@69 394 silk_encoder_control_FIX *psEncCtrl, /* I/O Pointer to Silk FIX encoder control struct */
Chris@69 395 const opus_int16 x16[], /* I Input signal */
Chris@69 396 opus_int condCoding /* I The type of conditional coding used so far for this frame */
Chris@69 397 )
Chris@69 398 {
Chris@69 399 opus_int32 TempGains_Q16[ MAX_NB_SUBFR ];
Chris@69 400 SideInfoIndices *psIndices_LBRR = &psEnc->sCmn.indices_LBRR[ psEnc->sCmn.nFramesEncoded ];
Chris@69 401 silk_nsq_state sNSQ_LBRR;
Chris@69 402
Chris@69 403 /*******************************************/
Chris@69 404 /* Control use of inband LBRR */
Chris@69 405 /*******************************************/
Chris@69 406 if( psEnc->sCmn.LBRR_enabled && psEnc->sCmn.speech_activity_Q8 > SILK_FIX_CONST( LBRR_SPEECH_ACTIVITY_THRES, 8 ) ) {
Chris@69 407 psEnc->sCmn.LBRR_flags[ psEnc->sCmn.nFramesEncoded ] = 1;
Chris@69 408
Chris@69 409 /* Copy noise shaping quantizer state and quantization indices from regular encoding */
Chris@69 410 silk_memcpy( &sNSQ_LBRR, &psEnc->sCmn.sNSQ, sizeof( silk_nsq_state ) );
Chris@69 411 silk_memcpy( psIndices_LBRR, &psEnc->sCmn.indices, sizeof( SideInfoIndices ) );
Chris@69 412
Chris@69 413 /* Save original gains */
Chris@69 414 silk_memcpy( TempGains_Q16, psEncCtrl->Gains_Q16, psEnc->sCmn.nb_subfr * sizeof( opus_int32 ) );
Chris@69 415
Chris@69 416 if( psEnc->sCmn.nFramesEncoded == 0 || psEnc->sCmn.LBRR_flags[ psEnc->sCmn.nFramesEncoded - 1 ] == 0 ) {
Chris@69 417 /* First frame in packet or previous frame not LBRR coded */
Chris@69 418 psEnc->sCmn.LBRRprevLastGainIndex = psEnc->sShape.LastGainIndex;
Chris@69 419
Chris@69 420 /* Increase Gains to get target LBRR rate */
Chris@69 421 psIndices_LBRR->GainsIndices[ 0 ] = psIndices_LBRR->GainsIndices[ 0 ] + psEnc->sCmn.LBRR_GainIncreases;
Chris@69 422 psIndices_LBRR->GainsIndices[ 0 ] = silk_min_int( psIndices_LBRR->GainsIndices[ 0 ], N_LEVELS_QGAIN - 1 );
Chris@69 423 }
Chris@69 424
Chris@69 425 /* Decode to get gains in sync with decoder */
Chris@69 426 /* Overwrite unquantized gains with quantized gains */
Chris@69 427 silk_gains_dequant( psEncCtrl->Gains_Q16, psIndices_LBRR->GainsIndices,
Chris@69 428 &psEnc->sCmn.LBRRprevLastGainIndex, condCoding == CODE_CONDITIONALLY, psEnc->sCmn.nb_subfr );
Chris@69 429
Chris@69 430 /*****************************************/
Chris@69 431 /* Noise shaping quantization */
Chris@69 432 /*****************************************/
Chris@69 433 if( psEnc->sCmn.nStatesDelayedDecision > 1 || psEnc->sCmn.warping_Q16 > 0 ) {
Chris@69 434 silk_NSQ_del_dec( &psEnc->sCmn, &sNSQ_LBRR, psIndices_LBRR, x16,
Chris@69 435 psEnc->sCmn.pulses_LBRR[ psEnc->sCmn.nFramesEncoded ], psEncCtrl->PredCoef_Q12[ 0 ], psEncCtrl->LTPCoef_Q14,
Chris@69 436 psEncCtrl->AR_Q13, psEncCtrl->HarmShapeGain_Q14, psEncCtrl->Tilt_Q14, psEncCtrl->LF_shp_Q14,
Chris@69 437 psEncCtrl->Gains_Q16, psEncCtrl->pitchL, psEncCtrl->Lambda_Q10, psEncCtrl->LTP_scale_Q14, psEnc->sCmn.arch );
Chris@69 438 } else {
Chris@69 439 silk_NSQ( &psEnc->sCmn, &sNSQ_LBRR, psIndices_LBRR, x16,
Chris@69 440 psEnc->sCmn.pulses_LBRR[ psEnc->sCmn.nFramesEncoded ], psEncCtrl->PredCoef_Q12[ 0 ], psEncCtrl->LTPCoef_Q14,
Chris@69 441 psEncCtrl->AR_Q13, psEncCtrl->HarmShapeGain_Q14, psEncCtrl->Tilt_Q14, psEncCtrl->LF_shp_Q14,
Chris@69 442 psEncCtrl->Gains_Q16, psEncCtrl->pitchL, psEncCtrl->Lambda_Q10, psEncCtrl->LTP_scale_Q14, psEnc->sCmn.arch );
Chris@69 443 }
Chris@69 444
Chris@69 445 /* Restore original gains */
Chris@69 446 silk_memcpy( psEncCtrl->Gains_Q16, TempGains_Q16, psEnc->sCmn.nb_subfr * sizeof( opus_int32 ) );
Chris@69 447 }
Chris@69 448 }