annotate wmavoice.c @ 12274:1d207bb5cd29 libavcodec

Use nested ifs instead of &&, which appears to not work with %ifidn (i.e. this construct was always enabled, even for <ssse3 versions).
author rbultje
date Mon, 26 Jul 2010 13:56:51 +0000
parents 7323559a53fd
children 2ba9068e748d
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
11123
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1 /*
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
2 * Windows Media Audio Voice decoder.
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
3 * Copyright (c) 2009 Ronald S. Bultje
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
4 *
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
5 * This file is part of FFmpeg.
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
6 *
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
7 * FFmpeg is free software; you can redistribute it and/or
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
8 * modify it under the terms of the GNU Lesser General Public
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
9 * License as published by the Free Software Foundation; either
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
10 * version 2.1 of the License, or (at your option) any later version.
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
11 *
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
12 * FFmpeg is distributed in the hope that it will be useful,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
13 * but WITHOUT ANY WARRANTY; without even the implied warranty of
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
14 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
15 * Lesser General Public License for more details.
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
16 *
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
17 * You should have received a copy of the GNU Lesser General Public
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
18 * License along with FFmpeg; if not, write to the Free Software
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
19 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
20 */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
21
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
22 /**
11644
7dd2a45249a9 Remove explicit filename from Doxygen @file commands.
diego
parents: 11560
diff changeset
23 * @file
11123
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
24 * @brief Windows Media Audio Voice compatible decoder
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
25 * @author Ronald S. Bultje <rsbultje@gmail.com>
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
26 */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
27
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
28 #include <math.h>
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
29 #include "avcodec.h"
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
30 #include "get_bits.h"
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
31 #include "put_bits.h"
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
32 #include "wmavoice_data.h"
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
33 #include "celp_math.h"
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
34 #include "celp_filters.h"
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
35 #include "acelp_vectors.h"
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
36 #include "acelp_filters.h"
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
37 #include "lsp.h"
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
38 #include "libavutil/lzo.h"
11653
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
39 #include "avfft.h"
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
40 #include "fft.h"
11123
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
41
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
42 #define MAX_BLOCKS 8 ///< maximum number of blocks per frame
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
43 #define MAX_LSPS 16 ///< maximum filter order
11653
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
44 #define MAX_LSPS_ALIGN16 16 ///< same as #MAX_LSPS; needs to be multiple
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
45 ///< of 16 for ASM input buffer alignment
11123
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
46 #define MAX_FRAMES 3 ///< maximum number of frames per superframe
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
47 #define MAX_FRAMESIZE 160 ///< maximum number of samples per frame
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
48 #define MAX_SIGNAL_HISTORY 416 ///< maximum excitation signal history
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
49 #define MAX_SFRAMESIZE (MAX_FRAMESIZE * MAX_FRAMES)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
50 ///< maximum number of samples per superframe
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
51 #define SFRAME_CACHE_MAXSIZE 256 ///< maximum cache size for frame data that
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
52 ///< was split over two packets
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
53 #define VLC_NBITS 6 ///< number of bits to read per VLC iteration
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
54
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
55 /**
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
56 * Frame type VLC coding.
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
57 */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
58 static VLC frame_type_vlc;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
59
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
60 /**
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
61 * Adaptive codebook types.
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
62 */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
63 enum {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
64 ACB_TYPE_NONE = 0, ///< no adaptive codebook (only hardcoded fixed)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
65 ACB_TYPE_ASYMMETRIC = 1, ///< adaptive codebook with per-frame pitch, which
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
66 ///< we interpolate to get a per-sample pitch.
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
67 ///< Signal is generated using an asymmetric sinc
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
68 ///< window function
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
69 ///< @note see #wmavoice_ipol1_coeffs
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
70 ACB_TYPE_HAMMING = 2 ///< Per-block pitch with signal generation using
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
71 ///< a Hamming sinc window function
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
72 ///< @note see #wmavoice_ipol2_coeffs
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
73 };
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
74
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
75 /**
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
76 * Fixed codebook types.
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
77 */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
78 enum {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
79 FCB_TYPE_SILENCE = 0, ///< comfort noise during silence
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
80 ///< generated from a hardcoded (fixed) codebook
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
81 ///< with per-frame (low) gain values
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
82 FCB_TYPE_HARDCODED = 1, ///< hardcoded (fixed) codebook with per-block
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
83 ///< gain values
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
84 FCB_TYPE_AW_PULSES = 2, ///< Pitch-adaptive window (AW) pulse signals,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
85 ///< used in particular for low-bitrate streams
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
86 FCB_TYPE_EXC_PULSES = 3, ///< Innovation (fixed) codebook pulse sets in
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
87 ///< combinations of either single pulses or
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
88 ///< pulse pairs
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
89 };
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
90
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
91 /**
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
92 * Description of frame types.
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
93 */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
94 static const struct frame_type_desc {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
95 uint8_t n_blocks; ///< amount of blocks per frame (each block
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
96 ///< (contains 160/#n_blocks samples)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
97 uint8_t log_n_blocks; ///< log2(#n_blocks)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
98 uint8_t acb_type; ///< Adaptive codebook type (ACB_TYPE_*)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
99 uint8_t fcb_type; ///< Fixed codebook type (FCB_TYPE_*)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
100 uint8_t dbl_pulses; ///< how many pulse vectors have pulse pairs
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
101 ///< (rather than just one single pulse)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
102 ///< only if #fcb_type == #FCB_TYPE_EXC_PULSES
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
103 uint16_t frame_size; ///< the amount of bits that make up the block
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
104 ///< data (per frame)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
105 } frame_descs[17] = {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
106 { 1, 0, ACB_TYPE_NONE, FCB_TYPE_SILENCE, 0, 0 },
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
107 { 2, 1, ACB_TYPE_NONE, FCB_TYPE_HARDCODED, 0, 28 },
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
108 { 2, 1, ACB_TYPE_ASYMMETRIC, FCB_TYPE_AW_PULSES, 0, 46 },
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
109 { 2, 1, ACB_TYPE_ASYMMETRIC, FCB_TYPE_EXC_PULSES, 2, 80 },
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
110 { 2, 1, ACB_TYPE_ASYMMETRIC, FCB_TYPE_EXC_PULSES, 5, 104 },
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
111 { 4, 2, ACB_TYPE_ASYMMETRIC, FCB_TYPE_EXC_PULSES, 0, 108 },
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
112 { 4, 2, ACB_TYPE_ASYMMETRIC, FCB_TYPE_EXC_PULSES, 2, 132 },
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
113 { 4, 2, ACB_TYPE_ASYMMETRIC, FCB_TYPE_EXC_PULSES, 5, 168 },
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
114 { 2, 1, ACB_TYPE_HAMMING, FCB_TYPE_EXC_PULSES, 0, 64 },
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
115 { 2, 1, ACB_TYPE_HAMMING, FCB_TYPE_EXC_PULSES, 2, 80 },
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
116 { 2, 1, ACB_TYPE_HAMMING, FCB_TYPE_EXC_PULSES, 5, 104 },
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
117 { 4, 2, ACB_TYPE_HAMMING, FCB_TYPE_EXC_PULSES, 0, 108 },
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
118 { 4, 2, ACB_TYPE_HAMMING, FCB_TYPE_EXC_PULSES, 2, 132 },
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
119 { 4, 2, ACB_TYPE_HAMMING, FCB_TYPE_EXC_PULSES, 5, 168 },
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
120 { 8, 3, ACB_TYPE_HAMMING, FCB_TYPE_EXC_PULSES, 0, 176 },
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
121 { 8, 3, ACB_TYPE_HAMMING, FCB_TYPE_EXC_PULSES, 2, 208 },
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
122 { 8, 3, ACB_TYPE_HAMMING, FCB_TYPE_EXC_PULSES, 5, 256 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
123 };
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
124
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
125 /**
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
126 * WMA Voice decoding context.
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
127 */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
128 typedef struct {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
129 /**
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
130 * @defgroup struct_global Global values
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
131 * Global values, specified in the stream header / extradata or used
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
132 * all over.
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
133 * @{
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
134 */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
135 GetBitContext gb; ///< packet bitreader. During decoder init,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
136 ///< it contains the extradata from the
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
137 ///< demuxer. During decoding, it contains
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
138 ///< packet data.
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
139 int8_t vbm_tree[25]; ///< converts VLC codes to frame type
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
140
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
141 int spillover_bitsize; ///< number of bits used to specify
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
142 ///< #spillover_nbits in the packet header
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
143 ///< = ceil(log2(ctx->block_align << 3))
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
144 int history_nsamples; ///< number of samples in history for signal
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
145 ///< prediction (through ACB)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
146
11653
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
147 /* postfilter specific values */
11123
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
148 int do_apf; ///< whether to apply the averaged
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
149 ///< projection filter (APF)
11653
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
150 int denoise_strength; ///< strength of denoising in Wiener filter
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
151 ///< [0-11]
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
152 int denoise_tilt_corr; ///< Whether to apply tilt correction to the
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
153 ///< Wiener filter coefficients (postfilter)
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
154 int dc_level; ///< Predicted amount of DC noise, based
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
155 ///< on which a DC removal filter is used
11123
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
156
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
157 int lsps; ///< number of LSPs per frame [10 or 16]
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
158 int lsp_q_mode; ///< defines quantizer defaults [0, 1]
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
159 int lsp_def_mode; ///< defines different sets of LSP defaults
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
160 ///< [0, 1]
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
161 int frame_lsp_bitsize; ///< size (in bits) of LSPs, when encoded
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
162 ///< per-frame (independent coding)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
163 int sframe_lsp_bitsize; ///< size (in bits) of LSPs, when encoded
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
164 ///< per superframe (residual coding)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
165
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
166 int min_pitch_val; ///< base value for pitch parsing code
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
167 int max_pitch_val; ///< max value + 1 for pitch parsing
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
168 int pitch_nbits; ///< number of bits used to specify the
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
169 ///< pitch value in the frame header
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
170 int block_pitch_nbits; ///< number of bits used to specify the
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
171 ///< first block's pitch value
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
172 int block_pitch_range; ///< range of the block pitch
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
173 int block_delta_pitch_nbits; ///< number of bits used to specify the
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
174 ///< delta pitch between this and the last
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
175 ///< block's pitch value, used in all but
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
176 ///< first block
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
177 int block_delta_pitch_hrange; ///< 1/2 range of the delta (full range is
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
178 ///< from -this to +this-1)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
179 uint16_t block_conv_table[4]; ///< boundaries for block pitch unit/scale
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
180 ///< conversion
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
181
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
182 /**
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
183 * @}
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
184 * @defgroup struct_packet Packet values
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
185 * Packet values, specified in the packet header or related to a packet.
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
186 * A packet is considered to be a single unit of data provided to this
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
187 * decoder by the demuxer.
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
188 * @{
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
189 */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
190 int spillover_nbits; ///< number of bits of the previous packet's
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
191 ///< last superframe preceeding this
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
192 ///< packet's first full superframe (useful
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
193 ///< for re-synchronization also)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
194 int has_residual_lsps; ///< if set, superframes contain one set of
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
195 ///< LSPs that cover all frames, encoded as
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
196 ///< independent and residual LSPs; if not
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
197 ///< set, each frame contains its own, fully
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
198 ///< independent, LSPs
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
199 int skip_bits_next; ///< number of bits to skip at the next call
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
200 ///< to #wmavoice_decode_packet() (since
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
201 ///< they're part of the previous superframe)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
202
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
203 uint8_t sframe_cache[SFRAME_CACHE_MAXSIZE + FF_INPUT_BUFFER_PADDING_SIZE];
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
204 ///< cache for superframe data split over
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
205 ///< multiple packets
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
206 int sframe_cache_size; ///< set to >0 if we have data from an
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
207 ///< (incomplete) superframe from a previous
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
208 ///< packet that spilled over in the current
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
209 ///< packet; specifies the amount of bits in
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
210 ///< #sframe_cache
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
211 PutBitContext pb; ///< bitstream writer for #sframe_cache
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
212
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
213 /**
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
214 * @}
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
215 * @defgroup struct_frame Frame and superframe values
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
216 * Superframe and frame data - these can change from frame to frame,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
217 * although some of them do in that case serve as a cache / history for
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
218 * the next frame or superframe.
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
219 * @{
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
220 */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
221 double prev_lsps[MAX_LSPS]; ///< LSPs of the last frame of the previous
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
222 ///< superframe
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
223 int last_pitch_val; ///< pitch value of the previous frame
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
224 int last_acb_type; ///< frame type [0-2] of the previous frame
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
225 int pitch_diff_sh16; ///< ((cur_pitch_val - #last_pitch_val)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
226 ///< << 16) / #MAX_FRAMESIZE
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
227 float silence_gain; ///< set for use in blocks if #ACB_TYPE_NONE
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
228
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
229 int aw_idx_is_ext; ///< whether the AW index was encoded in
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
230 ///< 8 bits (instead of 6)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
231 int aw_pulse_range; ///< the range over which #aw_pulse_set1()
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
232 ///< can apply the pulse, relative to the
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
233 ///< value in aw_first_pulse_off. The exact
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
234 ///< position of the first AW-pulse is within
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
235 ///< [pulse_off, pulse_off + this], and
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
236 ///< depends on bitstream values; [16 or 24]
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
237 int aw_n_pulses[2]; ///< number of AW-pulses in each block; note
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
238 ///< that this number can be negative (in
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
239 ///< which case it basically means "zero")
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
240 int aw_first_pulse_off[2]; ///< index of first sample to which to
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
241 ///< apply AW-pulses, or -0xff if unset
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
242 int aw_next_pulse_off_cache; ///< the position (relative to start of the
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
243 ///< second block) at which pulses should
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
244 ///< start to be positioned, serves as a
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
245 ///< cache for pitch-adaptive window pulses
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
246 ///< between blocks
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
247
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
248 int frame_cntr; ///< current frame index [0 - 0xFFFE]; is
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
249 ///< only used for comfort noise in #pRNG()
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
250 float gain_pred_err[6]; ///< cache for gain prediction
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
251 float excitation_history[MAX_SIGNAL_HISTORY];
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
252 ///< cache of the signal of previous
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
253 ///< superframes, used as a history for
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
254 ///< signal generation
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
255 float synth_history[MAX_LSPS]; ///< see #excitation_history
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
256 /**
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
257 * @}
11653
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
258 * @defgroup post_filter Postfilter values
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
259 * Varibales used for postfilter implementation, mostly history for
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
260 * smoothing and so on, and context variables for FFT/iFFT.
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
261 * @{
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
262 */
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
263 RDFTContext rdft, irdft; ///< contexts for FFT-calculation in the
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
264 ///< postfilter (for denoise filter)
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
265 DCTContext dct, dst; ///< contexts for phase shift (in Hilbert
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
266 ///< transform, part of postfilter)
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
267 float sin[511], cos[511]; ///< 8-bit cosine/sine windows over [-pi,pi]
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
268 ///< range
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
269 float postfilter_agc; ///< gain control memory, used in
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
270 ///< #adaptive_gain_control()
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
271 float dcf_mem[2]; ///< DC filter history
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
272 float zero_exc_pf[MAX_SIGNAL_HISTORY + MAX_SFRAMESIZE];
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
273 ///< zero filter output (i.e. excitation)
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
274 ///< by postfilter
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
275 float denoise_filter_cache[MAX_FRAMESIZE];
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
276 int denoise_filter_cache_size; ///< samples in #denoise_filter_cache
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
277 DECLARE_ALIGNED(16, float, tilted_lpcs_pf)[0x80];
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
278 ///< aligned buffer for LPC tilting
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
279 DECLARE_ALIGNED(16, float, denoise_coeffs_pf)[0x80];
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
280 ///< aligned buffer for denoise coefficients
11673
da59da296153 Fix buffer overrun (or, well, actually a typo, 80 should be 0x80...).
rbultje
parents: 11653
diff changeset
281 DECLARE_ALIGNED(16, float, synth_filter_out_buf)[0x80 + MAX_LSPS_ALIGN16];
11653
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
282 ///< aligned buffer for postfilter speech
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
283 ///< synthesis
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
284 /**
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
285 * @}
11123
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
286 */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
287 } WMAVoiceContext;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
288
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
289 /**
12024
fdafbcef52f5 Fix grammar errors in documentation
mru
parents: 12021
diff changeset
290 * Set up the variable bit mode (VBM) tree from container extradata.
11123
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
291 * @param gb bit I/O context.
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
292 * The bit context (s->gb) should be loaded with byte 23-46 of the
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
293 * container extradata (i.e. the ones containing the VBM tree).
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
294 * @param vbm_tree pointer to array to which the decoded VBM tree will be
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
295 * written.
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
296 * @return 0 on success, <0 on error.
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
297 */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
298 static av_cold int decode_vbmtree(GetBitContext *gb, int8_t vbm_tree[25])
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
299 {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
300 static const uint8_t bits[] = {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
301 2, 2, 2, 4, 4, 4,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
302 6, 6, 6, 8, 8, 8,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
303 10, 10, 10, 12, 12, 12,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
304 14, 14, 14, 14
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
305 };
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
306 static const uint16_t codes[] = {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
307 0x0000, 0x0001, 0x0002, // 00/01/10
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
308 0x000c, 0x000d, 0x000e, // 11+00/01/10
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
309 0x003c, 0x003d, 0x003e, // 1111+00/01/10
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
310 0x00fc, 0x00fd, 0x00fe, // 111111+00/01/10
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
311 0x03fc, 0x03fd, 0x03fe, // 11111111+00/01/10
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
312 0x0ffc, 0x0ffd, 0x0ffe, // 1111111111+00/01/10
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
313 0x3ffc, 0x3ffd, 0x3ffe, 0x3fff // 111111111111+xx
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
314 };
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
315 int cntr[8], n, res;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
316
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
317 memset(vbm_tree, 0xff, sizeof(vbm_tree));
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
318 memset(cntr, 0, sizeof(cntr));
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
319 for (n = 0; n < 17; n++) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
320 res = get_bits(gb, 3);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
321 if (cntr[res] > 3) // should be >= 3 + (res == 7))
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
322 return -1;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
323 vbm_tree[res * 3 + cntr[res]++] = n;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
324 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
325 INIT_VLC_STATIC(&frame_type_vlc, VLC_NBITS, sizeof(bits),
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
326 bits, 1, 1, codes, 2, 2, 132);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
327 return 0;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
328 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
329
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
330 /**
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
331 * Set up decoder with parameters from demuxer (extradata etc.).
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
332 */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
333 static av_cold int wmavoice_decode_init(AVCodecContext *ctx)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
334 {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
335 int n, flags, pitch_range, lsp16_flag;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
336 WMAVoiceContext *s = ctx->priv_data;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
337
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
338 /**
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
339 * Extradata layout:
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
340 * - byte 0-18: WMAPro-in-WMAVoice extradata (see wmaprodec.c),
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
341 * - byte 19-22: flags field (annoyingly in LE; see below for known
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
342 * values),
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
343 * - byte 23-46: variable bitmode tree (really just 17 * 3 bits,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
344 * rest is 0).
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
345 */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
346 if (ctx->extradata_size != 46) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
347 av_log(ctx, AV_LOG_ERROR,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
348 "Invalid extradata size %d (should be 46)\n",
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
349 ctx->extradata_size);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
350 return -1;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
351 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
352 flags = AV_RL32(ctx->extradata + 18);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
353 s->spillover_bitsize = 3 + av_ceil_log2(ctx->block_align);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
354 s->do_apf = flags & 0x1;
11653
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
355 if (s->do_apf) {
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
356 ff_rdft_init(&s->rdft, 7, DFT_R2C);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
357 ff_rdft_init(&s->irdft, 7, IDFT_C2R);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
358 ff_dct_init(&s->dct, 6, DCT_I);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
359 ff_dct_init(&s->dst, 6, DST_I);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
360
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
361 ff_sine_window_init(s->cos, 256);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
362 memcpy(&s->sin[255], s->cos, 256 * sizeof(s->cos[0]));
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
363 for (n = 0; n < 255; n++) {
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
364 s->sin[n] = -s->sin[510 - n];
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
365 s->cos[510 - n] = s->cos[n];
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
366 }
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
367 }
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
368 s->denoise_strength = (flags >> 2) & 0xF;
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
369 if (s->denoise_strength >= 12) {
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
370 av_log(ctx, AV_LOG_ERROR,
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
371 "Invalid denoise filter strength %d (max=11)\n",
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
372 s->denoise_strength);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
373 return -1;
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
374 }
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
375 s->denoise_tilt_corr = !!(flags & 0x40);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
376 s->dc_level = (flags >> 7) & 0xF;
11123
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
377 s->lsp_q_mode = !!(flags & 0x2000);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
378 s->lsp_def_mode = !!(flags & 0x4000);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
379 lsp16_flag = flags & 0x1000;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
380 if (lsp16_flag) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
381 s->lsps = 16;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
382 s->frame_lsp_bitsize = 34;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
383 s->sframe_lsp_bitsize = 60;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
384 } else {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
385 s->lsps = 10;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
386 s->frame_lsp_bitsize = 24;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
387 s->sframe_lsp_bitsize = 48;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
388 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
389 for (n = 0; n < s->lsps; n++)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
390 s->prev_lsps[n] = M_PI * (n + 1.0) / (s->lsps + 1.0);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
391
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
392 init_get_bits(&s->gb, ctx->extradata + 22, (ctx->extradata_size - 22) << 3);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
393 if (decode_vbmtree(&s->gb, s->vbm_tree) < 0) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
394 av_log(ctx, AV_LOG_ERROR, "Invalid VBM tree; broken extradata?\n");
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
395 return -1;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
396 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
397
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
398 s->min_pitch_val = ((ctx->sample_rate << 8) / 400 + 50) >> 8;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
399 s->max_pitch_val = ((ctx->sample_rate << 8) * 37 / 2000 + 50) >> 8;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
400 pitch_range = s->max_pitch_val - s->min_pitch_val;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
401 s->pitch_nbits = av_ceil_log2(pitch_range);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
402 s->last_pitch_val = 40;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
403 s->last_acb_type = ACB_TYPE_NONE;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
404 s->history_nsamples = s->max_pitch_val + 8;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
405
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
406 if (s->min_pitch_val < 1 || s->history_nsamples > MAX_SIGNAL_HISTORY) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
407 int min_sr = ((((1 << 8) - 50) * 400) + 0xFF) >> 8,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
408 max_sr = ((((MAX_SIGNAL_HISTORY - 8) << 8) + 205) * 2000 / 37) >> 8;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
409
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
410 av_log(ctx, AV_LOG_ERROR,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
411 "Unsupported samplerate %d (min=%d, max=%d)\n",
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
412 ctx->sample_rate, min_sr, max_sr); // 322-22097 Hz
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
413
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
414 return -1;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
415 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
416
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
417 s->block_conv_table[0] = s->min_pitch_val;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
418 s->block_conv_table[1] = (pitch_range * 25) >> 6;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
419 s->block_conv_table[2] = (pitch_range * 44) >> 6;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
420 s->block_conv_table[3] = s->max_pitch_val - 1;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
421 s->block_delta_pitch_hrange = (pitch_range >> 3) & ~0xF;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
422 s->block_delta_pitch_nbits = 1 + av_ceil_log2(s->block_delta_pitch_hrange);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
423 s->block_pitch_range = s->block_conv_table[2] +
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
424 s->block_conv_table[3] + 1 +
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
425 2 * (s->block_conv_table[1] - 2 * s->min_pitch_val);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
426 s->block_pitch_nbits = av_ceil_log2(s->block_pitch_range);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
427
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
428 ctx->sample_fmt = SAMPLE_FMT_FLT;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
429
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
430 return 0;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
431 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
432
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
433 /**
11653
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
434 * @defgroup postfilter Postfilter functions
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
435 * Postfilter functions (gain control, wiener denoise filter, DC filter,
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
436 * kalman smoothening, plus surrounding code to wrap it)
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
437 * @{
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
438 */
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
439 /**
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
440 * Adaptive gain control (as used in postfilter).
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
441 *
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
442 * Identical to #ff_adaptive_gain_control() in acelp_vectors.c, except
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
443 * that the energy here is calculated using sum(abs(...)), whereas the
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
444 * other codecs (e.g. AMR-NB, SIPRO) use sqrt(dotproduct(...)).
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
445 *
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
446 * @param out output buffer for filtered samples
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
447 * @param in input buffer containing the samples as they are after the
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
448 * postfilter steps so far
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
449 * @param speech_synth input buffer containing speech synth before postfilter
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
450 * @param size input buffer size
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
451 * @param alpha exponential filter factor
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
452 * @param gain_mem pointer to filter memory (single float)
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
453 */
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
454 static void adaptive_gain_control(float *out, const float *in,
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
455 const float *speech_synth,
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
456 int size, float alpha, float *gain_mem)
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
457 {
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
458 int i;
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
459 float speech_energy = 0.0, postfilter_energy = 0.0, gain_scale_factor;
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
460 float mem = *gain_mem;
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
461
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
462 for (i = 0; i < size; i++) {
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
463 speech_energy += fabsf(speech_synth[i]);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
464 postfilter_energy += fabsf(in[i]);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
465 }
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
466 gain_scale_factor = (1.0 - alpha) * speech_energy / postfilter_energy;
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
467
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
468 for (i = 0; i < size; i++) {
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
469 mem = alpha * mem + gain_scale_factor;
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
470 out[i] = in[i] * mem;
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
471 }
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
472
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
473 *gain_mem = mem;
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
474 }
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
475
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
476 /**
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
477 * Kalman smoothing function.
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
478 *
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
479 * This function looks back pitch +/- 3 samples back into history to find
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
480 * the best fitting curve (that one giving the optimal gain of the two
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
481 * signals, i.e. the highest dot product between the two), and then
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
482 * uses that signal history to smoothen the output of the speech synthesis
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
483 * filter.
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
484 *
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
485 * @param s WMA Voice decoding context
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
486 * @param pitch pitch of the speech signal
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
487 * @param in input speech signal
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
488 * @param out output pointer for smoothened signal
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
489 * @param size input/output buffer size
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
490 *
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
491 * @returns -1 if no smoothening took place, e.g. because no optimal
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
492 * fit could be found, or 0 on success.
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
493 */
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
494 static int kalman_smoothen(WMAVoiceContext *s, int pitch,
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
495 const float *in, float *out, int size)
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
496 {
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
497 int n;
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
498 float optimal_gain = 0, dot;
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
499 const float *ptr = &in[-FFMAX(s->min_pitch_val, pitch - 3)],
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
500 *end = &in[-FFMIN(s->max_pitch_val, pitch + 3)],
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
501 *best_hist_ptr;
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
502
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
503 /* find best fitting point in history */
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
504 do {
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
505 dot = ff_dot_productf(in, ptr, size);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
506 if (dot > optimal_gain) {
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
507 optimal_gain = dot;
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
508 best_hist_ptr = ptr;
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
509 }
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
510 } while (--ptr >= end);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
511
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
512 if (optimal_gain <= 0)
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
513 return -1;
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
514 dot = ff_dot_productf(best_hist_ptr, best_hist_ptr, size);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
515 if (dot <= 0) // would be 1.0
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
516 return -1;
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
517
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
518 if (optimal_gain <= dot) {
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
519 dot = dot / (dot + 0.6 * optimal_gain); // 0.625-1.000
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
520 } else
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
521 dot = 0.625;
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
522
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
523 /* actual smoothing */
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
524 for (n = 0; n < size; n++)
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
525 out[n] = best_hist_ptr[n] + dot * (in[n] - best_hist_ptr[n]);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
526
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
527 return 0;
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
528 }
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
529
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
530 /**
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
531 * Get the tilt factor of a formant filter from its transfer function
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
532 * @see #tilt_factor() in amrnbdec.c, which does essentially the same,
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
533 * but somehow (??) it does a speech synthesis filter in the
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
534 * middle, which is missing here
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
535 *
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
536 * @param lpcs LPC coefficients
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
537 * @param n_lpcs Size of LPC buffer
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
538 * @returns the tilt factor
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
539 */
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
540 static float tilt_factor(const float *lpcs, int n_lpcs)
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
541 {
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
542 float rh0, rh1;
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
543
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
544 rh0 = 1.0 + ff_dot_productf(lpcs, lpcs, n_lpcs);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
545 rh1 = lpcs[0] + ff_dot_productf(lpcs, &lpcs[1], n_lpcs - 1);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
546
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
547 return rh1 / rh0;
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
548 }
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
549
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
550 /**
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
551 * Derive denoise filter coefficients (in real domain) from the LPCs.
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
552 */
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
553 static void calc_input_response(WMAVoiceContext *s, float *lpcs,
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
554 int fcb_type, float *coeffs, int remainder)
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
555 {
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
556 float last_coeff, min = 15.0, max = -15.0;
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
557 float irange, angle_mul, gain_mul, range, sq;
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
558 int n, idx;
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
559
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
560 /* Create frequency power spectrum of speech input (i.e. RDFT of LPCs) */
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
561 ff_rdft_calc(&s->rdft, lpcs);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
562 #define log_range(var, assign) do { \
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
563 float tmp = log10f(assign); var = tmp; \
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
564 max = FFMAX(max, tmp); min = FFMIN(min, tmp); \
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
565 } while (0)
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
566 log_range(last_coeff, lpcs[1] * lpcs[1]);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
567 for (n = 1; n < 64; n++)
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
568 log_range(lpcs[n], lpcs[n * 2] * lpcs[n * 2] +
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
569 lpcs[n * 2 + 1] * lpcs[n * 2 + 1]);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
570 log_range(lpcs[0], lpcs[0] * lpcs[0]);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
571 #undef log_range
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
572 range = max - min;
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
573 lpcs[64] = last_coeff;
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
574
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
575 /* Now, use this spectrum to pick out these frequencies with higher
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
576 * (relative) power/energy (which we then take to be "not noise"),
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
577 * and set up a table (still in lpc[]) of (relative) gains per frequency.
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
578 * These frequencies will be maintained, while others ("noise") will be
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
579 * decreased in the filter output. */
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
580 irange = 64.0 / range; // so irange*(max-value) is in the range [0, 63]
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
581 gain_mul = range * (fcb_type == FCB_TYPE_HARDCODED ? (5.0 / 13.0) :
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
582 (5.0 / 14.7));
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
583 angle_mul = gain_mul * (8.0 * M_LN10 / M_PI);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
584 for (n = 0; n <= 64; n++) {
12273
7323559a53fd Rename pow variable to pwr.
rbultje
parents: 12150
diff changeset
585 float pwr;
11653
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
586
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
587 idx = FFMAX(0, lrint((max - lpcs[n]) * irange) - 1);
12273
7323559a53fd Rename pow variable to pwr.
rbultje
parents: 12150
diff changeset
588 pwr = wmavoice_denoise_power_table[s->denoise_strength][idx];
7323559a53fd Rename pow variable to pwr.
rbultje
parents: 12150
diff changeset
589 lpcs[n] = angle_mul * pwr;
11653
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
590
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
591 /* 70.57 =~ 1/log10(1.0331663) */
12273
7323559a53fd Rename pow variable to pwr.
rbultje
parents: 12150
diff changeset
592 idx = (pwr * gain_mul - 0.0295) * 70.570526123;
11653
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
593 if (idx > 127) { // fallback if index falls outside table range
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
594 coeffs[n] = wmavoice_energy_table[127] *
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
595 powf(1.0331663, idx - 127);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
596 } else
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
597 coeffs[n] = wmavoice_energy_table[FFMAX(0, idx)];
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
598 }
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
599
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
600 /* calculate the Hilbert transform of the gains, which we do (since this
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
601 * is a sinus input) by doing a phase shift (in theory, H(sin())=cos()).
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
602 * Hilbert_Transform(RDFT(x)) = Laplace_Transform(x), which calculates the
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
603 * "moment" of the LPCs in this filter. */
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
604 ff_dct_calc(&s->dct, lpcs);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
605 ff_dct_calc(&s->dst, lpcs);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
606
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
607 /* Split out the coefficient indexes into phase/magnitude pairs */
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
608 idx = 255 + av_clip(lpcs[64], -255, 255);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
609 coeffs[0] = coeffs[0] * s->cos[idx];
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
610 idx = 255 + av_clip(lpcs[64] - 2 * lpcs[63], -255, 255);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
611 last_coeff = coeffs[64] * s->cos[idx];
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
612 for (n = 63;; n--) {
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
613 idx = 255 + av_clip(-lpcs[64] - 2 * lpcs[n - 1], -255, 255);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
614 coeffs[n * 2 + 1] = coeffs[n] * s->sin[idx];
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
615 coeffs[n * 2] = coeffs[n] * s->cos[idx];
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
616
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
617 if (!--n) break;
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
618
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
619 idx = 255 + av_clip( lpcs[64] - 2 * lpcs[n - 1], -255, 255);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
620 coeffs[n * 2 + 1] = coeffs[n] * s->sin[idx];
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
621 coeffs[n * 2] = coeffs[n] * s->cos[idx];
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
622 }
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
623 coeffs[1] = last_coeff;
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
624
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
625 /* move into real domain */
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
626 ff_rdft_calc(&s->irdft, coeffs);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
627
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
628 /* tilt correction and normalize scale */
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
629 memset(&coeffs[remainder], 0, sizeof(coeffs[0]) * (128 - remainder));
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
630 if (s->denoise_tilt_corr) {
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
631 float tilt_mem = 0;
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
632
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
633 coeffs[remainder - 1] = 0;
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
634 ff_tilt_compensation(&tilt_mem,
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
635 -1.8 * tilt_factor(coeffs, remainder - 1),
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
636 coeffs, remainder);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
637 }
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
638 sq = (1.0 / 64.0) * sqrtf(1 / ff_dot_productf(coeffs, coeffs, remainder));
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
639 for (n = 0; n < remainder; n++)
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
640 coeffs[n] *= sq;
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
641 }
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
642
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
643 /**
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
644 * This function applies a Wiener filter on the (noisy) speech signal as
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
645 * a means to denoise it.
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
646 *
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
647 * - take RDFT of LPCs to get the power spectrum of the noise + speech;
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
648 * - using this power spectrum, calculate (for each frequency) the Wiener
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
649 * filter gain, which depends on the frequency power and desired level
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
650 * of noise subtraction (when set too high, this leads to artifacts)
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
651 * We can do this symmetrically over the X-axis (so 0-4kHz is the inverse
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
652 * of 4-8kHz);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
653 * - by doing a phase shift, calculate the Hilbert transform of this array
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
654 * of per-frequency filter-gains to get the filtering coefficients;
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
655 * - smoothen/normalize/de-tilt these filter coefficients as desired;
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
656 * - take RDFT of noisy sound, apply the coefficients and take its IRDFT
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
657 * to get the denoised speech signal;
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
658 * - the leftover (i.e. output of the IRDFT on denoised speech data beyond
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
659 * the frame boundary) are saved and applied to subsequent frames by an
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
660 * overlap-add method (otherwise you get clicking-artifacts).
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
661 *
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
662 * @param s WMA Voice decoding context
12114
b00c10643f24 Fix two doxy warnings.
rbultje
parents: 12024
diff changeset
663 * @param fcb_type Frame (codebook) type
11653
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
664 * @param synth_pf input: the noisy speech signal, output: denoised speech
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
665 * data; should be 16-byte aligned (for ASM purposes)
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
666 * @param size size of the speech data
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
667 * @param lpcs LPCs used to synthesize this frame's speech data
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
668 */
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
669 static void wiener_denoise(WMAVoiceContext *s, int fcb_type,
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
670 float *synth_pf, int size,
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
671 const float *lpcs)
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
672 {
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
673 int remainder, lim, n;
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
674
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
675 if (fcb_type != FCB_TYPE_SILENCE) {
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
676 float *tilted_lpcs = s->tilted_lpcs_pf,
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
677 *coeffs = s->denoise_coeffs_pf, tilt_mem = 0;
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
678
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
679 tilted_lpcs[0] = 1.0;
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
680 memcpy(&tilted_lpcs[1], lpcs, sizeof(lpcs[0]) * s->lsps);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
681 memset(&tilted_lpcs[s->lsps + 1], 0,
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
682 sizeof(tilted_lpcs[0]) * (128 - s->lsps - 1));
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
683 ff_tilt_compensation(&tilt_mem, 0.7 * tilt_factor(lpcs, s->lsps),
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
684 tilted_lpcs, s->lsps + 2);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
685
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
686 /* The IRDFT output (127 samples for 7-bit filter) beyond the frame
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
687 * size is applied to the next frame. All input beyond this is zero,
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
688 * and thus all output beyond this will go towards zero, hence we can
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
689 * limit to min(size-1, 127-size) as a performance consideration. */
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
690 remainder = FFMIN(127 - size, size - 1);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
691 calc_input_response(s, tilted_lpcs, fcb_type, coeffs, remainder);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
692
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
693 /* apply coefficients (in frequency spectrum domain), i.e. complex
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
694 * number multiplication */
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
695 memset(&synth_pf[size], 0, sizeof(synth_pf[0]) * (128 - size));
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
696 ff_rdft_calc(&s->rdft, synth_pf);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
697 ff_rdft_calc(&s->rdft, coeffs);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
698 synth_pf[0] *= coeffs[0];
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
699 synth_pf[1] *= coeffs[1];
11675
5c1363c233b8 Another buffer overflow, fixes issue1758.
rbultje
parents: 11673
diff changeset
700 for (n = 1; n < 64; n++) {
11653
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
701 float v1 = synth_pf[n * 2], v2 = synth_pf[n * 2 + 1];
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
702 synth_pf[n * 2] = v1 * coeffs[n * 2] - v2 * coeffs[n * 2 + 1];
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
703 synth_pf[n * 2 + 1] = v2 * coeffs[n * 2] + v1 * coeffs[n * 2 + 1];
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
704 }
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
705 ff_rdft_calc(&s->irdft, synth_pf);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
706 }
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
707
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
708 /* merge filter output with the history of previous runs */
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
709 if (s->denoise_filter_cache_size) {
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
710 lim = FFMIN(s->denoise_filter_cache_size, size);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
711 for (n = 0; n < lim; n++)
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
712 synth_pf[n] += s->denoise_filter_cache[n];
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
713 s->denoise_filter_cache_size -= lim;
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
714 memmove(s->denoise_filter_cache, &s->denoise_filter_cache[size],
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
715 sizeof(s->denoise_filter_cache[0]) * s->denoise_filter_cache_size);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
716 }
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
717
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
718 /* move remainder of filter output into a cache for future runs */
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
719 if (fcb_type != FCB_TYPE_SILENCE) {
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
720 lim = FFMIN(remainder, s->denoise_filter_cache_size);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
721 for (n = 0; n < lim; n++)
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
722 s->denoise_filter_cache[n] += synth_pf[size + n];
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
723 if (lim < remainder) {
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
724 memcpy(&s->denoise_filter_cache[lim], &synth_pf[size + lim],
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
725 sizeof(s->denoise_filter_cache[0]) * (remainder - lim));
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
726 s->denoise_filter_cache_size = remainder;
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
727 }
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
728 }
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
729 }
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
730
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
731 /**
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
732 * Averaging projection filter, the postfilter used in WMAVoice.
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
733 *
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
734 * This uses the following steps:
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
735 * - A zero-synthesis filter (generate excitation from synth signal)
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
736 * - Kalman smoothing on excitation, based on pitch
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
737 * - Re-synthesized smoothened output
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
738 * - Iterative Wiener denoise filter
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
739 * - Adaptive gain filter
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
740 * - DC filter
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
741 *
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
742 * @param s WMAVoice decoding context
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
743 * @param synth Speech synthesis output (before postfilter)
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
744 * @param samples Output buffer for filtered samples
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
745 * @param size Buffer size of synth & samples
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
746 * @param lpcs Generated LPCs used for speech synthesis
12114
b00c10643f24 Fix two doxy warnings.
rbultje
parents: 12024
diff changeset
747 * @param zero_exc_pf destination for zero synthesis filter (16-byte aligned)
11653
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
748 * @param fcb_type Frame type (silence, hardcoded, AW-pulses or FCB-pulses)
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
749 * @param pitch Pitch of the input signal
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
750 */
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
751 static void postfilter(WMAVoiceContext *s, const float *synth,
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
752 float *samples, int size,
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
753 const float *lpcs, float *zero_exc_pf,
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
754 int fcb_type, int pitch)
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
755 {
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
756 float synth_filter_in_buf[MAX_FRAMESIZE / 2],
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
757 *synth_pf = &s->synth_filter_out_buf[MAX_LSPS_ALIGN16],
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
758 *synth_filter_in = zero_exc_pf;
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
759
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
760 assert(size <= MAX_FRAMESIZE / 2);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
761
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
762 /* generate excitation from input signal */
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
763 ff_celp_lp_zero_synthesis_filterf(zero_exc_pf, lpcs, synth, size, s->lsps);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
764
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
765 if (fcb_type >= FCB_TYPE_AW_PULSES &&
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
766 !kalman_smoothen(s, pitch, zero_exc_pf, synth_filter_in_buf, size))
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
767 synth_filter_in = synth_filter_in_buf;
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
768
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
769 /* re-synthesize speech after smoothening, and keep history */
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
770 ff_celp_lp_synthesis_filterf(synth_pf, lpcs,
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
771 synth_filter_in, size, s->lsps);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
772 memcpy(&synth_pf[-s->lsps], &synth_pf[size - s->lsps],
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
773 sizeof(synth_pf[0]) * s->lsps);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
774
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
775 wiener_denoise(s, fcb_type, synth_pf, size, lpcs);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
776
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
777 adaptive_gain_control(samples, synth_pf, synth, size, 0.99,
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
778 &s->postfilter_agc);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
779
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
780 if (s->dc_level > 8) {
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
781 /* remove ultra-low frequency DC noise / highpass filter;
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
782 * coefficients are identical to those used in SIPR decoding,
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
783 * and very closely resemble those used in AMR-NB decoding. */
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
784 ff_acelp_apply_order_2_transfer_function(samples, samples,
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
785 (const float[2]) { -1.99997, 1.0 },
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
786 (const float[2]) { -1.9330735188, 0.93589198496 },
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
787 0.93980580475, s->dcf_mem, size);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
788 }
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
789 }
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
790 /**
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
791 * @}
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
792 */
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
793
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
794 /**
11123
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
795 * Dequantize LSPs
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
796 * @param lsps output pointer to the array that will hold the LSPs
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
797 * @param num number of LSPs to be dequantized
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
798 * @param values quantized values, contains n_stages values
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
799 * @param sizes range (i.e. max value) of each quantized value
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
800 * @param n_stages number of dequantization runs
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
801 * @param table dequantization table to be used
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
802 * @param mul_q LSF multiplier
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
803 * @param base_q base (lowest) LSF values
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
804 */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
805 static void dequant_lsps(double *lsps, int num,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
806 const uint16_t *values,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
807 const uint16_t *sizes,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
808 int n_stages, const uint8_t *table,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
809 const double *mul_q,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
810 const double *base_q)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
811 {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
812 int n, m;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
813
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
814 memset(lsps, 0, num * sizeof(*lsps));
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
815 for (n = 0; n < n_stages; n++) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
816 const uint8_t *t_off = &table[values[n] * num];
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
817 double base = base_q[n], mul = mul_q[n];
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
818
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
819 for (m = 0; m < num; m++)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
820 lsps[m] += base + mul * t_off[m];
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
821
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
822 table += sizes[n] * num;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
823 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
824 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
825
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
826 /**
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
827 * @defgroup lsp_dequant LSP dequantization routines
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
828 * LSP dequantization routines, for 10/16LSPs and independent/residual coding.
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
829 * @note we assume enough bits are available, caller should check.
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
830 * lsp10i() consumes 24 bits; lsp10r() consumes an additional 24 bits;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
831 * lsp16i() consumes 34 bits; lsp16r() consumes an additional 26 bits.
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
832 * @{
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
833 */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
834 /**
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
835 * Parse 10 independently-coded LSPs.
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
836 */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
837 static void dequant_lsp10i(GetBitContext *gb, double *lsps)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
838 {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
839 static const uint16_t vec_sizes[4] = { 256, 64, 32, 32 };
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
840 static const double mul_lsf[4] = {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
841 5.2187144800e-3, 1.4626986422e-3,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
842 9.6179549166e-4, 1.1325736225e-3
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
843 };
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
844 static const double base_lsf[4] = {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
845 M_PI * -2.15522e-1, M_PI * -6.1646e-2,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
846 M_PI * -3.3486e-2, M_PI * -5.7408e-2
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
847 };
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
848 uint16_t v[4];
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
849
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
850 v[0] = get_bits(gb, 8);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
851 v[1] = get_bits(gb, 6);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
852 v[2] = get_bits(gb, 5);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
853 v[3] = get_bits(gb, 5);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
854
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
855 dequant_lsps(lsps, 10, v, vec_sizes, 4, wmavoice_dq_lsp10i,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
856 mul_lsf, base_lsf);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
857 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
858
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
859 /**
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
860 * Parse 10 independently-coded LSPs, and then derive the tables to
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
861 * generate LSPs for the other frames from them (residual coding).
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
862 */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
863 static void dequant_lsp10r(GetBitContext *gb,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
864 double *i_lsps, const double *old,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
865 double *a1, double *a2, int q_mode)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
866 {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
867 static const uint16_t vec_sizes[3] = { 128, 64, 64 };
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
868 static const double mul_lsf[3] = {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
869 2.5807601174e-3, 1.2354460219e-3, 1.1763821673e-3
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
870 };
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
871 static const double base_lsf[3] = {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
872 M_PI * -1.07448e-1, M_PI * -5.2706e-2, M_PI * -5.1634e-2
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
873 };
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
874 const float (*ipol_tab)[2][10] = q_mode ?
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
875 wmavoice_lsp10_intercoeff_b : wmavoice_lsp10_intercoeff_a;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
876 uint16_t interpol, v[3];
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
877 int n;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
878
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
879 dequant_lsp10i(gb, i_lsps);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
880
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
881 interpol = get_bits(gb, 5);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
882 v[0] = get_bits(gb, 7);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
883 v[1] = get_bits(gb, 6);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
884 v[2] = get_bits(gb, 6);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
885
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
886 for (n = 0; n < 10; n++) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
887 double delta = old[n] - i_lsps[n];
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
888 a1[n] = ipol_tab[interpol][0][n] * delta + i_lsps[n];
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
889 a1[10 + n] = ipol_tab[interpol][1][n] * delta + i_lsps[n];
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
890 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
891
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
892 dequant_lsps(a2, 20, v, vec_sizes, 3, wmavoice_dq_lsp10r,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
893 mul_lsf, base_lsf);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
894 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
895
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
896 /**
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
897 * Parse 16 independently-coded LSPs.
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
898 */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
899 static void dequant_lsp16i(GetBitContext *gb, double *lsps)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
900 {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
901 static const uint16_t vec_sizes[5] = { 256, 64, 128, 64, 128 };
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
902 static const double mul_lsf[5] = {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
903 3.3439586280e-3, 6.9908173703e-4,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
904 3.3216608306e-3, 1.0334960326e-3,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
905 3.1899104283e-3
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
906 };
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
907 static const double base_lsf[5] = {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
908 M_PI * -1.27576e-1, M_PI * -2.4292e-2,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
909 M_PI * -1.28094e-1, M_PI * -3.2128e-2,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
910 M_PI * -1.29816e-1
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
911 };
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
912 uint16_t v[5];
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
913
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
914 v[0] = get_bits(gb, 8);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
915 v[1] = get_bits(gb, 6);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
916 v[2] = get_bits(gb, 7);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
917 v[3] = get_bits(gb, 6);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
918 v[4] = get_bits(gb, 7);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
919
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
920 dequant_lsps( lsps, 5, v, vec_sizes, 2,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
921 wmavoice_dq_lsp16i1, mul_lsf, base_lsf);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
922 dequant_lsps(&lsps[5], 5, &v[2], &vec_sizes[2], 2,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
923 wmavoice_dq_lsp16i2, &mul_lsf[2], &base_lsf[2]);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
924 dequant_lsps(&lsps[10], 6, &v[4], &vec_sizes[4], 1,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
925 wmavoice_dq_lsp16i3, &mul_lsf[4], &base_lsf[4]);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
926 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
927
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
928 /**
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
929 * Parse 16 independently-coded LSPs, and then derive the tables to
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
930 * generate LSPs for the other frames from them (residual coding).
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
931 */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
932 static void dequant_lsp16r(GetBitContext *gb,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
933 double *i_lsps, const double *old,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
934 double *a1, double *a2, int q_mode)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
935 {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
936 static const uint16_t vec_sizes[3] = { 128, 128, 128 };
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
937 static const double mul_lsf[3] = {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
938 1.2232979501e-3, 1.4062241527e-3, 1.6114744851e-3
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
939 };
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
940 static const double base_lsf[3] = {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
941 M_PI * -5.5830e-2, M_PI * -5.2908e-2, M_PI * -5.4776e-2
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
942 };
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
943 const float (*ipol_tab)[2][16] = q_mode ?
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
944 wmavoice_lsp16_intercoeff_b : wmavoice_lsp16_intercoeff_a;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
945 uint16_t interpol, v[3];
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
946 int n;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
947
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
948 dequant_lsp16i(gb, i_lsps);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
949
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
950 interpol = get_bits(gb, 5);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
951 v[0] = get_bits(gb, 7);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
952 v[1] = get_bits(gb, 7);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
953 v[2] = get_bits(gb, 7);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
954
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
955 for (n = 0; n < 16; n++) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
956 double delta = old[n] - i_lsps[n];
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
957 a1[n] = ipol_tab[interpol][0][n] * delta + i_lsps[n];
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
958 a1[16 + n] = ipol_tab[interpol][1][n] * delta + i_lsps[n];
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
959 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
960
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
961 dequant_lsps( a2, 10, v, vec_sizes, 1,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
962 wmavoice_dq_lsp16r1, mul_lsf, base_lsf);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
963 dequant_lsps(&a2[10], 10, &v[1], &vec_sizes[1], 1,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
964 wmavoice_dq_lsp16r2, &mul_lsf[1], &base_lsf[1]);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
965 dequant_lsps(&a2[20], 12, &v[2], &vec_sizes[2], 1,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
966 wmavoice_dq_lsp16r3, &mul_lsf[2], &base_lsf[2]);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
967 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
968
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
969 /**
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
970 * @}
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
971 * @defgroup aw Pitch-adaptive window coding functions
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
972 * The next few functions are for pitch-adaptive window coding.
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
973 * @{
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
974 */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
975 /**
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
976 * Parse the offset of the first pitch-adaptive window pulses, and
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
977 * the distribution of pulses between the two blocks in this frame.
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
978 * @param s WMA Voice decoding context private data
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
979 * @param gb bit I/O context
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
980 * @param pitch pitch for each block in this frame
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
981 */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
982 static void aw_parse_coords(WMAVoiceContext *s, GetBitContext *gb,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
983 const int *pitch)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
984 {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
985 static const int16_t start_offset[94] = {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
986 -11, -9, -7, -5, -3, -1, 1, 3, 5, 7, 9, 11,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
987 13, 15, 18, 17, 19, 20, 21, 22, 23, 24, 25, 26,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
988 27, 28, 29, 30, 31, 32, 33, 35, 37, 39, 41, 43,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
989 45, 47, 49, 51, 53, 55, 57, 59, 61, 63, 65, 67,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
990 69, 71, 73, 75, 77, 79, 81, 83, 85, 87, 89, 91,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
991 93, 95, 97, 99, 101, 103, 105, 107, 109, 111, 113, 115,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
992 117, 119, 121, 123, 125, 127, 129, 131, 133, 135, 137, 139,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
993 141, 143, 145, 147, 149, 151, 153, 155, 157, 159
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
994 };
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
995 int bits, offset;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
996
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
997 /* position of pulse */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
998 s->aw_idx_is_ext = 0;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
999 if ((bits = get_bits(gb, 6)) >= 54) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1000 s->aw_idx_is_ext = 1;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1001 bits += (bits - 54) * 3 + get_bits(gb, 2);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1002 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1003
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1004 /* for a repeated pulse at pulse_off with a pitch_lag of pitch[], count
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1005 * the distribution of the pulses in each block contained in this frame. */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1006 s->aw_pulse_range = FFMIN(pitch[0], pitch[1]) > 32 ? 24 : 16;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1007 for (offset = start_offset[bits]; offset < 0; offset += pitch[0]) ;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1008 s->aw_n_pulses[0] = (pitch[0] - 1 + MAX_FRAMESIZE / 2 - offset) / pitch[0];
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1009 s->aw_first_pulse_off[0] = offset - s->aw_pulse_range / 2;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1010 offset += s->aw_n_pulses[0] * pitch[0];
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1011 s->aw_n_pulses[1] = (pitch[1] - 1 + MAX_FRAMESIZE - offset) / pitch[1];
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1012 s->aw_first_pulse_off[1] = offset - (MAX_FRAMESIZE + s->aw_pulse_range) / 2;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1013
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1014 /* if continuing from a position before the block, reset position to
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1015 * start of block (when corrected for the range over which it can be
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1016 * spread in aw_pulse_set1()). */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1017 if (start_offset[bits] < MAX_FRAMESIZE / 2) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1018 while (s->aw_first_pulse_off[1] - pitch[1] + s->aw_pulse_range > 0)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1019 s->aw_first_pulse_off[1] -= pitch[1];
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1020 if (start_offset[bits] < 0)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1021 while (s->aw_first_pulse_off[0] - pitch[0] + s->aw_pulse_range > 0)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1022 s->aw_first_pulse_off[0] -= pitch[0];
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1023 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1024 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1025
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1026 /**
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1027 * Apply second set of pitch-adaptive window pulses.
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1028 * @param s WMA Voice decoding context private data
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1029 * @param gb bit I/O context
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1030 * @param block_idx block index in frame [0, 1]
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1031 * @param fcb structure containing fixed codebook vector info
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1032 */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1033 static void aw_pulse_set2(WMAVoiceContext *s, GetBitContext *gb,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1034 int block_idx, AMRFixed *fcb)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1035 {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1036 uint16_t use_mask[7]; // only 5 are used, rest is padding
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1037 /* in this function, idx is the index in the 80-bit (+ padding) use_mask
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1038 * bit-array. Since use_mask consists of 16-bit values, the lower 4 bits
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1039 * of idx are the position of the bit within a particular item in the
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1040 * array (0 being the most significant bit, and 15 being the least
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1041 * significant bit), and the remainder (>> 4) is the index in the
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1042 * use_mask[]-array. This is faster and uses less memory than using a
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1043 * 80-byte/80-int array. */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1044 int pulse_off = s->aw_first_pulse_off[block_idx],
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1045 pulse_start, n, idx, range, aidx, start_off = 0;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1046
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1047 /* set offset of first pulse to within this block */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1048 if (s->aw_n_pulses[block_idx] > 0)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1049 while (pulse_off + s->aw_pulse_range < 1)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1050 pulse_off += fcb->pitch_lag;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1051
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1052 /* find range per pulse */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1053 if (s->aw_n_pulses[0] > 0) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1054 if (block_idx == 0) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1055 range = 32;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1056 } else /* block_idx = 1 */ {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1057 range = 8;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1058 if (s->aw_n_pulses[block_idx] > 0)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1059 pulse_off = s->aw_next_pulse_off_cache;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1060 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1061 } else
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1062 range = 16;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1063 pulse_start = s->aw_n_pulses[block_idx] > 0 ? pulse_off - range / 2 : 0;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1064
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1065 /* aw_pulse_set1() already applies pulses around pulse_off (to be exactly,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1066 * in the range of [pulse_off, pulse_off + s->aw_pulse_range], and thus
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1067 * we exclude that range from being pulsed again in this function. */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1068 memset( use_mask, -1, 5 * sizeof(use_mask[0]));
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1069 memset(&use_mask[5], 0, 2 * sizeof(use_mask[0]));
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1070 if (s->aw_n_pulses[block_idx] > 0)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1071 for (idx = pulse_off; idx < MAX_FRAMESIZE / 2; idx += fcb->pitch_lag) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1072 int excl_range = s->aw_pulse_range; // always 16 or 24
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1073 uint16_t *use_mask_ptr = &use_mask[idx >> 4];
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1074 int first_sh = 16 - (idx & 15);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1075 *use_mask_ptr++ &= 0xFFFF << first_sh;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1076 excl_range -= first_sh;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1077 if (excl_range >= 16) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1078 *use_mask_ptr++ = 0;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1079 *use_mask_ptr &= 0xFFFF >> (excl_range - 16);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1080 } else
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1081 *use_mask_ptr &= 0xFFFF >> excl_range;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1082 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1083
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1084 /* find the 'aidx'th offset that is not excluded */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1085 aidx = get_bits(gb, s->aw_n_pulses[0] > 0 ? 5 - 2 * block_idx : 4);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1086 for (n = 0; n <= aidx; pulse_start++) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1087 for (idx = pulse_start; idx < 0; idx += fcb->pitch_lag) ;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1088 if (idx >= MAX_FRAMESIZE / 2) { // find from zero
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1089 if (use_mask[0]) idx = 0x0F;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1090 else if (use_mask[1]) idx = 0x1F;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1091 else if (use_mask[2]) idx = 0x2F;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1092 else if (use_mask[3]) idx = 0x3F;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1093 else if (use_mask[4]) idx = 0x4F;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1094 else return;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1095 idx -= av_log2_16bit(use_mask[idx >> 4]);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1096 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1097 if (use_mask[idx >> 4] & (0x8000 >> (idx & 15))) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1098 use_mask[idx >> 4] &= ~(0x8000 >> (idx & 15));
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1099 n++;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1100 start_off = idx;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1101 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1102 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1103
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1104 fcb->x[fcb->n] = start_off;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1105 fcb->y[fcb->n] = get_bits1(gb) ? -1.0 : 1.0;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1106 fcb->n++;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1107
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1108 /* set offset for next block, relative to start of that block */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1109 n = (MAX_FRAMESIZE / 2 - start_off) % fcb->pitch_lag;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1110 s->aw_next_pulse_off_cache = n ? fcb->pitch_lag - n : 0;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1111 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1112
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1113 /**
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1114 * Apply first set of pitch-adaptive window pulses.
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1115 * @param s WMA Voice decoding context private data
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1116 * @param gb bit I/O context
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1117 * @param block_idx block index in frame [0, 1]
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1118 * @param fcb storage location for fixed codebook pulse info
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1119 */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1120 static void aw_pulse_set1(WMAVoiceContext *s, GetBitContext *gb,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1121 int block_idx, AMRFixed *fcb)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1122 {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1123 int val = get_bits(gb, 12 - 2 * (s->aw_idx_is_ext && !block_idx));
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1124 float v;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1125
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1126 if (s->aw_n_pulses[block_idx] > 0) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1127 int n, v_mask, i_mask, sh, n_pulses;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1128
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1129 if (s->aw_pulse_range == 24) { // 3 pulses, 1:sign + 3:index each
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1130 n_pulses = 3;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1131 v_mask = 8;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1132 i_mask = 7;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1133 sh = 4;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1134 } else { // 4 pulses, 1:sign + 2:index each
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1135 n_pulses = 4;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1136 v_mask = 4;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1137 i_mask = 3;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1138 sh = 3;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1139 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1140
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1141 for (n = n_pulses - 1; n >= 0; n--, val >>= sh) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1142 fcb->y[fcb->n] = (val & v_mask) ? -1.0 : 1.0;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1143 fcb->x[fcb->n] = (val & i_mask) * n_pulses + n +
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1144 s->aw_first_pulse_off[block_idx];
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1145 while (fcb->x[fcb->n] < 0)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1146 fcb->x[fcb->n] += fcb->pitch_lag;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1147 if (fcb->x[fcb->n] < MAX_FRAMESIZE / 2)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1148 fcb->n++;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1149 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1150 } else {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1151 int num2 = (val & 0x1FF) >> 1, delta, idx;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1152
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1153 if (num2 < 1 * 79) { delta = 1; idx = num2 + 1; }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1154 else if (num2 < 2 * 78) { delta = 3; idx = num2 + 1 - 1 * 77; }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1155 else if (num2 < 3 * 77) { delta = 5; idx = num2 + 1 - 2 * 76; }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1156 else { delta = 7; idx = num2 + 1 - 3 * 75; }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1157 v = (val & 0x200) ? -1.0 : 1.0;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1158
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1159 fcb->no_repeat_mask |= 3 << fcb->n;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1160 fcb->x[fcb->n] = idx - delta;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1161 fcb->y[fcb->n] = v;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1162 fcb->x[fcb->n + 1] = idx;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1163 fcb->y[fcb->n + 1] = (val & 1) ? -v : v;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1164 fcb->n += 2;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1165 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1166 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1167
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1168 /**
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1169 * @}
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1170 *
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1171 * Generate a random number from frame_cntr and block_idx, which will lief
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1172 * in the range [0, 1000 - block_size] (so it can be used as an index in a
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1173 * table of size 1000 of which you want to read block_size entries).
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1174 *
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1175 * @param frame_cntr current frame number
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1176 * @param block_num current block index
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1177 * @param block_size amount of entries we want to read from a table
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1178 * that has 1000 entries
11556
b94e1810ce4c Replace @returns by @return.
benoit
parents: 11123
diff changeset
1179 * @return a (non-)random number in the [0, 1000 - block_size] range.
11123
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1180 */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1181 static int pRNG(int frame_cntr, int block_num, int block_size)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1182 {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1183 /* array to simplify the calculation of z:
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1184 * y = (x % 9) * 5 + 6;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1185 * z = (49995 * x) / y;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1186 * Since y only has 9 values, we can remove the division by using a
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1187 * LUT and using FASTDIV-style divisions. For each of the 9 values
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1188 * of y, we can rewrite z as:
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1189 * z = x * (49995 / y) + x * ((49995 % y) / y)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1190 * In this table, each col represents one possible value of y, the
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1191 * first number is 49995 / y, and the second is the FASTDIV variant
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1192 * of 49995 % y / y. */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1193 static const unsigned int div_tbl[9][2] = {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1194 { 8332, 3 * 715827883U }, // y = 6
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1195 { 4545, 0 * 390451573U }, // y = 11
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1196 { 3124, 11 * 268435456U }, // y = 16
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1197 { 2380, 15 * 204522253U }, // y = 21
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1198 { 1922, 23 * 165191050U }, // y = 26
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1199 { 1612, 23 * 138547333U }, // y = 31
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1200 { 1388, 27 * 119304648U }, // y = 36
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1201 { 1219, 16 * 104755300U }, // y = 41
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1202 { 1086, 39 * 93368855U } // y = 46
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1203 };
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1204 unsigned int z, y, x = MUL16(block_num, 1877) + frame_cntr;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1205 if (x >= 0xFFFF) x -= 0xFFFF; // max value of x is 8*1877+0xFFFE=0x13AA6,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1206 // so this is effectively a modulo (%)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1207 y = x - 9 * MULH(477218589, x); // x % 9
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1208 z = (uint16_t) (x * div_tbl[y][0] + UMULH(x, div_tbl[y][1]));
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1209 // z = x * 49995 / (y * 5 + 6)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1210 return z % (1000 - block_size);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1211 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1212
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1213 /**
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1214 * Parse hardcoded signal for a single block.
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1215 * @note see #synth_block().
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1216 */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1217 static void synth_block_hardcoded(WMAVoiceContext *s, GetBitContext *gb,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1218 int block_idx, int size,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1219 const struct frame_type_desc *frame_desc,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1220 float *excitation)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1221 {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1222 float gain;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1223 int n, r_idx;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1224
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1225 assert(size <= MAX_FRAMESIZE);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1226
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1227 /* Set the offset from which we start reading wmavoice_std_codebook */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1228 if (frame_desc->fcb_type == FCB_TYPE_SILENCE) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1229 r_idx = pRNG(s->frame_cntr, block_idx, size);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1230 gain = s->silence_gain;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1231 } else /* FCB_TYPE_HARDCODED */ {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1232 r_idx = get_bits(gb, 8);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1233 gain = wmavoice_gain_universal[get_bits(gb, 6)];
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1234 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1235
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1236 /* Clear gain prediction parameters */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1237 memset(s->gain_pred_err, 0, sizeof(s->gain_pred_err));
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1238
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1239 /* Apply gain to hardcoded codebook and use that as excitation signal */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1240 for (n = 0; n < size; n++)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1241 excitation[n] = wmavoice_std_codebook[r_idx + n] * gain;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1242 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1243
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1244 /**
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1245 * Parse FCB/ACB signal for a single block.
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1246 * @note see #synth_block().
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1247 */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1248 static void synth_block_fcb_acb(WMAVoiceContext *s, GetBitContext *gb,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1249 int block_idx, int size,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1250 int block_pitch_sh2,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1251 const struct frame_type_desc *frame_desc,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1252 float *excitation)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1253 {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1254 static const float gain_coeff[6] = {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1255 0.8169, -0.06545, 0.1726, 0.0185, -0.0359, 0.0458
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1256 };
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1257 float pulses[MAX_FRAMESIZE / 2], pred_err, acb_gain, fcb_gain;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1258 int n, idx, gain_weight;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1259 AMRFixed fcb;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1260
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1261 assert(size <= MAX_FRAMESIZE / 2);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1262 memset(pulses, 0, sizeof(*pulses) * size);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1263
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1264 fcb.pitch_lag = block_pitch_sh2 >> 2;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1265 fcb.pitch_fac = 1.0;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1266 fcb.no_repeat_mask = 0;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1267 fcb.n = 0;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1268
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1269 /* For the other frame types, this is where we apply the innovation
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1270 * (fixed) codebook pulses of the speech signal. */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1271 if (frame_desc->fcb_type == FCB_TYPE_AW_PULSES) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1272 aw_pulse_set1(s, gb, block_idx, &fcb);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1273 aw_pulse_set2(s, gb, block_idx, &fcb);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1274 } else /* FCB_TYPE_EXC_PULSES */ {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1275 int offset_nbits = 5 - frame_desc->log_n_blocks;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1276
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1277 fcb.no_repeat_mask = -1;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1278 /* similar to ff_decode_10_pulses_35bits(), but with single pulses
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1279 * (instead of double) for a subset of pulses */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1280 for (n = 0; n < 5; n++) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1281 float sign;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1282 int pos1, pos2;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1283
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1284 sign = get_bits1(gb) ? 1.0 : -1.0;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1285 pos1 = get_bits(gb, offset_nbits);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1286 fcb.x[fcb.n] = n + 5 * pos1;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1287 fcb.y[fcb.n++] = sign;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1288 if (n < frame_desc->dbl_pulses) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1289 pos2 = get_bits(gb, offset_nbits);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1290 fcb.x[fcb.n] = n + 5 * pos2;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1291 fcb.y[fcb.n++] = (pos1 < pos2) ? -sign : sign;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1292 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1293 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1294 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1295 ff_set_fixed_vector(pulses, &fcb, 1.0, size);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1296
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1297 /* Calculate gain for adaptive & fixed codebook signal.
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1298 * see ff_amr_set_fixed_gain(). */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1299 idx = get_bits(gb, 7);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1300 fcb_gain = expf(ff_dot_productf(s->gain_pred_err, gain_coeff, 6) -
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1301 5.2409161640 + wmavoice_gain_codebook_fcb[idx]);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1302 acb_gain = wmavoice_gain_codebook_acb[idx];
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1303 pred_err = av_clipf(wmavoice_gain_codebook_fcb[idx],
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1304 -2.9957322736 /* log(0.05) */,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1305 1.6094379124 /* log(5.0) */);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1306
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1307 gain_weight = 8 >> frame_desc->log_n_blocks;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1308 memmove(&s->gain_pred_err[gain_weight], s->gain_pred_err,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1309 sizeof(*s->gain_pred_err) * (6 - gain_weight));
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1310 for (n = 0; n < gain_weight; n++)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1311 s->gain_pred_err[n] = pred_err;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1312
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1313 /* Calculation of adaptive codebook */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1314 if (frame_desc->acb_type == ACB_TYPE_ASYMMETRIC) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1315 int len;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1316 for (n = 0; n < size; n += len) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1317 int next_idx_sh16;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1318 int abs_idx = block_idx * size + n;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1319 int pitch_sh16 = (s->last_pitch_val << 16) +
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1320 s->pitch_diff_sh16 * abs_idx;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1321 int pitch = (pitch_sh16 + 0x6FFF) >> 16;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1322 int idx_sh16 = ((pitch << 16) - pitch_sh16) * 8 + 0x58000;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1323 idx = idx_sh16 >> 16;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1324 if (s->pitch_diff_sh16) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1325 if (s->pitch_diff_sh16 > 0) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1326 next_idx_sh16 = (idx_sh16) &~ 0xFFFF;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1327 } else
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1328 next_idx_sh16 = (idx_sh16 + 0x10000) &~ 0xFFFF;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1329 len = av_clip((idx_sh16 - next_idx_sh16) / s->pitch_diff_sh16 / 8,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1330 1, size - n);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1331 } else
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1332 len = size;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1333
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1334 ff_acelp_interpolatef(&excitation[n], &excitation[n - pitch],
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1335 wmavoice_ipol1_coeffs, 17,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1336 idx, 9, len);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1337 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1338 } else /* ACB_TYPE_HAMMING */ {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1339 int block_pitch = block_pitch_sh2 >> 2;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1340 idx = block_pitch_sh2 & 3;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1341 if (idx) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1342 ff_acelp_interpolatef(excitation, &excitation[-block_pitch],
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1343 wmavoice_ipol2_coeffs, 4,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1344 idx, 8, size);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1345 } else
12021
18c23a632001 Silense one warning:
cehoyos
parents: 11675
diff changeset
1346 av_memcpy_backptr((uint8_t *) excitation, sizeof(float) * block_pitch,
11123
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1347 sizeof(float) * size);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1348 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1349
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1350 /* Interpolate ACB/FCB and use as excitation signal */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1351 ff_weighted_vector_sumf(excitation, excitation, pulses,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1352 acb_gain, fcb_gain, size);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1353 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1354
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1355 /**
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1356 * Parse data in a single block.
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1357 * @note we assume enough bits are available, caller should check.
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1358 *
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1359 * @param s WMA Voice decoding context private data
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1360 * @param gb bit I/O context
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1361 * @param block_idx index of the to-be-read block
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1362 * @param size amount of samples to be read in this block
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1363 * @param block_pitch_sh2 pitch for this block << 2
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1364 * @param lsps LSPs for (the end of) this frame
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1365 * @param prev_lsps LSPs for the last frame
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1366 * @param frame_desc frame type descriptor
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1367 * @param excitation target memory for the ACB+FCB interpolated signal
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1368 * @param synth target memory for the speech synthesis filter output
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1369 * @return 0 on success, <0 on error.
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1370 */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1371 static void synth_block(WMAVoiceContext *s, GetBitContext *gb,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1372 int block_idx, int size,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1373 int block_pitch_sh2,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1374 const double *lsps, const double *prev_lsps,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1375 const struct frame_type_desc *frame_desc,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1376 float *excitation, float *synth)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1377 {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1378 double i_lsps[MAX_LSPS];
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1379 float lpcs[MAX_LSPS];
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1380 float fac;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1381 int n;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1382
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1383 if (frame_desc->acb_type == ACB_TYPE_NONE)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1384 synth_block_hardcoded(s, gb, block_idx, size, frame_desc, excitation);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1385 else
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1386 synth_block_fcb_acb(s, gb, block_idx, size, block_pitch_sh2,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1387 frame_desc, excitation);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1388
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1389 /* convert interpolated LSPs to LPCs */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1390 fac = (block_idx + 0.5) / frame_desc->n_blocks;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1391 for (n = 0; n < s->lsps; n++) // LSF -> LSP
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1392 i_lsps[n] = cos(prev_lsps[n] + fac * (lsps[n] - prev_lsps[n]));
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1393 ff_acelp_lspd2lpc(i_lsps, lpcs, s->lsps >> 1);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1394
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1395 /* Speech synthesis */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1396 ff_celp_lp_synthesis_filterf(synth, lpcs, excitation, size, s->lsps);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1397 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1398
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1399 /**
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1400 * Synthesize output samples for a single frame.
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1401 * @note we assume enough bits are available, caller should check.
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1402 *
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1403 * @param ctx WMA Voice decoder context
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1404 * @param gb bit I/O context (s->gb or one for cross-packet superframes)
11653
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
1405 * @param frame_idx Frame number within superframe [0-2]
11123
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1406 * @param samples pointer to output sample buffer, has space for at least 160
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1407 * samples
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1408 * @param lsps LSP array
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1409 * @param prev_lsps array of previous frame's LSPs
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1410 * @param excitation target buffer for excitation signal
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1411 * @param synth target buffer for synthesized speech data
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1412 * @return 0 on success, <0 on error.
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1413 */
11653
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
1414 static int synth_frame(AVCodecContext *ctx, GetBitContext *gb, int frame_idx,
11123
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1415 float *samples,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1416 const double *lsps, const double *prev_lsps,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1417 float *excitation, float *synth)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1418 {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1419 WMAVoiceContext *s = ctx->priv_data;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1420 int n, n_blocks_x2, log_n_blocks_x2, cur_pitch_val;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1421 int pitch[MAX_BLOCKS], last_block_pitch;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1422
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1423 /* Parse frame type ("frame header"), see frame_descs */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1424 int bd_idx = s->vbm_tree[get_vlc2(gb, frame_type_vlc.table, 6, 3)],
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1425 block_nsamples = MAX_FRAMESIZE / frame_descs[bd_idx].n_blocks;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1426
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1427 if (bd_idx < 0) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1428 av_log(ctx, AV_LOG_ERROR,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1429 "Invalid frame type VLC code, skipping\n");
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1430 return -1;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1431 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1432
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1433 /* Pitch calculation for ACB_TYPE_ASYMMETRIC ("pitch-per-frame") */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1434 if (frame_descs[bd_idx].acb_type == ACB_TYPE_ASYMMETRIC) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1435 /* Pitch is provided per frame, which is interpreted as the pitch of
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1436 * the last sample of the last block of this frame. We can interpolate
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1437 * the pitch of other blocks (and even pitch-per-sample) by gradually
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1438 * incrementing/decrementing prev_frame_pitch to cur_pitch_val. */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1439 n_blocks_x2 = frame_descs[bd_idx].n_blocks << 1;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1440 log_n_blocks_x2 = frame_descs[bd_idx].log_n_blocks + 1;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1441 cur_pitch_val = s->min_pitch_val + get_bits(gb, s->pitch_nbits);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1442 cur_pitch_val = FFMIN(cur_pitch_val, s->max_pitch_val - 1);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1443 if (s->last_acb_type == ACB_TYPE_NONE ||
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1444 20 * abs(cur_pitch_val - s->last_pitch_val) >
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1445 (cur_pitch_val + s->last_pitch_val))
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1446 s->last_pitch_val = cur_pitch_val;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1447
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1448 /* pitch per block */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1449 for (n = 0; n < frame_descs[bd_idx].n_blocks; n++) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1450 int fac = n * 2 + 1;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1451
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1452 pitch[n] = (MUL16(fac, cur_pitch_val) +
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1453 MUL16((n_blocks_x2 - fac), s->last_pitch_val) +
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1454 frame_descs[bd_idx].n_blocks) >> log_n_blocks_x2;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1455 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1456
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1457 /* "pitch-diff-per-sample" for calculation of pitch per sample */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1458 s->pitch_diff_sh16 =
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1459 ((cur_pitch_val - s->last_pitch_val) << 16) / MAX_FRAMESIZE;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1460 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1461
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1462 /* Global gain (if silence) and pitch-adaptive window coordinates */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1463 switch (frame_descs[bd_idx].fcb_type) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1464 case FCB_TYPE_SILENCE:
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1465 s->silence_gain = wmavoice_gain_silence[get_bits(gb, 8)];
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1466 break;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1467 case FCB_TYPE_AW_PULSES:
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1468 aw_parse_coords(s, gb, pitch);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1469 break;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1470 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1471
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1472 for (n = 0; n < frame_descs[bd_idx].n_blocks; n++) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1473 int bl_pitch_sh2;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1474
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1475 /* Pitch calculation for ACB_TYPE_HAMMING ("pitch-per-block") */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1476 switch (frame_descs[bd_idx].acb_type) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1477 case ACB_TYPE_HAMMING: {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1478 /* Pitch is given per block. Per-block pitches are encoded as an
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1479 * absolute value for the first block, and then delta values
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1480 * relative to this value) for all subsequent blocks. The scale of
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1481 * this pitch value is semi-logaritmic compared to its use in the
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1482 * decoder, so we convert it to normal scale also. */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1483 int block_pitch,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1484 t1 = (s->block_conv_table[1] - s->block_conv_table[0]) << 2,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1485 t2 = (s->block_conv_table[2] - s->block_conv_table[1]) << 1,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1486 t3 = s->block_conv_table[3] - s->block_conv_table[2] + 1;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1487
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1488 if (n == 0) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1489 block_pitch = get_bits(gb, s->block_pitch_nbits);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1490 } else
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1491 block_pitch = last_block_pitch - s->block_delta_pitch_hrange +
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1492 get_bits(gb, s->block_delta_pitch_nbits);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1493 /* Convert last_ so that any next delta is within _range */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1494 last_block_pitch = av_clip(block_pitch,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1495 s->block_delta_pitch_hrange,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1496 s->block_pitch_range -
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1497 s->block_delta_pitch_hrange);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1498
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1499 /* Convert semi-log-style scale back to normal scale */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1500 if (block_pitch < t1) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1501 bl_pitch_sh2 = (s->block_conv_table[0] << 2) + block_pitch;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1502 } else {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1503 block_pitch -= t1;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1504 if (block_pitch < t2) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1505 bl_pitch_sh2 =
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1506 (s->block_conv_table[1] << 2) + (block_pitch << 1);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1507 } else {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1508 block_pitch -= t2;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1509 if (block_pitch < t3) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1510 bl_pitch_sh2 =
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1511 (s->block_conv_table[2] + block_pitch) << 2;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1512 } else
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1513 bl_pitch_sh2 = s->block_conv_table[3] << 2;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1514 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1515 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1516 pitch[n] = bl_pitch_sh2 >> 2;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1517 break;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1518 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1519
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1520 case ACB_TYPE_ASYMMETRIC: {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1521 bl_pitch_sh2 = pitch[n] << 2;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1522 break;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1523 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1524
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1525 default: // ACB_TYPE_NONE has no pitch
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1526 bl_pitch_sh2 = 0;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1527 break;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1528 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1529
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1530 synth_block(s, gb, n, block_nsamples, bl_pitch_sh2,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1531 lsps, prev_lsps, &frame_descs[bd_idx],
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1532 &excitation[n * block_nsamples],
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1533 &synth[n * block_nsamples]);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1534 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1535
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1536 /* Averaging projection filter, if applicable. Else, just copy samples
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1537 * from synthesis buffer */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1538 if (s->do_apf) {
11653
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
1539 double i_lsps[MAX_LSPS];
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
1540 float lpcs[MAX_LSPS];
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
1541
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
1542 for (n = 0; n < s->lsps; n++) // LSF -> LSP
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
1543 i_lsps[n] = cos(0.5 * (prev_lsps[n] + lsps[n]));
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
1544 ff_acelp_lspd2lpc(i_lsps, lpcs, s->lsps >> 1);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
1545 postfilter(s, synth, samples, 80, lpcs,
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
1546 &s->zero_exc_pf[s->history_nsamples + MAX_FRAMESIZE * frame_idx],
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
1547 frame_descs[bd_idx].fcb_type, pitch[0]);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
1548
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
1549 for (n = 0; n < s->lsps; n++) // LSF -> LSP
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
1550 i_lsps[n] = cos(lsps[n]);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
1551 ff_acelp_lspd2lpc(i_lsps, lpcs, s->lsps >> 1);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
1552 postfilter(s, &synth[80], &samples[80], 80, lpcs,
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
1553 &s->zero_exc_pf[s->history_nsamples + MAX_FRAMESIZE * frame_idx + 80],
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
1554 frame_descs[bd_idx].fcb_type, pitch[0]);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
1555 } else
11652
8b6f3d3b55cb Move clipping of audio samples (for those codecs outputting float) from decoder
rbultje
parents: 11644
diff changeset
1556 memcpy(samples, synth, 160 * sizeof(synth[0]));
11123
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1557
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1558 /* Cache values for next frame */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1559 s->frame_cntr++;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1560 if (s->frame_cntr >= 0xFFFF) s->frame_cntr -= 0xFFFF; // i.e. modulo (%)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1561 s->last_acb_type = frame_descs[bd_idx].acb_type;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1562 switch (frame_descs[bd_idx].acb_type) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1563 case ACB_TYPE_NONE:
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1564 s->last_pitch_val = 0;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1565 break;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1566 case ACB_TYPE_ASYMMETRIC:
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1567 s->last_pitch_val = cur_pitch_val;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1568 break;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1569 case ACB_TYPE_HAMMING:
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1570 s->last_pitch_val = pitch[frame_descs[bd_idx].n_blocks - 1];
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1571 break;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1572 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1573
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1574 return 0;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1575 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1576
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1577 /**
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1578 * Ensure minimum value for first item, maximum value for last value,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1579 * proper spacing between each value and proper ordering.
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1580 *
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1581 * @param lsps array of LSPs
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1582 * @param num size of LSP array
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1583 *
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1584 * @note basically a double version of #ff_acelp_reorder_lsf(), might be
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1585 * useful to put in a generic location later on. Parts are also
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1586 * present in #ff_set_min_dist_lsf() + #ff_sort_nearly_sorted_floats(),
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1587 * which is in float.
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1588 */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1589 static void stabilize_lsps(double *lsps, int num)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1590 {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1591 int n, m, l;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1592
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1593 /* set minimum value for first, maximum value for last and minimum
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1594 * spacing between LSF values.
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1595 * Very similar to ff_set_min_dist_lsf(), but in double. */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1596 lsps[0] = FFMAX(lsps[0], 0.0015 * M_PI);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1597 for (n = 1; n < num; n++)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1598 lsps[n] = FFMAX(lsps[n], lsps[n - 1] + 0.0125 * M_PI);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1599 lsps[num - 1] = FFMIN(lsps[num - 1], 0.9985 * M_PI);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1600
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1601 /* reorder (looks like one-time / non-recursed bubblesort).
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1602 * Very similar to ff_sort_nearly_sorted_floats(), but in double. */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1603 for (n = 1; n < num; n++) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1604 if (lsps[n] < lsps[n - 1]) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1605 for (m = 1; m < num; m++) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1606 double tmp = lsps[m];
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1607 for (l = m - 1; l >= 0; l--) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1608 if (lsps[l] <= tmp) break;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1609 lsps[l + 1] = lsps[l];
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1610 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1611 lsps[l + 1] = tmp;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1612 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1613 break;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1614 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1615 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1616 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1617
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1618 /**
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1619 * Test if there's enough bits to read 1 superframe.
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1620 *
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1621 * @param orig_gb bit I/O context used for reading. This function
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1622 * does not modify the state of the bitreader; it
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1623 * only uses it to copy the current stream position
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1624 * @param s WMA Voice decoding context private data
11556
b94e1810ce4c Replace @returns by @return.
benoit
parents: 11123
diff changeset
1625 * @return -1 if unsupported, 1 on not enough bits or 0 if OK.
11123
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1626 */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1627 static int check_bits_for_superframe(GetBitContext *orig_gb,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1628 WMAVoiceContext *s)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1629 {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1630 GetBitContext s_gb, *gb = &s_gb;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1631 int n, need_bits, bd_idx;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1632 const struct frame_type_desc *frame_desc;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1633
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1634 /* initialize a copy */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1635 init_get_bits(gb, orig_gb->buffer, orig_gb->size_in_bits);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1636 skip_bits_long(gb, get_bits_count(orig_gb));
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1637 assert(get_bits_left(gb) == get_bits_left(orig_gb));
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1638
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1639 /* superframe header */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1640 if (get_bits_left(gb) < 14)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1641 return 1;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1642 if (!get_bits1(gb))
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1643 return -1; // WMAPro-in-WMAVoice superframe
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1644 if (get_bits1(gb)) skip_bits(gb, 12); // number of samples in superframe
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1645 if (s->has_residual_lsps) { // residual LSPs (for all frames)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1646 if (get_bits_left(gb) < s->sframe_lsp_bitsize)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1647 return 1;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1648 skip_bits_long(gb, s->sframe_lsp_bitsize);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1649 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1650
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1651 /* frames */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1652 for (n = 0; n < MAX_FRAMES; n++) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1653 int aw_idx_is_ext = 0;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1654
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1655 if (!s->has_residual_lsps) { // independent LSPs (per-frame)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1656 if (get_bits_left(gb) < s->frame_lsp_bitsize) return 1;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1657 skip_bits_long(gb, s->frame_lsp_bitsize);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1658 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1659 bd_idx = s->vbm_tree[get_vlc2(gb, frame_type_vlc.table, 6, 3)];
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1660 if (bd_idx < 0)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1661 return -1; // invalid frame type VLC code
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1662 frame_desc = &frame_descs[bd_idx];
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1663 if (frame_desc->acb_type == ACB_TYPE_ASYMMETRIC) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1664 if (get_bits_left(gb) < s->pitch_nbits)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1665 return 1;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1666 skip_bits_long(gb, s->pitch_nbits);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1667 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1668 if (frame_desc->fcb_type == FCB_TYPE_SILENCE) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1669 skip_bits(gb, 8);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1670 } else if (frame_desc->fcb_type == FCB_TYPE_AW_PULSES) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1671 int tmp = get_bits(gb, 6);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1672 if (tmp >= 0x36) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1673 skip_bits(gb, 2);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1674 aw_idx_is_ext = 1;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1675 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1676 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1677
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1678 /* blocks */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1679 if (frame_desc->acb_type == ACB_TYPE_HAMMING) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1680 need_bits = s->block_pitch_nbits +
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1681 (frame_desc->n_blocks - 1) * s->block_delta_pitch_nbits;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1682 } else if (frame_desc->fcb_type == FCB_TYPE_AW_PULSES) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1683 need_bits = 2 * !aw_idx_is_ext;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1684 } else
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1685 need_bits = 0;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1686 need_bits += frame_desc->frame_size;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1687 if (get_bits_left(gb) < need_bits)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1688 return 1;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1689 skip_bits_long(gb, need_bits);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1690 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1691
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1692 return 0;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1693 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1694
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1695 /**
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1696 * Synthesize output samples for a single superframe. If we have any data
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1697 * cached in s->sframe_cache, that will be used instead of whatever is loaded
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1698 * in s->gb.
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1699 *
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1700 * WMA Voice superframes contain 3 frames, each containing 160 audio samples,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1701 * to give a total of 480 samples per frame. See #synth_frame() for frame
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1702 * parsing. In addition to 3 frames, superframes can also contain the LSPs
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1703 * (if these are globally specified for all frames (residually); they can
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1704 * also be specified individually per-frame. See the s->has_residual_lsps
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1705 * option), and can specify the number of samples encoded in this superframe
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1706 * (if less than 480), usually used to prevent blanks at track boundaries.
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1707 *
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1708 * @param ctx WMA Voice decoder context
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1709 * @param samples pointer to output buffer for voice samples
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1710 * @param data_size pointer containing the size of #samples on input, and the
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1711 * amount of #samples filled on output
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1712 * @return 0 on success, <0 on error or 1 if there was not enough data to
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1713 * fully parse the superframe
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1714 */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1715 static int synth_superframe(AVCodecContext *ctx,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1716 float *samples, int *data_size)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1717 {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1718 WMAVoiceContext *s = ctx->priv_data;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1719 GetBitContext *gb = &s->gb, s_gb;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1720 int n, res, n_samples = 480;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1721 double lsps[MAX_FRAMES][MAX_LSPS];
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1722 const double *mean_lsf = s->lsps == 16 ?
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1723 wmavoice_mean_lsf16[s->lsp_def_mode] : wmavoice_mean_lsf10[s->lsp_def_mode];
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1724 float excitation[MAX_SIGNAL_HISTORY + MAX_SFRAMESIZE + 12];
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1725 float synth[MAX_LSPS + MAX_SFRAMESIZE];
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1726
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1727 memcpy(synth, s->synth_history,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1728 s->lsps * sizeof(*synth));
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1729 memcpy(excitation, s->excitation_history,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1730 s->history_nsamples * sizeof(*excitation));
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1731
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1732 if (s->sframe_cache_size > 0) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1733 gb = &s_gb;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1734 init_get_bits(gb, s->sframe_cache, s->sframe_cache_size);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1735 s->sframe_cache_size = 0;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1736 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1737
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1738 if ((res = check_bits_for_superframe(gb, s)) == 1) return 1;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1739
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1740 /* First bit is speech/music bit, it differentiates between WMAVoice
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1741 * speech samples (the actual codec) and WMAVoice music samples, which
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1742 * are really WMAPro-in-WMAVoice-superframes. I've never seen those in
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1743 * the wild yet. */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1744 if (!get_bits1(gb)) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1745 av_log_missing_feature(ctx, "WMAPro-in-WMAVoice support", 1);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1746 return -1;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1747 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1748
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1749 /* (optional) nr. of samples in superframe; always <= 480 and >= 0 */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1750 if (get_bits1(gb)) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1751 if ((n_samples = get_bits(gb, 12)) > 480) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1752 av_log(ctx, AV_LOG_ERROR,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1753 "Superframe encodes >480 samples (%d), not allowed\n",
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1754 n_samples);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1755 return -1;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1756 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1757 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1758 /* Parse LSPs, if global for the superframe (can also be per-frame). */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1759 if (s->has_residual_lsps) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1760 double prev_lsps[MAX_LSPS], a1[MAX_LSPS * 2], a2[MAX_LSPS * 2];
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1761
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1762 for (n = 0; n < s->lsps; n++)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1763 prev_lsps[n] = s->prev_lsps[n] - mean_lsf[n];
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1764
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1765 if (s->lsps == 10) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1766 dequant_lsp10r(gb, lsps[2], prev_lsps, a1, a2, s->lsp_q_mode);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1767 } else /* s->lsps == 16 */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1768 dequant_lsp16r(gb, lsps[2], prev_lsps, a1, a2, s->lsp_q_mode);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1769
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1770 for (n = 0; n < s->lsps; n++) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1771 lsps[0][n] = mean_lsf[n] + (a1[n] - a2[n * 2]);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1772 lsps[1][n] = mean_lsf[n] + (a1[s->lsps + n] - a2[n * 2 + 1]);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1773 lsps[2][n] += mean_lsf[n];
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1774 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1775 for (n = 0; n < 3; n++)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1776 stabilize_lsps(lsps[n], s->lsps);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1777 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1778
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1779 /* Parse frames, optionally preceeded by per-frame (independent) LSPs. */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1780 for (n = 0; n < 3; n++) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1781 if (!s->has_residual_lsps) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1782 int m;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1783
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1784 if (s->lsps == 10) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1785 dequant_lsp10i(gb, lsps[n]);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1786 } else /* s->lsps == 16 */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1787 dequant_lsp16i(gb, lsps[n]);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1788
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1789 for (m = 0; m < s->lsps; m++)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1790 lsps[n][m] += mean_lsf[m];
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1791 stabilize_lsps(lsps[n], s->lsps);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1792 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1793
11653
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
1794 if ((res = synth_frame(ctx, gb, n,
11123
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1795 &samples[n * MAX_FRAMESIZE],
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1796 lsps[n], n == 0 ? s->prev_lsps : lsps[n - 1],
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1797 &excitation[s->history_nsamples + n * MAX_FRAMESIZE],
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1798 &synth[s->lsps + n * MAX_FRAMESIZE])))
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1799 return res;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1800 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1801
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1802 /* Statistics? FIXME - we don't check for length, a slight overrun
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1803 * will be caught by internal buffer padding, and anything else
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1804 * will be skipped, not read. */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1805 if (get_bits1(gb)) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1806 res = get_bits(gb, 4);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1807 skip_bits(gb, 10 * (res + 1));
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1808 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1809
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1810 /* Specify nr. of output samples */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1811 *data_size = n_samples * sizeof(float);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1812
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1813 /* Update history */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1814 memcpy(s->prev_lsps, lsps[2],
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1815 s->lsps * sizeof(*s->prev_lsps));
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1816 memcpy(s->synth_history, &synth[MAX_SFRAMESIZE],
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1817 s->lsps * sizeof(*synth));
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1818 memcpy(s->excitation_history, &excitation[MAX_SFRAMESIZE],
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1819 s->history_nsamples * sizeof(*excitation));
11653
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
1820 if (s->do_apf)
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
1821 memmove(s->zero_exc_pf, &s->zero_exc_pf[MAX_SFRAMESIZE],
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
1822 s->history_nsamples * sizeof(*s->zero_exc_pf));
11123
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1823
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1824 return 0;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1825 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1826
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1827 /**
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1828 * Parse the packet header at the start of each packet (input data to this
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1829 * decoder).
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1830 *
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1831 * @param s WMA Voice decoding context private data
11556
b94e1810ce4c Replace @returns by @return.
benoit
parents: 11123
diff changeset
1832 * @return 1 if not enough bits were available, or 0 on success.
11123
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1833 */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1834 static int parse_packet_header(WMAVoiceContext *s)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1835 {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1836 GetBitContext *gb = &s->gb;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1837 unsigned int res;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1838
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1839 if (get_bits_left(gb) < 11)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1840 return 1;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1841 skip_bits(gb, 4); // packet sequence number
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1842 s->has_residual_lsps = get_bits1(gb);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1843 do {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1844 res = get_bits(gb, 6); // number of superframes per packet
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1845 // (minus first one if there is spillover)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1846 if (get_bits_left(gb) < 6 * (res == 0x3F) + s->spillover_bitsize)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1847 return 1;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1848 } while (res == 0x3F);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1849 s->spillover_nbits = get_bits(gb, s->spillover_bitsize);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1850
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1851 return 0;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1852 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1853
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1854 /**
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1855 * Copy (unaligned) bits from gb/data/size to pb.
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1856 *
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1857 * @param pb target buffer to copy bits into
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1858 * @param data source buffer to copy bits from
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1859 * @param size size of the source data, in bytes
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1860 * @param gb bit I/O context specifying the current position in the source.
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1861 * data. This function might use this to align the bit position to
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1862 * a whole-byte boundary before calling #ff_copy_bits() on aligned
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1863 * source data
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1864 * @param nbits the amount of bits to copy from source to target
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1865 *
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1866 * @note after calling this function, the current position in the input bit
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1867 * I/O context is undefined.
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1868 */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1869 static void copy_bits(PutBitContext *pb,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1870 const uint8_t *data, int size,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1871 GetBitContext *gb, int nbits)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1872 {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1873 int rmn_bytes, rmn_bits;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1874
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1875 rmn_bits = rmn_bytes = get_bits_left(gb);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1876 if (rmn_bits < nbits)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1877 return;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1878 rmn_bits &= 7; rmn_bytes >>= 3;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1879 if ((rmn_bits = FFMIN(rmn_bits, nbits)) > 0)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1880 put_bits(pb, rmn_bits, get_bits(gb, rmn_bits));
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1881 ff_copy_bits(pb, data + size - rmn_bytes,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1882 FFMIN(nbits - rmn_bits, rmn_bytes << 3));
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1883 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1884
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1885 /**
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1886 * Packet decoding: a packet is anything that the (ASF) demuxer contains,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1887 * and we expect that the demuxer / application provides it to us as such
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1888 * (else you'll probably get garbage as output). Every packet has a size of
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1889 * ctx->block_align bytes, starts with a packet header (see
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1890 * #parse_packet_header()), and then a series of superframes. Superframe
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1891 * boundaries may exceed packets, i.e. superframes can split data over
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1892 * multiple (two) packets.
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1893 *
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1894 * For more information about frames, see #synth_superframe().
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1895 */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1896 static int wmavoice_decode_packet(AVCodecContext *ctx, void *data,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1897 int *data_size, AVPacket *avpkt)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1898 {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1899 WMAVoiceContext *s = ctx->priv_data;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1900 GetBitContext *gb = &s->gb;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1901 int size, res, pos;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1902
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1903 if (*data_size < 480 * sizeof(float)) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1904 av_log(ctx, AV_LOG_ERROR,
12150
30867f2c9009 Use correct length modifier for size comparison in printf expression, fixes:
diego
parents: 12114
diff changeset
1905 "Output buffer too small (%d given - %zu needed)\n",
11123
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1906 *data_size, 480 * sizeof(float));
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1907 return -1;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1908 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1909 *data_size = 0;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1910
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1911 /* Packets are sometimes a multiple of ctx->block_align, with a packet
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1912 * header at each ctx->block_align bytes. However, FFmpeg's ASF demuxer
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1913 * feeds us ASF packets, which may concatenate multiple "codec" packets
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1914 * in a single "muxer" packet, so we artificially emulate that by
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1915 * capping the packet size at ctx->block_align. */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1916 for (size = avpkt->size; size > ctx->block_align; size -= ctx->block_align);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1917 if (!size)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1918 return 0;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1919 init_get_bits(&s->gb, avpkt->data, size << 3);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1920
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1921 /* size == ctx->block_align is used to indicate whether we are dealing with
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1922 * a new packet or a packet of which we already read the packet header
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1923 * previously. */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1924 if (size == ctx->block_align) { // new packet header
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1925 if ((res = parse_packet_header(s)) < 0)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1926 return res;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1927
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1928 /* If the packet header specifies a s->spillover_nbits, then we want
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1929 * to push out all data of the previous packet (+ spillover) before
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1930 * continuing to parse new superframes in the current packet. */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1931 if (s->spillover_nbits > 0) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1932 if (s->sframe_cache_size > 0) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1933 int cnt = get_bits_count(gb);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1934 copy_bits(&s->pb, avpkt->data, size, gb, s->spillover_nbits);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1935 flush_put_bits(&s->pb);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1936 s->sframe_cache_size += s->spillover_nbits;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1937 if ((res = synth_superframe(ctx, data, data_size)) == 0 &&
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1938 *data_size > 0) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1939 cnt += s->spillover_nbits;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1940 s->skip_bits_next = cnt & 7;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1941 return cnt >> 3;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1942 } else
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1943 skip_bits_long (gb, s->spillover_nbits - cnt +
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1944 get_bits_count(gb)); // resync
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1945 } else
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1946 skip_bits_long(gb, s->spillover_nbits); // resync
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1947 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1948 } else if (s->skip_bits_next)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1949 skip_bits(gb, s->skip_bits_next);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1950
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1951 /* Try parsing superframes in current packet */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1952 s->sframe_cache_size = 0;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1953 s->skip_bits_next = 0;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1954 pos = get_bits_left(gb);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1955 if ((res = synth_superframe(ctx, data, data_size)) < 0) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1956 return res;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1957 } else if (*data_size > 0) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1958 int cnt = get_bits_count(gb);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1959 s->skip_bits_next = cnt & 7;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1960 return cnt >> 3;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1961 } else if ((s->sframe_cache_size = pos) > 0) {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1962 /* rewind bit reader to start of last (incomplete) superframe... */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1963 init_get_bits(gb, avpkt->data, size << 3);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1964 skip_bits_long(gb, (size << 3) - pos);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1965 assert(get_bits_left(gb) == pos);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1966
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1967 /* ...and cache it for spillover in next packet */
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1968 init_put_bits(&s->pb, s->sframe_cache, SFRAME_CACHE_MAXSIZE);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1969 copy_bits(&s->pb, avpkt->data, size, gb, s->sframe_cache_size);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1970 // FIXME bad - just copy bytes as whole and add use the
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1971 // skip_bits_next field
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1972 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1973
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1974 return size;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1975 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1976
11653
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
1977 static av_cold int wmavoice_decode_end(AVCodecContext *ctx)
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
1978 {
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
1979 WMAVoiceContext *s = ctx->priv_data;
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
1980
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
1981 if (s->do_apf) {
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
1982 ff_rdft_end(&s->rdft);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
1983 ff_rdft_end(&s->irdft);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
1984 ff_dct_end(&s->dct);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
1985 ff_dct_end(&s->dst);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
1986 }
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
1987
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
1988 return 0;
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
1989 }
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
1990
11123
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1991 static av_cold void wmavoice_flush(AVCodecContext *ctx)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1992 {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1993 WMAVoiceContext *s = ctx->priv_data;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1994 int n;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1995
11653
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
1996 s->postfilter_agc = 0;
11123
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1997 s->sframe_cache_size = 0;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1998 s->skip_bits_next = 0;
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
1999 for (n = 0; n < s->lsps; n++)
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
2000 s->prev_lsps[n] = M_PI * (n + 1.0) / (s->lsps + 1.0);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
2001 memset(s->excitation_history, 0,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
2002 sizeof(*s->excitation_history) * MAX_SIGNAL_HISTORY);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
2003 memset(s->synth_history, 0,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
2004 sizeof(*s->synth_history) * MAX_LSPS);
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
2005 memset(s->gain_pred_err, 0,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
2006 sizeof(s->gain_pred_err));
11653
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
2007
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
2008 if (s->do_apf) {
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
2009 memset(&s->synth_filter_out_buf[MAX_LSPS_ALIGN16 - s->lsps], 0,
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
2010 sizeof(*s->synth_filter_out_buf) * s->lsps);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
2011 memset(s->dcf_mem, 0,
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
2012 sizeof(*s->dcf_mem) * 2);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
2013 memset(s->zero_exc_pf, 0,
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
2014 sizeof(*s->zero_exc_pf) * s->history_nsamples);
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
2015 memset(s->denoise_filter_cache, 0, sizeof(s->denoise_filter_cache));
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
2016 }
11123
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
2017 }
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
2018
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
2019 AVCodec wmavoice_decoder = {
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
2020 "wmavoice",
11560
8a4984c5cacc Define AVMediaType enum, and use it instead of enum CodecType, which
stefano
parents: 11556
diff changeset
2021 AVMEDIA_TYPE_AUDIO,
11123
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
2022 CODEC_ID_WMAVOICE,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
2023 sizeof(WMAVoiceContext),
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
2024 wmavoice_decode_init,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
2025 NULL,
11653
28638e0d6e63 WMAVoice postfilter.
rbultje
parents: 11652
diff changeset
2026 wmavoice_decode_end,
11123
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
2027 wmavoice_decode_packet,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
2028 CODEC_CAP_SUBFRAMES,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
2029 .flush = wmavoice_flush,
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
2030 .long_name = NULL_IF_CONFIG_SMALL("Windows Media Audio Voice"),
d59349627f52 WMAVoice decoder.
rbultje
parents:
diff changeset
2031 };