libavcodec/vc1.c

   1 /*
   2  * VC-1 and WMV3 decoder common code
   3  * Copyright (c) 2006-2007 Konstantin Shishkov
   4  * Partly based on vc9.c (c) 2005 Anonymous, Alex Beregszaszi, Michael Niedermayer
   5  *
   6  * This file is part of Libav.
   7  *
   8  * Libav is free software; you can redistribute it and/or
   9  * modify it under the terms of the GNU Lesser General Public
  10  * License as published by the Free Software Foundation; either
  11  * version 2.1 of the License, or (at your option) any later version.
  12  *
  13  * Libav is distributed in the hope that it will be useful,
  14  * but WITHOUT ANY WARRANTY; without even the implied warranty of
  15  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  16  * Lesser General Public License for more details.
  17  *
  18  * You should have received a copy of the GNU Lesser General Public
  19  * License along with Libav; if not, write to the Free Software
  20  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
  21  */
  22
  23 /**
  24  * @file
  25  * VC-1 and WMV3 decoder common code
  26  *
  27  */
  28 #include "internal.h"
  29 #include "dsputil.h"
  30 #include "avcodec.h"
  31 #include "mpegvideo.h"
  32 #include "vc1.h"
  33 #include "vc1data.h"
  34 #include "msmpeg4data.h"
  35 #include "unary.h"
  36 #include "simple_idct.h"
  37
  38 #undef NDEBUG
  39 #include <assert.h>
  40
  41 /***********************************************************************/
  42 /**
  43  * @name VC-1 Bitplane decoding
  44  * @see 8.7, p56
  45  * @{
  46  */
  47
  48 /**
  49  * Imode types
  50  * @{
  51  */
  52 enum Imode {
  53     IMODE_RAW,
  54     IMODE_NORM2,
  55     IMODE_DIFF2,
  56     IMODE_NORM6,
  57     IMODE_DIFF6,
  58     IMODE_ROWSKIP,
  59     IMODE_COLSKIP
  60 };
  61 /** @} */ //imode defines
  62
  63 /** Decode rows by checking if they are skipped
  64  * @param plane Buffer to store decoded bits
  65  * @param[in] width Width of this buffer
  66  * @param[in] height Height of this buffer
  67  * @param[in] stride of this buffer
  68  */
  69 static void decode_rowskip(uint8_t* plane, int width, int height, int stride, GetBitContext *gb){
  70     int x, y;
  71
  72     for (y=0; y<height; y++){
  73         if (!get_bits1(gb)) //rowskip
  74             memset(plane, 0, width);
  75         else
  76             for (x=0; x<width; x++)
  77                 plane[x] = get_bits1(gb);
  78         plane += stride;
  79     }
  80 }
  81
  82 /** Decode columns by checking if they are skipped
  83  * @param plane Buffer to store decoded bits
  84  * @param[in] width Width of this buffer
  85  * @param[in] height Height of this buffer
  86  * @param[in] stride of this buffer
  87  * @todo FIXME: Optimize
  88  */
  89 static void decode_colskip(uint8_t* plane, int width, int height, int stride, GetBitContext *gb){
  90     int x, y;
  91
  92     for (x=0; x<width; x++){
  93         if (!get_bits1(gb)) //colskip
  94             for (y=0; y<height; y++)
  95                 plane[y*stride] = 0;
  96         else
  97             for (y=0; y<height; y++)
  98                 plane[y*stride] = get_bits1(gb);
  99         plane ++;
 100     }
 101 }
 102
 103 /** Decode a bitplane's bits
 104  * @param data bitplane where to store the decode bits
 105  * @param[out] raw_flag pointer to the flag indicating that this bitplane is not coded explicitly
 106  * @param v VC-1 context for bit reading and logging
 107  * @return Status
 108  * @todo FIXME: Optimize
 109  */
 110 static int bitplane_decoding(uint8_t* data, int *raw_flag, VC1Context *v)
 111 {
 112     GetBitContext *gb = &v->s.gb;
 113
 114     int imode, x, y, code, offset;
 115     uint8_t invert, *planep = data;
 116     int width, height, stride;
 117
 118     width = v->s.mb_width;
 119     height = v->s.mb_height;
 120     stride = v->s.mb_stride;
 121     invert = get_bits1(gb);
 122     imode = get_vlc2(gb, ff_vc1_imode_vlc.table, VC1_IMODE_VLC_BITS, 1);
 123
 124     *raw_flag = 0;
 125     switch (imode)
 126     {
 127     case IMODE_RAW:
 128         //Data is actually read in the MB layer (same for all tests == "raw")
 129         *raw_flag = 1; //invert ignored
 130         return invert;
 131     case IMODE_DIFF2:
 132     case IMODE_NORM2:
 133         if ((height * width) & 1)
 134         {
 135             *planep++ = get_bits1(gb);
 136             offset = 1;
 137         }
 138         else offset = 0;
 139         // decode bitplane as one long line
 140         for (y = offset; y < height * width; y += 2) {
 141             code = get_vlc2(gb, ff_vc1_norm2_vlc.table, VC1_NORM2_VLC_BITS, 1);
 142             *planep++ = code & 1;
 143             offset++;
 144             if(offset == width) {
 145                 offset = 0;
 146                 planep += stride - width;
 147             }
 148             *planep++ = code >> 1;
 149             offset++;
 150             if(offset == width) {
 151                 offset = 0;
 152                 planep += stride - width;
 153             }
 154         }
 155         break;
 156     case IMODE_DIFF6:
 157     case IMODE_NORM6:
 158         if(!(height % 3) && (width % 3)) { // use 2x3 decoding
 159             for(y = 0; y < height; y+= 3) {
 160                 for(x = width & 1; x < width; x += 2) {
 161                     code = get_vlc2(gb, ff_vc1_norm6_vlc.table, VC1_NORM6_VLC_BITS, 2);
 162                     if(code < 0){
 163                         av_log(v->s.avctx, AV_LOG_DEBUG, "invalid NORM-6 VLC\n");
 164                         return -1;
 165                     }
 166                     planep[x + 0] = (code >> 0) & 1;
 167                     planep[x + 1] = (code >> 1) & 1;
 168                     planep[x + 0 + stride] = (code >> 2) & 1;
 169                     planep[x + 1 + stride] = (code >> 3) & 1;
 170                     planep[x + 0 + stride * 2] = (code >> 4) & 1;
 171                     planep[x + 1 + stride * 2] = (code >> 5) & 1;
 172                 }
 173                 planep += stride * 3;
 174             }
 175             if(width & 1) decode_colskip(data, 1, height, stride, &v->s.gb);
 176         } else { // 3x2
 177             planep += (height & 1) * stride;
 178             for(y = height & 1; y < height; y += 2) {
 179                 for(x = width % 3; x < width; x += 3) {
 180                     code = get_vlc2(gb, ff_vc1_norm6_vlc.table, VC1_NORM6_VLC_BITS, 2);
 181                     if(code < 0){
 182                         av_log(v->s.avctx, AV_LOG_DEBUG, "invalid NORM-6 VLC\n");
 183                         return -1;
 184                     }
 185                     planep[x + 0] = (code >> 0) & 1;
 186                     planep[x + 1] = (code >> 1) & 1;
 187                     planep[x + 2] = (code >> 2) & 1;
 188                     planep[x + 0 + stride] = (code >> 3) & 1;
 189                     planep[x + 1 + stride] = (code >> 4) & 1;
 190                     planep[x + 2 + stride] = (code >> 5) & 1;
 191                 }
 192                 planep += stride * 2;
 193             }
 194             x = width % 3;
 195             if(x) decode_colskip(data  ,             x, height    , stride, &v->s.gb);
 196             if(height & 1) decode_rowskip(data+x, width - x, 1, stride, &v->s.gb);
 197         }
 198         break;
 199     case IMODE_ROWSKIP:
 200         decode_rowskip(data, width, height, stride, &v->s.gb);
 201         break;
 202     case IMODE_COLSKIP:
 203         decode_colskip(data, width, height, stride, &v->s.gb);
 204         break;
 205     default: break;
 206     }
 207
 208     /* Applying diff operator */
 209     if (imode == IMODE_DIFF2 || imode == IMODE_DIFF6)
 210     {
 211         planep = data;
 212         planep[0] ^= invert;
 213         for (x=1; x<width; x++)
 214             planep[x] ^= planep[x-1];
 215         for (y=1; y<height; y++)
 216         {
 217             planep += stride;
 218             planep[0] ^= planep[-stride];
 219             for (x=1; x<width; x++)
 220             {
 221                 if (planep[x-1] != planep[x-stride]) planep[x] ^= invert;
 222                 else                                 planep[x] ^= planep[x-1];
 223             }
 224         }
 225     }
 226     else if (invert)
 227     {
 228         planep = data;
 229         for (x=0; x<stride*height; x++) planep[x] = !planep[x]; //FIXME stride
 230     }
 231     return (imode<<1) + invert;
 232 }
 233
 234 /** @} */ //Bitplane group
 235
 236 /***********************************************************************/
 237 /** VOP Dquant decoding
 238  * @param v VC-1 Context
 239  */
 240 static int vop_dquant_decoding(VC1Context *v)
 241 {
 242     GetBitContext *gb = &v->s.gb;
 243     int pqdiff;
 244
 245     //variable size
 246     if (v->dquant == 2)
 247     {
 248         pqdiff = get_bits(gb, 3);
 249         if (pqdiff == 7) v->altpq = get_bits(gb, 5);
 250         else v->altpq = v->pq + pqdiff + 1;
 251     }
 252     else
 253     {
 254         v->dquantfrm = get_bits1(gb);
 255         if ( v->dquantfrm )
 256         {
 257             v->dqprofile = get_bits(gb, 2);
 258             switch (v->dqprofile)
 259             {
 260             case DQPROFILE_SINGLE_EDGE:
 261             case DQPROFILE_DOUBLE_EDGES:
 262                 v->dqsbedge = get_bits(gb, 2);
 263                 break;
 264             case DQPROFILE_ALL_MBS:
 265                 v->dqbilevel = get_bits1(gb);
 266                 if(!v->dqbilevel)
 267                     v->halfpq = 0;
 268             default: break; //Forbidden ?
 269             }
 270             if (v->dqbilevel || v->dqprofile != DQPROFILE_ALL_MBS)
 271             {
 272                 pqdiff = get_bits(gb, 3);
 273                 if (pqdiff == 7) v->altpq = get_bits(gb, 5);
 274                 else v->altpq = v->pq + pqdiff + 1;
 275             }
 276         }
 277     }
 278     return 0;
 279 }
 280
 281 static int decode_sequence_header_adv(VC1Context *v, GetBitContext *gb);
 282
 283 /**
 284  * Decode Simple/Main Profiles sequence header
 285  * @see Figure 7-8, p16-17
 286  * @param avctx Codec context
 287  * @param gb GetBit context initialized from Codec context extra_data
 288  * @return Status
 289  */
 290 int vc1_decode_sequence_header(AVCodecContext *avctx, VC1Context *v, GetBitContext *gb)
 291 {
 292     av_log(avctx, AV_LOG_DEBUG, "Header: %0X\n", show_bits(gb, 32));
 293     v->profile = get_bits(gb, 2);
 294     if (v->profile == PROFILE_COMPLEX)
 295     {
 296         av_log(avctx, AV_LOG_WARNING, "WMV3 Complex Profile is not fully supported\n");
 297     }
 298
 299     if (v->profile == PROFILE_ADVANCED)
 300     {
 301         v->zz_8x4 = ff_vc1_adv_progressive_8x4_zz;
 302         v->zz_4x8 = ff_vc1_adv_progressive_4x8_zz;
 303         return decode_sequence_header_adv(v, gb);
 304     }
 305     else
 306     {
 307         v->zz_8x4 = wmv2_scantableA;
 308         v->zz_4x8 = wmv2_scantableB;
 309         v->res_y411   = get_bits1(gb);
 310         v->res_sprite = get_bits1(gb);
 311         if (v->res_y411)
 312         {
 313             av_log(avctx, AV_LOG_ERROR,
 314                    "Old interlaced mode is not supported\n");
 315             return -1;
 316         }
 317     }
 318
 319     // (fps-2)/4 (->30)
 320     v->frmrtq_postproc = get_bits(gb, 3); //common
 321     // (bitrate-32kbps)/64kbps
 322     v->bitrtq_postproc = get_bits(gb, 5); //common
 323     v->s.loop_filter = get_bits1(gb); //common
 324     if(v->s.loop_filter == 1 && v->profile == PROFILE_SIMPLE)
 325     {
 326         av_log(avctx, AV_LOG_ERROR,
 327                "LOOPFILTER shall not be enabled in Simple Profile\n");
 328     }
 329     if(v->s.avctx->skip_loop_filter >= AVDISCARD_ALL)
 330         v->s.loop_filter = 0;
 331
 332     v->res_x8 = get_bits1(gb); //reserved
 333     v->multires = get_bits1(gb);
 334     v->res_fasttx = get_bits1(gb);
 335     if (!v->res_fasttx)
 336     {
 337         v->vc1dsp.vc1_inv_trans_8x8 = ff_simple_idct_8;
 338         v->vc1dsp.vc1_inv_trans_8x4 = ff_simple_idct84_add;
 339         v->vc1dsp.vc1_inv_trans_4x8 = ff_simple_idct48_add;
 340         v->vc1dsp.vc1_inv_trans_4x4 = ff_simple_idct44_add;
 341         v->vc1dsp.vc1_inv_trans_8x8_dc = ff_simple_idct_add_8;
 342         v->vc1dsp.vc1_inv_trans_8x4_dc = ff_simple_idct84_add;
 343         v->vc1dsp.vc1_inv_trans_4x8_dc = ff_simple_idct48_add;
 344         v->vc1dsp.vc1_inv_trans_4x4_dc = ff_simple_idct44_add;
 345     }
 346
 347     v->fastuvmc =  get_bits1(gb); //common
 348     if (!v->profile && !v->fastuvmc)
 349     {
 350         av_log(avctx, AV_LOG_ERROR,
 351                "FASTUVMC unavailable in Simple Profile\n");
 352         return -1;
 353     }
 354     v->extended_mv =  get_bits1(gb); //common
 355     if (!v->profile && v->extended_mv)
 356     {
 357         av_log(avctx, AV_LOG_ERROR,
 358                "Extended MVs unavailable in Simple Profile\n");
 359         return -1;
 360     }
 361     v->dquant =  get_bits(gb, 2); //common
 362     v->vstransform =  get_bits1(gb); //common
 363
 364     v->res_transtab = get_bits1(gb);
 365     if (v->res_transtab)
 366     {
 367         av_log(avctx, AV_LOG_ERROR,
 368                "1 for reserved RES_TRANSTAB is forbidden\n");
 369         return -1;
 370     }
 371
 372     v->overlap = get_bits1(gb); //common
 373
 374     v->s.resync_marker = get_bits1(gb);
 375     v->rangered = get_bits1(gb);
 376     if (v->rangered && v->profile == PROFILE_SIMPLE)
 377     {
 378         av_log(avctx, AV_LOG_INFO,
 379                "RANGERED should be set to 0 in Simple Profile\n");
 380     }
 381
 382     v->s.max_b_frames = avctx->max_b_frames = get_bits(gb, 3); //common
 383     v->quantizer_mode = get_bits(gb, 2); //common
 384
 385     v->finterpflag = get_bits1(gb); //common
 386
 387     if (v->res_sprite) {
 388         v->s.avctx->width  = v->s.avctx->coded_width  = get_bits(gb, 11);
 389         v->s.avctx->height = v->s.avctx->coded_height = get_bits(gb, 11);
 390         skip_bits(gb, 5); //frame rate
 391         v->res_x8 = get_bits1(gb);
 392         if (get_bits1(gb)) { // something to do with DC VLC selection
 393             av_log(avctx, AV_LOG_ERROR, "Unsupported sprite feature\n");
 394             return -1;
 395         }
 396         skip_bits(gb, 3); //slice code
 397         v->res_rtm_flag = 0;
 398     } else {
 399         v->res_rtm_flag = get_bits1(gb); //reserved
 400     }
 401     if (!v->res_rtm_flag)
 402     {
 403 //            av_log(avctx, AV_LOG_ERROR,
 404 //                   "0 for reserved RES_RTM_FLAG is forbidden\n");
 405         av_log(avctx, AV_LOG_ERROR,
 406                "Old WMV3 version detected, some frames may be decoded incorrectly\n");
 407         //return -1;
 408     }
 409     //TODO: figure out what they mean (always 0x402F)
 410     if(!v->res_fasttx) skip_bits(gb, 16);
 411     av_log(avctx, AV_LOG_DEBUG,
 412                "Profile %i:\nfrmrtq_postproc=%i, bitrtq_postproc=%i\n"
 413                "LoopFilter=%i, MultiRes=%i, FastUVMC=%i, Extended MV=%i\n"
 414                "Rangered=%i, VSTransform=%i, Overlap=%i, SyncMarker=%i\n"
 415                "DQuant=%i, Quantizer mode=%i, Max B frames=%i\n",
 416                v->profile, v->frmrtq_postproc, v->bitrtq_postproc,
 417                v->s.loop_filter, v->multires, v->fastuvmc, v->extended_mv,
 418                v->rangered, v->vstransform, v->overlap, v->s.resync_marker,
 419                v->dquant, v->quantizer_mode, avctx->max_b_frames
 420                );
 421     return 0;
 422 }
 423
 424 static int decode_sequence_header_adv(VC1Context *v, GetBitContext *gb)
 425 {
 426     v->res_rtm_flag = 1;
 427     v->level = get_bits(gb, 3);
 428     if(v->level >= 5)
 429     {
 430         av_log(v->s.avctx, AV_LOG_ERROR, "Reserved LEVEL %i\n",v->level);
 431     }
 432     v->chromaformat = get_bits(gb, 2);
 433     if (v->chromaformat != 1)
 434     {
 435         av_log(v->s.avctx, AV_LOG_ERROR,
 436                "Only 4:2:0 chroma format supported\n");
 437         return -1;
 438     }
 439
 440     // (fps-2)/4 (->30)
 441     v->frmrtq_postproc = get_bits(gb, 3); //common
 442     // (bitrate-32kbps)/64kbps
 443     v->bitrtq_postproc = get_bits(gb, 5); //common
 444     v->postprocflag = get_bits1(gb); //common
 445
 446     v->s.avctx->coded_width = (get_bits(gb, 12) + 1) << 1;
 447     v->s.avctx->coded_height = (get_bits(gb, 12) + 1) << 1;
 448     v->s.avctx->width = v->s.avctx->coded_width;
 449     v->s.avctx->height = v->s.avctx->coded_height;
 450     v->broadcast = get_bits1(gb);
 451     v->interlace = get_bits1(gb);
 452     v->tfcntrflag = get_bits1(gb);
 453     v->finterpflag = get_bits1(gb);
 454     skip_bits1(gb); // reserved
 455
 456     v->s.h_edge_pos = v->s.avctx->coded_width;
 457     v->s.v_edge_pos = v->s.avctx->coded_height;
 458
 459     av_log(v->s.avctx, AV_LOG_DEBUG,
 460                "Advanced Profile level %i:\nfrmrtq_postproc=%i, bitrtq_postproc=%i\n"
 461                "LoopFilter=%i, ChromaFormat=%i, Pulldown=%i, Interlace: %i\n"
 462                "TFCTRflag=%i, FINTERPflag=%i\n",
 463                v->level, v->frmrtq_postproc, v->bitrtq_postproc,
 464                v->s.loop_filter, v->chromaformat, v->broadcast, v->interlace,
 465                v->tfcntrflag, v->finterpflag
 466                );
 467
 468     v->psf = get_bits1(gb);
 469     if(v->psf) { //PsF, 6.1.13
 470         av_log(v->s.avctx, AV_LOG_ERROR, "Progressive Segmented Frame mode: not supported (yet)\n");
 471         return -1;
 472     }
 473     v->s.max_b_frames = v->s.avctx->max_b_frames = 7;
 474     if(get_bits1(gb)) { //Display Info - decoding is not affected by it
 475         int w, h, ar = 0;
 476         av_log(v->s.avctx, AV_LOG_DEBUG, "Display extended info:\n");
 477         v->s.avctx->width  = w = get_bits(gb, 14) + 1;
 478         v->s.avctx->height = h = get_bits(gb, 14) + 1;
 479         av_log(v->s.avctx, AV_LOG_DEBUG, "Display dimensions: %ix%i\n", w, h);
 480         if(get_bits1(gb))
 481             ar = get_bits(gb, 4);
 482         if(ar && ar < 14){
 483             v->s.avctx->sample_aspect_ratio = ff_vc1_pixel_aspect[ar];
 484         }else if(ar == 15){
 485             w = get_bits(gb, 8) + 1;
 486             h = get_bits(gb, 8) + 1;
 487             v->s.avctx->sample_aspect_ratio = (AVRational){w, h};
 488         }
 489         av_log(v->s.avctx, AV_LOG_DEBUG, "Aspect: %i:%i\n", v->s.avctx->sample_aspect_ratio.num, v->s.avctx->sample_aspect_ratio.den);
 490
 491         if(get_bits1(gb)){ //framerate stuff
 492             if(get_bits1(gb)) {
 493                 v->s.avctx->time_base.num = 32;
 494                 v->s.avctx->time_base.den = get_bits(gb, 16) + 1;
 495             } else {
 496                 int nr, dr;
 497                 nr = get_bits(gb, 8);
 498                 dr = get_bits(gb, 4);
 499                 if(nr && nr < 8 && dr && dr < 3){
 500                     v->s.avctx->time_base.num = ff_vc1_fps_dr[dr - 1];
 501                     v->s.avctx->time_base.den = ff_vc1_fps_nr[nr - 1] * 1000;
 502                 }
 503             }
 504         }
 505
 506         if(get_bits1(gb)){
 507             v->color_prim = get_bits(gb, 8);
 508             v->transfer_char = get_bits(gb, 8);
 509             v->matrix_coef = get_bits(gb, 8);
 510         }
 511     }
 512
 513     v->hrd_param_flag = get_bits1(gb);
 514     if(v->hrd_param_flag) {
 515         int i;
 516         v->hrd_num_leaky_buckets = get_bits(gb, 5);
 517         skip_bits(gb, 4); //bitrate exponent
 518         skip_bits(gb, 4); //buffer size exponent
 519         for(i = 0; i < v->hrd_num_leaky_buckets; i++) {
 520             skip_bits(gb, 16); //hrd_rate[n]
 521             skip_bits(gb, 16); //hrd_buffer[n]
 522         }
 523     }
 524     return 0;
 525 }
 526
 527 int vc1_decode_entry_point(AVCodecContext *avctx, VC1Context *v, GetBitContext *gb)
 528 {
 529     int i;
 530
 531     av_log(avctx, AV_LOG_DEBUG, "Entry point: %08X\n", show_bits_long(gb, 32));
 532     v->broken_link = get_bits1(gb);
 533     v->closed_entry = get_bits1(gb);
 534     v->panscanflag = get_bits1(gb);
 535     v->refdist_flag = get_bits1(gb);
 536     v->s.loop_filter = get_bits1(gb);
 537     v->fastuvmc = get_bits1(gb);
 538     v->extended_mv = get_bits1(gb);
 539     v->dquant = get_bits(gb, 2);
 540     v->vstransform = get_bits1(gb);
 541     v->overlap = get_bits1(gb);
 542     v->quantizer_mode = get_bits(gb, 2);
 543
 544     if(v->hrd_param_flag){
 545         for(i = 0; i < v->hrd_num_leaky_buckets; i++) {
 546             skip_bits(gb, 8); //hrd_full[n]
 547         }
 548     }
 549
 550     if(get_bits1(gb)){
 551         avctx->coded_width = (get_bits(gb, 12)+1)<<1;
 552         avctx->coded_height = (get_bits(gb, 12)+1)<<1;
 553     }
 554     if(v->extended_mv)
 555         v->extended_dmv = get_bits1(gb);
 556     if((v->range_mapy_flag = get_bits1(gb))) {
 557         av_log(avctx, AV_LOG_ERROR, "Luma scaling is not supported, expect wrong picture\n");
 558         v->range_mapy = get_bits(gb, 3);
 559     }
 560     if((v->range_mapuv_flag = get_bits1(gb))) {
 561         av_log(avctx, AV_LOG_ERROR, "Chroma scaling is not supported, expect wrong picture\n");
 562         v->range_mapuv = get_bits(gb, 3);
 563     }
 564
 565     av_log(avctx, AV_LOG_DEBUG, "Entry point info:\n"
 566         "BrokenLink=%i, ClosedEntry=%i, PanscanFlag=%i\n"
 567         "RefDist=%i, Postproc=%i, FastUVMC=%i, ExtMV=%i\n"
 568         "DQuant=%i, VSTransform=%i, Overlap=%i, Qmode=%i\n",
 569         v->broken_link, v->closed_entry, v->panscanflag, v->refdist_flag, v->s.loop_filter,
 570         v->fastuvmc, v->extended_mv, v->dquant, v->vstransform, v->overlap, v->quantizer_mode);
 571
 572     return 0;
 573 }
 574
 575 int vc1_parse_frame_header(VC1Context *v, GetBitContext* gb)
 576 {
 577     int pqindex, lowquant, status;
 578
 579     if(v->finterpflag) v->interpfrm = get_bits1(gb);
 580     skip_bits(gb, 2); //framecnt unused
 581     v->rangeredfrm = 0;
 582     if (v->rangered) v->rangeredfrm = get_bits1(gb);
 583     v->s.pict_type = get_bits1(gb);
 584     if (v->s.avctx->max_b_frames) {
 585         if (!v->s.pict_type) {
 586             if (get_bits1(gb)) v->s.pict_type = AV_PICTURE_TYPE_I;
 587             else v->s.pict_type = AV_PICTURE_TYPE_B;
 588         } else v->s.pict_type = AV_PICTURE_TYPE_P;
 589     } else v->s.pict_type = v->s.pict_type ? AV_PICTURE_TYPE_P : AV_PICTURE_TYPE_I;
 590
 591     v->bi_type = 0;
 592     if(v->s.pict_type == AV_PICTURE_TYPE_B) {
 593         v->bfraction_lut_index = get_vlc2(gb, ff_vc1_bfraction_vlc.table, VC1_BFRACTION_VLC_BITS, 1);
 594         v->bfraction = ff_vc1_bfraction_lut[v->bfraction_lut_index];
 595         if(v->bfraction == 0) {
 596             v->s.pict_type = AV_PICTURE_TYPE_BI;
 597         }
 598     }
 599     if(v->s.pict_type == AV_PICTURE_TYPE_I || v->s.pict_type == AV_PICTURE_TYPE_BI)
 600         skip_bits(gb, 7); // skip buffer fullness
 601
 602     if(v->parse_only)
 603         return 0;
 604
 605     /* calculate RND */
 606     if(v->s.pict_type == AV_PICTURE_TYPE_I || v->s.pict_type == AV_PICTURE_TYPE_BI)
 607         v->rnd = 1;
 608     if(v->s.pict_type == AV_PICTURE_TYPE_P)
 609         v->rnd ^= 1;
 610
 611     /* Quantizer stuff */
 612     pqindex = get_bits(gb, 5);
 613     if(!pqindex) return -1;
 614     if (v->quantizer_mode == QUANT_FRAME_IMPLICIT)
 615         v->pq = ff_vc1_pquant_table[0][pqindex];
 616     else
 617         v->pq = ff_vc1_pquant_table[1][pqindex];
 618
 619     v->pquantizer = 1;
 620     if (v->quantizer_mode == QUANT_FRAME_IMPLICIT)
 621         v->pquantizer = pqindex < 9;
 622     if (v->quantizer_mode == QUANT_NON_UNIFORM)
 623         v->pquantizer = 0;
 624     v->pqindex = pqindex;
 625     if (pqindex < 9) v->halfpq = get_bits1(gb);
 626     else v->halfpq = 0;
 627     if (v->quantizer_mode == QUANT_FRAME_EXPLICIT)
 628         v->pquantizer = get_bits1(gb);
 629     v->dquantfrm = 0;
 630     if (v->extended_mv == 1) v->mvrange = get_unary(gb, 0, 3);
 631     v->k_x = v->mvrange + 9 + (v->mvrange >> 1); //k_x can be 9 10 12 13
 632     v->k_y = v->mvrange + 8; //k_y can be 8 9 10 11
 633     v->range_x = 1 << (v->k_x - 1);
 634     v->range_y = 1 << (v->k_y - 1);
 635     if (v->multires && v->s.pict_type != AV_PICTURE_TYPE_B) v->respic = get_bits(gb, 2);
 636
 637     if(v->res_x8 && (v->s.pict_type == AV_PICTURE_TYPE_I || v->s.pict_type == AV_PICTURE_TYPE_BI)){
 638         v->x8_type = get_bits1(gb);
 639     }else v->x8_type = 0;
 640 //av_log(v->s.avctx, AV_LOG_INFO, "%c Frame: QP=[%i]%i (+%i/2) %i\n",
 641 //        (v->s.pict_type == AV_PICTURE_TYPE_P) ? 'P' : ((v->s.pict_type == AV_PICTURE_TYPE_I) ? 'I' : 'B'), pqindex, v->pq, v->halfpq, v->rangeredfrm);
 642
 643     if(v->s.pict_type == AV_PICTURE_TYPE_I || v->s.pict_type == AV_PICTURE_TYPE_P) v->use_ic = 0;
 644
 645     switch(v->s.pict_type) {
 646     case AV_PICTURE_TYPE_P:
 647         if (v->pq < 5) v->tt_index = 0;
 648         else if(v->pq < 13) v->tt_index = 1;
 649         else v->tt_index = 2;
 650
 651         lowquant = (v->pq > 12) ? 0 : 1;
 652         v->mv_mode = ff_vc1_mv_pmode_table[lowquant][get_unary(gb, 1, 4)];
 653         if (v->mv_mode == MV_PMODE_INTENSITY_COMP)
 654         {
 655             int scale, shift, i;
 656             v->mv_mode2 = ff_vc1_mv_pmode_table2[lowquant][get_unary(gb, 1, 3)];
 657             v->lumscale = get_bits(gb, 6);
 658             v->lumshift = get_bits(gb, 6);
 659             v->use_ic = 1;
 660             /* fill lookup tables for intensity compensation */
 661             if(!v->lumscale) {
 662                 scale = -64;
 663                 shift = (255 - v->lumshift * 2) << 6;
 664                 if(v->lumshift > 31)
 665                     shift += 128 << 6;
 666             } else {
 667                 scale = v->lumscale + 32;
 668                 if(v->lumshift > 31)
 669                     shift = (v->lumshift - 64) << 6;
 670                 else
 671                     shift = v->lumshift << 6;
 672             }
 673             for(i = 0; i < 256; i++) {
 674                 v->luty[i] = av_clip_uint8((scale * i + shift + 32) >> 6);
 675                 v->lutuv[i] = av_clip_uint8((scale * (i - 128) + 128*64 + 32) >> 6);
 676             }
 677         }
 678         if(v->mv_mode == MV_PMODE_1MV_HPEL || v->mv_mode == MV_PMODE_1MV_HPEL_BILIN)
 679             v->s.quarter_sample = 0;
 680         else if(v->mv_mode == MV_PMODE_INTENSITY_COMP) {
 681             if(v->mv_mode2 == MV_PMODE_1MV_HPEL || v->mv_mode2 == MV_PMODE_1MV_HPEL_BILIN)
 682                 v->s.quarter_sample = 0;
 683             else
 684                 v->s.quarter_sample = 1;
 685         } else
 686             v->s.quarter_sample = 1;
 687         v->s.mspel = !(v->mv_mode == MV_PMODE_1MV_HPEL_BILIN || (v->mv_mode == MV_PMODE_INTENSITY_COMP && v->mv_mode2 == MV_PMODE_1MV_HPEL_BILIN));
 688
 689         if ((v->mv_mode == MV_PMODE_INTENSITY_COMP &&
 690                  v->mv_mode2 == MV_PMODE_MIXED_MV)
 691                 || v->mv_mode == MV_PMODE_MIXED_MV)
 692         {
 693             status = bitplane_decoding(v->mv_type_mb_plane, &v->mv_type_is_raw, v);
 694             if (status < 0) return -1;
 695             av_log(v->s.avctx, AV_LOG_DEBUG, "MB MV Type plane encoding: "
 696                    "Imode: %i, Invert: %i\n", status>>1, status&1);
 697         } else {
 698             v->mv_type_is_raw = 0;
 699             memset(v->mv_type_mb_plane, 0, v->s.mb_stride * v->s.mb_height);
 700         }
 701         status = bitplane_decoding(v->s.mbskip_table, &v->skip_is_raw, v);
 702         if (status < 0) return -1;
 703         av_log(v->s.avctx, AV_LOG_DEBUG, "MB Skip plane encoding: "
 704                "Imode: %i, Invert: %i\n", status>>1, status&1);
 705
 706         /* Hopefully this is correct for P frames */
 707         v->s.mv_table_index = get_bits(gb, 2); //but using ff_vc1_ tables
 708         v->cbpcy_vlc = &ff_vc1_cbpcy_p_vlc[get_bits(gb, 2)];
 709
 710         if (v->dquant)
 711         {
 712             av_log(v->s.avctx, AV_LOG_DEBUG, "VOP DQuant info\n");
 713             vop_dquant_decoding(v);
 714         }
 715
 716         v->ttfrm = 0; //FIXME Is that so ?
 717         if (v->vstransform)
 718         {
 719             v->ttmbf = get_bits1(gb);
 720             if (v->ttmbf)
 721             {
 722                 v->ttfrm = ff_vc1_ttfrm_to_tt[get_bits(gb, 2)];
 723             }
 724         } else {
 725             v->ttmbf = 1;
 726             v->ttfrm = TT_8X8;
 727         }
 728         break;
 729     case AV_PICTURE_TYPE_B:
 730         if (v->pq < 5) v->tt_index = 0;
 731         else if(v->pq < 13) v->tt_index = 1;
 732         else v->tt_index = 2;
 733
 734         v->mv_mode = get_bits1(gb) ? MV_PMODE_1MV : MV_PMODE_1MV_HPEL_BILIN;
 735         v->s.quarter_sample = (v->mv_mode == MV_PMODE_1MV);
 736         v->s.mspel = v->s.quarter_sample;
 737
 738         status = bitplane_decoding(v->direct_mb_plane, &v->dmb_is_raw, v);
 739         if (status < 0) return -1;
 740         av_log(v->s.avctx, AV_LOG_DEBUG, "MB Direct Type plane encoding: "
 741                "Imode: %i, Invert: %i\n", status>>1, status&1);
 742         status = bitplane_decoding(v->s.mbskip_table, &v->skip_is_raw, v);
 743         if (status < 0) return -1;
 744         av_log(v->s.avctx, AV_LOG_DEBUG, "MB Skip plane encoding: "
 745                "Imode: %i, Invert: %i\n", status>>1, status&1);
 746
 747         v->s.mv_table_index = get_bits(gb, 2);
 748         v->cbpcy_vlc = &ff_vc1_cbpcy_p_vlc[get_bits(gb, 2)];
 749
 750         if (v->dquant)
 751         {
 752             av_log(v->s.avctx, AV_LOG_DEBUG, "VOP DQuant info\n");
 753             vop_dquant_decoding(v);
 754         }
 755
 756         v->ttfrm = 0;
 757         if (v->vstransform)
 758         {
 759             v->ttmbf = get_bits1(gb);
 760             if (v->ttmbf)
 761             {
 762                 v->ttfrm = ff_vc1_ttfrm_to_tt[get_bits(gb, 2)];
 763             }
 764         } else {
 765             v->ttmbf = 1;
 766             v->ttfrm = TT_8X8;
 767         }
 768         break;
 769     }
 770
 771     if(!v->x8_type)
 772     {
 773         /* AC Syntax */
 774         v->c_ac_table_index = decode012(gb);
 775         if (v->s.pict_type == AV_PICTURE_TYPE_I || v->s.pict_type == AV_PICTURE_TYPE_BI)
 776         {
 777             v->y_ac_table_index = decode012(gb);
 778         }
 779         /* DC Syntax */
 780         v->s.dc_table_index = get_bits1(gb);
 781     }
 782
 783     if(v->s.pict_type == AV_PICTURE_TYPE_BI) {
 784         v->s.pict_type = AV_PICTURE_TYPE_B;
 785         v->bi_type = 1;
 786     }
 787     return 0;
 788 }
 789
 790 int vc1_parse_frame_header_adv(VC1Context *v, GetBitContext* gb)
 791 {
 792     int pqindex, lowquant;
 793     int status;
 794
 795     v->p_frame_skipped = 0;
 796
 797     if(v->interlace){
 798         v->fcm = decode012(gb);
 799         if(v->fcm){
 800             if(!v->warn_interlaced++)
 801                 av_log(v->s.avctx, AV_LOG_ERROR, "Interlaced frames/fields support is not implemented\n");
 802             return -1;
 803         }
 804     }
 805     switch(get_unary(gb, 0, 4)) {
 806     case 0:
 807         v->s.pict_type = AV_PICTURE_TYPE_P;
 808         break;
 809     case 1:
 810         v->s.pict_type = AV_PICTURE_TYPE_B;
 811         break;
 812     case 2:
 813         v->s.pict_type = AV_PICTURE_TYPE_I;
 814         break;
 815     case 3:
 816         v->s.pict_type = AV_PICTURE_TYPE_BI;
 817         break;
 818     case 4:
 819         v->s.pict_type = AV_PICTURE_TYPE_P; // skipped pic
 820         v->p_frame_skipped = 1;
 821         return 0;
 822     }
 823     if(v->tfcntrflag)
 824         skip_bits(gb, 8);
 825     if(v->broadcast) {
 826         if(!v->interlace || v->psf) {
 827             v->rptfrm = get_bits(gb, 2);
 828         } else {
 829             v->tff = get_bits1(gb);
 830             v->rptfrm = get_bits1(gb);
 831         }
 832     }
 833     if(v->panscanflag) {
 834         av_log_missing_feature(v->s.avctx, "Pan-scan", 0);
 835         //...
 836     }
 837     v->rnd = get_bits1(gb);
 838     if(v->interlace)
 839         v->uvsamp = get_bits1(gb);
 840     if(v->finterpflag) v->interpfrm = get_bits1(gb);
 841     if(v->s.pict_type == AV_PICTURE_TYPE_B) {
 842         v->bfraction_lut_index = get_vlc2(gb, ff_vc1_bfraction_vlc.table, VC1_BFRACTION_VLC_BITS, 1);
 843         v->bfraction = ff_vc1_bfraction_lut[v->bfraction_lut_index];
 844         if(v->bfraction == 0) {
 845             v->s.pict_type = AV_PICTURE_TYPE_BI; /* XXX: should not happen here */
 846         }
 847     }
 848     pqindex = get_bits(gb, 5);
 849     if(!pqindex) return -1;
 850     v->pqindex = pqindex;
 851     if (v->quantizer_mode == QUANT_FRAME_IMPLICIT)
 852         v->pq = ff_vc1_pquant_table[0][pqindex];
 853     else
 854         v->pq = ff_vc1_pquant_table[1][pqindex];
 855
 856     v->pquantizer = 1;
 857     if (v->quantizer_mode == QUANT_FRAME_IMPLICIT)
 858         v->pquantizer = pqindex < 9;
 859     if (v->quantizer_mode == QUANT_NON_UNIFORM)
 860         v->pquantizer = 0;
 861     v->pqindex = pqindex;
 862     if (pqindex < 9) v->halfpq = get_bits1(gb);
 863     else v->halfpq = 0;
 864     if (v->quantizer_mode == QUANT_FRAME_EXPLICIT)
 865         v->pquantizer = get_bits1(gb);
 866     if(v->postprocflag)
 867         v->postproc = get_bits(gb, 2);
 868
 869     if(v->s.pict_type == AV_PICTURE_TYPE_I || v->s.pict_type == AV_PICTURE_TYPE_P) v->use_ic = 0;
 870
 871     if(v->parse_only)
 872         return 0;
 873
 874     switch(v->s.pict_type) {
 875     case AV_PICTURE_TYPE_I:
 876     case AV_PICTURE_TYPE_BI:
 877         status = bitplane_decoding(v->acpred_plane, &v->acpred_is_raw, v);
 878         if (status < 0) return -1;
 879         av_log(v->s.avctx, AV_LOG_DEBUG, "ACPRED plane encoding: "
 880                 "Imode: %i, Invert: %i\n", status>>1, status&1);
 881         v->condover = CONDOVER_NONE;
 882         if(v->overlap && v->pq <= 8) {
 883             v->condover = decode012(gb);
 884             if(v->condover == CONDOVER_SELECT) {
 885                 status = bitplane_decoding(v->over_flags_plane, &v->overflg_is_raw, v);
 886                 if (status < 0) return -1;
 887                 av_log(v->s.avctx, AV_LOG_DEBUG, "CONDOVER plane encoding: "
 888                         "Imode: %i, Invert: %i\n", status>>1, status&1);
 889             }
 890         }
 891         break;
 892     case AV_PICTURE_TYPE_P:
 893         if (v->extended_mv) v->mvrange = get_unary(gb, 0, 3);
 894         else v->mvrange = 0;
 895         v->k_x = v->mvrange + 9 + (v->mvrange >> 1); //k_x can be 9 10 12 13
 896         v->k_y = v->mvrange + 8; //k_y can be 8 9 10 11
 897         v->range_x = 1 << (v->k_x - 1);
 898         v->range_y = 1 << (v->k_y - 1);
 899
 900         if (v->pq < 5) v->tt_index = 0;
 901         else if(v->pq < 13) v->tt_index = 1;
 902         else v->tt_index = 2;
 903
 904         lowquant = (v->pq > 12) ? 0 : 1;
 905         v->mv_mode = ff_vc1_mv_pmode_table[lowquant][get_unary(gb, 1, 4)];
 906         if (v->mv_mode == MV_PMODE_INTENSITY_COMP)
 907         {
 908             int scale, shift, i;
 909             v->mv_mode2 = ff_vc1_mv_pmode_table2[lowquant][get_unary(gb, 1, 3)];
 910             v->lumscale = get_bits(gb, 6);
 911             v->lumshift = get_bits(gb, 6);
 912             /* fill lookup tables for intensity compensation */
 913             if(!v->lumscale) {
 914                 scale = -64;
 915                 shift = (255 - v->lumshift * 2) << 6;
 916                 if(v->lumshift > 31)
 917                     shift += 128 << 6;
 918             } else {
 919                 scale = v->lumscale + 32;
 920                 if(v->lumshift > 31)
 921                     shift = (v->lumshift - 64) << 6;
 922                 else
 923                     shift = v->lumshift << 6;
 924             }
 925             for(i = 0; i < 256; i++) {
 926                 v->luty[i] = av_clip_uint8((scale * i + shift + 32) >> 6);
 927                 v->lutuv[i] = av_clip_uint8((scale * (i - 128) + 128*64 + 32) >> 6);
 928             }
 929             v->use_ic = 1;
 930         }
 931         if(v->mv_mode == MV_PMODE_1MV_HPEL || v->mv_mode == MV_PMODE_1MV_HPEL_BILIN)
 932             v->s.quarter_sample = 0;
 933         else if(v->mv_mode == MV_PMODE_INTENSITY_COMP) {
 934             if(v->mv_mode2 == MV_PMODE_1MV_HPEL || v->mv_mode2 == MV_PMODE_1MV_HPEL_BILIN)
 935                 v->s.quarter_sample = 0;
 936             else
 937                 v->s.quarter_sample = 1;
 938         } else
 939             v->s.quarter_sample = 1;
 940         v->s.mspel = !(v->mv_mode == MV_PMODE_1MV_HPEL_BILIN || (v->mv_mode == MV_PMODE_INTENSITY_COMP && v->mv_mode2 == MV_PMODE_1MV_HPEL_BILIN));
 941
 942         if ((v->mv_mode == MV_PMODE_INTENSITY_COMP &&
 943                  v->mv_mode2 == MV_PMODE_MIXED_MV)
 944                 || v->mv_mode == MV_PMODE_MIXED_MV)
 945         {
 946             status = bitplane_decoding(v->mv_type_mb_plane, &v->mv_type_is_raw, v);
 947             if (status < 0) return -1;
 948             av_log(v->s.avctx, AV_LOG_DEBUG, "MB MV Type plane encoding: "
 949                    "Imode: %i, Invert: %i\n", status>>1, status&1);
 950         } else {
 951             v->mv_type_is_raw = 0;
 952             memset(v->mv_type_mb_plane, 0, v->s.mb_stride * v->s.mb_height);
 953         }
 954         status = bitplane_decoding(v->s.mbskip_table, &v->skip_is_raw, v);
 955         if (status < 0) return -1;
 956         av_log(v->s.avctx, AV_LOG_DEBUG, "MB Skip plane encoding: "
 957                "Imode: %i, Invert: %i\n", status>>1, status&1);
 958
 959         /* Hopefully this is correct for P frames */
 960         v->s.mv_table_index = get_bits(gb, 2); //but using ff_vc1_ tables
 961         v->cbpcy_vlc = &ff_vc1_cbpcy_p_vlc[get_bits(gb, 2)];
 962         if (v->dquant)
 963         {
 964             av_log(v->s.avctx, AV_LOG_DEBUG, "VOP DQuant info\n");
 965             vop_dquant_decoding(v);
 966         }
 967
 968         v->ttfrm = 0; //FIXME Is that so ?
 969         if (v->vstransform)
 970         {
 971             v->ttmbf = get_bits1(gb);
 972             if (v->ttmbf)
 973             {
 974                 v->ttfrm = ff_vc1_ttfrm_to_tt[get_bits(gb, 2)];
 975             }
 976         } else {
 977             v->ttmbf = 1;
 978             v->ttfrm = TT_8X8;
 979         }
 980         break;
 981     case AV_PICTURE_TYPE_B:
 982         if (v->extended_mv) v->mvrange = get_unary(gb, 0, 3);
 983         else v->mvrange = 0;
 984         v->k_x = v->mvrange + 9 + (v->mvrange >> 1); //k_x can be 9 10 12 13
 985         v->k_y = v->mvrange + 8; //k_y can be 8 9 10 11
 986         v->range_x = 1 << (v->k_x - 1);
 987         v->range_y = 1 << (v->k_y - 1);
 988
 989         if (v->pq < 5) v->tt_index = 0;
 990         else if(v->pq < 13) v->tt_index = 1;
 991         else v->tt_index = 2;
 992
 993         v->mv_mode = get_bits1(gb) ? MV_PMODE_1MV : MV_PMODE_1MV_HPEL_BILIN;
 994         v->s.quarter_sample = (v->mv_mode == MV_PMODE_1MV);
 995         v->s.mspel = v->s.quarter_sample;
 996
 997         status = bitplane_decoding(v->direct_mb_plane, &v->dmb_is_raw, v);
 998         if (status < 0) return -1;
 999         av_log(v->s.avctx, AV_LOG_DEBUG, "MB Direct Type plane encoding: "
1000                "Imode: %i, Invert: %i\n", status>>1, status&1);
1001         status = bitplane_decoding(v->s.mbskip_table, &v->skip_is_raw, v);
1002         if (status < 0) return -1;
1003         av_log(v->s.avctx, AV_LOG_DEBUG, "MB Skip plane encoding: "
1004                "Imode: %i, Invert: %i\n", status>>1, status&1);
1005
1006         v->s.mv_table_index = get_bits(gb, 2);
1007         v->cbpcy_vlc = &ff_vc1_cbpcy_p_vlc[get_bits(gb, 2)];
1008
1009         if (v->dquant)
1010         {
1011             av_log(v->s.avctx, AV_LOG_DEBUG, "VOP DQuant info\n");
1012             vop_dquant_decoding(v);
1013         }
1014
1015         v->ttfrm = 0;
1016         if (v->vstransform)
1017         {
1018             v->ttmbf = get_bits1(gb);
1019             if (v->ttmbf)
1020             {
1021                 v->ttfrm = ff_vc1_ttfrm_to_tt[get_bits(gb, 2)];
1022             }
1023         } else {
1024             v->ttmbf = 1;
1025             v->ttfrm = TT_8X8;
1026         }
1027         break;
1028     }
1029
1030     /* AC Syntax */
1031     v->c_ac_table_index = decode012(gb);
1032     if (v->s.pict_type == AV_PICTURE_TYPE_I || v->s.pict_type == AV_PICTURE_TYPE_BI)
1033     {
1034         v->y_ac_table_index = decode012(gb);
1035     }
1036     /* DC Syntax */
1037     v->s.dc_table_index = get_bits1(gb);
1038     if ((v->s.pict_type == AV_PICTURE_TYPE_I || v->s.pict_type == AV_PICTURE_TYPE_BI) && v->dquant) {
1039         av_log(v->s.avctx, AV_LOG_DEBUG, "VOP DQuant info\n");
1040         vop_dquant_decoding(v);
1041     }
1042
1043     v->bi_type = 0;
1044     if(v->s.pict_type == AV_PICTURE_TYPE_BI) {
1045         v->s.pict_type = AV_PICTURE_TYPE_B;
1046         v->bi_type = 1;
1047     }
1048     return 0;
1049 }