FFmpeg coverage

Directory: ../../../ffmpeg/
File: src/libavcodec/mlpdec.c
Date: 2024-04-19 17:50:32
Exec Total Coverage
Lines: 480 678 70.8%
Functions: 18 19 94.7%
Branches: 290 459 63.2%
Line Branch Exec Source
1 /*
2 * MLP decoder
3 * Copyright (c) 2007-2008 Ian Caulfield
4 *
5 * This file is part of FFmpeg.
6 *
7 * FFmpeg is free software; you can redistribute it and/or
8 * modify it under the terms of the GNU Lesser General Public
9 * License as published by the Free Software Foundation; either
10 * version 2.1 of the License, or (at your option) any later version.
11 *
12 * FFmpeg is distributed in the hope that it will be useful,
13 * but WITHOUT ANY WARRANTY; without even the implied warranty of
14 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
15 * Lesser General Public License for more details.
16 *
17 * You should have received a copy of the GNU Lesser General Public
18 * License along with FFmpeg; if not, write to the Free Software
19 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
20 */
21
22 /**
23 * @file
24 * MLP decoder
25 */
26
27 #include "config_components.h"
28
29 #include <stdint.h>
30
31 #include "avcodec.h"
32 #include "libavutil/internal.h"
33 #include "libavutil/intreadwrite.h"
34 #include "libavutil/channel_layout.h"
35 #include "libavutil/mem_internal.h"
36 #include "libavutil/thread.h"
37 #include "libavutil/opt.h"
38 #include "codec_internal.h"
39 #include "decode.h"
40 #include "get_bits.h"
41 #include "mlp_parse.h"
42 #include "mlpdsp.h"
43 #include "mlp.h"
44 #include "config.h"
45 #include "profiles.h"
46
47 /** number of bits used for VLC lookup - longest Huffman code is 9 */
48 #if ARCH_ARM
49 #define VLC_BITS 5
50 #define VLC_STATIC_SIZE 64
51 #else
52 #define VLC_BITS 9
53 #define VLC_STATIC_SIZE 512
54 #endif
55
56 typedef struct SubStream {
57 /// Set if a valid restart header has been read. Otherwise the substream cannot be decoded.
58 uint8_t restart_seen;
59 /// Set if end of stream is encountered
60 uint8_t end_of_stream;
61
62 //@{
63 /** restart header data */
64 /// The type of noise to be used in the rematrix stage.
65 uint16_t noise_type;
66
67 /// The index of the first channel coded in this substream.
68 uint8_t min_channel;
69 /// The index of the last channel coded in this substream.
70 uint8_t max_channel;
71 /// The coded channels mask in this substream.
72 uint64_t coded_channels;
73 /// The number of channels input into the rematrix stage.
74 uint8_t max_matrix_channel;
75 /// For each channel output by the matrix, the output channel to map it to
76 uint8_t ch_assign[MAX_CHANNELS];
77 /// The channel layout for this substream
78 uint64_t mask;
79 /// The matrix encoding mode for this substream
80 enum AVMatrixEncoding matrix_encoding;
81 enum AVMatrixEncoding prev_matrix_encoding;
82
83 /// Channel coding parameters for channels in the substream
84 ChannelParams channel_params[MAX_CHANNELS];
85
86 /// The left shift applied to random noise in 0x31ea substreams.
87 uint8_t noise_shift;
88 /// The current seed value for the pseudorandom noise generator(s).
89 uint32_t noisegen_seed;
90
91 /// Set if the substream contains extra info to check the size of VLC blocks.
92 uint8_t data_check_present;
93
94 /// Bitmask of which parameter sets are conveyed in a decoding parameter block.
95 uint8_t param_presence_flags;
96 //@}
97
98 //@{
99 /** matrix data */
100
101 /// Number of matrices to be applied.
102 uint8_t num_primitive_matrices;
103
104 /// matrix output channel
105 uint8_t matrix_out_ch[MAX_MATRICES];
106
107 /// Whether the LSBs of the matrix output are encoded in the bitstream.
108 uint8_t lsb_bypass[MAX_MATRICES];
109 /// Matrix coefficients, stored as 2.14 fixed point.
110 DECLARE_ALIGNED(32, int32_t, matrix_coeff)[MAX_MATRICES][MAX_CHANNELS];
111 /// Left shift to apply to noise values in 0x31eb substreams.
112 uint8_t matrix_noise_shift[MAX_MATRICES];
113 //@}
114
115 /// Left shift to apply to Huffman-decoded residuals.
116 uint8_t quant_step_size[MAX_CHANNELS];
117
118 /// number of PCM samples in current audio block
119 uint16_t blocksize;
120 /// Number of PCM samples decoded so far in this frame.
121 uint16_t blockpos;
122
123 /// Left shift to apply to decoded PCM values to get final 24-bit output.
124 int8_t output_shift[MAX_CHANNELS];
125
126 /// Running XOR of all output samples.
127 int32_t lossless_check_data;
128
129 } SubStream;
130
131 typedef struct MLPDecodeContext {
132 const AVClass *class;
133 AVCodecContext *avctx;
134
135 AVChannelLayout downmix_layout;
136
137 /// Current access unit being read has a major sync.
138 int is_major_sync_unit;
139
140 /// Size of the major sync unit, in bytes
141 int major_sync_header_size;
142
143 /// Set if a valid major sync block has been read. Otherwise no decoding is possible.
144 uint8_t params_valid;
145
146 /// Number of substreams contained within this stream.
147 uint8_t num_substreams;
148
149 /// Which substream of substreams carry 16-channel presentation
150 uint8_t extended_substream_info;
151
152 /// Which substream of substreams carry 2/6/8-channel presentation
153 uint8_t substream_info;
154
155 /// Index of the last substream to decode - further substreams are skipped.
156 uint8_t max_decoded_substream;
157
158 /// Stream needs channel reordering to comply with FFmpeg's channel order
159 uint8_t needs_reordering;
160
161 /// number of PCM samples contained in each frame
162 int access_unit_size;
163 /// next power of two above the number of samples in each frame
164 int access_unit_size_pow2;
165
166 SubStream substream[MAX_SUBSTREAMS];
167
168 int matrix_changed;
169 int filter_changed[MAX_CHANNELS][NUM_FILTERS];
170
171 int8_t noise_buffer[MAX_BLOCKSIZE_POW2];
172 int8_t bypassed_lsbs[MAX_BLOCKSIZE][MAX_CHANNELS];
173 DECLARE_ALIGNED(32, int32_t, sample_buffer)[MAX_BLOCKSIZE][MAX_CHANNELS];
174
175 MLPDSPContext dsp;
176 int32_t (*pack_output)(int32_t lossless_check_data,
177 uint16_t blockpos,
178 int32_t (*sample_buffer)[MAX_CHANNELS],
179 void *data,
180 uint8_t *ch_assign,
181 int8_t *output_shift,
182 uint8_t max_matrix_channel,
183 int is32);
184 } MLPDecodeContext;
185
186 static const enum AVChannel thd_channel_order[] = {
187 AV_CHAN_FRONT_LEFT, AV_CHAN_FRONT_RIGHT, // LR
188 AV_CHAN_FRONT_CENTER, // C
189 AV_CHAN_LOW_FREQUENCY, // LFE
190 AV_CHAN_SIDE_LEFT, AV_CHAN_SIDE_RIGHT, // LRs
191 AV_CHAN_TOP_FRONT_LEFT, AV_CHAN_TOP_FRONT_RIGHT, // LRvh
192 AV_CHAN_FRONT_LEFT_OF_CENTER, AV_CHAN_FRONT_RIGHT_OF_CENTER, // LRc
193 AV_CHAN_BACK_LEFT, AV_CHAN_BACK_RIGHT, // LRrs
194 AV_CHAN_BACK_CENTER, // Cs
195 AV_CHAN_TOP_CENTER, // Ts
196 AV_CHAN_SURROUND_DIRECT_LEFT, AV_CHAN_SURROUND_DIRECT_RIGHT, // LRsd
197 AV_CHAN_WIDE_LEFT, AV_CHAN_WIDE_RIGHT, // LRw
198 AV_CHAN_TOP_FRONT_CENTER, // Cvh
199 AV_CHAN_LOW_FREQUENCY_2, // LFE2
200 };
201
202 2220 static int mlp_channel_layout_subset(AVChannelLayout *layout, uint64_t mask)
203 {
204
3/4
✓ Branch 1 taken 215 times.
✓ Branch 2 taken 2005 times.
✓ Branch 3 taken 215 times.
✗ Branch 4 not taken.
 2435 return av_channel_layout_check(layout) &&
205 215 av_channel_layout_subset(layout, mask) ==
206 215 av_channel_layout_subset(layout, UINT64_MAX);
207 }
208
209 3958 static enum AVChannel thd_channel_layout_extract_channel(uint64_t channel_layout,
210 int index)
211 {
212 int i;
213
214
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 3958 times.
 3958 if (av_popcount64(channel_layout) <= index)
215 ✗ return AV_CHAN_NONE;
216
217
1/2
✓ Branch 0 taken 11285 times.
✗ Branch 1 not taken.
 11285 for (i = 0; i < FF_ARRAY_ELEMS(thd_channel_order); i++)
218
4/4
✓ Branch 0 taken 11237 times.
✓ Branch 1 taken 48 times.
✓ Branch 2 taken 3958 times.
✓ Branch 3 taken 7279 times.
 11285 if (channel_layout & (1ULL << thd_channel_order[i]) && !index--)
219 3958 return thd_channel_order[i];
220 ✗ return AV_CHAN_NONE;
221 }
222
223 static VLC huff_vlc[3];
224
225 /** Initialize static data, constant between all invocations of the codec. */
226
227 9 static av_cold void init_static(void)
228 {
229
2/2
✓ Branch 0 taken 27 times.
✓ Branch 1 taken 9 times.
 36 for (int i = 0; i < 3; i++) {
230 static VLCElem vlc_buf[3 * VLC_STATIC_SIZE];
231 27 huff_vlc[i].table = &vlc_buf[i * VLC_STATIC_SIZE];
232 27 huff_vlc[i].table_allocated = VLC_STATIC_SIZE;
233 27 vlc_init(&huff_vlc[i], VLC_BITS, 18,
234 &ff_mlp_huffman_tables[i][0][1], 2, 1,
235 &ff_mlp_huffman_tables[i][0][0], 2, 1, VLC_INIT_USE_STATIC);
236 }
237
238 9 ff_mlp_init_crc();
239 9 }
240
241 16834 static inline int32_t calculate_sign_huff(MLPDecodeContext *m,
242 unsigned int substr, unsigned int ch)
243 {
244 16834 SubStream *s = &m->substream[substr];
245 16834 ChannelParams *cp = &s->channel_params[ch];
246 16834 int lsb_bits = cp->huff_lsbs - s->quant_step_size[ch];
247
2/2
✓ Branch 0 taken 11899 times.
✓ Branch 1 taken 4935 times.
 16834 int sign_shift = lsb_bits + (cp->codebook ? 2 - cp->codebook : -1);
248 16834 int32_t sign_huff_offset = cp->huff_offset;
249
250
2/2
✓ Branch 0 taken 11899 times.
✓ Branch 1 taken 4935 times.
 16834 if (cp->codebook > 0)
251 11899 sign_huff_offset -= 7 << lsb_bits;
252
253
2/2
✓ Branch 0 taken 15281 times.
✓ Branch 1 taken 1553 times.
 16834 if (sign_shift >= 0)
254 15281 sign_huff_offset -= 1 << sign_shift;
255
256 16834 return sign_huff_offset;
257 }
258
259 /** Read a sample, consisting of either, both or neither of entropy-coded MSBs
260 * and plain LSBs. */
261
262 1121160 static inline int read_huff_channels(MLPDecodeContext *m, GetBitContext *gbp,
263 unsigned int substr, unsigned int pos)
264 {
265 1121160 SubStream *s = &m->substream[substr];
266 unsigned int mat, channel;
267
268
2/2
✓ Branch 0 taken 2054520 times.
✓ Branch 1 taken 1121160 times.
 3175680 for (mat = 0; mat < s->num_primitive_matrices; mat++)
269
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 2054520 times.
 2054520 if (s->lsb_bypass[mat])
270 ✗ m->bypassed_lsbs[pos + s->blockpos][mat] = get_bits1(gbp);
271
272
2/2
✓ Branch 0 taken 2734080 times.
✓ Branch 1 taken 1121160 times.
 3855240 for (channel = s->min_channel; channel <= s->max_channel; channel++) {
273 2734080 ChannelParams *cp = &s->channel_params[channel];
274 2734080 int codebook = cp->codebook;
275 2734080 int quant_step_size = s->quant_step_size[channel];
276 2734080 int lsb_bits = cp->huff_lsbs - quant_step_size;
277 2734080 int result = 0;
278
279
2/2
✓ Branch 0 taken 2490712 times.
✓ Branch 1 taken 243368 times.
 2734080 if (codebook > 0)
280 2490712 result = get_vlc2(gbp, huff_vlc[codebook-1].table,
281 VLC_BITS, (9 + VLC_BITS - 1) / VLC_BITS);
282
283
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 2734080 times.
 2734080 if (result < 0)
284 ✗ return AVERROR_INVALIDDATA;
285
286
2/2
✓ Branch 0 taken 2549944 times.
✓ Branch 1 taken 184136 times.
 2734080 if (lsb_bits > 0)
287 2549944 result = (result << lsb_bits) + get_bits_long(gbp, lsb_bits);
288
289 2734080 result += cp->sign_huff_offset;
290 2734080 result *= 1 << quant_step_size;
291
292 2734080 m->sample_buffer[pos + s->blockpos][channel] = result;
293 }
294
295 1121160 return 0;
296 }
297
298 15 static av_cold int mlp_decode_init(AVCodecContext *avctx)
299 {
300 static AVOnce init_static_once = AV_ONCE_INIT;
301 15 MLPDecodeContext *m = avctx->priv_data;
302 int substr;
303
304 15 m->avctx = avctx;
305
2/2
✓ Branch 0 taken 60 times.
✓ Branch 1 taken 15 times.
 75 for (substr = 0; substr < MAX_SUBSTREAMS; substr++)
306 60 m->substream[substr].lossless_check_data = 0xffffffff;
307 15 ff_mlpdsp_init(&m->dsp);
308
309
2/2
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 13 times.
 15 if (m->downmix_layout.nb_channels) {
310
1/4
✗ Branch 1 not taken.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
 2 if (!av_channel_layout_compare(&m->downmix_layout, &(AVChannelLayout)AV_CHANNEL_LAYOUT_STEREO) ||
311 ✗ !av_channel_layout_compare(&m->downmix_layout, &(AVChannelLayout)AV_CHANNEL_LAYOUT_STEREO_DOWNMIX)) {
312 2 av_channel_layout_uninit(&avctx->ch_layout);
313 2 avctx->ch_layout = (AVChannelLayout)AV_CHANNEL_LAYOUT_STEREO;
314 ✗ } else if (!av_channel_layout_compare(&m->downmix_layout, &(AVChannelLayout)AV_CHANNEL_LAYOUT_5POINT0)) {
315 ✗ av_channel_layout_uninit(&avctx->ch_layout);
316 ✗ avctx->ch_layout = (AVChannelLayout)AV_CHANNEL_LAYOUT_5POINT0;
317 ✗ } else if (!av_channel_layout_compare(&m->downmix_layout, &(AVChannelLayout)AV_CHANNEL_LAYOUT_5POINT1)) {
318 ✗ av_channel_layout_uninit(&avctx->ch_layout);
319 ✗ avctx->ch_layout = (AVChannelLayout)AV_CHANNEL_LAYOUT_5POINT1;
320 }
321 else
322 ✗ av_log(avctx, AV_LOG_WARNING, "Invalid downmix layout\n");
323 }
324
325 15 ff_thread_once(&init_static_once, init_static);
326
327 15 return 0;
328 }
329
330 /** Read a major sync info header - contains high level information about
331 * the stream - sample rate, channel arrangement etc. Most of this
332 * information is not actually necessary for decoding, only for playback.
333 */
334
335 1774 static int read_major_sync(MLPDecodeContext *m, GetBitContext *gb)
336 {
337 MLPHeaderInfo mh;
338 int substr, ret;
339
340
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 1774 times.
 1774 if ((ret = ff_mlp_read_major_sync(m->avctx, &mh, gb)) != 0)
341 ✗ return ret;
342
343
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 1774 times.
 1774 if (mh.group1_bits == 0) {
344 ✗ av_log(m->avctx, AV_LOG_ERROR, "invalid/unknown bits per sample\n");
345 ✗ return AVERROR_INVALIDDATA;
346 }
347
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 1774 times.
 1774 if (mh.group2_bits > mh.group1_bits) {
348 ✗ av_log(m->avctx, AV_LOG_ERROR,
349 "Channel group 2 cannot have more bits per sample than group 1.\n");
350 ✗ return AVERROR_INVALIDDATA;
351 }
352
353
1/4
✗ Branch 0 not taken.
✓ Branch 1 taken 1774 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
 1774 if (mh.group2_samplerate && mh.group2_samplerate != mh.group1_samplerate) {
354 ✗ av_log(m->avctx, AV_LOG_ERROR,
355 "Channel groups with differing sample rates are not currently supported.\n");
356 ✗ return AVERROR_INVALIDDATA;
357 }
358
359
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 1774 times.
 1774 if (mh.group1_samplerate == 0) {
360 ✗ av_log(m->avctx, AV_LOG_ERROR, "invalid/unknown sampling rate\n");
361 ✗ return AVERROR_INVALIDDATA;
362 }
363
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 1774 times.
 1774 if (mh.group1_samplerate > MAX_SAMPLERATE) {
364 ✗ av_log(m->avctx, AV_LOG_ERROR,
365 "Sampling rate %d is greater than the supported maximum (%d).\n",
366 mh.group1_samplerate, MAX_SAMPLERATE);
367 ✗ return AVERROR_INVALIDDATA;
368 }
369
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 1774 times.
 1774 if (mh.access_unit_size > MAX_BLOCKSIZE) {
370 ✗ av_log(m->avctx, AV_LOG_ERROR,
371 "Block size %d is greater than the supported maximum (%d).\n",
372 mh.access_unit_size, MAX_BLOCKSIZE);
373 ✗ return AVERROR_INVALIDDATA;
374 }
375
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 1774 times.
 1774 if (mh.access_unit_size_pow2 > MAX_BLOCKSIZE_POW2) {
376 ✗ av_log(m->avctx, AV_LOG_ERROR,
377 "Block size pow2 %d is greater than the supported maximum (%d).\n",
378 mh.access_unit_size_pow2, MAX_BLOCKSIZE_POW2);
379 ✗ return AVERROR_INVALIDDATA;
380 }
381
382
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 1774 times.
 1774 if (mh.num_substreams == 0)
383 ✗ return AVERROR_INVALIDDATA;
384
3/4
✓ Branch 0 taken 1117 times.
✓ Branch 1 taken 657 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 1117 times.
 1774 if (m->avctx->codec_id == AV_CODEC_ID_MLP && mh.num_substreams > 2) {
385 ✗ av_log(m->avctx, AV_LOG_ERROR, "MLP only supports up to 2 substreams.\n");
386 ✗ return AVERROR_INVALIDDATA;
387 }
388
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 1774 times.
 1774 if (mh.num_substreams > MAX_SUBSTREAMS) {
389 ✗ avpriv_request_sample(m->avctx,
390 "%d substreams (more than the "
391 "maximum supported by the decoder)",
392 mh.num_substreams);
393 ✗ return AVERROR_PATCHWELCOME;
394 }
395
396 1774 m->major_sync_header_size = mh.header_size;
397
398 1774 m->access_unit_size = mh.access_unit_size;
399 1774 m->access_unit_size_pow2 = mh.access_unit_size_pow2;
400
401 1774 m->num_substreams = mh.num_substreams;
402 1774 m->extended_substream_info = mh.extended_substream_info;
403 1774 m->substream_info = mh.substream_info;
404
405 /* If there is a 4th substream and the MSB of substream_info is set,
406 * there is a 16-channel spatial presentation (Atmos in TrueHD).
407 */
408
2/2
✓ Branch 0 taken 657 times.
✓ Branch 1 taken 1117 times.
 1774 if (m->avctx->codec_id == AV_CODEC_ID_TRUEHD
409
3/4
✓ Branch 0 taken 4 times.
✓ Branch 1 taken 653 times.
✓ Branch 2 taken 4 times.
✗ Branch 3 not taken.
 657 && m->num_substreams == 4 && m->substream_info >> 7 == 1) {
410 4 m->avctx->profile = AV_PROFILE_TRUEHD_ATMOS;
411 }
412
413 /* limit to decoding 3 substreams, as the 4th is used by Dolby Atmos for non-audio data */
414 1774 m->max_decoded_substream = FFMIN(m->num_substreams - 1, 2);
415
416 1774 m->avctx->sample_rate = mh.group1_samplerate;
417 1774 m->avctx->frame_size = mh.access_unit_size;
418
419 1774 m->avctx->bits_per_raw_sample = mh.group1_bits;
420
2/2
✓ Branch 0 taken 657 times.
✓ Branch 1 taken 1117 times.
 1774 if (mh.group1_bits > 16)
421 657 m->avctx->sample_fmt = AV_SAMPLE_FMT_S32;
422 else
423 1117 m->avctx->sample_fmt = AV_SAMPLE_FMT_S16;
424 3548 m->pack_output = m->dsp.mlp_select_pack_output(m->substream[m->max_decoded_substream].ch_assign,
425 1774 m->substream[m->max_decoded_substream].output_shift,
426 1774 m->substream[m->max_decoded_substream].max_matrix_channel,
427 1774 m->avctx->sample_fmt == AV_SAMPLE_FMT_S32);
428
429 1774 m->params_valid = 1;
430
2/2
✓ Branch 0 taken 7096 times.
✓ Branch 1 taken 1774 times.
 8870 for (substr = 0; substr < MAX_SUBSTREAMS; substr++)
431 7096 m->substream[substr].restart_seen = 0;
432
433 /* Set the layout for each substream. When there's more than one, the first
434 * substream is Stereo. Subsequent substreams' layouts are indicated in the
435 * major sync. */
436
2/2
✓ Branch 0 taken 1117 times.
✓ Branch 1 taken 657 times.
 1774 if (m->avctx->codec_id == AV_CODEC_ID_MLP) {
437
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 1117 times.
 1117 if (mh.stream_type != SYNC_MLP) {
438 ✗ avpriv_request_sample(m->avctx,
439 "unexpected stream_type %X in MLP",
440 mh.stream_type);
441 ✗ return AVERROR_PATCHWELCOME;
442 }
443
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 1117 times.
 1117 if ((substr = (mh.num_substreams > 1)))
444 ✗ m->substream[0].mask = AV_CH_LAYOUT_STEREO;
445 1117 m->substream[substr].mask = mh.channel_layout_mlp;
446 } else {
447
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 657 times.
 657 if (mh.stream_type != SYNC_TRUEHD) {
448 ✗ avpriv_request_sample(m->avctx,
449 "unexpected stream_type %X in !MLP",
450 mh.stream_type);
451 ✗ return AVERROR_PATCHWELCOME;
452 }
453 657 m->substream[1].mask = mh.channel_layout_thd_stream1;
454
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 657 times.
 657 if (mh.channels_thd_stream1 == 2 &&
455 ✗ mh.channels_thd_stream2 == 2 &&
456 ✗ m->avctx->ch_layout.nb_channels == 2)
457 ✗ m->substream[0].mask = AV_CH_LAYOUT_STEREO;
458
1/2
✓ Branch 0 taken 657 times.
✗ Branch 1 not taken.
 657 if ((substr = (mh.num_substreams > 1)))
459 657 m->substream[0].mask = AV_CH_LAYOUT_STEREO;
460
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 657 times.
 657 if (mh.num_substreams == 1 &&
461 ✗ mh.channels_thd_stream1 == 1 &&
462 ✗ mh.channels_thd_stream2 == 1 &&
463 ✗ m->avctx->ch_layout.nb_channels == 1)
464 ✗ m->substream[0].mask = AV_CH_LAYOUT_MONO;
465
2/2
✓ Branch 0 taken 4 times.
✓ Branch 1 taken 653 times.
 657 if (mh.num_substreams > 2)
466
1/2
✓ Branch 0 taken 4 times.
✗ Branch 1 not taken.
 4 if (mh.channel_layout_thd_stream2)
467 4 m->substream[2].mask = mh.channel_layout_thd_stream2;
468 else
469 ✗ m->substream[2].mask = mh.channel_layout_thd_stream1;
470
2/2
✓ Branch 0 taken 435 times.
✓ Branch 1 taken 222 times.
 657 if (m->avctx->ch_layout.nb_channels > 2)
471
2/2
✓ Branch 0 taken 4 times.
✓ Branch 1 taken 431 times.
 435 if (mh.num_substreams > 2)
472 4 m->substream[1].mask = mh.channel_layout_thd_stream1;
473 else
474 431 m->substream[mh.num_substreams > 1].mask = mh.channel_layout_thd_stream2;
475 }
476
477
1/4
✗ Branch 0 not taken.
✓ Branch 1 taken 1774 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
 1774 m->needs_reordering = mh.channel_arrangement >= 18 && mh.channel_arrangement <= 20;
478
479 /* Parse the TrueHD decoder channel modifiers and set each substream's
480 * AVMatrixEncoding accordingly.
481 *
482 * The meaning of the modifiers depends on the channel layout:
483 *
484 * - THD_CH_MODIFIER_LTRT, THD_CH_MODIFIER_LBINRBIN only apply to 2-channel
485 *
486 * - THD_CH_MODIFIER_MONO applies to 1-channel or 2-channel (dual mono)
487 *
488 * - THD_CH_MODIFIER_SURROUNDEX, THD_CH_MODIFIER_NOTSURROUNDEX only apply to
489 * layouts with an Ls/Rs channel pair
490 */
491
2/2
✓ Branch 0 taken 7096 times.
✓ Branch 1 taken 1774 times.
 8870 for (substr = 0; substr < MAX_SUBSTREAMS; substr++)
492 7096 m->substream[substr].matrix_encoding = AV_MATRIX_ENCODING_NONE;
493
2/2
✓ Branch 0 taken 657 times.
✓ Branch 1 taken 1117 times.
 1774 if (m->avctx->codec_id == AV_CODEC_ID_TRUEHD) {
494
2/2
✓ Branch 0 taken 4 times.
✓ Branch 1 taken 653 times.
 657 if (mh.num_substreams > 2 &&
495
1/2
✓ Branch 0 taken 4 times.
✗ Branch 1 not taken.
 4 mh.channel_layout_thd_stream2 & AV_CH_SIDE_LEFT &&
496
1/2
✓ Branch 0 taken 4 times.
✗ Branch 1 not taken.
 4 mh.channel_layout_thd_stream2 & AV_CH_SIDE_RIGHT &&
497
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 4 times.
 4 mh.channel_modifier_thd_stream2 == THD_CH_MODIFIER_SURROUNDEX)
498 ✗ m->substream[2].matrix_encoding = AV_MATRIX_ENCODING_DOLBYEX;
499
500
1/2
✓ Branch 0 taken 657 times.
✗ Branch 1 not taken.
 657 if (mh.num_substreams > 1 &&
501
2/2
✓ Branch 0 taken 649 times.
✓ Branch 1 taken 8 times.
 657 mh.channel_layout_thd_stream1 & AV_CH_SIDE_LEFT &&
502
1/2
✓ Branch 0 taken 649 times.
✗ Branch 1 not taken.
 649 mh.channel_layout_thd_stream1 & AV_CH_SIDE_RIGHT &&
503
2/2
✓ Branch 0 taken 4 times.
✓ Branch 1 taken 645 times.
 649 mh.channel_modifier_thd_stream1 == THD_CH_MODIFIER_SURROUNDEX)
504 4 m->substream[1].matrix_encoding = AV_MATRIX_ENCODING_DOLBYEX;
505
506
1/2
✓ Branch 0 taken 657 times.
✗ Branch 1 not taken.
 657 if (mh.num_substreams > 0)
507
2/3
✓ Branch 0 taken 4 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 653 times.
 657 switch (mh.channel_modifier_thd_stream0) {
508 4 case THD_CH_MODIFIER_LTRT:
509 4 m->substream[0].matrix_encoding = AV_MATRIX_ENCODING_DOLBY;
510 4 break;
511 ✗ case THD_CH_MODIFIER_LBINRBIN:
512 ✗ m->substream[0].matrix_encoding = AV_MATRIX_ENCODING_DOLBYHEADPHONE;
513 ✗ break;
514 653 default:
515 653 break;
516 }
517 }
518
519 1774 return 0;
520 }
521
522 /** Read a restart header from a block in a substream. This contains parameters
523 * required to decode the audio that do not change very often. Generally
524 * (always) present only in blocks following a major sync. */
525
526 2220 static int read_restart_header(MLPDecodeContext *m, GetBitContext *gbp,
527 const uint8_t *buf, unsigned int substr)
528 {
529 2220 SubStream *s = &m->substream[substr];
530 unsigned int ch;
531 int sync_word, tmp;
532 uint8_t checksum;
533 uint8_t lossless_check;
534 2220 int start_count = get_bits_count(gbp);
535 int min_channel, max_channel, max_matrix_channel, noise_type;
536 4440 const int std_max_matrix_channel = m->avctx->codec_id == AV_CODEC_ID_MLP
537 ? MAX_MATRIX_CHANNEL_MLP
538
2/2
✓ Branch 0 taken 1117 times.
✓ Branch 1 taken 1103 times.
 2220 : MAX_MATRIX_CHANNEL_TRUEHD;
539
540 2220 sync_word = get_bits(gbp, 13);
541
542
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 2220 times.
 2220 if (sync_word != 0x31ea >> 1) {
543 ✗ av_log(m->avctx, AV_LOG_ERROR,
544 "restart header sync incorrect (got 0x%04x)\n", sync_word);
545 ✗ return AVERROR_INVALIDDATA;
546 }
547
548 2220 noise_type = get_bits1(gbp);
549
550
3/4
✓ Branch 0 taken 1117 times.
✓ Branch 1 taken 1103 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 1117 times.
 2220 if (m->avctx->codec_id == AV_CODEC_ID_MLP && noise_type) {
551 ✗ av_log(m->avctx, AV_LOG_ERROR, "MLP must have 0x31ea sync word.\n");
552 ✗ return AVERROR_INVALIDDATA;
553 }
554
555 2220 skip_bits(gbp, 16); /* Output timestamp */
556
557 2220 min_channel = get_bits(gbp, 4);
558 2220 max_channel = get_bits(gbp, 4);
559 2220 max_matrix_channel = get_bits(gbp, 4);
560
561
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 2220 times.
 2220 if (max_matrix_channel > std_max_matrix_channel) {
562 ✗ av_log(m->avctx, AV_LOG_ERROR,
563 "Max matrix channel cannot be greater than %d.\n",
564 std_max_matrix_channel);
565 ✗ return AVERROR_INVALIDDATA;
566 }
567
568 /* This should happen for TrueHD streams with >6 channels and MLP's noise
569 * type. It is not yet known if this is allowed. */
570
3/4
✓ Branch 0 taken 4 times.
✓ Branch 1 taken 2216 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 4 times.
 2220 if (max_matrix_channel > MAX_MATRIX_CHANNEL_MLP && !noise_type) {
571 ✗ avpriv_request_sample(m->avctx,
572 "%d channels (more than the "
573 "maximum supported by the decoder)",
574 max_channel + 2);
575 ✗ return AVERROR_PATCHWELCOME;
576 }
577
578
2/4
✓ Branch 0 taken 2220 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 2220 times.
 2220 if (max_channel + 1 > MAX_CHANNELS || max_channel + 1 < min_channel)
579 ✗ return AVERROR_INVALIDDATA;
580
581 2220 s->min_channel = min_channel;
582 2220 s->max_channel = max_channel;
583 2220 s->coded_channels = ((1LL << (max_channel - min_channel + 1)) - 1) << min_channel;
584 2220 s->max_matrix_channel = max_matrix_channel;
585 2220 s->noise_type = noise_type;
586
587
2/2
✓ Branch 1 taken 215 times.
✓ Branch 2 taken 2005 times.
 2220 if (mlp_channel_layout_subset(&m->downmix_layout, s->mask) &&
588
1/2
✓ Branch 0 taken 215 times.
✗ Branch 1 not taken.
 215 m->max_decoded_substream > substr) {
589 215 av_log(m->avctx, AV_LOG_DEBUG,
590 "Extracting %d-channel downmix (0x%"PRIx64") from substream %d. "
591 "Further substreams will be skipped.\n",
592 215 s->max_channel + 1, s->mask, substr);
593 215 m->max_decoded_substream = substr;
594 }
595
596 2220 s->noise_shift = get_bits(gbp, 4);
597 2220 s->noisegen_seed = get_bits(gbp, 23);
598
599 2220 skip_bits(gbp, 19);
600
601 2220 s->data_check_present = get_bits1(gbp);
602 2220 lossless_check = get_bits(gbp, 8);
603
2/2
✓ Branch 0 taken 1774 times.
✓ Branch 1 taken 446 times.
 2220 if (substr == m->max_decoded_substream
604
2/2
✓ Branch 0 taken 1759 times.
✓ Branch 1 taken 15 times.
 1774 && s->lossless_check_data != 0xffffffff) {
605 1759 tmp = xor_32_to_8(s->lossless_check_data);
606
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 1759 times.
 1759 if (tmp != lossless_check)
607 ✗ av_log(m->avctx, AV_LOG_WARNING,
608 "Lossless check failed - expected %02x, calculated %02x.\n",
609 lossless_check, tmp);
610 }
611
612 2220 skip_bits(gbp, 16);
613
614 2220 memset(s->ch_assign, 0, sizeof(s->ch_assign));
615
616
2/2
✓ Branch 0 taken 6192 times.
✓ Branch 1 taken 2220 times.
 8412 for (ch = 0; ch <= s->max_matrix_channel; ch++) {
617 6192 int ch_assign = get_bits(gbp, 6);
618
2/2
✓ Branch 0 taken 3958 times.
✓ Branch 1 taken 2234 times.
 6192 if (m->avctx->codec_id == AV_CODEC_ID_TRUEHD) {
619 AVChannelLayout l;
620 3958 enum AVChannel channel = thd_channel_layout_extract_channel(s->mask, ch_assign);
621
622 3958 av_channel_layout_from_mask(&l, s->mask);
623 3958 ch_assign = av_channel_layout_index_from_channel(&l, channel);
624 }
625
2/4
✓ Branch 0 taken 6192 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 6192 times.
 6192 if (ch_assign < 0 || ch_assign > s->max_matrix_channel) {
626 ✗ avpriv_request_sample(m->avctx,
627 "Assignment of matrix channel %d to invalid output channel %d",
628 ch, ch_assign);
629 ✗ return AVERROR_PATCHWELCOME;
630 }
631 6192 s->ch_assign[ch_assign] = ch;
632 }
633
634 2220 checksum = ff_mlp_restart_checksum(buf, get_bits_count(gbp) - start_count);
635
636
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 2220 times.
 2220 if (checksum != get_bits(gbp, 8))
637 ✗ av_log(m->avctx, AV_LOG_ERROR, "restart header checksum error\n");
638
639 /* Set default decoding parameters. */
640 2220 s->param_presence_flags = 0xff;
641 2220 s->num_primitive_matrices = 0;
642 2220 s->blocksize = 8;
643 2220 s->lossless_check_data = 0;
644
645 2220 memset(s->output_shift , 0, sizeof(s->output_shift ));
646 2220 memset(s->quant_step_size, 0, sizeof(s->quant_step_size));
647
648
2/2
✓ Branch 0 taken 5292 times.
✓ Branch 1 taken 2220 times.
 7512 for (ch = s->min_channel; ch <= s->max_channel; ch++) {
649 5292 ChannelParams *cp = &s->channel_params[ch];
650 5292 cp->filter_params[FIR].order = 0;
651 5292 cp->filter_params[IIR].order = 0;
652 5292 cp->filter_params[FIR].shift = 0;
653 5292 cp->filter_params[IIR].shift = 0;
654
655 /* Default audio coding is 24-bit raw PCM. */
656 5292 cp->huff_offset = 0;
657 5292 cp->sign_huff_offset = -(1 << 23);
658 5292 cp->codebook = 0;
659 5292 cp->huff_lsbs = 24;
660 }
661
662
2/2
✓ Branch 0 taken 1774 times.
✓ Branch 1 taken 446 times.
 2220 if (substr == m->max_decoded_substream) {
663 1774 av_channel_layout_uninit(&m->avctx->ch_layout);
664 1774 av_channel_layout_from_mask(&m->avctx->ch_layout, s->mask);
665 3548 m->pack_output = m->dsp.mlp_select_pack_output(s->ch_assign,
666 1774 s->output_shift,
667 1774 s->max_matrix_channel,
668 1774 m->avctx->sample_fmt == AV_SAMPLE_FMT_S32);
669
670
3/4
✓ Branch 0 taken 1117 times.
✓ Branch 1 taken 657 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 1117 times.
 1774 if (m->avctx->codec_id == AV_CODEC_ID_MLP && m->needs_reordering) {
671 ✗ if (s->mask == (AV_CH_LAYOUT_QUAD|AV_CH_LOW_FREQUENCY) ||
672 ✗ s->mask == AV_CH_LAYOUT_5POINT0_BACK) {
673 ✗ int i = s->ch_assign[4];
674 ✗ s->ch_assign[4] = s->ch_assign[3];
675 ✗ s->ch_assign[3] = s->ch_assign[2];
676 ✗ s->ch_assign[2] = i;
677 ✗ } else if (s->mask == AV_CH_LAYOUT_5POINT1_BACK) {
678 ✗ FFSWAP(int, s->ch_assign[2], s->ch_assign[4]);
679 ✗ FFSWAP(int, s->ch_assign[3], s->ch_assign[5]);
680 }
681 }
682
683 }
684
685 2220 return 0;
686 }
687
688 /** Read parameters for one of the prediction filters. */
689
690 5829 static int read_filter_params(MLPDecodeContext *m, GetBitContext *gbp,
691 unsigned int substr, unsigned int channel,
692 unsigned int filter)
693 {
694 5829 SubStream *s = &m->substream[substr];
695 5829 FilterParams *fp = &s->channel_params[channel].filter_params[filter];
696
2/2
✓ Branch 0 taken 525 times.
✓ Branch 1 taken 5304 times.
 5829 const int max_order = filter ? MAX_IIR_ORDER : MAX_FIR_ORDER;
697
2/2
✓ Branch 0 taken 525 times.
✓ Branch 1 taken 5304 times.
 5829 const char fchar = filter ? 'I' : 'F';
698 int i, order;
699
700 // Filter is 0 for FIR, 1 for IIR.
701
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 5829 times.
 5829 av_assert0(filter < 2);
702
703
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 5829 times.
 5829 if (m->filter_changed[channel][filter]++ > 1) {
704 ✗ av_log(m->avctx, AV_LOG_ERROR, "Filters may change only once per access unit.\n");
705 ✗ return AVERROR_INVALIDDATA;
706 }
707
708 5829 order = get_bits(gbp, 4);
709
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 5829 times.
 5829 if (order > max_order) {
710 ✗ av_log(m->avctx, AV_LOG_ERROR,
711 "%cIR filter order %d is greater than maximum %d.\n",
712 fchar, order, max_order);
713 ✗ return AVERROR_INVALIDDATA;
714 }
715 5829 fp->order = order;
716
717
2/2
✓ Branch 0 taken 5824 times.
✓ Branch 1 taken 5 times.
 5829 if (order > 0) {
718 5824 int32_t *fcoeff = s->channel_params[channel].coeff[filter];
719 int coeff_bits, coeff_shift;
720
721 5824 fp->shift = get_bits(gbp, 4);
722
723 5824 coeff_bits = get_bits(gbp, 5);
724 5824 coeff_shift = get_bits(gbp, 3);
725
2/4
✓ Branch 0 taken 5824 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 5824 times.
 5824 if (coeff_bits < 1 || coeff_bits > 16) {
726 ✗ av_log(m->avctx, AV_LOG_ERROR,
727 "%cIR filter coeff_bits must be between 1 and 16.\n",
728 fchar);
729 ✗ return AVERROR_INVALIDDATA;
730 }
731
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 5824 times.
 5824 if (coeff_bits + coeff_shift > 16) {
732 ✗ av_log(m->avctx, AV_LOG_ERROR,
733 "Sum of coeff_bits and coeff_shift for %cIR filter must be 16 or less.\n",
734 fchar);
735 ✗ return AVERROR_INVALIDDATA;
736 }
737
738
2/2
✓ Branch 0 taken 24259 times.
✓ Branch 1 taken 5824 times.
 30083 for (i = 0; i < order; i++)
739 24259 fcoeff[i] = get_sbits(gbp, coeff_bits) * (1 << coeff_shift);
740
741
2/2
✓ Branch 1 taken 525 times.
✓ Branch 2 taken 5299 times.
 5824 if (get_bits1(gbp)) {
742 int state_bits, state_shift;
743
744
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 525 times.
 525 if (filter == FIR) {
745 ✗ av_log(m->avctx, AV_LOG_ERROR,
746 "FIR filter has state data specified.\n");
747 ✗ return AVERROR_INVALIDDATA;
748 }
749
750 525 state_bits = get_bits(gbp, 4);
751 525 state_shift = get_bits(gbp, 4);
752
753 /* TODO: Check validity of state data. */
754
755
2/2
✓ Branch 0 taken 972 times.
✓ Branch 1 taken 525 times.
 1497 for (i = 0; i < order; i++)
756
1/2
✓ Branch 0 taken 972 times.
✗ Branch 1 not taken.
 972 fp->state[i] = state_bits ? get_sbits(gbp, state_bits) * (1 << state_shift) : 0;
757 }
758 }
759
760 5829 return 0;
761 }
762
763 /** Read parameters for primitive matrices. */
764
765 1927 static int read_matrix_params(MLPDecodeContext *m, unsigned int substr, GetBitContext *gbp)
766 {
767 1927 SubStream *s = &m->substream[substr];
768 unsigned int mat, ch;
769 3854 const int max_primitive_matrices = m->avctx->codec_id == AV_CODEC_ID_MLP
770 ? MAX_MATRICES_MLP
771
2/2
✓ Branch 0 taken 832 times.
✓ Branch 1 taken 1095 times.
 1927 : MAX_MATRICES_TRUEHD;
772
773
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 1927 times.
 1927 if (m->matrix_changed++ > 1) {
774 ✗ av_log(m->avctx, AV_LOG_ERROR, "Matrices may change only once per access unit.\n");
775 ✗ return AVERROR_INVALIDDATA;
776 }
777
778 1927 s->num_primitive_matrices = get_bits(gbp, 4);
779
780
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 1927 times.
 1927 if (s->num_primitive_matrices > max_primitive_matrices) {
781 ✗ av_log(m->avctx, AV_LOG_ERROR,
782 "Number of primitive matrices cannot be greater than %d.\n",
783 max_primitive_matrices);
784 ✗ goto error;
785 }
786
787
2/2
✓ Branch 0 taken 3494 times.
✓ Branch 1 taken 1927 times.
 5421 for (mat = 0; mat < s->num_primitive_matrices; mat++) {
788 int frac_bits, max_chan;
789 3494 s->matrix_out_ch[mat] = get_bits(gbp, 4);
790 3494 frac_bits = get_bits(gbp, 4);
791 3494 s->lsb_bypass [mat] = get_bits1(gbp);
792
793
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 3494 times.
 3494 if (s->matrix_out_ch[mat] > s->max_matrix_channel) {
794 ✗ av_log(m->avctx, AV_LOG_ERROR,
795 "Invalid channel %d specified as output from matrix.\n",
796 ✗ s->matrix_out_ch[mat]);
797 ✗ goto error;
798 }
799
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 3494 times.
 3494 if (frac_bits > 14) {
800 ✗ av_log(m->avctx, AV_LOG_ERROR,
801 "Too many fractional bits specified.\n");
802 ✗ goto error;
803 }
804
805 3494 max_chan = s->max_matrix_channel;
806
2/2
✓ Branch 0 taken 2146 times.
✓ Branch 1 taken 1348 times.
 3494 if (!s->noise_type)
807 2146 max_chan+=2;
808
809
2/2
✓ Branch 0 taken 16736 times.
✓ Branch 1 taken 3494 times.
 20230 for (ch = 0; ch <= max_chan; ch++) {
810 16736 int coeff_val = 0;
811
2/2
✓ Branch 1 taken 9693 times.
✓ Branch 2 taken 7043 times.
 16736 if (get_bits1(gbp))
812 9693 coeff_val = get_sbits(gbp, frac_bits + 2);
813
814 16736 s->matrix_coeff[mat][ch] = coeff_val * (1 << (14 - frac_bits));
815 }
816
817
2/2
✓ Branch 0 taken 1348 times.
✓ Branch 1 taken 2146 times.
 3494 if (s->noise_type)
818 1348 s->matrix_noise_shift[mat] = get_bits(gbp, 4);
819 else
820 2146 s->matrix_noise_shift[mat] = 0;
821 }
822
823 1927 return 0;
824 ✗ error:
825 ✗ s->num_primitive_matrices = 0;
826 ✗ memset(s->matrix_out_ch, 0, sizeof(s->matrix_out_ch));
827
828 ✗ return AVERROR_INVALIDDATA;
829 }
830
831 /** Read channel parameters. */
832
833 15911 static int read_channel_params(MLPDecodeContext *m, unsigned int substr,
834 GetBitContext *gbp, unsigned int ch)
835 {
836 15911 SubStream *s = &m->substream[substr];
837 15911 ChannelParams *cp = &s->channel_params[ch];
838 15911 FilterParams *fir = &cp->filter_params[FIR];
839 15911 FilterParams *iir = &cp->filter_params[IIR];
840 int ret;
841
842
1/2
✓ Branch 0 taken 15911 times.
✗ Branch 1 not taken.
 15911 if (s->param_presence_flags & PARAM_FIR)
843
2/2
✓ Branch 1 taken 5304 times.
✓ Branch 2 taken 10607 times.
 15911 if (get_bits1(gbp))
844
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 5304 times.
 5304 if ((ret = read_filter_params(m, gbp, substr, ch, FIR)) < 0)
845 ✗ return ret;
846
847
1/2
✓ Branch 0 taken 15911 times.
✗ Branch 1 not taken.
 15911 if (s->param_presence_flags & PARAM_IIR)
848
2/2
✓ Branch 1 taken 525 times.
✓ Branch 2 taken 15386 times.
 15911 if (get_bits1(gbp))
849
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 525 times.
 525 if ((ret = read_filter_params(m, gbp, substr, ch, IIR)) < 0)
850 ✗ return ret;
851
852
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 15911 times.
 15911 if (fir->order + iir->order > 8) {
853 ✗ av_log(m->avctx, AV_LOG_ERROR, "Total filter orders too high.\n");
854 ✗ return AVERROR_INVALIDDATA;
855 }
856
857
4/4
✓ Branch 0 taken 10614 times.
✓ Branch 1 taken 5297 times.
✓ Branch 2 taken 990 times.
✓ Branch 3 taken 9624 times.
 15911 if (fir->order && iir->order &&
858
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 990 times.
 990 fir->shift != iir->shift) {
859 ✗ av_log(m->avctx, AV_LOG_ERROR,
860 "FIR and IIR filters must use the same precision.\n");
861 ✗ return AVERROR_INVALIDDATA;
862 }
863 /* The FIR and IIR filters must have the same precision.
864 * To simplify the filtering code, only the precision of the
865 * FIR filter is considered. If only the IIR filter is employed,
866 * the FIR filter precision is set to that of the IIR filter, so
867 * that the filtering code can use it. */
868
3/4
✓ Branch 0 taken 5297 times.
✓ Branch 1 taken 10614 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 5297 times.
 15911 if (!fir->order && iir->order)
869 ✗ fir->shift = iir->shift;
870
871
1/2
✓ Branch 0 taken 15911 times.
✗ Branch 1 not taken.
 15911 if (s->param_presence_flags & PARAM_HUFFOFFSET)
872
2/2
✓ Branch 1 taken 3516 times.
✓ Branch 2 taken 12395 times.
 15911 if (get_bits1(gbp))
873 3516 cp->huff_offset = get_sbits(gbp, 15);
874
875 15911 cp->codebook = get_bits(gbp, 2);
876 15911 cp->huff_lsbs = get_bits(gbp, 5);
877
878
3/4
✓ Branch 0 taken 11899 times.
✓ Branch 1 taken 4012 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 11899 times.
 15911 if (cp->codebook > 0 && cp->huff_lsbs > 24) {
879 ✗ av_log(m->avctx, AV_LOG_ERROR, "Invalid huff_lsbs.\n");
880 ✗ cp->huff_lsbs = 0;
881 ✗ return AVERROR_INVALIDDATA;
882 }
883
884 15911 return 0;
885 }
886
887 /** Read decoding parameters that change more often than those in the restart
888 * header. */
889
890 16036 static int read_decoding_params(MLPDecodeContext *m, GetBitContext *gbp,
891 unsigned int substr)
892 {
893 16036 SubStream *s = &m->substream[substr];
894 unsigned int ch;
895 16036 int ret = 0;
896 16036 unsigned recompute_sho = 0;
897
898
1/2
✓ Branch 0 taken 16036 times.
✗ Branch 1 not taken.
 16036 if (s->param_presence_flags & PARAM_PRESENCE)
899
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 16036 times.
 16036 if (get_bits1(gbp))
900 ✗ s->param_presence_flags = get_bits(gbp, 8);
901
902
1/2
✓ Branch 0 taken 16036 times.
✗ Branch 1 not taken.
 16036 if (s->param_presence_flags & PARAM_BLOCKSIZE)
903
2/2
✓ Branch 1 taken 4422 times.
✓ Branch 2 taken 11614 times.
 16036 if (get_bits1(gbp)) {
904 4422 s->blocksize = get_bits(gbp, 9);
905
2/4
✓ Branch 0 taken 4422 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 4422 times.
 4422 if (s->blocksize < 8 || s->blocksize > m->access_unit_size) {
906 ✗ av_log(m->avctx, AV_LOG_ERROR, "Invalid blocksize.\n");
907 ✗ s->blocksize = 0;
908 ✗ return AVERROR_INVALIDDATA;
909 }
910 }
911
912
1/2
✓ Branch 0 taken 16036 times.
✗ Branch 1 not taken.
 16036 if (s->param_presence_flags & PARAM_MATRIX)
913
2/2
✓ Branch 1 taken 1927 times.
✓ Branch 2 taken 14109 times.
 16036 if (get_bits1(gbp))
914
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 1927 times.
 1927 if ((ret = read_matrix_params(m, substr, gbp)) < 0)
915 ✗ return ret;
916
917
1/2
✓ Branch 0 taken 16036 times.
✗ Branch 1 not taken.
 16036 if (s->param_presence_flags & PARAM_OUTSHIFT)
918
2/2
✓ Branch 1 taken 1403 times.
✓ Branch 2 taken 14633 times.
 16036 if (get_bits1(gbp)) {
919
2/2
✓ Branch 0 taken 4558 times.
✓ Branch 1 taken 1403 times.
 5961 for (ch = 0; ch <= s->max_matrix_channel; ch++) {
920 4558 s->output_shift[ch] = get_sbits(gbp, 4);
921
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 4558 times.
 4558 if (s->output_shift[ch] < 0) {
922 ✗ avpriv_request_sample(m->avctx, "Negative output_shift");
923 ✗ s->output_shift[ch] = 0;
924 }
925 }
926
2/2
✓ Branch 0 taken 1371 times.
✓ Branch 1 taken 32 times.
 1403 if (substr == m->max_decoded_substream)
927 1371 m->pack_output = m->dsp.mlp_select_pack_output(s->ch_assign,
928 1371 s->output_shift,
929 1371 s->max_matrix_channel,
930 1371 m->avctx->sample_fmt == AV_SAMPLE_FMT_S32);
931 }
932
933
1/2
✓ Branch 0 taken 16036 times.
✗ Branch 1 not taken.
 16036 if (s->param_presence_flags & PARAM_QUANTSTEP)
934
2/2
✓ Branch 1 taken 2236 times.
✓ Branch 2 taken 13800 times.
 16036 if (get_bits1(gbp))
935
2/2
✓ Branch 0 taken 6192 times.
✓ Branch 1 taken 2236 times.
 8428 for (ch = 0; ch <= s->max_channel; ch++) {
936 6192 s->quant_step_size[ch] = get_bits(gbp, 4);
937
938 6192 recompute_sho |= 1<<ch;
939 }
940
941
2/2
✓ Branch 0 taken 38080 times.
✓ Branch 1 taken 16036 times.
 54116 for (ch = s->min_channel; ch <= s->max_channel; ch++)
942
2/2
✓ Branch 1 taken 15911 times.
✓ Branch 2 taken 22169 times.
 38080 if (get_bits1(gbp)) {
943 15911 recompute_sho |= 1<<ch;
944
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 15911 times.
 15911 if ((ret = read_channel_params(m, substr, gbp, ch)) < 0)
945 ✗ goto fail;
946 }
947
948
949 16036 fail:
950
2/2
✓ Branch 0 taken 44382 times.
✓ Branch 1 taken 16036 times.
 60418 for (ch = 0; ch <= s->max_channel; ch++) {
951
2/2
✓ Branch 0 taken 16834 times.
✓ Branch 1 taken 27548 times.
 44382 if (recompute_sho & (1<<ch)) {
952 16834 ChannelParams *cp = &s->channel_params[ch];
953
954
3/4
✓ Branch 0 taken 11899 times.
✓ Branch 1 taken 4935 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 11899 times.
 16834 if (cp->codebook > 0 && cp->huff_lsbs < s->quant_step_size[ch]) {
955 ✗ if (ret >= 0) {
956 ✗ av_log(m->avctx, AV_LOG_ERROR, "quant_step_size larger than huff_lsbs\n");
957 ✗ ret = AVERROR_INVALIDDATA;
958 }
959 ✗ s->quant_step_size[ch] = 0;
960 }
961
962 16834 cp->sign_huff_offset = calculate_sign_huff(m, substr, ch);
963 }
964 }
965 16036 return ret;
966 }
967
968 #define MSB_MASK(bits) (-(1 << (bits)))
969
970 /** Generate PCM samples using the prediction filters and residual values
971 * read from the data stream, and update the filter state. */
972
973 73644 static void filter_channel(MLPDecodeContext *m, unsigned int substr,
974 unsigned int channel)
975 {
976 73644 SubStream *s = &m->substream[substr];
977 73644 const int32_t *fircoeff = s->channel_params[channel].coeff[FIR];
978 int32_t state_buffer[NUM_FILTERS][MAX_BLOCKSIZE + MAX_FIR_ORDER];
979 73644 int32_t *firbuf = state_buffer[FIR] + MAX_BLOCKSIZE;
980 73644 int32_t *iirbuf = state_buffer[IIR] + MAX_BLOCKSIZE;
981 73644 FilterParams *fir = &s->channel_params[channel].filter_params[FIR];
982 73644 FilterParams *iir = &s->channel_params[channel].filter_params[IIR];
983 73644 unsigned int filter_shift = fir->shift;
984 73644 int32_t mask = MSB_MASK(s->quant_step_size[channel]);
985
986 73644 memcpy(firbuf, fir->state, MAX_FIR_ORDER * sizeof(int32_t));
987 73644 memcpy(iirbuf, iir->state, MAX_IIR_ORDER * sizeof(int32_t));
988
989 73644 m->dsp.mlp_filter_channel(firbuf, fircoeff,
990 73644 fir->order, iir->order,
991 73644 filter_shift, mask, s->blocksize,
992 73644 &m->sample_buffer[s->blockpos][channel]);
993
994 73644 memcpy(fir->state, firbuf - s->blocksize, MAX_FIR_ORDER * sizeof(int32_t));
995 73644 memcpy(iir->state, iirbuf - s->blocksize, MAX_IIR_ORDER * sizeof(int32_t));
996 73644 }
997
998 /** Read a block of PCM residual data (or actual if no filtering active). */
999
1000 30249 static int read_block_data(MLPDecodeContext *m, GetBitContext *gbp,
1001 unsigned int substr)
1002 {
1003 30249 SubStream *s = &m->substream[substr];
1004 30249 unsigned int i, ch, expected_stream_pos = 0;
1005 int ret;
1006
1007
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 30249 times.
 30249 if (s->data_check_present) {
1008 ✗ expected_stream_pos = get_bits_count(gbp);
1009 ✗ expected_stream_pos += get_bits(gbp, 16);
1010 ✗ avpriv_request_sample(m->avctx,
1011 "Substreams with VLC block size check info");
1012 }
1013
1014
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 30249 times.
 30249 if (s->blockpos + s->blocksize > m->access_unit_size) {
1015 ✗ av_log(m->avctx, AV_LOG_ERROR, "too many audio samples in frame\n");
1016 ✗ return AVERROR_INVALIDDATA;
1017 }
1018
1019 30249 memset(&m->bypassed_lsbs[s->blockpos][0], 0,
1020 30249 s->blocksize * sizeof(m->bypassed_lsbs[0]));
1021
1022
2/2
✓ Branch 0 taken 1121160 times.
✓ Branch 1 taken 30249 times.
 1151409 for (i = 0; i < s->blocksize; i++)
1023
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 1121160 times.
 1121160 if ((ret = read_huff_channels(m, gbp, substr, i)) < 0)
1024 ✗ return ret;
1025
1026
2/2
✓ Branch 0 taken 73644 times.
✓ Branch 1 taken 30249 times.
 103893 for (ch = s->min_channel; ch <= s->max_channel; ch++)
1027 73644 filter_channel(m, substr, ch);
1028
1029 30249 s->blockpos += s->blocksize;
1030
1031
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 30249 times.
 30249 if (s->data_check_present) {
1032 ✗ if (get_bits_count(gbp) != expected_stream_pos)
1033 ✗ av_log(m->avctx, AV_LOG_ERROR, "block data length mismatch\n");
1034 ✗ skip_bits(gbp, 8);
1035 }
1036
1037 30249 return 0;
1038 }
1039
1040 /** Data table used for TrueHD noise generation function. */
1041
1042 static const int8_t noise_table[256] = {
1043 30, 51, 22, 54, 3, 7, -4, 38, 14, 55, 46, 81, 22, 58, -3, 2,
1044 52, 31, -7, 51, 15, 44, 74, 30, 85, -17, 10, 33, 18, 80, 28, 62,
1045 10, 32, 23, 69, 72, 26, 35, 17, 73, 60, 8, 56, 2, 6, -2, -5,
1046 51, 4, 11, 50, 66, 76, 21, 44, 33, 47, 1, 26, 64, 48, 57, 40,
1047 38, 16, -10, -28, 92, 22, -18, 29, -10, 5, -13, 49, 19, 24, 70, 34,
1048 61, 48, 30, 14, -6, 25, 58, 33, 42, 60, 67, 17, 54, 17, 22, 30,
1049 67, 44, -9, 50, -11, 43, 40, 32, 59, 82, 13, 49, -14, 55, 60, 36,
1050 48, 49, 31, 47, 15, 12, 4, 65, 1, 23, 29, 39, 45, -2, 84, 69,
1051 0, 72, 37, 57, 27, 41, -15, -16, 35, 31, 14, 61, 24, 0, 27, 24,
1052 16, 41, 55, 34, 53, 9, 56, 12, 25, 29, 53, 5, 20, -20, -8, 20,
1053 13, 28, -3, 78, 38, 16, 11, 62, 46, 29, 21, 24, 46, 65, 43, -23,
1054 89, 18, 74, 21, 38, -12, 19, 12, -19, 8, 15, 33, 4, 57, 9, -8,
1055 36, 35, 26, 28, 7, 83, 63, 79, 75, 11, 3, 87, 37, 47, 34, 40,
1056 39, 19, 20, 42, 27, 34, 39, 77, 13, 42, 59, 64, 45, -1, 32, 37,
1057 45, -5, 53, -6, 7, 36, 50, 23, 6, 32, 9, -21, 18, 71, 27, 52,
1058 -25, 31, 35, 42, -1, 68, 63, 52, 26, 43, 66, 37, 41, 25, 40, 70,
1059 };
1060
1061 /** Noise generation functions.
1062 * I'm not sure what these are for - they seem to be some kind of pseudorandom
1063 * sequence generators, used to generate noise data which is used when the
1064 * channels are rematrixed. I'm not sure if they provide a practical benefit
1065 * to compression, or just obfuscate the decoder. Are they for some kind of
1066 * dithering? */
1067
1068 /** Generate two channels of noise, used in the matrix when
1069 * restart sync word == 0x31ea. */
1070
1071 12380 static void generate_2_noise_channels(MLPDecodeContext *m, unsigned int substr)
1072 {
1073 12380 SubStream *s = &m->substream[substr];
1074 unsigned int i;
1075 12380 uint32_t seed = s->noisegen_seed;
1076 12380 unsigned int maxchan = s->max_matrix_channel;
1077
1078
2/2
✓ Branch 0 taken 495200 times.
✓ Branch 1 taken 12380 times.
 507580 for (i = 0; i < s->blockpos; i++) {
1079 495200 uint16_t seed_shr7 = seed >> 7;
1080 495200 m->sample_buffer[i][maxchan+1] = ((int8_t)(seed >> 15)) * (1 << s->noise_shift);
1081 495200 m->sample_buffer[i][maxchan+2] = ((int8_t) seed_shr7) * (1 << s->noise_shift);
1082
1083 495200 seed = (seed << 16) ^ seed_shr7 ^ (seed_shr7 << 5);
1084 }
1085
1086 12380 s->noisegen_seed = seed;
1087 12380 }
1088
1089 /** Generate a block of noise, used when restart sync word == 0x31eb. */
1090
1091 7759 static void fill_noise_buffer(MLPDecodeContext *m, unsigned int substr)
1092 {
1093 7759 SubStream *s = &m->substream[substr];
1094 unsigned int i;
1095 7759 uint32_t seed = s->noisegen_seed;
1096
1097
2/2
✓ Branch 0 taken 496576 times.
✓ Branch 1 taken 7759 times.
 504335 for (i = 0; i < m->access_unit_size_pow2; i++) {
1098 496576 uint8_t seed_shr15 = seed >> 15;
1099 496576 m->noise_buffer[i] = noise_table[seed_shr15];
1100 496576 seed = (seed << 8) ^ seed_shr15 ^ (seed_shr15 << 5);
1101 }
1102
1103 7759 s->noisegen_seed = seed;
1104 7759 }
1105
1106 /** Write the audio data into the output buffer. */
1107
1108 20139 static int output_data(MLPDecodeContext *m, unsigned int substr,
1109 AVFrame *frame, int *got_frame_ptr)
1110 {
1111 20139 AVCodecContext *avctx = m->avctx;
1112 20139 SubStream *s = &m->substream[substr];
1113 unsigned int mat;
1114 unsigned int maxchan;
1115 int ret;
1116 20139 int is32 = (m->avctx->sample_fmt == AV_SAMPLE_FMT_S32);
1117
1118
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 20139 times.
 20139 if (m->avctx->ch_layout.nb_channels != s->max_matrix_channel + 1) {
1119 ✗ av_log(m->avctx, AV_LOG_ERROR, "channel count mismatch\n");
1120 ✗ return AVERROR_INVALIDDATA;
1121 }
1122
1123
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 20139 times.
 20139 if (!s->blockpos) {
1124 ✗ av_log(avctx, AV_LOG_ERROR, "No samples to output.\n");
1125 ✗ return AVERROR_INVALIDDATA;
1126 }
1127
1128 20139 maxchan = s->max_matrix_channel;
1129
2/2
✓ Branch 0 taken 12380 times.
✓ Branch 1 taken 7759 times.
 20139 if (!s->noise_type) {
1130 12380 generate_2_noise_channels(m, substr);
1131 12380 maxchan += 2;
1132 } else {
1133 7759 fill_noise_buffer(m, substr);
1134 }
1135
1136 /* Apply the channel matrices in turn to reconstruct the original audio
1137 * samples. */
1138
2/2
✓ Branch 0 taken 35059 times.
✓ Branch 1 taken 20139 times.
 55198 for (mat = 0; mat < s->num_primitive_matrices; mat++) {
1139 35059 unsigned int dest_ch = s->matrix_out_ch[mat];
1140 35059 m->dsp.mlp_rematrix_channel(&m->sample_buffer[0][0],
1141 35059 s->matrix_coeff[mat],
1142 35059 &m->bypassed_lsbs[0][mat],
1143 35059 m->noise_buffer,
1144 35059 s->num_primitive_matrices - mat,
1145 dest_ch,
1146 35059 s->blockpos,
1147 maxchan,
1148 35059 s->matrix_noise_shift[mat],
1149 m->access_unit_size_pow2,
1150 35059 MSB_MASK(s->quant_step_size[dest_ch]));
1151 }
1152
1153 /* get output buffer */
1154 20139 frame->nb_samples = s->blockpos;
1155
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 20139 times.
 20139 if ((ret = ff_get_buffer(avctx, frame, 0)) < 0)
1156 ✗ return ret;
1157 40278 s->lossless_check_data = m->pack_output(s->lossless_check_data,
1158 20139 s->blockpos,
1159 20139 m->sample_buffer,
1160 20139 frame->data[0],
1161 20139 s->ch_assign,
1162 20139 s->output_shift,
1163 20139 s->max_matrix_channel,
1164 is32);
1165
1166 /* Update matrix encoding side data */
1167
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 20139 times.
 20139 if (s->matrix_encoding != s->prev_matrix_encoding) {
1168 ✗ if ((ret = ff_side_data_update_matrix_encoding(frame, s->matrix_encoding)) < 0)
1169 ✗ return ret;
1170
1171 ✗ s->prev_matrix_encoding = s->matrix_encoding;
1172 }
1173
1174 20139 *got_frame_ptr = 1;
1175
1176 20139 return 0;
1177 }
1178
1179 /** Read an access unit from the stream.
1180 * @return negative on error, 0 if not enough data is present in the input stream,
1181 * otherwise the number of bytes consumed. */
1182
1183 20139 static int read_access_unit(AVCodecContext *avctx, AVFrame *frame,
1184 int *got_frame_ptr, AVPacket *avpkt)
1185 {
1186 20139 const uint8_t *buf = avpkt->data;
1187 20139 int buf_size = avpkt->size;
1188 20139 MLPDecodeContext *m = avctx->priv_data;
1189 GetBitContext gb;
1190 unsigned int length, substr;
1191 unsigned int substream_start;
1192 20139 unsigned int header_size = 4;
1193 20139 unsigned int substr_header_size = 0;
1194 uint8_t substream_parity_present[MAX_SUBSTREAMS];
1195 uint16_t substream_data_len[MAX_SUBSTREAMS];
1196 uint8_t parity_bits;
1197 int ret;
1198
1199
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 20139 times.
 20139 if (buf_size < 4)
1200 ✗ return AVERROR_INVALIDDATA;
1201
1202 20139 length = (AV_RB16(buf) & 0xfff) * 2;
1203
1204
2/4
✓ Branch 0 taken 20139 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 20139 times.
 20139 if (length < 4 || length > buf_size)
1205 ✗ return AVERROR_INVALIDDATA;
1206
1207 20139 init_get_bits(&gb, (buf + 4), (length - 4) * 8);
1208
1209 20139 m->is_major_sync_unit = 0;
1210
2/2
✓ Branch 1 taken 1774 times.
✓ Branch 2 taken 18365 times.
 20139 if (show_bits_long(&gb, 31) == (0xf8726fba >> 1)) {
1211
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 1774 times.
 1774 if (read_major_sync(m, &gb) < 0)
1212 ✗ goto error;
1213 1774 m->is_major_sync_unit = 1;
1214 1774 header_size += m->major_sync_header_size;
1215 }
1216
1217
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 20139 times.
 20139 if (!m->params_valid) {
1218 ✗ av_log(m->avctx, AV_LOG_WARNING,
1219 "Stream parameters not seen; skipping frame.\n");
1220 ✗ *got_frame_ptr = 0;
1221 ✗ return length;
1222 }
1223
1224 20139 substream_start = 0;
1225
1226
2/2
✓ Branch 0 taken 31571 times.
✓ Branch 1 taken 20139 times.
 51710 for (substr = 0; substr < m->num_substreams; substr++) {
1227 int extraword_present, checkdata_present, end, nonrestart_substr;
1228
1229 31571 extraword_present = get_bits1(&gb);
1230 31571 nonrestart_substr = get_bits1(&gb);
1231 31571 checkdata_present = get_bits1(&gb);
1232 31571 skip_bits1(&gb);
1233
1234 31571 end = get_bits(&gb, 12) * 2;
1235
1236 31571 substr_header_size += 2;
1237
1238
2/2
✓ Branch 0 taken 176 times.
✓ Branch 1 taken 31395 times.
 31571 if (extraword_present) {
1239
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 176 times.
 176 if (m->avctx->codec_id == AV_CODEC_ID_MLP) {
1240 ✗ av_log(m->avctx, AV_LOG_ERROR, "There must be no extraword for MLP.\n");
1241 ✗ goto error;
1242 }
1243 176 skip_bits(&gb, 16);
1244 176 substr_header_size += 2;
1245 }
1246
1247
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 31571 times.
 31571 if (length < header_size + substr_header_size) {
1248 ✗ av_log(m->avctx, AV_LOG_ERROR, "Insufficient data for headers\n");
1249 ✗ goto error;
1250 }
1251
1252
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 31571 times.
 31571 if (!(nonrestart_substr ^ m->is_major_sync_unit)) {
1253 ✗ av_log(m->avctx, AV_LOG_ERROR, "Invalid nonrestart_substr.\n");
1254 ✗ goto error;
1255 }
1256
1257
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 31571 times.
 31571 if (end + header_size + substr_header_size > length) {
1258 ✗ av_log(m->avctx, AV_LOG_ERROR,
1259 "Indicated length of substream %d data goes off end of "
1260 "packet.\n", substr);
1261 ✗ end = length - header_size - substr_header_size;
1262 }
1263
1264
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 31571 times.
 31571 if (end < substream_start) {
1265 ✗ av_log(avctx, AV_LOG_ERROR,
1266 "Indicated end offset of substream %d data "
1267 "is smaller than calculated start offset.\n",
1268 substr);
1269 ✗ goto error;
1270 }
1271
1272
2/2
✓ Branch 0 taken 3327 times.
✓ Branch 1 taken 28244 times.
 31571 if (substr > m->max_decoded_substream)
1273 3327 continue;
1274
1275 28244 substream_parity_present[substr] = checkdata_present;
1276 28244 substream_data_len[substr] = end - substream_start;
1277 28244 substream_start = end;
1278 }
1279
1280 20139 parity_bits = ff_mlp_calculate_parity(buf, 4);
1281 20139 parity_bits ^= ff_mlp_calculate_parity(buf + header_size, substr_header_size);
1282
1283
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 20139 times.
 20139 if ((((parity_bits >> 4) ^ parity_bits) & 0xF) != 0xF) {
1284 ✗ av_log(avctx, AV_LOG_ERROR, "Parity check failed.\n");
1285 ✗ goto error;
1286 }
1287
1288 20139 buf += header_size + substr_header_size;
1289
1290
2/2
✓ Branch 0 taken 28029 times.
✓ Branch 1 taken 20139 times.
 48168 for (substr = 0; substr <= m->max_decoded_substream; substr++) {
1291 28029 SubStream *s = &m->substream[substr];
1292
1293 28029 init_get_bits(&gb, buf, substream_data_len[substr] * 8);
1294
1295 28029 m->matrix_changed = 0;
1296 28029 memset(m->filter_changed, 0, sizeof(m->filter_changed));
1297
1298 28029 s->blockpos = 0;
1299 do {
1300
2/2
✓ Branch 1 taken 16036 times.
✓ Branch 2 taken 14213 times.
 30249 if (get_bits1(&gb)) {
1301
2/2
✓ Branch 1 taken 2220 times.
✓ Branch 2 taken 13816 times.
 16036 if (get_bits1(&gb)) {
1302 /* A restart header should be present. */
1303
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 2220 times.
 2220 if (read_restart_header(m, &gb, buf, substr) < 0)
1304 ✗ goto next_substr;
1305 2220 s->restart_seen = 1;
1306 }
1307
1308
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 16036 times.
 16036 if (!s->restart_seen)
1309 ✗ goto next_substr;
1310
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 16036 times.
 16036 if (read_decoding_params(m, &gb, substr) < 0)
1311 ✗ goto next_substr;
1312 }
1313
1314
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 30249 times.
 30249 if (!s->restart_seen)
1315 ✗ goto next_substr;
1316
1317
2/2
✓ Branch 0 taken 14504 times.
✓ Branch 1 taken 15745 times.
 30249 if (((avctx->ch_layout.nb_channels == 6 &&
1318
1/2
✓ Branch 0 taken 14504 times.
✗ Branch 1 not taken.
 14504 ((m->substream_info >> 2) & 0x3) != 0x3) ||
1319
2/2
✓ Branch 0 taken 405 times.
✓ Branch 1 taken 29844 times.
 30249 (avctx->ch_layout.nb_channels == 8 &&
1320
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 405 times.
 405 ((m->substream_info >> 4) & 0x7) != 0x7 &&
1321 ✗ ((m->substream_info >> 4) & 0x7) != 0x6 &&
1322 ✗ ((m->substream_info >> 4) & 0x7) != 0x3)) &&
1323 ✗ substr > 0 && substr < m->max_decoded_substream &&
1324 ✗ (s->min_channel <= m->substream[substr - 1].max_channel)) {
1325 ✗ av_log(avctx, AV_LOG_DEBUG,
1326 "Previous substream(%d) channels overlaps current substream(%d) channels, skipping.\n",
1327 substr - 1, substr);
1328 ✗ goto next_substr;
1329 }
1330
1331
2/2
✓ Branch 0 taken 8336 times.
✓ Branch 1 taken 21913 times.
 30249 if (substr != m->max_decoded_substream &&
1332
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 8336 times.
 8336 ((s->coded_channels & m->substream[m->max_decoded_substream].coded_channels) != 0)) {
1333 ✗ av_log(avctx, AV_LOG_DEBUG,
1334 "Current substream(%d) channels overlaps final substream(%d) channels, skipping.\n",
1335 ✗ substr, m->max_decoded_substream);
1336 ✗ goto next_substr;
1337 }
1338
1339
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 30249 times.
 30249 if ((ret = read_block_data(m, &gb, substr)) < 0)
1340 ✗ return ret;
1341
1342
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 30249 times.
 30249 if (get_bits_count(&gb) >= substream_data_len[substr] * 8)
1343 ✗ goto substream_length_mismatch;
1344
1345
2/2
✓ Branch 1 taken 2220 times.
✓ Branch 2 taken 28029 times.
 30249 } while (!get_bits1(&gb));
1346
1347 28029 skip_bits(&gb, (-get_bits_count(&gb)) & 15);
1348
1349
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 28029 times.
 28029 if (substream_data_len[substr] * 8 - get_bits_count(&gb) >= 32) {
1350 int shorten_by;
1351
1352 ✗ if (get_bits(&gb, 16) != 0xD234)
1353 ✗ return AVERROR_INVALIDDATA;
1354
1355 ✗ shorten_by = get_bits(&gb, 16);
1356 ✗ if (m->avctx->codec_id == AV_CODEC_ID_TRUEHD && shorten_by & 0x2000)
1357 ✗ s->blockpos -= FFMIN(shorten_by & 0x1FFF, s->blockpos);
1358 ✗ else if (m->avctx->codec_id == AV_CODEC_ID_MLP && shorten_by != 0xD234)
1359 ✗ return AVERROR_INVALIDDATA;
1360
1361 ✗ av_log(m->avctx, AV_LOG_DEBUG, "End of stream indicated.\n");
1362 ✗ s->end_of_stream = 1;
1363 }
1364
1365
1/2
✓ Branch 0 taken 28029 times.
✗ Branch 1 not taken.
 28029 if (substream_parity_present[substr]) {
1366 uint8_t parity, checksum;
1367
1368
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 28029 times.
 28029 if (substream_data_len[substr] * 8 - get_bits_count(&gb) != 16)
1369 ✗ goto substream_length_mismatch;
1370
1371 28029 parity = ff_mlp_calculate_parity(buf, substream_data_len[substr] - 2);
1372 28029 checksum = ff_mlp_checksum8 (buf, substream_data_len[substr] - 2);
1373
1374
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 28029 times.
 28029 if ((get_bits(&gb, 8) ^ parity) != 0xa9 )
1375 ✗ av_log(m->avctx, AV_LOG_ERROR, "Substream %d parity check failed.\n", substr);
1376
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 28029 times.
 28029 if ( get_bits(&gb, 8) != checksum)
1377 ✗ av_log(m->avctx, AV_LOG_ERROR, "Substream %d checksum failed.\n" , substr);
1378 }
1379
1380
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 28029 times.
 28029 if (substream_data_len[substr] * 8 != get_bits_count(&gb))
1381 ✗ goto substream_length_mismatch;
1382
1383 28029 next_substr:
1384
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 28029 times.
 28029 if (!s->restart_seen)
1385 ✗ av_log(m->avctx, AV_LOG_ERROR,
1386 "No restart header present in substream %d.\n", substr);
1387
1388 28029 buf += substream_data_len[substr];
1389 }
1390
1391
1/2
✗ Branch 1 not taken.
✓ Branch 2 taken 20139 times.
 20139 if ((ret = output_data(m, m->max_decoded_substream, frame, got_frame_ptr)) < 0)
1392 ✗ return ret;
1393
1394
2/2
✓ Branch 0 taken 28029 times.
✓ Branch 1 taken 20139 times.
 48168 for (substr = 0; substr <= m->max_decoded_substream; substr++){
1395 28029 SubStream *s = &m->substream[substr];
1396
1397
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 28029 times.
 28029 if (s->end_of_stream) {
1398 ✗ s->lossless_check_data = 0xffffffff;
1399 ✗ s->end_of_stream = 0;
1400 ✗ m->params_valid = 0;
1401 }
1402 }
1403
1404 20139 return length;
1405
1406 ✗ substream_length_mismatch:
1407 ✗ av_log(m->avctx, AV_LOG_ERROR, "substream %d length mismatch\n", substr);
1408 ✗ return AVERROR_INVALIDDATA;
1409
1410 ✗ error:
1411 ✗ m->params_valid = 0;
1412 ✗ return AVERROR_INVALIDDATA;
1413 }
1414
1415 ✗ static void mlp_decode_flush(AVCodecContext *avctx)
1416 {
1417 ✗ MLPDecodeContext *m = avctx->priv_data;
1418
1419 ✗ m->params_valid = 0;
1420 ✗ for (int substr = 0; substr <= m->max_decoded_substream; substr++){
1421 ✗ SubStream *s = &m->substream[substr];
1422
1423 ✗ s->lossless_check_data = 0xffffffff;
1424 ✗ s->prev_matrix_encoding = 0;
1425 }
1426 ✗ }
1427
1428 #define OFFSET(x) offsetof(MLPDecodeContext, x)
1429 #define FLAGS (AV_OPT_FLAG_DECODING_PARAM | AV_OPT_FLAG_AUDIO_PARAM)
1430 static const AVOption options[] = {
1431 { "downmix", "Request a specific channel layout from the decoder", OFFSET(downmix_layout),
1432 AV_OPT_TYPE_CHLAYOUT, {.str = NULL}, .flags = FLAGS },
1433 { NULL },
1434 };
1435
1436 static const AVClass mlp_decoder_class = {
1437 .class_name = "MLP decoder",
1438 .item_name = av_default_item_name,
1439 .option = options,
1440 .version = LIBAVUTIL_VERSION_INT,
1441 };
1442
1443 static const AVClass truehd_decoder_class = {
1444 .class_name = "TrueHD decoder",
1445 .item_name = av_default_item_name,
1446 .option = options,
1447 .version = LIBAVUTIL_VERSION_INT,
1448 };
1449
1450 #if CONFIG_MLP_DECODER
1451 const FFCodec ff_mlp_decoder = {
1452 .p.name = "mlp",
1453 CODEC_LONG_NAME("MLP (Meridian Lossless Packing)"),
1454 .p.type = AVMEDIA_TYPE_AUDIO,
1455 .p.id = AV_CODEC_ID_MLP,
1456 .priv_data_size = sizeof(MLPDecodeContext),
1457 .p.priv_class = &mlp_decoder_class,
1458 .init = mlp_decode_init,
1459 FF_CODEC_DECODE_CB(read_access_unit),
1460 .flush = mlp_decode_flush,
1461 .p.capabilities = AV_CODEC_CAP_DR1 | AV_CODEC_CAP_CHANNEL_CONF,
1462 };
1463 #endif
1464 #if CONFIG_TRUEHD_DECODER
1465 const FFCodec ff_truehd_decoder = {
1466 .p.name = "truehd",
1467 CODEC_LONG_NAME("TrueHD"),
1468 .p.type = AVMEDIA_TYPE_AUDIO,
1469 .p.id = AV_CODEC_ID_TRUEHD,
1470 .priv_data_size = sizeof(MLPDecodeContext),
1471 .p.priv_class = &truehd_decoder_class,
1472 .init = mlp_decode_init,
1473 FF_CODEC_DECODE_CB(read_access_unit),
1474 .flush = mlp_decode_flush,
1475 .p.capabilities = AV_CODEC_CAP_DR1 | AV_CODEC_CAP_CHANNEL_CONF,
1476 .p.profiles = NULL_IF_CONFIG_SMALL(ff_truehd_profiles),
1477 };
1478 #endif /* CONFIG_TRUEHD_DECODER */
1479