1885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org/*********************************************************************** 2885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.orgCopyright (c) 2006-2011, Skype Limited. All rights reserved. 3885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.orgRedistribution and use in source and binary forms, with or without 4885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.orgmodification, are permitted provided that the following conditions 5885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.orgare met: 6885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org- Redistributions of source code must retain the above copyright notice, 7885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.orgthis list of conditions and the following disclaimer. 8885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org- Redistributions in binary form must reproduce the above copyright 9885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.orgnotice, this list of conditions and the following disclaimer in the 10885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.orgdocumentation and/or other materials provided with the distribution. 11e3ea049fcaee2247e45f0ce793d4313babb4ef69tlegrand@chromium.org- Neither the name of Internet Society, IETF or IETF Trust, nor the 12885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.orgnames of specific contributors, may be used to endorse or promote 13885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.orgproducts derived from this software without specific prior written 14885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.orgpermission. 15e3ea049fcaee2247e45f0ce793d4313babb4ef69tlegrand@chromium.orgTHIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" 16885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.orgAND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE 17885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.orgIMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE 18885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.orgARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE 19885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.orgLIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR 20885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.orgCONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF 21885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.orgSUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS 22885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.orgINTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN 23885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.orgCONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) 24885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.orgARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE 25885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.orgPOSSIBILITY OF SUCH DAMAGE. 26885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org***********************************************************************/ 27885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 28885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org#ifdef HAVE_CONFIG_H 29885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org#include "config.h" 30885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org#endif 31885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 32885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org#include "main.h" 33e3ea049fcaee2247e45f0ce793d4313babb4ef69tlegrand@chromium.org#include "stack_alloc.h" 34885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 35885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org/* Silk VAD noise level estimation */ 363c3902f0ac13428394f14f78f0fab05ef3468d69tlegrand@google.comstatic OPUS_INLINE void silk_VAD_GetNoiseLevels( 37885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org const opus_int32 pX[ VAD_N_BANDS ], /* I subband energies */ 38885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org silk_VAD_state *psSilk_VAD /* I/O Pointer to Silk VAD state */ 39885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org); 40885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 41885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org/**********************************/ 42885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org/* Initialization of the Silk VAD */ 43885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org/**********************************/ 44885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.orgopus_int silk_VAD_Init( /* O Return value, 0 if success */ 45885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org silk_VAD_state *psSilk_VAD /* I/O Pointer to Silk VAD state */ 46885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org) 47885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org{ 48885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org opus_int b, ret = 0; 49885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 50885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* reset state memory */ 51885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org silk_memset( psSilk_VAD, 0, sizeof( silk_VAD_state ) ); 52885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 53885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* init noise levels */ 54885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* Initialize array with approx pink noise levels (psd proportional to inverse of frequency) */ 55885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org for( b = 0; b < VAD_N_BANDS; b++ ) { 56885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org psSilk_VAD->NoiseLevelBias[ b ] = silk_max_32( silk_DIV32_16( VAD_NOISE_LEVELS_BIAS, b + 1 ), 1 ); 57885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org } 58885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 59885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* Initialize state */ 60885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org for( b = 0; b < VAD_N_BANDS; b++ ) { 61885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org psSilk_VAD->NL[ b ] = silk_MUL( 100, psSilk_VAD->NoiseLevelBias[ b ] ); 62885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org psSilk_VAD->inv_NL[ b ] = silk_DIV32( silk_int32_MAX, psSilk_VAD->NL[ b ] ); 63885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org } 64885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org psSilk_VAD->counter = 15; 65885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 66885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* init smoothed energy-to-noise ratio*/ 67885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org for( b = 0; b < VAD_N_BANDS; b++ ) { 68885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org psSilk_VAD->NrgRatioSmth_Q8[ b ] = 100 * 256; /* 100 * 256 --> 20 dB SNR */ 69885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org } 70885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 71885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org return( ret ); 72885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org} 73885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 74885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org/* Weighting factors for tilt measure */ 75885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.orgstatic const opus_int32 tiltWeights[ VAD_N_BANDS ] = { 30000, 6000, -12000, -12000 }; 76885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 77885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org/***************************************/ 78885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org/* Get the speech activity level in Q8 */ 79885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org/***************************************/ 80885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.orgopus_int silk_VAD_GetSA_Q8( /* O Return value, 0 if success */ 81885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org silk_encoder_state *psEncC, /* I/O Encoder state */ 82885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org const opus_int16 pIn[] /* I PCM input */ 83885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org) 84885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org{ 85885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org opus_int SA_Q15, pSNR_dB_Q7, input_tilt; 86e3ea049fcaee2247e45f0ce793d4313babb4ef69tlegrand@chromium.org opus_int decimated_framelength1, decimated_framelength2; 87e3ea049fcaee2247e45f0ce793d4313babb4ef69tlegrand@chromium.org opus_int decimated_framelength; 88e3ea049fcaee2247e45f0ce793d4313babb4ef69tlegrand@chromium.org opus_int dec_subframe_length, dec_subframe_offset, SNR_Q7, i, b, s; 89885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org opus_int32 sumSquared, smooth_coef_Q16; 90885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org opus_int16 HPstateTmp; 91e3ea049fcaee2247e45f0ce793d4313babb4ef69tlegrand@chromium.org VARDECL( opus_int16, X ); 92885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org opus_int32 Xnrg[ VAD_N_BANDS ]; 93885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org opus_int32 NrgToNoiseRatio_Q8[ VAD_N_BANDS ]; 94885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org opus_int32 speech_nrg, x_tmp; 95e3ea049fcaee2247e45f0ce793d4313babb4ef69tlegrand@chromium.org opus_int X_offset[ VAD_N_BANDS ]; 96885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org opus_int ret = 0; 97885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org silk_VAD_state *psSilk_VAD = &psEncC->sVAD; 98e3ea049fcaee2247e45f0ce793d4313babb4ef69tlegrand@chromium.org SAVE_STACK; 99885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 100885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* Safety checks */ 101885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org silk_assert( VAD_N_BANDS == 4 ); 102885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org silk_assert( MAX_FRAME_LENGTH >= psEncC->frame_length ); 103885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org silk_assert( psEncC->frame_length <= 512 ); 104885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org silk_assert( psEncC->frame_length == 8 * silk_RSHIFT( psEncC->frame_length, 3 ) ); 105885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 106885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /***********************/ 107885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* Filter and Decimate */ 108885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /***********************/ 109e3ea049fcaee2247e45f0ce793d4313babb4ef69tlegrand@chromium.org decimated_framelength1 = silk_RSHIFT( psEncC->frame_length, 1 ); 110e3ea049fcaee2247e45f0ce793d4313babb4ef69tlegrand@chromium.org decimated_framelength2 = silk_RSHIFT( psEncC->frame_length, 2 ); 111e3ea049fcaee2247e45f0ce793d4313babb4ef69tlegrand@chromium.org decimated_framelength = silk_RSHIFT( psEncC->frame_length, 3 ); 112e3ea049fcaee2247e45f0ce793d4313babb4ef69tlegrand@chromium.org /* Decimate into 4 bands: 113e3ea049fcaee2247e45f0ce793d4313babb4ef69tlegrand@chromium.org 0 L 3L L 3L 5L 114e3ea049fcaee2247e45f0ce793d4313babb4ef69tlegrand@chromium.org - -- - -- -- 115e3ea049fcaee2247e45f0ce793d4313babb4ef69tlegrand@chromium.org 8 8 2 4 4 116e3ea049fcaee2247e45f0ce793d4313babb4ef69tlegrand@chromium.org 117e3ea049fcaee2247e45f0ce793d4313babb4ef69tlegrand@chromium.org [0-1 kHz| temp. |1-2 kHz| 2-4 kHz | 4-8 kHz | 118e3ea049fcaee2247e45f0ce793d4313babb4ef69tlegrand@chromium.org 119e3ea049fcaee2247e45f0ce793d4313babb4ef69tlegrand@chromium.org They're arranged to allow the minimal ( frame_length / 4 ) extra 120e3ea049fcaee2247e45f0ce793d4313babb4ef69tlegrand@chromium.org scratch space during the downsampling process */ 121e3ea049fcaee2247e45f0ce793d4313babb4ef69tlegrand@chromium.org X_offset[ 0 ] = 0; 122e3ea049fcaee2247e45f0ce793d4313babb4ef69tlegrand@chromium.org X_offset[ 1 ] = decimated_framelength + decimated_framelength2; 123e3ea049fcaee2247e45f0ce793d4313babb4ef69tlegrand@chromium.org X_offset[ 2 ] = X_offset[ 1 ] + decimated_framelength; 124e3ea049fcaee2247e45f0ce793d4313babb4ef69tlegrand@chromium.org X_offset[ 3 ] = X_offset[ 2 ] + decimated_framelength2; 125e3ea049fcaee2247e45f0ce793d4313babb4ef69tlegrand@chromium.org ALLOC( X, X_offset[ 3 ] + decimated_framelength1, opus_int16 ); 126e3ea049fcaee2247e45f0ce793d4313babb4ef69tlegrand@chromium.org 127885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* 0-8 kHz to 0-4 kHz and 4-8 kHz */ 128e3ea049fcaee2247e45f0ce793d4313babb4ef69tlegrand@chromium.org silk_ana_filt_bank_1( pIn, &psSilk_VAD->AnaState[ 0 ], 129e3ea049fcaee2247e45f0ce793d4313babb4ef69tlegrand@chromium.org X, &X[ X_offset[ 3 ] ], psEncC->frame_length ); 130885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 131885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* 0-4 kHz to 0-2 kHz and 2-4 kHz */ 132e3ea049fcaee2247e45f0ce793d4313babb4ef69tlegrand@chromium.org silk_ana_filt_bank_1( X, &psSilk_VAD->AnaState1[ 0 ], 133e3ea049fcaee2247e45f0ce793d4313babb4ef69tlegrand@chromium.org X, &X[ X_offset[ 2 ] ], decimated_framelength1 ); 134885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 135885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* 0-2 kHz to 0-1 kHz and 1-2 kHz */ 136e3ea049fcaee2247e45f0ce793d4313babb4ef69tlegrand@chromium.org silk_ana_filt_bank_1( X, &psSilk_VAD->AnaState2[ 0 ], 137e3ea049fcaee2247e45f0ce793d4313babb4ef69tlegrand@chromium.org X, &X[ X_offset[ 1 ] ], decimated_framelength2 ); 138885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 139885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /*********************************************/ 140885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* HP filter on lowest band (differentiator) */ 141885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /*********************************************/ 142e3ea049fcaee2247e45f0ce793d4313babb4ef69tlegrand@chromium.org X[ decimated_framelength - 1 ] = silk_RSHIFT( X[ decimated_framelength - 1 ], 1 ); 143e3ea049fcaee2247e45f0ce793d4313babb4ef69tlegrand@chromium.org HPstateTmp = X[ decimated_framelength - 1 ]; 144885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org for( i = decimated_framelength - 1; i > 0; i-- ) { 145e3ea049fcaee2247e45f0ce793d4313babb4ef69tlegrand@chromium.org X[ i - 1 ] = silk_RSHIFT( X[ i - 1 ], 1 ); 146e3ea049fcaee2247e45f0ce793d4313babb4ef69tlegrand@chromium.org X[ i ] -= X[ i - 1 ]; 147885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org } 148e3ea049fcaee2247e45f0ce793d4313babb4ef69tlegrand@chromium.org X[ 0 ] -= psSilk_VAD->HPstate; 149885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org psSilk_VAD->HPstate = HPstateTmp; 150885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 151885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /*************************************/ 152885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* Calculate the energy in each band */ 153885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /*************************************/ 154885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org for( b = 0; b < VAD_N_BANDS; b++ ) { 155885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* Find the decimated framelength in the non-uniformly divided bands */ 156885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org decimated_framelength = silk_RSHIFT( psEncC->frame_length, silk_min_int( VAD_N_BANDS - b, VAD_N_BANDS - 1 ) ); 157885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 158885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* Split length into subframe lengths */ 159885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org dec_subframe_length = silk_RSHIFT( decimated_framelength, VAD_INTERNAL_SUBFRAMES_LOG2 ); 160885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org dec_subframe_offset = 0; 161885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 162885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* Compute energy per sub-frame */ 163885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* initialize with summed energy of last subframe */ 164885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org Xnrg[ b ] = psSilk_VAD->XnrgSubfr[ b ]; 165885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org for( s = 0; s < VAD_INTERNAL_SUBFRAMES; s++ ) { 166885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org sumSquared = 0; 167885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org for( i = 0; i < dec_subframe_length; i++ ) { 168885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* The energy will be less than dec_subframe_length * ( silk_int16_MIN / 8 ) ^ 2. */ 169885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* Therefore we can accumulate with no risk of overflow (unless dec_subframe_length > 128) */ 170e3ea049fcaee2247e45f0ce793d4313babb4ef69tlegrand@chromium.org x_tmp = silk_RSHIFT( 171e3ea049fcaee2247e45f0ce793d4313babb4ef69tlegrand@chromium.org X[ X_offset[ b ] + i + dec_subframe_offset ], 3 ); 172885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org sumSquared = silk_SMLABB( sumSquared, x_tmp, x_tmp ); 173885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 174885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* Safety check */ 175885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org silk_assert( sumSquared >= 0 ); 176885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org } 177885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 178885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* Add/saturate summed energy of current subframe */ 179885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org if( s < VAD_INTERNAL_SUBFRAMES - 1 ) { 180885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org Xnrg[ b ] = silk_ADD_POS_SAT32( Xnrg[ b ], sumSquared ); 181885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org } else { 182885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* Look-ahead subframe */ 183885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org Xnrg[ b ] = silk_ADD_POS_SAT32( Xnrg[ b ], silk_RSHIFT( sumSquared, 1 ) ); 184885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org } 185885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 186885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org dec_subframe_offset += dec_subframe_length; 187885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org } 188885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org psSilk_VAD->XnrgSubfr[ b ] = sumSquared; 189885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org } 190885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 191885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /********************/ 192885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* Noise estimation */ 193885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /********************/ 194885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org silk_VAD_GetNoiseLevels( &Xnrg[ 0 ], psSilk_VAD ); 195885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 196885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /***********************************************/ 197885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* Signal-plus-noise to noise ratio estimation */ 198885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /***********************************************/ 199885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org sumSquared = 0; 200885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org input_tilt = 0; 201885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org for( b = 0; b < VAD_N_BANDS; b++ ) { 202885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org speech_nrg = Xnrg[ b ] - psSilk_VAD->NL[ b ]; 203885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org if( speech_nrg > 0 ) { 204885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* Divide, with sufficient resolution */ 205885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org if( ( Xnrg[ b ] & 0xFF800000 ) == 0 ) { 206885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org NrgToNoiseRatio_Q8[ b ] = silk_DIV32( silk_LSHIFT( Xnrg[ b ], 8 ), psSilk_VAD->NL[ b ] + 1 ); 207885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org } else { 208885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org NrgToNoiseRatio_Q8[ b ] = silk_DIV32( Xnrg[ b ], silk_RSHIFT( psSilk_VAD->NL[ b ], 8 ) + 1 ); 209885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org } 210885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 211885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* Convert to log domain */ 212885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org SNR_Q7 = silk_lin2log( NrgToNoiseRatio_Q8[ b ] ) - 8 * 128; 213885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 214885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* Sum-of-squares */ 215885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org sumSquared = silk_SMLABB( sumSquared, SNR_Q7, SNR_Q7 ); /* Q14 */ 216885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 217885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* Tilt measure */ 218885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org if( speech_nrg < ( (opus_int32)1 << 20 ) ) { 219885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* Scale down SNR value for small subband speech energies */ 220885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org SNR_Q7 = silk_SMULWB( silk_LSHIFT( silk_SQRT_APPROX( speech_nrg ), 6 ), SNR_Q7 ); 221885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org } 222885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org input_tilt = silk_SMLAWB( input_tilt, tiltWeights[ b ], SNR_Q7 ); 223885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org } else { 224885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org NrgToNoiseRatio_Q8[ b ] = 256; 225885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org } 226885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org } 227885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 228885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* Mean-of-squares */ 229885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org sumSquared = silk_DIV32_16( sumSquared, VAD_N_BANDS ); /* Q14 */ 230885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 231885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* Root-mean-square approximation, scale to dBs, and write to output pointer */ 232885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org pSNR_dB_Q7 = (opus_int16)( 3 * silk_SQRT_APPROX( sumSquared ) ); /* Q7 */ 233885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 234885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /*********************************/ 235885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* Speech Probability Estimation */ 236885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /*********************************/ 237885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org SA_Q15 = silk_sigm_Q15( silk_SMULWB( VAD_SNR_FACTOR_Q16, pSNR_dB_Q7 ) - VAD_NEGATIVE_OFFSET_Q5 ); 238885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 239885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /**************************/ 240885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* Frequency Tilt Measure */ 241885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /**************************/ 242885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org psEncC->input_tilt_Q15 = silk_LSHIFT( silk_sigm_Q15( input_tilt ) - 16384, 1 ); 243885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 244885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /**************************************************/ 245885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* Scale the sigmoid output based on power levels */ 246885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /**************************************************/ 247885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org speech_nrg = 0; 248885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org for( b = 0; b < VAD_N_BANDS; b++ ) { 249885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* Accumulate signal-without-noise energies, higher frequency bands have more weight */ 250885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org speech_nrg += ( b + 1 ) * silk_RSHIFT( Xnrg[ b ] - psSilk_VAD->NL[ b ], 4 ); 251885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org } 252885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 253885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* Power scaling */ 254885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org if( speech_nrg <= 0 ) { 255885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org SA_Q15 = silk_RSHIFT( SA_Q15, 1 ); 256885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org } else if( speech_nrg < 32768 ) { 257885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org if( psEncC->frame_length == 10 * psEncC->fs_kHz ) { 258885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org speech_nrg = silk_LSHIFT_SAT32( speech_nrg, 16 ); 259885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org } else { 260885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org speech_nrg = silk_LSHIFT_SAT32( speech_nrg, 15 ); 261885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org } 262885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 263885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* square-root */ 264885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org speech_nrg = silk_SQRT_APPROX( speech_nrg ); 265885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org SA_Q15 = silk_SMULWB( 32768 + speech_nrg, SA_Q15 ); 266885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org } 267885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 268885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* Copy the resulting speech activity in Q8 */ 269885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org psEncC->speech_activity_Q8 = silk_min_int( silk_RSHIFT( SA_Q15, 7 ), silk_uint8_MAX ); 270885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 271885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /***********************************/ 272885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* Energy Level and SNR estimation */ 273885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /***********************************/ 274885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* Smoothing coefficient */ 275885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org smooth_coef_Q16 = silk_SMULWB( VAD_SNR_SMOOTH_COEF_Q18, silk_SMULWB( (opus_int32)SA_Q15, SA_Q15 ) ); 276885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 277885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org if( psEncC->frame_length == 10 * psEncC->fs_kHz ) { 278885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org smooth_coef_Q16 >>= 1; 279885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org } 280885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 281885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org for( b = 0; b < VAD_N_BANDS; b++ ) { 282885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* compute smoothed energy-to-noise ratio per band */ 283885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org psSilk_VAD->NrgRatioSmth_Q8[ b ] = silk_SMLAWB( psSilk_VAD->NrgRatioSmth_Q8[ b ], 284885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org NrgToNoiseRatio_Q8[ b ] - psSilk_VAD->NrgRatioSmth_Q8[ b ], smooth_coef_Q16 ); 285885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 286885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* signal to noise ratio in dB per band */ 287885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org SNR_Q7 = 3 * ( silk_lin2log( psSilk_VAD->NrgRatioSmth_Q8[b] ) - 8 * 128 ); 288885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* quality = sigmoid( 0.25 * ( SNR_dB - 16 ) ); */ 289885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org psEncC->input_quality_bands_Q15[ b ] = silk_sigm_Q15( silk_RSHIFT( SNR_Q7 - 16 * 128, 4 ) ); 290885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org } 291885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 292e3ea049fcaee2247e45f0ce793d4313babb4ef69tlegrand@chromium.org RESTORE_STACK; 293885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org return( ret ); 294885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org} 295885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 296885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org/**************************/ 297885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org/* Noise level estimation */ 298885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org/**************************/ 2993c3902f0ac13428394f14f78f0fab05ef3468d69tlegrand@google.comstatic OPUS_INLINE void silk_VAD_GetNoiseLevels( 300885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org const opus_int32 pX[ VAD_N_BANDS ], /* I subband energies */ 301885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org silk_VAD_state *psSilk_VAD /* I/O Pointer to Silk VAD state */ 302885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org) 303885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org{ 304885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org opus_int k; 305885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org opus_int32 nl, nrg, inv_nrg; 306885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org opus_int coef, min_coef; 307885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 308885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* Initially faster smoothing */ 309885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org if( psSilk_VAD->counter < 1000 ) { /* 1000 = 20 sec */ 310885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org min_coef = silk_DIV32_16( silk_int16_MAX, silk_RSHIFT( psSilk_VAD->counter, 4 ) + 1 ); 311885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org } else { 312885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org min_coef = 0; 313885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org } 314885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 315885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org for( k = 0; k < VAD_N_BANDS; k++ ) { 316885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* Get old noise level estimate for current band */ 317885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org nl = psSilk_VAD->NL[ k ]; 318885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org silk_assert( nl >= 0 ); 319885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 320885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* Add bias */ 321885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org nrg = silk_ADD_POS_SAT32( pX[ k ], psSilk_VAD->NoiseLevelBias[ k ] ); 322885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org silk_assert( nrg > 0 ); 323885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 324885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* Invert energies */ 325885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org inv_nrg = silk_DIV32( silk_int32_MAX, nrg ); 326885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org silk_assert( inv_nrg >= 0 ); 327885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 328885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* Less update when subband energy is high */ 329885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org if( nrg > silk_LSHIFT( nl, 3 ) ) { 330885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org coef = VAD_NOISE_LEVEL_SMOOTH_COEF_Q16 >> 3; 331885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org } else if( nrg < nl ) { 332885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org coef = VAD_NOISE_LEVEL_SMOOTH_COEF_Q16; 333885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org } else { 334885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org coef = silk_SMULWB( silk_SMULWW( inv_nrg, nl ), VAD_NOISE_LEVEL_SMOOTH_COEF_Q16 << 1 ); 335885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org } 336885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 337885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* Initially faster smoothing */ 338885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org coef = silk_max_int( coef, min_coef ); 339885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 340885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* Smooth inverse energies */ 341885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org psSilk_VAD->inv_NL[ k ] = silk_SMLAWB( psSilk_VAD->inv_NL[ k ], inv_nrg - psSilk_VAD->inv_NL[ k ], coef ); 342885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org silk_assert( psSilk_VAD->inv_NL[ k ] >= 0 ); 343885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 344885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* Compute noise level by inverting again */ 345885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org nl = silk_DIV32( silk_int32_MAX, psSilk_VAD->inv_NL[ k ] ); 346885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org silk_assert( nl >= 0 ); 347885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 348885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* Limit noise levels (guarantee 7 bits of head room) */ 349885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org nl = silk_min( nl, 0x00FFFFFF ); 350885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 351885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* Store as part of state */ 352885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org psSilk_VAD->NL[ k ] = nl; 353885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org } 354885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org 355885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org /* Increment frame counter */ 356885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org psSilk_VAD->counter++; 357885f2ff5a7a7d6a73432d26a6c0ae9147e6b452sergeyu@chromium.org} 358