annotate src/opus-1.3/silk/resampler_rom.c @ 70:9e21af8f0420

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