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 27 matching lines...) Expand all Loading... |
38 opus_int32 chirp_Q16 /* I Chirp factor (typica
lly in the range 0 to 1) */ | 38 opus_int32 chirp_Q16 /* I Chirp factor (typica
lly in the range 0 to 1) */ |
39 ) | 39 ) |
40 { | 40 { |
41 opus_int i; | 41 opus_int i; |
42 opus_int32 chirp_minus_one_Q16 = chirp_Q16 - 65536; | 42 opus_int32 chirp_minus_one_Q16 = chirp_Q16 - 65536; |
43 | 43 |
44 /* NB: Dont use silk_SMULWB, instead of silk_RSHIFT_ROUND( silk_MUL(), 16 ),
below. */ | 44 /* NB: Dont use silk_SMULWB, instead of silk_RSHIFT_ROUND( silk_MUL(), 16 ),
below. */ |
45 /* Bias in silk_SMULWB can lead to unstable filters
*/ | 45 /* Bias in silk_SMULWB can lead to unstable filters
*/ |
46 for( i = 0; i < d - 1; i++ ) { | 46 for( i = 0; i < d - 1; i++ ) { |
47 ar[ i ] = (opus_int16)silk_RSHIFT_ROUND( silk_MUL( chirp_Q16, ar[ i ]
), 16 ); | 47 ar[ i ] = (opus_int16)silk_RSHIFT_ROUND( silk_MUL( chirp_Q16, ar[ i ]
), 16 ); |
48 chirp_Q16 += silk_RSHIFT_ROUND( silk_MUL( chirp_Q16, chirp_mi
nus_one_Q16 ), 16 ); | 48 chirp_Q16 += silk_RSHIFT_ROUND( silk_MUL( chirp_Q16, chirp_m
inus_one_Q16 ), 16 ); |
49 } | 49 } |
50 ar[ d - 1 ] = (opus_int16)silk_RSHIFT_ROUND( silk_MUL( chirp_Q16, ar[ d - 1
] ), 16 ); | 50 ar[ d - 1 ] = (opus_int16)silk_RSHIFT_ROUND( silk_MUL( chirp_Q16, ar[ d - 1
] ), 16 ); |
51 } | 51 } |
OLD | NEW |