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 29 matching lines...) Expand all Loading... |
40 #define MIN_RSHIFTS -16 | 40 #define MIN_RSHIFTS -16 |
41 #define MAX_RSHIFTS (32 - QA) | 41 #define MAX_RSHIFTS (32 - QA) |
42 | 42 |
43 /* Compute reflection coefficients from input signal */ | 43 /* Compute reflection coefficients from input signal */ |
44 void silk_burg_modified( | 44 void silk_burg_modified( |
45 opus_int32 *res_nrg, /* O Residual energy
*/ | 45 opus_int32 *res_nrg, /* O Residual energy
*/ |
46 opus_int *res_nrg_Q, /* O Residual energy Q va
lue */ | 46 opus_int *res_nrg_Q, /* O Residual energy Q va
lue */ |
47 opus_int32 A_Q16[], /* O Prediction coefficie
nts (length order) */ | 47 opus_int32 A_Q16[], /* O Prediction coefficie
nts (length order) */ |
48 const opus_int16 x[], /* I Input signal, length
: nb_subfr * ( D + subfr_length ) */ | 48 const opus_int16 x[], /* I Input signal, length
: nb_subfr * ( D + subfr_length ) */ |
49 const opus_int32 minInvGain_Q30, /* I Inverse of max predi
ction gain */ | 49 const opus_int32 minInvGain_Q30, /* I Inverse of max predi
ction gain */ |
50 const opus_int subfr_length, /* I Input signal subfram
e length (incl. D preceeding samples) */ | 50 const opus_int subfr_length, /* I Input signal subfram
e length (incl. D preceding samples) */ |
51 const opus_int nb_subfr, /* I Number of subframes
stacked in x */ | 51 const opus_int nb_subfr, /* I Number of subframes
stacked in x */ |
52 const opus_int D /* I Order
*/ | 52 const opus_int D /* I Order
*/ |
53 ) | 53 ) |
54 { | 54 { |
55 opus_int k, n, s, lz, rshifts, rshifts_extra, reached_max_gain; | 55 opus_int k, n, s, lz, rshifts, rshifts_extra, reached_max_gain; |
56 opus_int32 C0, num, nrg, rc_Q31, invGain_Q30, Atmp_QA, Atmp1, tmp1, tm
p2, x1, x2; | 56 opus_int32 C0, num, nrg, rc_Q31, invGain_Q30, Atmp_QA, Atmp1, tmp1, tm
p2, x1, x2; |
57 const opus_int16 *x_ptr; | 57 const opus_int16 *x_ptr; |
58 opus_int32 C_first_row[ SILK_MAX_ORDER_LPC ]; | 58 opus_int32 C_first_row[ SILK_MAX_ORDER_LPC ]; |
59 opus_int32 C_last_row[ SILK_MAX_ORDER_LPC ]; | 59 opus_int32 C_last_row[ SILK_MAX_ORDER_LPC ]; |
60 opus_int32 Af_QA[ SILK_MAX_ORDER_LPC ]; | 60 opus_int32 Af_QA[ SILK_MAX_ORDER_LPC ]; |
(...skipping 170 matching lines...) Expand 10 before | Expand all | Expand 10 after Loading... |
231 CAf[ k ] = silk_ADD_LSHIFT32( tmp1, silk_SMMUL( tmp2, rc_Q31
), 1 ); /* Q( -rshifts ) */ | 231 CAf[ k ] = silk_ADD_LSHIFT32( tmp1, silk_SMMUL( tmp2, rc_Q31
), 1 ); /* Q( -rshifts ) */ |
232 CAb[ n - k + 1 ] = silk_ADD_LSHIFT32( tmp2, silk_SMMUL( tmp1, rc_Q31
), 1 ); /* Q( -rshifts ) */ | 232 CAb[ n - k + 1 ] = silk_ADD_LSHIFT32( tmp2, silk_SMMUL( tmp1, rc_Q31
), 1 ); /* Q( -rshifts ) */ |
233 } | 233 } |
234 } | 234 } |
235 | 235 |
236 if( reached_max_gain ) { | 236 if( reached_max_gain ) { |
237 for( k = 0; k < D; k++ ) { | 237 for( k = 0; k < D; k++ ) { |
238 /* Scale coefficients */ | 238 /* Scale coefficients */ |
239 A_Q16[ k ] = -silk_RSHIFT_ROUND( Af_QA[ k ], QA - 16 ); | 239 A_Q16[ k ] = -silk_RSHIFT_ROUND( Af_QA[ k ], QA - 16 ); |
240 } | 240 } |
241 /* Subtract energy of preceeding samples from C0 */ | 241 /* Subtract energy of preceding samples from C0 */ |
242 if( rshifts > 0 ) { | 242 if( rshifts > 0 ) { |
243 for( s = 0; s < nb_subfr; s++ ) { | 243 for( s = 0; s < nb_subfr; s++ ) { |
244 x_ptr = x + s * subfr_length; | 244 x_ptr = x + s * subfr_length; |
245 C0 -= (opus_int32)silk_RSHIFT64( silk_inner_prod16_aligned_64( x
_ptr, x_ptr, D ), rshifts ); | 245 C0 -= (opus_int32)silk_RSHIFT64( silk_inner_prod16_aligned_64( x
_ptr, x_ptr, D ), rshifts ); |
246 } | 246 } |
247 } else { | 247 } else { |
248 for( s = 0; s < nb_subfr; s++ ) { | 248 for( s = 0; s < nb_subfr; s++ ) { |
249 x_ptr = x + s * subfr_length; | 249 x_ptr = x + s * subfr_length; |
250 C0 -= silk_LSHIFT32( silk_inner_prod_aligned( x_ptr, x_ptr, D ),
-rshifts ); | 250 C0 -= silk_LSHIFT32( silk_inner_prod_aligned( x_ptr, x_ptr, D ),
-rshifts ); |
251 } | 251 } |
252 } | 252 } |
253 /* Approximate residual energy */ | 253 /* Approximate residual energy */ |
254 *res_nrg = silk_LSHIFT( silk_SMMUL( invGain_Q30, C0 ), 2 ); | 254 *res_nrg = silk_LSHIFT( silk_SMMUL( invGain_Q30, C0 ), 2 ); |
255 *res_nrg_Q = -rshifts; | 255 *res_nrg_Q = -rshifts; |
256 } else { | 256 } else { |
257 /* Return residual energy */ | 257 /* Return residual energy */ |
258 nrg = CAf[ 0 ];
/* Q( -rshifts ) */ | 258 nrg = CAf[ 0 ];
/* Q( -rshifts ) */ |
259 tmp1 = (opus_int32)1 << 16;
/* Q16 */ | 259 tmp1 = (opus_int32)1 << 16;
/* Q16 */ |
260 for( k = 0; k < D; k++ ) { | 260 for( k = 0; k < D; k++ ) { |
261 Atmp1 = silk_RSHIFT_ROUND( Af_QA[ k ], QA - 16 );
/* Q16 */ | 261 Atmp1 = silk_RSHIFT_ROUND( Af_QA[ k ], QA - 16 );
/* Q16 */ |
262 nrg = silk_SMLAWW( nrg, CAf[ k + 1 ], Atmp1 );
/* Q( -rshifts ) */ | 262 nrg = silk_SMLAWW( nrg, CAf[ k + 1 ], Atmp1 );
/* Q( -rshifts ) */ |
263 tmp1 = silk_SMLAWW( tmp1, Atmp1, Atmp1 );
/* Q16 */ | 263 tmp1 = silk_SMLAWW( tmp1, Atmp1, Atmp1 );
/* Q16 */ |
264 A_Q16[ k ] = -Atmp1; | 264 A_Q16[ k ] = -Atmp1; |
265 } | 265 } |
266 *res_nrg = silk_SMLAWW( nrg, silk_SMMUL( FIND_LPC_COND_FAC, C0 ), -tmp1
); /* Q( -rshifts ) */ | 266 *res_nrg = silk_SMLAWW( nrg, silk_SMMUL( FIND_LPC_COND_FAC, C0 ), -tmp1
); /* Q( -rshifts ) */ |
267 *res_nrg_Q = -rshifts; | 267 *res_nrg_Q = -rshifts; |
268 } | 268 } |
269 } | 269 } |
OLD | NEW |