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 26 matching lines...) Expand all Loading... |
37 /* NB! State is kept internally and the */ | 37 /* NB! State is kept internally and the */ |
38 /* filter always starts with zero state */ | 38 /* filter always starts with zero state */ |
39 /* first d output samples are set to zero */ | 39 /* first d output samples are set to zero */ |
40 /*******************************************/ | 40 /*******************************************/ |
41 | 41 |
42 void silk_LPC_analysis_filter( | 42 void silk_LPC_analysis_filter( |
43 opus_int16 *out, /* O Output signal
*/ | 43 opus_int16 *out, /* O Output signal
*/ |
44 const opus_int16 *in, /* I Input signal
*/ | 44 const opus_int16 *in, /* I Input signal
*/ |
45 const opus_int16 *B, /* I MA prediction coeffi
cients, Q12 [order] */ | 45 const opus_int16 *B, /* I MA prediction coeffi
cients, Q12 [order] */ |
46 const opus_int32 len, /* I Signal length
*/ | 46 const opus_int32 len, /* I Signal length
*/ |
47 const opus_int32 d /* I Filter order
*/ | 47 const opus_int32 d, /* I Filter order
*/ |
| 48 int arch /* I Run-time architectur
e */ |
48 ) | 49 ) |
49 { | 50 { |
50 opus_int j; | 51 opus_int j; |
51 #ifdef FIXED_POINT | 52 #ifdef FIXED_POINT |
52 opus_int16 mem[SILK_MAX_ORDER_LPC]; | 53 opus_int16 mem[SILK_MAX_ORDER_LPC]; |
53 opus_int16 num[SILK_MAX_ORDER_LPC]; | 54 opus_int16 num[SILK_MAX_ORDER_LPC]; |
54 #else | 55 #else |
55 int ix; | 56 int ix; |
56 opus_int32 out32_Q12, out32; | 57 opus_int32 out32_Q12, out32; |
57 const opus_int16 *in_ptr; | 58 const opus_int16 *in_ptr; |
58 #endif | 59 #endif |
59 | 60 |
60 silk_assert( d >= 6 ); | 61 silk_assert( d >= 6 ); |
61 silk_assert( (d & 1) == 0 ); | 62 silk_assert( (d & 1) == 0 ); |
62 silk_assert( d <= len ); | 63 silk_assert( d <= len ); |
63 | 64 |
64 #ifdef FIXED_POINT | 65 #ifdef FIXED_POINT |
65 silk_assert( d <= SILK_MAX_ORDER_LPC ); | 66 silk_assert( d <= SILK_MAX_ORDER_LPC ); |
66 for ( j = 0; j < d; j++ ) { | 67 for ( j = 0; j < d; j++ ) { |
67 num[ j ] = -B[ j ]; | 68 num[ j ] = -B[ j ]; |
68 } | 69 } |
69 for (j=0;j<d;j++) { | 70 for (j=0;j<d;j++) { |
70 mem[ j ] = in[ d - j - 1 ]; | 71 mem[ j ] = in[ d - j - 1 ]; |
71 } | 72 } |
72 celt_fir( in + d, num, out + d, len - d, d, mem ); | 73 celt_fir( in + d, num, out + d, len - d, d, mem, arch ); |
73 for ( j = 0; j < d; j++ ) { | 74 for ( j = 0; j < d; j++ ) { |
74 out[ j ] = 0; | 75 out[ j ] = 0; |
75 } | 76 } |
76 #else | 77 #else |
| 78 (void)arch; |
77 for( ix = d; ix < len; ix++ ) { | 79 for( ix = d; ix < len; ix++ ) { |
78 in_ptr = &in[ ix - 1 ]; | 80 in_ptr = &in[ ix - 1 ]; |
79 | 81 |
80 out32_Q12 = silk_SMULBB( in_ptr[ 0 ], B[ 0 ] ); | 82 out32_Q12 = silk_SMULBB( in_ptr[ 0 ], B[ 0 ] ); |
81 /* Allowing wrap around so that two wraps can cancel each other. The rar
e | 83 /* Allowing wrap around so that two wraps can cancel each other. The rar
e |
82 cases where the result wraps around can only be triggered by invalid
streams*/ | 84 cases where the result wraps around can only be triggered by invalid
streams*/ |
83 out32_Q12 = silk_SMLABB_ovflw( out32_Q12, in_ptr[ -1 ], B[ 1 ] ); | 85 out32_Q12 = silk_SMLABB_ovflw( out32_Q12, in_ptr[ -1 ], B[ 1 ] ); |
84 out32_Q12 = silk_SMLABB_ovflw( out32_Q12, in_ptr[ -2 ], B[ 2 ] ); | 86 out32_Q12 = silk_SMLABB_ovflw( out32_Q12, in_ptr[ -2 ], B[ 2 ] ); |
85 out32_Q12 = silk_SMLABB_ovflw( out32_Q12, in_ptr[ -3 ], B[ 3 ] ); | 87 out32_Q12 = silk_SMLABB_ovflw( out32_Q12, in_ptr[ -3 ], B[ 3 ] ); |
86 out32_Q12 = silk_SMLABB_ovflw( out32_Q12, in_ptr[ -4 ], B[ 4 ] ); | 88 out32_Q12 = silk_SMLABB_ovflw( out32_Q12, in_ptr[ -4 ], B[ 4 ] ); |
(...skipping 10 matching lines...) Expand all Loading... |
97 out32 = silk_RSHIFT_ROUND( out32_Q12, 12 ); | 99 out32 = silk_RSHIFT_ROUND( out32_Q12, 12 ); |
98 | 100 |
99 /* Saturate output */ | 101 /* Saturate output */ |
100 out[ ix ] = (opus_int16)silk_SAT16( out32 ); | 102 out[ ix ] = (opus_int16)silk_SAT16( out32 ); |
101 } | 103 } |
102 | 104 |
103 /* Set first d output samples to zero */ | 105 /* Set first d output samples to zero */ |
104 silk_memset( out, 0, d * sizeof( opus_int16 ) ); | 106 silk_memset( out, 0, d * sizeof( opus_int16 ) ); |
105 #endif | 107 #endif |
106 } | 108 } |
OLD | NEW |