floatdsp: move vector_fmul_add from dsputil to avfloatdsp.
[ffmpeg.git] / libavformat / mxfdec.c
1 /*
2  * MXF demuxer.
3  * Copyright (c) 2006 SmartJog S.A., Baptiste Coudurier <baptiste dot coudurier at smartjog dot com>
4  *
5  * This file is part of Libav.
6  *
7  * Libav 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  * Libav 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 Libav; if not, write to the Free Software
19  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
20  */
21
22 /*
23  * References
24  * SMPTE 336M KLV Data Encoding Protocol Using Key-Length-Value
25  * SMPTE 377M MXF File Format Specifications
26  * SMPTE 378M Operational Pattern 1a
27  * SMPTE 379M MXF Generic Container
28  * SMPTE 381M Mapping MPEG Streams into the MXF Generic Container
29  * SMPTE 382M Mapping AES3 and Broadcast Wave Audio into the MXF Generic Container
30  * SMPTE 383M Mapping DV-DIF Data to the MXF Generic Container
31  *
32  * Principle
33  * Search for Track numbers which will identify essence element KLV packets.
34  * Search for SourcePackage which define tracks which contains Track numbers.
35  * Material Package contains tracks with reference to SourcePackage tracks.
36  * Search for Descriptors (Picture, Sound) which contains codec info and parameters.
37  * Assign Descriptors to correct Tracks.
38  *
39  * Metadata reading functions read Local Tags, get InstanceUID(0x3C0A) then add MetaDataSet to MXFContext.
40  * Metadata parsing resolves Strong References to objects.
41  *
42  * Simple demuxer, only OP1A supported and some files might not work at all.
43  * Only tracks with associated descriptors will be decoded. "Highly Desirable" SMPTE 377M D.1
44  */
45
46 //#define DEBUG
47
48 #include "libavutil/aes.h"
49 #include "libavutil/mathematics.h"
50 #include "libavcodec/bytestream.h"
51 #include "avformat.h"
52 #include "internal.h"
53 #include "mxf.h"
54
55 typedef enum {
56     Header,
57     BodyPartition,
58     Footer
59 } MXFPartitionType;
60
61 typedef enum {
62     OP1a = 1,
63     OP1b,
64     OP1c,
65     OP2a,
66     OP2b,
67     OP2c,
68     OP3a,
69     OP3b,
70     OP3c,
71     OPAtom,
72     OPSonyOpt,  /* FATE sample, violates the spec in places */
73 } MXFOP;
74
75 typedef struct {
76     int closed;
77     int complete;
78     MXFPartitionType type;
79     uint64_t previous_partition;
80     int index_sid;
81     int body_sid;
82     int64_t this_partition;
83     int64_t essence_offset;         ///< absolute offset of essence
84     int64_t essence_length;
85     int32_t kag_size;
86     int64_t header_byte_count;
87     int64_t index_byte_count;
88     int pack_length;
89 } MXFPartition;
90
91 typedef struct {
92     UID uid;
93     enum MXFMetadataSetType type;
94     UID source_container_ul;
95 } MXFCryptoContext;
96
97 typedef struct {
98     UID uid;
99     enum MXFMetadataSetType type;
100     UID source_package_uid;
101     UID data_definition_ul;
102     int64_t duration;
103     int64_t start_position;
104     int source_track_id;
105 } MXFStructuralComponent;
106
107 typedef struct {
108     UID uid;
109     enum MXFMetadataSetType type;
110     UID data_definition_ul;
111     UID *structural_components_refs;
112     int structural_components_count;
113     int64_t duration;
114 } MXFSequence;
115
116 typedef struct {
117     UID uid;
118     enum MXFMetadataSetType type;
119     MXFSequence *sequence; /* mandatory, and only one */
120     UID sequence_ref;
121     int track_id;
122     uint8_t track_number[4];
123     AVRational edit_rate;
124     int intra_only;
125 } MXFTrack;
126
127 typedef struct {
128     UID uid;
129     enum MXFMetadataSetType type;
130     UID essence_container_ul;
131     UID essence_codec_ul;
132     AVRational sample_rate;
133     AVRational aspect_ratio;
134     int width;
135     int height; /* Field height, not frame height */
136     int frame_layout; /* See MXFFrameLayout enum */
137     int channels;
138     int bits_per_sample;
139     unsigned int component_depth;
140     unsigned int horiz_subsampling;
141     unsigned int vert_subsampling;
142     UID *sub_descriptors_refs;
143     int sub_descriptors_count;
144     int linked_track_id;
145     uint8_t *extradata;
146     int extradata_size;
147     enum AVPixelFormat pix_fmt;
148 } MXFDescriptor;
149
150 typedef struct {
151     UID uid;
152     enum MXFMetadataSetType type;
153     int edit_unit_byte_count;
154     int index_sid;
155     int body_sid;
156     AVRational index_edit_rate;
157     uint64_t index_start_position;
158     uint64_t index_duration;
159     int8_t *temporal_offset_entries;
160     int *flag_entries;
161     uint64_t *stream_offset_entries;
162     int nb_index_entries;
163 } MXFIndexTableSegment;
164
165 typedef struct {
166     UID uid;
167     enum MXFMetadataSetType type;
168     UID package_uid;
169     UID *tracks_refs;
170     int tracks_count;
171     MXFDescriptor *descriptor; /* only one */
172     UID descriptor_ref;
173 } MXFPackage;
174
175 typedef struct {
176     UID uid;
177     enum MXFMetadataSetType type;
178 } MXFMetadataSet;
179
180 /* decoded index table */
181 typedef struct {
182     int index_sid;
183     int body_sid;
184     int nb_ptses;               /* number of PTSes or total duration of index */
185     int64_t first_dts;          /* DTS = EditUnit + first_dts */
186     int64_t *ptses;             /* maps EditUnit -> PTS */
187     int nb_segments;
188     MXFIndexTableSegment **segments;    /* sorted by IndexStartPosition */
189     AVIndexEntry *fake_index;   /* used for calling ff_index_search_timestamp() */
190 } MXFIndexTable;
191
192 typedef struct {
193     MXFPartition *partitions;
194     unsigned partitions_count;
195     MXFOP op;
196     UID *packages_refs;
197     int packages_count;
198     MXFMetadataSet **metadata_sets;
199     int metadata_sets_count;
200     AVFormatContext *fc;
201     struct AVAES *aesc;
202     uint8_t *local_tags;
203     int local_tags_count;
204     uint64_t footer_partition;
205     KLVPacket current_klv_data;
206     int current_klv_index;
207     int run_in;
208     MXFPartition *current_partition;
209     int parsing_backward;
210     int64_t last_forward_tell;
211     int last_forward_partition;
212     int current_edit_unit;
213     int nb_index_tables;
214     MXFIndexTable *index_tables;
215     int edit_units_per_packet;      ///< how many edit units to read at a time (PCM, OPAtom)
216 } MXFContext;
217
218 enum MXFWrappingScheme {
219     Frame,
220     Clip,
221 };
222
223 /* NOTE: klv_offset is not set (-1) for local keys */
224 typedef int MXFMetadataReadFunc(void *arg, AVIOContext *pb, int tag, int size, UID uid, int64_t klv_offset);
225
226 typedef struct {
227     const UID key;
228     MXFMetadataReadFunc *read;
229     int ctx_size;
230     enum MXFMetadataSetType type;
231 } MXFMetadataReadTableEntry;
232
233 /* partial keys to match */
234 static const uint8_t mxf_header_partition_pack_key[]       = { 0x06,0x0e,0x2b,0x34,0x02,0x05,0x01,0x01,0x0d,0x01,0x02,0x01,0x01,0x02 };
235 static const uint8_t mxf_essence_element_key[]             = { 0x06,0x0e,0x2b,0x34,0x01,0x02,0x01,0x01,0x0d,0x01,0x03,0x01 };
236 static const uint8_t mxf_avid_essence_element_key[]        = { 0x06,0x0e,0x2b,0x34,0x01,0x02,0x01,0x01,0x0e,0x04,0x03,0x01 };
237 static const uint8_t mxf_system_item_key[]                 = { 0x06,0x0E,0x2B,0x34,0x02,0x05,0x01,0x01,0x0D,0x01,0x03,0x01,0x04 };
238 static const uint8_t mxf_klv_key[]                         = { 0x06,0x0e,0x2b,0x34 };
239 /* complete keys to match */
240 static const uint8_t mxf_crypto_source_container_ul[]      = { 0x06,0x0e,0x2b,0x34,0x01,0x01,0x01,0x09,0x06,0x01,0x01,0x02,0x02,0x00,0x00,0x00 };
241 static const uint8_t mxf_encrypted_triplet_key[]           = { 0x06,0x0e,0x2b,0x34,0x02,0x04,0x01,0x07,0x0d,0x01,0x03,0x01,0x02,0x7e,0x01,0x00 };
242 static const uint8_t mxf_encrypted_essence_container[]     = { 0x06,0x0e,0x2b,0x34,0x04,0x01,0x01,0x07,0x0d,0x01,0x03,0x01,0x02,0x0b,0x01,0x00 };
243 static const uint8_t mxf_sony_mpeg4_extradata[]            = { 0x06,0x0e,0x2b,0x34,0x04,0x01,0x01,0x01,0x0e,0x06,0x06,0x02,0x02,0x01,0x00,0x00 };
244
245 #define IS_KLV_KEY(x, y) (!memcmp(x, y, sizeof(y)))
246
247 static int64_t klv_decode_ber_length(AVIOContext *pb)
248 {
249     uint64_t size = avio_r8(pb);
250     if (size & 0x80) { /* long form */
251         int bytes_num = size & 0x7f;
252         /* SMPTE 379M 5.3.4 guarantee that bytes_num must not exceed 8 bytes */
253         if (bytes_num > 8)
254             return AVERROR_INVALIDDATA;
255         size = 0;
256         while (bytes_num--)
257             size = size << 8 | avio_r8(pb);
258     }
259     return size;
260 }
261
262 static int mxf_read_sync(AVIOContext *pb, const uint8_t *key, unsigned size)
263 {
264     int i, b;
265     for (i = 0; i < size && !pb->eof_reached; i++) {
266         b = avio_r8(pb);
267         if (b == key[0])
268             i = 0;
269         else if (b != key[i])
270             i = -1;
271     }
272     return i == size;
273 }
274
275 static int klv_read_packet(KLVPacket *klv, AVIOContext *pb)
276 {
277     if (!mxf_read_sync(pb, mxf_klv_key, 4))
278         return AVERROR_INVALIDDATA;
279     klv->offset = avio_tell(pb) - 4;
280     memcpy(klv->key, mxf_klv_key, 4);
281     avio_read(pb, klv->key + 4, 12);
282     klv->length = klv_decode_ber_length(pb);
283     return klv->length == -1 ? -1 : 0;
284 }
285
286 static int mxf_get_stream_index(AVFormatContext *s, KLVPacket *klv)
287 {
288     int i;
289
290     for (i = 0; i < s->nb_streams; i++) {
291         MXFTrack *track = s->streams[i]->priv_data;
292         /* SMPTE 379M 7.3 */
293         if (!memcmp(klv->key + sizeof(mxf_essence_element_key), track->track_number, sizeof(track->track_number)))
294             return i;
295     }
296     /* return 0 if only one stream, for OP Atom files with 0 as track number */
297     return s->nb_streams == 1 ? 0 : -1;
298 }
299
300 /* XXX: use AVBitStreamFilter */
301 static int mxf_get_d10_aes3_packet(AVIOContext *pb, AVStream *st, AVPacket *pkt, int64_t length)
302 {
303     const uint8_t *buf_ptr, *end_ptr;
304     uint8_t *data_ptr;
305     int i;
306
307     if (length > 61444) /* worst case PAL 1920 samples 8 channels */
308         return AVERROR_INVALIDDATA;
309     length = av_get_packet(pb, pkt, length);
310     if (length < 0)
311         return length;
312     data_ptr = pkt->data;
313     end_ptr = pkt->data + length;
314     buf_ptr = pkt->data + 4; /* skip SMPTE 331M header */
315     for (; buf_ptr + st->codec->channels*4 < end_ptr; ) {
316         for (i = 0; i < st->codec->channels; i++) {
317             uint32_t sample = bytestream_get_le32(&buf_ptr);
318             if (st->codec->bits_per_coded_sample == 24)
319                 bytestream_put_le24(&data_ptr, (sample >> 4) & 0xffffff);
320             else
321                 bytestream_put_le16(&data_ptr, (sample >> 12) & 0xffff);
322         }
323         buf_ptr += 32 - st->codec->channels*4; // always 8 channels stored SMPTE 331M
324     }
325     av_shrink_packet(pkt, data_ptr - pkt->data);
326     return 0;
327 }
328
329 static int mxf_decrypt_triplet(AVFormatContext *s, AVPacket *pkt, KLVPacket *klv)
330 {
331     static const uint8_t checkv[16] = {0x43, 0x48, 0x55, 0x4b, 0x43, 0x48, 0x55, 0x4b, 0x43, 0x48, 0x55, 0x4b, 0x43, 0x48, 0x55, 0x4b};
332     MXFContext *mxf = s->priv_data;
333     AVIOContext *pb = s->pb;
334     int64_t end = avio_tell(pb) + klv->length;
335     int64_t size;
336     uint64_t orig_size;
337     uint64_t plaintext_size;
338     uint8_t ivec[16];
339     uint8_t tmpbuf[16];
340     int index;
341
342     if (!mxf->aesc && s->key && s->keylen == 16) {
343         mxf->aesc = av_aes_alloc();
344         if (!mxf->aesc)
345             return AVERROR(ENOMEM);
346         av_aes_init(mxf->aesc, s->key, 128, 1);
347     }
348     // crypto context
349     avio_skip(pb, klv_decode_ber_length(pb));
350     // plaintext offset
351     klv_decode_ber_length(pb);
352     plaintext_size = avio_rb64(pb);
353     // source klv key
354     klv_decode_ber_length(pb);
355     avio_read(pb, klv->key, 16);
356     if (!IS_KLV_KEY(klv, mxf_essence_element_key))
357         return AVERROR_INVALIDDATA;
358     index = mxf_get_stream_index(s, klv);
359     if (index < 0)
360         return AVERROR_INVALIDDATA;
361     // source size
362     klv_decode_ber_length(pb);
363     orig_size = avio_rb64(pb);
364     if (orig_size < plaintext_size)
365         return AVERROR_INVALIDDATA;
366     // enc. code
367     size = klv_decode_ber_length(pb);
368     if (size < 32 || size - 32 < orig_size)
369         return AVERROR_INVALIDDATA;
370     avio_read(pb, ivec, 16);
371     avio_read(pb, tmpbuf, 16);
372     if (mxf->aesc)
373         av_aes_crypt(mxf->aesc, tmpbuf, tmpbuf, 1, ivec, 1);
374     if (memcmp(tmpbuf, checkv, 16))
375         av_log(s, AV_LOG_ERROR, "probably incorrect decryption key\n");
376     size -= 32;
377     size = av_get_packet(pb, pkt, size);
378     if (size < 0)
379         return size;
380     else if (size < plaintext_size)
381         return AVERROR_INVALIDDATA;
382     size -= plaintext_size;
383     if (mxf->aesc)
384         av_aes_crypt(mxf->aesc, &pkt->data[plaintext_size],
385                      &pkt->data[plaintext_size], size >> 4, ivec, 1);
386     av_shrink_packet(pkt, orig_size);
387     pkt->stream_index = index;
388     avio_skip(pb, end - avio_tell(pb));
389     return 0;
390 }
391
392 static int mxf_read_primer_pack(void *arg, AVIOContext *pb, int tag, int size, UID uid, int64_t klv_offset)
393 {
394     MXFContext *mxf = arg;
395     int item_num = avio_rb32(pb);
396     int item_len = avio_rb32(pb);
397
398     if (item_len != 18) {
399         av_log_ask_for_sample(pb, "unsupported primer pack item length %d\n",
400                               item_len);
401         return AVERROR_PATCHWELCOME;
402     }
403     if (item_num > UINT_MAX / item_len)
404         return AVERROR_INVALIDDATA;
405     mxf->local_tags_count = item_num;
406     mxf->local_tags = av_malloc(item_num*item_len);
407     if (!mxf->local_tags)
408         return AVERROR(ENOMEM);
409     avio_read(pb, mxf->local_tags, item_num*item_len);
410     return 0;
411 }
412
413 static int mxf_read_partition_pack(void *arg, AVIOContext *pb, int tag, int size, UID uid, int64_t klv_offset)
414 {
415     MXFContext *mxf = arg;
416     MXFPartition *partition, *tmp_part;
417     UID op;
418     uint64_t footer_partition;
419     uint32_t nb_essence_containers;
420
421     if (mxf->partitions_count+1 >= UINT_MAX / sizeof(*mxf->partitions))
422         return AVERROR(ENOMEM);
423
424     tmp_part = av_realloc(mxf->partitions, (mxf->partitions_count + 1) * sizeof(*mxf->partitions));
425     if (!tmp_part)
426         return AVERROR(ENOMEM);
427     mxf->partitions = tmp_part;
428
429     if (mxf->parsing_backward) {
430         /* insert the new partition pack in the middle
431          * this makes the entries in mxf->partitions sorted by offset */
432         memmove(&mxf->partitions[mxf->last_forward_partition+1],
433                 &mxf->partitions[mxf->last_forward_partition],
434                 (mxf->partitions_count - mxf->last_forward_partition)*sizeof(*mxf->partitions));
435         partition = mxf->current_partition = &mxf->partitions[mxf->last_forward_partition];
436     } else {
437         mxf->last_forward_partition++;
438         partition = mxf->current_partition = &mxf->partitions[mxf->partitions_count];
439     }
440
441     memset(partition, 0, sizeof(*partition));
442     mxf->partitions_count++;
443     partition->pack_length = avio_tell(pb) - klv_offset + size;
444
445     switch(uid[13]) {
446     case 2:
447         partition->type = Header;
448         break;
449     case 3:
450         partition->type = BodyPartition;
451         break;
452     case 4:
453         partition->type = Footer;
454         break;
455     default:
456         av_log(mxf->fc, AV_LOG_ERROR, "unknown partition type %i\n", uid[13]);
457         return AVERROR_INVALIDDATA;
458     }
459
460     /* consider both footers to be closed (there is only Footer and CompleteFooter) */
461     partition->closed = partition->type == Footer || !(uid[14] & 1);
462     partition->complete = uid[14] > 2;
463     avio_skip(pb, 4);
464     partition->kag_size = avio_rb32(pb);
465     partition->this_partition = avio_rb64(pb);
466     partition->previous_partition = avio_rb64(pb);
467     footer_partition = avio_rb64(pb);
468     partition->header_byte_count = avio_rb64(pb);
469     partition->index_byte_count = avio_rb64(pb);
470     partition->index_sid = avio_rb32(pb);
471     avio_skip(pb, 8);
472     partition->body_sid = avio_rb32(pb);
473     avio_read(pb, op, sizeof(UID));
474     nb_essence_containers = avio_rb32(pb);
475
476     /* some files don'thave FooterPartition set in every partition */
477     if (footer_partition) {
478         if (mxf->footer_partition && mxf->footer_partition != footer_partition) {
479             av_log(mxf->fc, AV_LOG_ERROR,
480                    "inconsistent FooterPartition value: %"PRIu64" != %"PRIu64"\n",
481                    mxf->footer_partition, footer_partition);
482         } else {
483             mxf->footer_partition = footer_partition;
484         }
485     }
486
487     av_dlog(mxf->fc,
488             "PartitionPack: ThisPartition = 0x%"PRIX64
489             ", PreviousPartition = 0x%"PRIX64", "
490             "FooterPartition = 0x%"PRIX64", IndexSID = %i, BodySID = %i\n",
491             partition->this_partition,
492             partition->previous_partition, footer_partition,
493             partition->index_sid, partition->body_sid);
494
495     /* sanity check PreviousPartition if set */
496     if (partition->previous_partition &&
497         mxf->run_in + partition->previous_partition >= klv_offset) {
498         av_log(mxf->fc, AV_LOG_ERROR,
499                "PreviousPartition points to this partition or forward\n");
500         return AVERROR_INVALIDDATA;
501     }
502
503     if      (op[12] == 1 && op[13] == 1) mxf->op = OP1a;
504     else if (op[12] == 1 && op[13] == 2) mxf->op = OP1b;
505     else if (op[12] == 1 && op[13] == 3) mxf->op = OP1c;
506     else if (op[12] == 2 && op[13] == 1) mxf->op = OP2a;
507     else if (op[12] == 2 && op[13] == 2) mxf->op = OP2b;
508     else if (op[12] == 2 && op[13] == 3) mxf->op = OP2c;
509     else if (op[12] == 3 && op[13] == 1) mxf->op = OP3a;
510     else if (op[12] == 3 && op[13] == 2) mxf->op = OP3b;
511     else if (op[12] == 3 && op[13] == 3) mxf->op = OP3c;
512     else if (op[12] == 64&& op[13] == 1) mxf->op = OPSonyOpt;
513     else if (op[12] == 0x10) {
514         /* SMPTE 390m: "There shall be exactly one essence container"
515          * The following block deals with files that violate this, namely:
516          * 2011_DCPTEST_24FPS.V.mxf - two ECs, OP1a
517          * abcdefghiv016f56415e.mxf - zero ECs, OPAtom, output by Avid AirSpeed */
518         if (nb_essence_containers != 1) {
519             MXFOP op = nb_essence_containers ? OP1a : OPAtom;
520
521             /* only nag once */
522             if (!mxf->op)
523                 av_log(mxf->fc, AV_LOG_WARNING,
524                        "\"OPAtom\" with %u ECs - assuming %s\n",
525                        nb_essence_containers,
526                        op == OP1a ? "OP1a" : "OPAtom");
527
528             mxf->op = op;
529         } else
530             mxf->op = OPAtom;
531     } else {
532         av_log(mxf->fc, AV_LOG_ERROR, "unknown operational pattern: %02xh %02xh - guessing OP1a\n", op[12], op[13]);
533         mxf->op = OP1a;
534     }
535
536     if (partition->kag_size <= 0 || partition->kag_size > (1 << 20)) {
537         av_log(mxf->fc, AV_LOG_WARNING, "invalid KAGSize %i - guessing ", partition->kag_size);
538
539         if (mxf->op == OPSonyOpt)
540             partition->kag_size = 512;
541         else
542             partition->kag_size = 1;
543
544         av_log(mxf->fc, AV_LOG_WARNING, "%i\n", partition->kag_size);
545     }
546
547     return 0;
548 }
549
550 static int mxf_add_metadata_set(MXFContext *mxf, void *metadata_set)
551 {
552     MXFMetadataSet **tmp;
553     if (mxf->metadata_sets_count+1 >= UINT_MAX / sizeof(*mxf->metadata_sets))
554         return AVERROR(ENOMEM);
555     tmp = av_realloc(mxf->metadata_sets, (mxf->metadata_sets_count + 1) * sizeof(*mxf->metadata_sets));
556     if (!tmp)
557         return AVERROR(ENOMEM);
558     mxf->metadata_sets = tmp;
559     mxf->metadata_sets[mxf->metadata_sets_count] = metadata_set;
560     mxf->metadata_sets_count++;
561     return 0;
562 }
563
564 static int mxf_read_cryptographic_context(void *arg, AVIOContext *pb, int tag, int size, UID uid, int64_t klv_offset)
565 {
566     MXFCryptoContext *cryptocontext = arg;
567     if (size != 16)
568         return AVERROR_INVALIDDATA;
569     if (IS_KLV_KEY(uid, mxf_crypto_source_container_ul))
570         avio_read(pb, cryptocontext->source_container_ul, 16);
571     return 0;
572 }
573
574 static int mxf_read_content_storage(void *arg, AVIOContext *pb, int tag, int size, UID uid, int64_t klv_offset)
575 {
576     MXFContext *mxf = arg;
577     switch (tag) {
578     case 0x1901:
579         mxf->packages_count = avio_rb32(pb);
580         if (mxf->packages_count >= UINT_MAX / sizeof(UID))
581             return AVERROR_INVALIDDATA;
582         mxf->packages_refs = av_malloc(mxf->packages_count * sizeof(UID));
583         if (!mxf->packages_refs)
584             return AVERROR(ENOMEM);
585         avio_skip(pb, 4); /* useless size of objects, always 16 according to specs */
586         avio_read(pb, (uint8_t *)mxf->packages_refs, mxf->packages_count * sizeof(UID));
587         break;
588     }
589     return 0;
590 }
591
592 static int mxf_read_source_clip(void *arg, AVIOContext *pb, int tag, int size, UID uid, int64_t klv_offset)
593 {
594     MXFStructuralComponent *source_clip = arg;
595     switch(tag) {
596     case 0x0202:
597         source_clip->duration = avio_rb64(pb);
598         break;
599     case 0x1201:
600         source_clip->start_position = avio_rb64(pb);
601         break;
602     case 0x1101:
603         /* UMID, only get last 16 bytes */
604         avio_skip(pb, 16);
605         avio_read(pb, source_clip->source_package_uid, 16);
606         break;
607     case 0x1102:
608         source_clip->source_track_id = avio_rb32(pb);
609         break;
610     }
611     return 0;
612 }
613
614 static int mxf_read_material_package(void *arg, AVIOContext *pb, int tag, int size, UID uid, int64_t klv_offset)
615 {
616     MXFPackage *package = arg;
617     switch(tag) {
618     case 0x4403:
619         package->tracks_count = avio_rb32(pb);
620         if (package->tracks_count >= UINT_MAX / sizeof(UID))
621             return AVERROR_INVALIDDATA;
622         package->tracks_refs = av_malloc(package->tracks_count * sizeof(UID));
623         if (!package->tracks_refs)
624             return AVERROR(ENOMEM);
625         avio_skip(pb, 4); /* useless size of objects, always 16 according to specs */
626         avio_read(pb, (uint8_t *)package->tracks_refs, package->tracks_count * sizeof(UID));
627         break;
628     }
629     return 0;
630 }
631
632 static int mxf_read_track(void *arg, AVIOContext *pb, int tag, int size, UID uid, int64_t klv_offset)
633 {
634     MXFTrack *track = arg;
635     switch(tag) {
636     case 0x4801:
637         track->track_id = avio_rb32(pb);
638         break;
639     case 0x4804:
640         avio_read(pb, track->track_number, 4);
641         break;
642     case 0x4B01:
643         track->edit_rate.num = avio_rb32(pb);
644         track->edit_rate.den = avio_rb32(pb);
645         break;
646     case 0x4803:
647         avio_read(pb, track->sequence_ref, 16);
648         break;
649     }
650     return 0;
651 }
652
653 static int mxf_read_sequence(void *arg, AVIOContext *pb, int tag, int size, UID uid, int64_t klv_offset)
654 {
655     MXFSequence *sequence = arg;
656     switch(tag) {
657     case 0x0202:
658         sequence->duration = avio_rb64(pb);
659         break;
660     case 0x0201:
661         avio_read(pb, sequence->data_definition_ul, 16);
662         break;
663     case 0x1001:
664         sequence->structural_components_count = avio_rb32(pb);
665         if (sequence->structural_components_count >= UINT_MAX / sizeof(UID))
666             return AVERROR_INVALIDDATA;
667         sequence->structural_components_refs = av_malloc(sequence->structural_components_count * sizeof(UID));
668         if (!sequence->structural_components_refs)
669             return AVERROR(ENOMEM);
670         avio_skip(pb, 4); /* useless size of objects, always 16 according to specs */
671         avio_read(pb, (uint8_t *)sequence->structural_components_refs, sequence->structural_components_count * sizeof(UID));
672         break;
673     }
674     return 0;
675 }
676
677 static int mxf_read_source_package(void *arg, AVIOContext *pb, int tag, int size, UID uid, int64_t klv_offset)
678 {
679     MXFPackage *package = arg;
680     switch(tag) {
681     case 0x4403:
682         package->tracks_count = avio_rb32(pb);
683         if (package->tracks_count >= UINT_MAX / sizeof(UID))
684             return AVERROR_INVALIDDATA;
685         package->tracks_refs = av_malloc(package->tracks_count * sizeof(UID));
686         if (!package->tracks_refs)
687             return AVERROR(ENOMEM);
688         avio_skip(pb, 4); /* useless size of objects, always 16 according to specs */
689         avio_read(pb, (uint8_t *)package->tracks_refs, package->tracks_count * sizeof(UID));
690         break;
691     case 0x4401:
692         /* UMID, only get last 16 bytes */
693         avio_skip(pb, 16);
694         avio_read(pb, package->package_uid, 16);
695         break;
696     case 0x4701:
697         avio_read(pb, package->descriptor_ref, 16);
698         break;
699     }
700     return 0;
701 }
702
703 static int mxf_read_index_entry_array(AVIOContext *pb, MXFIndexTableSegment *segment)
704 {
705     int i, length;
706
707     segment->nb_index_entries = avio_rb32(pb);
708     if (!segment->nb_index_entries)
709         return 0;
710     else if (segment->nb_index_entries < 0 ||
711              segment->nb_index_entries >
712              (INT_MAX / sizeof(*segment->stream_offset_entries)))
713         return AVERROR(ENOMEM);
714
715     length = avio_rb32(pb);
716
717     segment->temporal_offset_entries = av_mallocz(segment->nb_index_entries *
718                                  sizeof(*segment->temporal_offset_entries));
719     segment->flag_entries            = av_mallocz(segment->nb_index_entries *
720                                  sizeof(*segment->flag_entries));
721     segment->stream_offset_entries   = av_mallocz(segment->nb_index_entries *
722                                  sizeof(*segment->stream_offset_entries));
723
724     if (!segment->flag_entries || !segment->stream_offset_entries ||
725         !segment->temporal_offset_entries) {
726         av_freep(&segment->flag_entries);
727         av_freep(&segment->stream_offset_entries);
728         av_freep(&segment->temporal_offset_entries);
729         return AVERROR(ENOMEM);
730     }
731
732     for (i = 0; i < segment->nb_index_entries; i++) {
733         segment->temporal_offset_entries[i] = avio_r8(pb);
734         avio_r8(pb);                                        /* KeyFrameOffset */
735         segment->flag_entries[i] = avio_r8(pb);
736         segment->stream_offset_entries[i] = avio_rb64(pb);
737         avio_skip(pb, length - 11);
738     }
739     return 0;
740 }
741
742 static int mxf_read_index_table_segment(void *arg, AVIOContext *pb, int tag, int size, UID uid, int64_t klv_offset)
743 {
744     MXFIndexTableSegment *segment = arg;
745     switch(tag) {
746     case 0x3F05:
747         segment->edit_unit_byte_count = avio_rb32(pb);
748         av_dlog(NULL, "EditUnitByteCount %d\n", segment->edit_unit_byte_count);
749         break;
750     case 0x3F06:
751         segment->index_sid = avio_rb32(pb);
752         av_dlog(NULL, "IndexSID %d\n", segment->index_sid);
753         break;
754     case 0x3F07:
755         segment->body_sid = avio_rb32(pb);
756         av_dlog(NULL, "BodySID %d\n", segment->body_sid);
757         break;
758     case 0x3F0A:
759         av_dlog(NULL, "IndexEntryArray found\n");
760         return mxf_read_index_entry_array(pb, segment);
761     case 0x3F0B:
762         segment->index_edit_rate.num = avio_rb32(pb);
763         segment->index_edit_rate.den = avio_rb32(pb);
764         av_dlog(NULL, "IndexEditRate %d/%d\n", segment->index_edit_rate.num,
765                 segment->index_edit_rate.den);
766         break;
767     case 0x3F0C:
768         segment->index_start_position = avio_rb64(pb);
769         av_dlog(NULL, "IndexStartPosition %"PRId64"\n", segment->index_start_position);
770         break;
771     case 0x3F0D:
772         segment->index_duration = avio_rb64(pb);
773         av_dlog(NULL, "IndexDuration %"PRId64"\n", segment->index_duration);
774         break;
775     }
776     return 0;
777 }
778
779 static void mxf_read_pixel_layout(AVIOContext *pb, MXFDescriptor *descriptor)
780 {
781     int code, value, ofs = 0;
782     char layout[16] = {0};
783
784     do {
785         code = avio_r8(pb);
786         value = avio_r8(pb);
787         av_dlog(NULL, "pixel layout: code %#x\n", code);
788
789         if (ofs < 16) {
790             layout[ofs++] = code;
791             layout[ofs++] = value;
792         }
793     } while (code != 0); /* SMPTE 377M E.2.46 */
794
795     ff_mxf_decode_pixel_layout(layout, &descriptor->pix_fmt);
796 }
797
798 static int mxf_read_generic_descriptor(void *arg, AVIOContext *pb, int tag, int size, UID uid, int64_t klv_offset)
799 {
800     MXFDescriptor *descriptor = arg;
801     descriptor->pix_fmt = AV_PIX_FMT_NONE;
802     switch(tag) {
803     case 0x3F01:
804         descriptor->sub_descriptors_count = avio_rb32(pb);
805         if (descriptor->sub_descriptors_count >= UINT_MAX / sizeof(UID))
806             return AVERROR_INVALIDDATA;
807         descriptor->sub_descriptors_refs = av_malloc(descriptor->sub_descriptors_count * sizeof(UID));
808         if (!descriptor->sub_descriptors_refs)
809             return AVERROR(ENOMEM);
810         avio_skip(pb, 4); /* useless size of objects, always 16 according to specs */
811         avio_read(pb, (uint8_t *)descriptor->sub_descriptors_refs, descriptor->sub_descriptors_count * sizeof(UID));
812         break;
813     case 0x3004:
814         avio_read(pb, descriptor->essence_container_ul, 16);
815         break;
816     case 0x3006:
817         descriptor->linked_track_id = avio_rb32(pb);
818         break;
819     case 0x3201: /* PictureEssenceCoding */
820         avio_read(pb, descriptor->essence_codec_ul, 16);
821         break;
822     case 0x3203:
823         descriptor->width = avio_rb32(pb);
824         break;
825     case 0x3202:
826         descriptor->height = avio_rb32(pb);
827         break;
828     case 0x320C:
829         descriptor->frame_layout = avio_r8(pb);
830         break;
831     case 0x320E:
832         descriptor->aspect_ratio.num = avio_rb32(pb);
833         descriptor->aspect_ratio.den = avio_rb32(pb);
834         break;
835     case 0x3301:
836         descriptor->component_depth = avio_rb32(pb);
837         break;
838     case 0x3302:
839         descriptor->horiz_subsampling = avio_rb32(pb);
840         break;
841     case 0x3308:
842         descriptor->vert_subsampling = avio_rb32(pb);
843         break;
844     case 0x3D03:
845         descriptor->sample_rate.num = avio_rb32(pb);
846         descriptor->sample_rate.den = avio_rb32(pb);
847         break;
848     case 0x3D06: /* SoundEssenceCompression */
849         avio_read(pb, descriptor->essence_codec_ul, 16);
850         break;
851     case 0x3D07:
852         descriptor->channels = avio_rb32(pb);
853         break;
854     case 0x3D01:
855         descriptor->bits_per_sample = avio_rb32(pb);
856         break;
857     case 0x3401:
858         mxf_read_pixel_layout(pb, descriptor);
859         break;
860     default:
861         /* Private uid used by SONY C0023S01.mxf */
862         if (IS_KLV_KEY(uid, mxf_sony_mpeg4_extradata)) {
863             descriptor->extradata = av_malloc(size + FF_INPUT_BUFFER_PADDING_SIZE);
864             if (!descriptor->extradata)
865                 return AVERROR(ENOMEM);
866             descriptor->extradata_size = size;
867             avio_read(pb, descriptor->extradata, size);
868         }
869         break;
870     }
871     return 0;
872 }
873
874 /*
875  * Match an uid independently of the version byte and up to len common bytes
876  * Returns: boolean
877  */
878 static int mxf_match_uid(const UID key, const UID uid, int len)
879 {
880     int i;
881     for (i = 0; i < len; i++) {
882         if (i != 7 && key[i] != uid[i])
883             return 0;
884     }
885     return 1;
886 }
887
888 static const MXFCodecUL *mxf_get_codec_ul(const MXFCodecUL *uls, UID *uid)
889 {
890     while (uls->uid[0]) {
891         if(mxf_match_uid(uls->uid, *uid, uls->matching_len))
892             break;
893         uls++;
894     }
895     return uls;
896 }
897
898 static void *mxf_resolve_strong_ref(MXFContext *mxf, UID *strong_ref, enum MXFMetadataSetType type)
899 {
900     int i;
901
902     if (!strong_ref)
903         return NULL;
904     for (i = 0; i < mxf->metadata_sets_count; i++) {
905         if (!memcmp(*strong_ref, mxf->metadata_sets[i]->uid, 16) &&
906             (type == AnyType || mxf->metadata_sets[i]->type == type)) {
907             return mxf->metadata_sets[i];
908         }
909     }
910     return NULL;
911 }
912
913 static const MXFCodecUL mxf_picture_essence_container_uls[] = {
914     // video essence container uls
915     { { 0x06,0x0E,0x2B,0x34,0x04,0x01,0x01,0x02,0x0D,0x01,0x03,0x01,0x02,0x04,0x60,0x01 }, 14, AV_CODEC_ID_MPEG2VIDEO }, /* MPEG-ES Frame wrapped */
916     { { 0x06,0x0E,0x2B,0x34,0x04,0x01,0x01,0x01,0x0D,0x01,0x03,0x01,0x02,0x02,0x41,0x01 }, 14,    AV_CODEC_ID_DVVIDEO }, /* DV 625 25mbps */
917     { { 0x06,0x0E,0x2B,0x34,0x04,0x01,0x01,0x01,0x0D,0x01,0x03,0x01,0x02,0x05,0x00,0x00 }, 14,   AV_CODEC_ID_RAWVIDEO }, /* Uncompressed Picture */
918     { { 0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00 },  0,      AV_CODEC_ID_NONE },
919 };
920
921 /* EC ULs for intra-only formats */
922 static const MXFCodecUL mxf_intra_only_essence_container_uls[] = {
923     { { 0x06,0x0E,0x2B,0x34,0x04,0x01,0x01,0x01,0x0D,0x01,0x03,0x01,0x02,0x01,0x00,0x00 }, 14, AV_CODEC_ID_MPEG2VIDEO }, /* MXF-GC SMPTE D-10 Mappings */
924     { { 0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00 },  0,       AV_CODEC_ID_NONE },
925 };
926
927 /* intra-only PictureEssenceCoding ULs, where no corresponding EC UL exists */
928 static const MXFCodecUL mxf_intra_only_picture_essence_coding_uls[] = {
929     { { 0x06,0x0E,0x2B,0x34,0x04,0x01,0x01,0x0A,0x04,0x01,0x02,0x02,0x01,0x32,0x00,0x00 }, 14,       AV_CODEC_ID_H264 }, /* H.264/MPEG-4 AVC Intra Profiles */
930     { { 0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00 },  0,       AV_CODEC_ID_NONE },
931 };
932
933 static const MXFCodecUL mxf_sound_essence_container_uls[] = {
934     // sound essence container uls
935     { { 0x06,0x0E,0x2B,0x34,0x04,0x01,0x01,0x01,0x0D,0x01,0x03,0x01,0x02,0x06,0x01,0x00 }, 14, AV_CODEC_ID_PCM_S16LE }, /* BWF Frame wrapped */
936     { { 0x06,0x0E,0x2B,0x34,0x04,0x01,0x01,0x02,0x0D,0x01,0x03,0x01,0x02,0x04,0x40,0x01 }, 14,       AV_CODEC_ID_MP2 }, /* MPEG-ES Frame wrapped, 0x40 ??? stream id */
937     { { 0x06,0x0E,0x2B,0x34,0x04,0x01,0x01,0x01,0x0D,0x01,0x03,0x01,0x02,0x01,0x01,0x01 }, 14, AV_CODEC_ID_PCM_S16LE }, /* D-10 Mapping 50Mbps PAL Extended Template */
938     { { 0x06,0x0E,0x2B,0x34,0x01,0x01,0x01,0xFF,0x4B,0x46,0x41,0x41,0x00,0x0D,0x4D,0x4F }, 14, AV_CODEC_ID_PCM_S16LE }, /* 0001GL00.MXF.A1.mxf_opatom.mxf */
939     { { 0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00 },  0,      AV_CODEC_ID_NONE },
940 };
941
942 static int mxf_get_sorted_table_segments(MXFContext *mxf, int *nb_sorted_segments, MXFIndexTableSegment ***sorted_segments)
943 {
944     int i, j, nb_segments = 0;
945     MXFIndexTableSegment **unsorted_segments;
946     int last_body_sid = -1, last_index_sid = -1, last_index_start = -1;
947
948     /* count number of segments, allocate arrays and copy unsorted segments */
949     for (i = 0; i < mxf->metadata_sets_count; i++)
950         if (mxf->metadata_sets[i]->type == IndexTableSegment)
951             nb_segments++;
952
953     if (!nb_segments)
954         return AVERROR_INVALIDDATA;
955
956     *sorted_segments  = av_mallocz(nb_segments * sizeof(**sorted_segments));
957     unsorted_segments = av_mallocz(nb_segments * sizeof(*unsorted_segments));
958     if (!*sorted_segments || !unsorted_segments) {
959         av_freep(sorted_segments);
960         av_free(unsorted_segments);
961         return AVERROR(ENOMEM);
962     }
963
964     for (i = j = 0; i < mxf->metadata_sets_count; i++)
965         if (mxf->metadata_sets[i]->type == IndexTableSegment)
966             unsorted_segments[j++] = (MXFIndexTableSegment*)mxf->metadata_sets[i];
967
968     *nb_sorted_segments = 0;
969
970     /* sort segments by {BodySID, IndexSID, IndexStartPosition}, remove duplicates while we're at it */
971     for (i = 0; i < nb_segments; i++) {
972         int best = -1, best_body_sid = -1, best_index_sid = -1, best_index_start = -1;
973
974         for (j = 0; j < nb_segments; j++) {
975             MXFIndexTableSegment *s = unsorted_segments[j];
976
977             /* Require larger BosySID, IndexSID or IndexStartPosition then the previous entry. This removes duplicates.
978              * We want the smallest values for the keys than what we currently have, unless this is the first such entry this time around.
979              */
980             if ((i == 0     || s->body_sid > last_body_sid || s->index_sid > last_index_sid || s->index_start_position > last_index_start) &&
981                 (best == -1 || s->body_sid < best_body_sid || s->index_sid < best_index_sid || s->index_start_position < best_index_start)) {
982                 best             = j;
983                 best_body_sid    = s->body_sid;
984                 best_index_sid   = s->index_sid;
985                 best_index_start = s->index_start_position;
986             }
987         }
988
989         /* no suitable entry found -> we're done */
990         if (best == -1)
991             break;
992
993         (*sorted_segments)[(*nb_sorted_segments)++] = unsorted_segments[best];
994         last_body_sid    = best_body_sid;
995         last_index_sid   = best_index_sid;
996         last_index_start = best_index_start;
997     }
998
999     av_free(unsorted_segments);
1000
1001     return 0;
1002 }
1003
1004 /**
1005  * Computes the absolute file offset of the given essence container offset
1006  */
1007 static int mxf_absolute_bodysid_offset(MXFContext *mxf, int body_sid, int64_t offset, int64_t *offset_out)
1008 {
1009     int x;
1010     int64_t offset_in = offset;     /* for logging */
1011
1012     for (x = 0; x < mxf->partitions_count; x++) {
1013         MXFPartition *p = &mxf->partitions[x];
1014
1015         if (p->body_sid != body_sid)
1016             continue;
1017
1018         if (offset < p->essence_length || !p->essence_length) {
1019             *offset_out = p->essence_offset + offset;
1020             return 0;
1021         }
1022
1023         offset -= p->essence_length;
1024     }
1025
1026     av_log(mxf->fc, AV_LOG_ERROR,
1027            "failed to find absolute offset of %"PRIX64" in BodySID %i - partial file?\n",
1028            offset_in, body_sid);
1029
1030     return AVERROR_INVALIDDATA;
1031 }
1032
1033 /**
1034  * Returns the end position of the essence container with given BodySID, or zero if unknown
1035  */
1036 static int64_t mxf_essence_container_end(MXFContext *mxf, int body_sid)
1037 {
1038     int x;
1039     int64_t ret = 0;
1040
1041     for (x = 0; x < mxf->partitions_count; x++) {
1042         MXFPartition *p = &mxf->partitions[x];
1043
1044         if (p->body_sid != body_sid)
1045             continue;
1046
1047         if (!p->essence_length)
1048             return 0;
1049
1050         ret = p->essence_offset + p->essence_length;
1051     }
1052
1053     return ret;
1054 }
1055
1056 /* EditUnit -> absolute offset */
1057 static int mxf_edit_unit_absolute_offset(MXFContext *mxf, MXFIndexTable *index_table, int64_t edit_unit, int64_t *edit_unit_out, int64_t *offset_out, int nag)
1058 {
1059     int i;
1060     int64_t offset_temp = 0;
1061
1062     for (i = 0; i < index_table->nb_segments; i++) {
1063         MXFIndexTableSegment *s = index_table->segments[i];
1064
1065         edit_unit = FFMAX(edit_unit, s->index_start_position);  /* clamp if trying to seek before start */
1066
1067         if (edit_unit < s->index_start_position + s->index_duration) {
1068             int64_t index = edit_unit - s->index_start_position;
1069
1070             if (s->edit_unit_byte_count)
1071                 offset_temp += s->edit_unit_byte_count * index;
1072             else if (s->nb_index_entries) {
1073                 if (s->nb_index_entries == 2 * s->index_duration + 1)
1074                     index *= 2;     /* Avid index */
1075
1076                 if (index < 0 || index >= s->nb_index_entries) {
1077                     av_log(mxf->fc, AV_LOG_ERROR, "IndexSID %i segment at %"PRId64" IndexEntryArray too small\n",
1078                            index_table->index_sid, s->index_start_position);
1079                     return AVERROR_INVALIDDATA;
1080                 }
1081
1082                 offset_temp = s->stream_offset_entries[index];
1083             } else {
1084                 av_log(mxf->fc, AV_LOG_ERROR, "IndexSID %i segment at %"PRId64" missing EditUnitByteCount and IndexEntryArray\n",
1085                        index_table->index_sid, s->index_start_position);
1086                 return AVERROR_INVALIDDATA;
1087             }
1088
1089             if (edit_unit_out)
1090                 *edit_unit_out = edit_unit;
1091
1092             return mxf_absolute_bodysid_offset(mxf, index_table->body_sid, offset_temp, offset_out);
1093         } else {
1094             /* EditUnitByteCount == 0 for VBR indexes, which is fine since they use explicit StreamOffsets */
1095             offset_temp += s->edit_unit_byte_count * s->index_duration;
1096         }
1097     }
1098
1099     if (nag)
1100         av_log(mxf->fc, AV_LOG_ERROR, "failed to map EditUnit %"PRId64" in IndexSID %i to an offset\n", edit_unit, index_table->index_sid);
1101
1102     return AVERROR_INVALIDDATA;
1103 }
1104
1105 static int mxf_compute_ptses_fake_index(MXFContext *mxf, MXFIndexTable *index_table)
1106 {
1107     int i, j, x;
1108     int8_t max_temporal_offset = -128;
1109
1110     /* first compute how many entries we have */
1111     for (i = 0; i < index_table->nb_segments; i++) {
1112         MXFIndexTableSegment *s = index_table->segments[i];
1113
1114         if (!s->nb_index_entries) {
1115             index_table->nb_ptses = 0;
1116             return 0;                               /* no TemporalOffsets */
1117         }
1118
1119         index_table->nb_ptses += s->index_duration;
1120     }
1121
1122     /* paranoid check */
1123     if (index_table->nb_ptses <= 0)
1124         return 0;
1125
1126     if (index_table->nb_ptses > INT_MAX / sizeof(AVIndexEntry))
1127         return AVERROR(ENOMEM);
1128
1129     index_table->ptses      = av_mallocz(index_table->nb_ptses *
1130                                          sizeof(int64_t));
1131     index_table->fake_index = av_mallocz(index_table->nb_ptses *
1132                                          sizeof(AVIndexEntry));
1133     if (!index_table->ptses || !index_table->fake_index) {
1134         av_freep(&index_table->ptses);
1135         return AVERROR(ENOMEM);
1136     }
1137
1138     /* we may have a few bad TemporalOffsets
1139      * make sure the corresponding PTSes don't have the bogus value 0 */
1140     for (x = 0; x < index_table->nb_ptses; x++)
1141         index_table->ptses[x] = AV_NOPTS_VALUE;
1142
1143     /**
1144      * We have this:
1145      *
1146      * x  TemporalOffset
1147      * 0:  0
1148      * 1:  1
1149      * 2:  1
1150      * 3: -2
1151      * 4:  1
1152      * 5:  1
1153      * 6: -2
1154      *
1155      * We want to transform it into this:
1156      *
1157      * x  DTS PTS
1158      * 0: -1   0
1159      * 1:  0   3
1160      * 2:  1   1
1161      * 3:  2   2
1162      * 4:  3   6
1163      * 5:  4   4
1164      * 6:  5   5
1165      *
1166      * We do this by bucket sorting x by x+TemporalOffset[x] into mxf->ptses,
1167      * then settings mxf->first_dts = -max(TemporalOffset[x]).
1168      * The latter makes DTS <= PTS.
1169      */
1170     for (i = x = 0; i < index_table->nb_segments; i++) {
1171         MXFIndexTableSegment *s = index_table->segments[i];
1172         int index_delta = 1;
1173         int n = s->nb_index_entries;
1174
1175         if (s->nb_index_entries == 2 * s->index_duration + 1) {
1176             index_delta = 2;    /* Avid index */
1177             /* ignore the last entry - it's the size of the essence container */
1178             n--;
1179         }
1180
1181         for (j = 0; j < n; j += index_delta, x++) {
1182             int offset = s->temporal_offset_entries[j] / index_delta;
1183             int index  = x + offset;
1184
1185             if (x >= index_table->nb_ptses) {
1186                 av_log(mxf->fc, AV_LOG_ERROR,
1187                        "x >= nb_ptses - IndexEntryCount %i < IndexDuration %"PRId64"?\n",
1188                        s->nb_index_entries, s->index_duration);
1189                 break;
1190             }
1191
1192             index_table->fake_index[x].timestamp = x;
1193             index_table->fake_index[x].flags = !(s->flag_entries[j] & 0x30) ? AVINDEX_KEYFRAME : 0;
1194
1195             if (index < 0 || index >= index_table->nb_ptses) {
1196                 av_log(mxf->fc, AV_LOG_ERROR,
1197                        "index entry %i + TemporalOffset %i = %i, which is out of bounds\n",
1198                        x, offset, index);
1199                 continue;
1200             }
1201
1202             index_table->ptses[index] = x;
1203             max_temporal_offset = FFMAX(max_temporal_offset, offset);
1204         }
1205     }
1206
1207     index_table->first_dts = -max_temporal_offset;
1208
1209     return 0;
1210 }
1211
1212 /**
1213  * Sorts and collects index table segments into index tables.
1214  * Also computes PTSes if possible.
1215  */
1216 static int mxf_compute_index_tables(MXFContext *mxf)
1217 {
1218     int i, j, k, ret, nb_sorted_segments;
1219     MXFIndexTableSegment **sorted_segments = NULL;
1220
1221     if ((ret = mxf_get_sorted_table_segments(mxf, &nb_sorted_segments, &sorted_segments)) ||
1222         nb_sorted_segments <= 0) {
1223         av_log(mxf->fc, AV_LOG_WARNING, "broken or empty index\n");
1224         return 0;
1225     }
1226
1227     /* sanity check and count unique BodySIDs/IndexSIDs */
1228     for (i = 0; i < nb_sorted_segments; i++) {
1229         if (i == 0 || sorted_segments[i-1]->index_sid != sorted_segments[i]->index_sid)
1230             mxf->nb_index_tables++;
1231         else if (sorted_segments[i-1]->body_sid != sorted_segments[i]->body_sid) {
1232             av_log(mxf->fc, AV_LOG_ERROR, "found inconsistent BodySID\n");
1233             ret = AVERROR_INVALIDDATA;
1234             goto finish_decoding_index;
1235         }
1236     }
1237
1238     if (mxf->nb_index_tables > INT_MAX / sizeof(MXFIndexTable) ||
1239         !(mxf->index_tables = av_mallocz(mxf->nb_index_tables *
1240                                          sizeof(MXFIndexTable)))) {
1241         av_log(mxf->fc, AV_LOG_ERROR, "failed to allocate index tables\n");
1242         ret = AVERROR(ENOMEM);
1243         goto finish_decoding_index;
1244     }
1245
1246     /* distribute sorted segments to index tables */
1247     for (i = j = 0; i < nb_sorted_segments; i++) {
1248         if (i != 0 && sorted_segments[i-1]->index_sid != sorted_segments[i]->index_sid) {
1249             /* next IndexSID */
1250             j++;
1251         }
1252
1253         mxf->index_tables[j].nb_segments++;
1254     }
1255
1256     for (i = j = 0; j < mxf->nb_index_tables; i += mxf->index_tables[j++].nb_segments) {
1257         MXFIndexTable *t = &mxf->index_tables[j];
1258
1259         if (t->nb_segments >
1260             (INT_MAX / sizeof(MXFIndexTableSegment *)) ||
1261             !(t->segments = av_mallocz(t->nb_segments *
1262                                        sizeof(MXFIndexTableSegment*)))) {
1263             av_log(mxf->fc, AV_LOG_ERROR, "failed to allocate IndexTableSegment"
1264                    " pointer array\n");
1265             ret = AVERROR(ENOMEM);
1266             goto finish_decoding_index;
1267         }
1268
1269         if (sorted_segments[i]->index_start_position)
1270             av_log(mxf->fc, AV_LOG_WARNING, "IndexSID %i starts at EditUnit %"PRId64" - seeking may not work as expected\n",
1271                    sorted_segments[i]->index_sid, sorted_segments[i]->index_start_position);
1272
1273         memcpy(t->segments, &sorted_segments[i], t->nb_segments * sizeof(MXFIndexTableSegment*));
1274         t->index_sid = sorted_segments[i]->index_sid;
1275         t->body_sid = sorted_segments[i]->body_sid;
1276
1277         if ((ret = mxf_compute_ptses_fake_index(mxf, t)) < 0)
1278             goto finish_decoding_index;
1279
1280         /* fix zero IndexDurations */
1281         for (k = 0; k < t->nb_segments; k++) {
1282             if (t->segments[k]->index_duration)
1283                 continue;
1284
1285             if (t->nb_segments > 1)
1286                 av_log(mxf->fc, AV_LOG_WARNING, "IndexSID %i segment %i has zero IndexDuration and there's more than one segment\n",
1287                        t->index_sid, k);
1288
1289             if (mxf->fc->nb_streams <= 0) {
1290                 av_log(mxf->fc, AV_LOG_WARNING, "no streams?\n");
1291                 break;
1292             }
1293
1294             /* assume the first stream's duration is reasonable
1295              * leave index_duration = 0 on further segments in case we have any (unlikely)
1296              */
1297             t->segments[k]->index_duration = mxf->fc->streams[0]->duration;
1298             break;
1299         }
1300     }
1301
1302     ret = 0;
1303 finish_decoding_index:
1304     av_free(sorted_segments);
1305     return ret;
1306 }
1307
1308 static int mxf_is_intra_only(MXFDescriptor *d)
1309 {
1310     return mxf_get_codec_ul(mxf_intra_only_essence_container_uls,
1311                             &d->essence_container_ul)->id != AV_CODEC_ID_NONE ||
1312            mxf_get_codec_ul(mxf_intra_only_picture_essence_coding_uls,
1313                             &d->essence_codec_ul)->id     != AV_CODEC_ID_NONE;
1314 }
1315
1316 static int mxf_parse_structural_metadata(MXFContext *mxf)
1317 {
1318     MXFPackage *material_package = NULL;
1319     MXFPackage *temp_package = NULL;
1320     int i, j, k, ret;
1321
1322     av_dlog(mxf->fc, "metadata sets count %d\n", mxf->metadata_sets_count);
1323     /* TODO: handle multiple material packages (OP3x) */
1324     for (i = 0; i < mxf->packages_count; i++) {
1325         material_package = mxf_resolve_strong_ref(mxf, &mxf->packages_refs[i], MaterialPackage);
1326         if (material_package) break;
1327     }
1328     if (!material_package) {
1329         av_log(mxf->fc, AV_LOG_ERROR, "no material package found\n");
1330         return AVERROR_INVALIDDATA;
1331     }
1332
1333     for (i = 0; i < material_package->tracks_count; i++) {
1334         MXFPackage *source_package = NULL;
1335         MXFTrack *material_track = NULL;
1336         MXFTrack *source_track = NULL;
1337         MXFTrack *temp_track = NULL;
1338         MXFDescriptor *descriptor = NULL;
1339         MXFStructuralComponent *component = NULL;
1340         UID *essence_container_ul = NULL;
1341         const MXFCodecUL *codec_ul = NULL;
1342         const MXFCodecUL *container_ul = NULL;
1343         const MXFCodecUL *pix_fmt_ul = NULL;
1344         AVStream *st;
1345
1346         if (!(material_track = mxf_resolve_strong_ref(mxf, &material_package->tracks_refs[i], Track))) {
1347             av_log(mxf->fc, AV_LOG_ERROR, "could not resolve material track strong ref\n");
1348             continue;
1349         }
1350
1351         if (!(material_track->sequence = mxf_resolve_strong_ref(mxf, &material_track->sequence_ref, Sequence))) {
1352             av_log(mxf->fc, AV_LOG_ERROR, "could not resolve material track sequence strong ref\n");
1353             continue;
1354         }
1355
1356         /* TODO: handle multiple source clips */
1357         for (j = 0; j < material_track->sequence->structural_components_count; j++) {
1358             /* TODO: handle timecode component */
1359             component = mxf_resolve_strong_ref(mxf, &material_track->sequence->structural_components_refs[j], SourceClip);
1360             if (!component)
1361                 continue;
1362
1363             for (k = 0; k < mxf->packages_count; k++) {
1364                 temp_package = mxf_resolve_strong_ref(mxf, &mxf->packages_refs[k], SourcePackage);
1365                 if (!temp_package)
1366                     continue;
1367                 if (!memcmp(temp_package->package_uid, component->source_package_uid, 16)) {
1368                     source_package = temp_package;
1369                     break;
1370                 }
1371             }
1372             if (!source_package) {
1373                 av_dlog(mxf->fc, "material track %d: no corresponding source package found\n", material_track->track_id);
1374                 break;
1375             }
1376             for (k = 0; k < source_package->tracks_count; k++) {
1377                 if (!(temp_track = mxf_resolve_strong_ref(mxf, &source_package->tracks_refs[k], Track))) {
1378                     av_log(mxf->fc, AV_LOG_ERROR, "could not resolve source track strong ref\n");
1379                     ret = AVERROR_INVALIDDATA;
1380                     goto fail_and_free;
1381                 }
1382                 if (temp_track->track_id == component->source_track_id) {
1383                     source_track = temp_track;
1384                     break;
1385                 }
1386             }
1387             if (!source_track) {
1388                 av_log(mxf->fc, AV_LOG_ERROR, "material track %d: no corresponding source track found\n", material_track->track_id);
1389                 break;
1390             }
1391         }
1392         if (!source_track || !component)
1393             continue;
1394
1395         if (!(source_track->sequence = mxf_resolve_strong_ref(mxf, &source_track->sequence_ref, Sequence))) {
1396             av_log(mxf->fc, AV_LOG_ERROR, "could not resolve source track sequence strong ref\n");
1397             ret = AVERROR_INVALIDDATA;
1398             goto fail_and_free;
1399         }
1400
1401         /* 0001GL00.MXF.A1.mxf_opatom.mxf has the same SourcePackageID as 0001GL.MXF.V1.mxf_opatom.mxf
1402          * This would result in both files appearing to have two streams. Work around this by sanity checking DataDefinition */
1403         if (memcmp(material_track->sequence->data_definition_ul, source_track->sequence->data_definition_ul, 16)) {
1404             av_log(mxf->fc, AV_LOG_ERROR, "material track %d: DataDefinition mismatch\n", material_track->track_id);
1405             continue;
1406         }
1407
1408         st = avformat_new_stream(mxf->fc, NULL);
1409         if (!st) {
1410             av_log(mxf->fc, AV_LOG_ERROR, "could not allocate stream\n");
1411             ret = AVERROR(ENOMEM);
1412             goto fail_and_free;
1413         }
1414         st->id = source_track->track_id;
1415         st->priv_data = source_track;
1416         st->duration = component->duration;
1417         if (st->duration == -1)
1418             st->duration = AV_NOPTS_VALUE;
1419         st->start_time = component->start_position;
1420         avpriv_set_pts_info(st, 64, material_track->edit_rate.den, material_track->edit_rate.num);
1421
1422         PRINT_KEY(mxf->fc, "data definition   ul", source_track->sequence->data_definition_ul);
1423         codec_ul = mxf_get_codec_ul(ff_mxf_data_definition_uls, &source_track->sequence->data_definition_ul);
1424         st->codec->codec_type = codec_ul->id;
1425
1426         source_package->descriptor = mxf_resolve_strong_ref(mxf, &source_package->descriptor_ref, AnyType);
1427         if (source_package->descriptor) {
1428             if (source_package->descriptor->type == MultipleDescriptor) {
1429                 for (j = 0; j < source_package->descriptor->sub_descriptors_count; j++) {
1430                     MXFDescriptor *sub_descriptor = mxf_resolve_strong_ref(mxf, &source_package->descriptor->sub_descriptors_refs[j], Descriptor);
1431
1432                     if (!sub_descriptor) {
1433                         av_log(mxf->fc, AV_LOG_ERROR, "could not resolve sub descriptor strong ref\n");
1434                         continue;
1435                     }
1436                     if (sub_descriptor->linked_track_id == source_track->track_id) {
1437                         descriptor = sub_descriptor;
1438                         break;
1439                     }
1440                 }
1441             } else if (source_package->descriptor->type == Descriptor)
1442                 descriptor = source_package->descriptor;
1443         }
1444         if (!descriptor) {
1445             av_log(mxf->fc, AV_LOG_INFO, "source track %d: stream %d, no descriptor found\n", source_track->track_id, st->index);
1446             continue;
1447         }
1448         PRINT_KEY(mxf->fc, "essence codec     ul", descriptor->essence_codec_ul);
1449         PRINT_KEY(mxf->fc, "essence container ul", descriptor->essence_container_ul);
1450         essence_container_ul = &descriptor->essence_container_ul;
1451         /* HACK: replacing the original key with mxf_encrypted_essence_container
1452          * is not allowed according to s429-6, try to find correct information anyway */
1453         if (IS_KLV_KEY(essence_container_ul, mxf_encrypted_essence_container)) {
1454             av_log(mxf->fc, AV_LOG_INFO, "broken encrypted mxf file\n");
1455             for (k = 0; k < mxf->metadata_sets_count; k++) {
1456                 MXFMetadataSet *metadata = mxf->metadata_sets[k];
1457                 if (metadata->type == CryptoContext) {
1458                     essence_container_ul = &((MXFCryptoContext *)metadata)->source_container_ul;
1459                     break;
1460                 }
1461             }
1462         }
1463
1464         /* TODO: drop PictureEssenceCoding and SoundEssenceCompression, only check EssenceContainer */
1465         codec_ul = mxf_get_codec_ul(ff_mxf_codec_uls, &descriptor->essence_codec_ul);
1466         st->codec->codec_id = codec_ul->id;
1467         if (descriptor->extradata) {
1468             st->codec->extradata = descriptor->extradata;
1469             st->codec->extradata_size = descriptor->extradata_size;
1470         }
1471         if (st->codec->codec_type == AVMEDIA_TYPE_VIDEO) {
1472             source_track->intra_only = mxf_is_intra_only(descriptor);
1473             container_ul = mxf_get_codec_ul(mxf_picture_essence_container_uls, essence_container_ul);
1474             if (st->codec->codec_id == AV_CODEC_ID_NONE)
1475                 st->codec->codec_id = container_ul->id;
1476             st->codec->width = descriptor->width;
1477             /* Field height, not frame height */
1478             st->codec->height = descriptor->height;
1479             switch (descriptor->frame_layout) {
1480                 case SegmentedFrame:
1481                     /* This one is a weird layout I don't fully understand. */
1482                     av_log(mxf->fc, AV_LOG_INFO,
1483                            "SegmentedFrame layout isn't currently supported\n");
1484                     break;
1485                 case FullFrame:
1486                     break;
1487                 case OneField:
1488                     /* Every other line is stored and needs to be duplicated. */
1489                     av_log(mxf->fc, AV_LOG_INFO,
1490                            "OneField frame layout isn't currently supported\n");
1491                     break;
1492                     /* The correct thing to do here is fall through, but by
1493                      * breaking we might be able to decode some streams at half
1494                      * the vertical resolution, rather than not al all.
1495                      * It's also for compatibility with the old behavior. */
1496                 case SeparateFields:
1497                 case MixedFields:
1498                     /* Turn field height into frame height. */
1499                     st->codec->height *= 2;
1500                 default:
1501                     av_log(mxf->fc, AV_LOG_INFO,
1502                            "Unknown frame layout type: %d\n",
1503                            descriptor->frame_layout);
1504             }
1505             if (st->codec->codec_id == AV_CODEC_ID_RAWVIDEO) {
1506                 st->codec->pix_fmt = descriptor->pix_fmt;
1507                 if (st->codec->pix_fmt == AV_PIX_FMT_NONE) {
1508                     pix_fmt_ul = mxf_get_codec_ul(ff_mxf_pixel_format_uls,
1509                                                   &descriptor->essence_codec_ul);
1510                     st->codec->pix_fmt = pix_fmt_ul->id;
1511                     if (st->codec->pix_fmt == AV_PIX_FMT_NONE) {
1512                         /* support files created before RP224v10 by defaulting to UYVY422
1513                            if subsampling is 4:2:2 and component depth is 8-bit */
1514                         if (descriptor->horiz_subsampling == 2 &&
1515                             descriptor->vert_subsampling == 1 &&
1516                             descriptor->component_depth == 8) {
1517                             st->codec->pix_fmt = AV_PIX_FMT_UYVY422;
1518                         }
1519                     }
1520                 }
1521             }
1522             st->need_parsing = AVSTREAM_PARSE_HEADERS;
1523         } else if (st->codec->codec_type == AVMEDIA_TYPE_AUDIO) {
1524             container_ul = mxf_get_codec_ul(mxf_sound_essence_container_uls, essence_container_ul);
1525             if (st->codec->codec_id == AV_CODEC_ID_NONE)
1526                 st->codec->codec_id = container_ul->id;
1527             st->codec->channels = descriptor->channels;
1528             st->codec->bits_per_coded_sample = descriptor->bits_per_sample;
1529
1530             if (descriptor->sample_rate.den > 0)
1531                 st->codec->sample_rate = descriptor->sample_rate.num / descriptor->sample_rate.den;
1532
1533             /* TODO: implement AV_CODEC_ID_RAWAUDIO */
1534             if (st->codec->codec_id == AV_CODEC_ID_PCM_S16LE) {
1535                 if (descriptor->bits_per_sample > 16 && descriptor->bits_per_sample <= 24)
1536                     st->codec->codec_id = AV_CODEC_ID_PCM_S24LE;
1537                 else if (descriptor->bits_per_sample == 32)
1538                     st->codec->codec_id = AV_CODEC_ID_PCM_S32LE;
1539             } else if (st->codec->codec_id == AV_CODEC_ID_PCM_S16BE) {
1540                 if (descriptor->bits_per_sample > 16 && descriptor->bits_per_sample <= 24)
1541                     st->codec->codec_id = AV_CODEC_ID_PCM_S24BE;
1542                 else if (descriptor->bits_per_sample == 32)
1543                     st->codec->codec_id = AV_CODEC_ID_PCM_S32BE;
1544             } else if (st->codec->codec_id == AV_CODEC_ID_MP2) {
1545                 st->need_parsing = AVSTREAM_PARSE_FULL;
1546             }
1547         }
1548         if (st->codec->codec_type != AVMEDIA_TYPE_DATA && (*essence_container_ul)[15] > 0x01) {
1549             /* TODO: decode timestamps */
1550             st->need_parsing = AVSTREAM_PARSE_TIMESTAMPS;
1551         }
1552     }
1553
1554     ret = 0;
1555 fail_and_free:
1556     return ret;
1557 }
1558
1559 static const MXFMetadataReadTableEntry mxf_metadata_read_table[] = {
1560     { { 0x06,0x0E,0x2B,0x34,0x02,0x05,0x01,0x01,0x0d,0x01,0x02,0x01,0x01,0x05,0x01,0x00 }, mxf_read_primer_pack },
1561     { { 0x06,0x0E,0x2B,0x34,0x02,0x05,0x01,0x01,0x0d,0x01,0x02,0x01,0x01,0x02,0x01,0x00 }, mxf_read_partition_pack },
1562     { { 0x06,0x0E,0x2B,0x34,0x02,0x05,0x01,0x01,0x0d,0x01,0x02,0x01,0x01,0x02,0x02,0x00 }, mxf_read_partition_pack },
1563     { { 0x06,0x0E,0x2B,0x34,0x02,0x05,0x01,0x01,0x0d,0x01,0x02,0x01,0x01,0x02,0x03,0x00 }, mxf_read_partition_pack },
1564     { { 0x06,0x0E,0x2B,0x34,0x02,0x05,0x01,0x01,0x0d,0x01,0x02,0x01,0x01,0x02,0x04,0x00 }, mxf_read_partition_pack },
1565     { { 0x06,0x0E,0x2B,0x34,0x02,0x05,0x01,0x01,0x0d,0x01,0x02,0x01,0x01,0x03,0x01,0x00 }, mxf_read_partition_pack },
1566     { { 0x06,0x0E,0x2B,0x34,0x02,0x05,0x01,0x01,0x0d,0x01,0x02,0x01,0x01,0x03,0x02,0x00 }, mxf_read_partition_pack },
1567     { { 0x06,0x0E,0x2B,0x34,0x02,0x05,0x01,0x01,0x0d,0x01,0x02,0x01,0x01,0x03,0x03,0x00 }, mxf_read_partition_pack },
1568     { { 0x06,0x0E,0x2B,0x34,0x02,0x05,0x01,0x01,0x0d,0x01,0x02,0x01,0x01,0x03,0x04,0x00 }, mxf_read_partition_pack },
1569     { { 0x06,0x0E,0x2B,0x34,0x02,0x05,0x01,0x01,0x0d,0x01,0x02,0x01,0x01,0x04,0x02,0x00 }, mxf_read_partition_pack },
1570     { { 0x06,0x0E,0x2B,0x34,0x02,0x05,0x01,0x01,0x0d,0x01,0x02,0x01,0x01,0x04,0x04,0x00 }, mxf_read_partition_pack },
1571     { { 0x06,0x0E,0x2B,0x34,0x02,0x53,0x01,0x01,0x0d,0x01,0x01,0x01,0x01,0x01,0x18,0x00 }, mxf_read_content_storage, 0, AnyType },
1572     { { 0x06,0x0E,0x2B,0x34,0x02,0x53,0x01,0x01,0x0d,0x01,0x01,0x01,0x01,0x01,0x37,0x00 }, mxf_read_source_package, sizeof(MXFPackage), SourcePackage },
1573     { { 0x06,0x0E,0x2B,0x34,0x02,0x53,0x01,0x01,0x0d,0x01,0x01,0x01,0x01,0x01,0x36,0x00 }, mxf_read_material_package, sizeof(MXFPackage), MaterialPackage },
1574     { { 0x06,0x0E,0x2B,0x34,0x02,0x53,0x01,0x01,0x0d,0x01,0x01,0x01,0x01,0x01,0x0F,0x00 }, mxf_read_sequence, sizeof(MXFSequence), Sequence },
1575     { { 0x06,0x0E,0x2B,0x34,0x02,0x53,0x01,0x01,0x0d,0x01,0x01,0x01,0x01,0x01,0x11,0x00 }, mxf_read_source_clip, sizeof(MXFStructuralComponent), SourceClip },
1576     { { 0x06,0x0E,0x2B,0x34,0x02,0x53,0x01,0x01,0x0d,0x01,0x01,0x01,0x01,0x01,0x44,0x00 }, mxf_read_generic_descriptor, sizeof(MXFDescriptor), MultipleDescriptor },
1577     { { 0x06,0x0E,0x2B,0x34,0x02,0x53,0x01,0x01,0x0d,0x01,0x01,0x01,0x01,0x01,0x42,0x00 }, mxf_read_generic_descriptor, sizeof(MXFDescriptor), Descriptor }, /* Generic Sound */
1578     { { 0x06,0x0E,0x2B,0x34,0x02,0x53,0x01,0x01,0x0d,0x01,0x01,0x01,0x01,0x01,0x28,0x00 }, mxf_read_generic_descriptor, sizeof(MXFDescriptor), Descriptor }, /* CDCI */
1579     { { 0x06,0x0E,0x2B,0x34,0x02,0x53,0x01,0x01,0x0d,0x01,0x01,0x01,0x01,0x01,0x29,0x00 }, mxf_read_generic_descriptor, sizeof(MXFDescriptor), Descriptor }, /* RGBA */
1580     { { 0x06,0x0E,0x2B,0x34,0x02,0x53,0x01,0x01,0x0d,0x01,0x01,0x01,0x01,0x01,0x51,0x00 }, mxf_read_generic_descriptor, sizeof(MXFDescriptor), Descriptor }, /* MPEG 2 Video */
1581     { { 0x06,0x0E,0x2B,0x34,0x02,0x53,0x01,0x01,0x0d,0x01,0x01,0x01,0x01,0x01,0x48,0x00 }, mxf_read_generic_descriptor, sizeof(MXFDescriptor), Descriptor }, /* Wave */
1582     { { 0x06,0x0E,0x2B,0x34,0x02,0x53,0x01,0x01,0x0d,0x01,0x01,0x01,0x01,0x01,0x47,0x00 }, mxf_read_generic_descriptor, sizeof(MXFDescriptor), Descriptor }, /* AES3 */
1583     { { 0x06,0x0E,0x2B,0x34,0x02,0x53,0x01,0x01,0x0d,0x01,0x01,0x01,0x01,0x01,0x3A,0x00 }, mxf_read_track, sizeof(MXFTrack), Track }, /* Static Track */
1584     { { 0x06,0x0E,0x2B,0x34,0x02,0x53,0x01,0x01,0x0d,0x01,0x01,0x01,0x01,0x01,0x3B,0x00 }, mxf_read_track, sizeof(MXFTrack), Track }, /* Generic Track */
1585     { { 0x06,0x0E,0x2B,0x34,0x02,0x53,0x01,0x01,0x0d,0x01,0x04,0x01,0x02,0x02,0x00,0x00 }, mxf_read_cryptographic_context, sizeof(MXFCryptoContext), CryptoContext },
1586     { { 0x06,0x0E,0x2B,0x34,0x02,0x53,0x01,0x01,0x0d,0x01,0x02,0x01,0x01,0x10,0x01,0x00 }, mxf_read_index_table_segment, sizeof(MXFIndexTableSegment), IndexTableSegment },
1587     { { 0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00 }, NULL, 0, AnyType },
1588 };
1589
1590 static int mxf_read_local_tags(MXFContext *mxf, KLVPacket *klv, MXFMetadataReadFunc *read_child, int ctx_size, enum MXFMetadataSetType type)
1591 {
1592     AVIOContext *pb = mxf->fc->pb;
1593     MXFMetadataSet *ctx = ctx_size ? av_mallocz(ctx_size) : mxf;
1594     uint64_t klv_end = avio_tell(pb) + klv->length;
1595
1596     if (!ctx)
1597         return AVERROR(ENOMEM);
1598     while (avio_tell(pb) + 4 < klv_end && !pb->eof_reached) {
1599         int ret;
1600         int tag = avio_rb16(pb);
1601         int size = avio_rb16(pb); /* KLV specified by 0x53 */
1602         uint64_t next = avio_tell(pb) + size;
1603         UID uid = {0};
1604
1605         av_dlog(mxf->fc, "local tag %#04x size %d\n", tag, size);
1606         if (!size) { /* ignore empty tag, needed for some files with empty UMID tag */
1607             av_log(mxf->fc, AV_LOG_ERROR, "local tag %#04x with 0 size\n", tag);
1608             continue;
1609         }
1610         if (tag > 0x7FFF) { /* dynamic tag */
1611             int i;
1612             for (i = 0; i < mxf->local_tags_count; i++) {
1613                 int local_tag = AV_RB16(mxf->local_tags+i*18);
1614                 if (local_tag == tag) {
1615                     memcpy(uid, mxf->local_tags+i*18+2, 16);
1616                     av_dlog(mxf->fc, "local tag %#04x\n", local_tag);
1617                     PRINT_KEY(mxf->fc, "uid", uid);
1618                 }
1619             }
1620         }
1621         if (ctx_size && tag == 0x3C0A)
1622             avio_read(pb, ctx->uid, 16);
1623         else if ((ret = read_child(ctx, pb, tag, size, uid, -1)) < 0)
1624             return ret;
1625
1626         /* Accept the 64k local set limit being exceeded (Avid). Don't accept
1627          * it extending past the end of the KLV though (zzuf5.mxf). */
1628         if (avio_tell(pb) > klv_end) {
1629             av_log(mxf->fc, AV_LOG_ERROR,
1630                    "local tag %#04x extends past end of local set @ %#"PRIx64"\n",
1631                    tag, klv->offset);
1632             return AVERROR_INVALIDDATA;
1633         } else if (avio_tell(pb) <= next)   /* only seek forward, else this can loop for a long time */
1634             avio_seek(pb, next, SEEK_SET);
1635     }
1636     if (ctx_size) ctx->type = type;
1637     return ctx_size ? mxf_add_metadata_set(mxf, ctx) : 0;
1638 }
1639
1640 /**
1641  * Seeks to the previous partition, if possible
1642  * @return <= 0 if we should stop parsing, > 0 if we should keep going
1643  */
1644 static int mxf_seek_to_previous_partition(MXFContext *mxf)
1645 {
1646     AVIOContext *pb = mxf->fc->pb;
1647
1648     if (!mxf->current_partition ||
1649         mxf->run_in + mxf->current_partition->previous_partition <= mxf->last_forward_tell)
1650         return 0;   /* we've parsed all partitions */
1651
1652     /* seek to previous partition */
1653     avio_seek(pb, mxf->run_in + mxf->current_partition->previous_partition, SEEK_SET);
1654     mxf->current_partition = NULL;
1655
1656     av_dlog(mxf->fc, "seeking to previous partition\n");
1657
1658     return 1;
1659 }
1660
1661 /**
1662  * Called when essence is encountered
1663  * @return <= 0 if we should stop parsing, > 0 if we should keep going
1664  */
1665 static int mxf_parse_handle_essence(MXFContext *mxf)
1666 {
1667     AVIOContext *pb = mxf->fc->pb;
1668     int64_t ret;
1669
1670     if (mxf->parsing_backward) {
1671         return mxf_seek_to_previous_partition(mxf);
1672     } else {
1673         if (!mxf->footer_partition) {
1674             av_dlog(mxf->fc, "no footer\n");
1675             return 0;
1676         }
1677
1678         av_dlog(mxf->fc, "seeking to footer\n");
1679
1680         /* remember where we were so we don't end up seeking further back than this */
1681         mxf->last_forward_tell = avio_tell(pb);
1682
1683         if (!pb->seekable) {
1684             av_log(mxf->fc, AV_LOG_INFO, "file is not seekable - not parsing footer\n");
1685             return -1;
1686         }
1687
1688         /* seek to footer partition and parse backward */
1689         if ((ret = avio_seek(pb, mxf->run_in + mxf->footer_partition, SEEK_SET)) < 0) {
1690             av_log(mxf->fc, AV_LOG_ERROR, "failed to seek to footer @ 0x%"PRIx64" (%"PRId64") - partial file?\n",
1691                    mxf->run_in + mxf->footer_partition, ret);
1692             return ret;
1693         }
1694
1695         mxf->current_partition = NULL;
1696         mxf->parsing_backward = 1;
1697     }
1698
1699     return 1;
1700 }
1701
1702 /**
1703  * Called when the next partition or EOF is encountered
1704  * @return <= 0 if we should stop parsing, > 0 if we should keep going
1705  */
1706 static int mxf_parse_handle_partition_or_eof(MXFContext *mxf)
1707 {
1708     return mxf->parsing_backward ? mxf_seek_to_previous_partition(mxf) : 1;
1709 }
1710
1711 /**
1712  * Figure out the proper offset and length of the essence container
1713  * in each partition
1714  */
1715 static void mxf_compute_essence_containers(MXFContext *mxf)
1716 {
1717     int x;
1718
1719     /* everything is already correct */
1720     if (mxf->op == OPAtom)
1721         return;
1722
1723     for (x = 0; x < mxf->partitions_count; x++) {
1724         MXFPartition *p = &mxf->partitions[x];
1725
1726         if (!p->body_sid)
1727             continue;       /* BodySID == 0 -> no essence */
1728
1729         if (x >= mxf->partitions_count - 1)
1730             break;          /* last partition - can't compute length (and we don't need to) */
1731
1732         /* essence container spans to the next partition */
1733         p->essence_length = mxf->partitions[x+1].this_partition - p->essence_offset;
1734
1735         if (p->essence_length < 0) {
1736             /* next ThisPartition < essence_offset */
1737             p->essence_length = 0;
1738             av_log(mxf->fc, AV_LOG_ERROR,
1739                    "partition %i: bad ThisPartition = %"PRIX64"\n",
1740                    x+1, mxf->partitions[x+1].this_partition);
1741         }
1742     }
1743 }
1744
1745 static int64_t round_to_kag(int64_t position, int kag_size)
1746 {
1747     /* TODO: account for run-in? the spec isn't clear whether KAG should account for it */
1748     /* NOTE: kag_size may be any integer between 1 - 2^10 */
1749     int64_t ret = (position / kag_size) * kag_size;
1750     return ret == position ? ret : ret + kag_size;
1751 }
1752
1753 static inline void compute_partition_essence_offset(AVFormatContext *s,
1754                                                     MXFContext *mxf,
1755                                                     KLVPacket *klv)
1756 {
1757     MXFPartition *cur_part = mxf->current_partition;
1758     /* for OP1a we compute essence_offset
1759      * for OPAtom we point essence_offset after the KL
1760      *     (usually op1a_essence_offset + 20 or 25)
1761      * TODO: for OP1a we could eliminate this entire if statement, always
1762      *       stopping parsing at op1a_essence_offset
1763      *       for OPAtom we still need the actual essence_offset though
1764      *       (the KL's length can vary)
1765      */
1766     int64_t op1a_essence_offset =
1767         round_to_kag(cur_part->this_partition + cur_part->pack_length,
1768                      cur_part->kag_size) +
1769         round_to_kag(cur_part->header_byte_count, cur_part->kag_size) +
1770         round_to_kag(cur_part->index_byte_count, cur_part->kag_size);
1771
1772     if (mxf->op == OPAtom) {
1773         /* point essence_offset to the actual data
1774          * OPAtom has all the essence in one big KLV
1775          */
1776         cur_part->essence_offset = avio_tell(s->pb);
1777         cur_part->essence_length = klv->length;
1778     } else {
1779         /* NOTE: op1a_essence_offset may be less than to klv.offset
1780          * (C0023S01.mxf)  */
1781         cur_part->essence_offset = op1a_essence_offset;
1782     }
1783 }
1784
1785 static int is_pcm(enum AVCodecID codec_id)
1786 {
1787     /* we only care about "normal" PCM codecs until we get samples */
1788     return codec_id >= AV_CODEC_ID_PCM_S16LE && codec_id < AV_CODEC_ID_PCM_S24DAUD;
1789 }
1790
1791 /**
1792  * Deal with the case where for some audio atoms EditUnitByteCount is
1793  * very small (2, 4..). In those cases we should read more than one
1794  * sample per call to mxf_read_packet().
1795  */
1796 static void mxf_handle_small_eubc(AVFormatContext *s)
1797 {
1798     MXFContext *mxf = s->priv_data;
1799
1800     /* assuming non-OPAtom == frame wrapped
1801      * no sane writer would wrap 2 byte PCM packets with 20 byte headers.. */
1802     if (mxf->op != OPAtom)
1803         return;
1804
1805     /* expect PCM with exactly one index table segment and a small (< 32) EUBC */
1806     if (s->nb_streams != 1                                     ||
1807         s->streams[0]->codec->codec_type != AVMEDIA_TYPE_AUDIO ||
1808         !is_pcm(s->streams[0]->codec->codec_id)                ||
1809         mxf->nb_index_tables != 1                              ||
1810         mxf->index_tables[0].nb_segments != 1                  ||
1811         mxf->index_tables[0].segments[0]->edit_unit_byte_count >= 32)
1812         return;
1813
1814     /* arbitrarily default to 48 kHz PAL audio frame size */
1815     /* TODO: We could compute this from the ratio between the audio
1816      *       and video edit rates for 48 kHz NTSC we could use the
1817      *       1802-1802-1802-1802-1801 pattern. */
1818     mxf->edit_units_per_packet = 1920;
1819 }
1820
1821 static int mxf_read_header(AVFormatContext *s)
1822 {
1823     MXFContext *mxf = s->priv_data;
1824     KLVPacket klv;
1825     int64_t essence_offset = 0;
1826     int ret;
1827
1828     mxf->last_forward_tell = INT64_MAX;
1829     mxf->edit_units_per_packet = 1;
1830
1831     if (!mxf_read_sync(s->pb, mxf_header_partition_pack_key, 14)) {
1832         av_log(s, AV_LOG_ERROR, "could not find header partition pack key\n");
1833         return AVERROR_INVALIDDATA;
1834     }
1835     avio_seek(s->pb, -14, SEEK_CUR);
1836     mxf->fc = s;
1837     mxf->run_in = avio_tell(s->pb);
1838
1839     while (!s->pb->eof_reached) {
1840         const MXFMetadataReadTableEntry *metadata;
1841
1842         if (klv_read_packet(&klv, s->pb) < 0) {
1843             /* EOF - seek to previous partition or stop */
1844             if(mxf_parse_handle_partition_or_eof(mxf) <= 0)
1845                 break;
1846             else
1847                 continue;
1848         }
1849
1850         PRINT_KEY(s, "read header", klv.key);
1851         av_dlog(s, "size %"PRIu64" offset %#"PRIx64"\n", klv.length, klv.offset);
1852         if (IS_KLV_KEY(klv.key, mxf_encrypted_triplet_key) ||
1853             IS_KLV_KEY(klv.key, mxf_essence_element_key) ||
1854             IS_KLV_KEY(klv.key, mxf_avid_essence_element_key) ||
1855             IS_KLV_KEY(klv.key, mxf_system_item_key)) {
1856
1857             if (!mxf->current_partition) {
1858                 av_log(mxf->fc, AV_LOG_ERROR,
1859                        "found essence prior to first PartitionPack\n");
1860                 return AVERROR_INVALIDDATA;
1861             }
1862
1863             if (!mxf->current_partition->essence_offset) {
1864                 compute_partition_essence_offset(s, mxf, &klv);
1865             }
1866
1867             if (!essence_offset)
1868                 essence_offset = klv.offset;
1869
1870             /* seek to footer, previous partition or stop */
1871             if (mxf_parse_handle_essence(mxf) <= 0)
1872                 break;
1873             continue;
1874         } else if (!memcmp(klv.key, mxf_header_partition_pack_key, 13) &&
1875                    klv.key[13] >= 2 && klv.key[13] <= 4 && mxf->current_partition) {
1876             /* next partition pack - keep going, seek to previous partition or stop */
1877             if(mxf_parse_handle_partition_or_eof(mxf) <= 0)
1878                 break;
1879             else if (mxf->parsing_backward)
1880                 continue;
1881             /* we're still parsing forward. proceed to parsing this partition pack */
1882         }
1883
1884         for (metadata = mxf_metadata_read_table; metadata->read; metadata++) {
1885             if (IS_KLV_KEY(klv.key, metadata->key)) {
1886                 int res;
1887                 if (klv.key[5] == 0x53) {
1888                     res = mxf_read_local_tags(mxf, &klv, metadata->read, metadata->ctx_size, metadata->type);
1889                 } else {
1890                     uint64_t next = avio_tell(s->pb) + klv.length;
1891                     res = metadata->read(mxf, s->pb, 0, klv.length, klv.key, klv.offset);
1892
1893                     /* only seek forward, else this can loop for a long time */
1894                     if (avio_tell(s->pb) > next) {
1895                         av_log(s, AV_LOG_ERROR, "read past end of KLV @ %#"PRIx64"\n",
1896                                klv.offset);
1897                         return AVERROR_INVALIDDATA;
1898                     }
1899
1900                     avio_seek(s->pb, next, SEEK_SET);
1901                 }
1902                 if (res < 0) {
1903                     av_log(s, AV_LOG_ERROR, "error reading header metadata\n");
1904                     return res;
1905                 }
1906                 break;
1907             }
1908         }
1909         if (!metadata->read)
1910             avio_skip(s->pb, klv.length);
1911     }
1912     /* FIXME avoid seek */
1913     if (!essence_offset)  {
1914         av_log(s, AV_LOG_ERROR, "no essence\n");
1915         return AVERROR_INVALIDDATA;
1916     }
1917     avio_seek(s->pb, essence_offset, SEEK_SET);
1918
1919     mxf_compute_essence_containers(mxf);
1920
1921     /* we need to do this before computing the index tables
1922      * to be able to fill in zero IndexDurations with st->duration */
1923     if ((ret = mxf_parse_structural_metadata(mxf)) < 0)
1924         return ret;
1925
1926     if ((ret = mxf_compute_index_tables(mxf)) < 0)
1927         return ret;
1928
1929     if (mxf->nb_index_tables > 1) {
1930         /* TODO: look up which IndexSID to use via EssenceContainerData */
1931         av_log(mxf->fc, AV_LOG_INFO, "got %i index tables - only the first one (IndexSID %i) will be used\n",
1932                mxf->nb_index_tables, mxf->index_tables[0].index_sid);
1933     } else if (mxf->nb_index_tables == 0 && mxf->op == OPAtom) {
1934         av_log(mxf->fc, AV_LOG_ERROR, "cannot demux OPAtom without an index\n");
1935         return AVERROR_INVALIDDATA;
1936     }
1937
1938     mxf_handle_small_eubc(s);
1939
1940     return 0;
1941 }
1942
1943 /**
1944  * Sets mxf->current_edit_unit based on what offset we're currently at.
1945  * @return next_ofs if OK, <0 on error
1946  */
1947 static int64_t mxf_set_current_edit_unit(MXFContext *mxf, int64_t current_offset)
1948 {
1949     int64_t last_ofs = -1, next_ofs = -1;
1950     MXFIndexTable *t = &mxf->index_tables[0];
1951
1952     /* this is called from the OP1a demuxing logic, which means there
1953      * may be no index tables */
1954     if (mxf->nb_index_tables <= 0)
1955         return -1;
1956
1957     /* find mxf->current_edit_unit so that the next edit unit starts ahead
1958      * of current_offset */
1959     while (mxf->current_edit_unit >= 0) {
1960         if (mxf_edit_unit_absolute_offset(mxf, t, mxf->current_edit_unit + 1,
1961                                           NULL, &next_ofs, 0) < 0)
1962             return -1;
1963
1964         if (next_ofs <= last_ofs) {
1965             /* large next_ofs didn't change or current_edit_unit wrapped
1966              * around this fixes the infinite loop on zzuf3.mxf */
1967             av_log(mxf->fc, AV_LOG_ERROR,
1968                    "next_ofs didn't change. not deriving packet timestamps\n");
1969             return -1;
1970         }
1971
1972         if (next_ofs > current_offset)
1973             break;
1974
1975         last_ofs = next_ofs;
1976         mxf->current_edit_unit++;
1977     }
1978
1979     /* not checking mxf->current_edit_unit >= t->nb_ptses here since CBR files
1980      * may lack IndexEntryArrays */
1981     if (mxf->current_edit_unit < 0)
1982         return -1;
1983
1984     return next_ofs;
1985 }
1986
1987 static int mxf_read_packet_old(AVFormatContext *s, AVPacket *pkt)
1988 {
1989     KLVPacket klv;
1990     MXFContext *mxf = s->priv_data;
1991
1992     while (!s->pb->eof_reached) {
1993         if (klv_read_packet(&klv, s->pb) < 0)
1994             return -1;
1995         PRINT_KEY(s, "read packet", klv.key);
1996         av_dlog(s, "size %"PRIu64" offset %#"PRIx64"\n", klv.length, klv.offset);
1997         if (IS_KLV_KEY(klv.key, mxf_encrypted_triplet_key)) {
1998             int res = mxf_decrypt_triplet(s, pkt, &klv);
1999             if (res < 0) {
2000                 av_log(s, AV_LOG_ERROR, "invalid encoded triplet\n");
2001                 return -1;
2002             }
2003             return 0;
2004         }
2005         if (IS_KLV_KEY(klv.key, mxf_essence_element_key) ||
2006             IS_KLV_KEY(klv.key, mxf_avid_essence_element_key)) {
2007             int index = mxf_get_stream_index(s, &klv);
2008             int64_t next_ofs, next_klv;
2009             AVStream *st;
2010             MXFTrack *track;
2011
2012             if (index < 0) {
2013                 av_log(s, AV_LOG_ERROR, "error getting stream index %d\n", AV_RB32(klv.key+12));
2014                 goto skip;
2015             }
2016
2017             st = s->streams[index];
2018             track = st->priv_data;
2019
2020             if (s->streams[index]->discard == AVDISCARD_ALL)
2021                 goto skip;
2022
2023             next_klv = avio_tell(s->pb) + klv.length;
2024             next_ofs = mxf_set_current_edit_unit(mxf, klv.offset);
2025
2026             if (next_ofs >= 0 && next_klv > next_ofs) {
2027                 /* if this check is hit then it's possible OPAtom was treated
2028                  * as OP1a truncate the packet since it's probably very large
2029                  * (>2 GiB is common) */
2030                 av_log_ask_for_sample(s,
2031                                       "KLV for edit unit %i extends into next "
2032                                       "edit unit - OPAtom misinterpreted as "
2033                                       "OP1a?\n",
2034                                       mxf->current_edit_unit);
2035                 klv.length = next_ofs - avio_tell(s->pb);
2036             }
2037
2038             /* check for 8 channels AES3 element */
2039             if (klv.key[12] == 0x06 && klv.key[13] == 0x01 && klv.key[14] == 0x10) {
2040                 if (mxf_get_d10_aes3_packet(s->pb, s->streams[index], pkt, klv.length) < 0) {
2041                     av_log(s, AV_LOG_ERROR, "error reading D-10 aes3 frame\n");
2042                     return -1;
2043                 }
2044             } else {
2045                 int ret = av_get_packet(s->pb, pkt, klv.length);
2046                 if (ret < 0)
2047                     return ret;
2048             }
2049             pkt->stream_index = index;
2050             pkt->pos = klv.offset;
2051
2052             if (s->streams[index]->codec->codec_type == AVMEDIA_TYPE_VIDEO && next_ofs >= 0) {
2053                 /* mxf->current_edit_unit good - see if we have an
2054                  * index table to derive timestamps from */
2055                 MXFIndexTable *t = &mxf->index_tables[0];
2056
2057                 if (mxf->nb_index_tables >= 1 &&
2058                     mxf->current_edit_unit < t->nb_ptses) {
2059                     pkt->dts = mxf->current_edit_unit + t->first_dts;
2060                     pkt->pts = t->ptses[mxf->current_edit_unit];
2061                 } else if (track->intra_only) {
2062                     /* intra-only -> PTS = EditUnit.
2063                      * let utils.c figure out DTS since it can be
2064                      * < PTS if low_delay = 0 (Sony IMX30) */
2065                     pkt->pts = mxf->current_edit_unit;
2066                 }
2067             }
2068
2069             /* seek for truncated packets */
2070             avio_seek(s->pb, next_klv, SEEK_SET);
2071
2072             return 0;
2073         } else
2074         skip:
2075             avio_skip(s->pb, klv.length);
2076     }
2077     return AVERROR_EOF;
2078 }
2079
2080 static int mxf_read_packet(AVFormatContext *s, AVPacket *pkt)
2081 {
2082     MXFContext *mxf = s->priv_data;
2083     int ret, size;
2084     int64_t ret64, pos, next_pos;
2085     AVStream *st;
2086     MXFIndexTable *t;
2087     int edit_units;
2088
2089     if (mxf->op != OPAtom)
2090         return mxf_read_packet_old(s, pkt);
2091
2092     /* OPAtom - clip wrapped demuxing */
2093     /* NOTE: mxf_read_header() makes sure nb_index_tables > 0 for OPAtom */
2094     st = s->streams[0];
2095     t = &mxf->index_tables[0];
2096
2097     if (mxf->current_edit_unit >= st->duration)
2098         return AVERROR_EOF;
2099
2100     edit_units = FFMIN(mxf->edit_units_per_packet, st->duration - mxf->current_edit_unit);
2101
2102     if ((ret = mxf_edit_unit_absolute_offset(mxf, t, mxf->current_edit_unit, NULL, &pos, 1)) < 0)
2103         return ret;
2104
2105     /* compute size by finding the next edit unit or the end of the essence container
2106      * not pretty, but it works */
2107     if ((ret = mxf_edit_unit_absolute_offset(mxf, t, mxf->current_edit_unit + edit_units, NULL, &next_pos, 0)) < 0 &&
2108         (next_pos = mxf_essence_container_end(mxf, t->body_sid)) <= 0) {
2109         av_log(s, AV_LOG_ERROR, "unable to compute the size of the last packet\n");
2110         return AVERROR_INVALIDDATA;
2111     }
2112
2113     if ((size = next_pos - pos) <= 0) {
2114         av_log(s, AV_LOG_ERROR, "bad size: %i\n", size);
2115         return AVERROR_INVALIDDATA;
2116     }
2117
2118     if ((ret64 = avio_seek(s->pb, pos, SEEK_SET)) < 0)
2119         return ret64;
2120
2121         if ((ret = av_get_packet(s->pb, pkt, size)) != size)
2122             return ret < 0 ? ret : AVERROR_EOF;
2123
2124     if (st->codec->codec_type == AVMEDIA_TYPE_VIDEO && t->ptses &&
2125         mxf->current_edit_unit >= 0 && mxf->current_edit_unit < t->nb_ptses) {
2126         pkt->dts = mxf->current_edit_unit + t->first_dts;
2127         pkt->pts = t->ptses[mxf->current_edit_unit];
2128     }
2129
2130     pkt->stream_index = 0;
2131     mxf->current_edit_unit += edit_units;
2132
2133     return 0;
2134 }
2135
2136
2137 static int mxf_read_close(AVFormatContext *s)
2138 {
2139     MXFContext *mxf = s->priv_data;
2140     MXFIndexTableSegment *seg;
2141     int i;
2142
2143     av_freep(&mxf->packages_refs);
2144
2145     for (i = 0; i < s->nb_streams; i++)
2146         s->streams[i]->priv_data = NULL;
2147
2148     for (i = 0; i < mxf->metadata_sets_count; i++) {
2149         switch (mxf->metadata_sets[i]->type) {
2150         case MultipleDescriptor:
2151             av_freep(&((MXFDescriptor *)mxf->metadata_sets[i])->sub_descriptors_refs);
2152             break;
2153         case Sequence:
2154             av_freep(&((MXFSequence *)mxf->metadata_sets[i])->structural_components_refs);
2155             break;
2156         case SourcePackage:
2157         case MaterialPackage:
2158             av_freep(&((MXFPackage *)mxf->metadata_sets[i])->tracks_refs);
2159             break;
2160         case IndexTableSegment:
2161             seg = (MXFIndexTableSegment *)mxf->metadata_sets[i];
2162             av_freep(&seg->temporal_offset_entries);
2163             av_freep(&seg->flag_entries);
2164             av_freep(&seg->stream_offset_entries);
2165             break;
2166         default:
2167             break;
2168         }
2169         av_freep(&mxf->metadata_sets[i]);
2170     }
2171     av_freep(&mxf->partitions);
2172     av_freep(&mxf->metadata_sets);
2173     av_freep(&mxf->aesc);
2174     av_freep(&mxf->local_tags);
2175
2176     for (i = 0; i < mxf->nb_index_tables; i++) {
2177         av_freep(&mxf->index_tables[i].segments);
2178         av_freep(&mxf->index_tables[i].ptses);
2179         av_freep(&mxf->index_tables[i].fake_index);
2180     }
2181     av_freep(&mxf->index_tables);
2182
2183     return 0;
2184 }
2185
2186 static int mxf_probe(AVProbeData *p) {
2187     uint8_t *bufp = p->buf;
2188     uint8_t *end = p->buf + p->buf_size;
2189
2190     if (p->buf_size < sizeof(mxf_header_partition_pack_key))
2191         return 0;
2192
2193     /* Must skip Run-In Sequence and search for MXF header partition pack key SMPTE 377M 5.5 */
2194     end -= sizeof(mxf_header_partition_pack_key);
2195     for (; bufp < end; bufp++) {
2196         if (IS_KLV_KEY(bufp, mxf_header_partition_pack_key))
2197             return AVPROBE_SCORE_MAX;
2198     }
2199     return 0;
2200 }
2201
2202 /* rudimentary byte seek */
2203 /* XXX: use MXF Index */
2204 static int mxf_read_seek(AVFormatContext *s, int stream_index, int64_t sample_time, int flags)
2205 {
2206     AVStream *st = s->streams[stream_index];
2207     int64_t seconds;
2208     MXFContext* mxf = s->priv_data;
2209     int64_t seekpos;
2210     int ret;
2211     MXFIndexTable *t;
2212
2213     if (mxf->nb_index_tables <= 0) {
2214     if (!s->bit_rate)
2215         return AVERROR_INVALIDDATA;
2216     if (sample_time < 0)
2217         sample_time = 0;
2218     seconds = av_rescale(sample_time, st->time_base.num, st->time_base.den);
2219
2220     if ((ret = avio_seek(s->pb, (s->bit_rate * seconds) >> 3, SEEK_SET)) < 0)
2221         return ret;
2222     ff_update_cur_dts(s, st, sample_time);
2223     } else {
2224         t = &mxf->index_tables[0];
2225
2226         /* clamp above zero, else ff_index_search_timestamp() returns negative
2227          * this also means we allow seeking before the start */
2228         sample_time = FFMAX(sample_time, 0);
2229
2230         if (t->fake_index) {
2231             /* behave as if we have a proper index */
2232             if ((sample_time = ff_index_search_timestamp(t->fake_index, t->nb_ptses, sample_time, flags)) < 0)
2233                 return sample_time;
2234         } else {
2235             /* no IndexEntryArray (one or more CBR segments)
2236              * make sure we don't seek past the end */
2237             sample_time = FFMIN(sample_time, st->duration - 1);
2238         }
2239
2240         if ((ret = mxf_edit_unit_absolute_offset(mxf, t, sample_time, &sample_time, &seekpos, 1)) << 0)
2241             return ret;
2242
2243         ff_update_cur_dts(s, st, sample_time);
2244         mxf->current_edit_unit = sample_time;
2245         avio_seek(s->pb, seekpos, SEEK_SET);
2246     }
2247     return 0;
2248 }
2249
2250 AVInputFormat ff_mxf_demuxer = {
2251     .name           = "mxf",
2252     .long_name      = NULL_IF_CONFIG_SMALL("MXF (Material eXchange Format)"),
2253     .priv_data_size = sizeof(MXFContext),
2254     .read_probe     = mxf_probe,
2255     .read_header    = mxf_read_header,
2256     .read_packet    = mxf_read_packet,
2257     .read_close     = mxf_read_close,
2258     .read_seek      = mxf_read_seek,
2259 };