| OLD | NEW |
| 1 /*********************************************************************** | 1 /*********************************************************************** |
| 2 Copyright (c) 2006-2011, Skype Limited. All rights reserved. | 2 Copyright (c) 2006-2011, Skype Limited. All rights reserved. |
| 3 Redistribution and use in source and binary forms, with or without | 3 Redistribution and use in source and binary forms, with or without |
| 4 modification, are permitted provided that the following conditions | 4 modification, are permitted provided that the following conditions |
| 5 are met: | 5 are met: |
| 6 - Redistributions of source code must retain the above copyright notice, | 6 - Redistributions of source code must retain the above copyright notice, |
| 7 this list of conditions and the following disclaimer. | 7 this list of conditions and the following disclaimer. |
| 8 - Redistributions in binary form must reproduce the above copyright | 8 - Redistributions in binary form must reproduce the above copyright |
| 9 notice, this list of conditions and the following disclaimer in the | 9 notice, this list of conditions and the following disclaimer in the |
| 10 documentation and/or other materials provided with the distribution. | 10 documentation and/or other materials provided with the distribution. |
| (...skipping 55 matching lines...) Expand 10 before | Expand all | Expand 10 after Loading... |
| 66 /* Convert to NLSFs */ | 66 /* Convert to NLSFs */ |
| 67 silk_A2NLSF_FLP( NLSF_Q15, a_tmp, psEncC->predictLPCOrder ); | 67 silk_A2NLSF_FLP( NLSF_Q15, a_tmp, psEncC->predictLPCOrder ); |
| 68 | 68 |
| 69 /* Search over interpolation indices to find the one with lowest residua
l energy */ | 69 /* Search over interpolation indices to find the one with lowest residua
l energy */ |
| 70 res_nrg_2nd = silk_float_MAX; | 70 res_nrg_2nd = silk_float_MAX; |
| 71 for( k = 3; k >= 0; k-- ) { | 71 for( k = 3; k >= 0; k-- ) { |
| 72 /* Interpolate NLSFs for first half */ | 72 /* Interpolate NLSFs for first half */ |
| 73 silk_interpolate( NLSF0_Q15, psEncC->prev_NLSFq_Q15, NLSF_Q15, k, ps
EncC->predictLPCOrder ); | 73 silk_interpolate( NLSF0_Q15, psEncC->prev_NLSFq_Q15, NLSF_Q15, k, ps
EncC->predictLPCOrder ); |
| 74 | 74 |
| 75 /* Convert to LPC for residual energy evaluation */ | 75 /* Convert to LPC for residual energy evaluation */ |
| 76 silk_NLSF2A_FLP( a_tmp, NLSF0_Q15, psEncC->predictLPCOrder ); | 76 silk_NLSF2A_FLP( a_tmp, NLSF0_Q15, psEncC->predictLPCOrder, psEncC->
arch ); |
| 77 | 77 |
| 78 /* Calculate residual energy with LSF interpolation */ | 78 /* Calculate residual energy with LSF interpolation */ |
| 79 silk_LPC_analysis_filter_FLP( LPC_res, a_tmp, x, 2 * subfr_length, p
sEncC->predictLPCOrder ); | 79 silk_LPC_analysis_filter_FLP( LPC_res, a_tmp, x, 2 * subfr_length, p
sEncC->predictLPCOrder ); |
| 80 res_nrg_interp = (silk_float)( | 80 res_nrg_interp = (silk_float)( |
| 81 silk_energy_FLP( LPC_res + psEncC->predictLPCOrder,
subfr_length - psEncC->predictLPCOrder ) + | 81 silk_energy_FLP( LPC_res + psEncC->predictLPCOrder,
subfr_length - psEncC->predictLPCOrder ) + |
| 82 silk_energy_FLP( LPC_res + psEncC->predictLPCOrder + subfr_lengt
h, subfr_length - psEncC->predictLPCOrder ) ); | 82 silk_energy_FLP( LPC_res + psEncC->predictLPCOrder + subfr_lengt
h, subfr_length - psEncC->predictLPCOrder ) ); |
| 83 | 83 |
| 84 /* Determine whether current interpolated NLSFs are best so far */ | 84 /* Determine whether current interpolated NLSFs are best so far */ |
| 85 if( res_nrg_interp < res_nrg ) { | 85 if( res_nrg_interp < res_nrg ) { |
| 86 /* Interpolation has lower residual energy */ | 86 /* Interpolation has lower residual energy */ |
| 87 res_nrg = res_nrg_interp; | 87 res_nrg = res_nrg_interp; |
| 88 psEncC->indices.NLSFInterpCoef_Q2 = (opus_int8)k; | 88 psEncC->indices.NLSFInterpCoef_Q2 = (opus_int8)k; |
| 89 } else if( res_nrg_interp > res_nrg_2nd ) { | 89 } else if( res_nrg_interp > res_nrg_2nd ) { |
| 90 /* No reason to continue iterating - residual energies will cont
inue to climb */ | 90 /* No reason to continue iterating - residual energies will cont
inue to climb */ |
| 91 break; | 91 break; |
| 92 } | 92 } |
| 93 res_nrg_2nd = res_nrg_interp; | 93 res_nrg_2nd = res_nrg_interp; |
| 94 } | 94 } |
| 95 } | 95 } |
| 96 | 96 |
| 97 if( psEncC->indices.NLSFInterpCoef_Q2 == 4 ) { | 97 if( psEncC->indices.NLSFInterpCoef_Q2 == 4 ) { |
| 98 /* NLSF interpolation is currently inactive, calculate NLSFs from full f
rame AR coefficients */ | 98 /* NLSF interpolation is currently inactive, calculate NLSFs from full f
rame AR coefficients */ |
| 99 silk_A2NLSF_FLP( NLSF_Q15, a, psEncC->predictLPCOrder ); | 99 silk_A2NLSF_FLP( NLSF_Q15, a, psEncC->predictLPCOrder ); |
| 100 } | 100 } |
| 101 | 101 |
| 102 silk_assert( psEncC->indices.NLSFInterpCoef_Q2 == 4 || | 102 silk_assert( psEncC->indices.NLSFInterpCoef_Q2 == 4 || |
| 103 ( psEncC->useInterpolatedNLSFs && !psEncC->first_frame_after_reset && ps
EncC->nb_subfr == MAX_NB_SUBFR ) ); | 103 ( psEncC->useInterpolatedNLSFs && !psEncC->first_frame_after_reset && ps
EncC->nb_subfr == MAX_NB_SUBFR ) ); |
| 104 } | 104 } |
| OLD | NEW |