annotate src/opus-1.3/silk/decode_parameters.c @ 69:7aeed7906520

Add Opus sources and macOS builds
author Chris Cannam
date Wed, 23 Jan 2019 13:48:08 +0000
parents
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 "main.h"
Chris@69 33
Chris@69 34 /* Decode parameters from payload */
Chris@69 35 void silk_decode_parameters(
Chris@69 36 silk_decoder_state *psDec, /* I/O State */
Chris@69 37 silk_decoder_control *psDecCtrl, /* I/O Decoder control */
Chris@69 38 opus_int condCoding /* I The type of conditional coding to use */
Chris@69 39 )
Chris@69 40 {
Chris@69 41 opus_int i, k, Ix;
Chris@69 42 opus_int16 pNLSF_Q15[ MAX_LPC_ORDER ], pNLSF0_Q15[ MAX_LPC_ORDER ];
Chris@69 43 const opus_int8 *cbk_ptr_Q7;
Chris@69 44
Chris@69 45 /* Dequant Gains */
Chris@69 46 silk_gains_dequant( psDecCtrl->Gains_Q16, psDec->indices.GainsIndices,
Chris@69 47 &psDec->LastGainIndex, condCoding == CODE_CONDITIONALLY, psDec->nb_subfr );
Chris@69 48
Chris@69 49 /****************/
Chris@69 50 /* Decode NLSFs */
Chris@69 51 /****************/
Chris@69 52 silk_NLSF_decode( pNLSF_Q15, psDec->indices.NLSFIndices, psDec->psNLSF_CB );
Chris@69 53
Chris@69 54 /* Convert NLSF parameters to AR prediction filter coefficients */
Chris@69 55 silk_NLSF2A( psDecCtrl->PredCoef_Q12[ 1 ], pNLSF_Q15, psDec->LPC_order, psDec->arch );
Chris@69 56
Chris@69 57 /* If just reset, e.g., because internal Fs changed, do not allow interpolation */
Chris@69 58 /* improves the case of packet loss in the first frame after a switch */
Chris@69 59 if( psDec->first_frame_after_reset == 1 ) {
Chris@69 60 psDec->indices.NLSFInterpCoef_Q2 = 4;
Chris@69 61 }
Chris@69 62
Chris@69 63 if( psDec->indices.NLSFInterpCoef_Q2 < 4 ) {
Chris@69 64 /* Calculation of the interpolated NLSF0 vector from the interpolation factor, */
Chris@69 65 /* the previous NLSF1, and the current NLSF1 */
Chris@69 66 for( i = 0; i < psDec->LPC_order; i++ ) {
Chris@69 67 pNLSF0_Q15[ i ] = psDec->prevNLSF_Q15[ i ] + silk_RSHIFT( silk_MUL( psDec->indices.NLSFInterpCoef_Q2,
Chris@69 68 pNLSF_Q15[ i ] - psDec->prevNLSF_Q15[ i ] ), 2 );
Chris@69 69 }
Chris@69 70
Chris@69 71 /* Convert NLSF parameters to AR prediction filter coefficients */
Chris@69 72 silk_NLSF2A( psDecCtrl->PredCoef_Q12[ 0 ], pNLSF0_Q15, psDec->LPC_order, psDec->arch );
Chris@69 73 } else {
Chris@69 74 /* Copy LPC coefficients for first half from second half */
Chris@69 75 silk_memcpy( psDecCtrl->PredCoef_Q12[ 0 ], psDecCtrl->PredCoef_Q12[ 1 ], psDec->LPC_order * sizeof( opus_int16 ) );
Chris@69 76 }
Chris@69 77
Chris@69 78 silk_memcpy( psDec->prevNLSF_Q15, pNLSF_Q15, psDec->LPC_order * sizeof( opus_int16 ) );
Chris@69 79
Chris@69 80 /* After a packet loss do BWE of LPC coefs */
Chris@69 81 if( psDec->lossCnt ) {
Chris@69 82 silk_bwexpander( psDecCtrl->PredCoef_Q12[ 0 ], psDec->LPC_order, BWE_AFTER_LOSS_Q16 );
Chris@69 83 silk_bwexpander( psDecCtrl->PredCoef_Q12[ 1 ], psDec->LPC_order, BWE_AFTER_LOSS_Q16 );
Chris@69 84 }
Chris@69 85
Chris@69 86 if( psDec->indices.signalType == TYPE_VOICED ) {
Chris@69 87 /*********************/
Chris@69 88 /* Decode pitch lags */
Chris@69 89 /*********************/
Chris@69 90
Chris@69 91 /* Decode pitch values */
Chris@69 92 silk_decode_pitch( psDec->indices.lagIndex, psDec->indices.contourIndex, psDecCtrl->pitchL, psDec->fs_kHz, psDec->nb_subfr );
Chris@69 93
Chris@69 94 /* Decode Codebook Index */
Chris@69 95 cbk_ptr_Q7 = silk_LTP_vq_ptrs_Q7[ psDec->indices.PERIndex ]; /* set pointer to start of codebook */
Chris@69 96
Chris@69 97 for( k = 0; k < psDec->nb_subfr; k++ ) {
Chris@69 98 Ix = psDec->indices.LTPIndex[ k ];
Chris@69 99 for( i = 0; i < LTP_ORDER; i++ ) {
Chris@69 100 psDecCtrl->LTPCoef_Q14[ k * LTP_ORDER + i ] = silk_LSHIFT( cbk_ptr_Q7[ Ix * LTP_ORDER + i ], 7 );
Chris@69 101 }
Chris@69 102 }
Chris@69 103
Chris@69 104 /**********************/
Chris@69 105 /* Decode LTP scaling */
Chris@69 106 /**********************/
Chris@69 107 Ix = psDec->indices.LTP_scaleIndex;
Chris@69 108 psDecCtrl->LTP_scale_Q14 = silk_LTPScales_table_Q14[ Ix ];
Chris@69 109 } else {
Chris@69 110 silk_memset( psDecCtrl->pitchL, 0, psDec->nb_subfr * sizeof( opus_int ) );
Chris@69 111 silk_memset( psDecCtrl->LTPCoef_Q14, 0, LTP_ORDER * psDec->nb_subfr * sizeof( opus_int16 ) );
Chris@69 112 psDec->indices.PERIndex = 0;
Chris@69 113 psDecCtrl->LTP_scale_Q14 = 0;
Chris@69 114 }
Chris@69 115 }