annotate acelp_filters.c @ 10358:028a4deb56d3 libavcodec

Sun Raster RGB support
author pross
date Sun, 04 Oct 2009 08:42:23 +0000
parents cfb009484d2d
children 9f35b262d3f0
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
6772
695f4ab2db4c various filters for ACELP-based codecs
voroshil
parents:
diff changeset
1 /*
695f4ab2db4c various filters for ACELP-based codecs
voroshil
parents:
diff changeset
2 * various filters for ACELP-based codecs
695f4ab2db4c various filters for ACELP-based codecs
voroshil
parents:
diff changeset
3 *
695f4ab2db4c various filters for ACELP-based codecs
voroshil
parents:
diff changeset
4 * Copyright (c) 2008 Vladimir Voroshilov
695f4ab2db4c various filters for ACELP-based codecs
voroshil
parents:
diff changeset
5 *
695f4ab2db4c various filters for ACELP-based codecs
voroshil
parents:
diff changeset
6 * This file is part of FFmpeg.
695f4ab2db4c various filters for ACELP-based codecs
voroshil
parents:
diff changeset
7 *
695f4ab2db4c various filters for ACELP-based codecs
voroshil
parents:
diff changeset
8 * FFmpeg is free software; you can redistribute it and/or
695f4ab2db4c various filters for ACELP-based codecs
voroshil
parents:
diff changeset
9 * modify it under the terms of the GNU Lesser General Public
695f4ab2db4c various filters for ACELP-based codecs
voroshil
parents:
diff changeset
10 * License as published by the Free Software Foundation; either
695f4ab2db4c various filters for ACELP-based codecs
voroshil
parents:
diff changeset
11 * version 2.1 of the License, or (at your option) any later version.
695f4ab2db4c various filters for ACELP-based codecs
voroshil
parents:
diff changeset
12 *
695f4ab2db4c various filters for ACELP-based codecs
voroshil
parents:
diff changeset
13 * FFmpeg is distributed in the hope that it will be useful,
695f4ab2db4c various filters for ACELP-based codecs
voroshil
parents:
diff changeset
14 * but WITHOUT ANY WARRANTY; without even the implied warranty of
695f4ab2db4c various filters for ACELP-based codecs
voroshil
parents:
diff changeset
15 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
695f4ab2db4c various filters for ACELP-based codecs
voroshil
parents:
diff changeset
16 * Lesser General Public License for more details.
695f4ab2db4c various filters for ACELP-based codecs
voroshil
parents:
diff changeset
17 *
695f4ab2db4c various filters for ACELP-based codecs
voroshil
parents:
diff changeset
18 * You should have received a copy of the GNU Lesser General Public
695f4ab2db4c various filters for ACELP-based codecs
voroshil
parents:
diff changeset
19 * License along with FFmpeg; if not, write to the Free Software
695f4ab2db4c various filters for ACELP-based codecs
voroshil
parents:
diff changeset
20 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
695f4ab2db4c various filters for ACELP-based codecs
voroshil
parents:
diff changeset
21 */
695f4ab2db4c various filters for ACELP-based codecs
voroshil
parents:
diff changeset
22
695f4ab2db4c various filters for ACELP-based codecs
voroshil
parents:
diff changeset
23 #include <inttypes.h>
695f4ab2db4c various filters for ACELP-based codecs
voroshil
parents:
diff changeset
24
695f4ab2db4c various filters for ACELP-based codecs
voroshil
parents:
diff changeset
25 #include "avcodec.h"
695f4ab2db4c various filters for ACELP-based codecs
voroshil
parents:
diff changeset
26 #include "acelp_filters.h"
695f4ab2db4c various filters for ACELP-based codecs
voroshil
parents:
diff changeset
27
9792
05ea4942df9b cosmetics: K&R style reformatting
diego
parents: 8049
diff changeset
28 const int16_t ff_acelp_interp_filter[61] = { /* (0.15) */
6856
94465a2c3b34 Move pitch vector interpolation code to acelp_filters
voroshil
parents: 6775
diff changeset
29 29443, 28346, 25207, 20449, 14701, 8693,
94465a2c3b34 Move pitch vector interpolation code to acelp_filters
voroshil
parents: 6775
diff changeset
30 3143, -1352, -4402, -5865, -5850, -4673,
94465a2c3b34 Move pitch vector interpolation code to acelp_filters
voroshil
parents: 6775
diff changeset
31 -2783, -672, 1211, 2536, 3130, 2991,
94465a2c3b34 Move pitch vector interpolation code to acelp_filters
voroshil
parents: 6775
diff changeset
32 2259, 1170, 0, -1001, -1652, -1868,
94465a2c3b34 Move pitch vector interpolation code to acelp_filters
voroshil
parents: 6775
diff changeset
33 -1666, -1147, -464, 218, 756, 1060,
94465a2c3b34 Move pitch vector interpolation code to acelp_filters
voroshil
parents: 6775
diff changeset
34 1099, 904, 550, 135, -245, -514,
94465a2c3b34 Move pitch vector interpolation code to acelp_filters
voroshil
parents: 6775
diff changeset
35 -634, -602, -451, -231, 0, 191,
94465a2c3b34 Move pitch vector interpolation code to acelp_filters
voroshil
parents: 6775
diff changeset
36 308, 340, 296, 198, 78, -36,
94465a2c3b34 Move pitch vector interpolation code to acelp_filters
voroshil
parents: 6775
diff changeset
37 -120, -163, -165, -132, -79, -19,
94465a2c3b34 Move pitch vector interpolation code to acelp_filters
voroshil
parents: 6775
diff changeset
38 34, 73, 91, 89, 70, 38,
94465a2c3b34 Move pitch vector interpolation code to acelp_filters
voroshil
parents: 6775
diff changeset
39 0,
94465a2c3b34 Move pitch vector interpolation code to acelp_filters
voroshil
parents: 6775
diff changeset
40 };
94465a2c3b34 Move pitch vector interpolation code to acelp_filters
voroshil
parents: 6775
diff changeset
41
9792
05ea4942df9b cosmetics: K&R style reformatting
diego
parents: 8049
diff changeset
42 void ff_acelp_interpolate(int16_t* out, const int16_t* in,
05ea4942df9b cosmetics: K&R style reformatting
diego
parents: 8049
diff changeset
43 const int16_t* filter_coeffs, int precision,
05ea4942df9b cosmetics: K&R style reformatting
diego
parents: 8049
diff changeset
44 int frac_pos, int filter_length, int length)
6856
94465a2c3b34 Move pitch vector interpolation code to acelp_filters
voroshil
parents: 6775
diff changeset
45 {
94465a2c3b34 Move pitch vector interpolation code to acelp_filters
voroshil
parents: 6775
diff changeset
46 int n, i;
94465a2c3b34 Move pitch vector interpolation code to acelp_filters
voroshil
parents: 6775
diff changeset
47
94465a2c3b34 Move pitch vector interpolation code to acelp_filters
voroshil
parents: 6775
diff changeset
48 assert(pitch_delay_frac >= 0 && pitch_delay_frac < precision);
94465a2c3b34 Move pitch vector interpolation code to acelp_filters
voroshil
parents: 6775
diff changeset
49
9792
05ea4942df9b cosmetics: K&R style reformatting
diego
parents: 8049
diff changeset
50 for (n = 0; n < length; n++) {
6856
94465a2c3b34 Move pitch vector interpolation code to acelp_filters
voroshil
parents: 6775
diff changeset
51 int idx = 0;
94465a2c3b34 Move pitch vector interpolation code to acelp_filters
voroshil
parents: 6775
diff changeset
52 int v = 0x4000;
94465a2c3b34 Move pitch vector interpolation code to acelp_filters
voroshil
parents: 6775
diff changeset
53
9792
05ea4942df9b cosmetics: K&R style reformatting
diego
parents: 8049
diff changeset
54 for (i = 0; i < filter_length;) {
6856
94465a2c3b34 Move pitch vector interpolation code to acelp_filters
voroshil
parents: 6775
diff changeset
55
94465a2c3b34 Move pitch vector interpolation code to acelp_filters
voroshil
parents: 6775
diff changeset
56 /* The reference G.729 and AMR fixed point code performs clipping after
94465a2c3b34 Move pitch vector interpolation code to acelp_filters
voroshil
parents: 6775
diff changeset
57 each of the two following accumulations.
94465a2c3b34 Move pitch vector interpolation code to acelp_filters
voroshil
parents: 6775
diff changeset
58 Since clipping affects only the synthetic OVERFLOW test without
94465a2c3b34 Move pitch vector interpolation code to acelp_filters
voroshil
parents: 6775
diff changeset
59 causing an int type overflow, it was moved outside the loop. */
94465a2c3b34 Move pitch vector interpolation code to acelp_filters
voroshil
parents: 6775
diff changeset
60
94465a2c3b34 Move pitch vector interpolation code to acelp_filters
voroshil
parents: 6775
diff changeset
61 /* R(x):=ac_v[-k+x]
94465a2c3b34 Move pitch vector interpolation code to acelp_filters
voroshil
parents: 6775
diff changeset
62 v += R(n-i)*ff_acelp_interp_filter(t+6i)
94465a2c3b34 Move pitch vector interpolation code to acelp_filters
voroshil
parents: 6775
diff changeset
63 v += R(n+i+1)*ff_acelp_interp_filter(6-t+6i) */
94465a2c3b34 Move pitch vector interpolation code to acelp_filters
voroshil
parents: 6775
diff changeset
64
7649
8c5d7ede9c96 rename pitch_delay_frac in ff_acelp_interpolate()
michael
parents: 7641
diff changeset
65 v += in[n + i] * filter_coeffs[idx + frac_pos];
6856
94465a2c3b34 Move pitch vector interpolation code to acelp_filters
voroshil
parents: 6775
diff changeset
66 idx += precision;
94465a2c3b34 Move pitch vector interpolation code to acelp_filters
voroshil
parents: 6775
diff changeset
67 i++;
7649
8c5d7ede9c96 rename pitch_delay_frac in ff_acelp_interpolate()
michael
parents: 7641
diff changeset
68 v += in[n - i] * filter_coeffs[idx - frac_pos];
6856
94465a2c3b34 Move pitch vector interpolation code to acelp_filters
voroshil
parents: 6775
diff changeset
69 }
9792
05ea4942df9b cosmetics: K&R style reformatting
diego
parents: 8049
diff changeset
70 if (av_clip_int16(v >> 15) != (v >> 15))
7692
3ced4fb23342 Replace cliping in ff_acelp_interpolate() by a check&av_log, this should be the
michael
parents: 7653
diff changeset
71 av_log(NULL, AV_LOG_WARNING, "overflow that would need cliping in ff_acelp_interpolate()\n");
3ced4fb23342 Replace cliping in ff_acelp_interpolate() by a check&av_log, this should be the
michael
parents: 7653
diff changeset
72 out[n] = v >> 15;
6856
94465a2c3b34 Move pitch vector interpolation code to acelp_filters
voroshil
parents: 6775
diff changeset
73 }
94465a2c3b34 Move pitch vector interpolation code to acelp_filters
voroshil
parents: 6775
diff changeset
74 }
94465a2c3b34 Move pitch vector interpolation code to acelp_filters
voroshil
parents: 6775
diff changeset
75
6772
695f4ab2db4c various filters for ACELP-based codecs
voroshil
parents:
diff changeset
76
9792
05ea4942df9b cosmetics: K&R style reformatting
diego
parents: 8049
diff changeset
77 void ff_acelp_high_pass_filter(int16_t* out, int hpf_f[2],
05ea4942df9b cosmetics: K&R style reformatting
diego
parents: 8049
diff changeset
78 const int16_t* in, int length)
6772
695f4ab2db4c various filters for ACELP-based codecs
voroshil
parents:
diff changeset
79 {
695f4ab2db4c various filters for ACELP-based codecs
voroshil
parents:
diff changeset
80 int i;
695f4ab2db4c various filters for ACELP-based codecs
voroshil
parents:
diff changeset
81 int tmp;
695f4ab2db4c various filters for ACELP-based codecs
voroshil
parents:
diff changeset
82
9792
05ea4942df9b cosmetics: K&R style reformatting
diego
parents: 8049
diff changeset
83 for (i = 0; i < length; i++) {
05ea4942df9b cosmetics: K&R style reformatting
diego
parents: 8049
diff changeset
84 tmp = (hpf_f[0]* 15836LL) >> 13;
05ea4942df9b cosmetics: K&R style reformatting
diego
parents: 8049
diff changeset
85 tmp += (hpf_f[1]* -7667LL) >> 13;
7695
073847371624 (cosmetics) Remove wrong comments.
voroshil
parents: 7692
diff changeset
86 tmp += 7699 * (in[i] - 2*in[i-1] + in[i-2]);
6772
695f4ab2db4c various filters for ACELP-based codecs
voroshil
parents:
diff changeset
87
7696
76cf62255b28 (cosmetics) Describe for which tests clipping is required.
voroshil
parents: 7695
diff changeset
88 /* With "+0x800" rounding, clipping is needed
76cf62255b28 (cosmetics) Describe for which tests clipping is required.
voroshil
parents: 7695
diff changeset
89 for ALGTHM and SPEECH tests. */
7695
073847371624 (cosmetics) Remove wrong comments.
voroshil
parents: 7692
diff changeset
90 out[i] = av_clip_int16((tmp + 0x800) >> 12);
6772
695f4ab2db4c various filters for ACELP-based codecs
voroshil
parents:
diff changeset
91
695f4ab2db4c various filters for ACELP-based codecs
voroshil
parents:
diff changeset
92 hpf_f[1] = hpf_f[0];
695f4ab2db4c various filters for ACELP-based codecs
voroshil
parents:
diff changeset
93 hpf_f[0] = tmp;
695f4ab2db4c various filters for ACELP-based codecs
voroshil
parents:
diff changeset
94 }
695f4ab2db4c various filters for ACELP-based codecs
voroshil
parents:
diff changeset
95 }
10057
cfb009484d2d Add a function that can apply an order 2 rational transfer function in-place.
superdump
parents: 9792
diff changeset
96
cfb009484d2d Add a function that can apply an order 2 rational transfer function in-place.
superdump
parents: 9792
diff changeset
97 void ff_acelp_apply_order_2_transfer_function(float *buf,
cfb009484d2d Add a function that can apply an order 2 rational transfer function in-place.
superdump
parents: 9792
diff changeset
98 const float zero_coeffs[2],
cfb009484d2d Add a function that can apply an order 2 rational transfer function in-place.
superdump
parents: 9792
diff changeset
99 const float pole_coeffs[2],
cfb009484d2d Add a function that can apply an order 2 rational transfer function in-place.
superdump
parents: 9792
diff changeset
100 float gain, float mem[2], int n)
cfb009484d2d Add a function that can apply an order 2 rational transfer function in-place.
superdump
parents: 9792
diff changeset
101 {
cfb009484d2d Add a function that can apply an order 2 rational transfer function in-place.
superdump
parents: 9792
diff changeset
102 int i;
cfb009484d2d Add a function that can apply an order 2 rational transfer function in-place.
superdump
parents: 9792
diff changeset
103 float tmp;
cfb009484d2d Add a function that can apply an order 2 rational transfer function in-place.
superdump
parents: 9792
diff changeset
104
cfb009484d2d Add a function that can apply an order 2 rational transfer function in-place.
superdump
parents: 9792
diff changeset
105 for (i = 0; i < n; i++) {
cfb009484d2d Add a function that can apply an order 2 rational transfer function in-place.
superdump
parents: 9792
diff changeset
106 tmp = gain * buf[i] - pole_coeffs[0] * mem[0] - pole_coeffs[1] * mem[1];
cfb009484d2d Add a function that can apply an order 2 rational transfer function in-place.
superdump
parents: 9792
diff changeset
107 buf[i] = tmp + zero_coeffs[0] * mem[0] + zero_coeffs[1] * mem[1];
cfb009484d2d Add a function that can apply an order 2 rational transfer function in-place.
superdump
parents: 9792
diff changeset
108
cfb009484d2d Add a function that can apply an order 2 rational transfer function in-place.
superdump
parents: 9792
diff changeset
109 mem[1] = mem[0];
cfb009484d2d Add a function that can apply an order 2 rational transfer function in-place.
superdump
parents: 9792
diff changeset
110 mem[0] = tmp;
cfb009484d2d Add a function that can apply an order 2 rational transfer function in-place.
superdump
parents: 9792
diff changeset
111 }
cfb009484d2d Add a function that can apply an order 2 rational transfer function in-place.
superdump
parents: 9792
diff changeset
112 }