falkTX
/
FFmpeg
mirror of https://github.com/falkTX/FFmpeg.git
1
0
Fork 0
Code Issues 0 Releases 338 Wiki Activity
You can not select more than 25 topics Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
75430 Commits
32 Branches
283MB
Tree: ae82d4f0fa
Branches Tags
${ item.name }
Create branch ${ searchTerm }
from 'ae82d4f0fa'
${ noResults }
FFmpeg/libavcodec/vp56.c

843 lines
26KB
Raw Blame History 

  1. /*

  2.  * Copyright (C) 2006  Aurelien Jacobs <aurel@gnuage.org>

  3.  *

  4.  * This file is part of FFmpeg.

  5.  *

  6.  * FFmpeg is free software; you can redistribute it and/or

  7.  * modify it under the terms of the GNU Lesser General Public

  8.  * License as published by the Free Software Foundation; either

  9.  * version 2.1 of the License, or (at your option) any later version.

  10.  *

  11.  * FFmpeg is distributed in the hope that it will be useful,

  12.  * but WITHOUT ANY WARRANTY; without even the implied warranty of

  13.  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU

  14.  * Lesser General Public License for more details.

  15.  *

  16.  * You should have received a copy of the GNU Lesser General Public

  17.  * License along with FFmpeg; if not, write to the Free Software

  18.  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA

  19.  */

  20. 

  21. /**

  22.  * @file

  23.  * VP5 and VP6 compatible video decoder (common features)

  24.  */

  25. 

  26. #include "avcodec.h"

  27. #include "bytestream.h"

  28. #include "internal.h"

  29. #include "h264chroma.h"

  30. #include "vp56.h"

  31. #include "vp56data.h"

  32. 

  33. 

  34. void ff_vp56_init_dequant(VP56Context *s, int quantizer)

  35. {

  36.     s->quantizer = quantizer;

  37.     s->dequant_dc = ff_vp56_dc_dequant[quantizer] << 2;

  38.     s->dequant_ac = ff_vp56_ac_dequant[quantizer] << 2;

  39. }

  40. 

  41. static int vp56_get_vectors_predictors(VP56Context *s, int row, int col,

  42.                                        VP56Frame ref_frame)

  43. {

  44.     int nb_pred = 0;

  45.     VP56mv vect[2] = {{0,0}, {0,0}};

  46.     int pos, offset;

  47.     VP56mv mvp;

  48. 

  49.     for (pos=0; pos<12; pos++) {

  50.         mvp.x = col + ff_vp56_candidate_predictor_pos[pos][0];

  51.         mvp.y = row + ff_vp56_candidate_predictor_pos[pos][1];

  52.         if (mvp.x < 0 || mvp.x >= s->mb_width ||

  53.             mvp.y < 0 || mvp.y >= s->mb_height)

  54.             continue;

  55.         offset = mvp.x + s->mb_width*mvp.y;

  56. 

  57.         if (ff_vp56_reference_frame[s->macroblocks[offset].type] != ref_frame)

  58.             continue;

  59.         if ((s->macroblocks[offset].mv.x == vect[0].x &&

  60.              s->macroblocks[offset].mv.y == vect[0].y) ||

  61.             (s->macroblocks[offset].mv.x == 0 &&

  62.              s->macroblocks[offset].mv.y == 0))

  63.             continue;

  64. 

  65.         vect[nb_pred++] = s->macroblocks[offset].mv;

  66.         if (nb_pred > 1) {

  67.             nb_pred = -1;

  68.             break;

  69.         }

  70.         s->vector_candidate_pos = pos;

  71.     }

  72. 

  73.     s->vector_candidate[0] = vect[0];

  74.     s->vector_candidate[1] = vect[1];

  75. 

  76.     return nb_pred+1;

  77. }

  78. 

  79. static void vp56_parse_mb_type_models(VP56Context *s)

  80. {

  81.     VP56RangeCoder *c = &s->c;

  82.     VP56Model *model = s->modelp;

  83.     int i, ctx, type;

  84. 

  85.     for (ctx=0; ctx<3; ctx++) {

  86.         if (vp56_rac_get_prob_branchy(c, 174)) {

  87.             int idx = vp56_rac_gets(c, 4);

  88.             memcpy(model->mb_types_stats[ctx],

  89.                    ff_vp56_pre_def_mb_type_stats[idx][ctx],

  90.                    sizeof(model->mb_types_stats[ctx]));

  91.         }

  92.         if (vp56_rac_get_prob_branchy(c, 254)) {

  93.             for (type=0; type<10; type++) {

  94.                 for(i=0; i<2; i++) {

  95.                     if (vp56_rac_get_prob_branchy(c, 205)) {

  96.                         int delta, sign = vp56_rac_get(c);

  97. 

  98.                         delta = vp56_rac_get_tree(c, ff_vp56_pmbtm_tree,

  99.                                                   ff_vp56_mb_type_model_model);

  100.                         if (!delta)

  101.                             delta = 4 * vp56_rac_gets(c, 7);

  102.                         model->mb_types_stats[ctx][type][i] += (delta ^ -sign) + sign;

  103.                     }

  104.                 }

  105.             }

  106.         }

  107.     }

  108. 

  109.     /* compute MB type probability tables based on previous MB type */

  110.     for (ctx=0; ctx<3; ctx++) {

  111.         int p[10];

  112. 

  113.         for (type=0; type<10; type++)

  114.             p[type] = 100 * model->mb_types_stats[ctx][type][1];

  115. 

  116.         for (type=0; type<10; type++) {

  117.             int p02, p34, p0234, p17, p56, p89, p5689, p156789;

  118. 

  119.             /* conservative MB type probability */

  120.             model->mb_type[ctx][type][0] = 255 - (255 * model->mb_types_stats[ctx][type][0]) / (1 + model->mb_types_stats[ctx][type][0] + model->mb_types_stats[ctx][type][1]);

  121. 

  122.             p[type] = 0;    /* same MB type => weight is null */

  123. 

  124.             /* binary tree parsing probabilities */

  125.             p02 = p[0] + p[2];

  126.             p34 = p[3] + p[4];

  127.             p0234 = p02 + p34;

  128.             p17 = p[1] + p[7];

  129.             p56 = p[5] + p[6];

  130.             p89 = p[8] + p[9];

  131.             p5689 = p56 + p89;

  132.             p156789 = p17 + p5689;

  133. 

  134.             model->mb_type[ctx][type][1] = 1 + 255 * p0234/(1+p0234+p156789);

  135.             model->mb_type[ctx][type][2] = 1 + 255 * p02  / (1+p0234);

  136.             model->mb_type[ctx][type][3] = 1 + 255 * p17  / (1+p156789);

  137.             model->mb_type[ctx][type][4] = 1 + 255 * p[0] / (1+p02);

  138.             model->mb_type[ctx][type][5] = 1 + 255 * p[3] / (1+p34);

  139.             model->mb_type[ctx][type][6] = 1 + 255 * p[1] / (1+p17);

  140.             model->mb_type[ctx][type][7] = 1 + 255 * p56  / (1+p5689);

  141.             model->mb_type[ctx][type][8] = 1 + 255 * p[5] / (1+p56);

  142.             model->mb_type[ctx][type][9] = 1 + 255 * p[8] / (1+p89);

  143. 

  144.             /* restore initial value */

  145.             p[type] = 100 * model->mb_types_stats[ctx][type][1];

  146.         }

  147.     }

  148. }

  149. 

  150. static VP56mb vp56_parse_mb_type(VP56Context *s,

  151.                                  VP56mb prev_type, int ctx)

  152. {

  153.     uint8_t *mb_type_model = s->modelp->mb_type[ctx][prev_type];

  154.     VP56RangeCoder *c = &s->c;

  155. 

  156.     if (vp56_rac_get_prob_branchy(c, mb_type_model[0]))

  157.         return prev_type;

  158.     else

  159.         return vp56_rac_get_tree(c, ff_vp56_pmbt_tree, mb_type_model);

  160. }

  161. 

  162. static void vp56_decode_4mv(VP56Context *s, int row, int col)

  163. {

  164.     VP56mv mv = {0,0};

  165.     int type[4];

  166.     int b;

  167. 

  168.     /* parse each block type */

  169.     for (b=0; b<4; b++) {

  170.         type[b] = vp56_rac_gets(&s->c, 2);

  171.         if (type[b])

  172.             type[b]++;  /* only returns 0, 2, 3 or 4 (all INTER_PF) */

  173.     }

  174. 

  175.     /* get vectors */

  176.     for (b=0; b<4; b++) {

  177.         switch (type[b]) {

  178.             case VP56_MB_INTER_NOVEC_PF:

  179.                 s->mv[b] = (VP56mv) {0,0};

  180.                 break;

  181.             case VP56_MB_INTER_DELTA_PF:

  182.                 s->parse_vector_adjustment(s, &s->mv[b]);

  183.                 break;

  184.             case VP56_MB_INTER_V1_PF:

  185.                 s->mv[b] = s->vector_candidate[0];

  186.                 break;

  187.             case VP56_MB_INTER_V2_PF:

  188.                 s->mv[b] = s->vector_candidate[1];

  189.                 break;

  190.         }

  191.         mv.x += s->mv[b].x;

  192.         mv.y += s->mv[b].y;

  193.     }

  194. 

  195.     /* this is the one selected for the whole MB for prediction */

  196.     s->macroblocks[row * s->mb_width + col].mv = s->mv[3];

  197. 

  198.     /* chroma vectors are average luma vectors */

  199.     if (s->avctx->codec->id == AV_CODEC_ID_VP5) {

  200.         s->mv[4].x = s->mv[5].x = RSHIFT(mv.x,2);

  201.         s->mv[4].y = s->mv[5].y = RSHIFT(mv.y,2);

  202.     } else {

  203.         s->mv[4] = s->mv[5] = (VP56mv) {mv.x/4, mv.y/4};

  204.     }

  205. }

  206. 

  207. static VP56mb vp56_decode_mv(VP56Context *s, int row, int col)

  208. {

  209.     VP56mv *mv, vect = {0,0};

  210.     int ctx, b;

  211. 

  212.     ctx = vp56_get_vectors_predictors(s, row, col, VP56_FRAME_PREVIOUS);

  213.     s->mb_type = vp56_parse_mb_type(s, s->mb_type, ctx);

  214.     s->macroblocks[row * s->mb_width + col].type = s->mb_type;

  215. 

  216.     switch (s->mb_type) {

  217.         case VP56_MB_INTER_V1_PF:

  218.             mv = &s->vector_candidate[0];

  219.             break;

  220. 

  221.         case VP56_MB_INTER_V2_PF:

  222.             mv = &s->vector_candidate[1];

  223.             break;

  224. 

  225.         case VP56_MB_INTER_V1_GF:

  226.             vp56_get_vectors_predictors(s, row, col, VP56_FRAME_GOLDEN);

  227.             mv = &s->vector_candidate[0];

  228.             break;

  229. 

  230.         case VP56_MB_INTER_V2_GF:

  231.             vp56_get_vectors_predictors(s, row, col, VP56_FRAME_GOLDEN);

  232.             mv = &s->vector_candidate[1];

  233.             break;

  234. 

  235.         case VP56_MB_INTER_DELTA_PF:

  236.             s->parse_vector_adjustment(s, &vect);

  237.             mv = &vect;

  238.             break;

  239. 

  240.         case VP56_MB_INTER_DELTA_GF:

  241.             vp56_get_vectors_predictors(s, row, col, VP56_FRAME_GOLDEN);

  242.             s->parse_vector_adjustment(s, &vect);

  243.             mv = &vect;

  244.             break;

  245. 

  246.         case VP56_MB_INTER_4V:

  247.             vp56_decode_4mv(s, row, col);

  248.             return s->mb_type;

  249. 

  250.         default:

  251.             mv = &vect;

  252.             break;

  253.     }

  254. 

  255.     s->macroblocks[row*s->mb_width + col].mv = *mv;

  256. 

  257.     /* same vector for all blocks */

  258.     for (b=0; b<6; b++)

  259.         s->mv[b] = *mv;

  260. 

  261.     return s->mb_type;

  262. }

  263. 

  264. static VP56mb vp56_conceal_mv(VP56Context *s, int row, int col)

  265. {

  266.     VP56mv *mv, vect = {0,0};

  267.     int b;

  268. 

  269.     s->mb_type = VP56_MB_INTER_NOVEC_PF;

  270.     s->macroblocks[row * s->mb_width + col].type = s->mb_type;

  271. 

  272.     mv = &vect;

  273. 

  274.     s->macroblocks[row*s->mb_width + col].mv = *mv;

  275. 

  276.     /* same vector for all blocks */

  277.     for (b=0; b<6; b++)

  278.         s->mv[b] = *mv;

  279. 

  280.     return s->mb_type;

  281. }

  282. 

  283. static void vp56_add_predictors_dc(VP56Context *s, VP56Frame ref_frame)

  284. {

  285.     int idx = s->idct_scantable[0];

  286.     int b;

  287. 

  288.     for (b=0; b<6; b++) {

  289.         VP56RefDc *ab = &s->above_blocks[s->above_block_idx[b]];

  290.         VP56RefDc *lb = &s->left_block[ff_vp56_b6to4[b]];

  291.         int count = 0;

  292.         int dc = 0;

  293.         int i;

  294. 

  295.         if (ref_frame == lb->ref_frame) {

  296.             dc += lb->dc_coeff;

  297.             count++;

  298.         }

  299.         if (ref_frame == ab->ref_frame) {

  300.             dc += ab->dc_coeff;

  301.             count++;

  302.         }

  303.         if (s->avctx->codec->id == AV_CODEC_ID_VP5)

  304.             for (i=0; i<2; i++)

  305.                 if (count < 2 && ref_frame == ab[-1+2*i].ref_frame) {

  306.                     dc += ab[-1+2*i].dc_coeff;

  307.                     count++;

  308.                 }

  309.         if (count == 0)

  310.             dc = s->prev_dc[ff_vp56_b2p[b]][ref_frame];

  311.         else if (count == 2)

  312.             dc /= 2;

  313. 

  314.         s->block_coeff[b][idx] += dc;

  315.         s->prev_dc[ff_vp56_b2p[b]][ref_frame] = s->block_coeff[b][idx];

  316.         ab->dc_coeff = s->block_coeff[b][idx];

  317.         ab->ref_frame = ref_frame;

  318.         lb->dc_coeff = s->block_coeff[b][idx];

  319.         lb->ref_frame = ref_frame;

  320.         s->block_coeff[b][idx] *= s->dequant_dc;

  321.     }

  322. }

  323. 

  324. static void vp56_deblock_filter(VP56Context *s, uint8_t *yuv,

  325.                                 ptrdiff_t stride, int dx, int dy)

  326. {

  327.     int t = ff_vp56_filter_threshold[s->quantizer];

  328.     if (dx)  s->vp56dsp.edge_filter_hor(yuv +         10-dx , stride, t);

  329.     if (dy)  s->vp56dsp.edge_filter_ver(yuv + stride*(10-dy), stride, t);

  330. }

  331. 

  332. static void vp56_mc(VP56Context *s, int b, int plane, uint8_t *src,

  333.                     ptrdiff_t stride, int x, int y)

  334. {

  335.     uint8_t *dst = s->frames[VP56_FRAME_CURRENT]->data[plane] + s->block_offset[b];

  336.     uint8_t *src_block;

  337.     int src_offset;

  338.     int overlap_offset = 0;

  339.     int mask = s->vp56_coord_div[b] - 1;

  340.     int deblock_filtering = s->deblock_filtering;

  341.     int dx;

  342.     int dy;

  343. 

  344.     if (s->avctx->skip_loop_filter >= AVDISCARD_ALL ||

  345.         (s->avctx->skip_loop_filter >= AVDISCARD_NONKEY

  346.          && !s->frames[VP56_FRAME_CURRENT]->key_frame))

  347.         deblock_filtering = 0;

  348. 

  349.     dx = s->mv[b].x / s->vp56_coord_div[b];

  350.     dy = s->mv[b].y / s->vp56_coord_div[b];

  351. 

  352.     if (b >= 4) {

  353.         x /= 2;

  354.         y /= 2;

  355.     }

  356.     x += dx - 2;

  357.     y += dy - 2;

  358. 

  359.     if (x<0 || x+12>=s->plane_width[plane] ||

  360.         y<0 || y+12>=s->plane_height[plane]) {

  361.         s->vdsp.emulated_edge_mc(s->edge_emu_buffer,

  362.                                  src + s->block_offset[b] + (dy-2)*stride + (dx-2),

  363.                                  stride, stride,

  364.                                  12, 12, x, y,

  365.                                  s->plane_width[plane],

  366.                                  s->plane_height[plane]);

  367.         src_block = s->edge_emu_buffer;

  368.         src_offset = 2 + 2*stride;

  369.     } else if (deblock_filtering) {

  370.         /* only need a 12x12 block, but there is no such dsp function, */

  371.         /* so copy a 16x12 block */

  372.         s->hdsp.put_pixels_tab[0][0](s->edge_emu_buffer,

  373.                                      src + s->block_offset[b] + (dy-2)*stride + (dx-2),

  374.                                      stride, 12);

  375.         src_block = s->edge_emu_buffer;

  376.         src_offset = 2 + 2*stride;

  377.     } else {

  378.         src_block = src;

  379.         src_offset = s->block_offset[b] + dy*stride + dx;

  380.     }

  381. 

  382.     if (deblock_filtering)

  383.         vp56_deblock_filter(s, src_block, stride, dx&7, dy&7);

  384. 

  385.     if (s->mv[b].x & mask)

  386.         overlap_offset += (s->mv[b].x > 0) ? 1 : -1;

  387.     if (s->mv[b].y & mask)

  388.         overlap_offset += (s->mv[b].y > 0) ? stride : -stride;

  389. 

  390.     if (overlap_offset) {

  391.         if (s->filter)

  392.             s->filter(s, dst, src_block, src_offset, src_offset+overlap_offset,

  393.                       stride, s->mv[b], mask, s->filter_selection, b<4);

  394.         else

  395.             s->vp3dsp.put_no_rnd_pixels_l2(dst, src_block+src_offset,

  396.                                            src_block+src_offset+overlap_offset,

  397.                                            stride, 8);

  398.     } else {

  399.         s->hdsp.put_pixels_tab[1][0](dst, src_block+src_offset, stride, 8);

  400.     }

  401. }

  402. 

  403. static int vp56_decode_mb(VP56Context *s, int row, int col, int is_alpha)

  404. {

  405.     AVFrame *frame_current, *frame_ref;

  406.     VP56mb mb_type;

  407.     VP56Frame ref_frame;

  408.     int b, ab, b_max, plane, off;

  409.     int ret;

  410. 

  411.     if (s->frames[VP56_FRAME_CURRENT]->key_frame)

  412.         mb_type = VP56_MB_INTRA;

  413.     else

  414.         mb_type = vp56_decode_mv(s, row, col);

  415.     ref_frame = ff_vp56_reference_frame[mb_type];

  416. 

  417.     ret = s->parse_coeff(s);

  418.     if (ret < 0)

  419.         return ret;

  420. 

  421.     vp56_add_predictors_dc(s, ref_frame);

  422. 

  423.     frame_current = s->frames[VP56_FRAME_CURRENT];

  424.     frame_ref = s->frames[ref_frame];

  425.     if (mb_type != VP56_MB_INTRA && !frame_ref->data[0])

  426.         return 0;

  427. 

  428.     ab = 6*is_alpha;

  429.     b_max = 6 - 2*is_alpha;

  430. 

  431.     switch (mb_type) {

  432.         case VP56_MB_INTRA:

  433.             for (b=0; b<b_max; b++) {

  434.                 plane = ff_vp56_b2p[b+ab];

  435.                 s->vp3dsp.idct_put(frame_current->data[plane] + s->block_offset[b],

  436.                                 s->stride[plane], s->block_coeff[b]);

  437.             }

  438.             break;

  439. 

  440.         case VP56_MB_INTER_NOVEC_PF:

  441.         case VP56_MB_INTER_NOVEC_GF:

  442.             for (b=0; b<b_max; b++) {

  443.                 plane = ff_vp56_b2p[b+ab];

  444.                 off = s->block_offset[b];

  445.                 s->hdsp.put_pixels_tab[1][0](frame_current->data[plane] + off,

  446.                                              frame_ref->data[plane] + off,

  447.                                              s->stride[plane], 8);

  448.                 s->vp3dsp.idct_add(frame_current->data[plane] + off,

  449.                                 s->stride[plane], s->block_coeff[b]);

  450.             }

  451.             break;

  452. 

  453.         case VP56_MB_INTER_DELTA_PF:

  454.         case VP56_MB_INTER_V1_PF:

  455.         case VP56_MB_INTER_V2_PF:

  456.         case VP56_MB_INTER_DELTA_GF:

  457.         case VP56_MB_INTER_4V:

  458.         case VP56_MB_INTER_V1_GF:

  459.         case VP56_MB_INTER_V2_GF:

  460.             for (b=0; b<b_max; b++) {

  461.                 int x_off = b==1 || b==3 ? 8 : 0;

  462.                 int y_off = b==2 || b==3 ? 8 : 0;

  463.                 plane = ff_vp56_b2p[b+ab];

  464.                 vp56_mc(s, b, plane, frame_ref->data[plane], s->stride[plane],

  465.                         16*col+x_off, 16*row+y_off);

  466.                 s->vp3dsp.idct_add(frame_current->data[plane] + s->block_offset[b],

  467.                                 s->stride[plane], s->block_coeff[b]);

  468.             }

  469.             break;

  470.     }

  471. 

  472.     if (is_alpha) {

  473.         s->block_coeff[4][0] = 0;

  474.         s->block_coeff[5][0] = 0;

  475.     }

  476.     return 0;

  477. }

  478. 

  479. static int vp56_conceal_mb(VP56Context *s, int row, int col, int is_alpha)

  480. {

  481.     AVFrame *frame_current, *frame_ref;

  482.     VP56mb mb_type;

  483.     VP56Frame ref_frame;

  484.     int b, ab, b_max, plane, off;

  485. 

  486.     if (s->frames[VP56_FRAME_CURRENT]->key_frame)

  487.         mb_type = VP56_MB_INTRA;

  488.     else

  489.         mb_type = vp56_conceal_mv(s, row, col);

  490.     ref_frame = ff_vp56_reference_frame[mb_type];

  491. 

  492.     frame_current = s->frames[VP56_FRAME_CURRENT];

  493.     frame_ref = s->frames[ref_frame];

  494.     if (mb_type != VP56_MB_INTRA && !frame_ref->data[0])

  495.         return 0;

  496. 

  497.     ab = 6*is_alpha;

  498.     b_max = 6 - 2*is_alpha;

  499. 

  500.     switch (mb_type) {

  501.         case VP56_MB_INTRA:

  502.             for (b=0; b<b_max; b++) {

  503.                 plane = ff_vp56_b2p[b+ab];

  504.                 s->vp3dsp.idct_put(frame_current->data[plane] + s->block_offset[b],

  505.                                 s->stride[plane], s->block_coeff[b]);

  506.             }

  507.             break;

  508. 

  509.         case VP56_MB_INTER_NOVEC_PF:

  510.         case VP56_MB_INTER_NOVEC_GF:

  511.             for (b=0; b<b_max; b++) {

  512.                 plane = ff_vp56_b2p[b+ab];

  513.                 off = s->block_offset[b];

  514.                 s->hdsp.put_pixels_tab[1][0](frame_current->data[plane] + off,

  515.                                              frame_ref->data[plane] + off,

  516.                                              s->stride[plane], 8);

  517.                 s->vp3dsp.idct_add(frame_current->data[plane] + off,

  518.                                 s->stride[plane], s->block_coeff[b]);

  519.             }

  520.             break;

  521.     }

  522. 

  523.     if (is_alpha) {

  524.         s->block_coeff[4][0] = 0;

  525.         s->block_coeff[5][0] = 0;

  526.     }

  527.     return 0;

  528. }

  529. 

  530. static int vp56_size_changed(VP56Context *s)

  531. {

  532.     AVCodecContext *avctx = s->avctx;

  533.     int stride = s->frames[VP56_FRAME_CURRENT]->linesize[0];

  534.     int i;

  535. 

  536.     s->plane_width[0]  = s->plane_width[3]  = avctx->coded_width;

  537.     s->plane_width[1]  = s->plane_width[2]  = avctx->coded_width/2;

  538.     s->plane_height[0] = s->plane_height[3] = avctx->coded_height;

  539.     s->plane_height[1] = s->plane_height[2] = avctx->coded_height/2;

  540. 

  541.     for (i=0; i<4; i++)

  542.         s->stride[i] = s->flip * s->frames[VP56_FRAME_CURRENT]->linesize[i];

  543. 

  544.     s->mb_width  = (avctx->coded_width +15) / 16;

  545.     s->mb_height = (avctx->coded_height+15) / 16;

  546. 

  547.     if (s->mb_width > 1000 || s->mb_height > 1000) {

  548.         ff_set_dimensions(avctx, 0, 0);

  549.         av_log(avctx, AV_LOG_ERROR, "picture too big\n");

  550.         return AVERROR_INVALIDDATA;

  551.     }

  552. 

  553.     av_reallocp_array(&s->above_blocks, 4*s->mb_width+6,

  554.                       sizeof(*s->above_blocks));

  555.     av_reallocp_array(&s->macroblocks, s->mb_width*s->mb_height,

  556.                       sizeof(*s->macroblocks));

  557.     av_free(s->edge_emu_buffer_alloc);

  558.     s->edge_emu_buffer_alloc = av_malloc(16*stride);

  559.     s->edge_emu_buffer = s->edge_emu_buffer_alloc;

  560.     if (!s->above_blocks || !s->macroblocks || !s->edge_emu_buffer_alloc)

  561.         return AVERROR(ENOMEM);

  562.     if (s->flip < 0)

  563.         s->edge_emu_buffer += 15 * stride;

  564. 

  565.     if (s->alpha_context)

  566.         return vp56_size_changed(s->alpha_context);

  567. 

  568.     return 0;

  569. }

  570. 

  571. static int ff_vp56_decode_mbs(AVCodecContext *avctx, void *, int, int);

  572. 

  573. int ff_vp56_decode_frame(AVCodecContext *avctx, void *data, int *got_frame,

  574.                          AVPacket *avpkt)

  575. {

  576.     const uint8_t *buf = avpkt->data;

  577.     VP56Context *s = avctx->priv_data;

  578.     AVFrame *const p = s->frames[VP56_FRAME_CURRENT];

  579.     int remaining_buf_size = avpkt->size;

  580.     int av_uninit(alpha_offset);

  581.     int i, res;

  582.     int ret;

  583. 

  584.     if (s->has_alpha) {

  585.         if (remaining_buf_size < 3)

  586.             return AVERROR_INVALIDDATA;

  587.         alpha_offset = bytestream_get_be24(&buf);

  588.         remaining_buf_size -= 3;

  589.         if (remaining_buf_size < alpha_offset)

  590.             return AVERROR_INVALIDDATA;

  591.     }

  592. 

  593.     res = s->parse_header(s, buf, remaining_buf_size);

  594.     if (res < 0)

  595.         return res;

  596. 

  597.     if (res == VP56_SIZE_CHANGE) {

  598.         for (i = 0; i < 4; i++) {

  599.             av_frame_unref(s->frames[i]);

  600.             if (s->alpha_context)

  601.                 av_frame_unref(s->alpha_context->frames[i]);

  602.         }

  603.     }

  604. 

  605.     ret = ff_get_buffer(avctx, p, AV_GET_BUFFER_FLAG_REF);

  606.     if (ret < 0)

  607.         return ret;

  608. 

  609.     if (avctx->pix_fmt == AV_PIX_FMT_YUVA420P) {

  610.         av_frame_unref(s->alpha_context->frames[VP56_FRAME_CURRENT]);

  611.         if ((ret = av_frame_ref(s->alpha_context->frames[VP56_FRAME_CURRENT], p)) < 0) {

  612.             av_frame_unref(p);

  613.             return ret;

  614.         }

  615.     }

  616. 

  617.     if (res == VP56_SIZE_CHANGE) {

  618.         if (vp56_size_changed(s)) {

  619.             av_frame_unref(p);

  620.             return AVERROR_INVALIDDATA;

  621.         }

  622.     }

  623. 

  624.     if (avctx->pix_fmt == AV_PIX_FMT_YUVA420P) {

  625.         int bak_w = avctx->width;

  626.         int bak_h = avctx->height;

  627.         int bak_cw = avctx->coded_width;

  628.         int bak_ch = avctx->coded_height;

  629.         buf += alpha_offset;

  630.         remaining_buf_size -= alpha_offset;

  631. 

  632.         res = s->alpha_context->parse_header(s->alpha_context, buf, remaining_buf_size);

  633.         if (res != 0) {

  634.             if(res==VP56_SIZE_CHANGE) {

  635.                 av_log(avctx, AV_LOG_ERROR, "Alpha reconfiguration\n");

  636.                 avctx->width  = bak_w;

  637.                 avctx->height = bak_h;

  638.                 avctx->coded_width  = bak_cw;

  639.                 avctx->coded_height = bak_ch;

  640.             }

  641.             av_frame_unref(p);

  642.             return AVERROR_INVALIDDATA;

  643.         }

  644.     }

  645. 

  646.     avctx->execute2(avctx, ff_vp56_decode_mbs, 0, 0, (avctx->pix_fmt == AV_PIX_FMT_YUVA420P) + 1);

  647. 

  648.     if ((res = av_frame_ref(data, p)) < 0)

  649.         return res;

  650.     *got_frame = 1;

  651. 

  652.     return avpkt->size;

  653. }

  654. 

  655. static int ff_vp56_decode_mbs(AVCodecContext *avctx, void *data,

  656.                               int jobnr, int threadnr)

  657. {

  658.     VP56Context *s0 = avctx->priv_data;

  659.     int is_alpha = (jobnr == 1);

  660.     VP56Context *s = is_alpha ? s0->alpha_context : s0;

  661.     AVFrame *const p = s->frames[VP56_FRAME_CURRENT];

  662.     int mb_row, mb_col, mb_row_flip, mb_offset = 0;

  663.     int block, y, uv;

  664.     ptrdiff_t stride_y, stride_uv;

  665.     int res;

  666.     int damaged = 0;

  667. 

  668.     if (p->key_frame) {

  669.         p->pict_type = AV_PICTURE_TYPE_I;

  670.         s->default_models_init(s);

  671.         for (block=0; block<s->mb_height*s->mb_width; block++)

  672.             s->macroblocks[block].type = VP56_MB_INTRA;

  673.     } else {

  674.         p->pict_type = AV_PICTURE_TYPE_P;

  675.         vp56_parse_mb_type_models(s);

  676.         s->parse_vector_models(s);

  677.         s->mb_type = VP56_MB_INTER_NOVEC_PF;

  678.     }

  679. 

  680.     if (s->parse_coeff_models(s))

  681.         goto next;

  682. 

  683.     memset(s->prev_dc, 0, sizeof(s->prev_dc));

  684.     s->prev_dc[1][VP56_FRAME_CURRENT] = 128;

  685.     s->prev_dc[2][VP56_FRAME_CURRENT] = 128;

  686. 

  687.     for (block=0; block < 4*s->mb_width+6; block++) {

  688.         s->above_blocks[block].ref_frame = VP56_FRAME_NONE;

  689.         s->above_blocks[block].dc_coeff = 0;

  690.         s->above_blocks[block].not_null_dc = 0;

  691.     }

  692.     s->above_blocks[2*s->mb_width + 2].ref_frame = VP56_FRAME_CURRENT;

  693.     s->above_blocks[3*s->mb_width + 4].ref_frame = VP56_FRAME_CURRENT;

  694. 

  695.     stride_y  = p->linesize[0];

  696.     stride_uv = p->linesize[1];

  697. 

  698.     if (s->flip < 0)

  699.         mb_offset = 7;

  700. 

  701.     /* main macroblocks loop */

  702.     for (mb_row=0; mb_row<s->mb_height; mb_row++) {

  703.         if (s->flip < 0)

  704.             mb_row_flip = s->mb_height - mb_row - 1;

  705.         else

  706.             mb_row_flip = mb_row;

  707. 

  708.         for (block=0; block<4; block++) {

  709.             s->left_block[block].ref_frame = VP56_FRAME_NONE;

  710.             s->left_block[block].dc_coeff = 0;

  711.             s->left_block[block].not_null_dc = 0;

  712.         }

  713.         memset(s->coeff_ctx, 0, sizeof(s->coeff_ctx));

  714.         memset(s->coeff_ctx_last, 24, sizeof(s->coeff_ctx_last));

  715. 

  716.         s->above_block_idx[0] = 1;

  717.         s->above_block_idx[1] = 2;

  718.         s->above_block_idx[2] = 1;

  719.         s->above_block_idx[3] = 2;

  720.         s->above_block_idx[4] = 2*s->mb_width + 2 + 1;

  721.         s->above_block_idx[5] = 3*s->mb_width + 4 + 1;

  722. 

  723.         s->block_offset[s->frbi] = (mb_row_flip*16 + mb_offset) * stride_y;

  724.         s->block_offset[s->srbi] = s->block_offset[s->frbi] + 8*stride_y;

  725.         s->block_offset[1] = s->block_offset[0] + 8;

  726.         s->block_offset[3] = s->block_offset[2] + 8;

  727.         s->block_offset[4] = (mb_row_flip*8 + mb_offset) * stride_uv;

  728.         s->block_offset[5] = s->block_offset[4];

  729. 

  730.         for (mb_col=0; mb_col<s->mb_width; mb_col++) {

  731.             if (!damaged) {

  732.                 int ret = vp56_decode_mb(s, mb_row, mb_col, is_alpha);

  733.                 if (ret < 0)

  734.                     damaged = 1;

  735.             }

  736.             if (damaged)

  737.                 vp56_conceal_mb(s, mb_row, mb_col, is_alpha);

  738. 

  739.             for (y=0; y<4; y++) {

  740.                 s->above_block_idx[y] += 2;

  741.                 s->block_offset[y] += 16;

  742.             }

  743. 

  744.             for (uv=4; uv<6; uv++) {

  745.                 s->above_block_idx[uv] += 1;

  746.                 s->block_offset[uv] += 8;

  747.             }

  748.         }

  749.     }

  750. 

  751. next:

  752.     if (p->key_frame || s->golden_frame) {

  753.         av_frame_unref(s->frames[VP56_FRAME_GOLDEN]);

  754.         if ((res = av_frame_ref(s->frames[VP56_FRAME_GOLDEN], p)) < 0)

  755.             return res;

  756.     }

  757. 

  758.     av_frame_unref(s->frames[VP56_FRAME_PREVIOUS]);

  759.     FFSWAP(AVFrame *, s->frames[VP56_FRAME_CURRENT],

  760.                       s->frames[VP56_FRAME_PREVIOUS]);

  761.     return 0;

  762. }

  763. 

  764. av_cold int ff_vp56_init(AVCodecContext *avctx, int flip, int has_alpha)

  765. {

  766.     VP56Context *s = avctx->priv_data;

  767.     return ff_vp56_init_context(avctx, s, flip, has_alpha);

  768. }

  769. 

  770. av_cold int ff_vp56_init_context(AVCodecContext *avctx, VP56Context *s,

  771.                                   int flip, int has_alpha)

  772. {

  773.     int i;

  774. 

  775.     s->avctx = avctx;

  776.     avctx->pix_fmt = has_alpha ? AV_PIX_FMT_YUVA420P : AV_PIX_FMT_YUV420P;

  777.     if (avctx->skip_alpha) avctx->pix_fmt = AV_PIX_FMT_YUV420P;

  778. 

  779.     ff_h264chroma_init(&s->h264chroma, 8);

  780.     ff_hpeldsp_init(&s->hdsp, avctx->flags);

  781.     ff_videodsp_init(&s->vdsp, 8);

  782.     ff_vp3dsp_init(&s->vp3dsp, avctx->flags);

  783.     ff_vp56dsp_init(&s->vp56dsp, avctx->codec->id);

  784.     for (i = 0; i < 64; i++) {

  785. #define TRANSPOSE(x) (((x) >> 3) | (((x) & 7) << 3))

  786.         s->idct_scantable[i] = TRANSPOSE(ff_zigzag_direct[i]);

  787. #undef TRANSPOSE

  788.     }

  789. 

  790.     for (i = 0; i < FF_ARRAY_ELEMS(s->frames); i++) {

  791.         s->frames[i] = av_frame_alloc();

  792.         if (!s->frames[i]) {

  793.             ff_vp56_free(avctx);

  794.             return AVERROR(ENOMEM);

  795.         }

  796.     }

  797.     s->edge_emu_buffer_alloc = NULL;

  798. 

  799.     s->above_blocks = NULL;

  800.     s->macroblocks = NULL;

  801.     s->quantizer = -1;

  802.     s->deblock_filtering = 1;

  803.     s->golden_frame = 0;

  804. 

  805.     s->filter = NULL;

  806. 

  807.     s->has_alpha = has_alpha;

  808. 

  809.     s->modelp = &s->model;

  810. 

  811.     if (flip) {

  812.         s->flip = -1;

  813.         s->frbi = 2;

  814.         s->srbi = 0;

  815.     } else {

  816.         s->flip = 1;

  817.         s->frbi = 0;

  818.         s->srbi = 2;

  819.     }

  820. 

  821.     return 0;

  822. }

  823. 

  824. av_cold int ff_vp56_free(AVCodecContext *avctx)

  825. {

  826.     VP56Context *s = avctx->priv_data;

  827.     return ff_vp56_free_context(s);

  828. }

  829. 

  830. av_cold int ff_vp56_free_context(VP56Context *s)

  831. {

  832.     int i;

  833. 

  834.     av_freep(&s->above_blocks);

  835.     av_freep(&s->macroblocks);

  836.     av_freep(&s->edge_emu_buffer_alloc);

  837. 

  838.     for (i = 0; i < FF_ARRAY_ELEMS(s->frames); i++)

  839.         av_frame_free(&s->frames[i]);

  840. 

  841.     return 0;

  842. }