svq3: use hpeldsp instead of dsputil for half-pel functions.
[ffmpeg.git] / libavcodec / exr.c
1 /*
2  * OpenEXR (.exr) image decoder
3  * Copyright (c) 2009 Jimmy Christensen
4  *
5  * This file is part of FFmpeg.
6  *
7  * FFmpeg is free software; you can redistribute it and/or
8  * modify it under the terms of the GNU Lesser General Public
9  * License as published by the Free Software Foundation; either
10  * version 2.1 of the License, or (at your option) any later version.
11  *
12  * FFmpeg is distributed in the hope that it will be useful,
13  * but WITHOUT ANY WARRANTY; without even the implied warranty of
14  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
15  * Lesser General Public License for more details.
16  *
17  * You should have received a copy of the GNU Lesser General Public
18  * License along with FFmpeg; if not, write to the Free Software
19  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
20  */
21
22 /**
23  * @file
24  * OpenEXR decoder
25  * @author Jimmy Christensen
26  *
27  * For more information on the OpenEXR format, visit:
28  *  http://openexr.com/
29  *
30  * exr_flt2uint() and exr_halflt2uint() is credited to  Reimar Döffinger
31  */
32
33 #include <zlib.h>
34
35 #include "avcodec.h"
36 #include "bytestream.h"
37 #include "mathops.h"
38 #include "thread.h"
39 #include "libavutil/imgutils.h"
40 #include "libavutil/avassert.h"
41
42 enum ExrCompr {
43     EXR_RAW   = 0,
44     EXR_RLE   = 1,
45     EXR_ZIP1  = 2,
46     EXR_ZIP16 = 3,
47     EXR_PIZ   = 4,
48     EXR_PXR24 = 5,
49     EXR_B44   = 6,
50     EXR_B44A  = 7,
51 };
52
53 enum ExrPixelType {
54     EXR_UINT,
55     EXR_HALF,
56     EXR_FLOAT
57 };
58
59 typedef struct EXRChannel {
60     int               xsub, ysub;
61     enum ExrPixelType pixel_type;
62 } EXRChannel;
63
64 typedef struct EXRThreadData {
65     uint8_t *uncompressed_data;
66     int uncompressed_size;
67
68     uint8_t *tmp;
69     int tmp_size;
70 } EXRThreadData;
71
72 typedef struct EXRContext {
73     AVFrame *picture;
74     int compr;
75     enum ExrPixelType pixel_type;
76     int channel_offsets[4]; // 0 = red, 1 = green, 2 = blue and 3 = alpha
77     const AVPixFmtDescriptor *desc;
78
79     uint32_t xmax, xmin;
80     uint32_t ymax, ymin;
81     uint32_t xdelta, ydelta;
82
83     int ysize;
84
85     uint64_t scan_line_size;
86     int scan_lines_per_block;
87
88     const uint8_t *buf, *table;
89     int buf_size;
90
91     EXRChannel *channels;
92     int nb_channels;
93
94     EXRThreadData *thread_data;
95     int thread_data_size;
96 } EXRContext;
97
98 /**
99  * Converts from 32-bit float as uint32_t to uint16_t
100  *
101  * @param v 32-bit float
102  * @return normalized 16-bit unsigned int
103  */
104 static inline uint16_t exr_flt2uint(uint32_t v)
105 {
106     unsigned int exp = v >> 23;
107     // "HACK": negative values result in exp<  0, so clipping them to 0
108     // is also handled by this condition, avoids explicit check for sign bit.
109     if (exp<= 127 + 7 - 24) // we would shift out all bits anyway
110         return 0;
111     if (exp >= 127)
112         return 0xffff;
113     v &= 0x007fffff;
114     return (v + (1 << 23)) >> (127 + 7 - exp);
115 }
116
117 /**
118  * Converts from 16-bit float as uint16_t to uint16_t
119  *
120  * @param v 16-bit float
121  * @return normalized 16-bit unsigned int
122  */
123 static inline uint16_t exr_halflt2uint(uint16_t v)
124 {
125     unsigned exp = 14 - (v >> 10);
126     if (exp >= 14) {
127         if (exp == 14) return (v >> 9) & 1;
128         else           return (v & 0x8000) ? 0 : 0xffff;
129     }
130     v <<= 6;
131     return (v + (1 << 16)) >> (exp + 1);
132 }
133
134 /**
135  * Gets the size of the header variable
136  *
137  * @param **buf the current pointer location in the header where
138  * the variable data starts
139  * @param *buf_end pointer location of the end of the buffer
140  * @return size of variable data
141  */
142 static unsigned int get_header_variable_length(const uint8_t **buf,
143                                                const uint8_t *buf_end)
144 {
145     unsigned int variable_buffer_data_size = bytestream_get_le32(buf);
146     if (variable_buffer_data_size >= buf_end - *buf)
147         return 0;
148     return variable_buffer_data_size;
149 }
150
151 /**
152  * Checks if the variable name corresponds with it's data type
153  *
154  * @param *avctx the AVCodecContext
155  * @param **buf the current pointer location in the header where
156  * the variable name starts
157  * @param *buf_end pointer location of the end of the buffer
158  * @param *value_name name of the varible to check
159  * @param *value_type type of the varible to check
160  * @param minimum_length minimum length of the variable data
161  * @param variable_buffer_data_size variable length read from the header
162  * after it's checked
163  * @return negative if variable is invalid
164  */
165 static int check_header_variable(AVCodecContext *avctx,
166                                               const uint8_t **buf,
167                                               const uint8_t *buf_end,
168                                               const char *value_name,
169                                               const char *value_type,
170                                               unsigned int minimum_length,
171                                               unsigned int *variable_buffer_data_size)
172 {
173     if (buf_end - *buf >= minimum_length && !strcmp(*buf, value_name)) {
174         *buf += strlen(value_name)+1;
175         if (!strcmp(*buf, value_type)) {
176             *buf += strlen(value_type)+1;
177             *variable_buffer_data_size = get_header_variable_length(buf, buf_end);
178             if (!*variable_buffer_data_size)
179                 av_log(avctx, AV_LOG_ERROR, "Incomplete header\n");
180             return 1;
181         }
182         *buf -= strlen(value_name)+1;
183         av_log(avctx, AV_LOG_WARNING, "Unknown data type for header variable %s\n", value_name);
184     }
185     return -1;
186 }
187
188 static void predictor(uint8_t *src, int size)
189 {
190     uint8_t *t = src + 1;
191     uint8_t *stop = src + size;
192
193     while (t < stop) {
194         int d = (int)t[-1] + (int)t[0] - 128;
195         t[0] = d;
196         ++t;
197     }
198 }
199
200 static void reorder_pixels(uint8_t *src, uint8_t *dst, int size)
201 {
202     const int8_t *t1 = src;
203     const int8_t *t2 = src + (size + 1) / 2;
204     int8_t *s = dst;
205     int8_t *stop = s + size;
206
207     while (1) {
208         if (s < stop)
209             *(s++) = *(t1++);
210         else
211             break;
212
213         if (s < stop)
214             *(s++) = *(t2++);
215         else
216             break;
217     }
218 }
219
220 static int zip_uncompress(const uint8_t *src, int compressed_size,
221                           int uncompressed_size, EXRThreadData *td)
222 {
223     unsigned long dest_len = uncompressed_size;
224
225     if (uncompress(td->tmp, &dest_len, src, compressed_size) != Z_OK ||
226         dest_len != uncompressed_size)
227         return AVERROR(EINVAL);
228
229     predictor(td->tmp, uncompressed_size);
230     reorder_pixels(td->tmp, td->uncompressed_data, uncompressed_size);
231
232     return 0;
233 }
234
235 static int rle_uncompress(const uint8_t *src, int compressed_size,
236                           int uncompressed_size, EXRThreadData *td)
237 {
238     int8_t *d = (int8_t *)td->tmp;
239     const int8_t *s = (const int8_t *)src;
240     int ssize = compressed_size;
241     int dsize = uncompressed_size;
242     int8_t *dend = d + dsize;
243     int count;
244
245     while (ssize > 0) {
246         count = *s++;
247
248         if (count < 0) {
249             count = -count;
250
251             if ((dsize -= count    ) < 0 ||
252                 (ssize -= count + 1) < 0)
253                 return -1;
254
255             while (count--)
256                 *d++ = *s++;
257         } else {
258             count++;
259
260             if ((dsize -= count) < 0 ||
261                 (ssize -= 2    ) < 0)
262                 return -1;
263
264             while (count--)
265                 *d++ = *s;
266
267             s++;
268         }
269     }
270
271     if (dend != d)
272         return AVERROR_INVALIDDATA;
273
274     predictor(td->tmp, uncompressed_size);
275     reorder_pixels(td->tmp, td->uncompressed_data, uncompressed_size);
276
277     return 0;
278 }
279
280 static int pxr24_uncompress(EXRContext *s, const uint8_t *src,
281                             int compressed_size, int uncompressed_size,
282                             EXRThreadData *td)
283 {
284     unsigned long dest_len = uncompressed_size;
285     const uint8_t *in = td->tmp;
286     uint8_t *out;
287     int c, i, j;
288
289     if (uncompress(td->tmp, &dest_len, src, compressed_size) != Z_OK ||
290         dest_len != uncompressed_size)
291         return AVERROR(EINVAL);
292
293     out = td->uncompressed_data;
294     for (i = 0; i < s->ysize; i++) {
295         for (c = 0; c < s->nb_channels; c++) {
296             EXRChannel *channel = &s->channels[c];
297             const uint8_t *ptr[4];
298             uint32_t pixel = 0;
299
300             switch (channel->pixel_type) {
301             case EXR_FLOAT:
302                 ptr[0] = in;
303                 ptr[1] = ptr[0] + s->xdelta;
304                 ptr[2] = ptr[1] + s->xdelta;
305                 in = ptr[2] + s->xdelta;
306
307                 for (j = 0; j < s->xdelta; ++j) {
308                     uint32_t diff = (*(ptr[0]++) << 24) |
309                                     (*(ptr[1]++) << 16) |
310                                     (*(ptr[2]++) <<  8);
311                     pixel += diff;
312                     bytestream_put_le32(&out, pixel);
313                 }
314                 break;
315             case EXR_HALF:
316                 ptr[0] = in;
317                 ptr[1] = ptr[0] + s->xdelta;
318                 in = ptr[1] + s->xdelta;
319                 for (j = 0; j < s->xdelta; j++) {
320                     uint32_t diff = (*(ptr[0]++) << 8) | *(ptr[1]++);
321
322                     pixel += diff;
323                     bytestream_put_le16(&out, pixel);
324                 }
325                 break;
326             default:
327                 av_assert1(0);
328             }
329         }
330     }
331
332     return 0;
333 }
334
335 static int decode_block(AVCodecContext *avctx, void *tdata,
336                         int jobnr, int threadnr)
337 {
338     EXRContext *s = avctx->priv_data;
339     AVFrame *const p = s->picture;
340     EXRThreadData *td = &s->thread_data[threadnr];
341     const uint8_t *channel_buffer[4] = { 0 };
342     const uint8_t *buf = s->buf;
343     uint64_t line_offset, uncompressed_size;
344     uint32_t xdelta = s->xdelta;
345     uint16_t *ptr_x;
346     uint8_t *ptr;
347     int32_t data_size, line;
348     const uint8_t *src;
349     int axmax = (avctx->width - (s->xmax + 1)) * 2 * s->desc->nb_components;
350     int bxmin = s->xmin * 2 * s->desc->nb_components;
351     int ret, i, x, buf_size = s->buf_size;
352
353     line_offset = AV_RL64(s->table + jobnr * 8);
354     // Check if the buffer has the required bytes needed from the offset
355     if (line_offset > buf_size - 8)
356         return AVERROR_INVALIDDATA;
357
358     src = buf + line_offset + 8;
359     line = AV_RL32(src - 8);
360     if (line < s->ymin || line > s->ymax)
361         return AVERROR_INVALIDDATA;
362
363     data_size = AV_RL32(src - 4);
364     if (data_size <= 0 || data_size > buf_size)
365         return AVERROR_INVALIDDATA;
366
367     s->ysize = FFMIN(s->scan_lines_per_block, s->ymax - line + 1);
368     uncompressed_size = s->scan_line_size * s->ysize;
369     if ((s->compr == EXR_RAW && (data_size != uncompressed_size ||
370                                  line_offset > buf_size - uncompressed_size)) ||
371         (s->compr != EXR_RAW && (data_size > uncompressed_size ||
372                                  line_offset > buf_size - data_size))) {
373         return AVERROR_INVALIDDATA;
374     }
375
376     if (data_size < uncompressed_size) {
377         av_fast_padded_malloc(&td->uncompressed_data, &td->uncompressed_size, uncompressed_size);
378         av_fast_padded_malloc(&td->tmp, &td->tmp_size, uncompressed_size);
379         if (!td->uncompressed_data || !td->tmp)
380             return AVERROR(ENOMEM);
381
382         switch (s->compr) {
383         case EXR_ZIP1:
384         case EXR_ZIP16:
385             ret = zip_uncompress(src, data_size, uncompressed_size, td);
386             break;
387         case EXR_PXR24:
388             ret = pxr24_uncompress(s, src, data_size, uncompressed_size, td);
389             break;
390         case EXR_RLE:
391             ret = rle_uncompress(src, data_size, uncompressed_size, td);
392         }
393
394         src = td->uncompressed_data;
395     }
396
397     channel_buffer[0] = src + xdelta * s->channel_offsets[0];
398     channel_buffer[1] = src + xdelta * s->channel_offsets[1];
399     channel_buffer[2] = src + xdelta * s->channel_offsets[2];
400     if (s->channel_offsets[3] >= 0)
401         channel_buffer[3] = src + xdelta * s->channel_offsets[3];
402
403     ptr = p->data[0] + line * p->linesize[0];
404     for (i = 0; i < s->scan_lines_per_block && line + i <= s->ymax; i++, ptr += p->linesize[0]) {
405         const uint8_t *r, *g, *b, *a;
406
407         r = channel_buffer[0];
408         g = channel_buffer[1];
409         b = channel_buffer[2];
410         if (channel_buffer[3])
411             a = channel_buffer[3];
412
413         ptr_x = (uint16_t *)ptr;
414
415         // Zero out the start if xmin is not 0
416         memset(ptr_x, 0, bxmin);
417         ptr_x += s->xmin * s->desc->nb_components;
418         if (s->pixel_type == EXR_FLOAT) {
419             // 32-bit
420             for (x = 0; x < xdelta; x++) {
421                 *ptr_x++ = exr_flt2uint(bytestream_get_le32(&r));
422                 *ptr_x++ = exr_flt2uint(bytestream_get_le32(&g));
423                 *ptr_x++ = exr_flt2uint(bytestream_get_le32(&b));
424                 if (channel_buffer[3])
425                     *ptr_x++ = exr_flt2uint(bytestream_get_le32(&a));
426             }
427         } else {
428             // 16-bit
429             for (x = 0; x < xdelta; x++) {
430                 *ptr_x++ = exr_halflt2uint(bytestream_get_le16(&r));
431                 *ptr_x++ = exr_halflt2uint(bytestream_get_le16(&g));
432                 *ptr_x++ = exr_halflt2uint(bytestream_get_le16(&b));
433                 if (channel_buffer[3])
434                     *ptr_x++ = exr_halflt2uint(bytestream_get_le16(&a));
435             }
436         }
437
438         // Zero out the end if xmax+1 is not w
439         memset(ptr_x, 0, axmax);
440
441         channel_buffer[0] += s->scan_line_size;
442         channel_buffer[1] += s->scan_line_size;
443         channel_buffer[2] += s->scan_line_size;
444         if (channel_buffer[3])
445             channel_buffer[3] += s->scan_line_size;
446     }
447
448     return 0;
449 }
450
451 static int decode_frame(AVCodecContext *avctx,
452                         void *data,
453                         int *got_frame,
454                         AVPacket *avpkt)
455 {
456     const uint8_t *buf      = avpkt->data;
457     unsigned int   buf_size = avpkt->size;
458     const uint8_t *buf_end  = buf + buf_size;
459
460     EXRContext *const s = avctx->priv_data;
461     ThreadFrame frame = { .f = data };
462     AVFrame *picture  = data;
463     uint8_t *ptr;
464
465     int i, y, magic_number, version, flags, ret;
466     int w = 0;
467     int h = 0;
468
469     int out_line_size;
470     int scan_line_blocks;
471
472     unsigned int current_channel_offset = 0;
473
474     s->xmin = ~0;
475     s->xmax = ~0;
476     s->ymin = ~0;
477     s->ymax = ~0;
478     s->xdelta = ~0;
479     s->ydelta = ~0;
480     s->channel_offsets[0] = -1;
481     s->channel_offsets[1] = -1;
482     s->channel_offsets[2] = -1;
483     s->channel_offsets[3] = -1;
484     s->pixel_type = -1;
485     s->nb_channels = 0;
486     s->compr = -1;
487     s->buf = buf;
488     s->buf_size = buf_size;
489
490     if (buf_size < 10) {
491         av_log(avctx, AV_LOG_ERROR, "Too short header to parse\n");
492         return AVERROR_INVALIDDATA;
493     }
494
495     magic_number = bytestream_get_le32(&buf);
496     if (magic_number != 20000630) { // As per documentation of OpenEXR it's supposed to be int 20000630 little-endian
497         av_log(avctx, AV_LOG_ERROR, "Wrong magic number %d\n", magic_number);
498         return AVERROR_INVALIDDATA;
499     }
500
501     version = bytestream_get_byte(&buf);
502     if (version != 2) {
503         av_log(avctx, AV_LOG_ERROR, "Unsupported version %d\n", version);
504         return AVERROR_PATCHWELCOME;
505     }
506
507     flags = bytestream_get_le24(&buf);
508     if (flags & 0x2) {
509         av_log(avctx, AV_LOG_ERROR, "Tile based images are not supported\n");
510         return AVERROR_PATCHWELCOME;
511     }
512
513     // Parse the header
514     while (buf < buf_end && buf[0]) {
515         unsigned int variable_buffer_data_size;
516         // Process the channel list
517         if (check_header_variable(avctx, &buf, buf_end, "channels", "chlist", 38, &variable_buffer_data_size) >= 0) {
518             const uint8_t *channel_list_end;
519             if (!variable_buffer_data_size)
520                 return AVERROR_INVALIDDATA;
521
522             channel_list_end = buf + variable_buffer_data_size;
523             while (channel_list_end - buf >= 19) {
524                 EXRChannel *channel;
525                 int current_pixel_type = -1;
526                 int channel_index = -1;
527                 int xsub, ysub;
528
529                 if (!strcmp(buf, "R"))
530                     channel_index = 0;
531                 else if (!strcmp(buf, "G"))
532                     channel_index = 1;
533                 else if (!strcmp(buf, "B"))
534                     channel_index = 2;
535                 else if (!strcmp(buf, "A"))
536                     channel_index = 3;
537                 else
538                     av_log(avctx, AV_LOG_WARNING, "Unsupported channel %.256s\n", buf);
539
540                 while (bytestream_get_byte(&buf) && buf < channel_list_end)
541                     continue; /* skip */
542
543                 if (channel_list_end - * &buf < 4) {
544                     av_log(avctx, AV_LOG_ERROR, "Incomplete header\n");
545                     return AVERROR_INVALIDDATA;
546                 }
547
548                 current_pixel_type = bytestream_get_le32(&buf);
549                 if (current_pixel_type > 2) {
550                     av_log(avctx, AV_LOG_ERROR, "Unknown pixel type\n");
551                     return AVERROR_INVALIDDATA;
552                 }
553
554                 buf += 4;
555                 xsub = bytestream_get_le32(&buf);
556                 ysub = bytestream_get_le32(&buf);
557                 if (xsub != 1 || ysub != 1) {
558                     av_log(avctx, AV_LOG_ERROR, "Unsupported subsampling %dx%d\n", xsub, ysub);
559                     return AVERROR_PATCHWELCOME;
560                 }
561
562                 if (channel_index >= 0) {
563                     if (s->pixel_type != -1 && s->pixel_type != current_pixel_type) {
564                         av_log(avctx, AV_LOG_ERROR, "RGB channels not of the same depth\n");
565                         return AVERROR_INVALIDDATA;
566                     }
567                     s->pixel_type = current_pixel_type;
568                     s->channel_offsets[channel_index] = current_channel_offset;
569                 }
570
571                 s->channels = av_realloc_f(s->channels, ++s->nb_channels, sizeof(EXRChannel));
572                 if (!s->channels)
573                     return AVERROR(ENOMEM);
574                 channel = &s->channels[s->nb_channels - 1];
575                 channel->pixel_type = current_pixel_type;
576                 channel->xsub = xsub;
577                 channel->ysub = ysub;
578
579                 current_channel_offset += 1 << current_pixel_type;
580             }
581
582             /* Check if all channels are set with an offset or if the channels
583              * are causing an overflow  */
584
585             if (FFMIN3(s->channel_offsets[0],
586                        s->channel_offsets[1],
587                        s->channel_offsets[2]) < 0) {
588                 if (s->channel_offsets[0] < 0)
589                     av_log(avctx, AV_LOG_ERROR, "Missing red channel\n");
590                 if (s->channel_offsets[1] < 0)
591                     av_log(avctx, AV_LOG_ERROR, "Missing green channel\n");
592                 if (s->channel_offsets[2] < 0)
593                     av_log(avctx, AV_LOG_ERROR, "Missing blue channel\n");
594                 return AVERROR_INVALIDDATA;
595             }
596
597             buf = channel_list_end;
598             continue;
599         } else if (check_header_variable(avctx, &buf, buf_end, "dataWindow", "box2i", 31, &variable_buffer_data_size) >= 0) {
600             if (!variable_buffer_data_size)
601                 return AVERROR_INVALIDDATA;
602
603             s->xmin = AV_RL32(buf);
604             s->ymin = AV_RL32(buf + 4);
605             s->xmax = AV_RL32(buf + 8);
606             s->ymax = AV_RL32(buf + 12);
607             s->xdelta = (s->xmax - s->xmin) + 1;
608             s->ydelta = (s->ymax - s->ymin) + 1;
609
610             buf += variable_buffer_data_size;
611             continue;
612         } else if (check_header_variable(avctx, &buf, buf_end, "displayWindow", "box2i", 34, &variable_buffer_data_size) >= 0) {
613             if (!variable_buffer_data_size)
614                 return AVERROR_INVALIDDATA;
615
616             w = AV_RL32(buf + 8) + 1;
617             h = AV_RL32(buf + 12) + 1;
618
619             buf += variable_buffer_data_size;
620             continue;
621         } else if (check_header_variable(avctx, &buf, buf_end, "lineOrder", "lineOrder", 25, &variable_buffer_data_size) >= 0) {
622             if (!variable_buffer_data_size)
623                 return AVERROR_INVALIDDATA;
624
625             av_log(avctx, AV_LOG_DEBUG, "line order : %d\n", *buf);
626             if (*buf > 2) {
627                 av_log(avctx, AV_LOG_ERROR, "Unknown line order\n");
628                 return AVERROR_INVALIDDATA;
629             }
630
631             buf += variable_buffer_data_size;
632             continue;
633         } else if (check_header_variable(avctx, &buf, buf_end, "pixelAspectRatio", "float", 31, &variable_buffer_data_size) >= 0) {
634             if (!variable_buffer_data_size)
635                 return AVERROR_INVALIDDATA;
636
637             avctx->sample_aspect_ratio = av_d2q(av_int2float(AV_RL32(buf)), 255);
638
639             buf += variable_buffer_data_size;
640             continue;
641         } else if (check_header_variable(avctx, &buf, buf_end, "compression", "compression", 29, &variable_buffer_data_size) >= 0) {
642             if (!variable_buffer_data_size)
643                 return AVERROR_INVALIDDATA;
644
645             if (s->compr == -1)
646                 s->compr = *buf;
647             else
648                 av_log(avctx, AV_LOG_WARNING, "Found more than one compression attribute\n");
649
650             buf += variable_buffer_data_size;
651             continue;
652         }
653
654         // Check if there is enough bytes for a header
655         if (buf_end - buf <= 9) {
656             av_log(avctx, AV_LOG_ERROR, "Incomplete header\n");
657             return AVERROR_INVALIDDATA;
658         }
659
660         // Process unknown variables
661         for (i = 0; i < 2; i++) {
662             // Skip variable name/type
663             while (++buf < buf_end)
664                 if (buf[0] == 0x0)
665                     break;
666         }
667         buf++;
668         // Skip variable length
669         if (buf_end - buf >= 5) {
670             variable_buffer_data_size = get_header_variable_length(&buf, buf_end);
671             if (!variable_buffer_data_size) {
672                 av_log(avctx, AV_LOG_ERROR, "Incomplete header\n");
673                 return AVERROR_INVALIDDATA;
674             }
675             buf += variable_buffer_data_size;
676         }
677     }
678
679     if (s->compr == -1) {
680         av_log(avctx, AV_LOG_ERROR, "Missing compression attribute\n");
681         return AVERROR_INVALIDDATA;
682     }
683
684     if (buf >= buf_end) {
685         av_log(avctx, AV_LOG_ERROR, "Incomplete frame\n");
686         return AVERROR_INVALIDDATA;
687     }
688     buf++;
689
690     switch (s->pixel_type) {
691     case EXR_FLOAT:
692     case EXR_HALF:
693         if (s->channel_offsets[3] >= 0)
694             avctx->pix_fmt = AV_PIX_FMT_RGBA64;
695         else
696             avctx->pix_fmt = AV_PIX_FMT_RGB48;
697         break;
698     case EXR_UINT:
699         av_log_missing_feature(avctx, "32-bit unsigned int", 1);
700         return AVERROR_PATCHWELCOME;
701     default:
702         av_log(avctx, AV_LOG_ERROR, "Missing channel list\n");
703         return AVERROR_INVALIDDATA;
704     }
705
706     switch (s->compr) {
707     case EXR_RAW:
708     case EXR_RLE:
709     case EXR_ZIP1:
710         s->scan_lines_per_block = 1;
711         break;
712     case EXR_PXR24:
713     case EXR_ZIP16:
714         s->scan_lines_per_block = 16;
715         break;
716     default:
717         av_log(avctx, AV_LOG_ERROR, "Compression type %d is not supported\n", s->compr);
718         return AVERROR_PATCHWELCOME;
719     }
720
721     if (av_image_check_size(w, h, 0, avctx))
722         return AVERROR_INVALIDDATA;
723
724     // Verify the xmin, xmax, ymin, ymax and xdelta before setting the actual image size
725     if (s->xmin > s->xmax ||
726         s->ymin > s->ymax ||
727         s->xdelta != s->xmax - s->xmin + 1 ||
728         s->xmax >= w || s->ymax >= h) {
729         av_log(avctx, AV_LOG_ERROR, "Wrong sizing or missing size information\n");
730         return AVERROR_INVALIDDATA;
731     }
732
733     if (w != avctx->width || h != avctx->height) {
734         avcodec_set_dimensions(avctx, w, h);
735     }
736
737     s->desc = av_pix_fmt_desc_get(avctx->pix_fmt);
738     out_line_size = avctx->width * 2 * s->desc->nb_components;
739     s->scan_line_size = s->xdelta * current_channel_offset;
740     scan_line_blocks = (s->ydelta + s->scan_lines_per_block - 1) / s->scan_lines_per_block;
741
742     if (s->compr != EXR_RAW) {
743         size_t thread_data_size, prev_size;
744         EXRThreadData *m;
745
746         prev_size = s->thread_data_size;
747         if (av_size_mult(avctx->thread_count, sizeof(EXRThreadData), &thread_data_size))
748             return AVERROR(EINVAL);
749
750         m = av_fast_realloc(s->thread_data, &s->thread_data_size, thread_data_size);
751         if (!m)
752             return AVERROR(ENOMEM);
753         s->thread_data = m;
754         memset(s->thread_data + prev_size, 0, s->thread_data_size - prev_size);
755     }
756
757     if ((ret = ff_thread_get_buffer(avctx, &frame, 0)) < 0) {
758         av_log(avctx, AV_LOG_ERROR, "get_buffer() failed\n");
759         return ret;
760     }
761
762     if (buf_end - buf < scan_line_blocks * 8)
763         return AVERROR_INVALIDDATA;
764     s->table = buf;
765     ptr = picture->data[0];
766
767     // Zero out the start if ymin is not 0
768     for (y = 0; y < s->ymin; y++) {
769         memset(ptr, 0, out_line_size);
770         ptr += picture->linesize[0];
771     }
772
773     s->picture = picture;
774     avctx->execute2(avctx, decode_block, s->thread_data, NULL, scan_line_blocks);
775
776     // Zero out the end if ymax+1 is not h
777     for (y = s->ymax + 1; y < avctx->height; y++) {
778         memset(ptr, 0, out_line_size);
779         ptr += picture->linesize[0];
780     }
781
782     *got_frame = 1;
783
784     return buf_size;
785 }
786
787 static av_cold int decode_end(AVCodecContext *avctx)
788 {
789     EXRContext *s = avctx->priv_data;
790     int i;
791
792     for (i = 0; i < s->thread_data_size / sizeof(EXRThreadData); i++) {
793         EXRThreadData *td = &s->thread_data[i];
794         av_free(td->uncompressed_data);
795         av_free(td->tmp);
796     }
797
798     av_freep(&s->thread_data);
799     s->thread_data_size = 0;
800     av_freep(&s->channels);
801
802     return 0;
803 }
804
805 AVCodec ff_exr_decoder = {
806     .name               = "exr",
807     .type               = AVMEDIA_TYPE_VIDEO,
808     .id                 = AV_CODEC_ID_EXR,
809     .priv_data_size     = sizeof(EXRContext),
810     .close              = decode_end,
811     .decode             = decode_frame,
812     .capabilities       = CODEC_CAP_DR1 | CODEC_CAP_FRAME_THREADS | CODEC_CAP_SLICE_THREADS,
813     .long_name          = NULL_IF_CONFIG_SMALL("OpenEXR image"),
814 };