media/libopus/silk/decode_parameters.c

changeset 0
6474c204b198
     1.1 --- /dev/null	Thu Jan 01 00:00:00 1970 +0000
     1.2 +++ b/media/libopus/silk/decode_parameters.c	Wed Dec 31 06:09:35 2014 +0100
     1.3 @@ -0,0 +1,115 @@
     1.4 +/***********************************************************************
     1.5 +Copyright (c) 2006-2011, Skype Limited. All rights reserved.
     1.6 +Redistribution and use in source and binary forms, with or without
     1.7 +modification, are permitted provided that the following conditions
     1.8 +are met:
     1.9 +- Redistributions of source code must retain the above copyright notice,
    1.10 +this list of conditions and the following disclaimer.
    1.11 +- Redistributions in binary form must reproduce the above copyright
    1.12 +notice, this list of conditions and the following disclaimer in the
    1.13 +documentation and/or other materials provided with the distribution.
    1.14 +- Neither the name of Internet Society, IETF or IETF Trust, nor the
    1.15 +names of specific contributors, may be used to endorse or promote
    1.16 +products derived from this software without specific prior written
    1.17 +permission.
    1.18 +THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
    1.19 +AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
    1.20 +IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
    1.21 +ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
    1.22 +LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
    1.23 +CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
    1.24 +SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
    1.25 +INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
    1.26 +CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
    1.27 +ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
    1.28 +POSSIBILITY OF SUCH DAMAGE.
    1.29 +***********************************************************************/
    1.30 +
    1.31 +#ifdef HAVE_CONFIG_H
    1.32 +#include "config.h"
    1.33 +#endif
    1.34 +
    1.35 +#include "main.h"
    1.36 +
    1.37 +/* Decode parameters from payload */
    1.38 +void silk_decode_parameters(
    1.39 +    silk_decoder_state          *psDec,                         /* I/O  State                                       */
    1.40 +    silk_decoder_control        *psDecCtrl,                     /* I/O  Decoder control                             */
    1.41 +    opus_int                    condCoding                      /* I    The type of conditional coding to use       */
    1.42 +)
    1.43 +{
    1.44 +    opus_int   i, k, Ix;
    1.45 +    opus_int16 pNLSF_Q15[ MAX_LPC_ORDER ], pNLSF0_Q15[ MAX_LPC_ORDER ];
    1.46 +    const opus_int8 *cbk_ptr_Q7;
    1.47 +
    1.48 +    /* Dequant Gains */
    1.49 +    silk_gains_dequant( psDecCtrl->Gains_Q16, psDec->indices.GainsIndices,
    1.50 +        &psDec->LastGainIndex, condCoding == CODE_CONDITIONALLY, psDec->nb_subfr );
    1.51 +
    1.52 +    /****************/
    1.53 +    /* Decode NLSFs */
    1.54 +    /****************/
    1.55 +    silk_NLSF_decode( pNLSF_Q15, psDec->indices.NLSFIndices, psDec->psNLSF_CB );
    1.56 +
    1.57 +    /* Convert NLSF parameters to AR prediction filter coefficients */
    1.58 +    silk_NLSF2A( psDecCtrl->PredCoef_Q12[ 1 ], pNLSF_Q15, psDec->LPC_order );
    1.59 +
    1.60 +    /* If just reset, e.g., because internal Fs changed, do not allow interpolation */
    1.61 +    /* improves the case of packet loss in the first frame after a switch           */
    1.62 +    if( psDec->first_frame_after_reset == 1 ) {
    1.63 +        psDec->indices.NLSFInterpCoef_Q2 = 4;
    1.64 +    }
    1.65 +
    1.66 +    if( psDec->indices.NLSFInterpCoef_Q2 < 4 ) {
    1.67 +        /* Calculation of the interpolated NLSF0 vector from the interpolation factor, */
    1.68 +        /* the previous NLSF1, and the current NLSF1                                   */
    1.69 +        for( i = 0; i < psDec->LPC_order; i++ ) {
    1.70 +            pNLSF0_Q15[ i ] = psDec->prevNLSF_Q15[ i ] + silk_RSHIFT( silk_MUL( psDec->indices.NLSFInterpCoef_Q2,
    1.71 +                pNLSF_Q15[ i ] - psDec->prevNLSF_Q15[ i ] ), 2 );
    1.72 +        }
    1.73 +
    1.74 +        /* Convert NLSF parameters to AR prediction filter coefficients */
    1.75 +        silk_NLSF2A( psDecCtrl->PredCoef_Q12[ 0 ], pNLSF0_Q15, psDec->LPC_order );
    1.76 +    } else {
    1.77 +        /* Copy LPC coefficients for first half from second half */
    1.78 +        silk_memcpy( psDecCtrl->PredCoef_Q12[ 0 ], psDecCtrl->PredCoef_Q12[ 1 ], psDec->LPC_order * sizeof( opus_int16 ) );
    1.79 +    }
    1.80 +
    1.81 +    silk_memcpy( psDec->prevNLSF_Q15, pNLSF_Q15, psDec->LPC_order * sizeof( opus_int16 ) );
    1.82 +
    1.83 +    /* After a packet loss do BWE of LPC coefs */
    1.84 +    if( psDec->lossCnt ) {
    1.85 +        silk_bwexpander( psDecCtrl->PredCoef_Q12[ 0 ], psDec->LPC_order, BWE_AFTER_LOSS_Q16 );
    1.86 +        silk_bwexpander( psDecCtrl->PredCoef_Q12[ 1 ], psDec->LPC_order, BWE_AFTER_LOSS_Q16 );
    1.87 +    }
    1.88 +
    1.89 +    if( psDec->indices.signalType == TYPE_VOICED ) {
    1.90 +        /*********************/
    1.91 +        /* Decode pitch lags */
    1.92 +        /*********************/
    1.93 +
    1.94 +        /* Decode pitch values */
    1.95 +        silk_decode_pitch( psDec->indices.lagIndex, psDec->indices.contourIndex, psDecCtrl->pitchL, psDec->fs_kHz, psDec->nb_subfr );
    1.96 +
    1.97 +        /* Decode Codebook Index */
    1.98 +        cbk_ptr_Q7 = silk_LTP_vq_ptrs_Q7[ psDec->indices.PERIndex ]; /* set pointer to start of codebook */
    1.99 +
   1.100 +        for( k = 0; k < psDec->nb_subfr; k++ ) {
   1.101 +            Ix = psDec->indices.LTPIndex[ k ];
   1.102 +            for( i = 0; i < LTP_ORDER; i++ ) {
   1.103 +                psDecCtrl->LTPCoef_Q14[ k * LTP_ORDER + i ] = silk_LSHIFT( cbk_ptr_Q7[ Ix * LTP_ORDER + i ], 7 );
   1.104 +            }
   1.105 +        }
   1.106 +
   1.107 +        /**********************/
   1.108 +        /* Decode LTP scaling */
   1.109 +        /**********************/
   1.110 +        Ix = psDec->indices.LTP_scaleIndex;
   1.111 +        psDecCtrl->LTP_scale_Q14 = silk_LTPScales_table_Q14[ Ix ];
   1.112 +    } else {
   1.113 +        silk_memset( psDecCtrl->pitchL,      0,             psDec->nb_subfr * sizeof( opus_int   ) );
   1.114 +        silk_memset( psDecCtrl->LTPCoef_Q14, 0, LTP_ORDER * psDec->nb_subfr * sizeof( opus_int16 ) );
   1.115 +        psDec->indices.PERIndex  = 0;
   1.116 +        psDecCtrl->LTP_scale_Q14 = 0;
   1.117 +    }
   1.118 +}

mercurial