media/libopus/silk/decode_parameters.c

Thu, 22 Jan 2015 13:21:57 +0100

author
Michael Schloh von Bennewitz <michael@schloh.com>
date
Thu, 22 Jan 2015 13:21:57 +0100
branch
TOR_BUG_9701
changeset 15
b8a032363ba2
permissions
-rw-r--r--

Incorporate requested changes from Mozilla in review:
https://bugzilla.mozilla.org/show_bug.cgi?id=1123480#c6

michael@0 1 /***********************************************************************
michael@0 2 Copyright (c) 2006-2011, Skype Limited. All rights reserved.
michael@0 3 Redistribution and use in source and binary forms, with or without
michael@0 4 modification, are permitted provided that the following conditions
michael@0 5 are met:
michael@0 6 - Redistributions of source code must retain the above copyright notice,
michael@0 7 this list of conditions and the following disclaimer.
michael@0 8 - Redistributions in binary form must reproduce the above copyright
michael@0 9 notice, this list of conditions and the following disclaimer in the
michael@0 10 documentation and/or other materials provided with the distribution.
michael@0 11 - Neither the name of Internet Society, IETF or IETF Trust, nor the
michael@0 12 names of specific contributors, may be used to endorse or promote
michael@0 13 products derived from this software without specific prior written
michael@0 14 permission.
michael@0 15 THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
michael@0 16 AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
michael@0 17 IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
michael@0 18 ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
michael@0 19 LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
michael@0 20 CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
michael@0 21 SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
michael@0 22 INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
michael@0 23 CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
michael@0 24 ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
michael@0 25 POSSIBILITY OF SUCH DAMAGE.
michael@0 26 ***********************************************************************/
michael@0 27
michael@0 28 #ifdef HAVE_CONFIG_H
michael@0 29 #include "config.h"
michael@0 30 #endif
michael@0 31
michael@0 32 #include "main.h"
michael@0 33
michael@0 34 /* Decode parameters from payload */
michael@0 35 void silk_decode_parameters(
michael@0 36 silk_decoder_state *psDec, /* I/O State */
michael@0 37 silk_decoder_control *psDecCtrl, /* I/O Decoder control */
michael@0 38 opus_int condCoding /* I The type of conditional coding to use */
michael@0 39 )
michael@0 40 {
michael@0 41 opus_int i, k, Ix;
michael@0 42 opus_int16 pNLSF_Q15[ MAX_LPC_ORDER ], pNLSF0_Q15[ MAX_LPC_ORDER ];
michael@0 43 const opus_int8 *cbk_ptr_Q7;
michael@0 44
michael@0 45 /* Dequant Gains */
michael@0 46 silk_gains_dequant( psDecCtrl->Gains_Q16, psDec->indices.GainsIndices,
michael@0 47 &psDec->LastGainIndex, condCoding == CODE_CONDITIONALLY, psDec->nb_subfr );
michael@0 48
michael@0 49 /****************/
michael@0 50 /* Decode NLSFs */
michael@0 51 /****************/
michael@0 52 silk_NLSF_decode( pNLSF_Q15, psDec->indices.NLSFIndices, psDec->psNLSF_CB );
michael@0 53
michael@0 54 /* Convert NLSF parameters to AR prediction filter coefficients */
michael@0 55 silk_NLSF2A( psDecCtrl->PredCoef_Q12[ 1 ], pNLSF_Q15, psDec->LPC_order );
michael@0 56
michael@0 57 /* If just reset, e.g., because internal Fs changed, do not allow interpolation */
michael@0 58 /* improves the case of packet loss in the first frame after a switch */
michael@0 59 if( psDec->first_frame_after_reset == 1 ) {
michael@0 60 psDec->indices.NLSFInterpCoef_Q2 = 4;
michael@0 61 }
michael@0 62
michael@0 63 if( psDec->indices.NLSFInterpCoef_Q2 < 4 ) {
michael@0 64 /* Calculation of the interpolated NLSF0 vector from the interpolation factor, */
michael@0 65 /* the previous NLSF1, and the current NLSF1 */
michael@0 66 for( i = 0; i < psDec->LPC_order; i++ ) {
michael@0 67 pNLSF0_Q15[ i ] = psDec->prevNLSF_Q15[ i ] + silk_RSHIFT( silk_MUL( psDec->indices.NLSFInterpCoef_Q2,
michael@0 68 pNLSF_Q15[ i ] - psDec->prevNLSF_Q15[ i ] ), 2 );
michael@0 69 }
michael@0 70
michael@0 71 /* Convert NLSF parameters to AR prediction filter coefficients */
michael@0 72 silk_NLSF2A( psDecCtrl->PredCoef_Q12[ 0 ], pNLSF0_Q15, psDec->LPC_order );
michael@0 73 } else {
michael@0 74 /* Copy LPC coefficients for first half from second half */
michael@0 75 silk_memcpy( psDecCtrl->PredCoef_Q12[ 0 ], psDecCtrl->PredCoef_Q12[ 1 ], psDec->LPC_order * sizeof( opus_int16 ) );
michael@0 76 }
michael@0 77
michael@0 78 silk_memcpy( psDec->prevNLSF_Q15, pNLSF_Q15, psDec->LPC_order * sizeof( opus_int16 ) );
michael@0 79
michael@0 80 /* After a packet loss do BWE of LPC coefs */
michael@0 81 if( psDec->lossCnt ) {
michael@0 82 silk_bwexpander( psDecCtrl->PredCoef_Q12[ 0 ], psDec->LPC_order, BWE_AFTER_LOSS_Q16 );
michael@0 83 silk_bwexpander( psDecCtrl->PredCoef_Q12[ 1 ], psDec->LPC_order, BWE_AFTER_LOSS_Q16 );
michael@0 84 }
michael@0 85
michael@0 86 if( psDec->indices.signalType == TYPE_VOICED ) {
michael@0 87 /*********************/
michael@0 88 /* Decode pitch lags */
michael@0 89 /*********************/
michael@0 90
michael@0 91 /* Decode pitch values */
michael@0 92 silk_decode_pitch( psDec->indices.lagIndex, psDec->indices.contourIndex, psDecCtrl->pitchL, psDec->fs_kHz, psDec->nb_subfr );
michael@0 93
michael@0 94 /* Decode Codebook Index */
michael@0 95 cbk_ptr_Q7 = silk_LTP_vq_ptrs_Q7[ psDec->indices.PERIndex ]; /* set pointer to start of codebook */
michael@0 96
michael@0 97 for( k = 0; k < psDec->nb_subfr; k++ ) {
michael@0 98 Ix = psDec->indices.LTPIndex[ k ];
michael@0 99 for( i = 0; i < LTP_ORDER; i++ ) {
michael@0 100 psDecCtrl->LTPCoef_Q14[ k * LTP_ORDER + i ] = silk_LSHIFT( cbk_ptr_Q7[ Ix * LTP_ORDER + i ], 7 );
michael@0 101 }
michael@0 102 }
michael@0 103
michael@0 104 /**********************/
michael@0 105 /* Decode LTP scaling */
michael@0 106 /**********************/
michael@0 107 Ix = psDec->indices.LTP_scaleIndex;
michael@0 108 psDecCtrl->LTP_scale_Q14 = silk_LTPScales_table_Q14[ Ix ];
michael@0 109 } else {
michael@0 110 silk_memset( psDecCtrl->pitchL, 0, psDec->nb_subfr * sizeof( opus_int ) );
michael@0 111 silk_memset( psDecCtrl->LTPCoef_Q14, 0, LTP_ORDER * psDec->nb_subfr * sizeof( opus_int16 ) );
michael@0 112 psDec->indices.PERIndex = 0;
michael@0 113 psDecCtrl->LTP_scale_Q14 = 0;
michael@0 114 }
michael@0 115 }

mercurial