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.
88307 Commits
32 Branches
380MB
Tree: 8548434337
Branches Tags
${ item.name }
Create branch ${ searchTerm }
from '8548434337'
${ noResults }
FFmpeg/libavcodec/vc1_block.c

2990 lines
117KB
Raw Blame History 

  1. /*

  2.  * VC-1 and WMV3 decoder

  3.  * Copyright (c) 2011 Mashiat Sarker Shakkhar

  4.  * Copyright (c) 2006-2007 Konstantin Shishkov

  5.  * Partly based on vc9.c (c) 2005 Anonymous, Alex Beregszaszi, Michael Niedermayer

  6.  *

  7.  * This file is part of FFmpeg.

  8.  *

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

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

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

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

  13.  *

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

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

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

  17.  * Lesser General Public License for more details.

  18.  *

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

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

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

  22.  */

  23. 

  24. /**

  25.  * @file

  26.  * VC-1 and WMV3 block decoding routines

  27.  */

  28. 

  29. #include "avcodec.h"

  30. #include "mpegutils.h"

  31. #include "mpegvideo.h"

  32. #include "msmpeg4data.h"

  33. #include "unary.h"

  34. #include "vc1.h"

  35. #include "vc1_pred.h"

  36. #include "vc1acdata.h"

  37. #include "vc1data.h"

  38. 

  39. #define MB_INTRA_VLC_BITS 9

  40. #define DC_VLC_BITS 9

  41. 

  42. // offset tables for interlaced picture MVDATA decoding

  43. static const uint8_t offset_table[2][9] = {

  44.     {  0,  1,  2,  4,  8, 16, 32,  64, 128 },

  45.     {  0,  1,  3,  7, 15, 31, 63, 127, 255 },

  46. };

  47. 

  48. /***********************************************************************/

  49. /**

  50.  * @name VC-1 Bitplane decoding

  51.  * @see 8.7, p56

  52.  * @{

  53.  */

  54. 

  55. 

  56. static inline void init_block_index(VC1Context *v)

  57. {

  58.     MpegEncContext *s = &v->s;

  59.     ff_init_block_index(s);

  60.     if (v->field_mode && !(v->second_field ^ v->tff)) {

  61.         s->dest[0] += s->current_picture_ptr->f->linesize[0];

  62.         s->dest[1] += s->current_picture_ptr->f->linesize[1];

  63.         s->dest[2] += s->current_picture_ptr->f->linesize[2];

  64.     }

  65. }

  66. 

  67. /** @} */ //Bitplane group

  68. 

  69. static void vc1_put_signed_blocks_clamped(VC1Context *v)

  70. {

  71.     MpegEncContext *s = &v->s;

  72.     int topleft_mb_pos, top_mb_pos;

  73.     int stride_y, fieldtx = 0;

  74.     int v_dist;

  75. 

  76.     /* The put pixels loop is always one MB row behind the decoding loop,

  77.      * because we can only put pixels when overlap filtering is done, and

  78.      * for filtering of the bottom edge of a MB, we need the next MB row

  79.      * present as well.

  80.      * Within the row, the put pixels loop is also one MB col behind the

  81.      * decoding loop. The reason for this is again, because for filtering

  82.      * of the right MB edge, we need the next MB present. */

  83.     if (!s->first_slice_line) {

  84.         if (s->mb_x) {

  85.             topleft_mb_pos = (s->mb_y - 1) * s->mb_stride + s->mb_x - 1;

  86.             if (v->fcm == ILACE_FRAME)

  87.                 fieldtx = v->fieldtx_plane[topleft_mb_pos];

  88.             stride_y       = s->linesize << fieldtx;

  89.             v_dist         = (16 - fieldtx) >> (fieldtx == 0);

  90.             s->idsp.put_signed_pixels_clamped(v->block[v->topleft_blk_idx][0],

  91.                                               s->dest[0] - 16 * s->linesize - 16,

  92.                                               stride_y);

  93.             s->idsp.put_signed_pixels_clamped(v->block[v->topleft_blk_idx][1],

  94.                                               s->dest[0] - 16 * s->linesize - 8,

  95.                                               stride_y);

  96.             s->idsp.put_signed_pixels_clamped(v->block[v->topleft_blk_idx][2],

  97.                                               s->dest[0] - v_dist * s->linesize - 16,

  98.                                               stride_y);

  99.             s->idsp.put_signed_pixels_clamped(v->block[v->topleft_blk_idx][3],

  100.                                               s->dest[0] - v_dist * s->linesize - 8,

  101.                                               stride_y);

  102.             if (!CONFIG_GRAY || !(s->avctx->flags & AV_CODEC_FLAG_GRAY)) {

  103.             s->idsp.put_signed_pixels_clamped(v->block[v->topleft_blk_idx][4],

  104.                                               s->dest[1] - 8 * s->uvlinesize - 8,

  105.                                               s->uvlinesize);

  106.             s->idsp.put_signed_pixels_clamped(v->block[v->topleft_blk_idx][5],

  107.                                               s->dest[2] - 8 * s->uvlinesize - 8,

  108.                                               s->uvlinesize);

  109.             }

  110.         }

  111.         if (s->mb_x == s->mb_width - 1) {

  112.             top_mb_pos = (s->mb_y - 1) * s->mb_stride + s->mb_x;

  113.             if (v->fcm == ILACE_FRAME)

  114.                 fieldtx = v->fieldtx_plane[top_mb_pos];

  115.             stride_y   = s->linesize << fieldtx;

  116.             v_dist     = fieldtx ? 15 : 8;

  117.             s->idsp.put_signed_pixels_clamped(v->block[v->top_blk_idx][0],

  118.                                               s->dest[0] - 16 * s->linesize,

  119.                                               stride_y);

  120.             s->idsp.put_signed_pixels_clamped(v->block[v->top_blk_idx][1],

  121.                                               s->dest[0] - 16 * s->linesize + 8,

  122.                                               stride_y);

  123.             s->idsp.put_signed_pixels_clamped(v->block[v->top_blk_idx][2],

  124.                                               s->dest[0] - v_dist * s->linesize,

  125.                                               stride_y);

  126.             s->idsp.put_signed_pixels_clamped(v->block[v->top_blk_idx][3],

  127.                                               s->dest[0] - v_dist * s->linesize + 8,

  128.                                               stride_y);

  129.             if (!CONFIG_GRAY || !(s->avctx->flags & AV_CODEC_FLAG_GRAY)) {

  130.             s->idsp.put_signed_pixels_clamped(v->block[v->top_blk_idx][4],

  131.                                               s->dest[1] - 8 * s->uvlinesize,

  132.                                               s->uvlinesize);

  133.             s->idsp.put_signed_pixels_clamped(v->block[v->top_blk_idx][5],

  134.                                               s->dest[2] - 8 * s->uvlinesize,

  135.                                               s->uvlinesize);

  136.             }

  137.         }

  138.     }

  139. 

  140. #define inc_blk_idx(idx) do { \

  141.         idx++; \

  142.         if (idx >= v->n_allocated_blks) \

  143.             idx = 0; \

  144.     } while (0)

  145. 

  146.     inc_blk_idx(v->topleft_blk_idx);

  147.     inc_blk_idx(v->top_blk_idx);

  148.     inc_blk_idx(v->left_blk_idx);

  149.     inc_blk_idx(v->cur_blk_idx);

  150. }

  151. 

  152. /***********************************************************************/

  153. /**

  154.  * @name VC-1 Block-level functions

  155.  * @see 7.1.4, p91 and 8.1.1.7, p(1)04

  156.  * @{

  157.  */

  158. 

  159. /**

  160.  * @def GET_MQUANT

  161.  * @brief Get macroblock-level quantizer scale

  162.  */

  163. #define GET_MQUANT()                                           \

  164.     if (v->dquantfrm) {                                        \

  165.         int edges = 0;                                         \

  166.         if (v->dqprofile == DQPROFILE_ALL_MBS) {               \

  167.             if (v->dqbilevel) {                                \

  168.                 mquant = (get_bits1(gb)) ? v->altpq : v->pq;   \

  169.             } else {                                           \

  170.                 mqdiff = get_bits(gb, 3);                      \

  171.                 if (mqdiff != 7)                               \

  172.                     mquant = v->pq + mqdiff;                   \

  173.                 else                                           \

  174.                     mquant = get_bits(gb, 5);                  \

  175.             }                                                  \

  176.         }                                                      \

  177.         if (v->dqprofile == DQPROFILE_SINGLE_EDGE)             \

  178.             edges = 1 << v->dqsbedge;                          \

  179.         else if (v->dqprofile == DQPROFILE_DOUBLE_EDGES)       \

  180.             edges = (3 << v->dqsbedge) % 15;                   \

  181.         else if (v->dqprofile == DQPROFILE_FOUR_EDGES)         \

  182.             edges = 15;                                        \

  183.         if ((edges&1) && !s->mb_x)                             \

  184.             mquant = v->altpq;                                 \

  185.         if ((edges&2) && s->first_slice_line)                  \

  186.             mquant = v->altpq;                                 \

  187.         if ((edges&4) && s->mb_x == (s->mb_width - 1))         \

  188.             mquant = v->altpq;                                 \

  189.         if ((edges&8) && s->mb_y == (s->mb_height - 1))        \

  190.             mquant = v->altpq;                                 \

  191.         if (!mquant || mquant > 31) {                          \

  192.             av_log(v->s.avctx, AV_LOG_ERROR,                   \

  193.                    "Overriding invalid mquant %d\n", mquant);  \

  194.             mquant = 1;                                        \

  195.         }                                                      \

  196.     }

  197. 

  198. /**

  199.  * @def GET_MVDATA(_dmv_x, _dmv_y)

  200.  * @brief Get MV differentials

  201.  * @see MVDATA decoding from 8.3.5.2, p(1)20

  202.  * @param _dmv_x Horizontal differential for decoded MV

  203.  * @param _dmv_y Vertical differential for decoded MV

  204.  */

  205. #define GET_MVDATA(_dmv_x, _dmv_y)                                      \

  206.     index = 1 + get_vlc2(gb, ff_vc1_mv_diff_vlc[s->mv_table_index].table, \

  207.                          VC1_MV_DIFF_VLC_BITS, 2);                      \

  208.     if (index > 36) {                                                   \

  209.         mb_has_coeffs = 1;                                              \

  210.         index -= 37;                                                    \

  211.     } else                                                              \

  212.         mb_has_coeffs = 0;                                              \

  213.     s->mb_intra = 0;                                                    \

  214.     if (!index) {                                                       \

  215.         _dmv_x = _dmv_y = 0;                                            \

  216.     } else if (index == 35) {                                           \

  217.         _dmv_x = get_bits(gb, v->k_x - 1 + s->quarter_sample);          \

  218.         _dmv_y = get_bits(gb, v->k_y - 1 + s->quarter_sample);          \

  219.     } else if (index == 36) {                                           \

  220.         _dmv_x = 0;                                                     \

  221.         _dmv_y = 0;                                                     \

  222.         s->mb_intra = 1;                                                \

  223.     } else {                                                            \

  224.         index1 = index % 6;                                             \

  225.         _dmv_x = offset_table[1][index1];                               \

  226.         val = size_table[index1] - (!s->quarter_sample && index1 == 5); \

  227.         if (val > 0) {                                                  \

  228.             val = get_bits(gb, val);                                    \

  229.             sign = 0 - (val & 1);                                       \

  230.             _dmv_x = (sign ^ ((val >> 1) + _dmv_x)) - sign;             \

  231.         }                                                               \

  232.                                                                         \

  233.         index1 = index / 6;                                             \

  234.         _dmv_y = offset_table[1][index1];                               \

  235.         val = size_table[index1] - (!s->quarter_sample && index1 == 5); \

  236.         if (val > 0) {                                                  \

  237.             val = get_bits(gb, val);                                    \

  238.             sign = 0 - (val & 1);                                       \

  239.             _dmv_y = (sign ^ ((val >> 1) + _dmv_y)) - sign;             \

  240.         }                                                               \

  241.     }

  242. 

  243. static av_always_inline void get_mvdata_interlaced(VC1Context *v, int *dmv_x,

  244.                                                    int *dmv_y, int *pred_flag)

  245. {

  246.     int index, index1;

  247.     int extend_x, extend_y;

  248.     GetBitContext *gb = &v->s.gb;

  249.     int bits, esc;

  250.     int val, sign;

  251. 

  252.     if (v->numref) {

  253.         bits = VC1_2REF_MVDATA_VLC_BITS;

  254.         esc  = 125;

  255.     } else {

  256.         bits = VC1_1REF_MVDATA_VLC_BITS;

  257.         esc  = 71;

  258.     }

  259.     extend_x = v->dmvrange & 1;

  260.     extend_y = (v->dmvrange >> 1) & 1;

  261.     index = get_vlc2(gb, v->imv_vlc->table, bits, 3);

  262.     if (index == esc) {

  263.         *dmv_x = get_bits(gb, v->k_x);

  264.         *dmv_y = get_bits(gb, v->k_y);

  265.         if (v->numref) {

  266.             if (pred_flag)

  267.                 *pred_flag = *dmv_y & 1;

  268.             *dmv_y = (*dmv_y + (*dmv_y & 1)) >> 1;

  269.         }

  270.     }

  271.     else {

  272.         av_assert0(index < esc);

  273.         index1 = (index + 1) % 9;

  274.         if (index1 != 0) {

  275.             val    = get_bits(gb, index1 + extend_x);

  276.             sign   = 0 - (val & 1);

  277.             *dmv_x = (sign ^ ((val >> 1) + offset_table[extend_x][index1])) - sign;

  278.         } else

  279.             *dmv_x = 0;

  280.         index1 = (index + 1) / 9;

  281.         if (index1 > v->numref) {

  282.             val    = get_bits(gb, (index1 >> v->numref) + extend_y);

  283.             sign   = 0 - (val & 1);

  284.             *dmv_y = (sign ^ ((val >> 1) + offset_table[extend_y][index1 >> v->numref])) - sign;

  285.         } else

  286.             *dmv_y = 0;

  287.         if (v->numref && pred_flag)

  288.             *pred_flag = index1 & 1;

  289.     }

  290. }

  291. 

  292. /** Reconstruct motion vector for B-frame and do motion compensation

  293.  */

  294. static inline void vc1_b_mc(VC1Context *v, int dmv_x[2], int dmv_y[2],

  295.                             int direct, int mode)

  296. {

  297.     if (direct) {

  298.         ff_vc1_mc_1mv(v, 0);

  299.         ff_vc1_interp_mc(v);

  300.         return;

  301.     }

  302.     if (mode == BMV_TYPE_INTERPOLATED) {

  303.         ff_vc1_mc_1mv(v, 0);

  304.         ff_vc1_interp_mc(v);

  305.         return;

  306.     }

  307. 

  308.     ff_vc1_mc_1mv(v, (mode == BMV_TYPE_BACKWARD));

  309. }

  310. 

  311. /** Get predicted DC value for I-frames only

  312.  * prediction dir: left=0, top=1

  313.  * @param s MpegEncContext

  314.  * @param overlap flag indicating that overlap filtering is used

  315.  * @param pq integer part of picture quantizer

  316.  * @param[in] n block index in the current MB

  317.  * @param dc_val_ptr Pointer to DC predictor

  318.  * @param dir_ptr Prediction direction for use in AC prediction

  319.  */

  320. static inline int vc1_i_pred_dc(MpegEncContext *s, int overlap, int pq, int n,

  321.                                 int16_t **dc_val_ptr, int *dir_ptr)

  322. {

  323.     int a, b, c, wrap, pred, scale;

  324.     int16_t *dc_val;

  325.     static const uint16_t dcpred[32] = {

  326.         -1, 1024,  512,  341,  256,  205,  171,  146,  128,

  327.              114,  102,   93,   85,   79,   73,   68,   64,

  328.               60,   57,   54,   51,   49,   47,   45,   43,

  329.               41,   39,   38,   37,   35,   34,   33

  330.     };

  331. 

  332.     /* find prediction - wmv3_dc_scale always used here in fact */

  333.     if (n < 4) scale = s->y_dc_scale;

  334.     else       scale = s->c_dc_scale;

  335. 

  336.     wrap   = s->block_wrap[n];

  337.     dc_val = s->dc_val[0] + s->block_index[n];

  338. 

  339.     /* B A

  340.      * C X

  341.      */

  342.     c = dc_val[ - 1];

  343.     b = dc_val[ - 1 - wrap];

  344.     a = dc_val[ - wrap];

  345. 

  346.     if (pq < 9 || !overlap) {

  347.         /* Set outer values */

  348.         if (s->first_slice_line && (n != 2 && n != 3))

  349.             b = a = dcpred[scale];

  350.         if (s->mb_x == 0 && (n != 1 && n != 3))

  351.             b = c = dcpred[scale];

  352.     } else {

  353.         /* Set outer values */

  354.         if (s->first_slice_line && (n != 2 && n != 3))

  355.             b = a = 0;

  356.         if (s->mb_x == 0 && (n != 1 && n != 3))

  357.             b = c = 0;

  358.     }

  359. 

  360.     if (abs(a - b) <= abs(b - c)) {

  361.         pred     = c;

  362.         *dir_ptr = 1; // left

  363.     } else {

  364.         pred     = a;

  365.         *dir_ptr = 0; // top

  366.     }

  367. 

  368.     /* update predictor */

  369.     *dc_val_ptr = &dc_val[0];

  370.     return pred;

  371. }

  372. 

  373. 

  374. /** Get predicted DC value

  375.  * prediction dir: left=0, top=1

  376.  * @param s MpegEncContext

  377.  * @param overlap flag indicating that overlap filtering is used

  378.  * @param pq integer part of picture quantizer

  379.  * @param[in] n block index in the current MB

  380.  * @param a_avail flag indicating top block availability

  381.  * @param c_avail flag indicating left block availability

  382.  * @param dc_val_ptr Pointer to DC predictor

  383.  * @param dir_ptr Prediction direction for use in AC prediction

  384.  */

  385. static inline int ff_vc1_pred_dc(MpegEncContext *s, int overlap, int pq, int n,

  386.                               int a_avail, int c_avail,

  387.                               int16_t **dc_val_ptr, int *dir_ptr)

  388. {

  389.     int a, b, c, wrap, pred;

  390.     int16_t *dc_val;

  391.     int mb_pos = s->mb_x + s->mb_y * s->mb_stride;

  392.     int q1, q2 = 0;

  393.     int dqscale_index;

  394. 

  395.     /* scale predictors if needed */

  396.     q1 = s->current_picture.qscale_table[mb_pos];

  397.     dqscale_index = s->y_dc_scale_table[q1] - 1;

  398.     if (dqscale_index < 0)

  399.         return 0;

  400. 

  401.     wrap = s->block_wrap[n];

  402.     dc_val = s->dc_val[0] + s->block_index[n];

  403. 

  404.     /* B A

  405.      * C X

  406.      */

  407.     c = dc_val[ - 1];

  408.     b = dc_val[ - 1 - wrap];

  409.     a = dc_val[ - wrap];

  410. 

  411.     if (c_avail && (n != 1 && n != 3)) {

  412.         q2 = s->current_picture.qscale_table[mb_pos - 1];

  413.         if (q2 && q2 != q1)

  414.             c = (c * s->y_dc_scale_table[q2] * ff_vc1_dqscale[dqscale_index] + 0x20000) >> 18;

  415.     }

  416.     if (a_avail && (n != 2 && n != 3)) {

  417.         q2 = s->current_picture.qscale_table[mb_pos - s->mb_stride];

  418.         if (q2 && q2 != q1)

  419.             a = (a * s->y_dc_scale_table[q2] * ff_vc1_dqscale[dqscale_index] + 0x20000) >> 18;

  420.     }

  421.     if (a_avail && c_avail && (n != 3)) {

  422.         int off = mb_pos;

  423.         if (n != 1)

  424.             off--;

  425.         if (n != 2)

  426.             off -= s->mb_stride;

  427.         q2 = s->current_picture.qscale_table[off];

  428.         if (q2 && q2 != q1)

  429.             b = (b * s->y_dc_scale_table[q2] * ff_vc1_dqscale[dqscale_index] + 0x20000) >> 18;

  430.     }

  431. 

  432.     if (c_avail && (!a_avail || abs(a - b) <= abs(b - c))) {

  433.         pred     = c;

  434.         *dir_ptr = 1; // left

  435.     } else if (a_avail) {

  436.         pred     = a;

  437.         *dir_ptr = 0; // top

  438.     } else {

  439.         pred     = 0;

  440.         *dir_ptr = 1; // left

  441.     }

  442. 

  443.     /* update predictor */

  444.     *dc_val_ptr = &dc_val[0];

  445.     return pred;

  446. }

  447. 

  448. /** @} */ // Block group

  449. 

  450. /**

  451.  * @name VC1 Macroblock-level functions in Simple/Main Profiles

  452.  * @see 7.1.4, p91 and 8.1.1.7, p(1)04

  453.  * @{

  454.  */

  455. 

  456. static inline int vc1_coded_block_pred(MpegEncContext * s, int n,

  457.                                        uint8_t **coded_block_ptr)

  458. {

  459.     int xy, wrap, pred, a, b, c;

  460. 

  461.     xy   = s->block_index[n];

  462.     wrap = s->b8_stride;

  463. 

  464.     /* B C

  465.      * A X

  466.      */

  467.     a = s->coded_block[xy - 1       ];

  468.     b = s->coded_block[xy - 1 - wrap];

  469.     c = s->coded_block[xy     - wrap];

  470. 

  471.     if (b == c) {

  472.         pred = a;

  473.     } else {

  474.         pred = c;

  475.     }

  476. 

  477.     /* store value */

  478.     *coded_block_ptr = &s->coded_block[xy];

  479. 

  480.     return pred;

  481. }

  482. 

  483. /**

  484.  * Decode one AC coefficient

  485.  * @param v The VC1 context

  486.  * @param last Last coefficient

  487.  * @param skip How much zero coefficients to skip

  488.  * @param value Decoded AC coefficient value

  489.  * @param codingset set of VLC to decode data

  490.  * @see 8.1.3.4

  491.  */

  492. static int vc1_decode_ac_coeff(VC1Context *v, int *last, int *skip,

  493.                                 int *value, int codingset)

  494. {

  495.     GetBitContext *gb = &v->s.gb;

  496.     int index, run, level, lst, sign;

  497. 

  498.     index = get_vlc2(gb, ff_vc1_ac_coeff_table[codingset].table, AC_VLC_BITS, 3);

  499.     if (index < 0)

  500.         return index;

  501.     if (index != ff_vc1_ac_sizes[codingset] - 1) {

  502.         run   = vc1_index_decode_table[codingset][index][0];

  503.         level = vc1_index_decode_table[codingset][index][1];

  504.         lst   = index >= vc1_last_decode_table[codingset] || get_bits_left(gb) < 0;

  505.         sign  = get_bits1(gb);

  506.     } else {

  507.         int escape = decode210(gb);

  508.         if (escape != 2) {

  509.             index = get_vlc2(gb, ff_vc1_ac_coeff_table[codingset].table, AC_VLC_BITS, 3);

  510.             run   = vc1_index_decode_table[codingset][index][0];

  511.             level = vc1_index_decode_table[codingset][index][1];

  512.             lst   = index >= vc1_last_decode_table[codingset];

  513.             if (escape == 0) {

  514.                 if (lst)

  515.                     level += vc1_last_delta_level_table[codingset][run];

  516.                 else

  517.                     level += vc1_delta_level_table[codingset][run];

  518.             } else {

  519.                 if (lst)

  520.                     run += vc1_last_delta_run_table[codingset][level] + 1;

  521.                 else

  522.                     run += vc1_delta_run_table[codingset][level] + 1;

  523.             }

  524.             sign = get_bits1(gb);

  525.         } else {

  526.             lst = get_bits1(gb);

  527.             if (v->s.esc3_level_length == 0) {

  528.                 if (v->pq < 8 || v->dquantfrm) { // table 59

  529.                     v->s.esc3_level_length = get_bits(gb, 3);

  530.                     if (!v->s.esc3_level_length)

  531.                         v->s.esc3_level_length = get_bits(gb, 2) + 8;

  532.                 } else { // table 60

  533.                     v->s.esc3_level_length = get_unary(gb, 1, 6) + 2;

  534.                 }

  535.                 v->s.esc3_run_length = 3 + get_bits(gb, 2);

  536.             }

  537.             run   = get_bits(gb, v->s.esc3_run_length);

  538.             sign  = get_bits1(gb);

  539.             level = get_bits(gb, v->s.esc3_level_length);

  540.         }

  541.     }

  542. 

  543.     *last  = lst;

  544.     *skip  = run;

  545.     *value = (level ^ -sign) + sign;

  546. 

  547.     return 0;

  548. }

  549. 

  550. /** Decode intra block in intra frames - should be faster than decode_intra_block

  551.  * @param v VC1Context

  552.  * @param block block to decode

  553.  * @param[in] n subblock index

  554.  * @param coded are AC coeffs present or not

  555.  * @param codingset set of VLC to decode data

  556.  */

  557. static int vc1_decode_i_block(VC1Context *v, int16_t block[64], int n,

  558.                               int coded, int codingset)

  559. {

  560.     GetBitContext *gb = &v->s.gb;

  561.     MpegEncContext *s = &v->s;

  562.     int dc_pred_dir = 0; /* Direction of the DC prediction used */

  563.     int i;

  564.     int16_t *dc_val;

  565.     int16_t *ac_val, *ac_val2;

  566.     int dcdiff, scale;

  567. 

  568.     /* Get DC differential */

  569.     if (n < 4) {

  570.         dcdiff = get_vlc2(&s->gb, ff_msmp4_dc_luma_vlc[s->dc_table_index].table, DC_VLC_BITS, 3);

  571.     } else {

  572.         dcdiff = get_vlc2(&s->gb, ff_msmp4_dc_chroma_vlc[s->dc_table_index].table, DC_VLC_BITS, 3);

  573.     }

  574.     if (dcdiff < 0) {

  575.         av_log(s->avctx, AV_LOG_ERROR, "Illegal DC VLC\n");

  576.         return -1;

  577.     }

  578.     if (dcdiff) {

  579.         const int m = (v->pq == 1 || v->pq == 2) ? 3 - v->pq : 0;

  580.         if (dcdiff == 119 /* ESC index value */) {

  581.             dcdiff = get_bits(gb, 8 + m);

  582.         } else {

  583.             if (m)

  584.                 dcdiff = (dcdiff << m) + get_bits(gb, m) - ((1 << m) - 1);

  585.         }

  586.         if (get_bits1(gb))

  587.             dcdiff = -dcdiff;

  588.     }

  589. 

  590.     /* Prediction */

  591.     dcdiff += vc1_i_pred_dc(&v->s, v->overlap, v->pq, n, &dc_val, &dc_pred_dir);

  592.     *dc_val = dcdiff;

  593. 

  594.     /* Store the quantized DC coeff, used for prediction */

  595.     if (n < 4)

  596.         scale = s->y_dc_scale;

  597.     else

  598.         scale = s->c_dc_scale;

  599.     block[0] = dcdiff * scale;

  600. 

  601.     ac_val  = s->ac_val[0][s->block_index[n]];

  602.     ac_val2 = ac_val;

  603.     if (dc_pred_dir) // left

  604.         ac_val -= 16;

  605.     else // top

  606.         ac_val -= 16 * s->block_wrap[n];

  607. 

  608.     scale = v->pq * 2 + v->halfpq;

  609. 

  610.     //AC Decoding

  611.     i = !!coded;

  612. 

  613.     if (coded) {

  614.         int last = 0, skip, value;

  615.         const uint8_t *zz_table;

  616.         int k;

  617. 

  618.         if (v->s.ac_pred) {

  619.             if (!dc_pred_dir)

  620.                 zz_table = v->zz_8x8[2];

  621.             else

  622.                 zz_table = v->zz_8x8[3];

  623.         } else

  624.             zz_table = v->zz_8x8[1];

  625. 

  626.         while (!last) {

  627.             int ret = vc1_decode_ac_coeff(v, &last, &skip, &value, codingset);

  628.             if (ret < 0)

  629.                 return ret;

  630.             i += skip;

  631.             if (i > 63)

  632.                 break;

  633.             block[zz_table[i++]] = value;

  634.         }

  635. 

  636.         /* apply AC prediction if needed */

  637.         if (s->ac_pred) {

  638.             int sh;

  639.             if (dc_pred_dir) { // left

  640.                 sh = v->left_blk_sh;

  641.             } else { // top

  642.                 sh = v->top_blk_sh;

  643.                 ac_val += 8;

  644.             }

  645.             for (k = 1; k < 8; k++)

  646.                 block[k << sh] += ac_val[k];

  647.         }

  648.         /* save AC coeffs for further prediction */

  649.         for (k = 1; k < 8; k++) {

  650.             ac_val2[k]     = block[k << v->left_blk_sh];

  651.             ac_val2[k + 8] = block[k << v->top_blk_sh];

  652.         }

  653. 

  654.         /* scale AC coeffs */

  655.         for (k = 1; k < 64; k++)

  656.             if (block[k]) {

  657.                 block[k] *= scale;

  658.                 if (!v->pquantizer)

  659.                     block[k] += (block[k] < 0) ? -v->pq : v->pq;

  660.             }

  661. 

  662.     } else {

  663.         int k;

  664. 

  665.         memset(ac_val2, 0, 16 * 2);

  666. 

  667.         /* apply AC prediction if needed */

  668.         if (s->ac_pred) {

  669.             int sh;

  670.             if (dc_pred_dir) { //left

  671.                 sh = v->left_blk_sh;

  672.             } else { // top

  673.                 sh = v->top_blk_sh;

  674.                 ac_val  += 8;

  675.                 ac_val2 += 8;

  676.             }

  677.             memcpy(ac_val2, ac_val, 8 * 2);

  678.             for (k = 1; k < 8; k++) {

  679.                 block[k << sh] = ac_val[k] * scale;

  680.                 if (!v->pquantizer && block[k << sh])

  681.                     block[k << sh] += (block[k << sh] < 0) ? -v->pq : v->pq;

  682.             }

  683.         }

  684.     }

  685.     if (s->ac_pred) i = 63;

  686.     s->block_last_index[n] = i;

  687. 

  688.     return 0;

  689. }

  690. 

  691. /** Decode intra block in intra frames - should be faster than decode_intra_block

  692.  * @param v VC1Context

  693.  * @param block block to decode

  694.  * @param[in] n subblock number

  695.  * @param coded are AC coeffs present or not

  696.  * @param codingset set of VLC to decode data

  697.  * @param mquant quantizer value for this macroblock

  698.  */

  699. static int vc1_decode_i_block_adv(VC1Context *v, int16_t block[64], int n,

  700.                                   int coded, int codingset, int mquant)

  701. {

  702.     GetBitContext *gb = &v->s.gb;

  703.     MpegEncContext *s = &v->s;

  704.     int dc_pred_dir = 0; /* Direction of the DC prediction used */

  705.     int i;

  706.     int16_t *dc_val = NULL;

  707.     int16_t *ac_val, *ac_val2;

  708.     int dcdiff;

  709.     int a_avail = v->a_avail, c_avail = v->c_avail;

  710.     int use_pred = s->ac_pred;

  711.     int scale;

  712.     int q1, q2 = 0;

  713.     int mb_pos = s->mb_x + s->mb_y * s->mb_stride;

  714. 

  715.     /* Get DC differential */

  716.     if (n < 4) {

  717.         dcdiff = get_vlc2(&s->gb, ff_msmp4_dc_luma_vlc[s->dc_table_index].table, DC_VLC_BITS, 3);

  718.     } else {

  719.         dcdiff = get_vlc2(&s->gb, ff_msmp4_dc_chroma_vlc[s->dc_table_index].table, DC_VLC_BITS, 3);

  720.     }

  721.     if (dcdiff < 0) {

  722.         av_log(s->avctx, AV_LOG_ERROR, "Illegal DC VLC\n");

  723.         return -1;

  724.     }

  725.     if (dcdiff) {

  726.         const int m = (mquant == 1 || mquant == 2) ? 3 - mquant : 0;

  727.         if (dcdiff == 119 /* ESC index value */) {

  728.             dcdiff = get_bits(gb, 8 + m);

  729.         } else {

  730.             if (m)

  731.                 dcdiff = (dcdiff << m) + get_bits(gb, m) - ((1 << m) - 1);

  732.         }

  733.         if (get_bits1(gb))

  734.             dcdiff = -dcdiff;

  735.     }

  736. 

  737.     /* Prediction */

  738.     dcdiff += ff_vc1_pred_dc(&v->s, v->overlap, mquant, n, v->a_avail, v->c_avail, &dc_val, &dc_pred_dir);

  739.     *dc_val = dcdiff;

  740. 

  741.     /* Store the quantized DC coeff, used for prediction */

  742.     if (n < 4)

  743.         scale = s->y_dc_scale;

  744.     else

  745.         scale = s->c_dc_scale;

  746.     block[0] = dcdiff * scale;

  747. 

  748.     /* check if AC is needed at all */

  749.     if (!a_avail && !c_avail)

  750.         use_pred = 0;

  751. 

  752.     scale = mquant * 2 + ((mquant == v->pq) ? v->halfpq : 0);

  753. 

  754.     ac_val  = s->ac_val[0][s->block_index[n]];

  755.     ac_val2 = ac_val;

  756.     if (dc_pred_dir) // left

  757.         ac_val -= 16;

  758.     else // top

  759.         ac_val -= 16 * s->block_wrap[n];

  760. 

  761.     q1 = s->current_picture.qscale_table[mb_pos];

  762.     if (n == 3)

  763.         q2 = q1;

  764.     else if (dc_pred_dir) {

  765.         if (n == 1)

  766.             q2 = q1;

  767.         else if (c_avail && mb_pos)

  768.             q2 = s->current_picture.qscale_table[mb_pos - 1];

  769.     } else {

  770.         if (n == 2)

  771.             q2 = q1;

  772.         else if (a_avail && mb_pos >= s->mb_stride)

  773.             q2 = s->current_picture.qscale_table[mb_pos - s->mb_stride];

  774.     }

  775. 

  776.     //AC Decoding

  777.     i = 1;

  778. 

  779.     if (coded) {

  780.         int last = 0, skip, value;

  781.         const uint8_t *zz_table;

  782.         int k;

  783. 

  784.         if (v->s.ac_pred) {

  785.             if (!use_pred && v->fcm == ILACE_FRAME) {

  786.                 zz_table = v->zzi_8x8;

  787.             } else {

  788.                 if (!dc_pred_dir) // top

  789.                     zz_table = v->zz_8x8[2];

  790.                 else // left

  791.                     zz_table = v->zz_8x8[3];

  792.             }

  793.         } else {

  794.             if (v->fcm != ILACE_FRAME)

  795.                 zz_table = v->zz_8x8[1];

  796.             else

  797.                 zz_table = v->zzi_8x8;

  798.         }

  799. 

  800.         while (!last) {

  801.             int ret = vc1_decode_ac_coeff(v, &last, &skip, &value, codingset);

  802.             if (ret < 0)

  803.                 return ret;

  804.             i += skip;

  805.             if (i > 63)

  806.                 break;

  807.             block[zz_table[i++]] = value;

  808.         }

  809. 

  810.         /* apply AC prediction if needed */

  811.         if (use_pred) {

  812.             int sh;

  813.             if (dc_pred_dir) { // left

  814.                 sh = v->left_blk_sh;

  815.             } else { // top

  816.                 sh = v->top_blk_sh;

  817.                 ac_val += 8;

  818.             }

  819.             /* scale predictors if needed*/

  820.             if (q2 && q1 != q2) {

  821.                 q1 = q1 * 2 + ((q1 == v->pq) ? v->halfpq : 0) - 1;

  822.                 if (q1 < 1)

  823.                     return AVERROR_INVALIDDATA;

  824.                 q2 = q2 * 2 + ((q2 == v->pq) ? v->halfpq : 0) - 1;

  825.                 for (k = 1; k < 8; k++)

  826.                     block[k << sh] += (ac_val[k] * q2 * ff_vc1_dqscale[q1 - 1] + 0x20000) >> 18;

  827.             } else {

  828.                 for (k = 1; k < 8; k++)

  829.                     block[k << sh] += ac_val[k];

  830.             }

  831.         }

  832.         /* save AC coeffs for further prediction */

  833.         for (k = 1; k < 8; k++) {

  834.             ac_val2[k    ] = block[k << v->left_blk_sh];

  835.             ac_val2[k + 8] = block[k << v->top_blk_sh];

  836.         }

  837. 

  838.         /* scale AC coeffs */

  839.         for (k = 1; k < 64; k++)

  840.             if (block[k]) {

  841.                 block[k] *= scale;

  842.                 if (!v->pquantizer)

  843.                     block[k] += (block[k] < 0) ? -mquant : mquant;

  844.             }

  845. 

  846.     } else { // no AC coeffs

  847.         int k;

  848. 

  849.         memset(ac_val2, 0, 16 * 2);

  850. 

  851.         /* apply AC prediction if needed */

  852.         if (use_pred) {

  853.             int sh;

  854.             if (dc_pred_dir) { // left

  855.                 sh = v->left_blk_sh;

  856.             } else { // top

  857.                 sh = v->top_blk_sh;

  858.                 ac_val  += 8;

  859.                 ac_val2 += 8;

  860.             }

  861.             memcpy(ac_val2, ac_val, 8 * 2);

  862.             if (q2 && q1 != q2) {

  863.                 q1 = q1 * 2 + ((q1 == v->pq) ? v->halfpq : 0) - 1;

  864.                 q2 = q2 * 2 + ((q2 == v->pq) ? v->halfpq : 0) - 1;

  865.                 if (q1 < 1)

  866.                     return AVERROR_INVALIDDATA;

  867.                 for (k = 1; k < 8; k++)

  868.                     ac_val2[k] = (ac_val2[k] * q2 * ff_vc1_dqscale[q1 - 1] + 0x20000) >> 18;

  869.             }

  870.             for (k = 1; k < 8; k++) {

  871.                 block[k << sh] = ac_val2[k] * scale;

  872.                 if (!v->pquantizer && block[k << sh])

  873.                     block[k << sh] += (block[k << sh] < 0) ? -mquant : mquant;

  874.             }

  875.         }

  876.     }

  877.     if (use_pred) i = 63;

  878.     s->block_last_index[n] = i;

  879. 

  880.     return 0;

  881. }

  882. 

  883. /** Decode intra block in inter frames - more generic version than vc1_decode_i_block

  884.  * @param v VC1Context

  885.  * @param block block to decode

  886.  * @param[in] n subblock index

  887.  * @param coded are AC coeffs present or not

  888.  * @param mquant block quantizer

  889.  * @param codingset set of VLC to decode data

  890.  */

  891. static int vc1_decode_intra_block(VC1Context *v, int16_t block[64], int n,

  892.                                   int coded, int mquant, int codingset)

  893. {

  894.     GetBitContext *gb = &v->s.gb;

  895.     MpegEncContext *s = &v->s;

  896.     int dc_pred_dir = 0; /* Direction of the DC prediction used */

  897.     int i;

  898.     int16_t *dc_val = NULL;

  899.     int16_t *ac_val, *ac_val2;

  900.     int dcdiff;

  901.     int mb_pos = s->mb_x + s->mb_y * s->mb_stride;

  902.     int a_avail = v->a_avail, c_avail = v->c_avail;

  903.     int use_pred = s->ac_pred;

  904.     int scale;

  905.     int q1, q2 = 0;

  906. 

  907.     s->bdsp.clear_block(block);

  908. 

  909.     /* XXX: Guard against dumb values of mquant */

  910.     mquant = av_clip_uintp2(mquant, 5);

  911. 

  912.     /* Set DC scale - y and c use the same */

  913.     s->y_dc_scale = s->y_dc_scale_table[mquant];

  914.     s->c_dc_scale = s->c_dc_scale_table[mquant];

  915. 

  916.     /* Get DC differential */

  917.     if (n < 4) {

  918.         dcdiff = get_vlc2(&s->gb, ff_msmp4_dc_luma_vlc[s->dc_table_index].table, DC_VLC_BITS, 3);

  919.     } else {

  920.         dcdiff = get_vlc2(&s->gb, ff_msmp4_dc_chroma_vlc[s->dc_table_index].table, DC_VLC_BITS, 3);

  921.     }

  922.     if (dcdiff < 0) {

  923.         av_log(s->avctx, AV_LOG_ERROR, "Illegal DC VLC\n");

  924.         return -1;

  925.     }

  926.     if (dcdiff) {

  927.         const int m = (mquant == 1 || mquant == 2) ? 3 - mquant : 0;

  928.         if (dcdiff == 119 /* ESC index value */) {

  929.             dcdiff = get_bits(gb, 8 + m);

  930.         } else {

  931.             if (m)

  932.                 dcdiff = (dcdiff << m) + get_bits(gb, m) - ((1 << m) - 1);

  933.         }

  934.         if (get_bits1(gb))

  935.             dcdiff = -dcdiff;

  936.     }

  937. 

  938.     /* Prediction */

  939.     dcdiff += ff_vc1_pred_dc(&v->s, v->overlap, mquant, n, a_avail, c_avail, &dc_val, &dc_pred_dir);

  940.     *dc_val = dcdiff;

  941. 

  942.     /* Store the quantized DC coeff, used for prediction */

  943. 

  944.     if (n < 4) {

  945.         block[0] = dcdiff * s->y_dc_scale;

  946.     } else {

  947.         block[0] = dcdiff * s->c_dc_scale;

  948.     }

  949. 

  950.     //AC Decoding

  951.     i = 1;

  952. 

  953.     /* check if AC is needed at all and adjust direction if needed */

  954.     if (!a_avail) dc_pred_dir = 1;

  955.     if (!c_avail) dc_pred_dir = 0;

  956.     if (!a_avail && !c_avail) use_pred = 0;

  957.     ac_val = s->ac_val[0][s->block_index[n]];

  958.     ac_val2 = ac_val;

  959. 

  960.     scale = mquant * 2 + v->halfpq;

  961. 

  962.     if (dc_pred_dir) //left

  963.         ac_val -= 16;

  964.     else //top

  965.         ac_val -= 16 * s->block_wrap[n];

  966. 

  967.     q1 = s->current_picture.qscale_table[mb_pos];

  968.     if (dc_pred_dir && c_avail && mb_pos)

  969.         q2 = s->current_picture.qscale_table[mb_pos - 1];

  970.     if (!dc_pred_dir && a_avail && mb_pos >= s->mb_stride)

  971.         q2 = s->current_picture.qscale_table[mb_pos - s->mb_stride];

  972.     if (dc_pred_dir && n == 1)

  973.         q2 = q1;

  974.     if (!dc_pred_dir && n == 2)

  975.         q2 = q1;

  976.     if (n == 3) q2 = q1;

  977. 

  978.     if (coded) {

  979.         int last = 0, skip, value;

  980.         int k;

  981. 

  982.         while (!last) {

  983.             int ret = vc1_decode_ac_coeff(v, &last, &skip, &value, codingset);

  984.             if (ret < 0)

  985.                 return ret;

  986.             i += skip;

  987.             if (i > 63)

  988.                 break;

  989.             if (v->fcm == PROGRESSIVE)

  990.                 block[v->zz_8x8[0][i++]] = value;

  991.             else {

  992.                 if (use_pred && (v->fcm == ILACE_FRAME)) {

  993.                     if (!dc_pred_dir) // top

  994.                         block[v->zz_8x8[2][i++]] = value;

  995.                     else // left

  996.                         block[v->zz_8x8[3][i++]] = value;

  997.                 } else {

  998.                     block[v->zzi_8x8[i++]] = value;

  999.                 }

  1000.             }

  1001.         }

  1002. 

  1003.         /* apply AC prediction if needed */

  1004.         if (use_pred) {

  1005.             /* scale predictors if needed*/

  1006.             if (q2 && q1 != q2) {

  1007.                 q1 = q1 * 2 + ((q1 == v->pq) ? v->halfpq : 0) - 1;

  1008.                 q2 = q2 * 2 + ((q2 == v->pq) ? v->halfpq : 0) - 1;

  1009. 

  1010.                 if (q1 < 1)

  1011.                     return AVERROR_INVALIDDATA;

  1012.                 if (dc_pred_dir) { // left

  1013.                     for (k = 1; k < 8; k++)

  1014.                         block[k << v->left_blk_sh] += (ac_val[k] * q2 * ff_vc1_dqscale[q1 - 1] + 0x20000) >> 18;

  1015.                 } else { //top

  1016.                     for (k = 1; k < 8; k++)

  1017.                         block[k << v->top_blk_sh] += (ac_val[k + 8] * q2 * ff_vc1_dqscale[q1 - 1] + 0x20000) >> 18;

  1018.                 }

  1019.             } else {

  1020.                 if (dc_pred_dir) { // left

  1021.                     for (k = 1; k < 8; k++)

  1022.                         block[k << v->left_blk_sh] += ac_val[k];

  1023.                 } else { // top

  1024.                     for (k = 1; k < 8; k++)

  1025.                         block[k << v->top_blk_sh] += ac_val[k + 8];

  1026.                 }

  1027.             }

  1028.         }

  1029.         /* save AC coeffs for further prediction */

  1030.         for (k = 1; k < 8; k++) {

  1031.             ac_val2[k    ] = block[k << v->left_blk_sh];

  1032.             ac_val2[k + 8] = block[k << v->top_blk_sh];

  1033.         }

  1034. 

  1035.         /* scale AC coeffs */

  1036.         for (k = 1; k < 64; k++)

  1037.             if (block[k]) {

  1038.                 block[k] *= scale;

  1039.                 if (!v->pquantizer)

  1040.                     block[k] += (block[k] < 0) ? -mquant : mquant;

  1041.             }

  1042. 

  1043.         if (use_pred) i = 63;

  1044.     } else { // no AC coeffs

  1045.         int k;

  1046. 

  1047.         memset(ac_val2, 0, 16 * 2);

  1048.         if (dc_pred_dir) { // left

  1049.             if (use_pred) {

  1050.                 memcpy(ac_val2, ac_val, 8 * 2);

  1051.                 if (q2 && q1 != q2) {

  1052.                     q1 = q1 * 2 + ((q1 == v->pq) ? v->halfpq : 0) - 1;

  1053.                     q2 = q2 * 2 + ((q2 == v->pq) ? v->halfpq : 0) - 1;

  1054.                     if (q1 < 1)

  1055.                         return AVERROR_INVALIDDATA;

  1056.                     for (k = 1; k < 8; k++)

  1057.                         ac_val2[k] = (ac_val2[k] * q2 * ff_vc1_dqscale[q1 - 1] + 0x20000) >> 18;

  1058.                 }

  1059.             }

  1060.         } else { // top

  1061.             if (use_pred) {

  1062.                 memcpy(ac_val2 + 8, ac_val + 8, 8 * 2);

  1063.                 if (q2 && q1 != q2) {

  1064.                     q1 = q1 * 2 + ((q1 == v->pq) ? v->halfpq : 0) - 1;

  1065.                     q2 = q2 * 2 + ((q2 == v->pq) ? v->halfpq : 0) - 1;

  1066.                     if (q1 < 1)

  1067.                         return AVERROR_INVALIDDATA;

  1068.                     for (k = 1; k < 8; k++)

  1069.                         ac_val2[k + 8] = (ac_val2[k + 8] * q2 * ff_vc1_dqscale[q1 - 1] + 0x20000) >> 18;

  1070.                 }

  1071.             }

  1072.         }

  1073. 

  1074.         /* apply AC prediction if needed */

  1075.         if (use_pred) {

  1076.             if (dc_pred_dir) { // left

  1077.                 for (k = 1; k < 8; k++) {

  1078.                     block[k << v->left_blk_sh] = ac_val2[k] * scale;

  1079.                     if (!v->pquantizer && block[k << v->left_blk_sh])

  1080.                         block[k << v->left_blk_sh] += (block[k << v->left_blk_sh] < 0) ? -mquant : mquant;

  1081.                 }

  1082.             } else { // top

  1083.                 for (k = 1; k < 8; k++) {

  1084.                     block[k << v->top_blk_sh] = ac_val2[k + 8] * scale;

  1085.                     if (!v->pquantizer && block[k << v->top_blk_sh])

  1086.                         block[k << v->top_blk_sh] += (block[k << v->top_blk_sh] < 0) ? -mquant : mquant;

  1087.                 }

  1088.             }

  1089.             i = 63;

  1090.         }

  1091.     }

  1092.     s->block_last_index[n] = i;

  1093. 

  1094.     return 0;

  1095. }

  1096. 

  1097. /** Decode P block

  1098.  */

  1099. static int vc1_decode_p_block(VC1Context *v, int16_t block[64], int n,

  1100.                               int mquant, int ttmb, int first_block,

  1101.                               uint8_t *dst, int linesize, int skip_block,

  1102.                               int *ttmb_out)

  1103. {

  1104.     MpegEncContext *s = &v->s;

  1105.     GetBitContext *gb = &s->gb;

  1106.     int i, j;

  1107.     int subblkpat = 0;

  1108.     int scale, off, idx, last, skip, value;

  1109.     int ttblk = ttmb & 7;

  1110.     int pat = 0;

  1111. 

  1112.     s->bdsp.clear_block(block);

  1113. 

  1114.     if (ttmb == -1) {

  1115.         ttblk = ff_vc1_ttblk_to_tt[v->tt_index][get_vlc2(gb, ff_vc1_ttblk_vlc[v->tt_index].table, VC1_TTBLK_VLC_BITS, 1)];

  1116.     }

  1117.     if (ttblk == TT_4X4) {

  1118.         subblkpat = ~(get_vlc2(gb, ff_vc1_subblkpat_vlc[v->tt_index].table, VC1_SUBBLKPAT_VLC_BITS, 1) + 1);

  1119.     }

  1120.     if ((ttblk != TT_8X8 && ttblk != TT_4X4)

  1121.         && ((v->ttmbf || (ttmb != -1 && (ttmb & 8) && !first_block))

  1122.             || (!v->res_rtm_flag && !first_block))) {

  1123.         subblkpat = decode012(gb);

  1124.         if (subblkpat)

  1125.             subblkpat ^= 3; // swap decoded pattern bits

  1126.         if (ttblk == TT_8X4_TOP || ttblk == TT_8X4_BOTTOM)

  1127.             ttblk = TT_8X4;

  1128.         if (ttblk == TT_4X8_RIGHT || ttblk == TT_4X8_LEFT)

  1129.             ttblk = TT_4X8;

  1130.     }

  1131.     scale = 2 * mquant + ((v->pq == mquant) ? v->halfpq : 0);

  1132. 

  1133.     // convert transforms like 8X4_TOP to generic TT and SUBBLKPAT

  1134.     if (ttblk == TT_8X4_TOP || ttblk == TT_8X4_BOTTOM) {

  1135.         subblkpat = 2 - (ttblk == TT_8X4_TOP);

  1136.         ttblk     = TT_8X4;

  1137.     }

  1138.     if (ttblk == TT_4X8_RIGHT || ttblk == TT_4X8_LEFT) {

  1139.         subblkpat = 2 - (ttblk == TT_4X8_LEFT);

  1140.         ttblk     = TT_4X8;

  1141.     }

  1142.     switch (ttblk) {

  1143.     case TT_8X8:

  1144.         pat  = 0xF;

  1145.         i    = 0;

  1146.         last = 0;

  1147.         while (!last) {

  1148.             int ret = vc1_decode_ac_coeff(v, &last, &skip, &value, v->codingset2);

  1149.             if (ret < 0)

  1150.                 return ret;

  1151.             i += skip;

  1152.             if (i > 63)

  1153.                 break;

  1154.             if (!v->fcm)

  1155.                 idx = v->zz_8x8[0][i++];

  1156.             else

  1157.                 idx = v->zzi_8x8[i++];

  1158.             block[idx] = value * scale;

  1159.             if (!v->pquantizer)

  1160.                 block[idx] += (block[idx] < 0) ? -mquant : mquant;

  1161.         }

  1162.         if (!skip_block) {

  1163.             if (i == 1)

  1164.                 v->vc1dsp.vc1_inv_trans_8x8_dc(dst, linesize, block);

  1165.             else {

  1166.                 v->vc1dsp.vc1_inv_trans_8x8(block);

  1167.                 s->idsp.add_pixels_clamped(block, dst, linesize);

  1168.             }

  1169.         }

  1170.         break;

  1171.     case TT_4X4:

  1172.         pat = ~subblkpat & 0xF;

  1173.         for (j = 0; j < 4; j++) {

  1174.             last = subblkpat & (1 << (3 - j));

  1175.             i    = 0;

  1176.             off  = (j & 1) * 4 + (j & 2) * 16;

  1177.             while (!last) {

  1178.                 int ret = vc1_decode_ac_coeff(v, &last, &skip, &value, v->codingset2);

  1179.                 if (ret < 0)

  1180.                     return ret;

  1181.                 i += skip;

  1182.                 if (i > 15)

  1183.                     break;

  1184.                 if (!v->fcm)

  1185.                     idx = ff_vc1_simple_progressive_4x4_zz[i++];

  1186.                 else

  1187.                     idx = ff_vc1_adv_interlaced_4x4_zz[i++];

  1188.                 block[idx + off] = value * scale;

  1189.                 if (!v->pquantizer)

  1190.                     block[idx + off] += (block[idx + off] < 0) ? -mquant : mquant;

  1191.             }

  1192.             if (!(subblkpat & (1 << (3 - j))) && !skip_block) {

  1193.                 if (i == 1)

  1194.                     v->vc1dsp.vc1_inv_trans_4x4_dc(dst + (j & 1) * 4 + (j & 2) * 2 * linesize, linesize, block + off);

  1195.                 else

  1196.                     v->vc1dsp.vc1_inv_trans_4x4(dst + (j & 1) * 4 + (j & 2) *  2 * linesize, linesize, block + off);

  1197.             }

  1198.         }

  1199.         break;

  1200.     case TT_8X4:

  1201.         pat = ~((subblkpat & 2) * 6 + (subblkpat & 1) * 3) & 0xF;

  1202.         for (j = 0; j < 2; j++) {

  1203.             last = subblkpat & (1 << (1 - j));

  1204.             i    = 0;

  1205.             off  = j * 32;

  1206.             while (!last) {

  1207.                 int ret = vc1_decode_ac_coeff(v, &last, &skip, &value, v->codingset2);

  1208.                 if (ret < 0)

  1209.                     return ret;

  1210.                 i += skip;

  1211.                 if (i > 31)

  1212.                     break;

  1213.                 if (!v->fcm)

  1214.                     idx = v->zz_8x4[i++] + off;

  1215.                 else

  1216.                     idx = ff_vc1_adv_interlaced_8x4_zz[i++] + off;

  1217.                 block[idx] = value * scale;

  1218.                 if (!v->pquantizer)

  1219.                     block[idx] += (block[idx] < 0) ? -mquant : mquant;

  1220.             }

  1221.             if (!(subblkpat & (1 << (1 - j))) && !skip_block) {

  1222.                 if (i == 1)

  1223.                     v->vc1dsp.vc1_inv_trans_8x4_dc(dst + j * 4 * linesize, linesize, block + off);

  1224.                 else

  1225.                     v->vc1dsp.vc1_inv_trans_8x4(dst + j * 4 * linesize, linesize, block + off);

  1226.             }

  1227.         }

  1228.         break;

  1229.     case TT_4X8:

  1230.         pat = ~(subblkpat * 5) & 0xF;

  1231.         for (j = 0; j < 2; j++) {

  1232.             last = subblkpat & (1 << (1 - j));

  1233.             i    = 0;

  1234.             off  = j * 4;

  1235.             while (!last) {

  1236.                 int ret = vc1_decode_ac_coeff(v, &last, &skip, &value, v->codingset2);

  1237.                 if (ret < 0)

  1238.                     return ret;

  1239.                 i += skip;

  1240.                 if (i > 31)

  1241.                     break;

  1242.                 if (!v->fcm)

  1243.                     idx = v->zz_4x8[i++] + off;

  1244.                 else

  1245.                     idx = ff_vc1_adv_interlaced_4x8_zz[i++] + off;

  1246.                 block[idx] = value * scale;

  1247.                 if (!v->pquantizer)

  1248.                     block[idx] += (block[idx] < 0) ? -mquant : mquant;

  1249.             }

  1250.             if (!(subblkpat & (1 << (1 - j))) && !skip_block) {

  1251.                 if (i == 1)

  1252.                     v->vc1dsp.vc1_inv_trans_4x8_dc(dst + j * 4, linesize, block + off);

  1253.                 else

  1254.                     v->vc1dsp.vc1_inv_trans_4x8(dst + j*4, linesize, block + off);

  1255.             }

  1256.         }

  1257.         break;

  1258.     }

  1259.     if (ttmb_out)

  1260.         *ttmb_out |= ttblk << (n * 4);

  1261.     return pat;

  1262. }

  1263. 

  1264. /** @} */ // Macroblock group

  1265. 

  1266. static const uint8_t size_table[6] = { 0, 2, 3, 4,  5,  8 };

  1267. 

  1268. /** Decode one P-frame MB

  1269.  */

  1270. static int vc1_decode_p_mb(VC1Context *v)

  1271. {

  1272.     MpegEncContext *s = &v->s;

  1273.     GetBitContext *gb = &s->gb;

  1274.     int i, j;

  1275.     int mb_pos = s->mb_x + s->mb_y * s->mb_stride;

  1276.     int cbp; /* cbp decoding stuff */

  1277.     int mqdiff, mquant; /* MB quantization */

  1278.     int ttmb = v->ttfrm; /* MB Transform type */

  1279. 

  1280.     int mb_has_coeffs = 1; /* last_flag */

  1281.     int dmv_x, dmv_y; /* Differential MV components */

  1282.     int index, index1; /* LUT indexes */

  1283.     int val, sign; /* temp values */

  1284.     int first_block = 1;

  1285.     int dst_idx, off;

  1286.     int skipped, fourmv;

  1287.     int block_cbp = 0, pat, block_tt = 0, block_intra = 0;

  1288. 

  1289.     mquant = v->pq; /* lossy initialization */

  1290. 

  1291.     if (v->mv_type_is_raw)

  1292.         fourmv = get_bits1(gb);

  1293.     else

  1294.         fourmv = v->mv_type_mb_plane[mb_pos];

  1295.     if (v->skip_is_raw)

  1296.         skipped = get_bits1(gb);

  1297.     else

  1298.         skipped = v->s.mbskip_table[mb_pos];

  1299. 

  1300.     if (!fourmv) { /* 1MV mode */

  1301.         if (!skipped) {

  1302.             GET_MVDATA(dmv_x, dmv_y);

  1303. 

  1304.             if (s->mb_intra) {

  1305.                 s->current_picture.motion_val[1][s->block_index[0]][0] = 0;

  1306.                 s->current_picture.motion_val[1][s->block_index[0]][1] = 0;

  1307.             }

  1308.             s->current_picture.mb_type[mb_pos] = s->mb_intra ? MB_TYPE_INTRA : MB_TYPE_16x16;

  1309.             ff_vc1_pred_mv(v, 0, dmv_x, dmv_y, 1, v->range_x, v->range_y, v->mb_type[0], 0, 0);

  1310. 

  1311.             /* FIXME Set DC val for inter block ? */

  1312.             if (s->mb_intra && !mb_has_coeffs) {

  1313.                 GET_MQUANT();

  1314.                 s->ac_pred = get_bits1(gb);

  1315.                 cbp        = 0;

  1316.             } else if (mb_has_coeffs) {

  1317.                 if (s->mb_intra)

  1318.                     s->ac_pred = get_bits1(gb);

  1319.                 cbp = get_vlc2(&v->s.gb, v->cbpcy_vlc->table, VC1_CBPCY_P_VLC_BITS, 2);

  1320.                 GET_MQUANT();

  1321.             } else {

  1322.                 mquant = v->pq;

  1323.                 cbp    = 0;

  1324.             }

  1325.             s->current_picture.qscale_table[mb_pos] = mquant;

  1326. 

  1327.             if (!v->ttmbf && !s->mb_intra && mb_has_coeffs)

  1328.                 ttmb = get_vlc2(gb, ff_vc1_ttmb_vlc[v->tt_index].table,

  1329.                                 VC1_TTMB_VLC_BITS, 2);

  1330.             if (!s->mb_intra) ff_vc1_mc_1mv(v, 0);

  1331.             dst_idx = 0;

  1332.             for (i = 0; i < 6; i++) {

  1333.                 s->dc_val[0][s->block_index[i]] = 0;

  1334.                 dst_idx += i >> 2;

  1335.                 val = ((cbp >> (5 - i)) & 1);

  1336.                 off = (i & 4) ? 0 : ((i & 1) * 8 + (i & 2) * 4 * s->linesize);

  1337.                 v->mb_type[0][s->block_index[i]] = s->mb_intra;

  1338.                 if (s->mb_intra) {

  1339.                     /* check if prediction blocks A and C are available */

  1340.                     v->a_avail = v->c_avail = 0;

  1341.                     if (i == 2 || i == 3 || !s->first_slice_line)

  1342.                         v->a_avail = v->mb_type[0][s->block_index[i] - s->block_wrap[i]];

  1343.                     if (i == 1 || i == 3 || s->mb_x)

  1344.                         v->c_avail = v->mb_type[0][s->block_index[i] - 1];

  1345. 

  1346.                     vc1_decode_intra_block(v, s->block[i], i, val, mquant,

  1347.                                            (i & 4) ? v->codingset2 : v->codingset);

  1348.                     if (CONFIG_GRAY && (i > 3) && (s->avctx->flags & AV_CODEC_FLAG_GRAY))

  1349.                         continue;

  1350.                     v->vc1dsp.vc1_inv_trans_8x8(s->block[i]);

  1351.                     if (v->rangeredfrm)

  1352.                         for (j = 0; j < 64; j++)

  1353.                             s->block[i][j] <<= 1;

  1354.                     s->idsp.put_signed_pixels_clamped(s->block[i],

  1355.                                                       s->dest[dst_idx] + off,

  1356.                                                       i & 4 ? s->uvlinesize

  1357.                                                             : s->linesize);

  1358.                     if (v->pq >= 9 && v->overlap) {

  1359.                         if (v->c_avail)

  1360.                             v->vc1dsp.vc1_h_overlap(s->dest[dst_idx] + off, i & 4 ? s->uvlinesize : s->linesize);

  1361.                         if (v->a_avail)

  1362.                             v->vc1dsp.vc1_v_overlap(s->dest[dst_idx] + off, i & 4 ? s->uvlinesize : s->linesize);

  1363.                     }

  1364.                     block_cbp   |= 0xF << (i << 2);

  1365.                     block_intra |= 1 << i;

  1366.                 } else if (val) {

  1367.                     pat = vc1_decode_p_block(v, s->block[i], i, mquant, ttmb, first_block,

  1368.                                              s->dest[dst_idx] + off, (i & 4) ? s->uvlinesize : s->linesize,

  1369.                                              CONFIG_GRAY && (i & 4) && (s->avctx->flags & AV_CODEC_FLAG_GRAY), &block_tt);

  1370.                     block_cbp |= pat << (i << 2);

  1371.                     if (!v->ttmbf && ttmb < 8)

  1372.                         ttmb = -1;

  1373.                     first_block = 0;

  1374.                 }

  1375.             }

  1376.         } else { // skipped

  1377.             s->mb_intra = 0;

  1378.             for (i = 0; i < 6; i++) {

  1379.                 v->mb_type[0][s->block_index[i]] = 0;

  1380.                 s->dc_val[0][s->block_index[i]]  = 0;

  1381.             }

  1382.             s->current_picture.mb_type[mb_pos]      = MB_TYPE_SKIP;

  1383.             s->current_picture.qscale_table[mb_pos] = 0;

  1384.             ff_vc1_pred_mv(v, 0, 0, 0, 1, v->range_x, v->range_y, v->mb_type[0], 0, 0);

  1385.             ff_vc1_mc_1mv(v, 0);

  1386.         }

  1387.     } else { // 4MV mode

  1388.         if (!skipped /* unskipped MB */) {

  1389.             int intra_count = 0, coded_inter = 0;

  1390.             int is_intra[6], is_coded[6];

  1391.             /* Get CBPCY */

  1392.             cbp = get_vlc2(&v->s.gb, v->cbpcy_vlc->table, VC1_CBPCY_P_VLC_BITS, 2);

  1393.             for (i = 0; i < 6; i++) {

  1394.                 val = ((cbp >> (5 - i)) & 1);

  1395.                 s->dc_val[0][s->block_index[i]] = 0;

  1396.                 s->mb_intra                     = 0;

  1397.                 if (i < 4) {

  1398.                     dmv_x = dmv_y = 0;

  1399.                     s->mb_intra   = 0;

  1400.                     mb_has_coeffs = 0;

  1401.                     if (val) {

  1402.                         GET_MVDATA(dmv_x, dmv_y);

  1403.                     }

  1404.                     ff_vc1_pred_mv(v, i, dmv_x, dmv_y, 0, v->range_x, v->range_y, v->mb_type[0], 0, 0);

  1405.                     if (!s->mb_intra)

  1406.                         ff_vc1_mc_4mv_luma(v, i, 0, 0);

  1407.                     intra_count += s->mb_intra;

  1408.                     is_intra[i]  = s->mb_intra;

  1409.                     is_coded[i]  = mb_has_coeffs;

  1410.                 }

  1411.                 if (i & 4) {

  1412.                     is_intra[i] = (intra_count >= 3);

  1413.                     is_coded[i] = val;

  1414.                 }

  1415.                 if (i == 4)

  1416.                     ff_vc1_mc_4mv_chroma(v, 0);

  1417.                 v->mb_type[0][s->block_index[i]] = is_intra[i];

  1418.                 if (!coded_inter)

  1419.                     coded_inter = !is_intra[i] & is_coded[i];

  1420.             }

  1421.             // if there are no coded blocks then don't do anything more

  1422.             dst_idx = 0;

  1423.             if (!intra_count && !coded_inter)

  1424.                 goto end;

  1425.             GET_MQUANT();

  1426.             s->current_picture.qscale_table[mb_pos] = mquant;

  1427.             /* test if block is intra and has pred */

  1428.             {

  1429.                 int intrapred = 0;

  1430.                 for (i = 0; i < 6; i++)

  1431.                     if (is_intra[i]) {

  1432.                         if (((!s->first_slice_line || (i == 2 || i == 3)) && v->mb_type[0][s->block_index[i] - s->block_wrap[i]])

  1433.                             || ((s->mb_x || (i == 1 || i == 3)) && v->mb_type[0][s->block_index[i] - 1])) {

  1434.                             intrapred = 1;

  1435.                             break;

  1436.                         }

  1437.                     }

  1438.                 if (intrapred)

  1439.                     s->ac_pred = get_bits1(gb);

  1440.                 else

  1441.                     s->ac_pred = 0;

  1442.             }

  1443.             if (!v->ttmbf && coded_inter)

  1444.                 ttmb = get_vlc2(gb, ff_vc1_ttmb_vlc[v->tt_index].table, VC1_TTMB_VLC_BITS, 2);

  1445.             for (i = 0; i < 6; i++) {

  1446.                 dst_idx    += i >> 2;

  1447.                 off         = (i & 4) ? 0 : ((i & 1) * 8 + (i & 2) * 4 * s->linesize);

  1448.                 s->mb_intra = is_intra[i];

  1449.                 if (is_intra[i]) {

  1450.                     /* check if prediction blocks A and C are available */

  1451.                     v->a_avail = v->c_avail = 0;

  1452.                     if (i == 2 || i == 3 || !s->first_slice_line)

  1453.                         v->a_avail = v->mb_type[0][s->block_index[i] - s->block_wrap[i]];

  1454.                     if (i == 1 || i == 3 || s->mb_x)

  1455.                         v->c_avail = v->mb_type[0][s->block_index[i] - 1];

  1456. 

  1457.                     vc1_decode_intra_block(v, s->block[i], i, is_coded[i], mquant,

  1458.                                            (i & 4) ? v->codingset2 : v->codingset);

  1459.                     if (CONFIG_GRAY && (i > 3) && (s->avctx->flags & AV_CODEC_FLAG_GRAY))

  1460.                         continue;

  1461.                     v->vc1dsp.vc1_inv_trans_8x8(s->block[i]);

  1462.                     if (v->rangeredfrm)

  1463.                         for (j = 0; j < 64; j++)

  1464.                             s->block[i][j] <<= 1;

  1465.                     s->idsp.put_signed_pixels_clamped(s->block[i],

  1466.                                                       s->dest[dst_idx] + off,

  1467.                                                       (i & 4) ? s->uvlinesize

  1468.                                                               : s->linesize);

  1469.                     if (v->pq >= 9 && v->overlap) {

  1470.                         if (v->c_avail)

  1471.                             v->vc1dsp.vc1_h_overlap(s->dest[dst_idx] + off, i & 4 ? s->uvlinesize : s->linesize);

  1472.                         if (v->a_avail)

  1473.                             v->vc1dsp.vc1_v_overlap(s->dest[dst_idx] + off, i & 4 ? s->uvlinesize : s->linesize);

  1474.                     }

  1475.                     block_cbp   |= 0xF << (i << 2);

  1476.                     block_intra |= 1 << i;

  1477.                 } else if (is_coded[i]) {

  1478.                     pat = vc1_decode_p_block(v, s->block[i], i, mquant, ttmb,

  1479.                                              first_block, s->dest[dst_idx] + off,

  1480.                                              (i & 4) ? s->uvlinesize : s->linesize,

  1481.                                              CONFIG_GRAY && (i & 4) && (s->avctx->flags & AV_CODEC_FLAG_GRAY),

  1482.                                              &block_tt);

  1483.                     block_cbp |= pat << (i << 2);

  1484.                     if (!v->ttmbf && ttmb < 8)

  1485.                         ttmb = -1;

  1486.                     first_block = 0;

  1487.                 }

  1488.             }

  1489.         } else { // skipped MB

  1490.             s->mb_intra                               = 0;

  1491.             s->current_picture.qscale_table[mb_pos] = 0;

  1492.             for (i = 0; i < 6; i++) {

  1493.                 v->mb_type[0][s->block_index[i]] = 0;

  1494.                 s->dc_val[0][s->block_index[i]]  = 0;

  1495.             }

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

  1497.                 ff_vc1_pred_mv(v, i, 0, 0, 0, v->range_x, v->range_y, v->mb_type[0], 0, 0);

  1498.                 ff_vc1_mc_4mv_luma(v, i, 0, 0);

  1499.             }

  1500.             ff_vc1_mc_4mv_chroma(v, 0);

  1501.             s->current_picture.qscale_table[mb_pos] = 0;

  1502.         }

  1503.     }

  1504. end:

  1505.     v->cbp[s->mb_x]      = block_cbp;

  1506.     v->ttblk[s->mb_x]    = block_tt;

  1507.     v->is_intra[s->mb_x] = block_intra;

  1508. 

  1509.     return 0;

  1510. }

  1511. 

  1512. /* Decode one macroblock in an interlaced frame p picture */

  1513. 

  1514. static int vc1_decode_p_mb_intfr(VC1Context *v)

  1515. {

  1516.     MpegEncContext *s = &v->s;

  1517.     GetBitContext *gb = &s->gb;

  1518.     int i;

  1519.     int mb_pos = s->mb_x + s->mb_y * s->mb_stride;

  1520.     int cbp = 0; /* cbp decoding stuff */

  1521.     int mqdiff, mquant; /* MB quantization */

  1522.     int ttmb = v->ttfrm; /* MB Transform type */

  1523. 

  1524.     int mb_has_coeffs = 1; /* last_flag */

  1525.     int dmv_x, dmv_y; /* Differential MV components */

  1526.     int val; /* temp value */

  1527.     int first_block = 1;

  1528.     int dst_idx, off;

  1529.     int skipped, fourmv = 0, twomv = 0;

  1530.     int block_cbp = 0, pat, block_tt = 0;

  1531.     int idx_mbmode = 0, mvbp;

  1532.     int stride_y, fieldtx;

  1533. 

  1534.     mquant = v->pq; /* Lossy initialization */

  1535. 

  1536.     if (v->skip_is_raw)

  1537.         skipped = get_bits1(gb);

  1538.     else

  1539.         skipped = v->s.mbskip_table[mb_pos];

  1540.     if (!skipped) {

  1541.         if (v->fourmvswitch)

  1542.             idx_mbmode = get_vlc2(gb, v->mbmode_vlc->table, VC1_INTFR_4MV_MBMODE_VLC_BITS, 2); // try getting this done

  1543.         else

  1544.             idx_mbmode = get_vlc2(gb, v->mbmode_vlc->table, VC1_INTFR_NON4MV_MBMODE_VLC_BITS, 2); // in a single line

  1545.         switch (ff_vc1_mbmode_intfrp[v->fourmvswitch][idx_mbmode][0]) {

  1546.         /* store the motion vector type in a flag (useful later) */

  1547.         case MV_PMODE_INTFR_4MV:

  1548.             fourmv = 1;

  1549.             v->blk_mv_type[s->block_index[0]] = 0;

  1550.             v->blk_mv_type[s->block_index[1]] = 0;

  1551.             v->blk_mv_type[s->block_index[2]] = 0;

  1552.             v->blk_mv_type[s->block_index[3]] = 0;

  1553.             break;

  1554.         case MV_PMODE_INTFR_4MV_FIELD:

  1555.             fourmv = 1;

  1556.             v->blk_mv_type[s->block_index[0]] = 1;

  1557.             v->blk_mv_type[s->block_index[1]] = 1;

  1558.             v->blk_mv_type[s->block_index[2]] = 1;

  1559.             v->blk_mv_type[s->block_index[3]] = 1;

  1560.             break;

  1561.         case MV_PMODE_INTFR_2MV_FIELD:

  1562.             twomv = 1;

  1563.             v->blk_mv_type[s->block_index[0]] = 1;

  1564.             v->blk_mv_type[s->block_index[1]] = 1;

  1565.             v->blk_mv_type[s->block_index[2]] = 1;

  1566.             v->blk_mv_type[s->block_index[3]] = 1;

  1567.             break;

  1568.         case MV_PMODE_INTFR_1MV:

  1569.             v->blk_mv_type[s->block_index[0]] = 0;

  1570.             v->blk_mv_type[s->block_index[1]] = 0;

  1571.             v->blk_mv_type[s->block_index[2]] = 0;

  1572.             v->blk_mv_type[s->block_index[3]] = 0;

  1573.             break;

  1574.         }

  1575.         if (ff_vc1_mbmode_intfrp[v->fourmvswitch][idx_mbmode][0] == MV_PMODE_INTFR_INTRA) { // intra MB

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

  1577.                 s->current_picture.motion_val[1][s->block_index[i]][0] = 0;

  1578.                 s->current_picture.motion_val[1][s->block_index[i]][1] = 0;

  1579.             }

  1580.             v->is_intra[s->mb_x] = 0x3f; // Set the bitfield to all 1.

  1581.             s->mb_intra          = 1;

  1582.             s->current_picture.mb_type[mb_pos] = MB_TYPE_INTRA;

  1583.             fieldtx = v->fieldtx_plane[mb_pos] = get_bits1(gb);

  1584.             mb_has_coeffs = get_bits1(gb);

  1585.             if (mb_has_coeffs)

  1586.                 cbp = 1 + get_vlc2(&v->s.gb, v->cbpcy_vlc->table, VC1_CBPCY_P_VLC_BITS, 2);

  1587.             v->s.ac_pred = v->acpred_plane[mb_pos] = get_bits1(gb);

  1588.             GET_MQUANT();

  1589.             s->current_picture.qscale_table[mb_pos] = mquant;

  1590.             /* Set DC scale - y and c use the same (not sure if necessary here) */

  1591.             s->y_dc_scale = s->y_dc_scale_table[mquant];

  1592.             s->c_dc_scale = s->c_dc_scale_table[mquant];

  1593.             dst_idx = 0;

  1594.             for (i = 0; i < 6; i++) {

  1595.                 v->a_avail = v->c_avail          = 0;

  1596.                 v->mb_type[0][s->block_index[i]] = 1;

  1597.                 s->dc_val[0][s->block_index[i]]  = 0;

  1598.                 dst_idx += i >> 2;

  1599.                 val = ((cbp >> (5 - i)) & 1);

  1600.                 if (i == 2 || i == 3 || !s->first_slice_line)

  1601.                     v->a_avail = v->mb_type[0][s->block_index[i] - s->block_wrap[i]];

  1602.                 if (i == 1 || i == 3 || s->mb_x)

  1603.                     v->c_avail = v->mb_type[0][s->block_index[i] - 1];

  1604. 

  1605.                 vc1_decode_intra_block(v, s->block[i], i, val, mquant,

  1606.                                        (i & 4) ? v->codingset2 : v->codingset);

  1607.                 if (CONFIG_GRAY && (i > 3) && (s->avctx->flags & AV_CODEC_FLAG_GRAY))

  1608.                     continue;

  1609.                 v->vc1dsp.vc1_inv_trans_8x8(s->block[i]);

  1610.                 if (i < 4) {

  1611.                     stride_y = s->linesize << fieldtx;

  1612.                     off = (fieldtx) ? ((i & 1) * 8) + ((i & 2) >> 1) * s->linesize : (i & 1) * 8 + 4 * (i & 2) * s->linesize;

  1613.                 } else {

  1614.                     stride_y = s->uvlinesize;

  1615.                     off = 0;

  1616.                 }

  1617.                 s->idsp.put_signed_pixels_clamped(s->block[i],

  1618.                                                   s->dest[dst_idx] + off,

  1619.                                                   stride_y);

  1620.                 //TODO: loop filter

  1621.             }

  1622. 

  1623.         } else { // inter MB

  1624.             mb_has_coeffs = ff_vc1_mbmode_intfrp[v->fourmvswitch][idx_mbmode][3];

  1625.             if (mb_has_coeffs)

  1626.                 cbp = 1 + get_vlc2(&v->s.gb, v->cbpcy_vlc->table, VC1_CBPCY_P_VLC_BITS, 2);

  1627.             if (ff_vc1_mbmode_intfrp[v->fourmvswitch][idx_mbmode][0] == MV_PMODE_INTFR_2MV_FIELD) {

  1628.                 v->twomvbp = get_vlc2(gb, v->twomvbp_vlc->table, VC1_2MV_BLOCK_PATTERN_VLC_BITS, 1);

  1629.             } else {

  1630.                 if ((ff_vc1_mbmode_intfrp[v->fourmvswitch][idx_mbmode][0] == MV_PMODE_INTFR_4MV)

  1631.                     || (ff_vc1_mbmode_intfrp[v->fourmvswitch][idx_mbmode][0] == MV_PMODE_INTFR_4MV_FIELD)) {

  1632.                     v->fourmvbp = get_vlc2(gb, v->fourmvbp_vlc->table, VC1_4MV_BLOCK_PATTERN_VLC_BITS, 1);

  1633.                 }

  1634.             }

  1635.             s->mb_intra = v->is_intra[s->mb_x] = 0;

  1636.             for (i = 0; i < 6; i++)

  1637.                 v->mb_type[0][s->block_index[i]] = 0;

  1638.             fieldtx = v->fieldtx_plane[mb_pos] = ff_vc1_mbmode_intfrp[v->fourmvswitch][idx_mbmode][1];

  1639.             /* for all motion vector read MVDATA and motion compensate each block */

  1640.             dst_idx = 0;

  1641.             if (fourmv) {

  1642.                 mvbp = v->fourmvbp;

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

  1644.                     dmv_x = dmv_y = 0;

  1645.                     if (mvbp & (8 >> i))

  1646.                         get_mvdata_interlaced(v, &dmv_x, &dmv_y, 0);

  1647.                     ff_vc1_pred_mv_intfr(v, i, dmv_x, dmv_y, 0, v->range_x, v->range_y, v->mb_type[0], 0);

  1648.                     ff_vc1_mc_4mv_luma(v, i, 0, 0);

  1649.                 }

  1650.                 ff_vc1_mc_4mv_chroma4(v, 0, 0, 0);

  1651.             } else if (twomv) {

  1652.                 mvbp  = v->twomvbp;

  1653.                 dmv_x = dmv_y = 0;

  1654.                 if (mvbp & 2) {

  1655.                     get_mvdata_interlaced(v, &dmv_x, &dmv_y, 0);

  1656.                 }

  1657.                 ff_vc1_pred_mv_intfr(v, 0, dmv_x, dmv_y, 2, v->range_x, v->range_y, v->mb_type[0], 0);

  1658.                 ff_vc1_mc_4mv_luma(v, 0, 0, 0);

  1659.                 ff_vc1_mc_4mv_luma(v, 1, 0, 0);

  1660.                 dmv_x = dmv_y = 0;

  1661.                 if (mvbp & 1) {

  1662.                     get_mvdata_interlaced(v, &dmv_x, &dmv_y, 0);

  1663.                 }

  1664.                 ff_vc1_pred_mv_intfr(v, 2, dmv_x, dmv_y, 2, v->range_x, v->range_y, v->mb_type[0], 0);

  1665.                 ff_vc1_mc_4mv_luma(v, 2, 0, 0);

  1666.                 ff_vc1_mc_4mv_luma(v, 3, 0, 0);

  1667.                 ff_vc1_mc_4mv_chroma4(v, 0, 0, 0);

  1668.             } else {

  1669.                 mvbp = ff_vc1_mbmode_intfrp[v->fourmvswitch][idx_mbmode][2];

  1670.                 dmv_x = dmv_y = 0;

  1671.                 if (mvbp) {

  1672.                     get_mvdata_interlaced(v, &dmv_x, &dmv_y, 0);

  1673.                 }

  1674.                 ff_vc1_pred_mv_intfr(v, 0, dmv_x, dmv_y, 1, v->range_x, v->range_y, v->mb_type[0], 0);

  1675.                 ff_vc1_mc_1mv(v, 0);

  1676.             }

  1677.             if (cbp)

  1678.                 GET_MQUANT();  // p. 227

  1679.             s->current_picture.qscale_table[mb_pos] = mquant;

  1680.             if (!v->ttmbf && cbp)

  1681.                 ttmb = get_vlc2(gb, ff_vc1_ttmb_vlc[v->tt_index].table, VC1_TTMB_VLC_BITS, 2);

  1682.             for (i = 0; i < 6; i++) {

  1683.                 s->dc_val[0][s->block_index[i]] = 0;

  1684.                 dst_idx += i >> 2;

  1685.                 val = ((cbp >> (5 - i)) & 1);

  1686.                 if (!fieldtx)

  1687.                     off = (i & 4) ? 0 : ((i & 1) * 8 + (i & 2) * 4 * s->linesize);

  1688.                 else

  1689.                     off = (i & 4) ? 0 : ((i & 1) * 8 + ((i > 1) * s->linesize));

  1690.                 if (val) {

  1691.                     pat = vc1_decode_p_block(v, s->block[i], i, mquant, ttmb,

  1692.                                              first_block, s->dest[dst_idx] + off,

  1693.                                              (i & 4) ? s->uvlinesize : (s->linesize << fieldtx),

  1694.                                              CONFIG_GRAY && (i & 4) && (s->avctx->flags & AV_CODEC_FLAG_GRAY), &block_tt);

  1695.                     block_cbp |= pat << (i << 2);

  1696.                     if (!v->ttmbf && ttmb < 8)

  1697.                         ttmb = -1;

  1698.                     first_block = 0;

  1699.                 }

  1700.             }

  1701.         }

  1702.     } else { // skipped

  1703.         s->mb_intra = v->is_intra[s->mb_x] = 0;

  1704.         for (i = 0; i < 6; i++) {

  1705.             v->mb_type[0][s->block_index[i]] = 0;

  1706.             s->dc_val[0][s->block_index[i]] = 0;

  1707.         }

  1708.         s->current_picture.mb_type[mb_pos]      = MB_TYPE_SKIP;

  1709.         s->current_picture.qscale_table[mb_pos] = 0;

  1710.         v->blk_mv_type[s->block_index[0]] = 0;

  1711.         v->blk_mv_type[s->block_index[1]] = 0;

  1712.         v->blk_mv_type[s->block_index[2]] = 0;

  1713.         v->blk_mv_type[s->block_index[3]] = 0;

  1714.         ff_vc1_pred_mv_intfr(v, 0, 0, 0, 1, v->range_x, v->range_y, v->mb_type[0], 0);

  1715.         ff_vc1_mc_1mv(v, 0);

  1716.     }

  1717.     if (s->mb_x == s->mb_width - 1)

  1718.         memmove(v->is_intra_base, v->is_intra, sizeof(v->is_intra_base[0])*s->mb_stride);

  1719.     return 0;

  1720. }

  1721. 

  1722. static int vc1_decode_p_mb_intfi(VC1Context *v)

  1723. {

  1724.     MpegEncContext *s = &v->s;

  1725.     GetBitContext *gb = &s->gb;

  1726.     int i;

  1727.     int mb_pos = s->mb_x + s->mb_y * s->mb_stride;

  1728.     int cbp = 0; /* cbp decoding stuff */

  1729.     int mqdiff, mquant; /* MB quantization */

  1730.     int ttmb = v->ttfrm; /* MB Transform type */

  1731. 

  1732.     int mb_has_coeffs = 1; /* last_flag */

  1733.     int dmv_x, dmv_y; /* Differential MV components */

  1734.     int val; /* temp values */

  1735.     int first_block = 1;

  1736.     int dst_idx, off;

  1737.     int pred_flag = 0;

  1738.     int block_cbp = 0, pat, block_tt = 0;

  1739.     int idx_mbmode = 0;

  1740. 

  1741.     mquant = v->pq; /* Lossy initialization */

  1742. 

  1743.     idx_mbmode = get_vlc2(gb, v->mbmode_vlc->table, VC1_IF_MBMODE_VLC_BITS, 2);

  1744.     if (idx_mbmode <= 1) { // intra MB

  1745.         v->is_intra[s->mb_x] = 0x3f; // Set the bitfield to all 1.

  1746.         s->mb_intra          = 1;

  1747.         s->current_picture.motion_val[1][s->block_index[0] + v->blocks_off][0] = 0;

  1748.         s->current_picture.motion_val[1][s->block_index[0] + v->blocks_off][1] = 0;

  1749.         s->current_picture.mb_type[mb_pos + v->mb_off] = MB_TYPE_INTRA;

  1750.         GET_MQUANT();

  1751.         s->current_picture.qscale_table[mb_pos] = mquant;

  1752.         /* Set DC scale - y and c use the same (not sure if necessary here) */

  1753.         s->y_dc_scale = s->y_dc_scale_table[mquant];

  1754.         s->c_dc_scale = s->c_dc_scale_table[mquant];

  1755.         v->s.ac_pred  = v->acpred_plane[mb_pos] = get_bits1(gb);

  1756.         mb_has_coeffs = idx_mbmode & 1;

  1757.         if (mb_has_coeffs)

  1758.             cbp = 1 + get_vlc2(&v->s.gb, v->cbpcy_vlc->table, VC1_ICBPCY_VLC_BITS, 2);

  1759.         dst_idx = 0;

  1760.         for (i = 0; i < 6; i++) {

  1761.             v->a_avail = v->c_avail          = 0;

  1762.             v->mb_type[0][s->block_index[i]] = 1;

  1763.             s->dc_val[0][s->block_index[i]]  = 0;

  1764.             dst_idx += i >> 2;

  1765.             val = ((cbp >> (5 - i)) & 1);

  1766.             if (i == 2 || i == 3 || !s->first_slice_line)

  1767.                 v->a_avail = v->mb_type[0][s->block_index[i] - s->block_wrap[i]];

  1768.             if (i == 1 || i == 3 || s->mb_x)

  1769.                 v->c_avail = v->mb_type[0][s->block_index[i] - 1];

  1770. 

  1771.             vc1_decode_intra_block(v, s->block[i], i, val, mquant,

  1772.                                    (i & 4) ? v->codingset2 : v->codingset);

  1773.             if (CONFIG_GRAY && (i > 3) && (s->avctx->flags & AV_CODEC_FLAG_GRAY))

  1774.                 continue;

  1775.             v->vc1dsp.vc1_inv_trans_8x8(s->block[i]);

  1776.             off  = (i & 4) ? 0 : ((i & 1) * 8 + (i & 2) * 4 * s->linesize);

  1777.             s->idsp.put_signed_pixels_clamped(s->block[i],

  1778.                                               s->dest[dst_idx] + off,

  1779.                                               (i & 4) ? s->uvlinesize

  1780.                                                       : s->linesize);

  1781.             // TODO: loop filter

  1782.         }

  1783.     } else {

  1784.         s->mb_intra = v->is_intra[s->mb_x] = 0;

  1785.         s->current_picture.mb_type[mb_pos + v->mb_off] = MB_TYPE_16x16;

  1786.         for (i = 0; i < 6; i++)

  1787.             v->mb_type[0][s->block_index[i]] = 0;

  1788.         if (idx_mbmode <= 5) { // 1-MV

  1789.             dmv_x = dmv_y = pred_flag = 0;

  1790.             if (idx_mbmode & 1) {

  1791.                 get_mvdata_interlaced(v, &dmv_x, &dmv_y, &pred_flag);

  1792.             }

  1793.             ff_vc1_pred_mv(v, 0, dmv_x, dmv_y, 1, v->range_x, v->range_y, v->mb_type[0], pred_flag, 0);

  1794.             ff_vc1_mc_1mv(v, 0);

  1795.             mb_has_coeffs = !(idx_mbmode & 2);

  1796.         } else { // 4-MV

  1797.             v->fourmvbp = get_vlc2(gb, v->fourmvbp_vlc->table, VC1_4MV_BLOCK_PATTERN_VLC_BITS, 1);

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

  1799.                 dmv_x = dmv_y = pred_flag = 0;

  1800.                 if (v->fourmvbp & (8 >> i))

  1801.                     get_mvdata_interlaced(v, &dmv_x, &dmv_y, &pred_flag);

  1802.                 ff_vc1_pred_mv(v, i, dmv_x, dmv_y, 0, v->range_x, v->range_y, v->mb_type[0], pred_flag, 0);

  1803.                 ff_vc1_mc_4mv_luma(v, i, 0, 0);

  1804.             }

  1805.             ff_vc1_mc_4mv_chroma(v, 0);

  1806.             mb_has_coeffs = idx_mbmode & 1;

  1807.         }

  1808.         if (mb_has_coeffs)

  1809.             cbp = 1 + get_vlc2(&v->s.gb, v->cbpcy_vlc->table, VC1_CBPCY_P_VLC_BITS, 2);

  1810.         if (cbp) {

  1811.             GET_MQUANT();

  1812.         }

  1813.         s->current_picture.qscale_table[mb_pos] = mquant;

  1814.         if (!v->ttmbf && cbp) {

  1815.             ttmb = get_vlc2(gb, ff_vc1_ttmb_vlc[v->tt_index].table, VC1_TTMB_VLC_BITS, 2);

  1816.         }

  1817.         dst_idx = 0;

  1818.         for (i = 0; i < 6; i++) {

  1819.             s->dc_val[0][s->block_index[i]] = 0;

  1820.             dst_idx += i >> 2;

  1821.             val = ((cbp >> (5 - i)) & 1);

  1822.             off = (i & 4) ? 0 : (i & 1) * 8 + (i & 2) * 4 * s->linesize;

  1823.             if (val) {

  1824.                 pat = vc1_decode_p_block(v, s->block[i], i, mquant, ttmb,

  1825.                                          first_block, s->dest[dst_idx] + off,

  1826.                                          (i & 4) ? s->uvlinesize : s->linesize,

  1827.                                          CONFIG_GRAY && (i & 4) && (s->avctx->flags & AV_CODEC_FLAG_GRAY),

  1828.                                          &block_tt);

  1829.                 block_cbp |= pat << (i << 2);

  1830.                 if (!v->ttmbf && ttmb < 8)

  1831.                     ttmb = -1;

  1832.                 first_block = 0;

  1833.             }

  1834.         }

  1835.     }

  1836.     if (s->mb_x == s->mb_width - 1)

  1837.         memmove(v->is_intra_base, v->is_intra, sizeof(v->is_intra_base[0]) * s->mb_stride);

  1838.     return 0;

  1839. }

  1840. 

  1841. /** Decode one B-frame MB (in Main profile)

  1842.  */

  1843. static void vc1_decode_b_mb(VC1Context *v)

  1844. {

  1845.     MpegEncContext *s = &v->s;

  1846.     GetBitContext *gb = &s->gb;

  1847.     int i, j;

  1848.     int mb_pos = s->mb_x + s->mb_y * s->mb_stride;

  1849.     int cbp = 0; /* cbp decoding stuff */

  1850.     int mqdiff, mquant; /* MB quantization */

  1851.     int ttmb = v->ttfrm; /* MB Transform type */

  1852.     int mb_has_coeffs = 0; /* last_flag */

  1853.     int index, index1; /* LUT indexes */

  1854.     int val, sign; /* temp values */

  1855.     int first_block = 1;

  1856.     int dst_idx, off;

  1857.     int skipped, direct;

  1858.     int dmv_x[2], dmv_y[2];

  1859.     int bmvtype = BMV_TYPE_BACKWARD;

  1860. 

  1861.     mquant      = v->pq; /* lossy initialization */

  1862.     s->mb_intra = 0;

  1863. 

  1864.     if (v->dmb_is_raw)

  1865.         direct = get_bits1(gb);

  1866.     else

  1867.         direct = v->direct_mb_plane[mb_pos];

  1868.     if (v->skip_is_raw)

  1869.         skipped = get_bits1(gb);

  1870.     else

  1871.         skipped = v->s.mbskip_table[mb_pos];

  1872. 

  1873.     dmv_x[0] = dmv_x[1] = dmv_y[0] = dmv_y[1] = 0;

  1874.     for (i = 0; i < 6; i++) {

  1875.         v->mb_type[0][s->block_index[i]] = 0;

  1876.         s->dc_val[0][s->block_index[i]]  = 0;

  1877.     }

  1878.     s->current_picture.qscale_table[mb_pos] = 0;

  1879. 

  1880.     if (!direct) {

  1881.         if (!skipped) {

  1882.             GET_MVDATA(dmv_x[0], dmv_y[0]);

  1883.             dmv_x[1] = dmv_x[0];

  1884.             dmv_y[1] = dmv_y[0];

  1885.         }

  1886.         if (skipped || !s->mb_intra) {

  1887.             bmvtype = decode012(gb);

  1888.             switch (bmvtype) {

  1889.             case 0:

  1890.                 bmvtype = (v->bfraction >= (B_FRACTION_DEN/2)) ? BMV_TYPE_BACKWARD : BMV_TYPE_FORWARD;

  1891.                 break;

  1892.             case 1:

  1893.                 bmvtype = (v->bfraction >= (B_FRACTION_DEN/2)) ? BMV_TYPE_FORWARD : BMV_TYPE_BACKWARD;

  1894.                 break;

  1895.             case 2:

  1896.                 bmvtype  = BMV_TYPE_INTERPOLATED;

  1897.                 dmv_x[0] = dmv_y[0] = 0;

  1898.             }

  1899.         }

  1900.     }

  1901.     for (i = 0; i < 6; i++)

  1902.         v->mb_type[0][s->block_index[i]] = s->mb_intra;

  1903. 

  1904.     if (skipped) {

  1905.         if (direct)

  1906.             bmvtype = BMV_TYPE_INTERPOLATED;

  1907.         ff_vc1_pred_b_mv(v, dmv_x, dmv_y, direct, bmvtype);

  1908.         vc1_b_mc(v, dmv_x, dmv_y, direct, bmvtype);

  1909.         return;

  1910.     }

  1911.     if (direct) {

  1912.         cbp = get_vlc2(&v->s.gb, v->cbpcy_vlc->table, VC1_CBPCY_P_VLC_BITS, 2);

  1913.         GET_MQUANT();

  1914.         s->mb_intra = 0;

  1915.         s->current_picture.qscale_table[mb_pos] = mquant;

  1916.         if (!v->ttmbf)

  1917.             ttmb = get_vlc2(gb, ff_vc1_ttmb_vlc[v->tt_index].table, VC1_TTMB_VLC_BITS, 2);

  1918.         dmv_x[0] = dmv_y[0] = dmv_x[1] = dmv_y[1] = 0;

  1919.         ff_vc1_pred_b_mv(v, dmv_x, dmv_y, direct, bmvtype);

  1920.         vc1_b_mc(v, dmv_x, dmv_y, direct, bmvtype);

  1921.     } else {

  1922.         if (!mb_has_coeffs && !s->mb_intra) {

  1923.             /* no coded blocks - effectively skipped */

  1924.             ff_vc1_pred_b_mv(v, dmv_x, dmv_y, direct, bmvtype);

  1925.             vc1_b_mc(v, dmv_x, dmv_y, direct, bmvtype);

  1926.             return;

  1927.         }

  1928.         if (s->mb_intra && !mb_has_coeffs) {

  1929.             GET_MQUANT();

  1930.             s->current_picture.qscale_table[mb_pos] = mquant;

  1931.             s->ac_pred = get_bits1(gb);

  1932.             cbp = 0;

  1933.             ff_vc1_pred_b_mv(v, dmv_x, dmv_y, direct, bmvtype);

  1934.         } else {

  1935.             if (bmvtype == BMV_TYPE_INTERPOLATED) {

  1936.                 GET_MVDATA(dmv_x[0], dmv_y[0]);

  1937.                 if (!mb_has_coeffs) {

  1938.                     /* interpolated skipped block */

  1939.                     ff_vc1_pred_b_mv(v, dmv_x, dmv_y, direct, bmvtype);

  1940.                     vc1_b_mc(v, dmv_x, dmv_y, direct, bmvtype);

  1941.                     return;

  1942.                 }

  1943.             }

  1944.             ff_vc1_pred_b_mv(v, dmv_x, dmv_y, direct, bmvtype);

  1945.             if (!s->mb_intra) {

  1946.                 vc1_b_mc(v, dmv_x, dmv_y, direct, bmvtype);

  1947.             }

  1948.             if (s->mb_intra)

  1949.                 s->ac_pred = get_bits1(gb);

  1950.             cbp = get_vlc2(&v->s.gb, v->cbpcy_vlc->table, VC1_CBPCY_P_VLC_BITS, 2);

  1951.             GET_MQUANT();

  1952.             s->current_picture.qscale_table[mb_pos] = mquant;

  1953.             if (!v->ttmbf && !s->mb_intra && mb_has_coeffs)

  1954.                 ttmb = get_vlc2(gb, ff_vc1_ttmb_vlc[v->tt_index].table, VC1_TTMB_VLC_BITS, 2);

  1955.         }

  1956.     }

  1957.     dst_idx = 0;

  1958.     for (i = 0; i < 6; i++) {

  1959.         s->dc_val[0][s->block_index[i]] = 0;

  1960.         dst_idx += i >> 2;

  1961.         val = ((cbp >> (5 - i)) & 1);

  1962.         off = (i & 4) ? 0 : ((i & 1) * 8 + (i & 2) * 4 * s->linesize);

  1963.         v->mb_type[0][s->block_index[i]] = s->mb_intra;

  1964.         if (s->mb_intra) {

  1965.             /* check if prediction blocks A and C are available */

  1966.             v->a_avail = v->c_avail = 0;

  1967.             if (i == 2 || i == 3 || !s->first_slice_line)

  1968.                 v->a_avail = v->mb_type[0][s->block_index[i] - s->block_wrap[i]];

  1969.             if (i == 1 || i == 3 || s->mb_x)

  1970.                 v->c_avail = v->mb_type[0][s->block_index[i] - 1];

  1971. 

  1972.             vc1_decode_intra_block(v, s->block[i], i, val, mquant,

  1973.                                    (i & 4) ? v->codingset2 : v->codingset);

  1974.             if (CONFIG_GRAY && (i > 3) && (s->avctx->flags & AV_CODEC_FLAG_GRAY))

  1975.                 continue;

  1976.             v->vc1dsp.vc1_inv_trans_8x8(s->block[i]);

  1977.             if (v->rangeredfrm)

  1978.                 for (j = 0; j < 64; j++)

  1979.                     s->block[i][j] <<= 1;

  1980.             s->idsp.put_signed_pixels_clamped(s->block[i],

  1981.                                               s->dest[dst_idx] + off,

  1982.                                               i & 4 ? s->uvlinesize

  1983.                                                     : s->linesize);

  1984.         } else if (val) {

  1985.             vc1_decode_p_block(v, s->block[i], i, mquant, ttmb,

  1986.                                first_block, s->dest[dst_idx] + off,

  1987.                                (i & 4) ? s->uvlinesize : s->linesize,

  1988.                                CONFIG_GRAY && (i & 4) && (s->avctx->flags & AV_CODEC_FLAG_GRAY), NULL);

  1989.             if (!v->ttmbf && ttmb < 8)

  1990.                 ttmb = -1;

  1991.             first_block = 0;

  1992.         }

  1993.     }

  1994. }

  1995. 

  1996. /** Decode one B-frame MB (in interlaced field B picture)

  1997.  */

  1998. static void vc1_decode_b_mb_intfi(VC1Context *v)

  1999. {

  2000.     MpegEncContext *s = &v->s;

  2001.     GetBitContext *gb = &s->gb;

  2002.     int i, j;

  2003.     int mb_pos = s->mb_x + s->mb_y * s->mb_stride;

  2004.     int cbp = 0; /* cbp decoding stuff */

  2005.     int mqdiff, mquant; /* MB quantization */

  2006.     int ttmb = v->ttfrm; /* MB Transform type */

  2007.     int mb_has_coeffs = 0; /* last_flag */

  2008.     int val; /* temp value */

  2009.     int first_block = 1;

  2010.     int dst_idx, off;

  2011.     int fwd;

  2012.     int dmv_x[2], dmv_y[2], pred_flag[2];

  2013.     int bmvtype = BMV_TYPE_BACKWARD;

  2014.     int idx_mbmode;

  2015. 

  2016.     mquant      = v->pq; /* Lossy initialization */

  2017.     s->mb_intra = 0;

  2018. 

  2019.     idx_mbmode = get_vlc2(gb, v->mbmode_vlc->table, VC1_IF_MBMODE_VLC_BITS, 2);

  2020.     if (idx_mbmode <= 1) { // intra MB

  2021.         v->is_intra[s->mb_x] = 0x3f; // Set the bitfield to all 1.

  2022.         s->mb_intra          = 1;

  2023.         s->current_picture.motion_val[1][s->block_index[0]][0] = 0;

  2024.         s->current_picture.motion_val[1][s->block_index[0]][1] = 0;

  2025.         s->current_picture.mb_type[mb_pos + v->mb_off]         = MB_TYPE_INTRA;

  2026.         GET_MQUANT();

  2027.         s->current_picture.qscale_table[mb_pos] = mquant;

  2028.         /* Set DC scale - y and c use the same (not sure if necessary here) */

  2029.         s->y_dc_scale = s->y_dc_scale_table[mquant];

  2030.         s->c_dc_scale = s->c_dc_scale_table[mquant];

  2031.         v->s.ac_pred  = v->acpred_plane[mb_pos] = get_bits1(gb);

  2032.         mb_has_coeffs = idx_mbmode & 1;

  2033.         if (mb_has_coeffs)

  2034.             cbp = 1 + get_vlc2(&v->s.gb, v->cbpcy_vlc->table, VC1_ICBPCY_VLC_BITS, 2);

  2035.         dst_idx = 0;

  2036.         for (i = 0; i < 6; i++) {

  2037.             v->a_avail = v->c_avail          = 0;

  2038.             v->mb_type[0][s->block_index[i]] = 1;

  2039.             s->dc_val[0][s->block_index[i]]  = 0;

  2040.             dst_idx += i >> 2;

  2041.             val = ((cbp >> (5 - i)) & 1);

  2042.             if (i == 2 || i == 3 || !s->first_slice_line)

  2043.                 v->a_avail = v->mb_type[0][s->block_index[i] - s->block_wrap[i]];

  2044.             if (i == 1 || i == 3 || s->mb_x)

  2045.                 v->c_avail = v->mb_type[0][s->block_index[i] - 1];

  2046. 

  2047.             vc1_decode_intra_block(v, s->block[i], i, val, mquant,

  2048.                                    (i & 4) ? v->codingset2 : v->codingset);

  2049.             if (CONFIG_GRAY && (i > 3) && (s->avctx->flags & AV_CODEC_FLAG_GRAY))

  2050.                 continue;

  2051.             v->vc1dsp.vc1_inv_trans_8x8(s->block[i]);

  2052.             if (v->rangeredfrm)

  2053.                 for (j = 0; j < 64; j++)

  2054.                     s->block[i][j] <<= 1;

  2055.             off  = (i & 4) ? 0 : ((i & 1) * 8 + (i & 2) * 4 * s->linesize);

  2056.             s->idsp.put_signed_pixels_clamped(s->block[i],

  2057.                                               s->dest[dst_idx] + off,

  2058.                                               (i & 4) ? s->uvlinesize

  2059.                                                       : s->linesize);

  2060.             // TODO: yet to perform loop filter

  2061.         }

  2062.     } else {

  2063.         s->mb_intra = v->is_intra[s->mb_x] = 0;

  2064.         s->current_picture.mb_type[mb_pos + v->mb_off] = MB_TYPE_16x16;

  2065.         for (i = 0; i < 6; i++)

  2066.             v->mb_type[0][s->block_index[i]] = 0;

  2067.         if (v->fmb_is_raw)

  2068.             fwd = v->forward_mb_plane[mb_pos] = get_bits1(gb);

  2069.         else

  2070.             fwd = v->forward_mb_plane[mb_pos];

  2071.         if (idx_mbmode <= 5) { // 1-MV

  2072.             int interpmvp = 0;

  2073.             dmv_x[0]     = dmv_x[1] = dmv_y[0] = dmv_y[1] = 0;

  2074.             pred_flag[0] = pred_flag[1] = 0;

  2075.             if (fwd)

  2076.                 bmvtype = BMV_TYPE_FORWARD;

  2077.             else {

  2078.                 bmvtype = decode012(gb);

  2079.                 switch (bmvtype) {

  2080.                 case 0:

  2081.                     bmvtype = BMV_TYPE_BACKWARD;

  2082.                     break;

  2083.                 case 1:

  2084.                     bmvtype = BMV_TYPE_DIRECT;

  2085.                     break;

  2086.                 case 2:

  2087.                     bmvtype   = BMV_TYPE_INTERPOLATED;

  2088.                     interpmvp = get_bits1(gb);

  2089.                 }

  2090.             }

  2091.             v->bmvtype = bmvtype;

  2092.             if (bmvtype != BMV_TYPE_DIRECT && idx_mbmode & 1) {

  2093.                 get_mvdata_interlaced(v, &dmv_x[bmvtype == BMV_TYPE_BACKWARD], &dmv_y[bmvtype == BMV_TYPE_BACKWARD], &pred_flag[bmvtype == BMV_TYPE_BACKWARD]);

  2094.             }

  2095.             if (interpmvp) {

  2096.                 get_mvdata_interlaced(v, &dmv_x[1], &dmv_y[1], &pred_flag[1]);

  2097.             }

  2098.             if (bmvtype == BMV_TYPE_DIRECT) {

  2099.                 dmv_x[0] = dmv_y[0] = pred_flag[0] = 0;

  2100.                 dmv_x[1] = dmv_y[1] = pred_flag[0] = 0;

  2101.                 if (!s->next_picture_ptr->field_picture) {

  2102.                     av_log(s->avctx, AV_LOG_ERROR, "Mixed field/frame direct mode not supported\n");

  2103.                     return;

  2104.                 }

  2105.             }

  2106.             ff_vc1_pred_b_mv_intfi(v, 0, dmv_x, dmv_y, 1, pred_flag);

  2107.             vc1_b_mc(v, dmv_x, dmv_y, (bmvtype == BMV_TYPE_DIRECT), bmvtype);

  2108.             mb_has_coeffs = !(idx_mbmode & 2);

  2109.         } else { // 4-MV

  2110.             if (fwd)

  2111.                 bmvtype = BMV_TYPE_FORWARD;

  2112.             v->bmvtype  = bmvtype;

  2113.             v->fourmvbp = get_vlc2(gb, v->fourmvbp_vlc->table, VC1_4MV_BLOCK_PATTERN_VLC_BITS, 1);

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

  2115.                 dmv_x[0] = dmv_y[0] = pred_flag[0] = 0;

  2116.                 dmv_x[1] = dmv_y[1] = pred_flag[1] = 0;

  2117.                 if (v->fourmvbp & (8 >> i)) {

  2118.                     get_mvdata_interlaced(v, &dmv_x[bmvtype == BMV_TYPE_BACKWARD],

  2119.                                              &dmv_y[bmvtype == BMV_TYPE_BACKWARD],

  2120.                                          &pred_flag[bmvtype == BMV_TYPE_BACKWARD]);

  2121.                 }

  2122.                 ff_vc1_pred_b_mv_intfi(v, i, dmv_x, dmv_y, 0, pred_flag);

  2123.                 ff_vc1_mc_4mv_luma(v, i, bmvtype == BMV_TYPE_BACKWARD, 0);

  2124.             }

  2125.             ff_vc1_mc_4mv_chroma(v, bmvtype == BMV_TYPE_BACKWARD);

  2126.             mb_has_coeffs = idx_mbmode & 1;

  2127.         }

  2128.         if (mb_has_coeffs)

  2129.             cbp = 1 + get_vlc2(&v->s.gb, v->cbpcy_vlc->table, VC1_CBPCY_P_VLC_BITS, 2);

  2130.         if (cbp) {

  2131.             GET_MQUANT();

  2132.         }

  2133.         s->current_picture.qscale_table[mb_pos] = mquant;

  2134.         if (!v->ttmbf && cbp) {

  2135.             ttmb = get_vlc2(gb, ff_vc1_ttmb_vlc[v->tt_index].table, VC1_TTMB_VLC_BITS, 2);

  2136.         }

  2137.         dst_idx = 0;

  2138.         for (i = 0; i < 6; i++) {

  2139.             s->dc_val[0][s->block_index[i]] = 0;

  2140.             dst_idx += i >> 2;

  2141.             val = ((cbp >> (5 - i)) & 1);

  2142.             off = (i & 4) ? 0 : (i & 1) * 8 + (i & 2) * 4 * s->linesize;

  2143.             if (val) {

  2144.                 vc1_decode_p_block(v, s->block[i], i, mquant, ttmb,

  2145.                                    first_block, s->dest[dst_idx] + off,

  2146.                                    (i & 4) ? s->uvlinesize : s->linesize,

  2147.                                    CONFIG_GRAY && (i & 4) && (s->avctx->flags & AV_CODEC_FLAG_GRAY), NULL);

  2148.                 if (!v->ttmbf && ttmb < 8)

  2149.                     ttmb = -1;

  2150.                 first_block = 0;

  2151.             }

  2152.         }

  2153.     }

  2154. }

  2155. 

  2156. /** Decode one B-frame MB (in interlaced frame B picture)

  2157.  */

  2158. static int vc1_decode_b_mb_intfr(VC1Context *v)

  2159. {

  2160.     MpegEncContext *s = &v->s;

  2161.     GetBitContext *gb = &s->gb;

  2162.     int i, j;

  2163.     int mb_pos = s->mb_x + s->mb_y * s->mb_stride;

  2164.     int cbp = 0; /* cbp decoding stuff */

  2165.     int mqdiff, mquant; /* MB quantization */

  2166.     int ttmb = v->ttfrm; /* MB Transform type */

  2167.     int mvsw = 0; /* motion vector switch */

  2168.     int mb_has_coeffs = 1; /* last_flag */

  2169.     int dmv_x, dmv_y; /* Differential MV components */

  2170.     int val; /* temp value */

  2171.     int first_block = 1;

  2172.     int dst_idx, off;

  2173.     int skipped, direct, twomv = 0;

  2174.     int block_cbp = 0, pat, block_tt = 0;

  2175.     int idx_mbmode = 0, mvbp;

  2176.     int stride_y, fieldtx;

  2177.     int bmvtype = BMV_TYPE_BACKWARD;

  2178.     int dir, dir2;

  2179. 

  2180.     mquant = v->pq; /* Lossy initialization */

  2181.     s->mb_intra = 0;

  2182.     if (v->skip_is_raw)

  2183.         skipped = get_bits1(gb);

  2184.     else

  2185.         skipped = v->s.mbskip_table[mb_pos];

  2186. 

  2187.     if (!skipped) {

  2188.         idx_mbmode = get_vlc2(gb, v->mbmode_vlc->table, VC1_INTFR_NON4MV_MBMODE_VLC_BITS, 2);

  2189.         if (ff_vc1_mbmode_intfrp[0][idx_mbmode][0] == MV_PMODE_INTFR_2MV_FIELD) {

  2190.             twomv = 1;

  2191.             v->blk_mv_type[s->block_index[0]] = 1;

  2192.             v->blk_mv_type[s->block_index[1]] = 1;

  2193.             v->blk_mv_type[s->block_index[2]] = 1;

  2194.             v->blk_mv_type[s->block_index[3]] = 1;

  2195.         } else {

  2196.             v->blk_mv_type[s->block_index[0]] = 0;

  2197.             v->blk_mv_type[s->block_index[1]] = 0;

  2198.             v->blk_mv_type[s->block_index[2]] = 0;

  2199.             v->blk_mv_type[s->block_index[3]] = 0;

  2200.         }

  2201.     }

  2202. 

  2203.     if (v->dmb_is_raw)

  2204.         direct = get_bits1(gb);

  2205.     else

  2206.         direct = v->direct_mb_plane[mb_pos];

  2207. 

  2208.     if (direct) {

  2209.         if (s->next_picture_ptr->field_picture)

  2210.             av_log(s->avctx, AV_LOG_WARNING, "Mixed frame/field direct mode not supported\n");

  2211.         s->mv[0][0][0] = s->current_picture.motion_val[0][s->block_index[0]][0] = scale_mv(s->next_picture.motion_val[1][s->block_index[0]][0], v->bfraction, 0, s->quarter_sample);

  2212.         s->mv[0][0][1] = s->current_picture.motion_val[0][s->block_index[0]][1] = scale_mv(s->next_picture.motion_val[1][s->block_index[0]][1], v->bfraction, 0, s->quarter_sample);

  2213.         s->mv[1][0][0] = s->current_picture.motion_val[1][s->block_index[0]][0] = scale_mv(s->next_picture.motion_val[1][s->block_index[0]][0], v->bfraction, 1, s->quarter_sample);

  2214.         s->mv[1][0][1] = s->current_picture.motion_val[1][s->block_index[0]][1] = scale_mv(s->next_picture.motion_val[1][s->block_index[0]][1], v->bfraction, 1, s->quarter_sample);

  2215. 

  2216.         if (twomv) {

  2217.             s->mv[0][2][0] = s->current_picture.motion_val[0][s->block_index[2]][0] = scale_mv(s->next_picture.motion_val[1][s->block_index[2]][0], v->bfraction, 0, s->quarter_sample);

  2218.             s->mv[0][2][1] = s->current_picture.motion_val[0][s->block_index[2]][1] = scale_mv(s->next_picture.motion_val[1][s->block_index[2]][1], v->bfraction, 0, s->quarter_sample);

  2219.             s->mv[1][2][0] = s->current_picture.motion_val[1][s->block_index[2]][0] = scale_mv(s->next_picture.motion_val[1][s->block_index[2]][0], v->bfraction, 1, s->quarter_sample);

  2220.             s->mv[1][2][1] = s->current_picture.motion_val[1][s->block_index[2]][1] = scale_mv(s->next_picture.motion_val[1][s->block_index[2]][1], v->bfraction, 1, s->quarter_sample);

  2221. 

  2222.             for (i = 1; i < 4; i += 2) {

  2223.                 s->mv[0][i][0] = s->current_picture.motion_val[0][s->block_index[i]][0] = s->mv[0][i-1][0];

  2224.                 s->mv[0][i][1] = s->current_picture.motion_val[0][s->block_index[i]][1] = s->mv[0][i-1][1];

  2225.                 s->mv[1][i][0] = s->current_picture.motion_val[1][s->block_index[i]][0] = s->mv[1][i-1][0];

  2226.                 s->mv[1][i][1] = s->current_picture.motion_val[1][s->block_index[i]][1] = s->mv[1][i-1][1];

  2227.             }

  2228.         } else {

  2229.             for (i = 1; i < 4; i++) {

  2230.                 s->mv[0][i][0] = s->current_picture.motion_val[0][s->block_index[i]][0] = s->mv[0][0][0];

  2231.                 s->mv[0][i][1] = s->current_picture.motion_val[0][s->block_index[i]][1] = s->mv[0][0][1];

  2232.                 s->mv[1][i][0] = s->current_picture.motion_val[1][s->block_index[i]][0] = s->mv[1][0][0];

  2233.                 s->mv[1][i][1] = s->current_picture.motion_val[1][s->block_index[i]][1] = s->mv[1][0][1];

  2234.             }

  2235.         }

  2236.     }

  2237. 

  2238.     if (ff_vc1_mbmode_intfrp[0][idx_mbmode][0] == MV_PMODE_INTFR_INTRA) { // intra MB

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

  2240.             s->mv[0][i][0] = s->current_picture.motion_val[0][s->block_index[i]][0] = 0;

  2241.             s->mv[0][i][1] = s->current_picture.motion_val[0][s->block_index[i]][1] = 0;

  2242.             s->mv[1][i][0] = s->current_picture.motion_val[1][s->block_index[i]][0] = 0;

  2243.             s->mv[1][i][1] = s->current_picture.motion_val[1][s->block_index[i]][1] = 0;

  2244.         }

  2245.         v->is_intra[s->mb_x] = 0x3f; // Set the bitfield to all 1.

  2246.         s->mb_intra          = 1;

  2247.         s->current_picture.mb_type[mb_pos] = MB_TYPE_INTRA;

  2248.         fieldtx = v->fieldtx_plane[mb_pos] = get_bits1(gb);

  2249.         mb_has_coeffs = get_bits1(gb);

  2250.         if (mb_has_coeffs)

  2251.             cbp = 1 + get_vlc2(&v->s.gb, v->cbpcy_vlc->table, VC1_CBPCY_P_VLC_BITS, 2);

  2252.         v->s.ac_pred = v->acpred_plane[mb_pos] = get_bits1(gb);

  2253.         GET_MQUANT();

  2254.         s->current_picture.qscale_table[mb_pos] = mquant;

  2255.         /* Set DC scale - y and c use the same (not sure if necessary here) */

  2256.         s->y_dc_scale = s->y_dc_scale_table[mquant];

  2257.         s->c_dc_scale = s->c_dc_scale_table[mquant];

  2258.         dst_idx = 0;

  2259.         for (i = 0; i < 6; i++) {

  2260.             v->a_avail = v->c_avail          = 0;

  2261.             v->mb_type[0][s->block_index[i]] = 1;

  2262.             s->dc_val[0][s->block_index[i]]  = 0;

  2263.             dst_idx += i >> 2;

  2264.             val = ((cbp >> (5 - i)) & 1);

  2265.             if (i == 2 || i == 3 || !s->first_slice_line)

  2266.                 v->a_avail = v->mb_type[0][s->block_index[i] - s->block_wrap[i]];

  2267.             if (i == 1 || i == 3 || s->mb_x)

  2268.                 v->c_avail = v->mb_type[0][s->block_index[i] - 1];

  2269. 

  2270.             vc1_decode_intra_block(v, s->block[i], i, val, mquant,

  2271.                                    (i & 4) ? v->codingset2 : v->codingset);

  2272.             if (CONFIG_GRAY && i > 3 && (s->avctx->flags & AV_CODEC_FLAG_GRAY))

  2273.                 continue;

  2274.             v->vc1dsp.vc1_inv_trans_8x8(s->block[i]);

  2275.             if (i < 4) {

  2276.                 stride_y = s->linesize << fieldtx;

  2277.                 off = (fieldtx) ? ((i & 1) * 8) + ((i & 2) >> 1) * s->linesize : (i & 1) * 8 + 4 * (i & 2) * s->linesize;

  2278.             } else {

  2279.                 stride_y = s->uvlinesize;

  2280.                 off = 0;

  2281.             }

  2282.             s->idsp.put_signed_pixels_clamped(s->block[i],

  2283.                                               s->dest[dst_idx] + off,

  2284.                                               stride_y);

  2285.         }

  2286.     } else {

  2287.         s->mb_intra = v->is_intra[s->mb_x] = 0;

  2288.         if (!direct) {

  2289.             if (skipped || !s->mb_intra) {

  2290.                 bmvtype = decode012(gb);

  2291.                 switch (bmvtype) {

  2292.                 case 0:

  2293.                     bmvtype = (v->bfraction >= (B_FRACTION_DEN/2)) ? BMV_TYPE_BACKWARD : BMV_TYPE_FORWARD;

  2294.                     break;

  2295.                 case 1:

  2296.                     bmvtype = (v->bfraction >= (B_FRACTION_DEN/2)) ? BMV_TYPE_FORWARD : BMV_TYPE_BACKWARD;

  2297.                     break;

  2298.                 case 2:

  2299.                     bmvtype  = BMV_TYPE_INTERPOLATED;

  2300.                 }

  2301.             }

  2302. 

  2303.             if (twomv && bmvtype != BMV_TYPE_INTERPOLATED)

  2304.                 mvsw = get_bits1(gb);

  2305.         }

  2306. 

  2307.         if (!skipped) { // inter MB

  2308.             mb_has_coeffs = ff_vc1_mbmode_intfrp[0][idx_mbmode][3];

  2309.             if (mb_has_coeffs)

  2310.                 cbp = 1 + get_vlc2(&v->s.gb, v->cbpcy_vlc->table, VC1_CBPCY_P_VLC_BITS, 2);

  2311.             if (!direct) {

  2312.                 if (bmvtype == BMV_TYPE_INTERPOLATED && twomv) {

  2313.                     v->fourmvbp = get_vlc2(gb, v->fourmvbp_vlc->table, VC1_4MV_BLOCK_PATTERN_VLC_BITS, 1);

  2314.                 } else if (bmvtype == BMV_TYPE_INTERPOLATED || twomv) {

  2315.                     v->twomvbp = get_vlc2(gb, v->twomvbp_vlc->table, VC1_2MV_BLOCK_PATTERN_VLC_BITS, 1);

  2316.                 }

  2317.             }

  2318. 

  2319.             for (i = 0; i < 6; i++)

  2320.                 v->mb_type[0][s->block_index[i]] = 0;

  2321.             fieldtx = v->fieldtx_plane[mb_pos] = ff_vc1_mbmode_intfrp[0][idx_mbmode][1];

  2322.             /* for all motion vector read MVDATA and motion compensate each block */

  2323.             dst_idx = 0;

  2324.             if (direct) {

  2325.                 if (twomv) {

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

  2327.                         ff_vc1_mc_4mv_luma(v, i, 0, 0);

  2328.                         ff_vc1_mc_4mv_luma(v, i, 1, 1);

  2329.                     }

  2330.                     ff_vc1_mc_4mv_chroma4(v, 0, 0, 0);

  2331.                     ff_vc1_mc_4mv_chroma4(v, 1, 1, 1);

  2332.                 } else {

  2333.                     ff_vc1_mc_1mv(v, 0);

  2334.                     ff_vc1_interp_mc(v);

  2335.                 }

  2336.             } else if (twomv && bmvtype == BMV_TYPE_INTERPOLATED) {

  2337.                 mvbp = v->fourmvbp;

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

  2339.                     dir = i==1 || i==3;

  2340.                     dmv_x = dmv_y = 0;

  2341.                     val = ((mvbp >> (3 - i)) & 1);

  2342.                     if (val)

  2343.                         get_mvdata_interlaced(v, &dmv_x, &dmv_y, 0);

  2344.                     j = i > 1 ? 2 : 0;

  2345.                     ff_vc1_pred_mv_intfr(v, j, dmv_x, dmv_y, 2, v->range_x, v->range_y, v->mb_type[0], dir);

  2346.                     ff_vc1_mc_4mv_luma(v, j, dir, dir);

  2347.                     ff_vc1_mc_4mv_luma(v, j+1, dir, dir);

  2348.                 }

  2349. 

  2350.                 ff_vc1_mc_4mv_chroma4(v, 0, 0, 0);

  2351.                 ff_vc1_mc_4mv_chroma4(v, 1, 1, 1);

  2352.             } else if (bmvtype == BMV_TYPE_INTERPOLATED) {

  2353.                 mvbp = v->twomvbp;

  2354.                 dmv_x = dmv_y = 0;

  2355.                 if (mvbp & 2)

  2356.                     get_mvdata_interlaced(v, &dmv_x, &dmv_y, 0);

  2357. 

  2358.                 ff_vc1_pred_mv_intfr(v, 0, dmv_x, dmv_y, 1, v->range_x, v->range_y, v->mb_type[0], 0);

  2359.                 ff_vc1_mc_1mv(v, 0);

  2360. 

  2361.                 dmv_x = dmv_y = 0;

  2362.                 if (mvbp & 1)

  2363.                     get_mvdata_interlaced(v, &dmv_x, &dmv_y, 0);

  2364. 

  2365.                 ff_vc1_pred_mv_intfr(v, 0, dmv_x, dmv_y, 1, v->range_x, v->range_y, v->mb_type[0], 1);

  2366.                 ff_vc1_interp_mc(v);

  2367.             } else if (twomv) {

  2368.                 dir = bmvtype == BMV_TYPE_BACKWARD;

  2369.                 dir2 = dir;

  2370.                 if (mvsw)

  2371.                     dir2 = !dir;

  2372.                 mvbp = v->twomvbp;

  2373.                 dmv_x = dmv_y = 0;

  2374.                 if (mvbp & 2)

  2375.                     get_mvdata_interlaced(v, &dmv_x, &dmv_y, 0);

  2376.                 ff_vc1_pred_mv_intfr(v, 0, dmv_x, dmv_y, 2, v->range_x, v->range_y, v->mb_type[0], dir);

  2377. 

  2378.                 dmv_x = dmv_y = 0;

  2379.                 if (mvbp & 1)

  2380.                     get_mvdata_interlaced(v, &dmv_x, &dmv_y, 0);

  2381.                 ff_vc1_pred_mv_intfr(v, 2, dmv_x, dmv_y, 2, v->range_x, v->range_y, v->mb_type[0], dir2);

  2382. 

  2383.                 if (mvsw) {

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

  2385.                         s->mv[dir][i+2][0] = s->mv[dir][i][0] = s->current_picture.motion_val[dir][s->block_index[i+2]][0] = s->current_picture.motion_val[dir][s->block_index[i]][0];

  2386.                         s->mv[dir][i+2][1] = s->mv[dir][i][1] = s->current_picture.motion_val[dir][s->block_index[i+2]][1] = s->current_picture.motion_val[dir][s->block_index[i]][1];

  2387.                         s->mv[dir2][i+2][0] = s->mv[dir2][i][0] = s->current_picture.motion_val[dir2][s->block_index[i]][0] = s->current_picture.motion_val[dir2][s->block_index[i+2]][0];

  2388.                         s->mv[dir2][i+2][1] = s->mv[dir2][i][1] = s->current_picture.motion_val[dir2][s->block_index[i]][1] = s->current_picture.motion_val[dir2][s->block_index[i+2]][1];

  2389.                     }

  2390.                 } else {

  2391.                     ff_vc1_pred_mv_intfr(v, 0, 0, 0, 2, v->range_x, v->range_y, v->mb_type[0], !dir);

  2392.                     ff_vc1_pred_mv_intfr(v, 2, 0, 0, 2, v->range_x, v->range_y, v->mb_type[0], !dir);

  2393.                 }

  2394. 

  2395.                 ff_vc1_mc_4mv_luma(v, 0, dir, 0);

  2396.                 ff_vc1_mc_4mv_luma(v, 1, dir, 0);

  2397.                 ff_vc1_mc_4mv_luma(v, 2, dir2, 0);

  2398.                 ff_vc1_mc_4mv_luma(v, 3, dir2, 0);

  2399.                 ff_vc1_mc_4mv_chroma4(v, dir, dir2, 0);

  2400.             } else {

  2401.                 dir = bmvtype == BMV_TYPE_BACKWARD;

  2402. 

  2403.                 mvbp = ff_vc1_mbmode_intfrp[0][idx_mbmode][2];

  2404.                 dmv_x = dmv_y = 0;

  2405.                 if (mvbp)

  2406.                     get_mvdata_interlaced(v, &dmv_x, &dmv_y, 0);

  2407. 

  2408.                 ff_vc1_pred_mv_intfr(v, 0, dmv_x, dmv_y, 1, v->range_x, v->range_y, v->mb_type[0], dir);

  2409.                 v->blk_mv_type[s->block_index[0]] = 1;

  2410.                 v->blk_mv_type[s->block_index[1]] = 1;

  2411.                 v->blk_mv_type[s->block_index[2]] = 1;

  2412.                 v->blk_mv_type[s->block_index[3]] = 1;

  2413.                 ff_vc1_pred_mv_intfr(v, 0, 0, 0, 2, v->range_x, v->range_y, 0, !dir);

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

  2415.                     s->mv[!dir][i+2][0] = s->mv[!dir][i][0] = s->current_picture.motion_val[!dir][s->block_index[i+2]][0] = s->current_picture.motion_val[!dir][s->block_index[i]][0];

  2416.                     s->mv[!dir][i+2][1] = s->mv[!dir][i][1] = s->current_picture.motion_val[!dir][s->block_index[i+2]][1] = s->current_picture.motion_val[!dir][s->block_index[i]][1];

  2417.                 }

  2418.                 ff_vc1_mc_1mv(v, dir);

  2419.             }

  2420. 

  2421.             if (cbp)

  2422.                 GET_MQUANT();  // p. 227

  2423.             s->current_picture.qscale_table[mb_pos] = mquant;

  2424.             if (!v->ttmbf && cbp)

  2425.                 ttmb = get_vlc2(gb, ff_vc1_ttmb_vlc[v->tt_index].table, VC1_TTMB_VLC_BITS, 2);

  2426.             for (i = 0; i < 6; i++) {

  2427.                 s->dc_val[0][s->block_index[i]] = 0;

  2428.                 dst_idx += i >> 2;

  2429.                 val = ((cbp >> (5 - i)) & 1);

  2430.                 if (!fieldtx)

  2431.                     off = (i & 4) ? 0 : ((i & 1) * 8 + (i & 2) * 4 * s->linesize);

  2432.                 else

  2433.                     off = (i & 4) ? 0 : ((i & 1) * 8 + ((i > 1) * s->linesize));

  2434.                 if (val) {

  2435.                     pat = vc1_decode_p_block(v, s->block[i], i, mquant, ttmb,

  2436.                                              first_block, s->dest[dst_idx] + off,

  2437.                                              (i & 4) ? s->uvlinesize : (s->linesize << fieldtx),

  2438.                                              CONFIG_GRAY && (i & 4) && (s->avctx->flags & AV_CODEC_FLAG_GRAY), &block_tt);

  2439.                     block_cbp |= pat << (i << 2);

  2440.                     if (!v->ttmbf && ttmb < 8)

  2441.                         ttmb = -1;

  2442.                     first_block = 0;

  2443.                 }

  2444.             }

  2445. 

  2446.         } else { // skipped

  2447.             dir = 0;

  2448.             for (i = 0; i < 6; i++) {

  2449.                 v->mb_type[0][s->block_index[i]] = 0;

  2450.                 s->dc_val[0][s->block_index[i]] = 0;

  2451.             }

  2452.             s->current_picture.mb_type[mb_pos]      = MB_TYPE_SKIP;

  2453.             s->current_picture.qscale_table[mb_pos] = 0;

  2454.             v->blk_mv_type[s->block_index[0]] = 0;

  2455.             v->blk_mv_type[s->block_index[1]] = 0;

  2456.             v->blk_mv_type[s->block_index[2]] = 0;

  2457.             v->blk_mv_type[s->block_index[3]] = 0;

  2458. 

  2459.             if (!direct) {

  2460.                 if (bmvtype == BMV_TYPE_INTERPOLATED) {

  2461.                     ff_vc1_pred_mv_intfr(v, 0, 0, 0, 1, v->range_x, v->range_y, v->mb_type[0], 0);

  2462.                     ff_vc1_pred_mv_intfr(v, 0, 0, 0, 1, v->range_x, v->range_y, v->mb_type[0], 1);

  2463.                 } else {

  2464.                     dir = bmvtype == BMV_TYPE_BACKWARD;

  2465.                     ff_vc1_pred_mv_intfr(v, 0, 0, 0, 1, v->range_x, v->range_y, v->mb_type[0], dir);

  2466.                     if (mvsw) {

  2467.                         int dir2 = dir;

  2468.                         if (mvsw)

  2469.                             dir2 = !dir;

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

  2471.                             s->mv[dir][i+2][0] = s->mv[dir][i][0] = s->current_picture.motion_val[dir][s->block_index[i+2]][0] = s->current_picture.motion_val[dir][s->block_index[i]][0];

  2472.                             s->mv[dir][i+2][1] = s->mv[dir][i][1] = s->current_picture.motion_val[dir][s->block_index[i+2]][1] = s->current_picture.motion_val[dir][s->block_index[i]][1];

  2473.                             s->mv[dir2][i+2][0] = s->mv[dir2][i][0] = s->current_picture.motion_val[dir2][s->block_index[i]][0] = s->current_picture.motion_val[dir2][s->block_index[i+2]][0];

  2474.                             s->mv[dir2][i+2][1] = s->mv[dir2][i][1] = s->current_picture.motion_val[dir2][s->block_index[i]][1] = s->current_picture.motion_val[dir2][s->block_index[i+2]][1];

  2475.                         }

  2476.                     } else {

  2477.                         v->blk_mv_type[s->block_index[0]] = 1;

  2478.                         v->blk_mv_type[s->block_index[1]] = 1;

  2479.                         v->blk_mv_type[s->block_index[2]] = 1;

  2480.                         v->blk_mv_type[s->block_index[3]] = 1;

  2481.                         ff_vc1_pred_mv_intfr(v, 0, 0, 0, 2, v->range_x, v->range_y, 0, !dir);

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

  2483.                             s->mv[!dir][i+2][0] = s->mv[!dir][i][0] = s->current_picture.motion_val[!dir][s->block_index[i+2]][0] = s->current_picture.motion_val[!dir][s->block_index[i]][0];

  2484.                             s->mv[!dir][i+2][1] = s->mv[!dir][i][1] = s->current_picture.motion_val[!dir][s->block_index[i+2]][1] = s->current_picture.motion_val[!dir][s->block_index[i]][1];

  2485.                         }

  2486.                     }

  2487.                 }

  2488.             }

  2489. 

  2490.             ff_vc1_mc_1mv(v, dir);

  2491.             if (direct || bmvtype == BMV_TYPE_INTERPOLATED) {

  2492.                 ff_vc1_interp_mc(v);

  2493.             }

  2494.         }

  2495.     }

  2496.     if (s->mb_x == s->mb_width - 1)

  2497.         memmove(v->is_intra_base, v->is_intra, sizeof(v->is_intra_base[0]) * s->mb_stride);

  2498.     v->cbp[s->mb_x]      = block_cbp;

  2499.     v->ttblk[s->mb_x]    = block_tt;

  2500.     return 0;

  2501. }

  2502. 

  2503. /** Decode blocks of I-frame

  2504.  */

  2505. static void vc1_decode_i_blocks(VC1Context *v)

  2506. {

  2507.     int k, j;

  2508.     MpegEncContext *s = &v->s;

  2509.     int cbp, val;

  2510.     uint8_t *coded_val;

  2511.     int mb_pos;

  2512. 

  2513.     /* select coding mode used for VLC tables selection */

  2514.     switch (v->y_ac_table_index) {

  2515.     case 0:

  2516.         v->codingset = (v->pqindex <= 8) ? CS_HIGH_RATE_INTRA : CS_LOW_MOT_INTRA;

  2517.         break;

  2518.     case 1:

  2519.         v->codingset = CS_HIGH_MOT_INTRA;

  2520.         break;

  2521.     case 2:

  2522.         v->codingset = CS_MID_RATE_INTRA;

  2523.         break;

  2524.     }

  2525. 

  2526.     switch (v->c_ac_table_index) {

  2527.     case 0:

  2528.         v->codingset2 = (v->pqindex <= 8) ? CS_HIGH_RATE_INTER : CS_LOW_MOT_INTER;

  2529.         break;

  2530.     case 1:

  2531.         v->codingset2 = CS_HIGH_MOT_INTER;

  2532.         break;

  2533.     case 2:

  2534.         v->codingset2 = CS_MID_RATE_INTER;

  2535.         break;

  2536.     }

  2537. 

  2538.     /* Set DC scale - y and c use the same */

  2539.     s->y_dc_scale = s->y_dc_scale_table[v->pq];

  2540.     s->c_dc_scale = s->c_dc_scale_table[v->pq];

  2541. 

  2542.     //do frame decode

  2543.     s->mb_x = s->mb_y = 0;

  2544.     s->mb_intra         = 1;

  2545.     s->first_slice_line = 1;

  2546.     for (s->mb_y = 0; s->mb_y < s->end_mb_y; s->mb_y++) {

  2547.         s->mb_x = 0;

  2548.         init_block_index(v);

  2549.         for (; s->mb_x < v->end_mb_x; s->mb_x++) {

  2550.             uint8_t *dst[6];

  2551.             ff_update_block_index(s);

  2552.             dst[0] = s->dest[0];

  2553.             dst[1] = dst[0] + 8;

  2554.             dst[2] = s->dest[0] + s->linesize * 8;

  2555.             dst[3] = dst[2] + 8;

  2556.             dst[4] = s->dest[1];

  2557.             dst[5] = s->dest[2];

  2558.             s->bdsp.clear_blocks(s->block[0]);

  2559.             mb_pos = s->mb_x + s->mb_y * s->mb_width;

  2560.             s->current_picture.mb_type[mb_pos]                     = MB_TYPE_INTRA;

  2561.             s->current_picture.qscale_table[mb_pos]                = v->pq;

  2562.             s->current_picture.motion_val[1][s->block_index[0]][0] = 0;

  2563.             s->current_picture.motion_val[1][s->block_index[0]][1] = 0;

  2564. 

  2565.             // do actual MB decoding and displaying

  2566.             cbp = get_vlc2(&v->s.gb, ff_msmp4_mb_i_vlc.table, MB_INTRA_VLC_BITS, 2);

  2567.             v->s.ac_pred = get_bits1(&v->s.gb);

  2568. 

  2569.             for (k = 0; k < 6; k++) {

  2570.                 val = ((cbp >> (5 - k)) & 1);

  2571. 

  2572.                 if (k < 4) {

  2573.                     int pred   = vc1_coded_block_pred(&v->s, k, &coded_val);

  2574.                     val        = val ^ pred;

  2575.                     *coded_val = val;

  2576.                 }

  2577.                 cbp |= val << (5 - k);

  2578. 

  2579.                 vc1_decode_i_block(v, s->block[k], k, val, (k < 4) ? v->codingset : v->codingset2);

  2580. 

  2581.                 if (CONFIG_GRAY && k > 3 && (s->avctx->flags & AV_CODEC_FLAG_GRAY))

  2582.                     continue;

  2583.                 v->vc1dsp.vc1_inv_trans_8x8(s->block[k]);

  2584.                 if (v->pq >= 9 && v->overlap) {

  2585.                     if (v->rangeredfrm)

  2586.                         for (j = 0; j < 64; j++)

  2587.                             s->block[k][j] <<= 1;

  2588.                     s->idsp.put_signed_pixels_clamped(s->block[k], dst[k],

  2589.                                                       k & 4 ? s->uvlinesize

  2590.                                                             : s->linesize);

  2591.                 } else {

  2592.                     if (v->rangeredfrm)

  2593.                         for (j = 0; j < 64; j++)

  2594.                             s->block[k][j] = (s->block[k][j] - 64) << 1;

  2595.                     s->idsp.put_pixels_clamped(s->block[k], dst[k],

  2596.                                                k & 4 ? s->uvlinesize

  2597.                                                      : s->linesize);

  2598.                 }

  2599.             }

  2600. 

  2601.             if (v->pq >= 9 && v->overlap) {

  2602.                 if (s->mb_x) {

  2603.                     v->vc1dsp.vc1_h_overlap(s->dest[0], s->linesize);

  2604.                     v->vc1dsp.vc1_h_overlap(s->dest[0] + 8 * s->linesize, s->linesize);

  2605.                     if (!CONFIG_GRAY || !(s->avctx->flags & AV_CODEC_FLAG_GRAY)) {

  2606.                         v->vc1dsp.vc1_h_overlap(s->dest[1], s->uvlinesize);

  2607.                         v->vc1dsp.vc1_h_overlap(s->dest[2], s->uvlinesize);

  2608.                     }

  2609.                 }

  2610.                 v->vc1dsp.vc1_h_overlap(s->dest[0] + 8, s->linesize);

  2611.                 v->vc1dsp.vc1_h_overlap(s->dest[0] + 8 * s->linesize + 8, s->linesize);

  2612.                 if (!s->first_slice_line) {

  2613.                     v->vc1dsp.vc1_v_overlap(s->dest[0], s->linesize);

  2614.                     v->vc1dsp.vc1_v_overlap(s->dest[0] + 8, s->linesize);

  2615.                     if (!CONFIG_GRAY || !(s->avctx->flags & AV_CODEC_FLAG_GRAY)) {

  2616.                         v->vc1dsp.vc1_v_overlap(s->dest[1], s->uvlinesize);

  2617.                         v->vc1dsp.vc1_v_overlap(s->dest[2], s->uvlinesize);

  2618.                     }

  2619.                 }

  2620.                 v->vc1dsp.vc1_v_overlap(s->dest[0] + 8 * s->linesize, s->linesize);

  2621.                 v->vc1dsp.vc1_v_overlap(s->dest[0] + 8 * s->linesize + 8, s->linesize);

  2622.             }

  2623.             if (v->s.loop_filter)

  2624.                 ff_vc1_loop_filter_iblk(v, v->pq);

  2625. 

  2626.             if (get_bits_count(&s->gb) > v->bits) {

  2627.                 ff_er_add_slice(&s->er, 0, 0, s->mb_x, s->mb_y, ER_MB_ERROR);

  2628.                 av_log(s->avctx, AV_LOG_ERROR, "Bits overconsumption: %i > %i\n",

  2629.                        get_bits_count(&s->gb), v->bits);

  2630.                 return;

  2631.             }

  2632.         }

  2633.         if (!v->s.loop_filter)

  2634.             ff_mpeg_draw_horiz_band(s, s->mb_y * 16, 16);

  2635.         else if (s->mb_y)

  2636.             ff_mpeg_draw_horiz_band(s, (s->mb_y - 1) * 16, 16);

  2637. 

  2638.         s->first_slice_line = 0;

  2639.     }

  2640.     if (v->s.loop_filter)

  2641.         ff_mpeg_draw_horiz_band(s, (s->end_mb_y - 1) * 16, 16);

  2642. 

  2643.     /* This is intentionally mb_height and not end_mb_y - unlike in advanced

  2644.      * profile, these only differ are when decoding MSS2 rectangles. */

  2645.     ff_er_add_slice(&s->er, 0, 0, s->mb_width - 1, s->mb_height - 1, ER_MB_END);

  2646. }

  2647. 

  2648. /** Decode blocks of I-frame for advanced profile

  2649.  */

  2650. static void vc1_decode_i_blocks_adv(VC1Context *v)

  2651. {

  2652.     int k;

  2653.     MpegEncContext *s = &v->s;

  2654.     int cbp, val;

  2655.     uint8_t *coded_val;

  2656.     int mb_pos;

  2657.     int mquant = v->pq;

  2658.     int mqdiff;

  2659.     GetBitContext *gb = &s->gb;

  2660. 

  2661.     /* select coding mode used for VLC tables selection */

  2662.     switch (v->y_ac_table_index) {

  2663.     case 0:

  2664.         v->codingset = (v->pqindex <= 8) ? CS_HIGH_RATE_INTRA : CS_LOW_MOT_INTRA;

  2665.         break;

  2666.     case 1:

  2667.         v->codingset = CS_HIGH_MOT_INTRA;

  2668.         break;

  2669.     case 2:

  2670.         v->codingset = CS_MID_RATE_INTRA;

  2671.         break;

  2672.     }

  2673. 

  2674.     switch (v->c_ac_table_index) {

  2675.     case 0:

  2676.         v->codingset2 = (v->pqindex <= 8) ? CS_HIGH_RATE_INTER : CS_LOW_MOT_INTER;

  2677.         break;

  2678.     case 1:

  2679.         v->codingset2 = CS_HIGH_MOT_INTER;

  2680.         break;

  2681.     case 2:

  2682.         v->codingset2 = CS_MID_RATE_INTER;

  2683.         break;

  2684.     }

  2685. 

  2686.     // do frame decode

  2687.     s->mb_x             = s->mb_y = 0;

  2688.     s->mb_intra         = 1;

  2689.     s->first_slice_line = 1;

  2690.     s->mb_y             = s->start_mb_y;

  2691.     if (s->start_mb_y) {

  2692.         s->mb_x = 0;

  2693.         init_block_index(v);

  2694.         memset(&s->coded_block[s->block_index[0] - s->b8_stride], 0,

  2695.                (1 + s->b8_stride) * sizeof(*s->coded_block));

  2696.     }

  2697.     for (; s->mb_y < s->end_mb_y; s->mb_y++) {

  2698.         s->mb_x = 0;

  2699.         init_block_index(v);

  2700.         for (;s->mb_x < s->mb_width; s->mb_x++) {

  2701.             int16_t (*block)[64] = v->block[v->cur_blk_idx];

  2702.             ff_update_block_index(s);

  2703.             s->bdsp.clear_blocks(block[0]);

  2704.             mb_pos = s->mb_x + s->mb_y * s->mb_stride;

  2705.             s->current_picture.mb_type[mb_pos + v->mb_off]                         = MB_TYPE_INTRA;

  2706.             s->current_picture.motion_val[1][s->block_index[0] + v->blocks_off][0] = 0;

  2707.             s->current_picture.motion_val[1][s->block_index[0] + v->blocks_off][1] = 0;

  2708. 

  2709.             // do actual MB decoding and displaying

  2710.             if (v->fieldtx_is_raw)

  2711.                 v->fieldtx_plane[mb_pos] = get_bits1(&v->s.gb);

  2712.             cbp = get_vlc2(&v->s.gb, ff_msmp4_mb_i_vlc.table, MB_INTRA_VLC_BITS, 2);

  2713.             if (v->acpred_is_raw)

  2714.                 v->s.ac_pred = get_bits1(&v->s.gb);

  2715.             else

  2716.                 v->s.ac_pred = v->acpred_plane[mb_pos];

  2717. 

  2718.             if (v->condover == CONDOVER_SELECT && v->overflg_is_raw)

  2719.                 v->over_flags_plane[mb_pos] = get_bits1(&v->s.gb);

  2720. 

  2721.             GET_MQUANT();

  2722. 

  2723.             s->current_picture.qscale_table[mb_pos] = mquant;

  2724.             /* Set DC scale - y and c use the same */

  2725.             s->y_dc_scale = s->y_dc_scale_table[mquant];

  2726.             s->c_dc_scale = s->c_dc_scale_table[mquant];

  2727. 

  2728.             for (k = 0; k < 6; k++) {

  2729.                 val = ((cbp >> (5 - k)) & 1);

  2730. 

  2731.                 if (k < 4) {

  2732.                     int pred   = vc1_coded_block_pred(&v->s, k, &coded_val);

  2733.                     val        = val ^ pred;

  2734.                     *coded_val = val;

  2735.                 }

  2736.                 cbp |= val << (5 - k);

  2737. 

  2738.                 v->a_avail = !s->first_slice_line || (k == 2 || k == 3);

  2739.                 v->c_avail = !!s->mb_x || (k == 1 || k == 3);

  2740. 

  2741.                 vc1_decode_i_block_adv(v, block[k], k, val,

  2742.                                        (k < 4) ? v->codingset : v->codingset2, mquant);

  2743. 

  2744.                 if (CONFIG_GRAY && k > 3 && (s->avctx->flags & AV_CODEC_FLAG_GRAY))

  2745.                     continue;

  2746.                 v->vc1dsp.vc1_inv_trans_8x8(block[k]);

  2747.             }

  2748. 

  2749.             ff_vc1_smooth_overlap_filter_iblk(v);

  2750.             vc1_put_signed_blocks_clamped(v);

  2751.             if (v->s.loop_filter)

  2752.                 ff_vc1_loop_filter_iblk_delayed(v, v->pq);

  2753. 

  2754.             if (get_bits_count(&s->gb) > v->bits) {

  2755.                 // TODO: may need modification to handle slice coding

  2756.                 ff_er_add_slice(&s->er, 0, s->start_mb_y, s->mb_x, s->mb_y, ER_MB_ERROR);

  2757.                 av_log(s->avctx, AV_LOG_ERROR, "Bits overconsumption: %i > %i\n",

  2758.                        get_bits_count(&s->gb), v->bits);

  2759.                 return;

  2760.             }

  2761.         }

  2762.         if (!v->s.loop_filter)

  2763.             ff_mpeg_draw_horiz_band(s, s->mb_y * 16, 16);

  2764.         else if (s->mb_y)

  2765.             ff_mpeg_draw_horiz_band(s, (s->mb_y-1) * 16, 16);

  2766.         s->first_slice_line = 0;

  2767.     }

  2768. 

  2769.     /* raw bottom MB row */

  2770.     s->mb_x = 0;

  2771.     init_block_index(v);

  2772.     for (; s->mb_x < s->mb_width; s->mb_x++) {

  2773.         ff_update_block_index(s);

  2774.         vc1_put_signed_blocks_clamped(v);

  2775.         if (v->s.loop_filter)

  2776.             ff_vc1_loop_filter_iblk_delayed(v, v->pq);

  2777.     }

  2778.     if (v->s.loop_filter)

  2779.         ff_mpeg_draw_horiz_band(s, (s->end_mb_y - 1) * 16, 16);

  2780.     ff_er_add_slice(&s->er, 0, s->start_mb_y << v->field_mode, s->mb_width - 1,

  2781.                     (s->end_mb_y << v->field_mode) - 1, ER_MB_END);

  2782. }

  2783. 

  2784. static void vc1_decode_p_blocks(VC1Context *v)

  2785. {

  2786.     MpegEncContext *s = &v->s;

  2787.     int apply_loop_filter;

  2788. 

  2789.     /* select coding mode used for VLC tables selection */

  2790.     switch (v->c_ac_table_index) {

  2791.     case 0:

  2792.         v->codingset = (v->pqindex <= 8) ? CS_HIGH_RATE_INTRA : CS_LOW_MOT_INTRA;

  2793.         break;

  2794.     case 1:

  2795.         v->codingset = CS_HIGH_MOT_INTRA;

  2796.         break;

  2797.     case 2:

  2798.         v->codingset = CS_MID_RATE_INTRA;

  2799.         break;

  2800.     }

  2801. 

  2802.     switch (v->c_ac_table_index) {

  2803.     case 0:

  2804.         v->codingset2 = (v->pqindex <= 8) ? CS_HIGH_RATE_INTER : CS_LOW_MOT_INTER;

  2805.         break;

  2806.     case 1:

  2807.         v->codingset2 = CS_HIGH_MOT_INTER;

  2808.         break;

  2809.     case 2:

  2810.         v->codingset2 = CS_MID_RATE_INTER;

  2811.         break;

  2812.     }

  2813. 

  2814.     apply_loop_filter   = s->loop_filter && !(s->avctx->skip_loop_filter >= AVDISCARD_NONKEY) &&

  2815.                           v->fcm == PROGRESSIVE;

  2816.     s->first_slice_line = 1;

  2817.     memset(v->cbp_base, 0, sizeof(v->cbp_base[0])*2*s->mb_stride);

  2818.     for (s->mb_y = s->start_mb_y; s->mb_y < s->end_mb_y; s->mb_y++) {

  2819.         s->mb_x = 0;

  2820.         init_block_index(v);

  2821.         for (; s->mb_x < s->mb_width; s->mb_x++) {

  2822.             ff_update_block_index(s);

  2823. 

  2824.             if (v->fcm == ILACE_FIELD)

  2825.                 vc1_decode_p_mb_intfi(v);

  2826.             else if (v->fcm == ILACE_FRAME)

  2827.                 vc1_decode_p_mb_intfr(v);

  2828.             else vc1_decode_p_mb(v);

  2829.             if (s->mb_y != s->start_mb_y && apply_loop_filter)

  2830.                 ff_vc1_apply_p_loop_filter(v);

  2831.             if (get_bits_count(&s->gb) > v->bits || get_bits_count(&s->gb) < 0) {

  2832.                 // TODO: may need modification to handle slice coding

  2833.                 ff_er_add_slice(&s->er, 0, s->start_mb_y, s->mb_x, s->mb_y, ER_MB_ERROR);

  2834.                 av_log(s->avctx, AV_LOG_ERROR, "Bits overconsumption: %i > %i at %ix%i\n",

  2835.                        get_bits_count(&s->gb), v->bits, s->mb_x, s->mb_y);

  2836.                 return;

  2837.             }

  2838.         }

  2839.         memmove(v->cbp_base,      v->cbp,      sizeof(v->cbp_base[0])      * s->mb_stride);

  2840.         memmove(v->ttblk_base,    v->ttblk,    sizeof(v->ttblk_base[0])    * s->mb_stride);

  2841.         memmove(v->is_intra_base, v->is_intra, sizeof(v->is_intra_base[0]) * s->mb_stride);

  2842.         memmove(v->luma_mv_base,  v->luma_mv,  sizeof(v->luma_mv_base[0])  * s->mb_stride);

  2843.         if (s->mb_y != s->start_mb_y)

  2844.             ff_mpeg_draw_horiz_band(s, (s->mb_y - 1) * 16, 16);

  2845.         s->first_slice_line = 0;

  2846.     }

  2847.     if (apply_loop_filter) {

  2848.         s->mb_x = 0;

  2849.         init_block_index(v);

  2850.         for (; s->mb_x < s->mb_width; s->mb_x++) {

  2851.             ff_update_block_index(s);

  2852.             ff_vc1_apply_p_loop_filter(v);

  2853.         }

  2854.     }

  2855.     if (s->end_mb_y >= s->start_mb_y)

  2856.         ff_mpeg_draw_horiz_band(s, (s->end_mb_y - 1) * 16, 16);

  2857.     ff_er_add_slice(&s->er, 0, s->start_mb_y << v->field_mode, s->mb_width - 1,

  2858.                     (s->end_mb_y << v->field_mode) - 1, ER_MB_END);

  2859. }

  2860. 

  2861. static void vc1_decode_b_blocks(VC1Context *v)

  2862. {

  2863.     MpegEncContext *s = &v->s;

  2864. 

  2865.     /* select coding mode used for VLC tables selection */

  2866.     switch (v->c_ac_table_index) {

  2867.     case 0:

  2868.         v->codingset = (v->pqindex <= 8) ? CS_HIGH_RATE_INTRA : CS_LOW_MOT_INTRA;

  2869.         break;

  2870.     case 1:

  2871.         v->codingset = CS_HIGH_MOT_INTRA;

  2872.         break;

  2873.     case 2:

  2874.         v->codingset = CS_MID_RATE_INTRA;

  2875.         break;

  2876.     }

  2877. 

  2878.     switch (v->c_ac_table_index) {

  2879.     case 0:

  2880.         v->codingset2 = (v->pqindex <= 8) ? CS_HIGH_RATE_INTER : CS_LOW_MOT_INTER;

  2881.         break;

  2882.     case 1:

  2883.         v->codingset2 = CS_HIGH_MOT_INTER;

  2884.         break;

  2885.     case 2:

  2886.         v->codingset2 = CS_MID_RATE_INTER;

  2887.         break;

  2888.     }

  2889. 

  2890.     s->first_slice_line = 1;

  2891.     for (s->mb_y = s->start_mb_y; s->mb_y < s->end_mb_y; s->mb_y++) {

  2892.         s->mb_x = 0;

  2893.         init_block_index(v);

  2894.         for (; s->mb_x < s->mb_width; s->mb_x++) {

  2895.             ff_update_block_index(s);

  2896. 

  2897.             if (v->fcm == ILACE_FIELD)

  2898.                 vc1_decode_b_mb_intfi(v);

  2899.             else if (v->fcm == ILACE_FRAME)

  2900.                 vc1_decode_b_mb_intfr(v);

  2901.             else

  2902.                 vc1_decode_b_mb(v);

  2903.             if (get_bits_count(&s->gb) > v->bits || get_bits_count(&s->gb) < 0) {

  2904.                 // TODO: may need modification to handle slice coding

  2905.                 ff_er_add_slice(&s->er, 0, s->start_mb_y, s->mb_x, s->mb_y, ER_MB_ERROR);

  2906.                 av_log(s->avctx, AV_LOG_ERROR, "Bits overconsumption: %i > %i at %ix%i\n",

  2907.                        get_bits_count(&s->gb), v->bits, s->mb_x, s->mb_y);

  2908.                 return;

  2909.             }

  2910.             if (v->s.loop_filter)

  2911.                 ff_vc1_loop_filter_iblk(v, v->pq);

  2912.         }

  2913.         if (!v->s.loop_filter)

  2914.             ff_mpeg_draw_horiz_band(s, s->mb_y * 16, 16);

  2915.         else if (s->mb_y)

  2916.             ff_mpeg_draw_horiz_band(s, (s->mb_y - 1) * 16, 16);

  2917.         s->first_slice_line = 0;

  2918.     }

  2919.     if (v->s.loop_filter)

  2920.         ff_mpeg_draw_horiz_band(s, (s->end_mb_y - 1) * 16, 16);

  2921.     ff_er_add_slice(&s->er, 0, s->start_mb_y << v->field_mode, s->mb_width - 1,

  2922.                     (s->end_mb_y << v->field_mode) - 1, ER_MB_END);

  2923. }

  2924. 

  2925. static void vc1_decode_skip_blocks(VC1Context *v)

  2926. {

  2927.     MpegEncContext *s = &v->s;

  2928. 

  2929.     if (!v->s.last_picture.f->data[0])

  2930.         return;

  2931. 

  2932.     ff_er_add_slice(&s->er, 0, s->start_mb_y, s->mb_width - 1, s->end_mb_y - 1, ER_MB_END);

  2933.     s->first_slice_line = 1;

  2934.     for (s->mb_y = s->start_mb_y; s->mb_y < s->end_mb_y; s->mb_y++) {

  2935.         s->mb_x = 0;

  2936.         init_block_index(v);

  2937.         ff_update_block_index(s);

  2938.         memcpy(s->dest[0], s->last_picture.f->data[0] + s->mb_y * 16 * s->linesize,   s->linesize   * 16);

  2939.         memcpy(s->dest[1], s->last_picture.f->data[1] + s->mb_y *  8 * s->uvlinesize, s->uvlinesize *  8);

  2940.         memcpy(s->dest[2], s->last_picture.f->data[2] + s->mb_y *  8 * s->uvlinesize, s->uvlinesize *  8);

  2941.         ff_mpeg_draw_horiz_band(s, s->mb_y * 16, 16);

  2942.         s->first_slice_line = 0;

  2943.     }

  2944.     s->pict_type = AV_PICTURE_TYPE_P;

  2945. }

  2946. 

  2947. void ff_vc1_decode_blocks(VC1Context *v)

  2948. {

  2949. 

  2950.     v->s.esc3_level_length = 0;

  2951.     if (v->x8_type) {

  2952.         ff_intrax8_decode_picture(&v->x8, &v->s.current_picture,

  2953.                                   &v->s.gb, &v->s.mb_x, &v->s.mb_y,

  2954.                                   2 * v->pq + v->halfpq, v->pq * !v->pquantizer,

  2955.                                   v->s.loop_filter, v->s.low_delay);

  2956. 

  2957.         ff_er_add_slice(&v->s.er, 0, 0,

  2958.                         (v->s.mb_x >> 1) - 1, (v->s.mb_y >> 1) - 1,

  2959.                         ER_MB_END);

  2960.     } else {

  2961.         v->cur_blk_idx     =  0;

  2962.         v->left_blk_idx    = -1;

  2963.         v->topleft_blk_idx =  1;

  2964.         v->top_blk_idx     =  2;

  2965.         switch (v->s.pict_type) {

  2966.         case AV_PICTURE_TYPE_I:

  2967.             if (v->profile == PROFILE_ADVANCED)

  2968.                 vc1_decode_i_blocks_adv(v);

  2969.             else

  2970.                 vc1_decode_i_blocks(v);

  2971.             break;

  2972.         case AV_PICTURE_TYPE_P:

  2973.             if (v->p_frame_skipped)

  2974.                 vc1_decode_skip_blocks(v);

  2975.             else

  2976.                 vc1_decode_p_blocks(v);

  2977.             break;

  2978.         case AV_PICTURE_TYPE_B:

  2979.             if (v->bi_type) {

  2980.                 if (v->profile == PROFILE_ADVANCED)

  2981.                     vc1_decode_i_blocks_adv(v);

  2982.                 else

  2983.                     vc1_decode_i_blocks(v);

  2984.             } else

  2985.                 vc1_decode_b_blocks(v);

  2986.             break;

  2987.         }

  2988.     }

  2989. }