git.sesse.net Git - vlc/blob - src/video_parser/vpar_blocks.c

   1 /*****************************************************************************
   2  * vpar_blocks.c : blocks parsing
   3  *****************************************************************************
   4  * Copyright (C) 1999, 2000 VideoLAN
   5  *
   6  * Authors:
   7  *
   8  * This program is free software; you can redistribute it and/or modify
   9  * it under the terms of the GNU General Public License as published by
  10  * the Free Software Foundation; either version 2 of the License, or
  11  * (at your option) any later version.
  12  *
  13  * This program 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
  16  * GNU General Public License for more details.
  17  *
  18  * You should have received a copy of the GNU General Public License
  19  * along with this program; if not, write to the Free Software
  20  * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA  02111, USA.
  21  *****************************************************************************/
  22
  23 /*****************************************************************************
  24  * Preamble
  25  *****************************************************************************/
  26 #include "defs.h"
  27
  28 #include <sys/types.h>                        /* on BSD, uio.h needs types.h */
  29 #include <sys/uio.h>                                            /* "input.h" */
  30
  31 #include "config.h"
  32 #include "common.h"
  33 #include "threads.h"
  34 #include "mtime.h"
  35 #include "plugins.h"
  36
  37 #include "intf_msg.h"
  38 #include "debug.h"                 /* XXX?? temporaire, requis par netlist.h */
  39
  40 #include "input.h"
  41 #include "input_netlist.h"
  42 #include "decoder_fifo.h"
  43 #include "video.h"
  44 #include "video_output.h"
  45
  46 #include "vdec_idct.h"
  47 #include "video_decoder.h"
  48 #include "vdec_motion.h"
  49
  50 #include "vpar_blocks.h"
  51 #include "vpar_headers.h"
  52 #include "vpar_synchro.h"
  53 #include "video_parser.h"
  54 #include "video_fifo.h"
  55
  56 static int i_count = 0;
  57
  58 /*
  59  * Welcome to vpar_blocks.c ! Here's where the heavy processor-critical parsing
  60  * task is done. This file is divided in several parts :
  61  *  - Initialization of the lookup tables
  62  *  - Decoding of coded blocks
  63  *  - Decoding of motion vectors
  64  *  - Decoding of the other macroblock structures
  65  *  - Picture data parsing management (slices and error handling)
  66  * It's a pretty long file. Good luck and have a nice day.
  67  */
  68
  69
  70 /*
  71  * Initialization tables
  72  */
  73
  74     /* Table for coded_block_pattern resolution */
  75 static lookup_t     pl_coded_pattern[512] =
  76     { {MB_ERROR, 0}, {0, 9}, {39, 9}, {27, 9}, {59, 9}, {55, 9}, {47, 9}, {31, 9},
  77     {58, 8}, {58, 8}, {54, 8}, {54, 8}, {46, 8}, {46, 8}, {30, 8}, {30, 8},
  78     {57, 8}, {57, 8}, {53, 8}, {53, 8}, {45, 8}, {45, 8}, {29, 8}, {29, 8},
  79     {38, 8}, {38, 8}, {26, 8}, {26, 8}, {37, 8}, {37, 8}, {25, 8}, {25, 8},
  80     {43, 8}, {43, 8}, {23, 8}, {23, 8}, {51, 8}, {51, 8}, {15, 8}, {15, 8},
  81     {42, 8}, {42, 8}, {22, 8}, {22, 8}, {50, 8}, {50, 8}, {14, 8}, {14, 8},
  82     {41, 8}, {41, 8}, {21, 8}, {21, 8}, {49, 8}, {49, 8}, {13, 8}, {13, 8},
  83     {35, 8}, {35, 8}, {19, 8}, {19, 8}, {11, 8}, {11, 8}, {7, 8}, {7, 8},
  84     {34, 7}, {34, 7}, {34, 7}, {34, 7}, {18, 7}, {18, 7}, {18, 7}, {18, 7},
  85     {10, 7}, {10, 7}, {10, 7}, {10, 7}, {6, 7}, {6, 7}, {6, 7}, {6, 7},
  86     {33, 7}, {33, 7}, {33, 7}, {33, 7}, {17, 7}, {17, 7}, {17, 7}, {17, 7},
  87     {9, 7}, {9, 7}, {9, 7}, {9, 7}, {5, 7}, {5, 7}, {5, 7}, {5, 7},
  88     {63, 6}, {63, 6}, {63, 6}, {63, 6}, {63, 6}, {63, 6}, {63, 6}, {63, 6},
  89     {3, 6}, {3, 6}, {3, 6}, {3, 6}, {3, 6}, {3, 6}, {3, 6}, {3, 6},
  90     {36, 6}, {36, 6}, {36, 6}, {36, 6}, {36, 6}, {36, 6}, {36, 6}, {36, 6},
  91     {24, 6}, {24, 6}, {24, 6}, {24, 6}, {24, 6}, {24, 6}, {24, 6}, {24, 6},
  92     {62, 5}, {62, 5}, {62, 5}, {62, 5}, {62, 5}, {62, 5}, {62, 5}, {62, 5},
  93     {62, 5}, {62, 5}, {62, 5}, {62, 5}, {62, 5}, {62, 5}, {62, 5}, {62, 5},
  94     {2, 5}, {2, 5}, {2, 5}, {2, 5}, {2, 5}, {2, 5}, {2, 5}, {2, 5},
  95     {2, 5}, {2, 5}, {2, 5}, {2, 5}, {2, 5}, {2, 5}, {2, 5}, {2, 5},
  96     {61, 5}, {61, 5}, {61, 5}, {61, 5}, {61, 5}, {61, 5}, {61, 5}, {61, 5},
  97     {61, 5}, {61, 5}, {61, 5}, {61, 5}, {61, 5}, {61, 5}, {61, 5}, {61, 5},
  98     {1, 5}, {1, 5}, {1, 5}, {1, 5}, {1, 5}, {1, 5}, {1, 5}, {1, 5},
  99     {1, 5}, {1, 5}, {1, 5}, {1, 5}, {1, 5}, {1, 5}, {1, 5}, {1, 5},
 100     {56, 5}, {56, 5}, {56, 5}, {56, 5}, {56, 5}, {56, 5}, {56, 5}, {56, 5},
 101     {56, 5}, {56, 5}, {56, 5}, {56, 5}, {56, 5}, {56, 5}, {56, 5}, {56, 5},
 102     {52, 5}, {52, 5}, {52, 5}, {52, 5}, {52, 5}, {52, 5}, {52, 5}, {52, 5},
 103     {52, 5}, {52, 5}, {52, 5}, {52, 5}, {52, 5}, {52, 5}, {52, 5}, {52, 5},
 104     {44, 5}, {44, 5}, {44, 5}, {44, 5}, {44, 5}, {44, 5}, {44, 5}, {44, 5},
 105     {44, 5}, {44, 5}, {44, 5}, {44, 5}, {44, 5}, {44, 5}, {44, 5}, {44, 5},
 106     {28, 5}, {28, 5}, {28, 5}, {28, 5}, {28, 5}, {28, 5}, {28, 5}, {28, 5},
 107     {28, 5}, {28, 5}, {28, 5}, {28, 5}, {28, 5}, {28, 5}, {28, 5}, {28, 5},
 108     {40, 5}, {40, 5}, {40, 5}, {40, 5}, {40, 5}, {40, 5}, {40, 5}, {40, 5},
 109     {40, 5}, {40, 5}, {40, 5}, {40, 5}, {40, 5}, {40, 5}, {40, 5}, {40, 5},
 110     {20, 5}, {20, 5}, {20, 5}, {20, 5}, {20, 5}, {20, 5}, {20, 5}, {20, 5},
 111     {20, 5}, {20, 5}, {20, 5}, {20, 5}, {20, 5}, {20, 5}, {20, 5}, {20, 5},
 112     {48, 5}, {48, 5}, {48, 5}, {48, 5}, {48, 5}, {48, 5}, {48, 5}, {48, 5},
 113     {48, 5}, {48, 5}, {48, 5}, {48, 5}, {48, 5}, {48, 5}, {48, 5}, {48, 5},
 114     {12, 5}, {12, 5}, {12, 5}, {12, 5}, {12, 5}, {12, 5}, {12, 5}, {12, 5},
 115     {12, 5}, {12, 5}, {12, 5}, {12, 5}, {12, 5}, {12, 5}, {12, 5}, {12, 5},
 116     {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4},
 117     {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4},
 118     {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4},
 119     {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4},
 120     {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4},
 121     {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4},
 122     {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4},
 123     {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4},
 124     {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4},
 125     {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4},
 126     {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4},
 127     {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4},
 128     {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4},
 129     {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4},
 130     {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4},
 131     {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4},
 132     {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3},
 133     {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3},
 134     {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3},
 135     {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3},
 136     {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3},
 137     {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3},
 138     {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3},
 139     {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3} };
 140
 141     /* Tables for dc DCT coefficients
 142      * Tables are cut in two parts to reduce memory occupation
 143      */
 144
 145     /* Table B-12, dct_dc_size_luminance, codes 00XXX ... 11110 */
 146
 147 static lookup_t     pl_dct_dc_lum_init_table_1[32] =
 148     { {1, 2}, {1, 2}, {1, 2}, {1, 2}, {1, 2}, {1, 2}, {1, 2}, {1, 2},
 149       {2, 2}, {2, 2}, {2, 2}, {2, 2}, {2, 2}, {2, 2}, {2, 2}, {2, 2},
 150       {0, 3}, {0, 3}, {0, 3}, {0, 3}, {3, 3}, {3, 3}, {3, 3}, {3, 3},
 151       {4, 3}, {4, 3}, {4, 3}, {4, 3}, {5, 4}, {5, 4}, {6, 5}, {MB_ERROR, 0}
 152     };
 153
 154 static lookup_t     ppl_dct_dc_init_table_1[2][32] =
 155 {    { {1, 2}, {1, 2}, {1, 2}, {1, 2}, {1, 2}, {1, 2}, {1, 2}, {1, 2},
 156       {2, 2}, {2, 2}, {2, 2}, {2, 2}, {2, 2}, {2, 2}, {2, 2}, {2, 2},
 157       {0, 3}, {0, 3}, {0, 3}, {0, 3}, {3, 3}, {3, 3}, {3, 3}, {3, 3},
 158       {4, 3}, {4, 3}, {4, 3}, {4, 3}, {5, 4}, {5, 4}, {6, 5}, {MB_ERROR, 0}},
 159 { {0, 2}, {0, 2}, {0, 2}, {0, 2}, {0, 2}, {0, 2}, {0, 2}, {0, 2},
 160       {1, 2}, {1, 2}, {1, 2}, {1, 2}, {1, 2}, {1, 2}, {1, 2}, {1, 2},
 161       {2, 2}, {2, 2}, {2, 2}, {2, 2}, {2, 2}, {2, 2}, {2, 2}, {2, 2},
 162       {3, 3}, {3, 3}, {3, 3}, {3, 3}, {4, 4}, {4, 4}, {5, 5}, {MB_ERROR, 0}
 163     }
 164     };
 165
 166     /* Table B-12, dct_dc_size_luminance, codes 111110xxx ... 111111111 */
 167
 168 static lookup_t     pl_dct_dc_lum_init_table_2[32] =
 169     { {7, 6}, {7, 6}, {7, 6}, {7, 6}, {7, 6}, {7, 6}, {7, 6}, {7, 6},
 170       {8, 7}, {8, 7}, {8, 7}, {8, 7}, {9, 8}, {9, 8}, {10,9}, {11,9},
 171       {MB_ERROR, 0}, {MB_ERROR, 0}, {MB_ERROR, 0}, {MB_ERROR, 0},
 172       {MB_ERROR, 0}, {MB_ERROR, 0}, {MB_ERROR, 0}, {MB_ERROR, 0},
 173       {MB_ERROR, 0}, {MB_ERROR, 0}, {MB_ERROR, 0}, {MB_ERROR, 0},
 174       {MB_ERROR, 0}, {MB_ERROR, 0}, {MB_ERROR, 0}, {MB_ERROR, 0}
 175     };
 176
 177 static lookup_t     ppl_dct_dc_init_table_2[2][32] =
 178 {    { {7, 6}, {7, 6}, {7, 6}, {7, 6}, {7, 6}, {7, 6}, {7, 6}, {7, 6},
 179       {8, 7}, {8, 7}, {8, 7}, {8, 7}, {9, 8}, {9, 8}, {10,9}, {11,9},
 180       {MB_ERROR, 0}, {MB_ERROR, 0}, {MB_ERROR, 0}, {MB_ERROR, 0},
 181       {MB_ERROR, 0}, {MB_ERROR, 0}, {MB_ERROR, 0}, {MB_ERROR, 0},
 182       {MB_ERROR, 0}, {MB_ERROR, 0}, {MB_ERROR, 0}, {MB_ERROR, 0},
 183       {MB_ERROR, 0}, {MB_ERROR, 0}, {MB_ERROR, 0}, {MB_ERROR, 0}},
 184     { {6, 6}, {6, 6}, {6, 6}, {6, 6}, {6, 6}, {6, 6}, {6, 6}, {6, 6},
 185       {6, 6}, {6, 6}, {6, 6}, {6, 6}, {6, 6}, {6, 6}, {6, 6}, {6, 6},
 186       {7, 7}, {7, 7}, {7, 7}, {7, 7}, {7, 7}, {7, 7}, {7, 7}, {7, 7},
 187       {8, 8}, {8, 8}, {8, 8}, {8, 8}, {9, 9}, {9, 9}, {10,10}, {11,10}
 188     }
 189     };
 190
 191     /* Table B-13, dct_dc_size_chrominance, codes 00xxx ... 11110 */
 192 static lookup_t     pl_dct_dc_chrom_init_table_1[32] =
 193  { {0, 2}, {0, 2}, {0, 2}, {0, 2}, {0, 2}, {0, 2}, {0, 2}, {0, 2},
 194       {1, 2}, {1, 2}, {1, 2}, {1, 2}, {1, 2}, {1, 2}, {1, 2}, {1, 2},
 195       {2, 2}, {2, 2}, {2, 2}, {2, 2}, {2, 2}, {2, 2}, {2, 2}, {2, 2},
 196       {3, 3}, {3, 3}, {3, 3}, {3, 3}, {4, 4}, {4, 4}, {5, 5}, {MB_ERROR, 0}
 197     };
 198
 199
 200    /* Table B-13, dct_dc_size_chrominance, codes 111110xxxx ... 1111111111 */
 201 static lookup_t     pl_dct_dc_chrom_init_table_2[32] =
 202     { {6, 6}, {6, 6}, {6, 6}, {6, 6}, {6, 6}, {6, 6}, {6, 6}, {6, 6},
 203       {6, 6}, {6, 6}, {6, 6}, {6, 6}, {6, 6}, {6, 6}, {6, 6}, {6, 6},
 204       {7, 7}, {7, 7}, {7, 7}, {7, 7}, {7, 7}, {7, 7}, {7, 7}, {7, 7},
 205       {8, 8}, {8, 8}, {8, 8}, {8, 8}, {9, 9}, {9, 9}, {10,10}, {11,10}
 206     };
 207
 208
 209     /* Tables for ac DCT coefficients. There are cut in many parts to save space */
 210     /* Table B-14, DCT coefficients table zero,
 211      * codes 0100 ... 1xxx (used for first (DC) coefficient)
 212      */
 213 static dct_lookup_t pl_DCT_tab_dc[12] =
 214     {
 215         {0,2,4}, {2,1,4}, {1,1,3}, {1,1,3},
 216         {0,1,1}, {0,1,1}, {0,1,1}, {0,1,1},
 217         {0,1,1}, {0,1,1}, {0,1,1}, {0,1,1}
 218     };
 219
 220     /* Table B-14, DCT coefficients table zero,
 221      * codes 0100 ... 1xxx (used for all other coefficients)
 222      */
 223 static dct_lookup_t pl_DCT_tab_ac[12] =
 224     {
 225         {0,2,4},  {2,1,4},  {1,1,3},  {1,1,3},
 226         {DCT_EOB,0,2}, {DCT_EOB,0,2}, {DCT_EOB,0,2}, {DCT_EOB,0,2}, /* EOB */
 227         {0,1,2},  {0,1,2},  {0,1,2},  {0,1,2}
 228     };
 229
 230     /* Table B-14, DCT coefficients table zero,
 231      * codes 000001xx ... 00111xxx
 232      */
 233 static dct_lookup_t pl_DCT_tab0[60] =
 234     {
 235         {DCT_ESCAPE,0,6}, {DCT_ESCAPE,0,6}, {DCT_ESCAPE,0,6}, {DCT_ESCAPE,0,6},
 236         /* Escape */
 237         {2,2,7}, {2,2,7}, {9,1,7}, {9,1,7},
 238         {0,4,7}, {0,4,7}, {8,1,7}, {8,1,7},
 239         {7,1,6}, {7,1,6}, {7,1,6}, {7,1,6},
 240         {6,1,6}, {6,1,6}, {6,1,6}, {6,1,6},
 241         {1,2,6}, {1,2,6}, {1,2,6}, {1,2,6},
 242         {5,1,6}, {5,1,6}, {5,1,6}, {5,1,6},
 243         {13,1,8}, {0,6,8}, {12,1,8}, {11,1,8},
 244         {3,2,8}, {1,3,8}, {0,5,8}, {10,1,8},
 245         {0,3,5}, {0,3,5}, {0,3,5}, {0,3,5},
 246         {0,3,5}, {0,3,5}, {0,3,5}, {0,3,5},
 247         {4,1,5}, {4,1,5}, {4,1,5}, {4,1,5},
 248         {4,1,5}, {4,1,5}, {4,1,5}, {4,1,5},
 249         {3,1,5}, {3,1,5}, {3,1,5}, {3,1,5},
 250         {3,1,5}, {3,1,5}, {3,1,5}, {3,1,5}
 251     };
 252
 253     /* Table B-15, DCT coefficients table one,
 254      * codes 000001xx ... 11111111
 255      */
 256 static dct_lookup_t pl_DCT_tab0a[252] =
 257     {
 258         {65,0,6}, {65,0,6}, {65,0,6}, {65,0,6}, /* Escape */
 259         {7,1,7}, {7,1,7}, {8,1,7}, {8,1,7},
 260         {6,1,7}, {6,1,7}, {2,2,7}, {2,2,7},
 261         {0,7,6}, {0,7,6}, {0,7,6}, {0,7,6},
 262         {0,6,6}, {0,6,6}, {0,6,6}, {0,6,6},
 263         {4,1,6}, {4,1,6}, {4,1,6}, {4,1,6},
 264         {5,1,6}, {5,1,6}, {5,1,6}, {5,1,6},
 265         {1,5,8}, {11,1,8}, {0,11,8}, {0,10,8},
 266         {13,1,8}, {12,1,8}, {3,2,8}, {1,4,8},
 267         {2,1,5}, {2,1,5}, {2,1,5}, {2,1,5},
 268         {2,1,5}, {2,1,5}, {2,1,5}, {2,1,5},
 269         {1,2,5}, {1,2,5}, {1,2,5}, {1,2,5},
 270         {1,2,5}, {1,2,5}, {1,2,5}, {1,2,5},
 271         {3,1,5}, {3,1,5}, {3,1,5}, {3,1,5},
 272         {3,1,5}, {3,1,5}, {3,1,5}, {3,1,5},
 273         {1,1,3}, {1,1,3}, {1,1,3}, {1,1,3},
 274         {1,1,3}, {1,1,3}, {1,1,3}, {1,1,3},
 275         {1,1,3}, {1,1,3}, {1,1,3}, {1,1,3},
 276         {1,1,3}, {1,1,3}, {1,1,3}, {1,1,3},
 277         {1,1,3}, {1,1,3}, {1,1,3}, {1,1,3},
 278         {1,1,3}, {1,1,3}, {1,1,3}, {1,1,3},
 279         {1,1,3}, {1,1,3}, {1,1,3}, {1,1,3},
 280         {1,1,3}, {1,1,3}, {1,1,3}, {1,1,3},
 281         {64,0,4}, {64,0,4}, {64,0,4}, {64,0,4}, /* EOB */
 282         {64,0,4}, {64,0,4}, {64,0,4}, {64,0,4},
 283         {64,0,4}, {64,0,4}, {64,0,4}, {64,0,4},
 284         {64,0,4}, {64,0,4}, {64,0,4}, {64,0,4},
 285         {0,3,4}, {0,3,4}, {0,3,4}, {0,3,4},
 286         {0,3,4}, {0,3,4}, {0,3,4}, {0,3,4},
 287         {0,3,4}, {0,3,4}, {0,3,4}, {0,3,4},
 288         {0,3,4}, {0,3,4}, {0,3,4}, {0,3,4},
 289         {0,1,2}, {0,1,2}, {0,1,2}, {0,1,2},
 290         {0,1,2}, {0,1,2}, {0,1,2}, {0,1,2},
 291         {0,1,2}, {0,1,2}, {0,1,2}, {0,1,2},
 292         {0,1,2}, {0,1,2}, {0,1,2}, {0,1,2},
 293         {0,1,2}, {0,1,2}, {0,1,2}, {0,1,2},
 294         {0,1,2}, {0,1,2}, {0,1,2}, {0,1,2},
 295         {0,1,2}, {0,1,2}, {0,1,2}, {0,1,2},
 296         {0,1,2}, {0,1,2}, {0,1,2}, {0,1,2},
 297         {0,1,2}, {0,1,2}, {0,1,2}, {0,1,2},
 298         {0,1,2}, {0,1,2}, {0,1,2}, {0,1,2},
 299         {0,1,2}, {0,1,2}, {0,1,2}, {0,1,2},
 300         {0,1,2}, {0,1,2}, {0,1,2}, {0,1,2},
 301         {0,1,2}, {0,1,2}, {0,1,2}, {0,1,2},
 302         {0,1,2}, {0,1,2}, {0,1,2}, {0,1,2},
 303         {0,1,2}, {0,1,2}, {0,1,2}, {0,1,2},
 304         {0,1,2}, {0,1,2}, {0,1,2}, {0,1,2},
 305         {0,2,3}, {0,2,3}, {0,2,3}, {0,2,3},
 306         {0,2,3}, {0,2,3}, {0,2,3}, {0,2,3},
 307         {0,2,3}, {0,2,3}, {0,2,3}, {0,2,3},
 308         {0,2,3}, {0,2,3}, {0,2,3}, {0,2,3},
 309         {0,2,3}, {0,2,3}, {0,2,3}, {0,2,3},
 310         {0,2,3}, {0,2,3}, {0,2,3}, {0,2,3},
 311         {0,2,3}, {0,2,3}, {0,2,3}, {0,2,3},
 312         {0,2,3}, {0,2,3}, {0,2,3}, {0,2,3},
 313         {0,4,5}, {0,4,5}, {0,4,5}, {0,4,5},
 314         {0,4,5}, {0,4,5}, {0,4,5}, {0,4,5},
 315         {0,5,5}, {0,5,5}, {0,5,5}, {0,5,5},
 316         {0,5,5}, {0,5,5}, {0,5,5}, {0,5,5},
 317         {9,1,7}, {9,1,7}, {1,3,7}, {1,3,7},
 318         {10,1,7}, {10,1,7}, {0,8,7}, {0,8,7},
 319         {0,9,7}, {0,9,7}, {0,12,8}, {0,13,8},
 320         {2,3,8}, {4,2,8}, {0,14,8}, {0,15,8}
 321     };
 322
 323     /* Table B-14, DCT coefficients table zero,
 324      * codes 0000001000 ... 0000001111
 325      */
 326 static dct_lookup_t pl_DCT_tab1[8] =
 327     {
 328         {16,1,10}, {5,2,10}, {0,7,10}, {2,3,10},
 329         {1,4,10}, {15,1,10}, {14,1,10}, {4,2,10}
 330     };
 331
 332     /* Table B-15, DCT coefficients table one,
 333      * codes 000000100x ... 000000111x
 334      */
 335 static dct_lookup_t pl_DCT_tab1a[8] =
 336     {
 337         {5,2,9}, {5,2,9}, {14,1,9}, {14,1,9},
 338         {2,4,10}, {16,1,10}, {15,1,9}, {15,1,9}
 339     };
 340
 341     /* Table B-14/15, DCT coefficients table zero / one,
 342      * codes 000000010000 ... 000000011111
 343      */
 344 static dct_lookup_t pl_DCT_tab2[16] =
 345     {
 346         {0,11,12}, {8,2,12}, {4,3,12}, {0,10,12},
 347         {2,4,12}, {7,2,12}, {21,1,12}, {20,1,12},
 348         {0,9,12}, {19,1,12}, {18,1,12}, {1,5,12},
 349         {3,3,12}, {0,8,12}, {6,2,12}, {17,1,12}
 350     };
 351
 352     /* Table B-14/15, DCT coefficients table zero / one,
 353      * codes 0000000010000 ... 0000000011111
 354      */
 355 static dct_lookup_t pl_DCT_tab3[16] =
 356     {
 357         {10,2,13}, {9,2,13}, {5,3,13}, {3,4,13},
 358         {2,5,13}, {1,7,13}, {1,6,13}, {0,15,13},
 359         {0,14,13}, {0,13,13}, {0,12,13}, {26,1,13},
 360         {25,1,13}, {24,1,13}, {23,1,13}, {22,1,13}
 361     };
 362
 363     /* Table B-14/15, DCT coefficients table zero / one,
 364      * codes 00000000010000 ... 00000000011111
 365      */
 366 static dct_lookup_t pl_DCT_tab4[16] =
 367     {
 368         {0,31,14}, {0,30,14}, {0,29,14}, {0,28,14},
 369         {0,27,14}, {0,26,14}, {0,25,14}, {0,24,14},
 370         {0,23,14}, {0,22,14}, {0,21,14}, {0,20,14},
 371         {0,19,14}, {0,18,14}, {0,17,14}, {0,16,14}
 372     };
 373
 374     /* Table B-14/15, DCT coefficients table zero / one,
 375      *   codes 000000000010000 ... 000000000011111
 376      */
 377 static dct_lookup_t pl_DCT_tab5[16] =
 378     {
 379     {0,40,15}, {0,39,15}, {0,38,15}, {0,37,15},
 380     {0,36,15}, {0,35,15}, {0,34,15}, {0,33,15},
 381     {0,32,15}, {1,14,15}, {1,13,15}, {1,12,15},
 382     {1,11,15}, {1,10,15}, {1,9,15}, {1,8,15}
 383     };
 384
 385     /* Table B-14/15, DCT coefficients table zero / one,
 386      * codes 0000000000010000 ... 0000000000011111
 387      */
 388 static dct_lookup_t pl_DCT_tab6[16] =
 389     {
 390         {1,18,16}, {1,17,16}, {1,16,16}, {1,15,16},
 391         {6,3,16}, {16,2,16}, {15,2,16}, {14,2,16},
 392         {13,2,16}, {12,2,16}, {11,2,16}, {31,1,16},
 393         {30,1,16}, {29,1,16}, {28,1,16}, {27,1,16}
 394     };
 395
 396     /* Pointers on tables of dct coefficients */
 397 static dct_lookup_t * ppl_dct_tab1[2] = { pl_DCT_tab_ac, pl_DCT_tab0a };
 398
 399 static dct_lookup_t * ppl_dct_tab2[2] = { pl_DCT_tab_ac, pl_DCT_tab_dc };
 400
 401
 402     /* Lookup Table for the chromatic component */
 403 static int pi_cc_index[12] = { 0, 0, 0, 0, 1, 2, 1, 2, 1, 2 };
 404
 405
 406 /*
 407  * Initialization of lookup tables
 408  */
 409
 410 /*****************************************************************************
 411  * vpar_InitCrop : Initialize the crop table for saturation
 412  *                 (ISO/IEC 13818-2 section 7.4.3)
 413  *****************************************************************************/
 414 #if defined(MPEG2_COMPLIANT) && !defined(VDEC_DFT)
 415 void vpar_InitCrop( vpar_thread_t * p_vpar )
 416 {
 417     int i_dummy;
 418
 419     p_vpar->pi_crop = p_vpar->pi_crop_buf + 4096;
 420
 421     for( i_dummy = -4096; i_dummy < -2048; i_dummy++ )
 422     {
 423         p_vpar->pi_crop[i_dummy] = -2048;
 424     }
 425     for( ; i_dummy < 2047; i_dummy++ )
 426     {
 427         p_vpar->pi_crop[i_dummy] = i_dummy;
 428     }
 429     for( ; i_dummy < 4095; i_dummy++ )
 430     {
 431         p_vpar->pi_crop[i_dummy] = 2047;
 432     }
 433 }
 434 #endif
 435
 436 /*****************************************************************************
 437  * vpar_InitMbAddrInc : Initialize the lookup table for mb_addr_inc
 438  *****************************************************************************/
 439
 440 /* Function for filling up the lookup table for mb_addr_inc */
 441 static void FillMbAddrIncTable( vpar_thread_t * p_vpar,
 442                                     int i_start, int i_end, int i_step,
 443                                     int * pi_value, int i_length )
 444 {
 445     int i_pos, i_offset;
 446     for( i_pos = i_start ; i_pos < i_end ; i_pos += i_step )
 447     {
 448         for( i_offset = 0 ; i_offset < i_step ; i_offset ++ )
 449         {
 450             p_vpar->pl_mb_addr_inc[i_pos + i_offset].i_value = * pi_value;
 451             p_vpar->pl_mb_addr_inc[i_pos + i_offset].i_length = i_length;
 452         }
 453         (*pi_value)--;
 454     }
 455 }
 456
 457 /* Function that initialize the table using the last one */
 458 void vpar_InitMbAddrInc( vpar_thread_t * p_vpar )
 459 {
 460     int i_dummy;
 461     int i_value;
 462
 463     for( i_dummy = 0 ; i_dummy < 8 ; i_dummy++ )
 464     {
 465         p_vpar->pl_mb_addr_inc[i_dummy].i_value = MB_ERROR;
 466         p_vpar->pl_mb_addr_inc[i_dummy].i_length = 0;
 467     }
 468
 469     p_vpar->pl_mb_addr_inc[8].i_value = MB_ADDRINC_ESCAPE;
 470     p_vpar->pl_mb_addr_inc[8].i_length = 11;
 471
 472     for( i_dummy = 9 ; i_dummy < 15 ; i_dummy ++ )
 473     {
 474         p_vpar->pl_mb_addr_inc[i_dummy].i_value =  MB_ERROR;
 475         p_vpar->pl_mb_addr_inc[i_dummy].i_length = 0;
 476     }
 477
 478     p_vpar->pl_mb_addr_inc[15].i_value = MB_ADDRINC_STUFFING;
 479     p_vpar->pl_mb_addr_inc[15].i_length = 11;
 480
 481     for( i_dummy = 16; i_dummy < 24; i_dummy++ )
 482     {
 483         p_vpar->pl_mb_addr_inc[i_dummy].i_value =  MB_ERROR;
 484         p_vpar->pl_mb_addr_inc[i_dummy].i_length = 0;
 485     }
 486
 487     i_value = 33;
 488
 489     FillMbAddrIncTable( p_vpar, 24, 36, 1, &i_value, 11 );
 490     FillMbAddrIncTable( p_vpar, 36, 48, 2, &i_value, 10 );
 491     FillMbAddrIncTable( p_vpar, 48, 96, 8, &i_value, 8 );
 492     FillMbAddrIncTable( p_vpar, 96, 128, 16, &i_value, 7 );
 493     FillMbAddrIncTable( p_vpar, 128, 256, 64, &i_value, 5 );
 494     FillMbAddrIncTable( p_vpar, 256, 512, 128, &i_value, 4 );
 495     FillMbAddrIncTable( p_vpar, 512, 1024, 256, &i_value, 3 );
 496     FillMbAddrIncTable( p_vpar, 1024, 2048, 1024, &i_value, 1 );
 497 }
 498
 499 /*****************************************************************************
 500  * vpar_Init*MBType : Initialize lookup table for the Macroblock type
 501  *****************************************************************************/
 502
 503 /* Fonction for filling up the tables */
 504 static void FillMBType( vpar_thread_t * p_vpar,
 505                                    int           i_mb_type,
 506                                    int           i_start,
 507                                    int           i_end,
 508                                    int           i_value,
 509                                    int           i_length )
 510 {
 511     int i_dummy;
 512
 513     for( i_dummy = i_start ; i_dummy < i_end ; i_dummy++ )
 514     {
 515         p_vpar->ppl_mb_type[i_mb_type][i_dummy].i_value = i_value;
 516         p_vpar->ppl_mb_type[i_mb_type][i_dummy].i_length = i_length;
 517     }
 518 }
 519
 520 /* Fonction that fills the table for P MB_Type */
 521 void vpar_InitPMBType( vpar_thread_t * p_vpar )
 522 {
 523     FillMBType( p_vpar, 0, 32, 64, MB_MOTION_FORWARD|MB_PATTERN, 1 );
 524     FillMBType( p_vpar, 0, 16, 32, MB_PATTERN, 2 );
 525     FillMBType( p_vpar, 0, 8, 16, MB_MOTION_FORWARD, 3 );
 526     FillMBType( p_vpar, 0, 6, 8, MB_INTRA, 5 );
 527     FillMBType( p_vpar, 0, 4, 6, MB_QUANT|MB_MOTION_FORWARD|MB_PATTERN, 5 );
 528     FillMBType( p_vpar, 0, 2, 4, MB_QUANT|MB_PATTERN, 5 );
 529     p_vpar->ppl_mb_type[0][1].i_value = MB_QUANT|MB_INTRA;
 530     p_vpar->ppl_mb_type[0][1].i_length = 6;
 531     p_vpar->ppl_mb_type[0][0].i_value = MB_ERROR;
 532     p_vpar->ppl_mb_type[0][0].i_length = 0;
 533 }
 534
 535 /* Fonction that fills the table for B MB_Type */
 536 void vpar_InitBMBType( vpar_thread_t * p_vpar )
 537 {
 538     FillMBType( p_vpar, 1, 48, 64, MB_MOTION_FORWARD
 539                                   |MB_MOTION_BACKWARD|MB_PATTERN, 2 );
 540     FillMBType( p_vpar, 1, 32, 48, MB_MOTION_FORWARD|MB_MOTION_BACKWARD, 2 );
 541     FillMBType( p_vpar, 1, 24, 32, MB_MOTION_BACKWARD|MB_PATTERN, 3 );
 542     FillMBType( p_vpar, 1, 16, 24, MB_MOTION_BACKWARD, 3 );
 543     FillMBType( p_vpar, 1, 12, 16, MB_MOTION_FORWARD|MB_PATTERN, 4 );
 544     FillMBType( p_vpar, 1, 8, 12, MB_MOTION_FORWARD, 4 );
 545     FillMBType( p_vpar, 1, 6, 8, MB_INTRA, 5 );
 546     FillMBType( p_vpar, 1, 4, 6, MB_QUANT|MB_MOTION_FORWARD
 547                                 |MB_MOTION_BACKWARD|MB_PATTERN, 5 );
 548     p_vpar->ppl_mb_type[1][3].i_value = MB_QUANT|MB_MOTION_FORWARD|MB_PATTERN;
 549     p_vpar->ppl_mb_type[1][3].i_length = 6;
 550     p_vpar->ppl_mb_type[1][2].i_value = MB_QUANT|MB_MOTION_BACKWARD|MB_PATTERN;
 551     p_vpar->ppl_mb_type[1][2].i_length = 6;
 552     p_vpar->ppl_mb_type[1][1].i_value = MB_QUANT|MB_INTRA;
 553     p_vpar->ppl_mb_type[1][1].i_length = 6;
 554     p_vpar->ppl_mb_type[1][0].i_value =MB_ERROR;
 555     p_vpar->ppl_mb_type[1][0].i_length = 0;
 556 }
 557
 558
 559 /*****************************************************************************
 560  * vpar_InitDCTTables : Initialize tables giving the length of the dct
 561  *                      coefficient from the vlc code
 562  *****************************************************************************/
 563
 564 /* First fonction for filling the table */
 565 static void FillDCTTable( dct_lookup_t * p_tab_dest, dct_lookup_t * p_tab_src,
 566                                      int i_step, int i_nb_elem, int i_offset )
 567 {
 568     int i_dummy, i_dummy2;
 569
 570     for( i_dummy=0 ; i_dummy < i_nb_elem ; i_dummy++ )
 571     {
 572         for( i_dummy2=0 ; i_dummy2 < i_step ; i_dummy2++ )
 573         {
 574             p_tab_dest[(i_dummy+i_offset)*i_step+i_dummy2] = p_tab_src[i_dummy];
 575         }
 576     }
 577 }
 578
 579
 580 /* Fonction that actually fills the table or create the pointers */
 581 void vpar_InitDCTTables( vpar_thread_t * p_vpar )
 582 {
 583     /* Tables are cut in two parts to reduce memory occupation */
 584     p_vpar->pppl_dct_dc_size[0][0] = pl_dct_dc_lum_init_table_1;
 585     p_vpar->pppl_dct_dc_size[0][1] = pl_dct_dc_lum_init_table_2;
 586     p_vpar->pppl_dct_dc_size[1][0] = pl_dct_dc_chrom_init_table_1;
 587     p_vpar->pppl_dct_dc_size[1][1] = pl_dct_dc_chrom_init_table_2;
 588
 589     /* XXX?? MB_ERROR is replaced by 0 because if we use -1 we
 590      * can block in DecodeMPEG2Intra and others */
 591     memset( p_vpar->ppl_dct_coef[0], 0, 16 );
 592     memset( p_vpar->ppl_dct_coef[1], 0, 16 );
 593
 594     /* For table B14 & B15, we have a pointer to tables */
 595     /* We fill the table thanks to the fonction defined above */
 596     FillDCTTable( p_vpar->ppl_dct_coef[0], pl_DCT_tab0, 256, 60,  4 );
 597     FillDCTTable( p_vpar->ppl_dct_coef[0], pl_DCT_tab1,  64,  8,  8 );
 598     FillDCTTable( p_vpar->ppl_dct_coef[0], pl_DCT_tab2,  16, 16, 16 );
 599     FillDCTTable( p_vpar->ppl_dct_coef[0], pl_DCT_tab3,   8, 16, 16 );
 600     FillDCTTable( p_vpar->ppl_dct_coef[0], pl_DCT_tab4,   4, 16, 16 );
 601     FillDCTTable( p_vpar->ppl_dct_coef[0], pl_DCT_tab5,   2, 16, 16 );
 602     FillDCTTable( p_vpar->ppl_dct_coef[0], pl_DCT_tab6,   1, 16, 16 );
 603
 604     FillDCTTable( p_vpar->ppl_dct_coef[1], pl_DCT_tab0a, 256, 60, 4 );
 605     FillDCTTable( p_vpar->ppl_dct_coef[1], pl_DCT_tab1a,  64,  8,  8 );
 606     FillDCTTable( p_vpar->ppl_dct_coef[1], pl_DCT_tab2,   16, 16, 16 );
 607     FillDCTTable( p_vpar->ppl_dct_coef[1], pl_DCT_tab3,    8, 16, 16 );
 608     FillDCTTable( p_vpar->ppl_dct_coef[1], pl_DCT_tab4,    4, 16, 16 );
 609     FillDCTTable( p_vpar->ppl_dct_coef[1], pl_DCT_tab5,    2, 16, 16 );
 610     FillDCTTable( p_vpar->ppl_dct_coef[1], pl_DCT_tab6,    1, 16, 16 );
 611 }
 612
 613
 614 /*
 615  * Block parsing
 616  */
 617
 618 /*****************************************************************************
 619  * DecodeMPEG1NonIntra : decode MPEG-1 non-intra blocks
 620  *****************************************************************************/
 621 static void DecodeMPEG1NonIntra( vpar_thread_t * p_vpar,
 622                                             macroblock_t * p_mb, int i_b,
 623                                             int i_chroma_format )
 624 {
 625     int         i_parse;
 626     int         i_nc;
 627     int         i_cc;
 628     int         i_coef;
 629     int         i_code;
 630     int         i_length;
 631     int         i_pos;
 632     int         i_run;
 633     int         i_level;
 634     boolean_t   b_dc;
 635     boolean_t   b_sign;
 636     boolean_t   b_chroma;
 637
 638     /* Give the chromatic component (0, 1, 2) */
 639     i_cc = pi_cc_index[i_b];
 640
 641     /* Determine whether it is luminance or not (chrominance) */
 642     b_chroma = ( i_cc + 1 ) >> 1;
 643
 644     /* There should be no D picture in non-intra blocks */
 645     if( p_vpar->picture.i_coding_type == D_CODING_TYPE )
 646         intf_ErrMsg("vpar error : D-picture in non intra block\n");
 647
 648     /* Decoding of the AC coefficients */
 649
 650     i_nc = 0;
 651     i_coef = 0;
 652     b_sign = 0;
 653
 654     for( i_parse = 0; !p_vpar->b_die; i_parse++ )
 655     {
 656         i_code = ShowBits( &p_vpar->bit_stream, 16 );
 657         if( i_code >= 16384 )
 658         {
 659             b_dc = !((i_parse + 62)/63);
 660             i_run =     ppl_dct_tab2[b_dc][(i_code>>12)-4].i_run;
 661             i_level =   ppl_dct_tab2[b_dc][(i_code>>12)-4].i_level;
 662             i_length =  ppl_dct_tab2[b_dc][(i_code>>12)-4].i_length;
 663         }
 664         else
 665         {
 666             i_run =     p_vpar->ppl_dct_coef[0][i_code].i_run;
 667             i_length =  p_vpar->ppl_dct_coef[0][i_code].i_length;
 668             i_level =   p_vpar->ppl_dct_coef[0][i_code].i_level;
 669         }
 670
 671         RemoveBits( &p_vpar->bit_stream, i_length );
 672
 673         switch( i_run )
 674         {
 675             case DCT_ESCAPE:
 676                 i_run = GetBits( &p_vpar->bit_stream, 6 );
 677                 i_level = GetBits( &p_vpar->bit_stream, 8 );
 678                 if (i_level == 0)
 679                     i_level = GetBits( &p_vpar->bit_stream, 8 );
 680                 else if (i_level == 128)
 681                     i_level = GetBits( &p_vpar->bit_stream, 8 ) - 256;
 682                 else if (i_level > 128)
 683                     i_level -= 256;
 684
 685                 if (i_level<0)
 686                     i_level = -i_level;
 687
 688                 break;
 689             case DCT_EOB:
 690 #ifdef HAVE_MMX
 691                 /* The MMX IDCT has a precision problem with non-intra
 692                  * blocks. */
 693                 p_mb->ppi_blocks[i_b][0] += 4;
 694 #endif
 695                 if( i_nc <= 1 )
 696                 {
 697                     p_mb->pf_idct[i_b] = vdec_SparseIDCT;
 698                     p_mb->pi_sparse_pos[i_b] = i_coef;
 699                 }
 700                 else
 701                 {
 702                     p_mb->pf_idct[i_b] = vdec_IDCT;
 703                 }
 704                 return;
 705
 706                 break;
 707             default:
 708                 b_sign = GetBits( &p_vpar->bit_stream, 1 );
 709         }
 710         i_coef = i_parse;
 711         i_parse += i_run;
 712         i_nc ++;
 713
 714         if( i_parse >= 64 )
 715         {
 716             break;
 717         }
 718
 719         i_pos = pi_scan[p_vpar->picture.b_alternate_scan][i_parse];
 720         i_level = ( ((i_level << 1) + 1) * p_vpar->mb.i_quantizer_scale
 721                     * p_vpar->sequence.nonintra_quant.pi_matrix[i_pos] ) >> 4;
 722         p_mb->ppi_blocks[i_b][i_pos] = b_sign ? -i_level : i_level;
 723     }
 724
 725     intf_ErrMsg("vpar error: DCT coeff (non-intra) is out of bounds\n");
 726     p_vpar->picture.b_error = 1;
 727 }
 728
 729 /*****************************************************************************
 730  * DecodeMPEG1Intra : decode MPEG-1 intra blocks
 731  *****************************************************************************/
 732 static void DecodeMPEG1Intra( vpar_thread_t * p_vpar,
 733                                          macroblock_t * p_mb, int i_b ,
 734                                          int i_chroma_format )
 735 {
 736     int             i_parse;
 737     int             i_nc;
 738     int             i_cc;
 739     int             i_coef;
 740     int             i_code;
 741     int             i_length;
 742     int             i_pos;
 743     int             i_dct_dc_size;
 744     int             i_dct_dc_diff;
 745     int             i_run;
 746     int             i_level;
 747     boolean_t       b_sign;
 748     boolean_t       b_chroma;
 749
 750     /* Give the chromatic component (0, 1, 2) */
 751     i_cc = pi_cc_index[i_b];
 752
 753     /* Determine whether it is luminance or not (chrominance) */
 754     b_chroma = ( i_cc + 1 ) >> 1;
 755
 756     /* decode length */
 757     i_code = ShowBits(&p_vpar->bit_stream, 5);
 758
 759     if (i_code<31)
 760     {
 761         i_dct_dc_size = ppl_dct_dc_init_table_1[b_chroma][i_code].i_value;
 762         i_length = ppl_dct_dc_init_table_1[b_chroma][i_code].i_length;
 763         RemoveBits( &p_vpar->bit_stream, i_length);
 764     }
 765     else
 766     {
 767         i_code = ShowBits(&p_vpar->bit_stream, (9+b_chroma)) - (0x1f0 * (b_chroma + 1));
 768         i_dct_dc_size = ppl_dct_dc_init_table_2[b_chroma][i_code].i_value;
 769         i_length = ppl_dct_dc_init_table_2[b_chroma][i_code].i_length;
 770         RemoveBits( &p_vpar->bit_stream, i_length);
 771     }
 772
 773     if (i_dct_dc_size == 0)
 774         i_dct_dc_diff = 0;
 775     else
 776     {
 777         i_dct_dc_diff = GetBits( &p_vpar->bit_stream, i_dct_dc_size);
 778         if ((i_dct_dc_diff & (1<<(i_dct_dc_size-1))) == 0)
 779             i_dct_dc_diff -= (1<<i_dct_dc_size) - 1;
 780     }
 781
 782     /* Read the actual code with the good length */
 783     p_vpar->mb.pi_dc_dct_pred[i_cc] += i_dct_dc_diff;
 784
 785     p_mb->ppi_blocks[i_b][0] = p_vpar->mb.pi_dc_dct_pred[i_cc];
 786
 787     i_nc = ( p_vpar->mb.pi_dc_dct_pred[i_cc] != 0 );
 788
 789
 790     if( p_vpar->picture.i_coding_type == D_CODING_TYPE )
 791     {
 792         /* Remove end_of_macroblock (always 1, prevents startcode emulation)
 793          * ISO/IEC 11172-2 section 2.4.2.7 and 2.4.3.6 */
 794         RemoveBits( &p_vpar->bit_stream, 1 );
 795         /* D pictures do not have AC coefficients */
 796         return;
 797     }
 798
 799
 800     /* Decoding of the AC coefficients */
 801     i_coef = 0;
 802     i_nc = 0;
 803     b_sign = 0;
 804
 805     for( i_parse = 1; !p_vpar->b_die/*i_parse < 64*/; i_parse++ )
 806     {
 807         i_code = ShowBits( &p_vpar->bit_stream, 16 );
 808         /* We use 2 main tables for the coefficients */
 809         if( i_code >= 16384 )
 810         {
 811             i_run =     ppl_dct_tab1[0][(i_code>>12)-4].i_run;
 812             i_level =   ppl_dct_tab1[0][(i_code>>12)-4].i_level;
 813             i_length =  ppl_dct_tab1[0][(i_code>>12)-4].i_length;
 814         }
 815         else
 816         {
 817             i_run =     p_vpar->ppl_dct_coef[0][i_code].i_run;
 818             i_length =  p_vpar->ppl_dct_coef[0][i_code].i_length;
 819             i_level =   p_vpar->ppl_dct_coef[0][i_code].i_level;
 820         }
 821
 822         RemoveBits( &p_vpar->bit_stream, i_length );
 823
 824         switch( i_run )
 825         {
 826             case DCT_ESCAPE:
 827                 i_run = GetBits( &p_vpar->bit_stream, 6 );
 828                 i_level = GetBits( &p_vpar->bit_stream, 8 );
 829                 if (i_level == 0)
 830                     i_level = GetBits( &p_vpar->bit_stream, 8 );
 831                 else if (i_level == 128)
 832                     i_level = GetBits( &p_vpar->bit_stream, 8 ) - 256;
 833                 else if (i_level > 128)
 834                     i_level -= 256;
 835                 break;
 836             case DCT_EOB:
 837                 if( i_nc <= 1 )
 838                 {
 839                     p_mb->pf_idct[i_b] = vdec_SparseIDCT;
 840                     p_mb->pi_sparse_pos[i_b] = i_coef;
 841                 }
 842                 else
 843                 {
 844                     p_mb->pf_idct[i_b] = vdec_IDCT;
 845                 }
 846                 return;
 847
 848                 break;
 849             default:
 850                 b_sign = GetBits( &p_vpar->bit_stream, 1 );
 851         }
 852
 853         /* Prepare the next block */
 854         i_coef = i_parse;
 855         i_parse += i_run;
 856         i_nc ++;
 857
 858         if( i_parse >= 64 )
 859         {
 860             /* We have an error in the stream */
 861             break;
 862         }
 863
 864         /* Determine the position of the block in the frame */
 865         i_pos = pi_scan[p_vpar->picture.b_alternate_scan][i_parse];
 866         i_level = ( i_level *
 867                     p_vpar->mb.i_quantizer_scale *
 868                     p_vpar->sequence.intra_quant.pi_matrix[i_pos] ) >> 3;
 869         p_mb->ppi_blocks[i_b][i_pos] = b_sign ? -i_level : i_level;
 870     }
 871
 872     intf_ErrMsg("vpar error: DCT coeff (intra) is out of bounds\n");
 873     p_vpar->picture.b_error = 1;
 874 }
 875
 876 /*****************************************************************************
 877  * DecodeMPEG2NonIntra : decode MPEG-2 non-intra blocks
 878  *****************************************************************************/
 879 static void DecodeMPEG2NonIntra( vpar_thread_t * p_vpar,
 880                                             macroblock_t * p_mb, int i_b,
 881                                             int i_chroma_format )
 882 {
 883     int         i_parse;
 884     int         i_nc;
 885     int         i_cc;
 886     int         i_coef;
 887     int         i_code;
 888     int         i_length;
 889     int         i_pos;
 890     int         i_run;
 891     int         i_level;
 892     boolean_t   b_dc;
 893     boolean_t   b_sign;
 894     boolean_t   b_chroma;
 895     int *       pi_quant;
 896
 897     /* Give the chromatic component (0, 1, 2) */
 898     i_cc = pi_cc_index[i_b];
 899
 900     /* Determine whether it is luminance or not (chrominance) */
 901     b_chroma = ( i_cc + 1 ) >> 1;
 902
 903     /* Give a pointer to the quantization matrices for intra blocks */
 904     if( (i_chroma_format == CHROMA_420) || (!b_chroma) )
 905     {
 906         pi_quant = p_vpar->sequence.nonintra_quant.pi_matrix;
 907     }
 908     else
 909     {
 910         pi_quant = p_vpar->sequence.chroma_nonintra_quant.pi_matrix;
 911     }
 912
 913     /* Decoding of the AC coefficients */
 914
 915     i_nc = 0;
 916     i_coef = 0;
 917     for( i_parse = 0; !p_vpar->b_die; i_parse++ )
 918     {
 919         i_code = ShowBits( &p_vpar->bit_stream, 16 );
 920         if( i_code >= 16384 )
 921         {
 922             b_dc = (i_parse == 0);
 923             i_run =     ppl_dct_tab2[b_dc][(i_code>>12)-4].i_run;
 924             i_level =   ppl_dct_tab2[b_dc][(i_code>>12)-4].i_level;
 925             i_length =  ppl_dct_tab2[b_dc][(i_code>>12)-4].i_length;
 926         }
 927         else
 928         {
 929             i_run =     p_vpar->ppl_dct_coef[0][i_code].i_run;
 930             i_length =  p_vpar->ppl_dct_coef[0][i_code].i_length;
 931             i_level =   p_vpar->ppl_dct_coef[0][i_code].i_level;
 932         }
 933
 934
 935         RemoveBits( &p_vpar->bit_stream, i_length );
 936
 937         switch( i_run )
 938         {
 939             case DCT_ESCAPE:
 940                 i_run = GetBits( &p_vpar->bit_stream, 6 );
 941                 i_level = GetBits( &p_vpar->bit_stream, 12 );
 942                 i_level = (b_sign = ( i_level > 2047 )) ? 4096 - i_level
 943                                                         : i_level;
 944                 break;
 945             case DCT_EOB:
 946 #ifdef HAVE_MMX
 947                 /* The MMX IDCT has a precision problem with non-intra
 948                  * blocks. */
 949                 p_mb->ppi_blocks[i_b][0] += 4;
 950 #endif
 951                 if( i_nc <= 1 )
 952                 {
 953                     p_mb->pf_idct[i_b] = vdec_SparseIDCT;
 954                     p_mb->pi_sparse_pos[i_b] = i_coef;
 955                 }
 956                 else
 957                 {
 958                     p_mb->pf_idct[i_b] = vdec_IDCT;
 959                 }
 960                 return;
 961
 962                 break;
 963             default:
 964                 b_sign = GetBits( &p_vpar->bit_stream, 1 );
 965         }
 966         i_coef = i_parse;
 967         i_parse += i_run;
 968         i_nc ++;
 969
 970         if( i_parse >= 64 )
 971         {
 972             break;
 973         }
 974
 975         i_pos = pi_scan[p_vpar->picture.b_alternate_scan][i_parse];
 976         i_level = ( ((i_level << 1) + 1) * p_vpar->mb.i_quantizer_scale
 977                     * pi_quant[i_pos] ) >> 5;
 978         p_mb->ppi_blocks[i_b][i_pos] = b_sign ? -i_level : i_level;
 979     }
 980
 981     intf_ErrMsg("vpar error: DCT coeff (non-intra) is out of bounds\n");
 982     p_vpar->picture.b_error = 1;
 983 }
 984
 985 /*****************************************************************************
 986  * DecodeMPEG2Intra : decode MPEG-2 intra blocks
 987  *****************************************************************************/
 988 static void DecodeMPEG2Intra( vpar_thread_t * p_vpar,
 989                                          macroblock_t * p_mb, int i_b,
 990                                          int i_chroma_format )
 991 {
 992     int             i_parse;
 993     int             i_nc;
 994     int             i_cc;
 995     int             i_coef;
 996     int             i_code;
 997     int             i_length;
 998     int             i_pos;
 999     int             i_dct_dc_size;
1000     int             i_dct_dc_diff;
1001     int             i_run;
1002     int             i_level;
1003     boolean_t       b_vlc_intra;
1004     boolean_t       b_sign;
1005     boolean_t       b_chroma;
1006     int *           pi_quant;
1007
1008     /* Give the chromatic component (0, 1, 2) */
1009     i_cc = pi_cc_index[i_b];
1010
1011     /* Determine whether it is luminance or not (chrominance) */
1012     b_chroma = ( i_cc + 1 ) >> 1;
1013
1014     /* Give a pointer to the quantization matrices for intra blocks */
1015     if( (i_chroma_format == CHROMA_420) || (!b_chroma) )
1016     {
1017         pi_quant = p_vpar->sequence.intra_quant.pi_matrix;
1018     }
1019     else
1020     {
1021         pi_quant = p_vpar->sequence.chroma_intra_quant.pi_matrix;
1022     }
1023
1024     /* decode length */
1025     i_code = ShowBits(&p_vpar->bit_stream, 5);
1026
1027     if (i_code<31)
1028     {
1029         i_dct_dc_size = ppl_dct_dc_init_table_1[b_chroma][i_code].i_value;
1030         i_length = ppl_dct_dc_init_table_1[b_chroma][i_code].i_length;
1031         RemoveBits( &p_vpar->bit_stream, i_length);
1032     }
1033     else
1034     {
1035         i_code = ShowBits(&p_vpar->bit_stream, (9+b_chroma)) - (0x1f0 * (b_chroma + 1));
1036         i_dct_dc_size = ppl_dct_dc_init_table_2[b_chroma][i_code].i_value;
1037         i_length = ppl_dct_dc_init_table_2[b_chroma][i_code].i_length;
1038         RemoveBits( &p_vpar->bit_stream, i_length);
1039     }
1040
1041     if (i_dct_dc_size == 0)
1042         i_dct_dc_diff = 0;
1043     else
1044     {
1045         i_dct_dc_diff = GetBits( &p_vpar->bit_stream, i_dct_dc_size);
1046         if ((i_dct_dc_diff & (1<<(i_dct_dc_size-1))) == 0)
1047             i_dct_dc_diff -= (1<<i_dct_dc_size) - 1;
1048     }
1049
1050     /* Read the actual code with the good length */
1051     p_vpar->mb.pi_dc_dct_pred[i_cc] += i_dct_dc_diff;
1052
1053     p_mb->ppi_blocks[i_b][0] = ( p_vpar->mb.pi_dc_dct_pred[i_cc] <<
1054                                ( 3 - p_vpar->picture.i_intra_dc_precision ) );
1055
1056     i_nc = ( p_vpar->mb.pi_dc_dct_pred[i_cc] != 0 );
1057
1058
1059     /* Decoding of the AC coefficients */
1060
1061     i_coef = 0;
1062     b_vlc_intra = p_vpar->picture.b_intra_vlc_format;
1063     for( i_parse = 1; !p_vpar->b_die/*i_parse < 64*/; i_parse++ )
1064     {
1065         i_code = ShowBits( &p_vpar->bit_stream, 16 );
1066         /* We use 2 main tables for the coefficients */
1067         if( i_code >= 16384 )
1068         {
1069             i_run =     ppl_dct_tab1[b_vlc_intra][(i_code>>(12-(4*b_vlc_intra)))-4].i_run;
1070             i_level =   ppl_dct_tab1[b_vlc_intra][(i_code>>(12-(4*b_vlc_intra)))-4].i_level;
1071             i_length =  ppl_dct_tab1[b_vlc_intra][(i_code>>(12-(4*b_vlc_intra)))-4].i_length;
1072         }
1073         else
1074         {
1075             i_run =     p_vpar->ppl_dct_coef[b_vlc_intra][i_code].i_run;
1076             i_length =  p_vpar->ppl_dct_coef[b_vlc_intra][i_code].i_length;
1077             i_level =   p_vpar->ppl_dct_coef[b_vlc_intra][i_code].i_level;
1078         }
1079
1080         RemoveBits( &p_vpar->bit_stream, i_length );
1081
1082         switch( i_run )
1083         {
1084             case DCT_ESCAPE:
1085                 i_run = GetBits( &p_vpar->bit_stream, 6 );
1086                 i_level = GetBits( &p_vpar->bit_stream, 12 );
1087                 i_level = (b_sign = ( i_level > 2047 )) ? 4096 - i_level
1088                                                         : i_level;
1089                 break;
1090             case DCT_EOB:
1091                 if( i_nc <= 1 )
1092                 {
1093                     p_mb->pf_idct[i_b] = vdec_SparseIDCT;
1094                     p_mb->pi_sparse_pos[i_b] = i_coef;
1095                 }
1096                 else
1097                 {
1098                     p_mb->pf_idct[i_b] = vdec_IDCT;
1099                 }
1100                 return;
1101
1102                 break;
1103             default:
1104                 b_sign = GetBits( &p_vpar->bit_stream, 1 );
1105         }
1106
1107         /* Prepare the next block */
1108         i_coef = i_parse;
1109         i_parse += i_run;
1110         i_nc ++;
1111
1112         if( i_parse >= 64 )
1113         {
1114             /* We have an error in the stream */
1115             break;
1116         }
1117
1118         /* Determine the position of the block in the frame */
1119         i_pos = pi_scan[p_vpar->picture.b_alternate_scan][i_parse];
1120         i_level = ( i_level *
1121                     p_vpar->mb.i_quantizer_scale *
1122                     pi_quant[i_pos] ) >> 4;
1123         p_mb->ppi_blocks[i_b][i_pos] = b_sign ? -i_level : i_level;
1124     }
1125
1126     intf_ErrMsg("vpar error: DCT coeff (intra) is out of bounds\n");
1127     p_vpar->picture.b_error = 1;
1128 }
1129
1130
1131 /*
1132  * Motion vectors
1133  */
1134
1135 /****************************************************************************
1136  * MotionCode : Parse the next motion code
1137  ****************************************************************************/
1138 static int MotionCode( vpar_thread_t * p_vpar )
1139 {
1140     int i_code;
1141     static lookup_t pl_mv_tab0[8] =
1142         { {-1,0}, {3,3}, {2,2}, {2,2}, {1,1}, {1,1}, {1,1}, {1,1} };
1143     /* Table B-10, motion_code, codes 0000011 ... 000011x */
1144     static lookup_t pl_mv_tab1[8] =
1145         { {-1,0}, {-1,0}, {-1,0}, {7,6}, {6,6}, {5,6}, {4,5}, {4,5} };
1146     /* Table B-10, motion_code, codes 0000001100 ... 000001011x */
1147     static lookup_t pl_mv_tab2[12] = {
1148         {16,9}, {15,9}, {14,9}, {13,9},
1149         {12,9}, {11,9}, {10,8}, {10,8},
1150         {9,8},  {9,8},  {8,8},  {8,8} };
1151
1152     if( GetBits( &p_vpar->bit_stream, 1 ) )
1153     {
1154         return 0;
1155     }
1156     if( (i_code = ShowBits( &p_vpar->bit_stream, 9) ) >= 64 )
1157     {
1158         i_code >>= 6;
1159         RemoveBits( &p_vpar->bit_stream, pl_mv_tab0[i_code].i_length );
1160         return( GetBits( &p_vpar->bit_stream, 1 ) ?
1161             -pl_mv_tab0[i_code].i_value : pl_mv_tab0[i_code].i_value );
1162     }
1163
1164     if( i_code >= 24 )
1165     {
1166         i_code >>= 3;
1167         RemoveBits( &p_vpar->bit_stream, pl_mv_tab1[i_code].i_length );
1168         return( GetBits( &p_vpar->bit_stream, 1 ) ?
1169             -pl_mv_tab1[i_code].i_value : pl_mv_tab1[i_code].i_value );
1170     }
1171
1172     if( (i_code -= 12) < 0 )
1173     {
1174         p_vpar->picture.b_error = 1;
1175         intf_DbgMsg( "vpar debug: Invalid motion_vector code\n" );
1176         return 0;
1177     }
1178
1179     RemoveBits( &p_vpar->bit_stream, pl_mv_tab2[i_code].i_length );
1180     return( GetBits( &p_vpar->bit_stream, 1 ) ?
1181         -pl_mv_tab2[i_code].i_value : pl_mv_tab2[i_code].i_value );
1182 }
1183
1184 /****************************************************************************
1185  * DecodeMotionVector : Decode a motion_vector
1186  ****************************************************************************/
1187 static void DecodeMotionVector( int * pi_prediction, int i_r_size,
1188         int i_motion_code, int i_motion_residual, int i_full_pel )
1189 {
1190     int i_limit, i_vector;
1191
1192     /* ISO/IEC 13818-1 section 7.6.3.1 */
1193     i_limit = 16 << i_r_size;
1194     i_vector = *pi_prediction >> i_full_pel;
1195
1196     if( i_motion_code > 0 )
1197     {
1198         i_vector += ((i_motion_code-1) << i_r_size) + i_motion_residual + 1;
1199         if( i_vector >= i_limit )
1200             i_vector -= i_limit + i_limit;
1201     }
1202     else if( i_motion_code < 0 )
1203     {
1204         i_vector -= ((-i_motion_code-1) << i_r_size) + i_motion_residual + 1;
1205         if( i_vector < -i_limit )
1206             i_vector += i_limit + i_limit;
1207     }
1208     *pi_prediction = i_vector << i_full_pel;
1209 }
1210
1211 /****************************************************************************
1212  * MotionVector : Parse the next motion_vector field
1213  ****************************************************************************/
1214 static void MotionVector( vpar_thread_t * p_vpar,
1215                                      macroblock_t * p_mb, int i_r,
1216                                      int i_s, int i_full_pel, int i_structure )
1217 {
1218     int i_motion_code, i_motion_residual;
1219     int i_r_size;
1220     int pi_dm_vector[2];
1221
1222     i_r_size = p_vpar->picture.ppi_f_code[i_s][0] - 1;
1223     i_motion_code = MotionCode( p_vpar );
1224     i_motion_residual = (i_r_size != 0 && i_motion_code != 0) ?
1225                         GetBits( &p_vpar->bit_stream, i_r_size) : 0;
1226     DecodeMotionVector( &p_vpar->mb.pppi_pmv[i_r][i_s][0], i_r_size,
1227                         i_motion_code, i_motion_residual, i_full_pel );
1228     p_mb->pppi_motion_vectors[i_r][i_s][0] = p_vpar->mb.pppi_pmv[i_r][i_s][0];
1229
1230     if( p_vpar->mb.b_dmv )
1231     {
1232         if( GetBits(&p_vpar->bit_stream, 1) )
1233         {
1234             pi_dm_vector[0] = GetBits( &p_vpar->bit_stream, 1 ) ? -1 : 1;
1235         }
1236         else
1237         {
1238             pi_dm_vector[0] = 0;
1239         }
1240     }
1241
1242     i_r_size = p_vpar->picture.ppi_f_code[i_s][1]-1;
1243     i_motion_code = MotionCode( p_vpar );
1244     i_motion_residual = (i_r_size != 0 && i_motion_code != 0) ?
1245                         GetBits( &p_vpar->bit_stream, i_r_size) : 0;
1246
1247
1248     if( (p_vpar->mb.i_mv_format == MOTION_FIELD)
1249         && (i_structure == FRAME_STRUCTURE) )
1250     {
1251          p_vpar->mb.pppi_pmv[i_r][i_s][1] >>= 1;
1252     }
1253
1254     DecodeMotionVector( &p_vpar->mb.pppi_pmv[i_r][i_s][1], i_r_size,
1255                         i_motion_code, i_motion_residual, i_full_pel );
1256
1257     if( (p_vpar->mb.i_mv_format == MOTION_FIELD)
1258         && (i_structure == FRAME_STRUCTURE) )
1259          p_vpar->mb.pppi_pmv[i_r][i_s][1] <<= 1;
1260
1261     p_mb->pppi_motion_vectors[i_r][i_s][1] = p_vpar->mb.pppi_pmv[i_r][i_s][1];
1262
1263     if( p_vpar->mb.b_dmv )
1264     {
1265         if( GetBits(&p_vpar->bit_stream, 1) )
1266         {
1267             pi_dm_vector[1] = GetBits( &p_vpar->bit_stream, 1 ) ? -1 : 1;
1268         }
1269         else
1270         {
1271             pi_dm_vector[1] = 0;
1272         }
1273
1274         /* Dual Prime Arithmetic (ISO/IEC 13818-2 section 7.6.3.6). */
1275
1276 #define i_mv_x  p_mb->pppi_motion_vectors[0][0][0]
1277         if( i_structure == FRAME_STRUCTURE )
1278         {
1279 #define i_mv_y  (p_mb->pppi_motion_vectors[0][0][1] << 1)
1280             if( p_vpar->picture.b_top_field_first )
1281             {
1282                 /* vector for prediction of top field from bottom field */
1283                 p_mb->ppi_dmv[0][0] = ((i_mv_x + (i_mv_x > 0)) >> 1) + pi_dm_vector[0];
1284                 p_mb->ppi_dmv[0][1] = ((i_mv_y + (i_mv_y > 0)) >> 1) + pi_dm_vector[1] - 1;
1285
1286                 /* vector for prediction of bottom field from top field */
1287                 p_mb->ppi_dmv[1][0] = ((3*i_mv_x + (i_mv_x > 0)) >> 1) + pi_dm_vector[0];
1288                 p_mb->ppi_dmv[1][1] = ((3*i_mv_y + (i_mv_y > 0)) >> 1) + pi_dm_vector[1] + 1;
1289             }
1290             else
1291             {
1292                 /* vector for prediction of top field from bottom field */
1293                 p_mb->ppi_dmv[0][0] = ((3*i_mv_x + (i_mv_x > 0)) >> 1) + pi_dm_vector[0];
1294                 p_mb->ppi_dmv[0][1] = ((3*i_mv_y + (i_mv_y > 0)) >> 1) + pi_dm_vector[1] - 1;
1295
1296                 /* vector for prediction of bottom field from top field */
1297                 p_mb->ppi_dmv[1][0] = ((i_mv_x + (i_mv_x > 0)) >> 1) + pi_dm_vector[0];
1298                 p_mb->ppi_dmv[1][1] = ((i_mv_y + (i_mv_y > 0)) >> 1) + pi_dm_vector[1] + 1;
1299             }
1300 #undef i_mv_y
1301         }
1302         else
1303         {
1304 #define i_mv_y  p_mb->pppi_motion_vectors[0][0][1]
1305             /* vector for prediction from field of opposite 'parity' */
1306             p_mb->ppi_dmv[0][0] = ((i_mv_x + (i_mv_x > 0)) >> 1) + pi_dm_vector[0];
1307             p_mb->ppi_dmv[0][1] = ((i_mv_y + (i_mv_y > 0)) >> 1) + pi_dm_vector[1];
1308
1309             /* correct for vertical field shift */
1310             if( p_vpar->picture.i_structure == TOP_FIELD )
1311                 p_mb->ppi_dmv[0][1]--;
1312             else
1313                 p_mb->ppi_dmv[0][1]++;
1314 #undef i_mv_y
1315         }
1316 #undef i_mv_x
1317     }
1318 }
1319
1320 /*****************************************************************************
1321  * DecodeMVMPEG1 : Parse the next MPEG-1 motion vectors
1322  *****************************************************************************/
1323 static void DecodeMVMPEG1( vpar_thread_t * p_vpar,
1324                             macroblock_t * p_mb, int i_s, int i_structure )
1325 {
1326     MotionVector( p_vpar, p_mb, 0, i_s,
1327                   p_vpar->picture.pb_full_pel_vector[i_s], i_structure );
1328 }
1329
1330 /*****************************************************************************
1331  * DecodeMVMPEG2 : Parse the next MPEG-2 motion_vectors field
1332  *****************************************************************************/
1333 static void DecodeMVMPEG2( vpar_thread_t * p_vpar,
1334                             macroblock_t * p_mb, int i_s, int i_structure )
1335 {
1336     if( p_vpar->mb.i_mv_count == 1 )
1337     {
1338         if( p_vpar->mb.i_mv_format == MOTION_FIELD && !p_vpar->mb.b_dmv )
1339         {
1340             p_mb->ppi_field_select[0][i_s] = p_mb->ppi_field_select[1][i_s]
1341                                             = GetBits( &p_vpar->bit_stream, 1 );
1342         }
1343         MotionVector( p_vpar, p_mb, 0, i_s, 0, i_structure );
1344         p_vpar->mb.pppi_pmv[1][i_s][0] = p_vpar->mb.pppi_pmv[0][i_s][0];
1345         p_vpar->mb.pppi_pmv[1][i_s][1] = p_vpar->mb.pppi_pmv[0][i_s][1];
1346         p_mb->pppi_motion_vectors[1][i_s][0] = p_vpar->mb.pppi_pmv[0][i_s][0];
1347         p_mb->pppi_motion_vectors[1][i_s][1] = p_vpar->mb.pppi_pmv[0][i_s][1];
1348     }
1349     else
1350     {
1351         p_mb->ppi_field_select[0][i_s] = GetBits( &p_vpar->bit_stream, 1 );
1352         MotionVector( p_vpar, p_mb, 0, i_s, 0, i_structure );
1353         p_mb->ppi_field_select[1][i_s] = GetBits( &p_vpar->bit_stream, 1 );
1354         MotionVector( p_vpar, p_mb, 1, i_s, 0, i_structure );
1355     }
1356 }
1357
1358
1359 /*
1360  * Macroblock information structures
1361  */
1362
1363 /*****************************************************************************
1364  * MacroblockAddressIncrement : Get the macroblock_address_increment field
1365  *****************************************************************************/
1366 static int MacroblockAddressIncrement( vpar_thread_t * p_vpar )
1367 {
1368     int i_addr_inc = 0;
1369     /* Index in the lookup table mb_addr_inc */
1370     int    i_index = ShowBits( &p_vpar->bit_stream, 11 );
1371
1372     /* Test the presence of the escape character */
1373     while( i_index == 8 )
1374     {
1375         RemoveBits( &p_vpar->bit_stream, 11 );
1376         i_addr_inc += 33;
1377         i_index = ShowBits( &p_vpar->bit_stream, 11 );
1378     }
1379
1380     /* Affect the value from the lookup table */
1381     i_addr_inc += p_vpar->pl_mb_addr_inc[i_index].i_value;
1382
1383     /* Dump the good number of bits */
1384     RemoveBits( &p_vpar->bit_stream, p_vpar->pl_mb_addr_inc[i_index].i_length );
1385
1386     return i_addr_inc;
1387 }
1388
1389 /*****************************************************************************
1390  * IMBType : macroblock_type in I pictures
1391  *****************************************************************************/
1392 static int IMBType( vpar_thread_t * p_vpar )
1393 {
1394     /* Take two bits for testing */
1395     int                 i_type = ShowBits( &p_vpar->bit_stream, 2 );
1396
1397     /* Lookup table for macroblock_type */
1398     static lookup_t     pl_mb_Itype[4] = { {MB_ERROR, 0},
1399                                            {MB_QUANT|MB_INTRA, 2},
1400                                            {MB_INTRA, 1},
1401                                            {MB_INTRA, 1} };
1402     /* Dump the good number of bits */
1403     RemoveBits( &p_vpar->bit_stream, pl_mb_Itype[i_type].i_length );
1404     return pl_mb_Itype[i_type].i_value;
1405 }
1406
1407 /*****************************************************************************
1408  * PMBType : macroblock_type in P pictures
1409  *****************************************************************************/
1410 static int PMBType( vpar_thread_t * p_vpar )
1411 {
1412     /* Testing on 6 bits */
1413     int                i_type = ShowBits( &p_vpar->bit_stream, 6 );
1414
1415     /* Dump the good number of bits */
1416     RemoveBits( &p_vpar->bit_stream, p_vpar->ppl_mb_type[0][i_type].i_length );
1417     /* return the value from the lookup table for P type */
1418     return p_vpar->ppl_mb_type[0][i_type].i_value;
1419 }
1420
1421 /*****************************************************************************
1422  * BMBType : macroblock_type in B pictures
1423  *****************************************************************************/
1424 static int BMBType( vpar_thread_t * p_vpar )
1425 {
1426      /* Testing on 6 bits */
1427     int                i_type = ShowBits( &p_vpar->bit_stream, 6 );
1428
1429     /* Dump the good number of bits */
1430     RemoveBits( &p_vpar->bit_stream, p_vpar->ppl_mb_type[1][i_type].i_length );
1431
1432     /* return the value from the lookup table for B type */
1433     return p_vpar->ppl_mb_type[1][i_type].i_value;
1434 }
1435
1436 /*****************************************************************************
1437  * DMBType : macroblock_type in D pictures
1438  *****************************************************************************/
1439 static int DMBType( vpar_thread_t * p_vpar )
1440 {
1441     return GetBits( &p_vpar->bit_stream, 1 );
1442 }
1443
1444 /*****************************************************************************
1445  * CodedPattern420 : coded_block_pattern with 4:2:0 chroma
1446  *****************************************************************************/
1447 static int CodedPattern420( vpar_thread_t * p_vpar )
1448 {
1449     /* Take the max 9 bits length vlc code for testing */
1450     int      i_vlc = ShowBits( &p_vpar->bit_stream, 9 );
1451
1452     /* Trash the good number of bits read in the lookup table */
1453     RemoveBits( &p_vpar->bit_stream, pl_coded_pattern[i_vlc].i_length );
1454
1455     /* return the value from the vlc table */
1456     return pl_coded_pattern[i_vlc].i_value;
1457 }
1458
1459 /*****************************************************************************
1460  * CodedPattern422 : coded_block_pattern with 4:2:2 chroma
1461  *****************************************************************************/
1462 static int CodedPattern422( vpar_thread_t * p_vpar )
1463 {
1464     int      i_vlc = ShowBits( &p_vpar->bit_stream, 9 );
1465
1466     RemoveBits( &p_vpar->bit_stream, pl_coded_pattern[i_vlc].i_length );
1467
1468     /* Supplementary 2 bits long code for 4:2:2 format */
1469     return pl_coded_pattern[i_vlc].i_value |
1470            (GetBits( &p_vpar->bit_stream, 2 ) << 6);
1471 }
1472
1473 /*****************************************************************************
1474  * CodedPattern444 : coded_block_pattern with 4:4:4 chroma
1475  *****************************************************************************/
1476 static int CodedPattern444( vpar_thread_t * p_vpar )
1477 {
1478     int      i_vlc = ShowBits( &p_vpar->bit_stream, 9 );
1479
1480     RemoveBits( &p_vpar->bit_stream, pl_coded_pattern[i_vlc].i_length );
1481
1482     return pl_coded_pattern[i_vlc].i_value |
1483            (GetBits( &p_vpar->bit_stream, 6 ) << 6);
1484 }
1485
1486 /*****************************************************************************
1487  * InitMacroblock : Initialize macroblock values
1488  *****************************************************************************/
1489 static void InitMacroblock( vpar_thread_t * p_vpar,
1490                                        macroblock_t * p_mb, int i_coding_type,
1491                                        int i_chroma_format,
1492                                        int i_structure,
1493                                        boolean_t b_second_field )
1494 {
1495     p_mb->i_chroma_nb_blocks = 1 << i_chroma_format;
1496     p_mb->p_picture = p_vpar->picture.p_picture;
1497
1498     if( i_coding_type == B_CODING_TYPE )
1499         p_mb->p_backward = p_vpar->sequence.p_backward;
1500     else
1501         p_mb->p_backward = NULL;
1502     if( (i_coding_type == P_CODING_TYPE) || (i_coding_type == B_CODING_TYPE) )
1503         p_mb->p_forward = p_vpar->sequence.p_forward;
1504     else
1505         p_mb->p_forward = NULL;
1506
1507     p_mb->i_l_x = p_vpar->mb.i_l_x;
1508     p_mb->i_c_x = p_vpar->mb.i_c_x;
1509     p_mb->i_motion_l_y = p_vpar->mb.i_l_y;
1510     p_mb->i_motion_c_y = p_vpar->mb.i_c_y;
1511     if( (p_mb->b_motion_field = (i_structure == BOTTOM_FIELD)) )
1512     {
1513         p_mb->i_motion_l_y--;
1514         p_mb->i_motion_c_y--;
1515     }
1516     p_mb->i_addb_l_stride = (p_mb->i_l_stride = p_vpar->picture.i_l_stride) - 8;
1517     p_mb->i_addb_c_stride = (p_mb->i_c_stride = p_vpar->picture.i_c_stride) - 8;
1518     p_mb->b_P_second = ( b_second_field && i_coding_type == P_CODING_TYPE );
1519 }
1520
1521 /*****************************************************************************
1522  * UpdateContext : Update the p_vpar contextual values
1523  *****************************************************************************/
1524 static void UpdateContext( vpar_thread_t * p_vpar, int i_structure )
1525 {
1526     /* Update macroblock real position. */
1527     p_vpar->mb.i_l_x += 16;
1528     p_vpar->mb.i_l_y += (p_vpar->mb.i_l_x / p_vpar->sequence.i_width)
1529                         * (2 - (i_structure == FRAME_STRUCTURE)) * 16;
1530     p_vpar->mb.i_l_x %= p_vpar->sequence.i_width;
1531
1532     p_vpar->mb.i_c_x += p_vpar->sequence.i_chroma_mb_width;
1533     p_vpar->mb.i_c_y += (p_vpar->mb.i_c_x / p_vpar->sequence.i_chroma_width)
1534                         * (2 - (i_structure == FRAME_STRUCTURE))
1535                         * p_vpar->sequence.i_chroma_mb_height;
1536     p_vpar->mb.i_c_x %= p_vpar->sequence.i_chroma_width;
1537 }
1538
1539 /*****************************************************************************
1540  * SkippedMacroblock : Generate a skipped macroblock with NULL motion vector
1541  *****************************************************************************/
1542 static void SkippedMacroblock( vpar_thread_t * p_vpar, int i_mb,
1543                                           int i_mb_base, int i_coding_type,
1544                                           int i_chroma_format,
1545                                           int i_structure,
1546                                           boolean_t b_second_field )
1547 {
1548     macroblock_t *  p_mb;
1549
1550     static f_motion_t   pf_motion_skipped[4][4] =
1551     {
1552         {NULL, NULL, NULL, NULL},
1553         {NULL, vdec_MotionFieldField420, vdec_MotionFieldField420,
1554             vdec_MotionFrameFrame420},
1555         {NULL, vdec_MotionFieldField422, vdec_MotionFieldField422,
1556             vdec_MotionFrameFrame422},
1557         {NULL, vdec_MotionFieldField444, vdec_MotionFieldField444,
1558             vdec_MotionFrameFrame444},
1559     };
1560
1561     if( i_coding_type == I_CODING_TYPE )
1562     {
1563         intf_DbgMsg("vpar error: skipped macroblock in I-picture\n");
1564         p_vpar->picture.b_error = 1;
1565         return;
1566     }
1567
1568     if( (p_mb = vpar_NewMacroblock( &p_vpar->vfifo )) == NULL )
1569     {
1570         /* b_die == 1 */
1571         return;
1572     }
1573 #ifdef VDEC_SMP
1574     p_vpar->picture.pp_mb[i_mb_base + i_mb] = p_mb;
1575 #endif
1576
1577     InitMacroblock( p_vpar, p_mb, i_coding_type, i_chroma_format,
1578                     i_structure, b_second_field );
1579
1580     /* Motion type is picture structure. */
1581     p_mb->pf_motion = pf_motion_skipped[i_chroma_format]
1582                                        [i_structure];
1583     p_mb->i_mb_type = MB_MOTION_FORWARD;
1584     p_mb->i_coded_block_pattern = 0;
1585     memset( p_mb->pppi_motion_vectors, 0, 8*sizeof(int) );
1586
1587     /* Set the field we use for motion compensation */
1588     p_mb->ppi_field_select[0][0] = p_mb->ppi_field_select[0][1]
1589                                  = ( i_structure == BOTTOM_FIELD );
1590
1591     UpdateContext( p_vpar, i_structure );
1592
1593 #ifndef VDEC_SMP
1594     /* Decode the macroblock NOW ! */
1595     vdec_DecodeMacroblock( p_vpar->pp_vdec[0], p_mb );
1596 #endif
1597 }
1598
1599 /*****************************************************************************
1600  * MacroblockModes : Get the macroblock_modes structure
1601  *****************************************************************************/
1602 static void MacroblockModes( vpar_thread_t * p_vpar,
1603                                         macroblock_t * p_mb,
1604                                         int i_chroma_format,
1605                                         int i_coding_type,
1606                                         int i_structure )
1607 {
1608     static int          ppi_mv_count[2][4] = { {0, 1, 2, 1}, {0, 2, 1, 1} };
1609     static int          ppi_mv_format[2][4] = { {0, 1, 1, 1}, {0, 1, 2, 1} };
1610
1611     /* Get macroblock_type. */
1612     switch( i_coding_type )
1613     {
1614     case P_CODING_TYPE:
1615         p_mb->i_mb_type = PMBType( p_vpar );
1616         break;
1617     case B_CODING_TYPE:
1618         p_mb->i_mb_type = BMBType( p_vpar );
1619         break;
1620     case I_CODING_TYPE:
1621         p_mb->i_mb_type = IMBType( p_vpar );
1622         break;
1623     case D_CODING_TYPE:
1624         p_mb->i_mb_type = DMBType( p_vpar );
1625     }
1626
1627     /* SCALABILITY : warning, we don't know if spatial_temporal_weight_code
1628      * has to be dropped, take care if you use scalable streams. */
1629     /* RemoveBits( &p_vpar->bit_stream, 2 ); */
1630
1631     if( (i_coding_type == P_CODING_TYPE || i_coding_type == B_CODING_TYPE)
1632         && (p_mb->i_mb_type & (MB_MOTION_FORWARD | MB_MOTION_BACKWARD)) )
1633     {
1634         if( !(i_structure == FRAME_STRUCTURE
1635                && p_vpar->picture.b_frame_pred_frame_dct) )
1636         {
1637             p_vpar->mb.i_motion_type = GetBits( &p_vpar->bit_stream, 2 );
1638         }
1639         else
1640         {
1641             p_vpar->mb.i_motion_type = MOTION_FRAME;
1642         }
1643
1644         /* XXX?? */
1645         p_vpar->mb.i_mv_count = ppi_mv_count[i_structure == FRAME_STRUCTURE]
1646                                             [p_vpar->mb.i_motion_type];
1647         p_vpar->mb.i_mv_format = ppi_mv_format[i_structure == FRAME_STRUCTURE]
1648                                               [p_vpar->mb.i_motion_type];
1649         p_vpar->mb.b_dmv = p_vpar->mb.i_motion_type == MOTION_DMV;
1650     }
1651
1652     p_vpar->mb.b_dct_type = 0;
1653     if( (i_structure == FRAME_STRUCTURE) &&
1654         (!p_vpar->picture.b_frame_pred_frame_dct) &&
1655         (p_mb->i_mb_type & (MB_PATTERN|MB_INTRA)) )
1656     {
1657         if( (p_vpar->mb.b_dct_type = GetBits( &p_vpar->bit_stream, 1 )) )
1658         {
1659             /* The DCT is coded on fields. Jump one line between each
1660              * sample. */
1661             p_mb->i_addb_l_stride <<= 1;
1662             p_mb->i_addb_l_stride += 8;
1663             /* With CHROMA_420, the DCT is necessarily frame-coded. */
1664             if( i_chroma_format != CHROMA_420 )
1665             {
1666                 p_mb->i_addb_c_stride <<= 1;
1667                 p_mb->i_addb_c_stride += 8;
1668             }
1669         }
1670     }
1671 }
1672
1673 /*****************************************************************************
1674  * ParseMacroblock : Parse the next macroblock
1675  *****************************************************************************/
1676 #define PARSEERROR                                                      \
1677 if( p_vpar->picture.b_error )                                           \
1678 {                                                                       \
1679     /* Mark this block as skipped (better than green blocks), and       \
1680      * go to the next slice. */                                         \
1681     (*pi_mb_address)--;                                                 \
1682     vpar_DestroyMacroblock( &p_vpar->vfifo, p_mb );                     \
1683     return;                                                             \
1684 }
1685
1686
1687 #define PARSEBLOCKS( MPEG1FUNC, MPEG2FUNC )                             \
1688 {                                                                       \
1689     i_mask = 1 << (3 + (1 << i_chroma_format));                         \
1690                                                                         \
1691     /* luminance */                                                     \
1692     p_data1 = p_mb->p_picture->p_y                                      \
1693         + p_mb->i_l_x + p_vpar->mb.i_l_y*(p_vpar->sequence.i_width);    \
1694                                                                         \
1695     for( i_b = 0 ; i_b < 4 ; i_b++, i_mask >>= 1 )                      \
1696     {                                                                   \
1697         if( p_mb->i_coded_block_pattern & i_mask )                      \
1698         {                                                               \
1699             memset( p_mb->ppi_blocks[i_b], 0, 64*sizeof(dctelem_t) );   \
1700             if( b_mpeg2 )                                               \
1701                 MPEG2FUNC( p_vpar, p_mb, i_b, i_chroma_format );        \
1702             else                                                        \
1703                 MPEG1FUNC( p_vpar, p_mb, i_b, i_chroma_format );        \
1704                                                                         \
1705             /* Calculate block coordinates. */                          \
1706             p_mb->p_data[i_b] = p_data1                                 \
1707                                 + pi_y[p_vpar->mb.b_dct_type][i_b]      \
1708                                 * p_vpar->sequence.i_width              \
1709                                 + pi_x[i_b];                            \
1710                                                                         \
1711             PARSEERROR                                                  \
1712         }                                                               \
1713     }                                                                   \
1714                                                                         \
1715     /* chrominance */                                                   \
1716     p_data1 = p_mb->p_picture->p_u                                      \
1717               + p_mb->i_c_x                                             \
1718               + p_vpar->mb.i_c_y                                        \
1719                 * (p_vpar->sequence.i_chroma_width);                    \
1720     p_data2 = p_mb->p_picture->p_v                                      \
1721                + p_mb->i_c_x                                            \
1722                + p_vpar->mb.i_c_y                                       \
1723                 * (p_vpar->sequence.i_chroma_width);                    \
1724                                                                         \
1725     for( i_b = 4; i_b < 4 + (1 << i_chroma_format);                     \
1726          i_b++, i_mask >>= 1 )                                          \
1727     {                                                                   \
1728         yuv_data_t *    pp_data[2] = {p_data1, p_data2};                \
1729                                                                         \
1730         if( p_mb->i_coded_block_pattern & i_mask )                      \
1731         {                                                               \
1732             memset( p_mb->ppi_blocks[i_b], 0, 64*sizeof(dctelem_t) );   \
1733             if( b_mpeg2 )                                               \
1734                 MPEG2FUNC( p_vpar, p_mb, i_b, i_chroma_format );        \
1735             else                                                        \
1736                 MPEG1FUNC( p_vpar, p_mb, i_b, i_chroma_format );        \
1737                                                                         \
1738             /* Calculate block coordinates. */                          \
1739             p_mb->p_data[i_b] = pp_data[i_b & 1]                        \
1740                                  + pi_y[p_vpar->mb.b_dct_type][i_b]     \
1741                                    * p_vpar->sequence.i_chroma_width    \
1742                                  + pi_x[i_b];                           \
1743                                                                         \
1744             PARSEERROR                                                  \
1745         }                                                               \
1746     }                                                                   \
1747 }
1748
1749 static void ParseMacroblock(
1750                            vpar_thread_t * p_vpar,
1751                            int * pi_mb_address,     /* previous address to be
1752                                                      * used for mb_addr_incr */
1753                            int i_mb_previous,          /* actual previous mb */
1754                            int i_mb_base,     /* non-zero if field structure */
1755                            /* The following parameters are explicit in
1756                             * optimized routines : */
1757                            boolean_t b_mpeg2,             /* you know what ? */
1758                            int i_coding_type,                /* I, P, B or D */
1759                            int i_chroma_format,     /* 4:2:0, 4:2:2 or 4:4:4 */
1760                            int i_structure,    /* T(OP), B(OTTOM) or F(RAME) */
1761                            boolean_t b_second_field )     /* second field of a
1762                                                            * field picture   */
1763 {
1764     static f_motion_t   pppf_motion[4][2][4] =
1765       {
1766         { {NULL, NULL, NULL, NULL},
1767           {NULL, NULL, NULL, NULL}
1768         },
1769         { {NULL, vdec_MotionFieldField420, vdec_MotionField16x8420,
1770             vdec_MotionFieldDMV420},
1771           {NULL, vdec_MotionFrameField420, vdec_MotionFrameFrame420,
1772             vdec_MotionFrameDMV420}
1773         },
1774         { {NULL, vdec_MotionFieldField422, vdec_MotionField16x8422,
1775             vdec_MotionFieldDMV422},
1776           {NULL, vdec_MotionFrameField422, vdec_MotionFrameFrame422,
1777             vdec_MotionFrameDMV422}
1778         },
1779         { {NULL, vdec_MotionFieldField444, vdec_MotionField16x8444,
1780             vdec_MotionFieldDMV444},
1781           {NULL, vdec_MotionFrameField444, vdec_MotionFrameFrame444,
1782             vdec_MotionFrameDMV444}
1783         }
1784       };
1785     static int      pi_x[12] = {0,8,0,8,0,0,0,0,8,8,8,8};
1786     static int      pi_y[2][12] = { {0,0,8,8,0,0,8,8,0,0,8,8},
1787                                     {0,0,1,1,0,0,1,1,0,0,1,1} };
1788     int             i_mb, i_b, i_mask;
1789     int i_inc;
1790     macroblock_t *  p_mb;
1791     yuv_data_t *    p_data1;
1792     yuv_data_t *    p_data2;
1793
1794 i_count++;
1795
1796     i_inc = MacroblockAddressIncrement( p_vpar );
1797     *pi_mb_address += i_inc;
1798
1799     if( i_inc < 0 )
1800     {
1801         intf_ErrMsg( "vpar error: bad address increment (%d)\n", i_inc );
1802         p_vpar->picture.b_error = 1;
1803         return;
1804     }
1805
1806     if( *pi_mb_address - i_mb_previous - 1 )
1807     {
1808         /* Skipped macroblock (ISO/IEC 13818-2 7.6.6). */
1809
1810         /* Reset DC predictors (7.2.1). */
1811         p_vpar->mb.pi_dc_dct_pred[0] = p_vpar->mb.pi_dc_dct_pred[1]
1812             = p_vpar->mb.pi_dc_dct_pred[2]
1813             = 1 << (7 + p_vpar->picture.i_intra_dc_precision);
1814
1815         if( i_coding_type == P_CODING_TYPE )
1816         {
1817             /* Reset motion vector predictors (ISO/IEC 13818-2 7.6.3.4). */
1818             memset( p_vpar->mb.pppi_pmv, 0, 8*sizeof(int) );
1819         }
1820
1821         for( i_mb = i_mb_previous + 1; i_mb < *pi_mb_address; i_mb++ )
1822         {
1823             SkippedMacroblock( p_vpar, i_mb, i_mb_base, i_coding_type,
1824                                i_chroma_format, i_structure, b_second_field );
1825         }
1826     }
1827
1828     /* Get a macroblock structure. */
1829     if( (p_mb = vpar_NewMacroblock( &p_vpar->vfifo )) == NULL )
1830     {
1831         /* b_die == 1 */
1832         return;
1833     }
1834 #ifdef VDEC_SMP
1835     p_vpar->picture.pp_mb[i_mb_base + *pi_mb_address] = p_mb;
1836 #endif
1837
1838     InitMacroblock( p_vpar, p_mb, i_coding_type, i_chroma_format,
1839                     i_structure, b_second_field );
1840
1841     /* Parse off macroblock_modes structure. */
1842     MacroblockModes( p_vpar, p_mb, i_chroma_format, i_coding_type,
1843                      i_structure );
1844
1845     if( p_mb->i_mb_type & MB_QUANT )
1846     {
1847         LoadQuantizerScale( p_vpar );
1848     }
1849
1850     if( (i_coding_type == P_CODING_TYPE || i_coding_type == B_CODING_TYPE)
1851          && (p_mb->i_mb_type & MB_MOTION_FORWARD) )
1852     {
1853         if( b_mpeg2 )
1854             DecodeMVMPEG2( p_vpar, p_mb, 0, i_structure );
1855         else
1856             DecodeMVMPEG1( p_vpar, p_mb, 0, i_structure );
1857         PARSEERROR
1858     }
1859
1860     if( (i_coding_type == B_CODING_TYPE)
1861          && (p_mb->i_mb_type & MB_MOTION_BACKWARD) )
1862     {
1863         if( b_mpeg2 )
1864             DecodeMVMPEG2( p_vpar, p_mb, 1, i_structure );
1865         else
1866             DecodeMVMPEG1( p_vpar, p_mb, 1, i_structure );
1867         PARSEERROR
1868     }
1869
1870     if( i_coding_type == P_CODING_TYPE
1871          && !(p_mb->i_mb_type & (MB_MOTION_FORWARD|MB_INTRA)) )
1872     {
1873         /* Special No-MC macroblock in P pictures (7.6.3.5). */
1874         p_mb->i_mb_type |= MB_MOTION_FORWARD;
1875         memset( p_vpar->mb.pppi_pmv, 0, 8*sizeof(int) );
1876         memset( p_mb->pppi_motion_vectors, 0, 8*sizeof(int) );
1877         p_vpar->mb.i_motion_type = 1 + (i_structure == FRAME_STRUCTURE);
1878         p_mb->ppi_field_select[0][0] = (i_structure == BOTTOM_FIELD);
1879     }
1880
1881     if( (i_coding_type != I_CODING_TYPE) && !(p_mb->i_mb_type & MB_INTRA) )
1882     {
1883         /* Reset DC predictors (7.2.1). */
1884         p_vpar->mb.pi_dc_dct_pred[0] = p_vpar->mb.pi_dc_dct_pred[1]
1885             = p_vpar->mb.pi_dc_dct_pred[2]
1886             = 1 << (7 + p_vpar->picture.i_intra_dc_precision);
1887
1888         /* Motion function pointer. */
1889         p_mb->pf_motion = pppf_motion[i_chroma_format]
1890                                      [i_structure == FRAME_STRUCTURE]
1891                                      [p_vpar->mb.i_motion_type];
1892
1893         if( p_mb->i_mb_type & MB_PATTERN )
1894         {
1895             switch( i_chroma_format )
1896             {
1897             case CHROMA_420:
1898                 p_mb->i_coded_block_pattern = CodedPattern420( p_vpar );
1899                 break;
1900             case CHROMA_422:
1901                 p_mb->i_coded_block_pattern = CodedPattern422( p_vpar );
1902                 break;
1903             case CHROMA_444:
1904                 p_mb->i_coded_block_pattern = CodedPattern444( p_vpar );
1905             }
1906         }
1907         else
1908         {
1909             p_mb->i_coded_block_pattern = 0;
1910         }
1911
1912         /*
1913          * Effectively decode blocks.
1914          */
1915         PARSEBLOCKS( DecodeMPEG1NonIntra, DecodeMPEG2NonIntra )
1916     }
1917     else
1918     {
1919         if( !p_vpar->picture.b_concealment_mv )
1920         {
1921             /* Reset MV predictors. */
1922             memset( p_vpar->mb.pppi_pmv, 0, 8*sizeof(int) );
1923         }
1924         else
1925         {
1926             if( b_mpeg2 )
1927                 DecodeMVMPEG2( p_vpar, p_mb, 0, i_structure );
1928             else
1929                 DecodeMVMPEG1( p_vpar, p_mb, 0, i_structure );
1930             RemoveBits( &p_vpar->bit_stream, 1 );
1931         }
1932
1933         if( p_mb->i_mb_type & MB_PATTERN )
1934         {
1935             switch( i_chroma_format )
1936             {
1937             case CHROMA_420:
1938                 p_mb->i_coded_block_pattern = CodedPattern420( p_vpar );
1939                 break;
1940             case CHROMA_422:
1941                 p_mb->i_coded_block_pattern = CodedPattern422( p_vpar );
1942                 break;
1943             case CHROMA_444:
1944                 p_mb->i_coded_block_pattern = CodedPattern444( p_vpar );
1945             }
1946         }
1947         else
1948         {
1949             p_mb->i_coded_block_pattern =
1950                                 (1 << (4 + (1 << i_chroma_format))) - 1;
1951         }
1952
1953         /*
1954          * Effectively decode blocks.
1955          */
1956         PARSEBLOCKS( DecodeMPEG1Intra, DecodeMPEG2Intra )
1957     }
1958
1959     if( !p_vpar->picture.b_error )
1960     {
1961         UpdateContext( p_vpar, i_structure );
1962 #ifndef VDEC_SMP
1963         /* Decode the macroblock NOW ! */
1964         vdec_DecodeMacroblock( p_vpar->pp_vdec[0], p_mb );
1965 #endif
1966     }
1967     else
1968     {
1969         /* Mark this block as skipped (better than green blocks), and go
1970          * to the next slice. */
1971         (*pi_mb_address)--;
1972         vpar_DestroyMacroblock( &p_vpar->vfifo, p_mb );
1973     }
1974 }
1975
1976 /*
1977  * Picture data parsing management
1978  */
1979
1980 /*****************************************************************************
1981  * SliceHeader : Parse the next slice structure
1982  *****************************************************************************/
1983 static void SliceHeader( vpar_thread_t * p_vpar,
1984                                     int * pi_mb_address, int i_mb_base,
1985                                     u32 i_vert_code, boolean_t b_high,
1986                                     boolean_t b_dp_scalable,
1987                                     boolean_t b_mpeg2, int i_coding_type,
1988                                     int i_chroma_format, int i_structure,
1989                                     boolean_t b_second_field )
1990 {
1991     int                     i_mb_address_save = *pi_mb_address;
1992
1993     p_vpar->picture.b_error = 0;
1994
1995     if( b_high )
1996     {
1997         /* Picture with more than 2800 lines. */
1998         i_vert_code += GetBits( &p_vpar->bit_stream, 3 ) << 7;
1999     }
2000     if( b_dp_scalable )
2001     {
2002         /* DATA_PARTITIONING scalability. */
2003         RemoveBits( &p_vpar->bit_stream, 7 ); /* priority_breakpoint */
2004     }
2005
2006     LoadQuantizerScale( p_vpar );
2007
2008     if( GetBits( &p_vpar->bit_stream, 1 ) )
2009     {
2010         /* intra_slice, slice_id */
2011         RemoveBits( &p_vpar->bit_stream, 8 );
2012         /* extra_information_slice */
2013         while( GetBits( &p_vpar->bit_stream, 1 ) )
2014         {
2015             RemoveBits( &p_vpar->bit_stream, 8 );
2016         }
2017     }
2018     *pi_mb_address = (i_vert_code - 1)*p_vpar->sequence.i_mb_width;
2019
2020     if( *pi_mb_address < i_mb_address_save )
2021     {
2022         intf_ErrMsg( "vpar error: slices do not follow, maybe a PES has been trashed\n" );
2023         p_vpar->picture.b_error = 1;
2024         return;
2025     }
2026
2027     /* Reset DC coefficients predictors (ISO/IEC 13818-2 7.2.1). */
2028     p_vpar->mb.pi_dc_dct_pred[0] = p_vpar->mb.pi_dc_dct_pred[1]
2029         = p_vpar->mb.pi_dc_dct_pred[2]
2030         = 1 << (7 + p_vpar->picture.i_intra_dc_precision);
2031
2032     /* Reset motion vector predictors (ISO/IEC 13818-2 7.6.3.4). */
2033     memset( p_vpar->mb.pppi_pmv, 0, 8*sizeof(int) );
2034
2035     do
2036     {
2037         if( i_mb_address_save >= p_vpar->sequence.i_mb_size )
2038         {
2039             p_vpar->picture.b_error = 1;
2040             return;
2041         }
2042         ParseMacroblock( p_vpar, pi_mb_address, i_mb_address_save,
2043                          i_mb_base, b_mpeg2, i_coding_type,
2044                          i_chroma_format, i_structure,
2045                          b_second_field );
2046         i_mb_address_save = *pi_mb_address;
2047         if( p_vpar->picture.b_error )
2048         {
2049             return;
2050         }
2051     }
2052     while( ShowBits( &p_vpar->bit_stream, 23 ) && !p_vpar->b_die );
2053     NextStartCode( p_vpar );
2054 }
2055
2056 /*****************************************************************************
2057  * PictureData : Parse off all macroblocks (ISO/IEC 13818-2 6.2.3.7)
2058  *****************************************************************************/
2059 static void PictureData( vpar_thread_t * p_vpar, int i_mb_base,
2060                                     boolean_t b_high, boolean_t b_dp_scalable,
2061                                     boolean_t b_mpeg2, int i_coding_type,
2062                                     int i_chroma_format, int i_structure,
2063                                     boolean_t b_second_field )
2064 {
2065     int         i_mb_address = 0;
2066     u32         i_dummy;
2067
2068     NextStartCode( p_vpar );
2069     while( ((i_coding_type != I_CODING_TYPE && i_coding_type != D_CODING_TYPE)
2070              || !p_vpar->picture.b_error)
2071            && i_mb_address+i_mb_base < p_vpar->sequence.i_mb_size
2072            && !p_vpar->b_die )
2073     {
2074         if( ((i_dummy = ShowBits( &p_vpar->bit_stream, 32 ))
2075                  < SLICE_START_CODE_MIN) ||
2076             (i_dummy > SLICE_START_CODE_MAX) )
2077         {
2078             intf_DbgMsg("vpar debug: premature end of picture\n");
2079             p_vpar->picture.b_error = 1;
2080             break;
2081         }
2082         RemoveBits32( &p_vpar->bit_stream );
2083
2084         /* Decode slice data. */
2085         SliceHeader( p_vpar, &i_mb_address, i_mb_base, i_dummy & 255,
2086                      b_high, b_dp_scalable, b_mpeg2, i_coding_type,
2087                      i_chroma_format, i_structure, b_second_field );
2088     }
2089
2090     /* Try to recover from error. If we missed less than half the
2091      * number of macroblocks of the picture, mark the missed ones
2092      * as skipped. */
2093     if( (i_coding_type == P_CODING_TYPE || i_coding_type == B_CODING_TYPE)
2094         && p_vpar->picture.b_error &&
2095         ( (i_mb_address-i_mb_base) > (p_vpar->sequence.i_mb_size >> 1)
2096            || (i_structure != FRAME_STRUCTURE
2097                && (i_mb_address-i_mb_base) > (p_vpar->sequence.i_mb_size >> 2) ) ) )
2098     {
2099         int         i_mb;
2100
2101         p_vpar->picture.b_error = 0;
2102         for( i_mb = i_mb_address + 1;
2103              i_mb < (p_vpar->sequence.i_mb_size
2104                      << (i_structure != FRAME_STRUCTURE));
2105              i_mb++ )
2106         {
2107             SkippedMacroblock( p_vpar, i_mb, i_mb_base, i_coding_type,
2108                                i_chroma_format, i_structure, b_second_field );
2109         }
2110     }
2111 }
2112
2113 /*****************************************************************************
2114  * vpar_PictureDataVWXYZ : Parse the next macroblock ; specific functions
2115  *****************************************************************************
2116  * V = MPEG2 ?
2117  * W = coding type ?
2118  * X = chroma format ?
2119  * Y = structure ?
2120  * Z = second field ?
2121  *****************************************************************************/
2122 void vpar_PictureDataGENERIC( vpar_thread_t * p_vpar, int i_mb_base )
2123 {
2124     PictureData( p_vpar, i_mb_base, (p_vpar->sequence.i_height > 2800),
2125                  (p_vpar->sequence.i_scalable_mode == SC_DP),
2126                  p_vpar->sequence.b_mpeg2, p_vpar->picture.i_coding_type,
2127                  p_vpar->sequence.i_chroma_format,
2128                  p_vpar->picture.i_structure,
2129                  (p_vpar->picture.i_structure !=
2130                     p_vpar->picture.i_current_structure) );
2131 }
2132
2133 #if (VPAR_OPTIM_LEVEL > 0)
2134 /* Optimizations for frame pictures */
2135 void vpar_PictureData2I420F0( vpar_thread_t * p_vpar, int i_mb_base )
2136 {
2137     PictureData( p_vpar, 0, 0, 0, 1, I_CODING_TYPE, CHROMA_420,
2138                  FRAME_STRUCTURE, 0 );
2139 }
2140
2141 void vpar_PictureData2P420F0( vpar_thread_t * p_vpar, int i_mb_base )
2142 {
2143     PictureData( p_vpar, 0, 0, 0, 1, P_CODING_TYPE, CHROMA_420,
2144                  FRAME_STRUCTURE, 0 );
2145 }
2146
2147 void vpar_PictureData2B420F0( vpar_thread_t * p_vpar, int i_mb_base )
2148 {
2149     PictureData( p_vpar, 0, 0, 0, 1, B_CODING_TYPE, CHROMA_420,
2150                  FRAME_STRUCTURE, 0 );
2151 }
2152 #endif
2153
2154 #if (VPAR_OPTIM_LEVEL > 1)
2155 /* Optimizations for field pictures */
2156 void vpar_PictureData2I420TZ( vpar_thread_t * p_vpar, int i_mb_base )
2157 {
2158     PictureData( p_vpar, i_mb_base, 0, 0, 1, I_CODING_TYPE, CHROMA_420,
2159                  TOP_FIELD, (p_vpar->picture.i_structure !=
2160                                 p_vpar->picture.i_current_structure) );
2161 }
2162
2163 void vpar_PictureData2P420TZ( vpar_thread_t * p_vpar, int i_mb_base )
2164 {
2165     PictureData( p_vpar, i_mb_base, 0, 0, 1, P_CODING_TYPE, CHROMA_420,
2166                  TOP_FIELD, (p_vpar->picture.i_structure !=
2167                                 p_vpar->picture.i_current_structure) );
2168 }
2169
2170 void vpar_PictureData2B420TZ( vpar_thread_t * p_vpar, int i_mb_base )
2171 {
2172     PictureData( p_vpar, i_mb_base, 0, 0, 1, B_CODING_TYPE, CHROMA_420,
2173                  TOP_FIELD, (p_vpar->picture.i_structure !=
2174                                 p_vpar->picture.i_current_structure) );
2175 }
2176
2177 void vpar_PictureData2I420BZ( vpar_thread_t * p_vpar, int i_mb_base )
2178 {
2179     PictureData( p_vpar, i_mb_base, 0, 0, 1, I_CODING_TYPE, CHROMA_420,
2180                  BOTTOM_FIELD, (p_vpar->picture.i_structure !=
2181                                 p_vpar->picture.i_current_structure) );
2182 }
2183
2184 void vpar_PictureData2P420BZ( vpar_thread_t * p_vpar, int i_mb_base )
2185 {
2186     PictureData( p_vpar, i_mb_base, 0, 0, 1, P_CODING_TYPE, CHROMA_420,
2187                  BOTTOM_FIELD, (p_vpar->picture.i_structure !=
2188                                 p_vpar->picture.i_current_structure) );
2189 }
2190
2191 void vpar_PictureData2B420BZ( vpar_thread_t * p_vpar, int i_mb_base )
2192 {
2193     PictureData( p_vpar, i_mb_base, 0, 0, 1, B_CODING_TYPE, CHROMA_420,
2194                  BOTTOM_FIELD, (p_vpar->picture.i_structure !=
2195                                 p_vpar->picture.i_current_structure) );
2196 }
2197 #endif