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