dsputil: move a bink-only function to binkdsp
[ffmpeg.git] / libavcodec / bink.c
1 /*
2  * Bink video decoder
3  * Copyright (c) 2009 Konstantin Shishkov
4  * Copyright (C) 2011 Peter Ross <pross@xvid.org>
5  *
6  * This file is part of Libav.
7  *
8  * Libav is free software; you can redistribute it and/or
9  * modify it under the terms of the GNU Lesser General Public
10  * License as published by the Free Software Foundation; either
11  * version 2.1 of the License, or (at your option) any later version.
12  *
13  * Libav is distributed in the hope that it will be useful,
14  * but WITHOUT ANY WARRANTY; without even the implied warranty of
15  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
16  * Lesser General Public License for more details.
17  *
18  * You should have received a copy of the GNU Lesser General Public
19  * License along with Libav; if not, write to the Free Software
20  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
21  */
22
23 #include "libavutil/imgutils.h"
24 #include "avcodec.h"
25 #include "dsputil.h"
26 #include "binkdata.h"
27 #include "binkdsp.h"
28 #include "mathops.h"
29
30 #define ALT_BITSTREAM_READER_LE
31 #include "get_bits.h"
32
33 #define BINK_FLAG_ALPHA 0x00100000
34 #define BINK_FLAG_GRAY  0x00020000
35
36 static VLC bink_trees[16];
37
38 /**
39  * IDs for different data types used in old version of Bink video codec
40  */
41 enum OldSources {
42     BINKB_SRC_BLOCK_TYPES = 0, ///< 8x8 block types
43     BINKB_SRC_COLORS,          ///< pixel values used for different block types
44     BINKB_SRC_PATTERN,         ///< 8-bit values for 2-colour pattern fill
45     BINKB_SRC_X_OFF,           ///< X components of motion value
46     BINKB_SRC_Y_OFF,           ///< Y components of motion value
47     BINKB_SRC_INTRA_DC,        ///< DC values for intrablocks with DCT
48     BINKB_SRC_INTER_DC,        ///< DC values for interblocks with DCT
49     BINKB_SRC_INTRA_Q,         ///< quantizer values for intrablocks with DCT
50     BINKB_SRC_INTER_Q,         ///< quantizer values for interblocks with DCT
51     BINKB_SRC_INTER_COEFS,     ///< number of coefficients for residue blocks
52
53     BINKB_NB_SRC
54 };
55
56 static const int binkb_bundle_sizes[BINKB_NB_SRC] = {
57     4, 8, 8, 5, 5, 11, 11, 4, 4, 7
58 };
59
60 static const int binkb_bundle_signed[BINKB_NB_SRC] = {
61     0, 0, 0, 1, 1, 0, 1, 0, 0, 0
62 };
63
64 static int32_t binkb_intra_quant[16][64];
65 static int32_t binkb_inter_quant[16][64];
66
67 /**
68  * IDs for different data types used in Bink video codec
69  */
70 enum Sources {
71     BINK_SRC_BLOCK_TYPES = 0, ///< 8x8 block types
72     BINK_SRC_SUB_BLOCK_TYPES, ///< 16x16 block types (a subset of 8x8 block types)
73     BINK_SRC_COLORS,          ///< pixel values used for different block types
74     BINK_SRC_PATTERN,         ///< 8-bit values for 2-colour pattern fill
75     BINK_SRC_X_OFF,           ///< X components of motion value
76     BINK_SRC_Y_OFF,           ///< Y components of motion value
77     BINK_SRC_INTRA_DC,        ///< DC values for intrablocks with DCT
78     BINK_SRC_INTER_DC,        ///< DC values for interblocks with DCT
79     BINK_SRC_RUN,             ///< run lengths for special fill block
80
81     BINK_NB_SRC
82 };
83
84 /**
85  * data needed to decode 4-bit Huffman-coded value
86  */
87 typedef struct Tree {
88     int     vlc_num;  ///< tree number (in bink_trees[])
89     uint8_t syms[16]; ///< leaf value to symbol mapping
90 } Tree;
91
92 #define GET_HUFF(gb, tree)  (tree).syms[get_vlc2(gb, bink_trees[(tree).vlc_num].table,\
93                                                  bink_trees[(tree).vlc_num].bits, 1)]
94
95 /**
96  * data structure used for decoding single Bink data type
97  */
98 typedef struct Bundle {
99     int     len;       ///< length of number of entries to decode (in bits)
100     Tree    tree;      ///< Huffman tree-related data
101     uint8_t *data;     ///< buffer for decoded symbols
102     uint8_t *data_end; ///< buffer end
103     uint8_t *cur_dec;  ///< pointer to the not yet decoded part of the buffer
104     uint8_t *cur_ptr;  ///< pointer to the data that is not read from buffer yet
105 } Bundle;
106
107 /*
108  * Decoder context
109  */
110 typedef struct BinkContext {
111     AVCodecContext *avctx;
112     DSPContext     dsp;
113     BinkDSPContext bdsp;
114     AVFrame        pic, last;
115     int            version;              ///< internal Bink file version
116     int            has_alpha;
117     int            swap_planes;
118
119     Bundle         bundle[BINKB_NB_SRC]; ///< bundles for decoding all data types
120     Tree           col_high[16];         ///< trees for decoding high nibble in "colours" data type
121     int            col_lastval;          ///< value of last decoded high nibble in "colours" data type
122 } BinkContext;
123
124 /**
125  * Bink video block types
126  */
127 enum BlockTypes {
128     SKIP_BLOCK = 0, ///< skipped block
129     SCALED_BLOCK,   ///< block has size 16x16
130     MOTION_BLOCK,   ///< block is copied from previous frame with some offset
131     RUN_BLOCK,      ///< block is composed from runs of colours with custom scan order
132     RESIDUE_BLOCK,  ///< motion block with some difference added
133     INTRA_BLOCK,    ///< intra DCT block
134     FILL_BLOCK,     ///< block is filled with single colour
135     INTER_BLOCK,    ///< motion block with DCT applied to the difference
136     PATTERN_BLOCK,  ///< block is filled with two colours following custom pattern
137     RAW_BLOCK,      ///< uncoded 8x8 block
138 };
139
140 /**
141  * Initialize length length in all bundles.
142  *
143  * @param c     decoder context
144  * @param width plane width
145  * @param bw    plane width in 8x8 blocks
146  */
147 static void init_lengths(BinkContext *c, int width, int bw)
148 {
149     c->bundle[BINK_SRC_BLOCK_TYPES].len = av_log2((width >> 3) + 511) + 1;
150
151     c->bundle[BINK_SRC_SUB_BLOCK_TYPES].len = av_log2((width >> 4) + 511) + 1;
152
153     c->bundle[BINK_SRC_COLORS].len = av_log2(bw*64 + 511) + 1;
154
155     c->bundle[BINK_SRC_INTRA_DC].len =
156     c->bundle[BINK_SRC_INTER_DC].len =
157     c->bundle[BINK_SRC_X_OFF].len =
158     c->bundle[BINK_SRC_Y_OFF].len = av_log2((width >> 3) + 511) + 1;
159
160     c->bundle[BINK_SRC_PATTERN].len = av_log2((bw << 3) + 511) + 1;
161
162     c->bundle[BINK_SRC_RUN].len = av_log2(bw*48 + 511) + 1;
163 }
164
165 /**
166  * Allocate memory for bundles.
167  *
168  * @param c decoder context
169  */
170 static av_cold void init_bundles(BinkContext *c)
171 {
172     int bw, bh, blocks;
173     int i;
174
175     bw = (c->avctx->width  + 7) >> 3;
176     bh = (c->avctx->height + 7) >> 3;
177     blocks = bw * bh;
178
179     for (i = 0; i < BINKB_NB_SRC; i++) {
180         c->bundle[i].data = av_malloc(blocks * 64);
181         c->bundle[i].data_end = c->bundle[i].data + blocks * 64;
182     }
183 }
184
185 /**
186  * Free memory used by bundles.
187  *
188  * @param c decoder context
189  */
190 static av_cold void free_bundles(BinkContext *c)
191 {
192     int i;
193     for (i = 0; i < BINKB_NB_SRC; i++)
194         av_freep(&c->bundle[i].data);
195 }
196
197 /**
198  * Merge two consequent lists of equal size depending on bits read.
199  *
200  * @param gb   context for reading bits
201  * @param dst  buffer where merged list will be written to
202  * @param src  pointer to the head of the first list (the second lists starts at src+size)
203  * @param size input lists size
204  */
205 static void merge(GetBitContext *gb, uint8_t *dst, uint8_t *src, int size)
206 {
207     uint8_t *src2 = src + size;
208     int size2 = size;
209
210     do {
211         if (!get_bits1(gb)) {
212             *dst++ = *src++;
213             size--;
214         } else {
215             *dst++ = *src2++;
216             size2--;
217         }
218     } while (size && size2);
219
220     while (size--)
221         *dst++ = *src++;
222     while (size2--)
223         *dst++ = *src2++;
224 }
225
226 /**
227  * Read information about Huffman tree used to decode data.
228  *
229  * @param gb   context for reading bits
230  * @param tree pointer for storing tree data
231  */
232 static void read_tree(GetBitContext *gb, Tree *tree)
233 {
234     uint8_t tmp1[16], tmp2[16], *in = tmp1, *out = tmp2;
235     int i, t, len;
236
237     tree->vlc_num = get_bits(gb, 4);
238     if (!tree->vlc_num) {
239         for (i = 0; i < 16; i++)
240             tree->syms[i] = i;
241         return;
242     }
243     if (get_bits1(gb)) {
244         len = get_bits(gb, 3);
245         memset(tmp1, 0, sizeof(tmp1));
246         for (i = 0; i <= len; i++) {
247             tree->syms[i] = get_bits(gb, 4);
248             tmp1[tree->syms[i]] = 1;
249         }
250         for (i = 0; i < 16; i++)
251             if (!tmp1[i])
252                 tree->syms[++len] = i;
253     } else {
254         len = get_bits(gb, 2);
255         for (i = 0; i < 16; i++)
256             in[i] = i;
257         for (i = 0; i <= len; i++) {
258             int size = 1 << i;
259             for (t = 0; t < 16; t += size << 1)
260                 merge(gb, out + t, in + t, size);
261             FFSWAP(uint8_t*, in, out);
262         }
263         memcpy(tree->syms, in, 16);
264     }
265 }
266
267 /**
268  * Prepare bundle for decoding data.
269  *
270  * @param gb          context for reading bits
271  * @param c           decoder context
272  * @param bundle_num  number of the bundle to initialize
273  */
274 static void read_bundle(GetBitContext *gb, BinkContext *c, int bundle_num)
275 {
276     int i;
277
278     if (bundle_num == BINK_SRC_COLORS) {
279         for (i = 0; i < 16; i++)
280             read_tree(gb, &c->col_high[i]);
281         c->col_lastval = 0;
282     }
283     if (bundle_num != BINK_SRC_INTRA_DC && bundle_num != BINK_SRC_INTER_DC)
284         read_tree(gb, &c->bundle[bundle_num].tree);
285     c->bundle[bundle_num].cur_dec =
286     c->bundle[bundle_num].cur_ptr = c->bundle[bundle_num].data;
287 }
288
289 /**
290  * common check before starting decoding bundle data
291  *
292  * @param gb context for reading bits
293  * @param b  bundle
294  * @param t  variable where number of elements to decode will be stored
295  */
296 #define CHECK_READ_VAL(gb, b, t) \
297     if (!b->cur_dec || (b->cur_dec > b->cur_ptr)) \
298         return 0; \
299     t = get_bits(gb, b->len); \
300     if (!t) { \
301         b->cur_dec = NULL; \
302         return 0; \
303     } \
304
305 static int read_runs(AVCodecContext *avctx, GetBitContext *gb, Bundle *b)
306 {
307     int t, v;
308     const uint8_t *dec_end;
309
310     CHECK_READ_VAL(gb, b, t);
311     dec_end = b->cur_dec + t;
312     if (dec_end > b->data_end) {
313         av_log(avctx, AV_LOG_ERROR, "Run value went out of bounds\n");
314         return -1;
315     }
316     if (get_bits1(gb)) {
317         v = get_bits(gb, 4);
318         memset(b->cur_dec, v, t);
319         b->cur_dec += t;
320     } else {
321         while (b->cur_dec < dec_end)
322             *b->cur_dec++ = GET_HUFF(gb, b->tree);
323     }
324     return 0;
325 }
326
327 static int read_motion_values(AVCodecContext *avctx, GetBitContext *gb, Bundle *b)
328 {
329     int t, sign, v;
330     const uint8_t *dec_end;
331
332     CHECK_READ_VAL(gb, b, t);
333     dec_end = b->cur_dec + t;
334     if (dec_end > b->data_end) {
335         av_log(avctx, AV_LOG_ERROR, "Too many motion values\n");
336         return -1;
337     }
338     if (get_bits1(gb)) {
339         v = get_bits(gb, 4);
340         if (v) {
341             sign = -get_bits1(gb);
342             v = (v ^ sign) - sign;
343         }
344         memset(b->cur_dec, v, t);
345         b->cur_dec += t;
346     } else {
347         do {
348             v = GET_HUFF(gb, b->tree);
349             if (v) {
350                 sign = -get_bits1(gb);
351                 v = (v ^ sign) - sign;
352             }
353             *b->cur_dec++ = v;
354         } while (b->cur_dec < dec_end);
355     }
356     return 0;
357 }
358
359 static const uint8_t bink_rlelens[4] = { 4, 8, 12, 32 };
360
361 static int read_block_types(AVCodecContext *avctx, GetBitContext *gb, Bundle *b)
362 {
363     int t, v;
364     int last = 0;
365     const uint8_t *dec_end;
366
367     CHECK_READ_VAL(gb, b, t);
368     dec_end = b->cur_dec + t;
369     if (dec_end > b->data_end) {
370         av_log(avctx, AV_LOG_ERROR, "Too many block type values\n");
371         return -1;
372     }
373     if (get_bits1(gb)) {
374         v = get_bits(gb, 4);
375         memset(b->cur_dec, v, t);
376         b->cur_dec += t;
377     } else {
378         do {
379             v = GET_HUFF(gb, b->tree);
380             if (v < 12) {
381                 last = v;
382                 *b->cur_dec++ = v;
383             } else {
384                 int run = bink_rlelens[v - 12];
385
386                 memset(b->cur_dec, last, run);
387                 b->cur_dec += run;
388             }
389         } while (b->cur_dec < dec_end);
390     }
391     return 0;
392 }
393
394 static int read_patterns(AVCodecContext *avctx, GetBitContext *gb, Bundle *b)
395 {
396     int t, v;
397     const uint8_t *dec_end;
398
399     CHECK_READ_VAL(gb, b, t);
400     dec_end = b->cur_dec + t;
401     if (dec_end > b->data_end) {
402         av_log(avctx, AV_LOG_ERROR, "Too many pattern values\n");
403         return -1;
404     }
405     while (b->cur_dec < dec_end) {
406         v  = GET_HUFF(gb, b->tree);
407         v |= GET_HUFF(gb, b->tree) << 4;
408         *b->cur_dec++ = v;
409     }
410
411     return 0;
412 }
413
414 static int read_colors(GetBitContext *gb, Bundle *b, BinkContext *c)
415 {
416     int t, sign, v;
417     const uint8_t *dec_end;
418
419     CHECK_READ_VAL(gb, b, t);
420     dec_end = b->cur_dec + t;
421     if (dec_end > b->data_end) {
422         av_log(c->avctx, AV_LOG_ERROR, "Too many color values\n");
423         return -1;
424     }
425     if (get_bits1(gb)) {
426         c->col_lastval = GET_HUFF(gb, c->col_high[c->col_lastval]);
427         v = GET_HUFF(gb, b->tree);
428         v = (c->col_lastval << 4) | v;
429         if (c->version < 'i') {
430             sign = ((int8_t) v) >> 7;
431             v = ((v & 0x7F) ^ sign) - sign;
432             v += 0x80;
433         }
434         memset(b->cur_dec, v, t);
435         b->cur_dec += t;
436     } else {
437         while (b->cur_dec < dec_end) {
438             c->col_lastval = GET_HUFF(gb, c->col_high[c->col_lastval]);
439             v = GET_HUFF(gb, b->tree);
440             v = (c->col_lastval << 4) | v;
441             if (c->version < 'i') {
442                 sign = ((int8_t) v) >> 7;
443                 v = ((v & 0x7F) ^ sign) - sign;
444                 v += 0x80;
445             }
446             *b->cur_dec++ = v;
447         }
448     }
449     return 0;
450 }
451
452 /** number of bits used to store first DC value in bundle */
453 #define DC_START_BITS 11
454
455 static int read_dcs(AVCodecContext *avctx, GetBitContext *gb, Bundle *b,
456                     int start_bits, int has_sign)
457 {
458     int i, j, len, len2, bsize, sign, v, v2;
459     int16_t *dst = (int16_t*)b->cur_dec;
460
461     CHECK_READ_VAL(gb, b, len);
462     v = get_bits(gb, start_bits - has_sign);
463     if (v && has_sign) {
464         sign = -get_bits1(gb);
465         v = (v ^ sign) - sign;
466     }
467     *dst++ = v;
468     len--;
469     for (i = 0; i < len; i += 8) {
470         len2 = FFMIN(len - i, 8);
471         bsize = get_bits(gb, 4);
472         if (bsize) {
473             for (j = 0; j < len2; j++) {
474                 v2 = get_bits(gb, bsize);
475                 if (v2) {
476                     sign = -get_bits1(gb);
477                     v2 = (v2 ^ sign) - sign;
478                 }
479                 v += v2;
480                 *dst++ = v;
481                 if (v < -32768 || v > 32767) {
482                     av_log(avctx, AV_LOG_ERROR, "DC value went out of bounds: %d\n", v);
483                     return -1;
484                 }
485             }
486         } else {
487             for (j = 0; j < len2; j++)
488                 *dst++ = v;
489         }
490     }
491
492     b->cur_dec = (uint8_t*)dst;
493     return 0;
494 }
495
496 /**
497  * Retrieve next value from bundle.
498  *
499  * @param c      decoder context
500  * @param bundle bundle number
501  */
502 static inline int get_value(BinkContext *c, int bundle)
503 {
504     int ret;
505
506     if (bundle < BINK_SRC_X_OFF || bundle == BINK_SRC_RUN)
507         return *c->bundle[bundle].cur_ptr++;
508     if (bundle == BINK_SRC_X_OFF || bundle == BINK_SRC_Y_OFF)
509         return (int8_t)*c->bundle[bundle].cur_ptr++;
510     ret = *(int16_t*)c->bundle[bundle].cur_ptr;
511     c->bundle[bundle].cur_ptr += 2;
512     return ret;
513 }
514
515 static void binkb_init_bundle(BinkContext *c, int bundle_num)
516 {
517     c->bundle[bundle_num].cur_dec =
518     c->bundle[bundle_num].cur_ptr = c->bundle[bundle_num].data;
519     c->bundle[bundle_num].len = 13;
520 }
521
522 static void binkb_init_bundles(BinkContext *c)
523 {
524     int i;
525     for (i = 0; i < BINKB_NB_SRC; i++)
526         binkb_init_bundle(c, i);
527 }
528
529 static int binkb_read_bundle(BinkContext *c, GetBitContext *gb, int bundle_num)
530 {
531     const int bits = binkb_bundle_sizes[bundle_num];
532     const int mask = 1 << (bits - 1);
533     const int issigned = binkb_bundle_signed[bundle_num];
534     Bundle *b = &c->bundle[bundle_num];
535     int i, len;
536
537     CHECK_READ_VAL(gb, b, len);
538     if (bits <= 8) {
539         if (!issigned) {
540             for (i = 0; i < len; i++)
541                 *b->cur_dec++ = get_bits(gb, bits);
542         } else {
543             for (i = 0; i < len; i++)
544                 *b->cur_dec++ = get_bits(gb, bits) - mask;
545         }
546     } else {
547         int16_t *dst = (int16_t*)b->cur_dec;
548
549         if (!issigned) {
550             for (i = 0; i < len; i++)
551                 *dst++ = get_bits(gb, bits);
552         } else {
553             for (i = 0; i < len; i++)
554                 *dst++ = get_bits(gb, bits) - mask;
555         }
556         b->cur_dec = (uint8_t*)dst;
557     }
558     return 0;
559 }
560
561 static inline int binkb_get_value(BinkContext *c, int bundle_num)
562 {
563     int16_t ret;
564     const int bits = binkb_bundle_sizes[bundle_num];
565
566     if (bits <= 8) {
567         int val = *c->bundle[bundle_num].cur_ptr++;
568         return binkb_bundle_signed[bundle_num] ? (int8_t)val : val;
569     }
570     ret = *(int16_t*)c->bundle[bundle_num].cur_ptr;
571     c->bundle[bundle_num].cur_ptr += 2;
572     return ret;
573 }
574
575 /**
576  * Read 8x8 block of DCT coefficients.
577  *
578  * @param gb       context for reading bits
579  * @param block    place for storing coefficients
580  * @param scan     scan order table
581  * @param quant_matrices quantization matrices
582  * @return 0 for success, negative value in other cases
583  */
584 static int read_dct_coeffs(GetBitContext *gb, int32_t block[64], const uint8_t *scan,
585                            const int32_t quant_matrices[16][64], int q)
586 {
587     int coef_list[128];
588     int mode_list[128];
589     int i, t, mask, bits, ccoef, mode, sign;
590     int list_start = 64, list_end = 64, list_pos;
591     int coef_count = 0;
592     int coef_idx[64];
593     int quant_idx;
594     const int32_t *quant;
595
596     coef_list[list_end] = 4;  mode_list[list_end++] = 0;
597     coef_list[list_end] = 24; mode_list[list_end++] = 0;
598     coef_list[list_end] = 44; mode_list[list_end++] = 0;
599     coef_list[list_end] = 1;  mode_list[list_end++] = 3;
600     coef_list[list_end] = 2;  mode_list[list_end++] = 3;
601     coef_list[list_end] = 3;  mode_list[list_end++] = 3;
602
603     bits = get_bits(gb, 4) - 1;
604     for (mask = 1 << bits; bits >= 0; mask >>= 1, bits--) {
605         list_pos = list_start;
606         while (list_pos < list_end) {
607             if (!(mode_list[list_pos] | coef_list[list_pos]) || !get_bits1(gb)) {
608                 list_pos++;
609                 continue;
610             }
611             ccoef = coef_list[list_pos];
612             mode  = mode_list[list_pos];
613             switch (mode) {
614             case 0:
615                 coef_list[list_pos] = ccoef + 4;
616                 mode_list[list_pos] = 1;
617             case 2:
618                 if (mode == 2) {
619                     coef_list[list_pos]   = 0;
620                     mode_list[list_pos++] = 0;
621                 }
622                 for (i = 0; i < 4; i++, ccoef++) {
623                     if (get_bits1(gb)) {
624                         coef_list[--list_start] = ccoef;
625                         mode_list[  list_start] = 3;
626                     } else {
627                         int t;
628                         if (!bits) {
629                             t = 1 - (get_bits1(gb) << 1);
630                         } else {
631                             t = get_bits(gb, bits) | mask;
632                             sign = -get_bits1(gb);
633                             t = (t ^ sign) - sign;
634                         }
635                         block[scan[ccoef]] = t;
636                         coef_idx[coef_count++] = ccoef;
637                     }
638                 }
639                 break;
640             case 1:
641                 mode_list[list_pos] = 2;
642                 for (i = 0; i < 3; i++) {
643                     ccoef += 4;
644                     coef_list[list_end]   = ccoef;
645                     mode_list[list_end++] = 2;
646                 }
647                 break;
648             case 3:
649                 if (!bits) {
650                     t = 1 - (get_bits1(gb) << 1);
651                 } else {
652                     t = get_bits(gb, bits) | mask;
653                     sign = -get_bits1(gb);
654                     t = (t ^ sign) - sign;
655                 }
656                 block[scan[ccoef]] = t;
657                 coef_idx[coef_count++] = ccoef;
658                 coef_list[list_pos]   = 0;
659                 mode_list[list_pos++] = 0;
660                 break;
661             }
662         }
663     }
664
665     if (q == -1) {
666         quant_idx = get_bits(gb, 4);
667     } else {
668         quant_idx = q;
669     }
670
671     quant = quant_matrices[quant_idx];
672
673     block[0] = (block[0] * quant[0]) >> 11;
674     for (i = 0; i < coef_count; i++) {
675         int idx = coef_idx[i];
676         block[scan[idx]] = (block[scan[idx]] * quant[idx]) >> 11;
677     }
678
679     return 0;
680 }
681
682 /**
683  * Read 8x8 block with residue after motion compensation.
684  *
685  * @param gb          context for reading bits
686  * @param block       place to store read data
687  * @param masks_count number of masks to decode
688  * @return 0 on success, negative value in other cases
689  */
690 static int read_residue(GetBitContext *gb, DCTELEM block[64], int masks_count)
691 {
692     int coef_list[128];
693     int mode_list[128];
694     int i, sign, mask, ccoef, mode;
695     int list_start = 64, list_end = 64, list_pos;
696     int nz_coeff[64];
697     int nz_coeff_count = 0;
698
699     coef_list[list_end] =  4; mode_list[list_end++] = 0;
700     coef_list[list_end] = 24; mode_list[list_end++] = 0;
701     coef_list[list_end] = 44; mode_list[list_end++] = 0;
702     coef_list[list_end] =  0; mode_list[list_end++] = 2;
703
704     for (mask = 1 << get_bits(gb, 3); mask; mask >>= 1) {
705         for (i = 0; i < nz_coeff_count; i++) {
706             if (!get_bits1(gb))
707                 continue;
708             if (block[nz_coeff[i]] < 0)
709                 block[nz_coeff[i]] -= mask;
710             else
711                 block[nz_coeff[i]] += mask;
712             masks_count--;
713             if (masks_count < 0)
714                 return 0;
715         }
716         list_pos = list_start;
717         while (list_pos < list_end) {
718             if (!(coef_list[list_pos] | mode_list[list_pos]) || !get_bits1(gb)) {
719                 list_pos++;
720                 continue;
721             }
722             ccoef = coef_list[list_pos];
723             mode  = mode_list[list_pos];
724             switch (mode) {
725             case 0:
726                 coef_list[list_pos] = ccoef + 4;
727                 mode_list[list_pos] = 1;
728             case 2:
729                 if (mode == 2) {
730                     coef_list[list_pos]   = 0;
731                     mode_list[list_pos++] = 0;
732                 }
733                 for (i = 0; i < 4; i++, ccoef++) {
734                     if (get_bits1(gb)) {
735                         coef_list[--list_start] = ccoef;
736                         mode_list[  list_start] = 3;
737                     } else {
738                         nz_coeff[nz_coeff_count++] = bink_scan[ccoef];
739                         sign = -get_bits1(gb);
740                         block[bink_scan[ccoef]] = (mask ^ sign) - sign;
741                         masks_count--;
742                         if (masks_count < 0)
743                             return 0;
744                     }
745                 }
746                 break;
747             case 1:
748                 mode_list[list_pos] = 2;
749                 for (i = 0; i < 3; i++) {
750                     ccoef += 4;
751                     coef_list[list_end]   = ccoef;
752                     mode_list[list_end++] = 2;
753                 }
754                 break;
755             case 3:
756                 nz_coeff[nz_coeff_count++] = bink_scan[ccoef];
757                 sign = -get_bits1(gb);
758                 block[bink_scan[ccoef]] = (mask ^ sign) - sign;
759                 coef_list[list_pos]   = 0;
760                 mode_list[list_pos++] = 0;
761                 masks_count--;
762                 if (masks_count < 0)
763                     return 0;
764                 break;
765             }
766         }
767     }
768
769     return 0;
770 }
771
772 /**
773  * Copy 8x8 block from source to destination, where src and dst may be overlapped
774  */
775 static inline void put_pixels8x8_overlapped(uint8_t *dst, uint8_t *src, int stride)
776 {
777     uint8_t tmp[64];
778     int i;
779     for (i = 0; i < 8; i++)
780         memcpy(tmp + i*8, src + i*stride, 8);
781     for (i = 0; i < 8; i++)
782         memcpy(dst + i*stride, tmp + i*8, 8);
783 }
784
785 static int binkb_decode_plane(BinkContext *c, GetBitContext *gb, int plane_idx,
786                               int is_key, int is_chroma)
787 {
788     int blk;
789     int i, j, bx, by;
790     uint8_t *dst, *ref, *ref_start, *ref_end;
791     int v, col[2];
792     const uint8_t *scan;
793     int xoff, yoff;
794     LOCAL_ALIGNED_16(DCTELEM, block, [64]);
795     LOCAL_ALIGNED_16(int32_t, dctblock, [64]);
796     int coordmap[64];
797     int ybias = is_key ? -15 : 0;
798     int qp;
799
800     const int stride = c->pic.linesize[plane_idx];
801     int bw = is_chroma ? (c->avctx->width  + 15) >> 4 : (c->avctx->width  + 7) >> 3;
802     int bh = is_chroma ? (c->avctx->height + 15) >> 4 : (c->avctx->height + 7) >> 3;
803
804     binkb_init_bundles(c);
805     ref_start = c->pic.data[plane_idx];
806     ref_end   = c->pic.data[plane_idx] + (bh * c->pic.linesize[plane_idx] + bw) * 8;
807
808     for (i = 0; i < 64; i++)
809         coordmap[i] = (i & 7) + (i >> 3) * stride;
810
811     for (by = 0; by < bh; by++) {
812         for (i = 0; i < BINKB_NB_SRC; i++) {
813             if (binkb_read_bundle(c, gb, i) < 0)
814                 return -1;
815         }
816
817         dst  = c->pic.data[plane_idx]  + 8*by*stride;
818         for (bx = 0; bx < bw; bx++, dst += 8) {
819             blk = binkb_get_value(c, BINKB_SRC_BLOCK_TYPES);
820             switch (blk) {
821             case 0:
822                 break;
823             case 1:
824                 scan = bink_patterns[get_bits(gb, 4)];
825                 i = 0;
826                 do {
827                     int mode, run;
828
829                     mode = get_bits1(gb);
830                     run = get_bits(gb, binkb_runbits[i]) + 1;
831
832                     i += run;
833                     if (i > 64) {
834                         av_log(c->avctx, AV_LOG_ERROR, "Run went out of bounds\n");
835                         return -1;
836                     }
837                     if (mode) {
838                         v = binkb_get_value(c, BINKB_SRC_COLORS);
839                         for (j = 0; j < run; j++)
840                             dst[coordmap[*scan++]] = v;
841                     } else {
842                         for (j = 0; j < run; j++)
843                             dst[coordmap[*scan++]] = binkb_get_value(c, BINKB_SRC_COLORS);
844                     }
845                 } while (i < 63);
846                 if (i == 63)
847                     dst[coordmap[*scan++]] = binkb_get_value(c, BINKB_SRC_COLORS);
848                 break;
849             case 2:
850                 memset(dctblock, 0, sizeof(*dctblock) * 64);
851                 dctblock[0] = binkb_get_value(c, BINKB_SRC_INTRA_DC);
852                 qp = binkb_get_value(c, BINKB_SRC_INTRA_Q);
853                 read_dct_coeffs(gb, dctblock, bink_scan, binkb_intra_quant, qp);
854                 c->bdsp.idct_put(dst, stride, dctblock);
855                 break;
856             case 3:
857                 xoff = binkb_get_value(c, BINKB_SRC_X_OFF);
858                 yoff = binkb_get_value(c, BINKB_SRC_Y_OFF) + ybias;
859                 ref = dst + xoff + yoff * stride;
860                 if (ref < ref_start || ref + 8*stride > ref_end) {
861                     av_log(c->avctx, AV_LOG_WARNING, "Reference block is out of bounds\n");
862                 } else if (ref + 8*stride < dst || ref >= dst + 8*stride) {
863                     c->dsp.put_pixels_tab[1][0](dst, ref, stride, 8);
864                 } else {
865                     put_pixels8x8_overlapped(dst, ref, stride);
866                 }
867                 c->dsp.clear_block(block);
868                 v = binkb_get_value(c, BINKB_SRC_INTER_COEFS);
869                 read_residue(gb, block, v);
870                 c->dsp.add_pixels8(dst, block, stride);
871                 break;
872             case 4:
873                 xoff = binkb_get_value(c, BINKB_SRC_X_OFF);
874                 yoff = binkb_get_value(c, BINKB_SRC_Y_OFF) + ybias;
875                 ref = dst + xoff + yoff * stride;
876                 if (ref < ref_start || ref + 8 * stride > ref_end) {
877                     av_log(c->avctx, AV_LOG_WARNING, "Reference block is out of bounds\n");
878                 } else if (ref + 8*stride < dst || ref >= dst + 8*stride) {
879                     c->dsp.put_pixels_tab[1][0](dst, ref, stride, 8);
880                 } else {
881                     put_pixels8x8_overlapped(dst, ref, stride);
882                 }
883                 memset(dctblock, 0, sizeof(*dctblock) * 64);
884                 dctblock[0] = binkb_get_value(c, BINKB_SRC_INTER_DC);
885                 qp = binkb_get_value(c, BINKB_SRC_INTER_Q);
886                 read_dct_coeffs(gb, dctblock, bink_scan, binkb_inter_quant, qp);
887                 c->bdsp.idct_add(dst, stride, dctblock);
888                 break;
889             case 5:
890                 v = binkb_get_value(c, BINKB_SRC_COLORS);
891                 c->dsp.fill_block_tab[1](dst, v, stride, 8);
892                 break;
893             case 6:
894                 for (i = 0; i < 2; i++)
895                     col[i] = binkb_get_value(c, BINKB_SRC_COLORS);
896                 for (i = 0; i < 8; i++) {
897                     v = binkb_get_value(c, BINKB_SRC_PATTERN);
898                     for (j = 0; j < 8; j++, v >>= 1)
899                         dst[i*stride + j] = col[v & 1];
900                 }
901                 break;
902             case 7:
903                 xoff = binkb_get_value(c, BINKB_SRC_X_OFF);
904                 yoff = binkb_get_value(c, BINKB_SRC_Y_OFF) + ybias;
905                 ref = dst + xoff + yoff * stride;
906                 if (ref < ref_start || ref + 8 * stride > ref_end) {
907                     av_log(c->avctx, AV_LOG_WARNING, "Reference block is out of bounds\n");
908                 } else if (ref + 8*stride < dst || ref >= dst + 8*stride) {
909                     c->dsp.put_pixels_tab[1][0](dst, ref, stride, 8);
910                 } else {
911                     put_pixels8x8_overlapped(dst, ref, stride);
912                 }
913                 break;
914             case 8:
915                 for (i = 0; i < 8; i++)
916                     memcpy(dst + i*stride, c->bundle[BINKB_SRC_COLORS].cur_ptr + i*8, 8);
917                 c->bundle[BINKB_SRC_COLORS].cur_ptr += 64;
918                 break;
919             default:
920                 av_log(c->avctx, AV_LOG_ERROR, "Unknown block type %d\n", blk);
921                 return -1;
922             }
923         }
924     }
925     if (get_bits_count(gb) & 0x1F) //next plane data starts at 32-bit boundary
926         skip_bits_long(gb, 32 - (get_bits_count(gb) & 0x1F));
927
928     return 0;
929 }
930
931 static int bink_decode_plane(BinkContext *c, GetBitContext *gb, int plane_idx,
932                              int is_chroma)
933 {
934     int blk;
935     int i, j, bx, by;
936     uint8_t *dst, *prev, *ref, *ref_start, *ref_end;
937     int v, col[2];
938     const uint8_t *scan;
939     int xoff, yoff;
940     LOCAL_ALIGNED_16(DCTELEM, block, [64]);
941     LOCAL_ALIGNED_16(uint8_t, ublock, [64]);
942     LOCAL_ALIGNED_16(int32_t, dctblock, [64]);
943     int coordmap[64];
944
945     const int stride = c->pic.linesize[plane_idx];
946     int bw = is_chroma ? (c->avctx->width  + 15) >> 4 : (c->avctx->width  + 7) >> 3;
947     int bh = is_chroma ? (c->avctx->height + 15) >> 4 : (c->avctx->height + 7) >> 3;
948     int width = c->avctx->width >> is_chroma;
949
950     init_lengths(c, FFMAX(width, 8), bw);
951     for (i = 0; i < BINK_NB_SRC; i++)
952         read_bundle(gb, c, i);
953
954     ref_start = c->last.data[plane_idx];
955     ref_end   = c->last.data[plane_idx]
956                 + (bw - 1 + c->last.linesize[plane_idx] * (bh - 1)) * 8;
957
958     for (i = 0; i < 64; i++)
959         coordmap[i] = (i & 7) + (i >> 3) * stride;
960
961     for (by = 0; by < bh; by++) {
962         if (read_block_types(c->avctx, gb, &c->bundle[BINK_SRC_BLOCK_TYPES]) < 0)
963             return -1;
964         if (read_block_types(c->avctx, gb, &c->bundle[BINK_SRC_SUB_BLOCK_TYPES]) < 0)
965             return -1;
966         if (read_colors(gb, &c->bundle[BINK_SRC_COLORS], c) < 0)
967             return -1;
968         if (read_patterns(c->avctx, gb, &c->bundle[BINK_SRC_PATTERN]) < 0)
969             return -1;
970         if (read_motion_values(c->avctx, gb, &c->bundle[BINK_SRC_X_OFF]) < 0)
971             return -1;
972         if (read_motion_values(c->avctx, gb, &c->bundle[BINK_SRC_Y_OFF]) < 0)
973             return -1;
974         if (read_dcs(c->avctx, gb, &c->bundle[BINK_SRC_INTRA_DC], DC_START_BITS, 0) < 0)
975             return -1;
976         if (read_dcs(c->avctx, gb, &c->bundle[BINK_SRC_INTER_DC], DC_START_BITS, 1) < 0)
977             return -1;
978         if (read_runs(c->avctx, gb, &c->bundle[BINK_SRC_RUN]) < 0)
979             return -1;
980
981         if (by == bh)
982             break;
983         dst  = c->pic.data[plane_idx]  + 8*by*stride;
984         prev = c->last.data[plane_idx] + 8*by*stride;
985         for (bx = 0; bx < bw; bx++, dst += 8, prev += 8) {
986             blk = get_value(c, BINK_SRC_BLOCK_TYPES);
987             // 16x16 block type on odd line means part of the already decoded block, so skip it
988             if ((by & 1) && blk == SCALED_BLOCK) {
989                 bx++;
990                 dst  += 8;
991                 prev += 8;
992                 continue;
993             }
994             switch (blk) {
995             case SKIP_BLOCK:
996                 c->dsp.put_pixels_tab[1][0](dst, prev, stride, 8);
997                 break;
998             case SCALED_BLOCK:
999                 blk = get_value(c, BINK_SRC_SUB_BLOCK_TYPES);
1000                 switch (blk) {
1001                 case RUN_BLOCK:
1002                     scan = bink_patterns[get_bits(gb, 4)];
1003                     i = 0;
1004                     do {
1005                         int run = get_value(c, BINK_SRC_RUN) + 1;
1006
1007                         i += run;
1008                         if (i > 64) {
1009                             av_log(c->avctx, AV_LOG_ERROR, "Run went out of bounds\n");
1010                             return -1;
1011                         }
1012                         if (get_bits1(gb)) {
1013                             v = get_value(c, BINK_SRC_COLORS);
1014                             for (j = 0; j < run; j++)
1015                                 ublock[*scan++] = v;
1016                         } else {
1017                             for (j = 0; j < run; j++)
1018                                 ublock[*scan++] = get_value(c, BINK_SRC_COLORS);
1019                         }
1020                     } while (i < 63);
1021                     if (i == 63)
1022                         ublock[*scan++] = get_value(c, BINK_SRC_COLORS);
1023                     break;
1024                 case INTRA_BLOCK:
1025                     memset(dctblock, 0, sizeof(*dctblock) * 64);
1026                     dctblock[0] = get_value(c, BINK_SRC_INTRA_DC);
1027                     read_dct_coeffs(gb, dctblock, bink_scan, bink_intra_quant, -1);
1028                     c->bdsp.idct_put(ublock, 8, dctblock);
1029                     break;
1030                 case FILL_BLOCK:
1031                     v = get_value(c, BINK_SRC_COLORS);
1032                     c->dsp.fill_block_tab[0](dst, v, stride, 16);
1033                     break;
1034                 case PATTERN_BLOCK:
1035                     for (i = 0; i < 2; i++)
1036                         col[i] = get_value(c, BINK_SRC_COLORS);
1037                     for (j = 0; j < 8; j++) {
1038                         v = get_value(c, BINK_SRC_PATTERN);
1039                         for (i = 0; i < 8; i++, v >>= 1)
1040                             ublock[i + j*8] = col[v & 1];
1041                     }
1042                     break;
1043                 case RAW_BLOCK:
1044                     for (j = 0; j < 8; j++)
1045                         for (i = 0; i < 8; i++)
1046                             ublock[i + j*8] = get_value(c, BINK_SRC_COLORS);
1047                     break;
1048                 default:
1049                     av_log(c->avctx, AV_LOG_ERROR, "Incorrect 16x16 block type %d\n", blk);
1050                     return -1;
1051                 }
1052                 if (blk != FILL_BLOCK)
1053                 c->bdsp.scale_block(ublock, dst, stride);
1054                 bx++;
1055                 dst  += 8;
1056                 prev += 8;
1057                 break;
1058             case MOTION_BLOCK:
1059                 xoff = get_value(c, BINK_SRC_X_OFF);
1060                 yoff = get_value(c, BINK_SRC_Y_OFF);
1061                 ref = prev + xoff + yoff * stride;
1062                 if (ref < ref_start || ref > ref_end) {
1063                     av_log(c->avctx, AV_LOG_ERROR, "Copy out of bounds @%d, %d\n",
1064                            bx*8 + xoff, by*8 + yoff);
1065                     return -1;
1066                 }
1067                 c->dsp.put_pixels_tab[1][0](dst, ref, stride, 8);
1068                 break;
1069             case RUN_BLOCK:
1070                 scan = bink_patterns[get_bits(gb, 4)];
1071                 i = 0;
1072                 do {
1073                     int run = get_value(c, BINK_SRC_RUN) + 1;
1074
1075                     i += run;
1076                     if (i > 64) {
1077                         av_log(c->avctx, AV_LOG_ERROR, "Run went out of bounds\n");
1078                         return -1;
1079                     }
1080                     if (get_bits1(gb)) {
1081                         v = get_value(c, BINK_SRC_COLORS);
1082                         for (j = 0; j < run; j++)
1083                             dst[coordmap[*scan++]] = v;
1084                     } else {
1085                         for (j = 0; j < run; j++)
1086                             dst[coordmap[*scan++]] = get_value(c, BINK_SRC_COLORS);
1087                     }
1088                 } while (i < 63);
1089                 if (i == 63)
1090                     dst[coordmap[*scan++]] = get_value(c, BINK_SRC_COLORS);
1091                 break;
1092             case RESIDUE_BLOCK:
1093                 xoff = get_value(c, BINK_SRC_X_OFF);
1094                 yoff = get_value(c, BINK_SRC_Y_OFF);
1095                 ref = prev + xoff + yoff * stride;
1096                 if (ref < ref_start || ref > ref_end) {
1097                     av_log(c->avctx, AV_LOG_ERROR, "Copy out of bounds @%d, %d\n",
1098                            bx*8 + xoff, by*8 + yoff);
1099                     return -1;
1100                 }
1101                 c->dsp.put_pixels_tab[1][0](dst, ref, stride, 8);
1102                 c->dsp.clear_block(block);
1103                 v = get_bits(gb, 7);
1104                 read_residue(gb, block, v);
1105                 c->dsp.add_pixels8(dst, block, stride);
1106                 break;
1107             case INTRA_BLOCK:
1108                 memset(dctblock, 0, sizeof(*dctblock) * 64);
1109                 dctblock[0] = get_value(c, BINK_SRC_INTRA_DC);
1110                 read_dct_coeffs(gb, dctblock, bink_scan, bink_intra_quant, -1);
1111                 c->bdsp.idct_put(dst, stride, dctblock);
1112                 break;
1113             case FILL_BLOCK:
1114                 v = get_value(c, BINK_SRC_COLORS);
1115                 c->dsp.fill_block_tab[1](dst, v, stride, 8);
1116                 break;
1117             case INTER_BLOCK:
1118                 xoff = get_value(c, BINK_SRC_X_OFF);
1119                 yoff = get_value(c, BINK_SRC_Y_OFF);
1120                 ref = prev + xoff + yoff * stride;
1121                 c->dsp.put_pixels_tab[1][0](dst, ref, stride, 8);
1122                 memset(dctblock, 0, sizeof(*dctblock) * 64);
1123                 dctblock[0] = get_value(c, BINK_SRC_INTER_DC);
1124                 read_dct_coeffs(gb, dctblock, bink_scan, bink_inter_quant, -1);
1125                 c->bdsp.idct_add(dst, stride, dctblock);
1126                 break;
1127             case PATTERN_BLOCK:
1128                 for (i = 0; i < 2; i++)
1129                     col[i] = get_value(c, BINK_SRC_COLORS);
1130                 for (i = 0; i < 8; i++) {
1131                     v = get_value(c, BINK_SRC_PATTERN);
1132                     for (j = 0; j < 8; j++, v >>= 1)
1133                         dst[i*stride + j] = col[v & 1];
1134                 }
1135                 break;
1136             case RAW_BLOCK:
1137                 for (i = 0; i < 8; i++)
1138                     memcpy(dst + i*stride, c->bundle[BINK_SRC_COLORS].cur_ptr + i*8, 8);
1139                 c->bundle[BINK_SRC_COLORS].cur_ptr += 64;
1140                 break;
1141             default:
1142                 av_log(c->avctx, AV_LOG_ERROR, "Unknown block type %d\n", blk);
1143                 return -1;
1144             }
1145         }
1146     }
1147     if (get_bits_count(gb) & 0x1F) //next plane data starts at 32-bit boundary
1148         skip_bits_long(gb, 32 - (get_bits_count(gb) & 0x1F));
1149
1150     return 0;
1151 }
1152
1153 static int decode_frame(AVCodecContext *avctx, void *data, int *data_size, AVPacket *pkt)
1154 {
1155     BinkContext * const c = avctx->priv_data;
1156     GetBitContext gb;
1157     int plane, plane_idx;
1158     int bits_count = pkt->size << 3;
1159
1160     if (c->version > 'b') {
1161         if(c->pic.data[0])
1162             avctx->release_buffer(avctx, &c->pic);
1163
1164         if(avctx->get_buffer(avctx, &c->pic) < 0){
1165             av_log(avctx, AV_LOG_ERROR, "get_buffer() failed\n");
1166             return -1;
1167         }
1168     } else {
1169         if(avctx->reget_buffer(avctx, &c->pic) < 0){
1170             av_log(avctx, AV_LOG_ERROR, "reget_buffer() failed\n");
1171             return -1;
1172         }
1173     }
1174
1175     init_get_bits(&gb, pkt->data, bits_count);
1176     if (c->has_alpha) {
1177         if (c->version >= 'i')
1178             skip_bits_long(&gb, 32);
1179         if (bink_decode_plane(c, &gb, 3, 0) < 0)
1180             return -1;
1181     }
1182     if (c->version >= 'i')
1183         skip_bits_long(&gb, 32);
1184
1185     for (plane = 0; plane < 3; plane++) {
1186         plane_idx = (!plane || !c->swap_planes) ? plane : (plane ^ 3);
1187
1188         if (c->version > 'b') {
1189             if (bink_decode_plane(c, &gb, plane_idx, !!plane) < 0)
1190                 return -1;
1191         } else {
1192             if (binkb_decode_plane(c, &gb, plane_idx, !pkt->pts, !!plane) < 0)
1193                 return -1;
1194         }
1195         if (get_bits_count(&gb) >= bits_count)
1196             break;
1197     }
1198     emms_c();
1199
1200     *data_size = sizeof(AVFrame);
1201     *(AVFrame*)data = c->pic;
1202
1203     if (c->version > 'b')
1204         FFSWAP(AVFrame, c->pic, c->last);
1205
1206     /* always report that the buffer was completely consumed */
1207     return pkt->size;
1208 }
1209
1210 /**
1211  * Caclulate quantization tables for version b
1212  */
1213 static av_cold void binkb_calc_quant(void)
1214 {
1215     uint8_t inv_bink_scan[64];
1216     double s[64];
1217     int i, j;
1218
1219     for (j = 0; j < 8; j++) {
1220         for (i = 0; i < 8; i++) {
1221             if (j && j != 4)
1222                if (i && i != 4)
1223                    s[j*8 + i] = cos(j * M_PI/16.0) * cos(i * M_PI/16.0) * 2.0;
1224                else
1225                    s[j*8 + i] = cos(j * M_PI/16.0) * sqrt(2.0);
1226             else
1227                if (i && i != 4)
1228                    s[j*8 + i] = cos(i * M_PI/16.0) * sqrt(2.0);
1229                else
1230                    s[j*8 + i] = 1.0;
1231         }
1232     }
1233
1234     for (i = 0; i < 64; i++)
1235         inv_bink_scan[bink_scan[i]] = i;
1236
1237     for (j = 0; j < 16; j++) {
1238         for (i = 0; i < 64; i++) {
1239             int k = inv_bink_scan[i];
1240             if (s[i] == 1.0) {
1241                 binkb_intra_quant[j][k] = (1L << 12) * binkb_intra_seed[i] *
1242                                           binkb_num[j]/binkb_den[j];
1243                 binkb_inter_quant[j][k] = (1L << 12) * binkb_inter_seed[i] *
1244                                           binkb_num[j]/binkb_den[j];
1245             } else {
1246                 binkb_intra_quant[j][k] = (1L << 12) * binkb_intra_seed[i] * s[i] *
1247                                           binkb_num[j]/(double)binkb_den[j];
1248                 binkb_inter_quant[j][k] = (1L << 12) * binkb_inter_seed[i] * s[i] *
1249                                           binkb_num[j]/(double)binkb_den[j];
1250             }
1251         }
1252     }
1253 }
1254
1255 static av_cold int decode_init(AVCodecContext *avctx)
1256 {
1257     BinkContext * const c = avctx->priv_data;
1258     static VLC_TYPE table[16 * 128][2];
1259     static int binkb_initialised = 0;
1260     int i;
1261     int flags;
1262
1263     c->version = avctx->codec_tag >> 24;
1264     if (avctx->extradata_size < 4) {
1265         av_log(avctx, AV_LOG_ERROR, "Extradata missing or too short\n");
1266         return -1;
1267     }
1268     flags = AV_RL32(avctx->extradata);
1269     c->has_alpha = flags & BINK_FLAG_ALPHA;
1270     c->swap_planes = c->version >= 'h';
1271     if (!bink_trees[15].table) {
1272         for (i = 0; i < 16; i++) {
1273             const int maxbits = bink_tree_lens[i][15];
1274             bink_trees[i].table = table + i*128;
1275             bink_trees[i].table_allocated = 1 << maxbits;
1276             init_vlc(&bink_trees[i], maxbits, 16,
1277                      bink_tree_lens[i], 1, 1,
1278                      bink_tree_bits[i], 1, 1, INIT_VLC_USE_NEW_STATIC | INIT_VLC_LE);
1279         }
1280     }
1281     c->avctx = avctx;
1282
1283     c->pic.data[0] = NULL;
1284
1285     if (av_image_check_size(avctx->width, avctx->height, 0, avctx) < 0) {
1286         return 1;
1287     }
1288
1289     avctx->pix_fmt = c->has_alpha ? PIX_FMT_YUVA420P : PIX_FMT_YUV420P;
1290
1291     avctx->idct_algo = FF_IDCT_BINK;
1292     dsputil_init(&c->dsp, avctx);
1293     ff_binkdsp_init(&c->bdsp);
1294
1295     init_bundles(c);
1296
1297     if (c->version == 'b') {
1298         if (!binkb_initialised) {
1299             binkb_calc_quant();
1300             binkb_initialised = 1;
1301         }
1302     }
1303
1304     return 0;
1305 }
1306
1307 static av_cold int decode_end(AVCodecContext *avctx)
1308 {
1309     BinkContext * const c = avctx->priv_data;
1310
1311     if (c->pic.data[0])
1312         avctx->release_buffer(avctx, &c->pic);
1313     if (c->last.data[0])
1314         avctx->release_buffer(avctx, &c->last);
1315
1316     free_bundles(c);
1317     return 0;
1318 }
1319
1320 AVCodec ff_bink_decoder = {
1321     "binkvideo",
1322     AVMEDIA_TYPE_VIDEO,
1323     CODEC_ID_BINKVIDEO,
1324     sizeof(BinkContext),
1325     decode_init,
1326     NULL,
1327     decode_end,
1328     decode_frame,
1329     .long_name = NULL_IF_CONFIG_SMALL("Bink video"),
1330 };