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 FFmpeg.
   7  *
   8  * FFmpeg is free software; you can redistribute it and/or
   9  * modify it under the terms of the GNU Lesser General Public
  10  * License as published by the Free Software Foundation; either
  11  * version 2.1 of the License, or (at your option) any later version.
  12  *
  13  * FFmpeg is distributed in the hope that it will be useful,
  14  * but WITHOUT ANY WARRANTY; without even the implied warranty of
  15  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  16  * Lesser General Public License for more details.
  17  *
  18  * You should have received a copy of the GNU Lesser General Public
  19  * License along with FFmpeg; if not, write to the Free Software
  20  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
  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         if (v->res_sprite) {
 318             av_log(avctx, AV_LOG_ERROR, "WMVP is not fully supported\n");
 319         }
 320     }
 321
 322     // (fps-2)/4 (->30)
 323     v->frmrtq_postproc = get_bits(gb, 3); //common
 324     // (bitrate-32kbps)/64kbps
 325     v->bitrtq_postproc = get_bits(gb, 5); //common
 326     v->s.loop_filter = get_bits1(gb); //common
 327     if(v->s.loop_filter == 1 && v->profile == PROFILE_SIMPLE)
 328     {
 329         av_log(avctx, AV_LOG_ERROR,
 330                "LOOPFILTER shall not be enabled in Simple Profile\n");
 331     }
 332     if(v->s.avctx->skip_loop_filter >= AVDISCARD_ALL)
 333         v->s.loop_filter = 0;
 334
 335     v->res_x8 = get_bits1(gb); //reserved
 336     v->multires = get_bits1(gb);
 337     v->res_fasttx = get_bits1(gb);
 338     if (!v->res_fasttx)
 339     {
 340         v->vc1dsp.vc1_inv_trans_8x8 = ff_simple_idct_8;
 341         v->vc1dsp.vc1_inv_trans_8x4 = ff_simple_idct84_add;
 342         v->vc1dsp.vc1_inv_trans_4x8 = ff_simple_idct48_add;
 343         v->vc1dsp.vc1_inv_trans_4x4 = ff_simple_idct44_add;
 344         v->vc1dsp.vc1_inv_trans_8x8_dc = ff_simple_idct_add_8;
 345         v->vc1dsp.vc1_inv_trans_8x4_dc = ff_simple_idct84_add;
 346         v->vc1dsp.vc1_inv_trans_4x8_dc = ff_simple_idct48_add;
 347         v->vc1dsp.vc1_inv_trans_4x4_dc = ff_simple_idct44_add;
 348     }
 349
 350     v->fastuvmc =  get_bits1(gb); //common
 351     if (!v->profile && !v->fastuvmc)
 352     {
 353         av_log(avctx, AV_LOG_ERROR,
 354                "FASTUVMC unavailable in Simple Profile\n");
 355         return -1;
 356     }
 357     v->extended_mv =  get_bits1(gb); //common
 358     if (!v->profile && v->extended_mv)
 359     {
 360         av_log(avctx, AV_LOG_ERROR,
 361                "Extended MVs unavailable in Simple Profile\n");
 362         return -1;
 363     }
 364     v->dquant =  get_bits(gb, 2); //common
 365     v->vstransform =  get_bits1(gb); //common
 366
 367     v->res_transtab = get_bits1(gb);
 368     if (v->res_transtab)
 369     {
 370         av_log(avctx, AV_LOG_ERROR,
 371                "1 for reserved RES_TRANSTAB is forbidden\n");
 372         return -1;
 373     }
 374
 375     v->overlap = get_bits1(gb); //common
 376
 377     v->s.resync_marker = get_bits1(gb);
 378     v->rangered = get_bits1(gb);
 379     if (v->rangered && v->profile == PROFILE_SIMPLE)
 380     {
 381         av_log(avctx, AV_LOG_INFO,
 382                "RANGERED should be set to 0 in Simple Profile\n");
 383     }
 384
 385     v->s.max_b_frames = avctx->max_b_frames = get_bits(gb, 3); //common
 386     v->quantizer_mode = get_bits(gb, 2); //common
 387
 388     v->finterpflag = get_bits1(gb); //common
 389
 390     if (v->res_sprite) {
 391         v->s.avctx->width  = v->s.avctx->coded_width  = get_bits(gb, 11);
 392         v->s.avctx->height = v->s.avctx->coded_height = get_bits(gb, 11);
 393         skip_bits(gb, 5); //frame rate
 394         v->res_x8 = get_bits1(gb);
 395         if (get_bits1(gb)) { // something to do with DC VLC selection
 396             av_log(avctx, AV_LOG_ERROR, "Unsupported sprite feature\n");
 397             return -1;
 398         }
 399         skip_bits(gb, 3); //slice code
 400         v->res_rtm_flag = 0;
 401     } else {
 402         v->res_rtm_flag = get_bits1(gb); //reserved
 403     }
 404     if (!v->res_rtm_flag)
 405     {
 406 //            av_log(avctx, AV_LOG_ERROR,
 407 //                   "0 for reserved RES_RTM_FLAG is forbidden\n");
 408         av_log(avctx, AV_LOG_ERROR,
 409                "Old WMV3 version detected, some frames may be decoded incorrectly\n");
 410         //return -1;
 411     }
 412     //TODO: figure out what they mean (always 0x402F)
 413     if(!v->res_fasttx) skip_bits(gb, 16);
 414     av_log(avctx, AV_LOG_DEBUG,
 415                "Profile %i:\nfrmrtq_postproc=%i, bitrtq_postproc=%i\n"
 416                "LoopFilter=%i, MultiRes=%i, FastUVMC=%i, Extended MV=%i\n"
 417                "Rangered=%i, VSTransform=%i, Overlap=%i, SyncMarker=%i\n"
 418                "DQuant=%i, Quantizer mode=%i, Max B frames=%i\n",
 419                v->profile, v->frmrtq_postproc, v->bitrtq_postproc,
 420                v->s.loop_filter, v->multires, v->fastuvmc, v->extended_mv,
 421                v->rangered, v->vstransform, v->overlap, v->s.resync_marker,
 422                v->dquant, v->quantizer_mode, avctx->max_b_frames
 423                );
 424     return 0;
 425 }
 426
 427 static int decode_sequence_header_adv(VC1Context *v, GetBitContext *gb)
 428 {
 429     v->res_rtm_flag = 1;
 430     v->level = get_bits(gb, 3);
 431     if(v->level >= 5)
 432     {
 433         av_log(v->s.avctx, AV_LOG_ERROR, "Reserved LEVEL %i\n",v->level);
 434     }
 435     v->chromaformat = get_bits(gb, 2);
 436     if (v->chromaformat != 1)
 437     {
 438         av_log(v->s.avctx, AV_LOG_ERROR,
 439                "Only 4:2:0 chroma format supported\n");
 440         return -1;
 441     }
 442
 443     // (fps-2)/4 (->30)
 444     v->frmrtq_postproc = get_bits(gb, 3); //common
 445     // (bitrate-32kbps)/64kbps
 446     v->bitrtq_postproc = get_bits(gb, 5); //common
 447     v->postprocflag = get_bits1(gb); //common
 448
 449     v->s.avctx->coded_width = (get_bits(gb, 12) + 1) << 1;
 450     v->s.avctx->coded_height = (get_bits(gb, 12) + 1) << 1;
 451     v->s.avctx->width = v->s.avctx->coded_width;
 452     v->s.avctx->height = v->s.avctx->coded_height;
 453     v->broadcast = get_bits1(gb);
 454     v->interlace = get_bits1(gb);
 455     v->tfcntrflag = get_bits1(gb);
 456     v->finterpflag = get_bits1(gb);
 457     skip_bits1(gb); // reserved
 458
 459     v->s.h_edge_pos = v->s.avctx->coded_width;
 460     v->s.v_edge_pos = v->s.avctx->coded_height;
 461
 462     av_log(v->s.avctx, AV_LOG_DEBUG,
 463                "Advanced Profile level %i:\nfrmrtq_postproc=%i, bitrtq_postproc=%i\n"
 464                "LoopFilter=%i, ChromaFormat=%i, Pulldown=%i, Interlace: %i\n"
 465                "TFCTRflag=%i, FINTERPflag=%i\n",
 466                v->level, v->frmrtq_postproc, v->bitrtq_postproc,
 467                v->s.loop_filter, v->chromaformat, v->broadcast, v->interlace,
 468                v->tfcntrflag, v->finterpflag
 469                );
 470
 471     v->psf = get_bits1(gb);
 472     if(v->psf) { //PsF, 6.1.13
 473         av_log(v->s.avctx, AV_LOG_ERROR, "Progressive Segmented Frame mode: not supported (yet)\n");
 474         return -1;
 475     }
 476     v->s.max_b_frames = v->s.avctx->max_b_frames = 7;
 477     if(get_bits1(gb)) { //Display Info - decoding is not affected by it
 478         int w, h, ar = 0;
 479         av_log(v->s.avctx, AV_LOG_DEBUG, "Display extended info:\n");
 480         v->s.avctx->width  = w = get_bits(gb, 14) + 1;
 481         v->s.avctx->height = h = get_bits(gb, 14) + 1;
 482         av_log(v->s.avctx, AV_LOG_DEBUG, "Display dimensions: %ix%i\n", w, h);
 483         if(get_bits1(gb))
 484             ar = get_bits(gb, 4);
 485         if(ar && ar < 14){
 486             v->s.avctx->sample_aspect_ratio = ff_vc1_pixel_aspect[ar];
 487         }else if(ar == 15){
 488             w = get_bits(gb, 8);
 489             h = get_bits(gb, 8);
 490             v->s.avctx->sample_aspect_ratio = (AVRational){w, h};
 491         }
 492         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);
 493
 494         if(get_bits1(gb)){ //framerate stuff
 495             if(get_bits1(gb)) {
 496                 v->s.avctx->time_base.num = 32;
 497                 v->s.avctx->time_base.den = get_bits(gb, 16) + 1;
 498             } else {
 499                 int nr, dr;
 500                 nr = get_bits(gb, 8);
 501                 dr = get_bits(gb, 4);
 502                 if(nr && nr < 8 && dr && dr < 3){
 503                     v->s.avctx->time_base.num = ff_vc1_fps_dr[dr - 1];
 504                     v->s.avctx->time_base.den = ff_vc1_fps_nr[nr - 1] * 1000;
 505                 }
 506             }
 507         }
 508
 509         if(get_bits1(gb)){
 510             v->color_prim = get_bits(gb, 8);
 511             v->transfer_char = get_bits(gb, 8);
 512             v->matrix_coef = get_bits(gb, 8);
 513         }
 514     }
 515
 516     v->hrd_param_flag = get_bits1(gb);
 517     if(v->hrd_param_flag) {
 518         int i;
 519         v->hrd_num_leaky_buckets = get_bits(gb, 5);
 520         skip_bits(gb, 4); //bitrate exponent
 521         skip_bits(gb, 4); //buffer size exponent
 522         for(i = 0; i < v->hrd_num_leaky_buckets; i++) {
 523             skip_bits(gb, 16); //hrd_rate[n]
 524             skip_bits(gb, 16); //hrd_buffer[n]
 525         }
 526     }
 527     return 0;
 528 }
 529
 530 int vc1_decode_entry_point(AVCodecContext *avctx, VC1Context *v, GetBitContext *gb)
 531 {
 532     int i;
 533
 534     av_log(avctx, AV_LOG_DEBUG, "Entry point: %08X\n", show_bits_long(gb, 32));
 535     v->broken_link = get_bits1(gb);
 536     v->closed_entry = get_bits1(gb);
 537     v->panscanflag = get_bits1(gb);
 538     v->refdist_flag = get_bits1(gb);
 539     v->s.loop_filter = get_bits1(gb);
 540     v->fastuvmc = get_bits1(gb);
 541     v->extended_mv = get_bits1(gb);
 542     v->dquant = get_bits(gb, 2);
 543     v->vstransform = get_bits1(gb);
 544     v->overlap = get_bits1(gb);
 545     v->quantizer_mode = get_bits(gb, 2);
 546
 547     if(v->hrd_param_flag){
 548         for(i = 0; i < v->hrd_num_leaky_buckets; i++) {
 549             skip_bits(gb, 8); //hrd_full[n]
 550         }
 551     }
 552
 553     if(get_bits1(gb)){
 554         avctx->coded_width = (get_bits(gb, 12)+1)<<1;
 555         avctx->coded_height = (get_bits(gb, 12)+1)<<1;
 556     }
 557     if(v->extended_mv)
 558         v->extended_dmv = get_bits1(gb);
 559     if((v->range_mapy_flag = get_bits1(gb))) {
 560         av_log(avctx, AV_LOG_ERROR, "Luma scaling is not supported, expect wrong picture\n");
 561         v->range_mapy = get_bits(gb, 3);
 562     }
 563     if((v->range_mapuv_flag = get_bits1(gb))) {
 564         av_log(avctx, AV_LOG_ERROR, "Chroma scaling is not supported, expect wrong picture\n");
 565         v->range_mapuv = get_bits(gb, 3);
 566     }
 567
 568     av_log(avctx, AV_LOG_DEBUG, "Entry point info:\n"
 569         "BrokenLink=%i, ClosedEntry=%i, PanscanFlag=%i\n"
 570         "RefDist=%i, Postproc=%i, FastUVMC=%i, ExtMV=%i\n"
 571         "DQuant=%i, VSTransform=%i, Overlap=%i, Qmode=%i\n",
 572         v->broken_link, v->closed_entry, v->panscanflag, v->refdist_flag, v->s.loop_filter,
 573         v->fastuvmc, v->extended_mv, v->dquant, v->vstransform, v->overlap, v->quantizer_mode);
 574
 575     return 0;
 576 }
 577
 578 int vc1_parse_frame_header(VC1Context *v, GetBitContext* gb)
 579 {
 580     int pqindex, lowquant, status;
 581
 582     if(v->finterpflag) v->interpfrm = get_bits1(gb);
 583     skip_bits(gb, 2); //framecnt unused
 584     v->rangeredfrm = 0;
 585     if (v->rangered) v->rangeredfrm = get_bits1(gb);
 586     v->s.pict_type = get_bits1(gb);
 587     if (v->s.avctx->max_b_frames) {
 588         if (!v->s.pict_type) {
 589             if (get_bits1(gb)) v->s.pict_type = AV_PICTURE_TYPE_I;
 590             else v->s.pict_type = AV_PICTURE_TYPE_B;
 591         } else v->s.pict_type = AV_PICTURE_TYPE_P;
 592     } else v->s.pict_type = v->s.pict_type ? AV_PICTURE_TYPE_P : AV_PICTURE_TYPE_I;
 593
 594     v->bi_type = 0;
 595     if(v->s.pict_type == AV_PICTURE_TYPE_B) {
 596         v->bfraction_lut_index = get_vlc2(gb, ff_vc1_bfraction_vlc.table, VC1_BFRACTION_VLC_BITS, 1);
 597         v->bfraction = ff_vc1_bfraction_lut[v->bfraction_lut_index];
 598         if(v->bfraction == 0) {
 599             v->s.pict_type = AV_PICTURE_TYPE_BI;
 600         }
 601     }
 602     if(v->s.pict_type == AV_PICTURE_TYPE_I || v->s.pict_type == AV_PICTURE_TYPE_BI)
 603         skip_bits(gb, 7); // skip buffer fullness
 604
 605     if(v->parse_only)
 606         return 0;
 607
 608     /* calculate RND */
 609     if(v->s.pict_type == AV_PICTURE_TYPE_I || v->s.pict_type == AV_PICTURE_TYPE_BI)
 610         v->rnd = 1;
 611     if(v->s.pict_type == AV_PICTURE_TYPE_P)
 612         v->rnd ^= 1;
 613
 614     /* Quantizer stuff */
 615     pqindex = get_bits(gb, 5);
 616     if(!pqindex) return -1;
 617     if (v->quantizer_mode == QUANT_FRAME_IMPLICIT)
 618         v->pq = ff_vc1_pquant_table[0][pqindex];
 619     else
 620         v->pq = ff_vc1_pquant_table[1][pqindex];
 621
 622     v->pquantizer = 1;
 623     if (v->quantizer_mode == QUANT_FRAME_IMPLICIT)
 624         v->pquantizer = pqindex < 9;
 625     if (v->quantizer_mode == QUANT_NON_UNIFORM)
 626         v->pquantizer = 0;
 627     v->pqindex = pqindex;
 628     if (pqindex < 9) v->halfpq = get_bits1(gb);
 629     else v->halfpq = 0;
 630     if (v->quantizer_mode == QUANT_FRAME_EXPLICIT)
 631         v->pquantizer = get_bits1(gb);
 632     v->dquantfrm = 0;
 633     if (v->extended_mv == 1) v->mvrange = get_unary(gb, 0, 3);
 634     v->k_x = v->mvrange + 9 + (v->mvrange >> 1); //k_x can be 9 10 12 13
 635     v->k_y = v->mvrange + 8; //k_y can be 8 9 10 11
 636     v->range_x = 1 << (v->k_x - 1);
 637     v->range_y = 1 << (v->k_y - 1);
 638     if (v->multires && v->s.pict_type != AV_PICTURE_TYPE_B) v->respic = get_bits(gb, 2);
 639
 640     if(v->res_x8 && (v->s.pict_type == AV_PICTURE_TYPE_I || v->s.pict_type == AV_PICTURE_TYPE_BI)){
 641         v->x8_type = get_bits1(gb);
 642     }else v->x8_type = 0;
 643 //av_log(v->s.avctx, AV_LOG_INFO, "%c Frame: QP=[%i]%i (+%i/2) %i\n",
 644 //        (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);
 645
 646     if(v->s.pict_type == AV_PICTURE_TYPE_I || v->s.pict_type == AV_PICTURE_TYPE_P) v->use_ic = 0;
 647
 648     switch(v->s.pict_type) {
 649     case AV_PICTURE_TYPE_P:
 650         if (v->pq < 5) v->tt_index = 0;
 651         else if(v->pq < 13) v->tt_index = 1;
 652         else v->tt_index = 2;
 653
 654         lowquant = (v->pq > 12) ? 0 : 1;
 655         v->mv_mode = ff_vc1_mv_pmode_table[lowquant][get_unary(gb, 1, 4)];
 656         if (v->mv_mode == MV_PMODE_INTENSITY_COMP)
 657         {
 658             int scale, shift, i;
 659             v->mv_mode2 = ff_vc1_mv_pmode_table2[lowquant][get_unary(gb, 1, 3)];
 660             v->lumscale = get_bits(gb, 6);
 661             v->lumshift = get_bits(gb, 6);
 662             v->use_ic = 1;
 663             /* fill lookup tables for intensity compensation */
 664             if(!v->lumscale) {
 665                 scale = -64;
 666                 shift = (255 - v->lumshift * 2) << 6;
 667                 if(v->lumshift > 31)
 668                     shift += 128 << 6;
 669             } else {
 670                 scale = v->lumscale + 32;
 671                 if(v->lumshift > 31)
 672                     shift = (v->lumshift - 64) << 6;
 673                 else
 674                     shift = v->lumshift << 6;
 675             }
 676             for(i = 0; i < 256; i++) {
 677                 v->luty[i] = av_clip_uint8((scale * i + shift + 32) >> 6);
 678                 v->lutuv[i] = av_clip_uint8((scale * (i - 128) + 128*64 + 32) >> 6);
 679             }
 680         }
 681         if(v->mv_mode == MV_PMODE_1MV_HPEL || v->mv_mode == MV_PMODE_1MV_HPEL_BILIN)
 682             v->s.quarter_sample = 0;
 683         else if(v->mv_mode == MV_PMODE_INTENSITY_COMP) {
 684             if(v->mv_mode2 == MV_PMODE_1MV_HPEL || v->mv_mode2 == MV_PMODE_1MV_HPEL_BILIN)
 685                 v->s.quarter_sample = 0;
 686             else
 687                 v->s.quarter_sample = 1;
 688         } else
 689             v->s.quarter_sample = 1;
 690         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));
 691
 692         if ((v->mv_mode == MV_PMODE_INTENSITY_COMP &&
 693                  v->mv_mode2 == MV_PMODE_MIXED_MV)
 694                 || v->mv_mode == MV_PMODE_MIXED_MV)
 695         {
 696             status = bitplane_decoding(v->mv_type_mb_plane, &v->mv_type_is_raw, v);
 697             if (status < 0) return -1;
 698             av_log(v->s.avctx, AV_LOG_DEBUG, "MB MV Type plane encoding: "
 699                    "Imode: %i, Invert: %i\n", status>>1, status&1);
 700         } else {
 701             v->mv_type_is_raw = 0;
 702             memset(v->mv_type_mb_plane, 0, v->s.mb_stride * v->s.mb_height);
 703         }
 704         status = bitplane_decoding(v->s.mbskip_table, &v->skip_is_raw, v);
 705         if (status < 0) return -1;
 706         av_log(v->s.avctx, AV_LOG_DEBUG, "MB Skip plane encoding: "
 707                "Imode: %i, Invert: %i\n", status>>1, status&1);
 708
 709         /* Hopefully this is correct for P frames */
 710         v->s.mv_table_index = get_bits(gb, 2); //but using ff_vc1_ tables
 711         v->cbpcy_vlc = &ff_vc1_cbpcy_p_vlc[get_bits(gb, 2)];
 712
 713         if (v->dquant)
 714         {
 715             av_log(v->s.avctx, AV_LOG_DEBUG, "VOP DQuant info\n");
 716             vop_dquant_decoding(v);
 717         }
 718
 719         v->ttfrm = 0; //FIXME Is that so ?
 720         if (v->vstransform)
 721         {
 722             v->ttmbf = get_bits1(gb);
 723             if (v->ttmbf)
 724             {
 725                 v->ttfrm = ff_vc1_ttfrm_to_tt[get_bits(gb, 2)];
 726             }
 727         } else {
 728             v->ttmbf = 1;
 729             v->ttfrm = TT_8X8;
 730         }
 731         break;
 732     case AV_PICTURE_TYPE_B:
 733         if (v->pq < 5) v->tt_index = 0;
 734         else if(v->pq < 13) v->tt_index = 1;
 735         else v->tt_index = 2;
 736
 737         v->mv_mode = get_bits1(gb) ? MV_PMODE_1MV : MV_PMODE_1MV_HPEL_BILIN;
 738         v->s.quarter_sample = (v->mv_mode == MV_PMODE_1MV);
 739         v->s.mspel = v->s.quarter_sample;
 740
 741         status = bitplane_decoding(v->direct_mb_plane, &v->dmb_is_raw, v);
 742         if (status < 0) return -1;
 743         av_log(v->s.avctx, AV_LOG_DEBUG, "MB Direct Type plane encoding: "
 744                "Imode: %i, Invert: %i\n", status>>1, status&1);
 745         status = bitplane_decoding(v->s.mbskip_table, &v->skip_is_raw, v);
 746         if (status < 0) return -1;
 747         av_log(v->s.avctx, AV_LOG_DEBUG, "MB Skip plane encoding: "
 748                "Imode: %i, Invert: %i\n", status>>1, status&1);
 749
 750         v->s.mv_table_index = get_bits(gb, 2);
 751         v->cbpcy_vlc = &ff_vc1_cbpcy_p_vlc[get_bits(gb, 2)];
 752
 753         if (v->dquant)
 754         {
 755             av_log(v->s.avctx, AV_LOG_DEBUG, "VOP DQuant info\n");
 756             vop_dquant_decoding(v);
 757         }
 758
 759         v->ttfrm = 0;
 760         if (v->vstransform)
 761         {
 762             v->ttmbf = get_bits1(gb);
 763             if (v->ttmbf)
 764             {
 765                 v->ttfrm = ff_vc1_ttfrm_to_tt[get_bits(gb, 2)];
 766             }
 767         } else {
 768             v->ttmbf = 1;
 769             v->ttfrm = TT_8X8;
 770         }
 771         break;
 772     }
 773
 774     if(!v->x8_type)
 775     {
 776         /* AC Syntax */
 777         v->c_ac_table_index = decode012(gb);
 778         if (v->s.pict_type == AV_PICTURE_TYPE_I || v->s.pict_type == AV_PICTURE_TYPE_BI)
 779         {
 780             v->y_ac_table_index = decode012(gb);
 781         }
 782         /* DC Syntax */
 783         v->s.dc_table_index = get_bits1(gb);
 784     }
 785
 786     if(v->s.pict_type == AV_PICTURE_TYPE_BI) {
 787         v->s.pict_type = AV_PICTURE_TYPE_B;
 788         v->bi_type = 1;
 789     }
 790     return 0;
 791 }
 792
 793 int vc1_parse_frame_header_adv(VC1Context *v, GetBitContext* gb)
 794 {
 795     int pqindex, lowquant;
 796     int status;
 797
 798     v->p_frame_skipped = 0;
 799
 800     if(v->interlace){
 801         v->fcm = decode012(gb);
 802         if(v->fcm){
 803             if(!v->warn_interlaced++)
 804                 av_log(v->s.avctx, AV_LOG_ERROR, "Interlaced frames/fields support is not implemented\n");
 805             return -1;
 806         }
 807     }
 808     switch(get_unary(gb, 0, 4)) {
 809     case 0:
 810         v->s.pict_type = AV_PICTURE_TYPE_P;
 811         break;
 812     case 1:
 813         v->s.pict_type = AV_PICTURE_TYPE_B;
 814         break;
 815     case 2:
 816         v->s.pict_type = AV_PICTURE_TYPE_I;
 817         break;
 818     case 3:
 819         v->s.pict_type = AV_PICTURE_TYPE_BI;
 820         break;
 821     case 4:
 822         v->s.pict_type = AV_PICTURE_TYPE_P; // skipped pic
 823         v->p_frame_skipped = 1;
 824         return 0;
 825     }
 826     if(v->tfcntrflag)
 827         skip_bits(gb, 8);
 828     if(v->broadcast) {
 829         if(!v->interlace || v->psf) {
 830             v->rptfrm = get_bits(gb, 2);
 831         } else {
 832             v->tff = get_bits1(gb);
 833             v->rff = get_bits1(gb);
 834         }
 835     }
 836     if(v->panscanflag) {
 837         av_log_missing_feature(v->s.avctx, "Pan-scan", 0);
 838         //...
 839     }
 840     v->rnd = get_bits1(gb);
 841     if(v->interlace)
 842         v->uvsamp = get_bits1(gb);
 843     if(v->finterpflag) v->interpfrm = get_bits1(gb);
 844     if(v->s.pict_type == AV_PICTURE_TYPE_B) {
 845         v->bfraction_lut_index = get_vlc2(gb, ff_vc1_bfraction_vlc.table, VC1_BFRACTION_VLC_BITS, 1);
 846         v->bfraction = ff_vc1_bfraction_lut[v->bfraction_lut_index];
 847         if(v->bfraction == 0) {
 848             v->s.pict_type = AV_PICTURE_TYPE_BI; /* XXX: should not happen here */
 849         }
 850     }
 851     pqindex = get_bits(gb, 5);
 852     if(!pqindex) return -1;
 853     v->pqindex = pqindex;
 854     if (v->quantizer_mode == QUANT_FRAME_IMPLICIT)
 855         v->pq = ff_vc1_pquant_table[0][pqindex];
 856     else
 857         v->pq = ff_vc1_pquant_table[1][pqindex];
 858
 859     v->pquantizer = 1;
 860     if (v->quantizer_mode == QUANT_FRAME_IMPLICIT)
 861         v->pquantizer = pqindex < 9;
 862     if (v->quantizer_mode == QUANT_NON_UNIFORM)
 863         v->pquantizer = 0;
 864     v->pqindex = pqindex;
 865     if (pqindex < 9) v->halfpq = get_bits1(gb);
 866     else v->halfpq = 0;
 867     if (v->quantizer_mode == QUANT_FRAME_EXPLICIT)
 868         v->pquantizer = get_bits1(gb);
 869     if(v->postprocflag)
 870         v->postproc = get_bits(gb, 2);
 871
 872     if(v->s.pict_type == AV_PICTURE_TYPE_I || v->s.pict_type == AV_PICTURE_TYPE_P) v->use_ic = 0;
 873
 874     if(v->parse_only)
 875         return 0;
 876
 877     switch(v->s.pict_type) {
 878     case AV_PICTURE_TYPE_I:
 879     case AV_PICTURE_TYPE_BI:
 880         status = bitplane_decoding(v->acpred_plane, &v->acpred_is_raw, v);
 881         if (status < 0) return -1;
 882         av_log(v->s.avctx, AV_LOG_DEBUG, "ACPRED plane encoding: "
 883                 "Imode: %i, Invert: %i\n", status>>1, status&1);
 884         v->condover = CONDOVER_NONE;
 885         if(v->overlap && v->pq <= 8) {
 886             v->condover = decode012(gb);
 887             if(v->condover == CONDOVER_SELECT) {
 888                 status = bitplane_decoding(v->over_flags_plane, &v->overflg_is_raw, v);
 889                 if (status < 0) return -1;
 890                 av_log(v->s.avctx, AV_LOG_DEBUG, "CONDOVER plane encoding: "
 891                         "Imode: %i, Invert: %i\n", status>>1, status&1);
 892             }
 893         }
 894         break;
 895     case AV_PICTURE_TYPE_P:
 896         if (v->extended_mv) v->mvrange = get_unary(gb, 0, 3);
 897         else v->mvrange = 0;
 898         v->k_x = v->mvrange + 9 + (v->mvrange >> 1); //k_x can be 9 10 12 13
 899         v->k_y = v->mvrange + 8; //k_y can be 8 9 10 11
 900         v->range_x = 1 << (v->k_x - 1);
 901         v->range_y = 1 << (v->k_y - 1);
 902
 903         if (v->pq < 5) v->tt_index = 0;
 904         else if(v->pq < 13) v->tt_index = 1;
 905         else v->tt_index = 2;
 906
 907         lowquant = (v->pq > 12) ? 0 : 1;
 908         v->mv_mode = ff_vc1_mv_pmode_table[lowquant][get_unary(gb, 1, 4)];
 909         if (v->mv_mode == MV_PMODE_INTENSITY_COMP)
 910         {
 911             int scale, shift, i;
 912             v->mv_mode2 = ff_vc1_mv_pmode_table2[lowquant][get_unary(gb, 1, 3)];
 913             v->lumscale = get_bits(gb, 6);
 914             v->lumshift = get_bits(gb, 6);
 915             /* fill lookup tables for intensity compensation */
 916             if(!v->lumscale) {
 917                 scale = -64;
 918                 shift = (255 - v->lumshift * 2) << 6;
 919                 if(v->lumshift > 31)
 920                     shift += 128 << 6;
 921             } else {
 922                 scale = v->lumscale + 32;
 923                 if(v->lumshift > 31)
 924                     shift = (v->lumshift - 64) << 6;
 925                 else
 926                     shift = v->lumshift << 6;
 927             }
 928             for(i = 0; i < 256; i++) {
 929                 v->luty[i] = av_clip_uint8((scale * i + shift + 32) >> 6);
 930                 v->lutuv[i] = av_clip_uint8((scale * (i - 128) + 128*64 + 32) >> 6);
 931             }
 932             v->use_ic = 1;
 933         }
 934         if(v->mv_mode == MV_PMODE_1MV_HPEL || v->mv_mode == MV_PMODE_1MV_HPEL_BILIN)
 935             v->s.quarter_sample = 0;
 936         else if(v->mv_mode == MV_PMODE_INTENSITY_COMP) {
 937             if(v->mv_mode2 == MV_PMODE_1MV_HPEL || v->mv_mode2 == MV_PMODE_1MV_HPEL_BILIN)
 938                 v->s.quarter_sample = 0;
 939             else
 940                 v->s.quarter_sample = 1;
 941         } else
 942             v->s.quarter_sample = 1;
 943         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));
 944
 945         if ((v->mv_mode == MV_PMODE_INTENSITY_COMP &&
 946                  v->mv_mode2 == MV_PMODE_MIXED_MV)
 947                 || v->mv_mode == MV_PMODE_MIXED_MV)
 948         {
 949             status = bitplane_decoding(v->mv_type_mb_plane, &v->mv_type_is_raw, v);
 950             if (status < 0) return -1;
 951             av_log(v->s.avctx, AV_LOG_DEBUG, "MB MV Type plane encoding: "
 952                    "Imode: %i, Invert: %i\n", status>>1, status&1);
 953         } else {
 954             v->mv_type_is_raw = 0;
 955             memset(v->mv_type_mb_plane, 0, v->s.mb_stride * v->s.mb_height);
 956         }
 957         status = bitplane_decoding(v->s.mbskip_table, &v->skip_is_raw, v);
 958         if (status < 0) return -1;
 959         av_log(v->s.avctx, AV_LOG_DEBUG, "MB Skip plane encoding: "
 960                "Imode: %i, Invert: %i\n", status>>1, status&1);
 961
 962         /* Hopefully this is correct for P frames */
 963         v->s.mv_table_index = get_bits(gb, 2); //but using ff_vc1_ tables
 964         v->cbpcy_vlc = &ff_vc1_cbpcy_p_vlc[get_bits(gb, 2)];
 965         if (v->dquant)
 966         {
 967             av_log(v->s.avctx, AV_LOG_DEBUG, "VOP DQuant info\n");
 968             vop_dquant_decoding(v);
 969         }
 970
 971         v->ttfrm = 0; //FIXME Is that so ?
 972         if (v->vstransform)
 973         {
 974             v->ttmbf = get_bits1(gb);
 975             if (v->ttmbf)
 976             {
 977                 v->ttfrm = ff_vc1_ttfrm_to_tt[get_bits(gb, 2)];
 978             }
 979         } else {
 980             v->ttmbf = 1;
 981             v->ttfrm = TT_8X8;
 982         }
 983         break;
 984     case AV_PICTURE_TYPE_B:
 985         if (v->extended_mv) v->mvrange = get_unary(gb, 0, 3);
 986         else v->mvrange = 0;
 987         v->k_x = v->mvrange + 9 + (v->mvrange >> 1); //k_x can be 9 10 12 13
 988         v->k_y = v->mvrange + 8; //k_y can be 8 9 10 11
 989         v->range_x = 1 << (v->k_x - 1);
 990         v->range_y = 1 << (v->k_y - 1);
 991
 992         if (v->pq < 5) v->tt_index = 0;
 993         else if(v->pq < 13) v->tt_index = 1;
 994         else v->tt_index = 2;
 995
 996         v->mv_mode = get_bits1(gb) ? MV_PMODE_1MV : MV_PMODE_1MV_HPEL_BILIN;
 997         v->s.quarter_sample = (v->mv_mode == MV_PMODE_1MV);
 998         v->s.mspel = v->s.quarter_sample;
 999
1000         status = bitplane_decoding(v->direct_mb_plane, &v->dmb_is_raw, v);
1001         if (status < 0) return -1;
1002         av_log(v->s.avctx, AV_LOG_DEBUG, "MB Direct Type plane encoding: "
1003                "Imode: %i, Invert: %i\n", status>>1, status&1);
1004         status = bitplane_decoding(v->s.mbskip_table, &v->skip_is_raw, v);
1005         if (status < 0) return -1;
1006         av_log(v->s.avctx, AV_LOG_DEBUG, "MB Skip plane encoding: "
1007                "Imode: %i, Invert: %i\n", status>>1, status&1);
1008
1009         v->s.mv_table_index = get_bits(gb, 2);
1010         v->cbpcy_vlc = &ff_vc1_cbpcy_p_vlc[get_bits(gb, 2)];
1011
1012         if (v->dquant)
1013         {
1014             av_log(v->s.avctx, AV_LOG_DEBUG, "VOP DQuant info\n");
1015             vop_dquant_decoding(v);
1016         }
1017
1018         v->ttfrm = 0;
1019         if (v->vstransform)
1020         {
1021             v->ttmbf = get_bits1(gb);
1022             if (v->ttmbf)
1023             {
1024                 v->ttfrm = ff_vc1_ttfrm_to_tt[get_bits(gb, 2)];
1025             }
1026         } else {
1027             v->ttmbf = 1;
1028             v->ttfrm = TT_8X8;
1029         }
1030         break;
1031     }
1032
1033     /* AC Syntax */
1034     v->c_ac_table_index = decode012(gb);
1035     if (v->s.pict_type == AV_PICTURE_TYPE_I || v->s.pict_type == AV_PICTURE_TYPE_BI)
1036     {
1037         v->y_ac_table_index = decode012(gb);
1038     }
1039     /* DC Syntax */
1040     v->s.dc_table_index = get_bits1(gb);
1041     if ((v->s.pict_type == AV_PICTURE_TYPE_I || v->s.pict_type == AV_PICTURE_TYPE_BI) && v->dquant) {
1042         av_log(v->s.avctx, AV_LOG_DEBUG, "VOP DQuant info\n");
1043         vop_dquant_decoding(v);
1044     }
1045
1046     v->bi_type = 0;
1047     if(v->s.pict_type == AV_PICTURE_TYPE_BI) {
1048         v->s.pict_type = AV_PICTURE_TYPE_B;
1049         v->bi_type = 1;
1050     }
1051     return 0;
1052 }