|
1 /*********************************************************************** |
|
2 Copyright (c) 2006-2011, Skype Limited. All rights reserved. |
|
3 Redistribution and use in source and binary forms, with or without |
|
4 modification, are permitted provided that the following conditions |
|
5 are met: |
|
6 - Redistributions of source code must retain the above copyright notice, |
|
7 this list of conditions and the following disclaimer. |
|
8 - Redistributions in binary form must reproduce the above copyright |
|
9 notice, this list of conditions and the following disclaimer in the |
|
10 documentation and/or other materials provided with the distribution. |
|
11 - Neither the name of Internet Society, IETF or IETF Trust, nor the |
|
12 names of specific contributors, may be used to endorse or promote |
|
13 products derived from this software without specific prior written |
|
14 permission. |
|
15 THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" |
|
16 AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE |
|
17 IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE |
|
18 ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE |
|
19 LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR |
|
20 CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF |
|
21 SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS |
|
22 INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN |
|
23 CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) |
|
24 ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE |
|
25 POSSIBILITY OF SUCH DAMAGE. |
|
26 ***********************************************************************/ |
|
27 |
|
28 #ifdef HAVE_CONFIG_H |
|
29 #include "config.h" |
|
30 #endif |
|
31 |
|
32 /* Filter coefficients for IIR/FIR polyphase resampling * |
|
33 * Total size: 179 Words (358 Bytes) */ |
|
34 |
|
35 #include "resampler_private.h" |
|
36 |
|
37 /* Matlab code for the notch filter coefficients: */ |
|
38 /* B = [1, 0.147, 1]; A = [1, 0.107, 0.89]; G = 0.93; freqz(G * B, A, 2^14, 16e3); axis([0, 8000, -10, 1]) */ |
|
39 /* fprintf('\t%6d, %6d, %6d, %6d\n', round(B(2)*2^16), round(-A(2)*2^16), round((1-A(3))*2^16), round(G*2^15)) */ |
|
40 /* const opus_int16 silk_resampler_up2_hq_notch[ 4 ] = { 9634, -7012, 7209, 30474 }; */ |
|
41 |
|
42 /* Tables with IIR and FIR coefficients for fractional downsamplers (123 Words) */ |
|
43 silk_DWORD_ALIGN const opus_int16 silk_Resampler_3_4_COEFS[ 2 + 3 * RESAMPLER_DOWN_ORDER_FIR0 / 2 ] = { |
|
44 -20694, -13867, |
|
45 -49, 64, 17, -157, 353, -496, 163, 11047, 22205, |
|
46 -39, 6, 91, -170, 186, 23, -896, 6336, 19928, |
|
47 -19, -36, 102, -89, -24, 328, -951, 2568, 15909, |
|
48 }; |
|
49 |
|
50 silk_DWORD_ALIGN const opus_int16 silk_Resampler_2_3_COEFS[ 2 + 2 * RESAMPLER_DOWN_ORDER_FIR0 / 2 ] = { |
|
51 -14457, -14019, |
|
52 64, 128, -122, 36, 310, -768, 584, 9267, 17733, |
|
53 12, 128, 18, -142, 288, -117, -865, 4123, 14459, |
|
54 }; |
|
55 |
|
56 silk_DWORD_ALIGN const opus_int16 silk_Resampler_1_2_COEFS[ 2 + RESAMPLER_DOWN_ORDER_FIR1 / 2 ] = { |
|
57 616, -14323, |
|
58 -10, 39, 58, -46, -84, 120, 184, -315, -541, 1284, 5380, 9024, |
|
59 }; |
|
60 |
|
61 silk_DWORD_ALIGN const opus_int16 silk_Resampler_1_3_COEFS[ 2 + RESAMPLER_DOWN_ORDER_FIR2 / 2 ] = { |
|
62 16102, -15162, |
|
63 -13, 0, 20, 26, 5, -31, -43, -4, 65, 90, 7, -157, -248, -44, 593, 1583, 2612, 3271, |
|
64 }; |
|
65 |
|
66 silk_DWORD_ALIGN const opus_int16 silk_Resampler_1_4_COEFS[ 2 + RESAMPLER_DOWN_ORDER_FIR2 / 2 ] = { |
|
67 22500, -15099, |
|
68 3, -14, -20, -15, 2, 25, 37, 25, -16, -71, -107, -79, 50, 292, 623, 982, 1288, 1464, |
|
69 }; |
|
70 |
|
71 silk_DWORD_ALIGN const opus_int16 silk_Resampler_1_6_COEFS[ 2 + RESAMPLER_DOWN_ORDER_FIR2 / 2 ] = { |
|
72 27540, -15257, |
|
73 17, 12, 8, 1, -10, -22, -30, -32, -22, 3, 44, 100, 168, 243, 317, 381, 429, 455, |
|
74 }; |
|
75 |
|
76 silk_DWORD_ALIGN const opus_int16 silk_Resampler_2_3_COEFS_LQ[ 2 + 2 * 2 ] = { |
|
77 -2797, -6507, |
|
78 4697, 10739, |
|
79 1567, 8276, |
|
80 }; |
|
81 |
|
82 /* Table with interplation fractions of 1/24, 3/24, 5/24, ... , 23/24 : 23/24 (46 Words) */ |
|
83 silk_DWORD_ALIGN const opus_int16 silk_resampler_frac_FIR_12[ 12 ][ RESAMPLER_ORDER_FIR_12 / 2 ] = { |
|
84 { 189, -600, 617, 30567 }, |
|
85 { 117, -159, -1070, 29704 }, |
|
86 { 52, 221, -2392, 28276 }, |
|
87 { -4, 529, -3350, 26341 }, |
|
88 { -48, 758, -3956, 23973 }, |
|
89 { -80, 905, -4235, 21254 }, |
|
90 { -99, 972, -4222, 18278 }, |
|
91 { -107, 967, -3957, 15143 }, |
|
92 { -103, 896, -3487, 11950 }, |
|
93 { -91, 773, -2865, 8798 }, |
|
94 { -71, 611, -2143, 5784 }, |
|
95 { -46, 425, -1375, 2996 }, |
|
96 }; |