annotate ra288.c @ 7434:6a029d8dcfa8 libavcodec

Cosmetics: remove braces and useless newline
author vitor
date Mon, 28 Jul 2008 04:07:18 +0000
parents 21ecdb66d5b0
children 014780c09494
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
1304
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
1 /*
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
2 * RealAudio 2.0 (28.8K)
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
3 * Copyright (c) 2003 the ffmpeg project
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
4 *
3947
c8c591fe26f8 Change license headers to say 'FFmpeg' instead of 'this program/this library'
diego
parents: 3036
diff changeset
5 * This file is part of FFmpeg.
c8c591fe26f8 Change license headers to say 'FFmpeg' instead of 'this program/this library'
diego
parents: 3036
diff changeset
6 *
c8c591fe26f8 Change license headers to say 'FFmpeg' instead of 'this program/this library'
diego
parents: 3036
diff changeset
7 * FFmpeg is free software; you can redistribute it and/or
1304
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
8 * modify it under the terms of the GNU Lesser General Public
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
9 * License as published by the Free Software Foundation; either
3947
c8c591fe26f8 Change license headers to say 'FFmpeg' instead of 'this program/this library'
diego
parents: 3036
diff changeset
10 * version 2.1 of the License, or (at your option) any later version.
1304
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
11 *
3947
c8c591fe26f8 Change license headers to say 'FFmpeg' instead of 'this program/this library'
diego
parents: 3036
diff changeset
12 * FFmpeg is distributed in the hope that it will be useful,
1304
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
13 * but WITHOUT ANY WARRANTY; without even the implied warranty of
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
14 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
15 * Lesser General Public License for more details.
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
16 *
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
17 * You should have received a copy of the GNU Lesser General Public
3947
c8c591fe26f8 Change license headers to say 'FFmpeg' instead of 'this program/this library'
diego
parents: 3036
diff changeset
18 * License along with FFmpeg; if not, write to the Free Software
3036
0b546eab515d Update licensing information: The FSF changed postal address.
diego
parents: 2979
diff changeset
19 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
1304
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
20 */
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
21
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
22 #include "avcodec.h"
7171
9176fe6f17c9 Rewrite unpack() using the bitstream reader
vitor
parents: 7169
diff changeset
23 #define ALT_BITSTREAM_READER_LE
9176fe6f17c9 Rewrite unpack() using the bitstream reader
vitor
parents: 7169
diff changeset
24 #include "bitstream.h"
1335
b4a72edb3a71 moved the tables into header files (and applied the 'static' patch). Nick: why do you like mergeing tables and code into one file, so making it unusable big?
al3x
parents: 1316
diff changeset
25 #include "ra288.h"
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2956
diff changeset
26
1304
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
27 typedef struct {
7428
e8eca33561c7 Rename two context buffers: pr{1,2} -> {sp,gain}_lpc. Also add a few doxy
vitor
parents: 7426
diff changeset
28 float sp_lpc[36]; ///< LPC coefficients for speech data (spec: A)
e8eca33561c7 Rename two context buffers: pr{1,2} -> {sp,gain}_lpc. Also add a few doxy
vitor
parents: 7426
diff changeset
29 float gain_lpc[10]; ///< LPC coefficients for gain (spec: GB)
7213
cf82f8a56940 Remove phasep context var, it is just phase*5
vitor
parents: 7212
diff changeset
30 int phase;
1304
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
31
7418
83d005a098b5 Rename st{1,2}{a,b} context variables to more meaningful names: {sp,gain}_{rec,hist}
vitor
parents: 7416
diff changeset
32 float sp_hist[111]; ///< Speech data history (spec: SB)
83d005a098b5 Rename st{1,2}{a,b} context variables to more meaningful names: {sp,gain}_{rec,hist}
vitor
parents: 7416
diff changeset
33
83d005a098b5 Rename st{1,2}{a,b} context variables to more meaningful names: {sp,gain}_{rec,hist}
vitor
parents: 7416
diff changeset
34 /** Speech part of the gain autocorrelation (spec: REXP) */
83d005a098b5 Rename st{1,2}{a,b} context variables to more meaningful names: {sp,gain}_{rec,hist}
vitor
parents: 7416
diff changeset
35 float sp_rec[37];
83d005a098b5 Rename st{1,2}{a,b} context variables to more meaningful names: {sp,gain}_{rec,hist}
vitor
parents: 7416
diff changeset
36
83d005a098b5 Rename st{1,2}{a,b} context variables to more meaningful names: {sp,gain}_{rec,hist}
vitor
parents: 7416
diff changeset
37 float gain_hist[38]; ///< Log-gain history (spec: SBLG)
83d005a098b5 Rename st{1,2}{a,b} context variables to more meaningful names: {sp,gain}_{rec,hist}
vitor
parents: 7416
diff changeset
38
83d005a098b5 Rename st{1,2}{a,b} context variables to more meaningful names: {sp,gain}_{rec,hist}
vitor
parents: 7416
diff changeset
39 /** Recursive part of the gain autocorrelation (spec: REXPLG) */
83d005a098b5 Rename st{1,2}{a,b} context variables to more meaningful names: {sp,gain}_{rec,hist}
vitor
parents: 7416
diff changeset
40 float gain_rec[11];
83d005a098b5 Rename st{1,2}{a,b} context variables to more meaningful names: {sp,gain}_{rec,hist}
vitor
parents: 7416
diff changeset
41
7167
9b93a3cf6440 Reindent the whole file
vitor
parents: 7040
diff changeset
42 float sb[41];
9b93a3cf6440 Reindent the whole file
vitor
parents: 7040
diff changeset
43 float lhist[10];
7411
d41948bf80a1 Give the context structure of ra288 a better name: s/Real288_internal/RA288Context/
vitor
parents: 7375
diff changeset
44 } RA288Context;
1304
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
45
7294
bf3952d5e76c Declare function parameters as const
vitor
parents: 7275
diff changeset
46 static inline float scalar_product_float(const float * v1, const float * v2,
bf3952d5e76c Declare function parameters as const
vitor
parents: 7275
diff changeset
47 int size)
7201
4d441356fa23 Create a function for float dot product instead of duplicating all over the file
vitor
parents: 7195
diff changeset
48 {
4d441356fa23 Create a function for float dot product instead of duplicating all over the file
vitor
parents: 7195
diff changeset
49 float res = 0.;
4d441356fa23 Create a function for float dot product instead of duplicating all over the file
vitor
parents: 7195
diff changeset
50
4d441356fa23 Create a function for float dot product instead of duplicating all over the file
vitor
parents: 7195
diff changeset
51 while (size--)
4d441356fa23 Create a function for float dot product instead of duplicating all over the file
vitor
parents: 7195
diff changeset
52 res += *v1++ * *v2++;
4d441356fa23 Create a function for float dot product instead of duplicating all over the file
vitor
parents: 7195
diff changeset
53
4d441356fa23 Create a function for float dot product instead of duplicating all over the file
vitor
parents: 7195
diff changeset
54 return res;
4d441356fa23 Create a function for float dot product instead of duplicating all over the file
vitor
parents: 7195
diff changeset
55 }
4d441356fa23 Create a function for float dot product instead of duplicating all over the file
vitor
parents: 7195
diff changeset
56
7415
649496faa1d2 Move colmult() function to the beginning of file to group DSP-related functions.
vitor
parents: 7414
diff changeset
57 static void colmult(float *tgt, const float *m1, const float *m2, int n)
649496faa1d2 Move colmult() function to the beginning of file to group DSP-related functions.
vitor
parents: 7414
diff changeset
58 {
649496faa1d2 Move colmult() function to the beginning of file to group DSP-related functions.
vitor
parents: 7414
diff changeset
59 while (n--)
7416
a65363f8ab4b Cosmetics: remove useless parenthesis
vitor
parents: 7415
diff changeset
60 *tgt++ = *m1++ * *m2++;
7415
649496faa1d2 Move colmult() function to the beginning of file to group DSP-related functions.
vitor
parents: 7414
diff changeset
61 }
649496faa1d2 Move colmult() function to the beginning of file to group DSP-related functions.
vitor
parents: 7414
diff changeset
62
1304
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
63 /* Decode and produce output */
7412
51da3d8200d0 Give the context variables of ra288 a better name: s/glob/ractx/
vitor
parents: 7411
diff changeset
64 static void decode(RA288Context *ractx, float gain, int cb_coef)
1304
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
65 {
7195
5e5c1178ae9a Write for loops in a more standard way. In my opinion, they are much more readable now.
vitor
parents: 7193
diff changeset
66 int x, y;
7268
7e51fc7805ac Do not declare as double a var that only stores a float
vitor
parents: 7267
diff changeset
67 double sumsum;
7e51fc7805ac Do not declare as double a var that only stores a float
vitor
parents: 7267
diff changeset
68 float sum, buffer[5];
7167
9b93a3cf6440 Reindent the whole file
vitor
parents: 7040
diff changeset
69
7412
51da3d8200d0 Give the context variables of ra288 a better name: s/glob/ractx/
vitor
parents: 7411
diff changeset
70 memmove(ractx->sb + 5, ractx->sb, 36 * sizeof(*ractx->sb));
7167
9b93a3cf6440 Reindent the whole file
vitor
parents: 7040
diff changeset
71
7202
2a4ac127112c Cosmetics: remove braces
vitor
parents: 7201
diff changeset
72 for (x=4; x >= 0; x--)
7428
e8eca33561c7 Rename two context buffers: pr{1,2} -> {sp,gain}_lpc. Also add a few doxy
vitor
parents: 7426
diff changeset
73 ractx->sb[x] = -scalar_product_float(ractx->sb + x + 1,
e8eca33561c7 Rename two context buffers: pr{1,2} -> {sp,gain}_lpc. Also add a few doxy
vitor
parents: 7426
diff changeset
74 ractx->sp_lpc, 36);
7167
9b93a3cf6440 Reindent the whole file
vitor
parents: 7040
diff changeset
75
7431
da24b77f46f6 Add a few comments pointing to the G.728 specification
vitor
parents: 7430
diff changeset
76 /* block 46 of G.728 spec */
7428
e8eca33561c7 Rename two context buffers: pr{1,2} -> {sp,gain}_lpc. Also add a few doxy
vitor
parents: 7426
diff changeset
77 sum = 32. - scalar_product_float(ractx->gain_lpc, ractx->lhist, 10);
7167
9b93a3cf6440 Reindent the whole file
vitor
parents: 7040
diff changeset
78
7431
da24b77f46f6 Add a few comments pointing to the G.728 specification
vitor
parents: 7430
diff changeset
79 /* block 47 of G.728 spec */
7269
04ecb6594999 Another av_clipf() usage
vitor
parents: 7268
diff changeset
80 sum = av_clipf(sum, 0, 60);
1304
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
81
7431
da24b77f46f6 Add a few comments pointing to the G.728 specification
vitor
parents: 7430
diff changeset
82 /* block 48 of G.728 spec */
7192
9517beb3b491 Simplify
vitor
parents: 7191
diff changeset
83 sumsum = exp(sum * 0.1151292546497) * gain; /* pow(10.0,sum/20)*f */
7167
9b93a3cf6440 Reindent the whole file
vitor
parents: 7040
diff changeset
84
7210
43048f7a1d7c Cosmetics: remove braces
vitor
parents: 7209
diff changeset
85 for (x=0; x < 5; x++)
7191
98c5f3839e9b Remove unneeded var
vitor
parents: 7190
diff changeset
86 buffer[x] = codetable[cb_coef][x] * sumsum;
1304
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
87
7275
c5bcb9e4c5b2 Revert r14218. FFMAX is a macro, so
vitor
parents: 7270
diff changeset
88 sum = scalar_product_float(buffer, buffer, 5) / 5;
c5bcb9e4c5b2 Revert r14218. FFMAX is a macro, so
vitor
parents: 7270
diff changeset
89
c5bcb9e4c5b2 Revert r14218. FFMAX is a macro, so
vitor
parents: 7270
diff changeset
90 sum = FFMAX(sum, 1);
7167
9b93a3cf6440 Reindent the whole file
vitor
parents: 7040
diff changeset
91
9b93a3cf6440 Reindent the whole file
vitor
parents: 7040
diff changeset
92 /* shift and store */
7412
51da3d8200d0 Give the context variables of ra288 a better name: s/glob/ractx/
vitor
parents: 7411
diff changeset
93 memmove(ractx->lhist, ractx->lhist - 1, 10 * sizeof(*ractx->lhist));
7167
9b93a3cf6440 Reindent the whole file
vitor
parents: 7040
diff changeset
94
7430
4ebfb53ff2df Remove the history buffer from the context. It can easily be evaluated
vitor
parents: 7428
diff changeset
95 *ractx->lhist = 10 * log10(sum) - 32;
1304
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
96
7167
9b93a3cf6440 Reindent the whole file
vitor
parents: 7040
diff changeset
97 for (x=1; x < 5; x++)
7195
5e5c1178ae9a Write for loops in a more standard way. In my opinion, they are much more readable now.
vitor
parents: 7193
diff changeset
98 for (y=x-1; y >= 0; y--)
7428
e8eca33561c7 Rename two context buffers: pr{1,2} -> {sp,gain}_lpc. Also add a few doxy
vitor
parents: 7426
diff changeset
99 buffer[x] -= ractx->sp_lpc[x-y-1] * buffer[y];
1304
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
100
7167
9b93a3cf6440 Reindent the whole file
vitor
parents: 7040
diff changeset
101 /* output */
7434
6a029d8dcfa8 Cosmetics: remove braces and useless newline
vitor
parents: 7433
diff changeset
102 for (x=0; x < 5; x++)
6a029d8dcfa8 Cosmetics: remove braces and useless newline
vitor
parents: 7433
diff changeset
103 ractx->sb[4-x] = av_clipf(ractx->sb[4-x] + buffer[x], -4095, 4095);
1304
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
104 }
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
105
7369
b0d8a3985f96 Rename pred() to eval_lpc_coeffs() and add a doxy comment
vitor
parents: 7368
diff changeset
106 /**
b0d8a3985f96 Rename pred() to eval_lpc_coeffs() and add a doxy comment
vitor
parents: 7368
diff changeset
107 * Converts autocorrelation coefficients to LPC coefficients using the
b0d8a3985f96 Rename pred() to eval_lpc_coeffs() and add a doxy comment
vitor
parents: 7368
diff changeset
108 * Levinson-Durbin algorithm. See blocks 37 and 50 of the G.728 specification.
b0d8a3985f96 Rename pred() to eval_lpc_coeffs() and add a doxy comment
vitor
parents: 7368
diff changeset
109 *
7370
61df1ddb1001 Follow FFmpeg convention of returning negative values on error in eval_lpc_coeffs()
vitor
parents: 7369
diff changeset
110 * @return 0 if success, -1 if fail
7369
b0d8a3985f96 Rename pred() to eval_lpc_coeffs() and add a doxy comment
vitor
parents: 7368
diff changeset
111 */
b0d8a3985f96 Rename pred() to eval_lpc_coeffs() and add a doxy comment
vitor
parents: 7368
diff changeset
112 static int eval_lpc_coeffs(const float *in, float *tgt, int n)
1304
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
113 {
7167
9b93a3cf6440 Reindent the whole file
vitor
parents: 7040
diff changeset
114 int x, y;
9b93a3cf6440 Reindent the whole file
vitor
parents: 7040
diff changeset
115 double f0, f1, f2;
1304
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
116
7167
9b93a3cf6440 Reindent the whole file
vitor
parents: 7040
diff changeset
117 if (in[n] == 0)
7370
61df1ddb1001 Follow FFmpeg convention of returning negative values on error in eval_lpc_coeffs()
vitor
parents: 7369
diff changeset
118 return -1;
1304
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
119
7167
9b93a3cf6440 Reindent the whole file
vitor
parents: 7040
diff changeset
120 if ((f0 = *in) <= 0)
7370
61df1ddb1001 Follow FFmpeg convention of returning negative values on error in eval_lpc_coeffs()
vitor
parents: 7369
diff changeset
121 return -1;
7167
9b93a3cf6440 Reindent the whole file
vitor
parents: 7040
diff changeset
122
7296
6221f1f81f9b Simplify
vitor
parents: 7295
diff changeset
123 in--; // To avoid a -1 subtraction in the inner loop
6221f1f81f9b Simplify
vitor
parents: 7295
diff changeset
124
7295
f3923fea3a8c Simplify
vitor
parents: 7294
diff changeset
125 for (x=1; x <= n; x++) {
7296
6221f1f81f9b Simplify
vitor
parents: 7295
diff changeset
126 f1 = in[x+1];
7195
5e5c1178ae9a Write for loops in a more standard way. In my opinion, they are much more readable now.
vitor
parents: 7193
diff changeset
127
5e5c1178ae9a Write for loops in a more standard way. In my opinion, they are much more readable now.
vitor
parents: 7193
diff changeset
128 for (y=0; y < x - 1; y++)
7296
6221f1f81f9b Simplify
vitor
parents: 7295
diff changeset
129 f1 += in[x-y]*tgt[y];
1304
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
130
7299
1b5bc92244de Simplify: do not overuse pointer aritmetic
vitor
parents: 7298
diff changeset
131 tgt[x-1] = f2 = -f1/f0;
1b5bc92244de Simplify: do not overuse pointer aritmetic
vitor
parents: 7298
diff changeset
132 for (y=0; y < x >> 1; y++) {
1b5bc92244de Simplify: do not overuse pointer aritmetic
vitor
parents: 7298
diff changeset
133 float temp = tgt[y] + tgt[x-y-2]*f2;
1b5bc92244de Simplify: do not overuse pointer aritmetic
vitor
parents: 7298
diff changeset
134 tgt[x-y-2] += tgt[y]*f2;
1b5bc92244de Simplify: do not overuse pointer aritmetic
vitor
parents: 7298
diff changeset
135 tgt[y] = temp;
7167
9b93a3cf6440 Reindent the whole file
vitor
parents: 7040
diff changeset
136 }
9b93a3cf6440 Reindent the whole file
vitor
parents: 7040
diff changeset
137 if ((f0 += f1*f2) < 0)
7370
61df1ddb1001 Follow FFmpeg convention of returning negative values on error in eval_lpc_coeffs()
vitor
parents: 7369
diff changeset
138 return -1;
1304
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
139 }
7295
f3923fea3a8c Simplify
vitor
parents: 7294
diff changeset
140
7370
61df1ddb1001 Follow FFmpeg convention of returning negative values on error in eval_lpc_coeffs()
vitor
parents: 7369
diff changeset
141 return 0;
1304
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
142 }
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
143
7294
bf3952d5e76c Declare function parameters as const
vitor
parents: 7275
diff changeset
144 static void prodsum(float *tgt, const float *src, int len, int n)
7169
f12513065d85 Move function to avoid forward declaration
vitor
parents: 7168
diff changeset
145 {
7201
4d441356fa23 Create a function for float dot product instead of duplicating all over the file
vitor
parents: 7195
diff changeset
146 for (; n >= 0; n--)
4d441356fa23 Create a function for float dot product instead of duplicating all over the file
vitor
parents: 7195
diff changeset
147 tgt[n] = scalar_product_float(src, src - n, len);
7169
f12513065d85 Move function to avoid forward declaration
vitor
parents: 7168
diff changeset
148
f12513065d85 Move function to avoid forward declaration
vitor
parents: 7168
diff changeset
149 }
f12513065d85 Move function to avoid forward declaration
vitor
parents: 7168
diff changeset
150
7367
f8e1a1ef4a56 Add comments to do_hybrid_window()
vitor
parents: 7366
diff changeset
151 /**
f8e1a1ef4a56 Add comments to do_hybrid_window()
vitor
parents: 7366
diff changeset
152 * Hybrid window filtering. See blocks 36 and 49 of the G.728 specification.
f8e1a1ef4a56 Add comments to do_hybrid_window()
vitor
parents: 7366
diff changeset
153 *
7433
21ecdb66d5b0 Change the way the input is passed to do_hybrid_filter(). Before, in[0] was
vitor
parents: 7432
diff changeset
154 * @note This function is slightly different from that described in the spec.
21ecdb66d5b0 Change the way the input is passed to do_hybrid_filter(). Before, in[0] was
vitor
parents: 7432
diff changeset
155 * It expects in[0] to be the newest sample and in[n-1] to be the oldest
21ecdb66d5b0 Change the way the input is passed to do_hybrid_filter(). Before, in[0] was
vitor
parents: 7432
diff changeset
156 * one stored. The spec has in the more ordinary way (in[0] the oldest
21ecdb66d5b0 Change the way the input is passed to do_hybrid_filter(). Before, in[0] was
vitor
parents: 7432
diff changeset
157 * and in[n-1] the newest).
21ecdb66d5b0 Change the way the input is passed to do_hybrid_filter(). Before, in[0] was
vitor
parents: 7432
diff changeset
158 *
7367
f8e1a1ef4a56 Add comments to do_hybrid_window()
vitor
parents: 7366
diff changeset
159 * @param order the order of the filter
f8e1a1ef4a56 Add comments to do_hybrid_window()
vitor
parents: 7366
diff changeset
160 * @param n the length of the input
7375
981b14d8b0e9 Cosmetics: grammar nits
vitor
parents: 7372
diff changeset
161 * @param non_rec the number of non-recursive samples
7367
f8e1a1ef4a56 Add comments to do_hybrid_window()
vitor
parents: 7366
diff changeset
162 * @param out the filter output
f8e1a1ef4a56 Add comments to do_hybrid_window()
vitor
parents: 7366
diff changeset
163 * @param in pointer to the input of the filter
f8e1a1ef4a56 Add comments to do_hybrid_window()
vitor
parents: 7366
diff changeset
164 * @param hist pointer to the input history of the filter. It is updated by
f8e1a1ef4a56 Add comments to do_hybrid_window()
vitor
parents: 7366
diff changeset
165 * this function.
7375
981b14d8b0e9 Cosmetics: grammar nits
vitor
parents: 7372
diff changeset
166 * @param out pointer to the non-recursive part of the output
7367
f8e1a1ef4a56 Add comments to do_hybrid_window()
vitor
parents: 7366
diff changeset
167 * @param out2 pointer to the recursive part of the output
f8e1a1ef4a56 Add comments to do_hybrid_window()
vitor
parents: 7366
diff changeset
168 * @param window pointer to the windowing function table
f8e1a1ef4a56 Add comments to do_hybrid_window()
vitor
parents: 7366
diff changeset
169 */
7366
d74804a3a18b Give parameters of do_hybrid_window() more meaningful names
vitor
parents: 7365
diff changeset
170 static void do_hybrid_window(int order, int n, int non_rec, const float *in,
d74804a3a18b Give parameters of do_hybrid_window() more meaningful names
vitor
parents: 7365
diff changeset
171 float *out, float *hist, float *out2,
d74804a3a18b Give parameters of do_hybrid_window() more meaningful names
vitor
parents: 7365
diff changeset
172 const float *window)
1304
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
173 {
7167
9b93a3cf6440 Reindent the whole file
vitor
parents: 7040
diff changeset
174 unsigned int x;
9b93a3cf6440 Reindent the whole file
vitor
parents: 7040
diff changeset
175 float buffer1[37];
9b93a3cf6440 Reindent the whole file
vitor
parents: 7040
diff changeset
176 float buffer2[37];
9b93a3cf6440 Reindent the whole file
vitor
parents: 7040
diff changeset
177 float work[111];
1304
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
178
7367
f8e1a1ef4a56 Add comments to do_hybrid_window()
vitor
parents: 7366
diff changeset
179 /* update history */
7366
d74804a3a18b Give parameters of do_hybrid_window() more meaningful names
vitor
parents: 7365
diff changeset
180 memmove(hist , hist + n, (order + non_rec)*sizeof(*hist));
7433
21ecdb66d5b0 Change the way the input is passed to do_hybrid_filter(). Before, in[0] was
vitor
parents: 7432
diff changeset
181
21ecdb66d5b0 Change the way the input is passed to do_hybrid_filter(). Before, in[0] was
vitor
parents: 7432
diff changeset
182 for (x=0; x < n; x++)
21ecdb66d5b0 Change the way the input is passed to do_hybrid_filter(). Before, in[0] was
vitor
parents: 7432
diff changeset
183 hist[order + non_rec + x] = in[n-x-1];
7363
d346c4da00bd Simplify co(), use memcpy/memmove and colmult() when useful.
vitor
parents: 7322
diff changeset
184
7366
d74804a3a18b Give parameters of do_hybrid_window() more meaningful names
vitor
parents: 7365
diff changeset
185 colmult(work, window, hist, order + n + non_rec);
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2956
diff changeset
186
7366
d74804a3a18b Give parameters of do_hybrid_window() more meaningful names
vitor
parents: 7365
diff changeset
187 prodsum(buffer1, work + order , n , order);
d74804a3a18b Give parameters of do_hybrid_window() more meaningful names
vitor
parents: 7365
diff changeset
188 prodsum(buffer2, work + order + n, non_rec, order);
1304
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
189
7366
d74804a3a18b Give parameters of do_hybrid_window() more meaningful names
vitor
parents: 7365
diff changeset
190 for (x=0; x <= order; x++) {
d74804a3a18b Give parameters of do_hybrid_window() more meaningful names
vitor
parents: 7365
diff changeset
191 out2[x] = out2[x] * 0.5625 + buffer1[x];
d74804a3a18b Give parameters of do_hybrid_window() more meaningful names
vitor
parents: 7365
diff changeset
192 out [x] = out2[x] + buffer2[x];
7167
9b93a3cf6440 Reindent the whole file
vitor
parents: 7040
diff changeset
193 }
7365
031d4866b192 Rename co() function to the much more descriptive name of
vitor
parents: 7364
diff changeset
194
031d4866b192 Rename co() function to the much more descriptive name of
vitor
parents: 7364
diff changeset
195 /* Multiply by the white noise correcting factor (WNCF) */
031d4866b192 Rename co() function to the much more descriptive name of
vitor
parents: 7364
diff changeset
196 *out *= 257./256.;
1304
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
197 }
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
198
7372
e97d0795ee70 Rename update() function to backward_filter() and add a doxy comment
vitor
parents: 7371
diff changeset
199 /**
e97d0795ee70 Rename update() function to backward_filter() and add a doxy comment
vitor
parents: 7371
diff changeset
200 * Backward synthesis filter. Find the LPC coefficients from past speech data.
e97d0795ee70 Rename update() function to backward_filter() and add a doxy comment
vitor
parents: 7371
diff changeset
201 */
7412
51da3d8200d0 Give the context variables of ra288 a better name: s/glob/ractx/
vitor
parents: 7411
diff changeset
202 static void backward_filter(RA288Context *ractx)
1304
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
203 {
7428
e8eca33561c7 Rename two context buffers: pr{1,2} -> {sp,gain}_lpc. Also add a few doxy
vitor
parents: 7426
diff changeset
204 float temp1[37]; // RTMP in the spec
e8eca33561c7 Rename two context buffers: pr{1,2} -> {sp,gain}_lpc. Also add a few doxy
vitor
parents: 7426
diff changeset
205 float temp2[11]; // GPTPMP in the spec
7430
4ebfb53ff2df Remove the history buffer from the context. It can easily be evaluated
vitor
parents: 7428
diff changeset
206
7433
21ecdb66d5b0 Change the way the input is passed to do_hybrid_filter(). Before, in[0] was
vitor
parents: 7432
diff changeset
207 do_hybrid_window(36, 40, 35, ractx->sb, temp1, ractx->sp_hist,
7428
e8eca33561c7 Rename two context buffers: pr{1,2} -> {sp,gain}_lpc. Also add a few doxy
vitor
parents: 7426
diff changeset
208 ractx->sp_rec, syn_window);
1304
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
209
7428
e8eca33561c7 Rename two context buffers: pr{1,2} -> {sp,gain}_lpc. Also add a few doxy
vitor
parents: 7426
diff changeset
210 if (!eval_lpc_coeffs(temp1, ractx->sp_lpc, 36))
e8eca33561c7 Rename two context buffers: pr{1,2} -> {sp,gain}_lpc. Also add a few doxy
vitor
parents: 7426
diff changeset
211 colmult(ractx->sp_lpc, ractx->sp_lpc, syn_bw_tab, 36);
7169
f12513065d85 Move function to avoid forward declaration
vitor
parents: 7168
diff changeset
212
7433
21ecdb66d5b0 Change the way the input is passed to do_hybrid_filter(). Before, in[0] was
vitor
parents: 7432
diff changeset
213 do_hybrid_window(10, 8, 20, ractx->lhist, temp2, ractx->gain_hist,
7428
e8eca33561c7 Rename two context buffers: pr{1,2} -> {sp,gain}_lpc. Also add a few doxy
vitor
parents: 7426
diff changeset
214 ractx->gain_rec, gain_window);
7169
f12513065d85 Move function to avoid forward declaration
vitor
parents: 7168
diff changeset
215
7428
e8eca33561c7 Rename two context buffers: pr{1,2} -> {sp,gain}_lpc. Also add a few doxy
vitor
parents: 7426
diff changeset
216 if (!eval_lpc_coeffs(temp2, ractx->gain_lpc, 10))
e8eca33561c7 Rename two context buffers: pr{1,2} -> {sp,gain}_lpc. Also add a few doxy
vitor
parents: 7426
diff changeset
217 colmult(ractx->gain_lpc, ractx->gain_lpc, gain_bw_tab, 10);
1304
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
218 }
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
219
7174
7c4589349c01 Remove useless wrapper around ra288_decode_frame()
vitor
parents: 7173
diff changeset
220 /* Decode a block (celp) */
7c4589349c01 Remove useless wrapper around ra288_decode_frame()
vitor
parents: 7173
diff changeset
221 static int ra288_decode_frame(AVCodecContext * avctx, void *data,
7c4589349c01 Remove useless wrapper around ra288_decode_frame()
vitor
parents: 7173
diff changeset
222 int *data_size, const uint8_t * buf,
7c4589349c01 Remove useless wrapper around ra288_decode_frame()
vitor
parents: 7173
diff changeset
223 int buf_size)
1313
d6da658fe6c6 finally works now
nickols_k
parents: 1308
diff changeset
224 {
7174
7c4589349c01 Remove useless wrapper around ra288_decode_frame()
vitor
parents: 7173
diff changeset
225 int16_t *out = data;
7167
9b93a3cf6440 Reindent the whole file
vitor
parents: 7040
diff changeset
226 int x, y;
7412
51da3d8200d0 Give the context variables of ra288 a better name: s/glob/ractx/
vitor
parents: 7411
diff changeset
227 RA288Context *ractx = avctx->priv_data;
7173
bc2fd265f52b Remove unpack() function, read the bitstream as needed
vitor
parents: 7172
diff changeset
228 GetBitContext gb;
7167
9b93a3cf6440 Reindent the whole file
vitor
parents: 7040
diff changeset
229
7174
7c4589349c01 Remove useless wrapper around ra288_decode_frame()
vitor
parents: 7173
diff changeset
230 if (buf_size < avctx->block_align) {
7c4589349c01 Remove useless wrapper around ra288_decode_frame()
vitor
parents: 7173
diff changeset
231 av_log(avctx, AV_LOG_ERROR,
7c4589349c01 Remove useless wrapper around ra288_decode_frame()
vitor
parents: 7173
diff changeset
232 "Error! Input buffer is too small [%d<%d]\n",
7c4589349c01 Remove useless wrapper around ra288_decode_frame()
vitor
parents: 7173
diff changeset
233 buf_size, avctx->block_align);
7c4589349c01 Remove useless wrapper around ra288_decode_frame()
vitor
parents: 7173
diff changeset
234 return 0;
7c4589349c01 Remove useless wrapper around ra288_decode_frame()
vitor
parents: 7173
diff changeset
235 }
7c4589349c01 Remove useless wrapper around ra288_decode_frame()
vitor
parents: 7173
diff changeset
236
7c4589349c01 Remove useless wrapper around ra288_decode_frame()
vitor
parents: 7173
diff changeset
237 init_get_bits(&gb, buf, avctx->block_align * 8);
1313
d6da658fe6c6 finally works now
nickols_k
parents: 1308
diff changeset
238
7167
9b93a3cf6440 Reindent the whole file
vitor
parents: 7040
diff changeset
239 for (x=0; x < 32; x++) {
7192
9517beb3b491 Simplify
vitor
parents: 7191
diff changeset
240 float gain = amptable[get_bits(&gb, 3)];
7173
bc2fd265f52b Remove unpack() function, read the bitstream as needed
vitor
parents: 7172
diff changeset
241 int cb_coef = get_bits(&gb, 6 + (x&1));
7425
ed615bbf4f8d Change slightly the meaning of RA288Context.phase. This
vitor
parents: 7419
diff changeset
242 ractx->phase = (x + 4) & 7;
7412
51da3d8200d0 Give the context variables of ra288 a better name: s/glob/ractx/
vitor
parents: 7411
diff changeset
243 decode(ractx, gain, cb_coef);
7167
9b93a3cf6440 Reindent the whole file
vitor
parents: 7040
diff changeset
244
7195
5e5c1178ae9a Write for loops in a more standard way. In my opinion, they are much more readable now.
vitor
parents: 7193
diff changeset
245 for (y=0; y < 5; y++)
7432
158d1acfe3bf Remove RA288Context.output buffer. This buffer is just RA288Context.sb
vitor
parents: 7431
diff changeset
246 *(out++) = 8 * ractx->sb[4 - y];
7167
9b93a3cf6440 Reindent the whole file
vitor
parents: 7040
diff changeset
247
7425
ed615bbf4f8d Change slightly the meaning of RA288Context.phase. This
vitor
parents: 7419
diff changeset
248 if (ractx->phase == 7)
7412
51da3d8200d0 Give the context variables of ra288 a better name: s/glob/ractx/
vitor
parents: 7411
diff changeset
249 backward_filter(ractx);
7167
9b93a3cf6440 Reindent the whole file
vitor
parents: 7040
diff changeset
250 }
9b93a3cf6440 Reindent the whole file
vitor
parents: 7040
diff changeset
251
7174
7c4589349c01 Remove useless wrapper around ra288_decode_frame()
vitor
parents: 7173
diff changeset
252 *data_size = (char *)out - (char *)data;
2956
5f51b1e0bed6 Cook compatibe decoder, patch by Benjamin Larsson
rtognimp
parents: 2288
diff changeset
253 return avctx->block_align;
1304
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
254 }
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
255
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
256 AVCodec ra_288_decoder =
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
257 {
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
258 "real_288",
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
259 CODEC_TYPE_AUDIO,
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
260 CODEC_ID_RA_288,
7411
d41948bf80a1 Give the context structure of ra288 a better name: s/Real288_internal/RA288Context/
vitor
parents: 7375
diff changeset
261 sizeof(RA288Context),
7168
c2691f22708b Remove useless initialization
vitor
parents: 7167
diff changeset
262 NULL,
1304
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
263 NULL,
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
264 NULL,
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
265 ra288_decode_frame,
7040
e943e1409077 Make AVCodec long_names definition conditional depending on CONFIG_SMALL.
stefano
parents: 6712
diff changeset
266 .long_name = NULL_IF_CONFIG_SMALL("RealAudio 2.0 (28.8K)"),
1304
e8543aab0cc9 RealAudio 1.0 (14_4) and 2.0 (28_8) native decoders. Author unknown, code from a mplayerhq (originally from public domain player for Amiga
nickols_k
parents:
diff changeset
267 };