2 * Amuse Graphics Movie decoder
4 * Copyright (c) 2018 Paul B Mahol
6 * This file is part of FFmpeg.
8 * FFmpeg 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.
13 * FFmpeg 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.
18 * You should have received a copy of the GNU Lesser General Public
19 * License along with FFmpeg; if not, write to the Free Software
20 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
27 #define BITSTREAM_READER_LE
30 #include "bytestream.h"
31 #include "copy_block.h"
36 static const uint8_t unscaled_luma[64] = {
37 16, 11, 10, 16, 24, 40, 51, 61, 12, 12, 14, 19,
38 26, 58, 60, 55, 14, 13, 16, 24, 40, 57, 69, 56,
39 14, 17, 22, 29, 51, 87, 80, 62, 18, 22, 37, 56,
40 68,109,103, 77, 24, 35, 55, 64, 81,104,113, 92,
41 49, 64, 78, 87,103,121,120,101, 72, 92, 95, 98,
45 static const uint8_t unscaled_chroma[64] = {
46 17, 18, 24, 47, 99, 99, 99, 99, 18, 21, 26, 66,
47 99, 99, 99, 99, 24, 26, 56, 99, 99, 99, 99, 99,
48 47, 66, 99, 99, 99, 99, 99, 99, 99, 99, 99, 99,
49 99, 99, 99, 99, 99, 99, 99, 99, 99, 99, 99, 99,
50 99, 99, 99, 99, 99, 99, 99, 99, 99, 99, 99, 99,
54 typedef struct MotionVector {
58 typedef struct AGMContext {
60 AVCodecContext *avctx;
77 unsigned padded_output_size;
80 MotionVector *mvectors;
81 unsigned mvectors_size;
87 int luma_quant_matrix[64];
88 int chroma_quant_matrix[64];
91 DECLARE_ALIGNED(32, int16_t, block)[64];
94 unsigned wblocks_size;
102 static int read_code(GetBitContext *gb, int *oskip, int *level, int *map, int mode)
104 int len = 0, skip = 0, max;
106 if (show_bits(gb, 2)) {
107 switch (show_bits(gb, 4)) {
147 len = ((show_bits(gb, 5) & 0x10) | 0xA0) >> 4;
151 return AVERROR_INVALIDDATA;
155 *level = get_bits(gb, len);
158 max = 1 << (len - 1);
160 *level = -(max + *level);
161 } else if (show_bits(gb, 3) & 4) {
164 if (show_bits(gb, 4)) {
165 if (show_bits(gb, 4) == 1) {
167 *oskip = get_bits(gb, 16);
169 *oskip = get_bits(gb, 4);
173 *oskip = get_bits(gb, 10);
175 } else if (mode == 0) {
176 *oskip = get_bits(gb, 10);
182 *oskip = get_bits(gb, 4);
191 static int decode_intra_blocks(AGMContext *s, GetBitContext *gb,
192 const int *quant_matrix, int *skip, int *dc_level)
194 const uint8_t *scantable = s->scantable.permutated;
195 int level, ret, map = 0;
197 memset(s->wblocks, 0, s->wblocks_size);
199 for (int i = 0; i < 64; i++) {
200 int16_t *block = s->wblocks + scantable[i];
202 for (int j = 0; j < s->blocks_w;) {
206 rskip = FFMIN(*skip, s->blocks_w - j);
209 for (int k = 0; k < rskip; k++)
210 block[64 * k] = *dc_level * quant_matrix[0];
215 ret = read_code(gb, skip, &level, &map, s->flags & 1);
222 block[0] = (i == 0 ? *dc_level : level) * quant_matrix[i];
232 static int decode_inter_blocks(AGMContext *s, GetBitContext *gb,
233 const int *quant_matrix, int *skip,
236 const uint8_t *scantable = s->scantable.permutated;
239 memset(s->wblocks, 0, s->wblocks_size);
240 memset(s->map, 0, s->map_size);
242 for (int i = 0; i < 64; i++) {
243 int16_t *block = s->wblocks + scantable[i];
245 for (int j = 0; j < s->blocks_w;) {
249 rskip = FFMIN(*skip, s->blocks_w - j);
254 ret = read_code(gb, skip, &level, &map[j], s->flags & 1);
258 block[0] = level * quant_matrix[i];
268 static int decode_intra_block(AGMContext *s, GetBitContext *gb,
269 const int *quant_matrix, int *skip, int *dc_level)
271 const uint8_t *scantable = s->scantable.permutated;
272 const int offset = s->plus ? 0 : 1024;
273 int16_t *block = s->block;
274 int level, ret, map = 0;
276 memset(block, 0, sizeof(s->block));
281 ret = read_code(gb, skip, &level, &map, s->flags & 1);
286 block[scantable[0]] = offset + *dc_level * quant_matrix[0];
288 for (int i = 1; i < 64;) {
292 rskip = FFMIN(*skip, 64 - i);
296 ret = read_code(gb, skip, &level, &map, s->flags & 1);
300 block[scantable[i]] = level * quant_matrix[i];
308 static int decode_intra_plane(AGMContext *s, GetBitContext *gb, int size,
309 const int *quant_matrix, AVFrame *frame,
312 int ret, skip = 0, dc_level = 0;
313 const int offset = s->plus ? 0 : 1024;
315 if ((ret = init_get_bits8(gb, s->gbyte.buffer, size)) < 0)
319 av_fast_padded_malloc(&s->wblocks, &s->wblocks_size,
320 64 * s->blocks_w * sizeof(*s->wblocks));
322 return AVERROR(ENOMEM);
324 for (int y = 0; y < s->blocks_h; y++) {
325 ret = decode_intra_blocks(s, gb, quant_matrix, &skip, &dc_level);
329 for (int x = 0; x < s->blocks_w; x++) {
330 s->wblocks[64 * x] += offset;
331 s->idsp.idct_put(frame->data[plane] + (s->blocks_h - 1 - y) * 8 * frame->linesize[plane] + x * 8,
332 frame->linesize[plane], s->wblocks + 64 * x);
336 for (int y = 0; y < s->blocks_h; y++) {
337 for (int x = 0; x < s->blocks_w; x++) {
338 ret = decode_intra_block(s, gb, quant_matrix, &skip, &dc_level);
342 s->idsp.idct_put(frame->data[plane] + (s->blocks_h - 1 - y) * 8 * frame->linesize[plane] + x * 8,
343 frame->linesize[plane], s->block);
349 if (get_bits_left(gb) < 0)
350 av_log(s->avctx, AV_LOG_WARNING, "overread\n");
351 if (get_bits_left(gb) > 0)
352 av_log(s->avctx, AV_LOG_WARNING, "underread: %d\n", get_bits_left(gb));
357 static int decode_inter_block(AGMContext *s, GetBitContext *gb,
358 const int *quant_matrix, int *skip,
361 const uint8_t *scantable = s->scantable.permutated;
362 int16_t *block = s->block;
365 memset(block, 0, sizeof(s->block));
367 for (int i = 0; i < 64;) {
371 rskip = FFMIN(*skip, 64 - i);
375 ret = read_code(gb, skip, &level, map, s->flags & 1);
379 block[scantable[i]] = level * quant_matrix[i];
387 static int decode_inter_plane(AGMContext *s, GetBitContext *gb, int size,
388 const int *quant_matrix, AVFrame *frame,
389 AVFrame *prev, int plane)
393 if ((ret = init_get_bits8(gb, s->gbyte.buffer, size)) < 0)
397 av_fast_padded_malloc(&s->wblocks, &s->wblocks_size,
398 64 * s->blocks_w * sizeof(*s->wblocks));
400 return AVERROR(ENOMEM);
402 av_fast_padded_malloc(&s->map, &s->map_size,
403 s->blocks_w * sizeof(*s->map));
405 return AVERROR(ENOMEM);
407 for (int y = 0; y < s->blocks_h; y++) {
408 ret = decode_inter_blocks(s, gb, quant_matrix, &skip, s->map);
412 for (int x = 0; x < s->blocks_w; x++) {
413 int shift = plane == 0;
414 int mvpos = (y >> shift) * (s->blocks_w >> shift) + (x >> shift);
415 int orig_mv_x = s->mvectors[mvpos].x;
416 int mv_x = s->mvectors[mvpos].x / (1 + !shift);
417 int mv_y = s->mvectors[mvpos].y / (1 + !shift);
418 int h = s->avctx->coded_height >> !shift;
419 int w = s->avctx->coded_width >> !shift;
422 if (orig_mv_x >= -32) {
423 if (y * 8 + mv_y < 0 || y * 8 + mv_y >= h ||
424 x * 8 + mv_x < 0 || x * 8 + mv_x >= w)
425 return AVERROR_INVALIDDATA;
427 copy_block8(frame->data[plane] + (s->blocks_h - 1 - y) * 8 * frame->linesize[plane] + x * 8,
428 prev->data[plane] + ((s->blocks_h - 1 - y) * 8 - mv_y) * prev->linesize[plane] + (x * 8 + mv_x),
429 frame->linesize[plane], prev->linesize[plane], 8);
431 s->idsp.idct(s->wblocks + x * 64);
432 for (int i = 0; i < 64; i++)
433 s->wblocks[i + x * 64] = (s->wblocks[i + x * 64] + 1) & 0xFFFC;
434 s->idsp.add_pixels_clamped(&s->wblocks[x*64], frame->data[plane] + (s->blocks_h - 1 - y) * 8 * frame->linesize[plane] + x * 8,
435 frame->linesize[plane]);
438 s->idsp.idct_put(frame->data[plane] + (s->blocks_h - 1 - y) * 8 * frame->linesize[plane] + x * 8,
439 frame->linesize[plane], s->wblocks + x * 64);
443 } else if (s->flags & 2) {
444 for (int y = 0; y < s->blocks_h; y++) {
445 for (int x = 0; x < s->blocks_w; x++) {
446 int shift = plane == 0;
447 int mvpos = (y >> shift) * (s->blocks_w >> shift) + (x >> shift);
448 int orig_mv_x = s->mvectors[mvpos].x;
449 int mv_x = s->mvectors[mvpos].x / (1 + !shift);
450 int mv_y = s->mvectors[mvpos].y / (1 + !shift);
451 int h = s->avctx->coded_height >> !shift;
452 int w = s->avctx->coded_width >> !shift;
455 ret = decode_inter_block(s, gb, quant_matrix, &skip, &map);
459 if (orig_mv_x >= -32) {
460 if (y * 8 + mv_y < 0 || y * 8 + mv_y >= h ||
461 x * 8 + mv_x < 0 || x * 8 + mv_x >= w)
462 return AVERROR_INVALIDDATA;
464 copy_block8(frame->data[plane] + (s->blocks_h - 1 - y) * 8 * frame->linesize[plane] + x * 8,
465 prev->data[plane] + ((s->blocks_h - 1 - y) * 8 - mv_y) * prev->linesize[plane] + (x * 8 + mv_x),
466 frame->linesize[plane], prev->linesize[plane], 8);
468 s->idsp.idct(s->block);
469 for (int i = 0; i < 64; i++)
470 s->block[i] = (s->block[i] + 1) & 0xFFFC;
471 s->idsp.add_pixels_clamped(s->block, frame->data[plane] + (s->blocks_h - 1 - y) * 8 * frame->linesize[plane] + x * 8,
472 frame->linesize[plane]);
475 s->idsp.idct_put(frame->data[plane] + (s->blocks_h - 1 - y) * 8 * frame->linesize[plane] + x * 8,
476 frame->linesize[plane], s->block);
480 } else if (s->flags & 1) {
481 av_fast_padded_malloc(&s->wblocks, &s->wblocks_size,
482 64 * s->blocks_w * sizeof(*s->wblocks));
484 return AVERROR(ENOMEM);
486 av_fast_padded_malloc(&s->map, &s->map_size,
487 s->blocks_w * sizeof(*s->map));
489 return AVERROR(ENOMEM);
491 for (int y = 0; y < s->blocks_h; y++) {
492 ret = decode_inter_blocks(s, gb, quant_matrix, &skip, s->map);
496 for (int x = 0; x < s->blocks_w; x++) {
499 s->idsp.idct_add(frame->data[plane] + (s->blocks_h - 1 - y) * 8 * frame->linesize[plane] + x * 8,
500 frame->linesize[plane], s->wblocks + 64 * x);
504 for (int y = 0; y < s->blocks_h; y++) {
505 for (int x = 0; x < s->blocks_w; x++) {
508 ret = decode_inter_block(s, gb, quant_matrix, &skip, &map);
514 s->idsp.idct_add(frame->data[plane] + (s->blocks_h - 1 - y) * 8 * frame->linesize[plane] + x * 8,
515 frame->linesize[plane], s->block);
521 if (get_bits_left(gb) < 0)
522 av_log(s->avctx, AV_LOG_WARNING, "overread\n");
523 if (get_bits_left(gb) > 0)
524 av_log(s->avctx, AV_LOG_WARNING, "underread: %d\n", get_bits_left(gb));
529 static void compute_quant_matrix(AGMContext *s, double qscale)
531 int luma[64], chroma[64];
532 double f = 1.0 - fabs(qscale);
534 if (!s->key_frame && (s->flags & 2)) {
536 for (int i = 0; i < 64; i++) {
537 luma[i] = FFMAX(1, 16 * f);
538 chroma[i] = FFMAX(1, 16 * f);
541 for (int i = 0; i < 64; i++) {
542 luma[i] = FFMAX(1, 16 - qscale * 32);
543 chroma[i] = FFMAX(1, 16 - qscale * 32);
548 for (int i = 0; i < 64; i++) {
549 luma[i] = FFMAX(1, unscaled_luma [(i & 7) * 8 + (i >> 3)] * f);
550 chroma[i] = FFMAX(1, unscaled_chroma[(i & 7) * 8 + (i >> 3)] * f);
553 for (int i = 0; i < 64; i++) {
554 luma[i] = FFMAX(1, 255.0 - (255 - unscaled_luma [(i & 7) * 8 + (i >> 3)]) * f);
555 chroma[i] = FFMAX(1, 255.0 - (255 - unscaled_chroma[(i & 7) * 8 + (i >> 3)]) * f);
560 for (int i = 0; i < 64; i++) {
561 int pos = ff_zigzag_direct[i];
563 s->luma_quant_matrix[i] = luma[pos] * ((pos / 8) & 1 ? -1 : 1);
564 s->chroma_quant_matrix[i] = chroma[pos] * ((pos / 8) & 1 ? -1 : 1);
568 static int decode_raw_intra_rgb(AVCodecContext *avctx, GetByteContext *gbyte, AVFrame *frame)
570 uint8_t *dst = frame->data[0] + (avctx->height - 1) * frame->linesize[0];
571 uint8_t r = 0, g = 0, b = 0;
573 for (int y = 0; y < avctx->height; y++) {
574 for (int x = 0; x < avctx->width; x++) {
575 dst[x*3+0] = bytestream2_get_byte(gbyte) + r;
577 dst[x*3+1] = bytestream2_get_byte(gbyte) + g;
579 dst[x*3+2] = bytestream2_get_byte(gbyte) + b;
582 dst -= frame->linesize[0];
588 static int fill_pixels(uint8_t **y0, uint8_t **y1,
589 uint8_t **u, uint8_t **v,
590 int ylinesize, int ulinesize, int vlinesize,
592 int *nx, int *ny, int *np, int w, int h)
594 uint8_t *y0dst = *y0;
595 uint8_t *y1dst = *y1;
598 int x = *nx, y = *ny, pos = *np;
601 y0dst[2*x+0] += fill[0];
602 y0dst[2*x+1] += fill[1];
603 y1dst[2*x+0] += fill[2];
604 y1dst[2*x+1] += fill[3];
606 } else if (pos == 1) {
615 y0dst -= 2*ylinesize;
616 y1dst -= 2*ylinesize;
620 y0dst[2*x+0] += fill[2];
621 y0dst[2*x+1] += fill[3];
623 } else if (pos == 2) {
624 y1dst[2*x+0] += fill[0];
625 y1dst[2*x+1] += fill[1];
634 y0dst -= 2*ylinesize;
635 y1dst -= 2*ylinesize;
653 static int decode_runlen_rgb(AVCodecContext *avctx, GetByteContext *gbyte, AVFrame *frame)
655 uint8_t *dst = frame->data[0] + (avctx->height - 1) * frame->linesize[0];
656 int runlen, y = 0, x = 0;
660 while (bytestream2_get_bytes_left(gbyte) > 0) {
661 code = bytestream2_peek_le32(gbyte);
662 runlen = code & 0xFFFFFF;
664 if (code >> 24 == 0x77) {
665 bytestream2_skip(gbyte, 4);
667 for (int i = 0; i < 4; i++)
668 fill[i] = bytestream2_get_byte(gbyte);
673 for (int i = 0; i < 4; i++) {
676 if (x >= frame->width * 3) {
679 dst -= frame->linesize[0];
680 if (y >= frame->height)
686 for (int i = 0; i < 4; i++)
687 fill[i] = bytestream2_get_byte(gbyte);
689 for (int i = 0; i < 4; i++) {
692 if (x >= frame->width * 3) {
695 dst -= frame->linesize[0];
696 if (y >= frame->height)
706 static int decode_runlen(AVCodecContext *avctx, GetByteContext *gbyte, AVFrame *frame)
708 uint8_t *y0dst = frame->data[0] + (avctx->height - 1) * frame->linesize[0];
709 uint8_t *y1dst = y0dst - frame->linesize[0];
710 uint8_t *udst = frame->data[1] + ((avctx->height >> 1) - 1) * frame->linesize[1];
711 uint8_t *vdst = frame->data[2] + ((avctx->height >> 1) - 1) * frame->linesize[2];
712 int runlen, y = 0, x = 0, pos = 0;
716 while (bytestream2_get_bytes_left(gbyte) > 0) {
717 code = bytestream2_peek_le32(gbyte);
718 runlen = code & 0xFFFFFF;
720 if (code >> 24 == 0x77) {
721 bytestream2_skip(gbyte, 4);
723 for (int i = 0; i < 4; i++)
724 fill[i] = bytestream2_get_byte(gbyte);
729 if (fill_pixels(&y0dst, &y1dst, &udst, &vdst,
739 for (int i = 0; i < 4; i++)
740 fill[i] = bytestream2_get_byte(gbyte);
742 if (fill_pixels(&y0dst, &y1dst, &udst, &vdst,
756 static int decode_raw_intra(AVCodecContext *avctx, GetByteContext *gbyte, AVFrame *frame)
758 uint8_t *y0dst = frame->data[0] + (avctx->height - 1) * frame->linesize[0];
759 uint8_t *y1dst = y0dst - frame->linesize[0];
760 uint8_t *udst = frame->data[1] + ((avctx->height >> 1) - 1) * frame->linesize[1];
761 uint8_t *vdst = frame->data[2] + ((avctx->height >> 1) - 1) * frame->linesize[2];
762 uint8_t ly0 = 0, ly1 = 0, ly2 = 0, ly3 = 0, lu = 0, lv = 0;
764 for (int y = 0; y < avctx->height / 2; y++) {
765 for (int x = 0; x < avctx->width / 2; x++) {
766 y0dst[x*2+0] = bytestream2_get_byte(gbyte) + ly0;
768 y0dst[x*2+1] = bytestream2_get_byte(gbyte) + ly1;
770 y1dst[x*2+0] = bytestream2_get_byte(gbyte) + ly2;
772 y1dst[x*2+1] = bytestream2_get_byte(gbyte) + ly3;
774 udst[x] = bytestream2_get_byte(gbyte) + lu;
776 vdst[x] = bytestream2_get_byte(gbyte) + lv;
780 y0dst -= 2*frame->linesize[0];
781 y1dst -= 2*frame->linesize[0];
782 udst -= frame->linesize[1];
783 vdst -= frame->linesize[2];
789 static int decode_intra(AVCodecContext *avctx, GetBitContext *gb, AVFrame *frame)
791 AGMContext *s = avctx->priv_data;
794 compute_quant_matrix(s, (2 * s->compression - 100) / 100.0);
796 s->blocks_w = avctx->coded_width >> 3;
797 s->blocks_h = avctx->coded_height >> 3;
799 ret = decode_intra_plane(s, gb, s->size[0], s->luma_quant_matrix, frame, 0);
803 bytestream2_skip(&s->gbyte, s->size[0]);
805 s->blocks_w = avctx->coded_width >> 4;
806 s->blocks_h = avctx->coded_height >> 4;
808 ret = decode_intra_plane(s, gb, s->size[1], s->chroma_quant_matrix, frame, 2);
812 bytestream2_skip(&s->gbyte, s->size[1]);
814 s->blocks_w = avctx->coded_width >> 4;
815 s->blocks_h = avctx->coded_height >> 4;
817 ret = decode_intra_plane(s, gb, s->size[2], s->chroma_quant_matrix, frame, 1);
824 static int decode_motion_vectors(AVCodecContext *avctx, GetBitContext *gb)
826 AGMContext *s = avctx->priv_data;
827 int nb_mvs = ((avctx->height + 15) >> 4) * ((avctx->width + 15) >> 4);
828 int ret, skip = 0, value, map;
830 av_fast_padded_malloc(&s->mvectors, &s->mvectors_size,
831 nb_mvs * sizeof(*s->mvectors));
833 return AVERROR(ENOMEM);
835 if ((ret = init_get_bits8(gb, s->gbyte.buffer, bytestream2_get_bytes_left(&s->gbyte) -
836 (s->size[0] + s->size[1] + s->size[2]))) < 0)
839 memset(s->mvectors, 0, sizeof(*s->mvectors) * nb_mvs);
841 for (int i = 0; i < nb_mvs; i++) {
842 ret = read_code(gb, &skip, &value, &map, 1);
845 s->mvectors[i].x = value;
849 for (int i = 0; i < nb_mvs; i++) {
850 ret = read_code(gb, &skip, &value, &map, 1);
853 s->mvectors[i].y = value;
857 if (get_bits_left(gb) <= 0)
858 return AVERROR_INVALIDDATA;
859 skip = (get_bits_count(gb) >> 3) + 1;
860 bytestream2_skip(&s->gbyte, skip);
865 static int decode_inter(AVCodecContext *avctx, GetBitContext *gb,
866 AVFrame *frame, AVFrame *prev)
868 AGMContext *s = avctx->priv_data;
871 compute_quant_matrix(s, (2 * s->compression - 100) / 100.0);
874 ret = decode_motion_vectors(avctx, gb);
879 s->blocks_w = avctx->coded_width >> 3;
880 s->blocks_h = avctx->coded_height >> 3;
882 ret = decode_inter_plane(s, gb, s->size[0], s->luma_quant_matrix, frame, prev, 0);
886 bytestream2_skip(&s->gbyte, s->size[0]);
888 s->blocks_w = avctx->coded_width >> 4;
889 s->blocks_h = avctx->coded_height >> 4;
891 ret = decode_inter_plane(s, gb, s->size[1], s->chroma_quant_matrix, frame, prev, 2);
895 bytestream2_skip(&s->gbyte, s->size[1]);
897 s->blocks_w = avctx->coded_width >> 4;
898 s->blocks_h = avctx->coded_height >> 4;
900 ret = decode_inter_plane(s, gb, s->size[2], s->chroma_quant_matrix, frame, prev, 1);
907 typedef struct Node {
912 static void get_tree_codes(uint32_t *codes, Node *nodes, int idx, uint32_t pfx, int bitpos)
914 if (idx < 256 && idx >= 0) {
917 get_tree_codes(codes, nodes, nodes[idx].child[0], pfx + (0 << bitpos), bitpos + 1);
918 get_tree_codes(codes, nodes, nodes[idx].child[1], pfx + (1 << bitpos), bitpos + 1);
922 static void make_new_tree(const uint8_t *bitlens, uint32_t *codes)
924 int zlcount = 0, curlen, idx, nindex, last, llast;
925 int blcounts[32] = { 0 };
931 for (int i = 0; i < 256; i++) {
932 int bitlen = bitlens[i];
933 int blcount = blcounts[bitlen];
935 zlcount += bitlen < 1;
936 syms[(bitlen << 8) + blcount] = i;
940 for (int i = 0; i < 512; i++) {
941 nodes[i].child[0] = -1;
942 nodes[i].child[1] = -1;
945 for (int i = 0; i < 256; i++) {
946 node_idx[i] = 257 + i;;
954 for (curlen = 1; curlen < 32; curlen++) {
955 if (blcounts[curlen] > 0) {
956 int max_zlcount = zlcount + blcounts[curlen];
958 for (int i = 0; zlcount < 256 && zlcount < max_zlcount; zlcount++, i++) {
959 int p = node_idx[nindex - 1 + 512];
960 int ch = syms[256 * curlen + i];
962 if (nodes[p].child[0] == -1) {
963 nodes[p].child[0] = ch;
965 nodes[p].child[1] = ch;
968 nodes[ch].parent = p;
977 p = node_idx[nindex - 1 + 512];
979 if (nodes[p].child[0] == -1) {
980 nodes[p].child[0] = ch;
982 nodes[p].child[1] = ch;
986 nodes[ch].parent = p;
991 for (int i = 0; i < idx; i++)
992 node_idx[512 + i] = old_idx[i];
1000 get_tree_codes(codes, nodes, 256, 0, 0);
1003 static int build_huff(const uint8_t *bitlen, VLC *vlc)
1005 uint32_t new_codes[256];
1007 uint8_t symbols[256];
1008 uint32_t codes[256];
1011 make_new_tree(bitlen, new_codes);
1013 for (int i = 0; i < 256; i++) {
1015 bits[nb_codes] = bitlen[i];
1016 codes[nb_codes] = new_codes[i];
1017 symbols[nb_codes] = i;
1023 return ff_init_vlc_sparse(vlc, 13, nb_codes,
1030 static int decode_huffman2(AVCodecContext *avctx, int header, int size)
1032 AGMContext *s = avctx->priv_data;
1033 GetBitContext *gb = &s->gb;
1037 if ((ret = init_get_bits8(gb, s->gbyte.buffer,
1038 bytestream2_get_bytes_left(&s->gbyte))) < 0)
1041 s->output_size = get_bits_long(gb, 32);
1043 av_fast_padded_malloc(&s->output, &s->padded_output_size, s->output_size);
1045 return AVERROR(ENOMEM);
1047 x = get_bits(gb, 1);
1048 len = 4 + get_bits(gb, 1);
1051 int count = get_bits(gb, 3) + 1;
1053 for (int i = 0; i < count; i++)
1054 cb[i] = get_bits(gb, len);
1056 for (int i = 0; i < 256; i++) {
1057 int idx = get_bits(gb, 3);
1061 for (int i = 0; i < 256; i++)
1062 lens[i] = get_bits(gb, len);
1065 if ((ret = build_huff(lens, &s->vlc)) < 0)
1069 while (get_bits_left(gb) > 0 && x < s->output_size) {
1070 int val = get_vlc2(gb, s->vlc.table, s->vlc.bits, 3);
1072 return AVERROR_INVALIDDATA;
1073 s->output[x++] = val;
1079 static int decode_frame(AVCodecContext *avctx, void *data,
1080 int *got_frame, AVPacket *avpkt)
1082 AGMContext *s = avctx->priv_data;
1083 GetBitContext *gb = &s->gb;
1084 GetByteContext *gbyte = &s->gbyte;
1085 AVFrame *frame = data;
1086 int w, h, width, height, header;
1087 unsigned compressed_size;
1094 bytestream2_init(gbyte, avpkt->data, avpkt->size);
1096 header = bytestream2_get_le32(gbyte);
1097 s->fflags = bytestream2_get_le32(gbyte);
1098 s->bitstream_size = s->fflags & 0x1FFFFFFF;
1100 av_log(avctx, AV_LOG_DEBUG, "fflags: %X\n", s->fflags);
1101 if (avpkt->size < s->bitstream_size + 8)
1102 return AVERROR_INVALIDDATA;
1104 s->key_frame = (avpkt->flags & AV_PKT_FLAG_KEY);
1105 frame->key_frame = s->key_frame;
1106 frame->pict_type = s->key_frame ? AV_PICTURE_TYPE_I : AV_PICTURE_TYPE_P;
1109 if (avctx->codec_tag == MKTAG('A', 'G', 'M', '0') ||
1110 avctx->codec_tag == MKTAG('A', 'G', 'M', '1'))
1111 return AVERROR_PATCHWELCOME;
1113 ret = decode_huffman2(avctx, header, (avpkt->size - s->bitstream_size) - 8);
1116 bytestream2_init(gbyte, s->output, s->output_size);
1117 } else if (!s->dct) {
1118 bytestream2_skip(gbyte, 4);
1123 w = bytestream2_get_le32(gbyte);
1124 h = bytestream2_get_le32(gbyte);
1125 if (w == INT32_MIN || h == INT32_MIN)
1126 return AVERROR_INVALIDDATA;
1136 width = avctx->width;
1137 height = avctx->height;
1138 if (w < width || h < height || w & 7 || h & 7)
1139 return AVERROR_INVALIDDATA;
1141 ret = ff_set_dimensions(avctx, w, h);
1144 avctx->width = width;
1145 avctx->height = height;
1147 s->compression = bytestream2_get_le32(gbyte);
1148 if (s->compression < 0 || s->compression > 100)
1149 return AVERROR_INVALIDDATA;
1151 for (int i = 0; i < 3; i++)
1152 s->size[i] = bytestream2_get_le32(gbyte);
1154 compressed_size = s->output_size;
1157 compressed_size = avpkt->size;
1160 if (s->size[0] < 0 || s->size[1] < 0 || s->size[2] < 0 ||
1161 skip + s->size[0] + s->size[1] + s->size[2] > compressed_size) {
1162 return AVERROR_INVALIDDATA;
1166 if ((ret = ff_get_buffer(avctx, frame, AV_GET_BUFFER_FLAG_REF)) < 0)
1169 if (frame->key_frame) {
1170 if (!s->dct && !s->rgb)
1171 ret = decode_raw_intra(avctx, gbyte, frame);
1172 else if (!s->dct && s->rgb)
1173 ret = decode_raw_intra_rgb(avctx, gbyte, frame);
1175 ret = decode_intra(avctx, gb, frame);
1177 if (!s->prev_frame->data[0]) {
1178 av_log(avctx, AV_LOG_ERROR, "Missing reference frame.\n");
1179 return AVERROR_INVALIDDATA;
1182 if (!(s->flags & 2)) {
1183 ret = av_frame_copy(frame, s->prev_frame);
1189 ret = decode_inter(avctx, gb, frame, s->prev_frame);
1190 } else if (!s->dct && !s->rgb) {
1191 ret = decode_runlen(avctx, gbyte, frame);
1193 ret = decode_runlen_rgb(avctx, gbyte, frame);
1199 av_frame_unref(s->prev_frame);
1200 if ((ret = av_frame_ref(s->prev_frame, frame)) < 0)
1203 frame->crop_top = avctx->coded_height - avctx->height;
1204 frame->crop_left = avctx->coded_width - avctx->width;
1211 static av_cold int decode_init(AVCodecContext *avctx)
1213 AGMContext *s = avctx->priv_data;
1215 s->rgb = avctx->codec_tag == MKTAG('A', 'G', 'M', '4');
1216 avctx->pix_fmt = s->rgb ? AV_PIX_FMT_BGR24 : AV_PIX_FMT_YUV420P;
1218 s->plus = avctx->codec_tag == MKTAG('A', 'G', 'M', '3') ||
1219 avctx->codec_tag == MKTAG('A', 'G', 'M', '7');
1221 s->dct = avctx->codec_tag != MKTAG('A', 'G', 'M', '4') &&
1222 avctx->codec_tag != MKTAG('A', 'G', 'M', '5');
1224 avctx->idct_algo = FF_IDCT_SIMPLE;
1225 ff_idctdsp_init(&s->idsp, avctx);
1226 ff_init_scantable(s->idsp.idct_permutation, &s->scantable, ff_zigzag_direct);
1228 s->prev_frame = av_frame_alloc();
1230 return AVERROR(ENOMEM);
1235 static void decode_flush(AVCodecContext *avctx)
1237 AGMContext *s = avctx->priv_data;
1239 av_frame_unref(s->prev_frame);
1242 static av_cold int decode_close(AVCodecContext *avctx)
1244 AGMContext *s = avctx->priv_data;
1246 ff_free_vlc(&s->vlc);
1247 av_frame_free(&s->prev_frame);
1248 av_freep(&s->mvectors);
1249 s->mvectors_size = 0;
1250 av_freep(&s->wblocks);
1251 s->wblocks_size = 0;
1252 av_freep(&s->output);
1253 s->padded_output_size = 0;
1260 AVCodec ff_agm_decoder = {
1262 .long_name = NULL_IF_CONFIG_SMALL("Amuse Graphics Movie"),
1263 .type = AVMEDIA_TYPE_VIDEO,
1264 .id = AV_CODEC_ID_AGM,
1265 .priv_data_size = sizeof(AGMContext),
1266 .init = decode_init,
1267 .close = decode_close,
1268 .decode = decode_frame,
1269 .flush = decode_flush,
1270 .capabilities = AV_CODEC_CAP_DR1,
1271 .caps_internal = FF_CODEC_CAP_INIT_THREADSAFE |
1272 FF_CODEC_CAP_INIT_CLEANUP |
1273 FF_CODEC_CAP_EXPORTS_CROPPING,