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 23 matching lines...) Expand all Loading... |
34 /* Find LPC and LTP coefficients */ | 34 /* Find LPC and LTP coefficients */ |
35 void silk_find_pred_coefs_FLP( | 35 void silk_find_pred_coefs_FLP( |
36 silk_encoder_state_FLP *psEnc, /* I/O
Encoder state FLP */ | 36 silk_encoder_state_FLP *psEnc, /* I/O
Encoder state FLP */ |
37 silk_encoder_control_FLP *psEncCtrl, /* I/O
Encoder control FLP */ | 37 silk_encoder_control_FLP *psEncCtrl, /* I/O
Encoder control FLP */ |
38 const silk_float res_pitch[], /* I
Residual from pitch analysis */ | 38 const silk_float res_pitch[], /* I
Residual from pitch analysis */ |
39 const silk_float x[], /* I
Speech signal */ | 39 const silk_float x[], /* I
Speech signal */ |
40 opus_int condCoding /* I
The type of conditional coding to use */ | 40 opus_int condCoding /* I
The type of conditional coding to use */ |
41 ) | 41 ) |
42 { | 42 { |
43 opus_int i; | 43 opus_int i; |
44 silk_float WLTP[ MAX_NB_SUBFR * LTP_ORDER * LTP_ORDER ]; | 44 silk_float XXLTP[ MAX_NB_SUBFR * LTP_ORDER * LTP_ORDER ]; |
45 silk_float invGains[ MAX_NB_SUBFR ], Wght[ MAX_NB_SUBFR ]; | 45 silk_float xXLTP[ MAX_NB_SUBFR * LTP_ORDER ]; |
| 46 silk_float invGains[ MAX_NB_SUBFR ]; |
46 opus_int16 NLSF_Q15[ MAX_LPC_ORDER ]; | 47 opus_int16 NLSF_Q15[ MAX_LPC_ORDER ]; |
47 const silk_float *x_ptr; | 48 const silk_float *x_ptr; |
48 silk_float *x_pre_ptr, LPC_in_pre[ MAX_NB_SUBFR * MAX_LPC_ORDER + MAX_
FRAME_LENGTH ]; | 49 silk_float *x_pre_ptr, LPC_in_pre[ MAX_NB_SUBFR * MAX_LPC_ORDER + MAX_
FRAME_LENGTH ]; |
49 silk_float minInvGain; | 50 silk_float minInvGain; |
50 | 51 |
51 /* Weighting for weighted least squares */ | 52 /* Weighting for weighted least squares */ |
52 for( i = 0; i < psEnc->sCmn.nb_subfr; i++ ) { | 53 for( i = 0; i < psEnc->sCmn.nb_subfr; i++ ) { |
53 silk_assert( psEncCtrl->Gains[ i ] > 0.0f ); | 54 silk_assert( psEncCtrl->Gains[ i ] > 0.0f ); |
54 invGains[ i ] = 1.0f / psEncCtrl->Gains[ i ]; | 55 invGains[ i ] = 1.0f / psEncCtrl->Gains[ i ]; |
55 Wght[ i ] = invGains[ i ] * invGains[ i ]; | |
56 } | 56 } |
57 | 57 |
58 if( psEnc->sCmn.indices.signalType == TYPE_VOICED ) { | 58 if( psEnc->sCmn.indices.signalType == TYPE_VOICED ) { |
59 /**********/ | 59 /**********/ |
60 /* VOICED */ | 60 /* VOICED */ |
61 /**********/ | 61 /**********/ |
62 silk_assert( psEnc->sCmn.ltp_mem_length - psEnc->sCmn.predictLPCOrder >=
psEncCtrl->pitchL[ 0 ] + LTP_ORDER / 2 ); | 62 silk_assert( psEnc->sCmn.ltp_mem_length - psEnc->sCmn.predictLPCOrder >=
psEncCtrl->pitchL[ 0 ] + LTP_ORDER / 2 ); |
63 | 63 |
64 /* LTP analysis */ | 64 /* LTP analysis */ |
65 silk_find_LTP_FLP( psEncCtrl->LTPCoef, WLTP, &psEncCtrl->LTPredCodGain,
res_pitch, | 65 silk_find_LTP_FLP( XXLTP, xXLTP, res_pitch, psEncCtrl->pitchL, psEnc->sC
mn.subfr_length, psEnc->sCmn.nb_subfr ); |
66 psEncCtrl->pitchL, Wght, psEnc->sCmn.subfr_length, psEnc->sCmn.nb_su
bfr, psEnc->sCmn.ltp_mem_length ); | |
67 | 66 |
68 /* Quantize LTP gain parameters */ | 67 /* Quantize LTP gain parameters */ |
69 silk_quant_LTP_gains_FLP( psEncCtrl->LTPCoef, psEnc->sCmn.indices.LTPInd
ex, &psEnc->sCmn.indices.PERIndex, | 68 silk_quant_LTP_gains_FLP( psEncCtrl->LTPCoef, psEnc->sCmn.indices.LTPInd
ex, &psEnc->sCmn.indices.PERIndex, |
70 &psEnc->sCmn.sum_log_gain_Q7, WLTP, psEnc->sCmn.mu_LTP_Q9, psEnc->sC
mn.LTPQuantLowComplexity, psEnc->sCmn.nb_subfr, | 69 &psEnc->sCmn.sum_log_gain_Q7, &psEncCtrl->LTPredCodGain, XXLTP, xXLT
P, psEnc->sCmn.subfr_length, psEnc->sCmn.nb_subfr, psEnc->sCmn.arch ); |
71 psEnc->sCmn.arch ); | |
72 | 70 |
73 /* Control LTP scaling */ | 71 /* Control LTP scaling */ |
74 silk_LTP_scale_ctrl_FLP( psEnc, psEncCtrl, condCoding ); | 72 silk_LTP_scale_ctrl_FLP( psEnc, psEncCtrl, condCoding ); |
75 | 73 |
76 /* Create LTP residual */ | 74 /* Create LTP residual */ |
77 silk_LTP_analysis_filter_FLP( LPC_in_pre, x - psEnc->sCmn.predictLPCOrde
r, psEncCtrl->LTPCoef, | 75 silk_LTP_analysis_filter_FLP( LPC_in_pre, x - psEnc->sCmn.predictLPCOrde
r, psEncCtrl->LTPCoef, |
78 psEncCtrl->pitchL, invGains, psEnc->sCmn.subfr_length, psEnc->sCmn.n
b_subfr, psEnc->sCmn.predictLPCOrder ); | 76 psEncCtrl->pitchL, invGains, psEnc->sCmn.subfr_length, psEnc->sCmn.n
b_subfr, psEnc->sCmn.predictLPCOrder ); |
79 } else { | 77 } else { |
80 /************/ | 78 /************/ |
81 /* UNVOICED */ | 79 /* UNVOICED */ |
(...skipping 27 matching lines...) Expand all Loading... |
109 silk_process_NLSFs_FLP( &psEnc->sCmn, psEncCtrl->PredCoef, NLSF_Q15, psEnc->
sCmn.prev_NLSFq_Q15 ); | 107 silk_process_NLSFs_FLP( &psEnc->sCmn, psEncCtrl->PredCoef, NLSF_Q15, psEnc->
sCmn.prev_NLSFq_Q15 ); |
110 | 108 |
111 /* Calculate residual energy using quantized LPC coefficients */ | 109 /* Calculate residual energy using quantized LPC coefficients */ |
112 silk_residual_energy_FLP( psEncCtrl->ResNrg, LPC_in_pre, psEncCtrl->PredCoef
, psEncCtrl->Gains, | 110 silk_residual_energy_FLP( psEncCtrl->ResNrg, LPC_in_pre, psEncCtrl->PredCoef
, psEncCtrl->Gains, |
113 psEnc->sCmn.subfr_length, psEnc->sCmn.nb_subfr, psEnc->sCmn.predictLPCOr
der ); | 111 psEnc->sCmn.subfr_length, psEnc->sCmn.nb_subfr, psEnc->sCmn.predictLPCOr
der ); |
114 | 112 |
115 /* Copy to prediction struct for use in next frame for interpolation */ | 113 /* Copy to prediction struct for use in next frame for interpolation */ |
116 silk_memcpy( psEnc->sCmn.prev_NLSFq_Q15, NLSF_Q15, sizeof( psEnc->sCmn.prev_
NLSFq_Q15 ) ); | 114 silk_memcpy( psEnc->sCmn.prev_NLSFq_Q15, NLSF_Q15, sizeof( psEnc->sCmn.prev_
NLSFq_Q15 ) ); |
117 } | 115 } |
118 | 116 |
OLD | NEW |