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

   1 /*****************************************************************************
   2  * vpar_blocks.c : blocks parsing
   3  * (c)1999 VideoLAN
   4  *****************************************************************************/
   5
   6 /*****************************************************************************
   7  * Preamble
   8  *****************************************************************************/
   9 #include <errno.h>
  10 #include <stdlib.h>
  11 #include <stdio.h>
  12 #include <unistd.h>
  13 #include <string.h>
  14 #include <sys/uio.h>
  15 #include <X11/Xlib.h>
  16 #include <X11/extensions/XShm.h>
  17
  18 #include "config.h"
  19 #include "common.h"
  20 #include "mtime.h"
  21 #include "vlc_thread.h"
  22
  23 #include "intf_msg.h"
  24 #include "debug.h"                    /* ?? temporaire, requis par netlist.h */
  25
  26 #include "input.h"
  27 #include "input_netlist.h"
  28 #include "decoder_fifo.h"
  29 #include "video.h"
  30 #include "video_output.h"
  31
  32 #include "vdec_idct.h"
  33 #include "video_decoder.h"
  34 #include "vdec_motion.h"
  35
  36 #include "vpar_blocks.h"
  37 #include "vpar_headers.h"
  38 #include "video_fifo.h"
  39 #include "vpar_synchro.h"
  40 #include "video_parser.h"
  41
  42 /*
  43  * Local prototypes
  44  */
  45 static __inline__ void InitMacroblock( vpar_thread_t * p_vpar,
  46                                        macroblock_t * p_mb );
  47 static __inline__ int MacroblockAddressIncrement( vpar_thread_t * p_vpar );
  48 static __inline__ void MacroblockModes( vpar_thread_t * p_vpar,
  49                                         macroblock_t * p_mb );
  50 typedef void (*f_decode_block_t)( vpar_thread_t *, macroblock_t *, int );
  51 static void vpar_DecodeMPEG1Non( vpar_thread_t * p_vpar, macroblock_t * p_mb, int i_b );
  52 static void vpar_DecodeMPEG1Intra( vpar_thread_t * p_vpar, macroblock_t * p_mb, int i_b );
  53 static void vpar_DecodeMPEG2Non( vpar_thread_t * p_vpar, macroblock_t * p_mb, int i_b );
  54 static void vpar_DecodeMPEG2Intra( vpar_thread_t * p_vpar, macroblock_t * p_mb, int i_b );
  55
  56 /*
  57  * Initialisation tables
  58  */
  59     /* Table for coded_block_pattern resolution */
  60 lookup_t     pl_coded_pattern_init_table[512] =
  61     { {MB_ERROR, 0}, {0, 9}, {39, 9}, {27, 9}, {59, 9}, {55, 9}, {47, 9}, {31, 9},
  62     {58, 8}, {58, 8}, {54, 8}, {54, 8}, {46, 8}, {46, 8}, {30, 8}, {30, 8},
  63     {57, 8}, {57, 8}, {53, 8}, {53, 8}, {45, 8}, {45, 8}, {29, 8}, {29, 8},
  64     {38, 8}, {38, 8}, {26, 8}, {26, 8}, {37, 8}, {37, 8}, {25, 8}, {25, 8},
  65     {43, 8}, {43, 8}, {23, 8}, {23, 8}, {51, 8}, {51, 8}, {15, 8}, {15, 8},
  66     {42, 8}, {42, 8}, {22, 8}, {22, 8}, {50, 8}, {50, 8}, {14, 8}, {14, 8},
  67     {41, 8}, {41, 8}, {21, 8}, {21, 8}, {49, 8}, {49, 8}, {13, 8}, {13, 8},
  68     {35, 8}, {35, 8}, {19, 8}, {19, 8}, {11, 8}, {11, 8}, {7, 8}, {7, 8},
  69     {34, 7}, {34, 7}, {34, 7}, {34, 7}, {18, 7}, {18, 7}, {18, 7}, {18, 7},
  70     {10, 7}, {10, 7}, {10, 7}, {10, 7}, {6, 7}, {6, 7}, {6, 7}, {6, 7},
  71     {33, 7}, {33, 7}, {33, 7}, {33, 7}, {17, 7}, {17, 7}, {17, 7}, {17, 7},
  72     {9, 7}, {9, 7}, {9, 7}, {9, 7}, {5, 7}, {5, 7}, {5, 7}, {5, 7},
  73     {63, 6}, {63, 6}, {63, 6}, {63, 6}, {63, 6}, {63, 6}, {63, 6}, {63, 6},
  74     {3, 6}, {3, 6}, {3, 6}, {3, 6}, {3, 6}, {3, 6}, {3, 6}, {3, 6},
  75     {36, 6}, {36, 6}, {36, 6}, {36, 6}, {36, 6}, {36, 6}, {36, 6}, {36, 6},
  76     {24, 6}, {24, 6}, {24, 6}, {24, 6}, {24, 6}, {24, 6}, {24, 6}, {24, 6},
  77     {62, 5}, {62, 5}, {62, 5}, {62, 5}, {62, 5}, {62, 5}, {62, 5}, {62, 5},
  78     {62, 5}, {62, 5}, {62, 5}, {62, 5}, {62, 5}, {62, 5}, {62, 5}, {62, 5},
  79     {2, 5}, {2, 5}, {2, 5}, {2, 5}, {2, 5}, {2, 5}, {2, 5}, {2, 5},
  80     {2, 5}, {2, 5}, {2, 5}, {2, 5}, {2, 5}, {2, 5}, {2, 5}, {2, 5},
  81     {61, 5}, {61, 5}, {61, 5}, {61, 5}, {61, 5}, {61, 5}, {61, 5}, {61, 5},
  82     {61, 5}, {61, 5}, {61, 5}, {61, 5}, {61, 5}, {61, 5}, {61, 5}, {61, 5},
  83     {1, 5}, {1, 5}, {1, 5}, {1, 5}, {1, 5}, {1, 5}, {1, 5}, {1, 5},
  84     {1, 5}, {1, 5}, {1, 5}, {1, 5}, {1, 5}, {1, 5}, {1, 5}, {1, 5},
  85     {56, 5}, {56, 5}, {56, 5}, {56, 5}, {56, 5}, {56, 5}, {56, 5}, {56, 5},
  86     {56, 5}, {56, 5}, {56, 5}, {56, 5}, {56, 5}, {56, 5}, {56, 5}, {56, 5},
  87     {52, 5}, {52, 5}, {52, 5}, {52, 5}, {52, 5}, {52, 5}, {52, 5}, {52, 5},
  88     {52, 5}, {52, 5}, {52, 5}, {52, 5}, {52, 5}, {52, 5}, {52, 5}, {52, 5},
  89     {44, 5}, {44, 5}, {44, 5}, {44, 5}, {44, 5}, {44, 5}, {44, 5}, {44, 5},
  90     {44, 5}, {44, 5}, {44, 5}, {44, 5}, {44, 5}, {44, 5}, {44, 5}, {44, 5},
  91     {28, 5}, {28, 5}, {28, 5}, {28, 5}, {28, 5}, {28, 5}, {28, 5}, {28, 5},
  92     {28, 5}, {28, 5}, {28, 5}, {28, 5}, {28, 5}, {28, 5}, {28, 5}, {28, 5},
  93     {40, 5}, {40, 5}, {40, 5}, {40, 5}, {40, 5}, {40, 5}, {40, 5}, {40, 5},
  94     {40, 5}, {40, 5}, {40, 5}, {40, 5}, {40, 5}, {40, 5}, {40, 5}, {40, 5},
  95     {20, 5}, {20, 5}, {20, 5}, {20, 5}, {20, 5}, {20, 5}, {20, 5}, {20, 5},
  96     {20, 5}, {20, 5}, {20, 5}, {20, 5}, {20, 5}, {20, 5}, {20, 5}, {20, 5},
  97     {48, 5}, {48, 5}, {48, 5}, {48, 5}, {48, 5}, {48, 5}, {48, 5}, {48, 5},
  98     {48, 5}, {48, 5}, {48, 5}, {48, 5}, {48, 5}, {48, 5}, {48, 5}, {48, 5},
  99     {12, 5}, {12, 5}, {12, 5}, {12, 5}, {12, 5}, {12, 5}, {12, 5}, {12, 5},
 100     {12, 5}, {12, 5}, {12, 5}, {12, 5}, {12, 5}, {12, 5}, {12, 5}, {12, 5},
 101     {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4},
 102     {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4},
 103     {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4},
 104     {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4},
 105     {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4},
 106     {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4},
 107     {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4},
 108     {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4},
 109     {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4},
 110     {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4},
 111     {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4},
 112     {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4},
 113     {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4},
 114     {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4},
 115     {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4},
 116     {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4},
 117     {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3},
 118     {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3},
 119     {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3},
 120     {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3},
 121     {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3},
 122     {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3},
 123     {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3},
 124     {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3} };
 125
 126     /* Tables for dc DCT coefficients
 127      * Tables are cut in two parts to reduce memory occupation
 128      */
 129
 130     /* Table B-12, dct_dc_size_luminance, codes 00xxx ... 11110 */
 131 lookup_t     pl_dct_dc_lum_init_table_1[32] =
 132     { {1, 2}, {1, 2}, {1, 2}, {1, 2}, {1, 2}, {1, 2}, {1, 2}, {1, 2},
 133       {2, 2}, {2, 2}, {2, 2}, {2, 2}, {2, 2}, {2, 2}, {2, 2}, {2, 2},
 134       {0, 3}, {0, 3}, {0, 3}, {0, 3}, {3, 3}, {3, 3}, {3, 3}, {3, 3},
 135       {4, 3}, {4, 3}, {4, 3}, {4, 3}, {5, 4}, {5, 4}, {6, 5}, {MB_ERROR, 0}
 136     };
 137
 138     /* Table B-12, dct_dc_size_luminance, codes 111110xxx ... 111111111 */
 139 lookup_t     pl_dct_dc_lum_init_table_2[32] =
 140     { {7, 6}, {7, 6}, {7, 6}, {7, 6}, {7, 6}, {7, 6}, {7, 6}, {7, 6},
 141       {8, 7}, {8, 7}, {8, 7}, {8, 7}, {9, 8}, {9, 8}, {10,9}, {11,9},
 142       {MB_ERROR, 0}, {MB_ERROR, 0}, {MB_ERROR, 0}, {MB_ERROR, 0},
 143       {MB_ERROR, 0}, {MB_ERROR, 0}, {MB_ERROR, 0}, {MB_ERROR, 0},
 144       {MB_ERROR, 0}, {MB_ERROR, 0}, {MB_ERROR, 0}, {MB_ERROR, 0},
 145       {MB_ERROR, 0}, {MB_ERROR, 0}, {MB_ERROR, 0}, {MB_ERROR, 0}
 146     };
 147
 148     /* Table B-13, dct_dc_size_chrominance, codes 00xxx ... 11110 */
 149 lookup_t     pl_dct_dc_chrom_init_table_1[32] =
 150     { {0, 2}, {0, 2}, {0, 2}, {0, 2}, {0, 2}, {0, 2}, {0, 2}, {0, 2},
 151       {1, 2}, {1, 2}, {1, 2}, {1, 2}, {1, 2}, {1, 2}, {1, 2}, {1, 2},
 152       {2, 2}, {2, 2}, {2, 2}, {2, 2}, {2, 2}, {2, 2}, {2, 2}, {2, 2},
 153       {3, 3}, {3, 3}, {3, 3}, {3, 3}, {4, 4}, {4, 4}, {5, 5}, {MB_ERROR, 0}
 154     };
 155
 156     /* Table B-13, dct_dc_size_chrominance, codes 111110xxxx ... 1111111111 */
 157 lookup_t     pl_dct_dc_chrom_init_table_2[32] =
 158     { {6, 6}, {6, 6}, {6, 6}, {6, 6}, {6, 6}, {6, 6}, {6, 6}, {6, 6},
 159       {6, 6}, {6, 6}, {6, 6}, {6, 6}, {6, 6}, {6, 6}, {6, 6}, {6, 6},
 160       {7, 7}, {7, 7}, {7, 7}, {7, 7}, {7, 7}, {7, 7}, {7, 7}, {7, 7},
 161       {8, 8}, {8, 8}, {8, 8}, {8, 8}, {9, 9}, {9, 9}, {10,10}, {11,10}
 162     };
 163
 164     /*
 165      * Structure to store the tables B14 & B15
 166      * Is constructed from the tables below
 167      */
 168     dct_lookup_t            ppl_dct_coef[2][16384];
 169
 170
 171     /* Tables for ac DCT coefficients. There are cut in many parts to save space */
 172     /* Table B-14, DCT coefficients table zero,
 173      * codes 0100 ... 1xxx (used for first (DC) coefficient)
 174      */
 175 dct_lookup_t pl_DCT_tab_dc[12] =
 176     {
 177         {0,2,4}, {2,1,4}, {1,1,3}, {1,1,3},
 178         {0,1,1}, {0,1,1}, {0,1,1}, {0,1,1},
 179         {0,1,1}, {0,1,1}, {0,1,1}, {0,1,1}
 180     };
 181
 182     /* Table B-14, DCT coefficients table zero,
 183      * codes 0100 ... 1xxx (used for all other coefficients)
 184      */
 185 dct_lookup_t pl_DCT_tab_ac[12] =
 186     {
 187         {0,2,4},  {2,1,4},  {1,1,3},  {1,1,3},
 188         {DCT_EOB,0,2}, {DCT_EOB,0,2}, {DCT_EOB,0,2}, {DCT_EOB,0,2}, /* EOB */
 189         {0,1,2},  {0,1,2},  {0,1,2},  {0,1,2}
 190     };
 191
 192     /* Table B-14, DCT coefficients table zero,
 193      * codes 000001xx ... 00111xxx
 194      */
 195 dct_lookup_t pl_DCT_tab0[60] =
 196     {
 197         {DCT_ESCAPE,0,6}, {DCT_ESCAPE,0,6}, {DCT_ESCAPE,0,6}, {DCT_ESCAPE,0,6},
 198         /* Escape */
 199         {2,2,7}, {2,2,7}, {9,1,7}, {9,1,7},
 200         {0,4,7}, {0,4,7}, {8,1,7}, {8,1,7},
 201         {7,1,6}, {7,1,6}, {7,1,6}, {7,1,6},
 202         {6,1,6}, {6,1,6}, {6,1,6}, {6,1,6},
 203         {1,2,6}, {1,2,6}, {1,2,6}, {1,2,6},
 204         {5,1,6}, {5,1,6}, {5,1,6}, {5,1,6},
 205         {13,1,8}, {0,6,8}, {12,1,8}, {11,1,8},
 206         {3,2,8}, {1,3,8}, {0,5,8}, {10,1,8},
 207         {0,3,5}, {0,3,5}, {0,3,5}, {0,3,5},
 208         {0,3,5}, {0,3,5}, {0,3,5}, {0,3,5},
 209         {4,1,5}, {4,1,5}, {4,1,5}, {4,1,5},
 210         {4,1,5}, {4,1,5}, {4,1,5}, {4,1,5},
 211         {3,1,5}, {3,1,5}, {3,1,5}, {3,1,5},
 212         {3,1,5}, {3,1,5}, {3,1,5}, {3,1,5}
 213     };
 214
 215     /* Table B-15, DCT coefficients table one,
 216      * codes 000001xx ... 11111111
 217      */
 218 dct_lookup_t pl_DCT_tab0a[252] =
 219     {
 220         {65,0,6}, {65,0,6}, {65,0,6}, {65,0,6}, /* Escape */
 221         {7,1,7}, {7,1,7}, {8,1,7}, {8,1,7},
 222         {6,1,7}, {6,1,7}, {2,2,7}, {2,2,7},
 223         {0,7,6}, {0,7,6}, {0,7,6}, {0,7,6},
 224         {0,6,6}, {0,6,6}, {0,6,6}, {0,6,6},
 225         {4,1,6}, {4,1,6}, {4,1,6}, {4,1,6},
 226         {5,1,6}, {5,1,6}, {5,1,6}, {5,1,6},
 227         {1,5,8}, {11,1,8}, {0,11,8}, {0,10,8},
 228         {13,1,8}, {12,1,8}, {3,2,8}, {1,4,8},
 229         {2,1,5}, {2,1,5}, {2,1,5}, {2,1,5},
 230         {2,1,5}, {2,1,5}, {2,1,5}, {2,1,5},
 231         {1,2,5}, {1,2,5}, {1,2,5}, {1,2,5},
 232         {1,2,5}, {1,2,5}, {1,2,5}, {1,2,5},
 233         {3,1,5}, {3,1,5}, {3,1,5}, {3,1,5},
 234         {3,1,5}, {3,1,5}, {3,1,5}, {3,1,5},
 235         {1,1,3}, {1,1,3}, {1,1,3}, {1,1,3},
 236         {1,1,3}, {1,1,3}, {1,1,3}, {1,1,3},
 237         {1,1,3}, {1,1,3}, {1,1,3}, {1,1,3},
 238         {1,1,3}, {1,1,3}, {1,1,3}, {1,1,3},
 239         {1,1,3}, {1,1,3}, {1,1,3}, {1,1,3},
 240         {1,1,3}, {1,1,3}, {1,1,3}, {1,1,3},
 241         {1,1,3}, {1,1,3}, {1,1,3}, {1,1,3},
 242         {1,1,3}, {1,1,3}, {1,1,3}, {1,1,3},
 243         {64,0,4}, {64,0,4}, {64,0,4}, {64,0,4}, /* EOB */
 244         {64,0,4}, {64,0,4}, {64,0,4}, {64,0,4},
 245         {64,0,4}, {64,0,4}, {64,0,4}, {64,0,4},
 246         {64,0,4}, {64,0,4}, {64,0,4}, {64,0,4},
 247         {0,3,4}, {0,3,4}, {0,3,4}, {0,3,4},
 248         {0,3,4}, {0,3,4}, {0,3,4}, {0,3,4},
 249         {0,3,4}, {0,3,4}, {0,3,4}, {0,3,4},
 250         {0,3,4}, {0,3,4}, {0,3,4}, {0,3,4},
 251         {0,1,2}, {0,1,2}, {0,1,2}, {0,1,2},
 252         {0,1,2}, {0,1,2}, {0,1,2}, {0,1,2},
 253         {0,1,2}, {0,1,2}, {0,1,2}, {0,1,2},
 254         {0,1,2}, {0,1,2}, {0,1,2}, {0,1,2},
 255         {0,1,2}, {0,1,2}, {0,1,2}, {0,1,2},
 256         {0,1,2}, {0,1,2}, {0,1,2}, {0,1,2},
 257         {0,1,2}, {0,1,2}, {0,1,2}, {0,1,2},
 258         {0,1,2}, {0,1,2}, {0,1,2}, {0,1,2},
 259         {0,1,2}, {0,1,2}, {0,1,2}, {0,1,2},
 260         {0,1,2}, {0,1,2}, {0,1,2}, {0,1,2},
 261         {0,1,2}, {0,1,2}, {0,1,2}, {0,1,2},
 262         {0,1,2}, {0,1,2}, {0,1,2}, {0,1,2},
 263         {0,1,2}, {0,1,2}, {0,1,2}, {0,1,2},
 264         {0,1,2}, {0,1,2}, {0,1,2}, {0,1,2},
 265         {0,1,2}, {0,1,2}, {0,1,2}, {0,1,2},
 266         {0,1,2}, {0,1,2}, {0,1,2}, {0,1,2},
 267         {0,2,3}, {0,2,3}, {0,2,3}, {0,2,3},
 268         {0,2,3}, {0,2,3}, {0,2,3}, {0,2,3},
 269         {0,2,3}, {0,2,3}, {0,2,3}, {0,2,3},
 270         {0,2,3}, {0,2,3}, {0,2,3}, {0,2,3},
 271         {0,2,3}, {0,2,3}, {0,2,3}, {0,2,3},
 272         {0,2,3}, {0,2,3}, {0,2,3}, {0,2,3},
 273         {0,2,3}, {0,2,3}, {0,2,3}, {0,2,3},
 274         {0,2,3}, {0,2,3}, {0,2,3}, {0,2,3},
 275         {0,4,5}, {0,4,5}, {0,4,5}, {0,4,5},
 276         {0,4,5}, {0,4,5}, {0,4,5}, {0,4,5},
 277         {0,5,5}, {0,5,5}, {0,5,5}, {0,5,5},
 278         {0,5,5}, {0,5,5}, {0,5,5}, {0,5,5},
 279         {9,1,7}, {9,1,7}, {1,3,7}, {1,3,7},
 280         {10,1,7}, {10,1,7}, {0,8,7}, {0,8,7},
 281         {0,9,7}, {0,9,7}, {0,12,8}, {0,13,8},
 282         {2,3,8}, {4,2,8}, {0,14,8}, {0,15,8}
 283     };
 284
 285
 286
 287
 288     /* Table B-14, DCT coefficients table zero,
 289      * codes 0000001000 ... 0000001111
 290      */
 291 dct_lookup_t pl_DCT_tab1[8] =
 292     {
 293         {16,1,10}, {5,2,10}, {0,7,10}, {2,3,10},
 294         {1,4,10}, {15,1,10}, {14,1,10}, {4,2,10}
 295     };
 296
 297     /* Table B-15, DCT coefficients table one,
 298      * codes 000000100x ... 000000111x
 299      */
 300 dct_lookup_t pl_DCT_tab1a[8] =
 301     {
 302         {5,2,9}, {5,2,9}, {14,1,9}, {14,1,9},
 303         {2,4,10}, {16,1,10}, {15,1,9}, {15,1,9}
 304     };
 305
 306     /* Table B-14/15, DCT coefficients table zero / one,
 307      * codes 000000010000 ... 000000011111
 308      */
 309 dct_lookup_t pl_DCT_tab2[16] =
 310     {
 311         {0,11,12}, {8,2,12}, {4,3,12}, {0,10,12},
 312         {2,4,12}, {7,2,12}, {21,1,12}, {20,1,12},
 313         {0,9,12}, {19,1,12}, {18,1,12}, {1,5,12},
 314         {3,3,12}, {0,8,12}, {6,2,12}, {17,1,12}
 315     };
 316
 317     /* Table B-14/15, DCT coefficients table zero / one,
 318      * codes 0000000010000 ... 0000000011111
 319      */
 320 dct_lookup_t pl_DCT_tab3[16] =
 321     {
 322         {10,2,13}, {9,2,13}, {5,3,13}, {3,4,13},
 323         {2,5,13}, {1,7,13}, {1,6,13}, {0,15,13},
 324         {0,14,13}, {0,13,13}, {0,12,13}, {26,1,13},
 325         {25,1,13}, {24,1,13}, {23,1,13}, {22,1,13}
 326     };
 327
 328     /* Table B-14/15, DCT coefficients table zero / one,
 329      * codes 00000000010000 ... 00000000011111
 330      */
 331 dct_lookup_t pl_DCT_tab4[16] =
 332     {
 333         {0,31,14}, {0,30,14}, {0,29,14}, {0,28,14},
 334         {0,27,14}, {0,26,14}, {0,25,14}, {0,24,14},
 335         {0,23,14}, {0,22,14}, {0,21,14}, {0,20,14},
 336         {0,19,14}, {0,18,14}, {0,17,14}, {0,16,14}
 337     };
 338
 339     /* Table B-14/15, DCT coefficients table zero / one,
 340      *   codes 000000000010000 ... 000000000011111
 341      */
 342 dct_lookup_t pl_DCT_tab5[16] =
 343     {
 344     {0,40,15}, {0,39,15}, {0,38,15}, {0,37,15},
 345     {0,36,15}, {0,35,15}, {0,34,15}, {0,33,15},
 346     {0,32,15}, {1,14,15}, {1,13,15}, {1,12,15},
 347     {1,11,15}, {1,10,15}, {1,9,15}, {1,8,15}
 348     };
 349
 350     /* Table B-14/15, DCT coefficients table zero / one,
 351      * codes 0000000000010000 ... 0000000000011111
 352      */
 353 dct_lookup_t pl_DCT_tab6[16] =
 354     {
 355         {1,18,16}, {1,17,16}, {1,16,16}, {1,15,16},
 356         {6,3,16}, {16,2,16}, {15,2,16}, {14,2,16},
 357         {13,2,16}, {12,2,16}, {11,2,16}, {31,1,16},
 358         {30,1,16}, {29,1,16}, {28,1,16}, {27,1,16}
 359     };
 360 /*
 361  * Initialization of lookup tables
 362  */
 363
 364 /*****************************************************************************
 365  * vpar_InitCrop : Initialize the crop table for saturation
 366  *                 (ISO/IEC 13818-2 section 7.4.3)
 367  *****************************************************************************/
 368 #if defined(MPEG2_COMPLIANT) && !defined(VDEC_DFT)
 369 void vpar_InitCrop( vpar_thread_t * p_vpar )
 370 {
 371     int i_dummy;
 372
 373     p_vpar->pi_crop = p_vpar->pi_crop_buf + 4096;
 374
 375     for( i_dummy = -4096; i_dummy < -2048; i_dummy++ )
 376     {
 377         p_vpar->pi_crop[i_dummy] = -2048;
 378     }
 379     for( ; i_dummy < 2047; i_dummy++ )
 380     {
 381         p_vpar->pi_crop[i_dummy] = i_dummy;
 382     }
 383     for( ; i_dummy < 4095; i_dummy++ )
 384     {
 385         p_vpar->pi_crop[i_dummy] = 2047;
 386     }
 387 }
 388 #endif
 389
 390 /*****************************************************************************
 391  * InitMbAddrInc : Initialize the lookup table for mb_addr_inc
 392  *****************************************************************************/
 393
 394 /* Function for filling up the lookup table for mb_addr_inc */
 395 void __inline__ FillMbAddrIncTable( vpar_thread_t * p_vpar,
 396                                     int i_start, int i_end, int i_step,
 397                                     int * pi_value, int i_length )
 398 {
 399     int i_pos, i_offset;
 400     for( i_pos = i_start ; i_pos < i_end ; i_pos += i_step )
 401     {
 402         for( i_offset = 0 ; i_offset < i_step ; i_offset ++ )
 403         {
 404             p_vpar->pl_mb_addr_inc[i_pos + i_offset].i_value = * pi_value;
 405             p_vpar->pl_mb_addr_inc[i_pos + i_offset].i_length = i_length;
 406         }
 407         (*pi_value)--;
 408     }
 409 }
 410
 411 /* Function that initialize the table using the last one */
 412 void vpar_InitMbAddrInc( vpar_thread_t * p_vpar )
 413 {
 414     int i_dummy;
 415     int i_value;
 416
 417     for( i_dummy = 0 ; i_dummy < 8 ; i_dummy++ )
 418     {
 419         p_vpar->pl_mb_addr_inc[i_dummy].i_value = MB_ERROR;
 420         p_vpar->pl_mb_addr_inc[i_dummy].i_length = 0;
 421     }
 422
 423     p_vpar->pl_mb_addr_inc[8].i_value = MB_ADDRINC_ESCAPE;
 424     p_vpar->pl_mb_addr_inc[8].i_length = 11;
 425
 426     for( i_dummy = 9 ; i_dummy < 15 ; i_dummy ++ )
 427     {
 428         p_vpar->pl_mb_addr_inc[i_dummy].i_value =  MB_ERROR;
 429         p_vpar->pl_mb_addr_inc[i_dummy].i_length = 0;
 430     }
 431
 432     p_vpar->pl_mb_addr_inc[15].i_value = MB_ADDRINC_STUFFING;
 433     p_vpar->pl_mb_addr_inc[15].i_length = 11;
 434
 435     for( i_dummy = 16; i_dummy < 24; i_dummy++ )
 436     {
 437         p_vpar->pl_mb_addr_inc[i_dummy].i_value =  MB_ERROR;
 438         p_vpar->pl_mb_addr_inc[i_dummy].i_length = 0;
 439     }
 440
 441     i_value = 33;
 442
 443     FillMbAddrIncTable( p_vpar, 24, 36, 1, &i_value, 11 );
 444     FillMbAddrIncTable( p_vpar, 36, 48, 2, &i_value, 10 );
 445     FillMbAddrIncTable( p_vpar, 48, 96, 8, &i_value, 8 );
 446     FillMbAddrIncTable( p_vpar, 96, 128, 16, &i_value, 7 );
 447     FillMbAddrIncTable( p_vpar, 128, 256, 64, &i_value, 5 );
 448     FillMbAddrIncTable( p_vpar, 256, 512, 128, &i_value, 4 );
 449     FillMbAddrIncTable( p_vpar, 512, 1024, 256, &i_value, 3 );
 450     FillMbAddrIncTable( p_vpar, 1024, 2048, 1024, &i_value, 1 );
 451 }
 452
 453 /*****************************************************************************
 454  * Init*MBType : Initialize lookup table for the Macroblock type
 455  *****************************************************************************/
 456
 457 /* Fonction for filling up the tables */
 458 static void __inline__ FillMBType( vpar_thread_t * p_vpar,
 459                                    int           i_mb_type,
 460                                    int           i_start,
 461                                    int           i_end,
 462                                    int           i_value,
 463                                    int           i_length )
 464 {
 465     int i_dummy;
 466
 467     for( i_dummy = i_start ; i_dummy < i_end ; i_dummy++ )
 468     {
 469         p_vpar->ppl_mb_type[i_mb_type][i_dummy].i_value = i_value;
 470         p_vpar->ppl_mb_type[i_mb_type][i_dummy].i_length = i_length;
 471     }
 472 }
 473
 474 /* Fonction that fills the table for P MB_Type */
 475 void vpar_InitPMBType( vpar_thread_t * p_vpar )
 476 {
 477     FillMBType( p_vpar, 0, 32, 64, MB_MOTION_FORWARD|MB_PATTERN, 1 );
 478     FillMBType( p_vpar, 0, 16, 32, MB_PATTERN, 2 );
 479     FillMBType( p_vpar, 0, 8, 16, MB_MOTION_FORWARD, 3 );
 480     FillMBType( p_vpar, 0, 6, 8, MB_INTRA, 5 );
 481     FillMBType( p_vpar, 0, 4, 6, MB_QUANT|MB_MOTION_FORWARD|MB_PATTERN, 5 );
 482     FillMBType( p_vpar, 0, 2, 4, MB_QUANT|MB_PATTERN, 5 );
 483     p_vpar->ppl_mb_type[0][1].i_value = MB_QUANT|MB_INTRA;
 484     p_vpar->ppl_mb_type[0][1].i_length = 6;
 485     p_vpar->ppl_mb_type[0][0].i_value = MB_ERROR;
 486     p_vpar->ppl_mb_type[0][0].i_length = 0;
 487 }
 488
 489 /* Fonction that fills the table for B MB_Type */
 490 void vpar_InitBMBType( vpar_thread_t * p_vpar )
 491 {
 492     FillMBType( p_vpar, 1, 48, 64, MB_MOTION_FORWARD
 493                                   |MB_MOTION_BACKWARD|MB_PATTERN, 2 );
 494     FillMBType( p_vpar, 1, 32, 48, MB_MOTION_FORWARD|MB_MOTION_BACKWARD, 2 );
 495     FillMBType( p_vpar, 1, 24, 32, MB_MOTION_BACKWARD|MB_PATTERN, 3 );
 496     FillMBType( p_vpar, 1, 16, 24, MB_MOTION_BACKWARD, 3 );
 497     FillMBType( p_vpar, 1, 12, 16, MB_MOTION_FORWARD|MB_PATTERN, 4 );
 498     FillMBType( p_vpar, 1, 8, 12, MB_MOTION_FORWARD, 4 );
 499     FillMBType( p_vpar, 1, 6, 8, MB_INTRA, 5 );
 500     FillMBType( p_vpar, 1, 4, 6, MB_QUANT|MB_MOTION_FORWARD
 501                                 |MB_MOTION_BACKWARD|MB_PATTERN, 5 );
 502     p_vpar->ppl_mb_type[1][3].i_value = MB_QUANT|MB_MOTION_FORWARD|MB_PATTERN;
 503     p_vpar->ppl_mb_type[1][3].i_length = 6;
 504     p_vpar->ppl_mb_type[1][2].i_value = MB_QUANT|MB_MOTION_BACKWARD|MB_PATTERN;
 505     p_vpar->ppl_mb_type[1][2].i_length = 6;
 506     p_vpar->ppl_mb_type[1][1].i_value = MB_QUANT|MB_INTRA;
 507     p_vpar->ppl_mb_type[1][1].i_length = 6;
 508     p_vpar->ppl_mb_type[1][0].i_value =MB_ERROR;
 509     p_vpar->ppl_mb_type[1][0].i_length = 0;
 510 }
 511
 512 /*****************************************************************************
 513  * InitCodedPattern : Initialize the lookup table for decoding
 514  *                    coded block pattern
 515  *****************************************************************************/
 516 void vpar_InitCodedPattern( vpar_thread_t * p_vpar )
 517 {
 518     p_vpar->pl_coded_pattern = (lookup_t*) pl_coded_pattern_init_table;
 519 }
 520
 521 /*****************************************************************************
 522  * InitDCT : Initialize tables giving the length of the dct coefficient
 523  *           from the vlc code
 524  *****************************************************************************/
 525
 526 /* First fonction for filling the table */
 527 static void __inline__ FillDCTTable( dct_lookup_t * p_tab_dest, dct_lookup_t * p_tab_src,
 528                                      int i_step, int i_nb_elem, int i_offset )
 529 {
 530     int i_dummy, i_dummy2;
 531
 532     for( i_dummy=0 ; i_dummy < i_nb_elem ; i_dummy++ )
 533     {
 534         for( i_dummy2=0 ; i_dummy2 < i_step ; i_dummy2++ )
 535         {
 536             p_tab_dest[(i_dummy+i_offset)*i_step+i_dummy2] = p_tab_src[i_dummy];
 537         }
 538     }
 539 }
 540
 541
 542 /* Fonction that actually fills the table or create the pointers */
 543 void vpar_InitDCTTables( vpar_thread_t * p_vpar )
 544 {
 545     /* Tables are cut in two parts to reduce memory occupation */
 546     p_vpar->pppl_dct_dc_size[0][0] = pl_dct_dc_lum_init_table_1;
 547     p_vpar->pppl_dct_dc_size[0][1] = pl_dct_dc_lum_init_table_2;
 548     p_vpar->pppl_dct_dc_size[1][0] = pl_dct_dc_chrom_init_table_1;
 549     p_vpar->pppl_dct_dc_size[1][1] = pl_dct_dc_chrom_init_table_2;
 550
 551     memset( ppl_dct_coef[0], MB_ERROR, 16 );
 552     memset( ppl_dct_coef[1], MB_ERROR, 16 );
 553
 554     /* For table B14 & B15, we have a pointer to tables */
 555     /* We fill the table thanks to the fonction defined above */
 556     FillDCTTable( ppl_dct_coef[0], pl_DCT_tab0, 256, 60,  4 );
 557     FillDCTTable( ppl_dct_coef[0], pl_DCT_tab1,  64,  8,  8 );
 558     FillDCTTable( ppl_dct_coef[0], pl_DCT_tab2,  16, 16, 16 );
 559     FillDCTTable( ppl_dct_coef[0], pl_DCT_tab3,   8, 16, 16 );
 560     FillDCTTable( ppl_dct_coef[0], pl_DCT_tab4,   4, 16, 16 );
 561     FillDCTTable( ppl_dct_coef[0], pl_DCT_tab5,   2, 16, 16 );
 562     FillDCTTable( ppl_dct_coef[0], pl_DCT_tab6,   1, 16, 16 );
 563
 564     FillDCTTable( ppl_dct_coef[1], pl_DCT_tab0a, 256, 60, 4 );
 565     FillDCTTable( ppl_dct_coef[1], pl_DCT_tab1a,  64,  8,  8 );
 566     FillDCTTable( ppl_dct_coef[1], pl_DCT_tab2,   16, 16, 16 );
 567     FillDCTTable( ppl_dct_coef[1], pl_DCT_tab3,    8, 16, 16 );
 568     FillDCTTable( ppl_dct_coef[1], pl_DCT_tab4,    4, 16, 16 );
 569     FillDCTTable( ppl_dct_coef[1], pl_DCT_tab5,    2, 16, 16 );
 570     FillDCTTable( ppl_dct_coef[1], pl_DCT_tab6,    1, 16, 16 );
 571 }
 572
 573 /*
 574  * Macroblock parsing functions
 575  */
 576
 577 /*****************************************************************************
 578  * vpar_ParseMacroblock : Parse the next macroblock
 579  *****************************************************************************/
 580 void vpar_ParseMacroblock( vpar_thread_t * p_vpar, int * pi_mb_address,
 581                            int i_mb_previous, int i_mb_base )
 582 {
 583     static f_addb_t ppf_addb_intra[2] = {vdec_AddBlock, vdec_CopyBlock};
 584     static f_decode_block_t pppf_decode_block[2][2] =
 585                 { {vpar_DecodeMPEG1Non, vpar_DecodeMPEG1Intra},
 586                   {vpar_DecodeMPEG2Non, vpar_DecodeMPEG2Intra} };
 587     static int      pi_x[12] = {0,8,0,8,0,0,0,0,8,8,8,8};
 588     static int      pi_y[2][12] = { {0,0,8,8,0,0,8,8,0,0,8,8},
 589                                     {0,0,1,1,0,0,1,1,0,0,1,1} };
 590     static int      pi_chroma_hor[4] = { 0, 1, 1, 0 };
 591     static int      pi_chroma_ver[4] = { 0, 1, 0, 0 };
 592
 593     int             i_mb, i_b, i_mask;
 594     macroblock_t *  p_mb;
 595     f_addb_t        pf_addb;
 596     yuv_data_t *    p_data1;
 597     yuv_data_t *    p_data2;
 598
 599     /************* DEBUG *************/
 600     static int i_count = 0;
 601     int i_inc;
 602     boolean_t b_stop = 0;
 603
 604     i_inc = MacroblockAddressIncrement( p_vpar );
 605
 606     *pi_mb_address += i_inc;
 607
 608     //*pi_mb_address += MacroblockAddressIncrement( p_vpar );
 609     b_stop = i_inc > 1;
 610     fprintf( stderr, "inc : %d (%d)\n", *pi_mb_address, i_inc );
 611 if( 0  )//i_count > 8 )
 612 {
 613     exit(0);
 614 }
 615 i_count++;
 616
 617 if( 0 )
 618     //i_count == 249)
 619     // i_count != *pi_mb_address)
 620     //b_stop )
 621 {
 622     fprintf( stderr, "i_count = %d (%d)\n", i_count, i_inc );
 623      fprintf( stderr, "%x", GetBits( &p_vpar->bit_stream, 16 ) );
 624      fprintf( stderr, "%x ", GetBits( &p_vpar->bit_stream, 16 ) );
 625      fprintf( stderr, "%x", GetBits( &p_vpar->bit_stream, 16 ) );
 626      fprintf( stderr, "%x\n", GetBits( &p_vpar->bit_stream, 16 ) );
 627      fprintf( stderr, "%x", GetBits( &p_vpar->bit_stream, 16 ) );
 628      fprintf( stderr, "%x ", GetBits( &p_vpar->bit_stream, 16 ) );
 629      fprintf( stderr, "%x", GetBits( &p_vpar->bit_stream, 16 ) );
 630      fprintf( stderr, "%x\n", GetBits( &p_vpar->bit_stream, 16 ) );
 631      exit(0);
 632 }
 633
 634     for( i_mb = i_mb_previous; i_mb < *pi_mb_address; i_mb++ )
 635     {
 636         /* Skipped macroblock (ISO/IEC 13818-2 7.6.6). */
 637         static int          pi_dc_dct_reinit[4] = {128,256,512,1024};
 638         static f_motion_t   pf_motion_skipped[4] = {NULL, vdec_MotionField,
 639                                 vdec_MotionField, vdec_MotionFrame};
 640
 641         /* Reset DC predictors (7.2.1). */
 642         p_vpar->slice.pi_dc_dct_pred[0] = p_vpar->slice.pi_dc_dct_pred[1]
 643             = p_vpar->slice.pi_dc_dct_pred[2]
 644             = pi_dc_dct_reinit[p_vpar->picture.i_intra_dc_precision];
 645
 646         if( p_vpar->picture.i_coding_type == P_CODING_TYPE )
 647         {
 648             /* Reset motion vector predictors (ISO/IEC 13818-2 7.6.3.4). */
 649             bzero( p_vpar->slice.pppi_pmv, 8*sizeof(int) );
 650         }
 651
 652         if( (p_mb = p_vpar->picture.pp_mb[i_mb_base + i_mb] =
 653              vpar_NewMacroblock( &p_vpar->vfifo )) == NULL )
 654         {
 655             p_vpar->picture.b_error = 1;
 656             intf_ErrMsg("vpar error: macroblock list is empty !\n");
 657             return;
 658         }
 659
 660         InitMacroblock( p_vpar, p_mb );
 661
 662         /* No IDCT nor AddBlock. */
 663         for( i_b = 0; i_b < 12; i_b++ )
 664         {
 665             p_mb->pf_idct[i_b] = vdec_DummyIDCT;
 666             p_mb->pf_addb[i_b] = vdec_DummyBlock;
 667         }
 668
 669         /* Motion type is picture structure. */
 670         p_mb->pf_motion = pf_motion_skipped[p_vpar->picture.i_structure];
 671
 672         /* Set the field we use for motion compensation */
 673         p_mb->ppi_field_select[0][0] = p_mb->ppi_field_select[0][1]
 674                                      = ( p_vpar->picture.i_current_structure == BOTTOM_FIELD );
 675     }
 676 //fprintf(stderr, "MB1\n");
 677     /* Get a macroblock structure. */
 678     if( (p_mb = p_vpar->picture.pp_mb[i_mb_base + *pi_mb_address] =
 679          vpar_NewMacroblock( &p_vpar->vfifo )) == NULL )
 680     {
 681         p_vpar->picture.b_error = 1;
 682         intf_ErrMsg("vpar error: macroblock list is empty !\n");
 683         return;
 684     }
 685
 686     InitMacroblock( p_vpar, p_mb );
 687
 688     /* Parse off macroblock_modes structure. */
 689     MacroblockModes( p_vpar, p_mb );
 690
 691     if( p_vpar->mb.i_mb_type & MB_QUANT )
 692     {
 693         LoadQuantizerScale( p_vpar );
 694     }
 695
 696     if( p_vpar->mb.i_mb_type & MB_MOTION_FORWARD )
 697     {
 698 fprintf( stderr, "motion !\n" );
 699         (*p_vpar->sequence.pf_decode_mv)( p_vpar, p_mb, 0 );
 700     }
 701
 702     if( p_vpar->mb.i_mb_type & MB_MOTION_BACKWARD )
 703     {
 704 fprintf( stderr, "motion2 !\n" );
 705         (*p_vpar->sequence.pf_decode_mv)( p_vpar, p_mb, 1 );
 706     }
 707
 708     if( p_vpar->picture.b_concealment_mv && (p_vpar->mb.i_mb_type & MB_INTRA) )
 709     {
 710         RemoveBits( &p_vpar->bit_stream, 1 );
 711     }
 712
 713     if( p_vpar->mb.i_mb_type & MB_PATTERN )
 714     {
 715         (*p_vpar->sequence.pf_decode_pattern)( p_vpar );
 716     }
 717     else
 718     {
 719         int     pi_coded_block_pattern[2] = {0,
 720                     (1 << (4+p_vpar->sequence.i_chroma_nb_blocks)) - 1};
 721         p_vpar->mb.i_coded_block_pattern = pi_coded_block_pattern
 722                                     [p_vpar->mb.i_mb_type & MB_INTRA];
 723     }
 724
 725     pf_addb = ppf_addb_intra[p_vpar->mb.i_mb_type & MB_INTRA];
 726
 727     /*
 728      * Effectively decode blocks.
 729      */
 730 //fprintf(stderr, "MB2\n");
 731     i_mask = 1 << (3 + p_vpar->sequence.i_chroma_nb_blocks);
 732
 733     /* luminance */
 734     p_data1 = p_mb->p_picture->p_y
 735               + p_mb->i_l_x + p_mb->i_l_y*(p_vpar->sequence.i_width);
 736
 737     for( i_b = 0 ; i_b < 4 ; i_b++, i_mask >>= 1 )
 738     {
 739         if( p_vpar->mb.i_coded_block_pattern & i_mask )
 740         {
 741             memset( p_mb->ppi_blocks[i_b], 0, 64*sizeof(elem_t) );
 742             (*pppf_decode_block[p_vpar->sequence.b_mpeg2]
 743                                [p_vpar->mb.i_mb_type & MB_INTRA])
 744                 ( p_vpar, p_mb, i_b );
 745
 746             /* decode_block has already set pf_idct and pi_sparse_pos. */
 747             p_mb->pf_addb[i_b] = pf_addb;
 748
 749             /* Calculate block coordinates. */
 750             p_mb->p_data[i_b] = p_data1
 751                                 + pi_y[p_vpar->mb.b_dct_type][i_b]
 752                                 * p_vpar->sequence.i_width;
 753         }
 754         else
 755         {
 756             /* Block not coded, so no IDCT, nor AddBlock */
 757             p_mb->pf_addb[i_b] = vdec_DummyBlock;
 758             p_mb->pf_idct[i_b] = vdec_DummyIDCT;
 759         }
 760  if( 0 )
 761     //i_count == 482 && i_b == 2 )
 762     // i_count != *pi_mb_address)
 763     //b_stop )
 764 {
 765     fprintf( stderr, "i_count = %d (%d)\n", i_count, i_inc );
 766      fprintf( stderr, "%4x", GetBits( &p_vpar->bit_stream, 16 ) );
 767      fprintf( stderr, "%4x ", GetBits( &p_vpar->bit_stream, 16 ) );
 768      fprintf( stderr, "%4x", GetBits( &p_vpar->bit_stream, 16 ) );
 769      fprintf( stderr, "%4x\n", GetBits( &p_vpar->bit_stream, 16 ) );
 770      fprintf( stderr, "%4x", GetBits( &p_vpar->bit_stream, 16 ) );
 771      fprintf( stderr, "%4x ", GetBits( &p_vpar->bit_stream, 16 ) );
 772      fprintf( stderr, "%4x", GetBits( &p_vpar->bit_stream, 16 ) );
 773      fprintf( stderr, "%4x\n", GetBits( &p_vpar->bit_stream, 16 ) );
 774      fprintf( stderr, "%4x", GetBits( &p_vpar->bit_stream, 16 ) );
 775      fprintf( stderr, "%4x ", GetBits( &p_vpar->bit_stream, 16 ) );
 776      fprintf( stderr, "%4x", GetBits( &p_vpar->bit_stream, 16 ) );
 777      fprintf( stderr, "%4x\n", GetBits( &p_vpar->bit_stream, 16 ) );
 778      fprintf( stderr, "%4x", GetBits( &p_vpar->bit_stream, 16 ) );
 779      fprintf( stderr, "%4x ", GetBits( &p_vpar->bit_stream, 16 ) );
 780      fprintf( stderr, "%4x", GetBits( &p_vpar->bit_stream, 16 ) );
 781      fprintf( stderr, "%4x\n", GetBits( &p_vpar->bit_stream, 16 ) );
 782      fprintf( stderr, "%4x", GetBits( &p_vpar->bit_stream, 16 ) );
 783      fprintf( stderr, "%4x ", GetBits( &p_vpar->bit_stream, 16 ) );
 784      fprintf( stderr, "%4x", GetBits( &p_vpar->bit_stream, 16 ) );
 785      fprintf( stderr, "%4x\n", GetBits( &p_vpar->bit_stream, 16 ) );
 786      exit(0);
 787 }
 788     }
 789
 790     /* chrominance */
 791     p_data1 = p_mb->p_picture->p_u
 792               + (p_mb->i_c_x >> pi_chroma_hor[p_vpar->sequence.i_chroma_format])
 793               + (p_mb->i_c_y >> pi_chroma_ver[p_vpar->sequence.i_chroma_format])
 794                 * (p_vpar->sequence.i_chroma_width);
 795     p_data2 = p_mb->p_picture->p_v
 796               + (p_mb->i_c_x >> pi_chroma_hor[p_vpar->sequence.i_chroma_format])
 797               + (p_mb->i_c_y >> pi_chroma_ver[p_vpar->sequence.i_chroma_format])
 798                 * (p_vpar->sequence.i_chroma_width);
 799
 800     for( i_b = 4; i_b < 4 + p_vpar->sequence.i_chroma_nb_blocks;
 801          i_b++, i_mask >>= 1 )
 802     {
 803         yuv_data_t *    pp_data[2] = {p_data1, p_data2};
 804
 805         if( p_vpar->mb.i_coded_block_pattern & i_mask )
 806         {
 807             memset( p_mb->ppi_blocks[i_b], 0, 64*sizeof(elem_t) );
 808             (*pppf_decode_block[p_vpar->sequence.b_mpeg2]
 809                                [p_vpar->mb.i_mb_type & MB_INTRA])
 810                 ( p_vpar, p_mb, i_b );
 811
 812             /* decode_block has already set pf_idct and pi_sparse_pos. */
 813             p_mb->pf_addb[i_b] = pf_addb;
 814
 815             /* Calculate block coordinates. */
 816             p_mb->p_data[i_b] = pp_data[i_b & 1]
 817                                  + pi_y[p_vpar->mb.b_dct_type][i_b]
 818                                    * p_vpar->sequence.i_chroma_width;
 819         }
 820         else
 821         {
 822             /* Block not coded, so no IDCT, nor AddBlock */
 823             p_mb->pf_addb[i_b] = vdec_DummyBlock;
 824             p_mb->pf_idct[i_b] = vdec_DummyIDCT;
 825         }
 826     }
 827
 828     if( !( p_vpar->mb.i_mb_type & MB_INTRA ) )
 829     {
 830         static int          pi_dc_dct_reinit[4] = {128,256,512,1024};
 831
 832         /* Reset DC predictors (7.2.1). */
 833         p_vpar->slice.pi_dc_dct_pred[0] = p_vpar->slice.pi_dc_dct_pred[1]
 834             = p_vpar->slice.pi_dc_dct_pred[2]
 835             = pi_dc_dct_reinit[p_vpar->picture.i_intra_dc_precision];
 836     }
 837     else if( !p_vpar->picture.b_concealment_mv )
 838     {
 839         /* Reset MV predictors. */
 840         bzero( p_vpar->slice.pppi_pmv, 8*sizeof(int) );
 841     }
 842  if( 0 )
 843     //i_count == 249)
 844     // i_count != *pi_mb_address)
 845     //b_stop )
 846 {
 847     fprintf( stderr, "i_count = %d (%d)\n", i_count, i_inc );
 848      fprintf( stderr, "%x", GetBits( &p_vpar->bit_stream, 16 ) );
 849      fprintf( stderr, "%x ", GetBits( &p_vpar->bit_stream, 16 ) );
 850      fprintf( stderr, "%x", GetBits( &p_vpar->bit_stream, 16 ) );
 851      fprintf( stderr, "%x\n", GetBits( &p_vpar->bit_stream, 16 ) );
 852      fprintf( stderr, "%x", GetBits( &p_vpar->bit_stream, 16 ) );
 853      fprintf( stderr, "%x ", GetBits( &p_vpar->bit_stream, 16 ) );
 854      fprintf( stderr, "%x", GetBits( &p_vpar->bit_stream, 16 ) );
 855      fprintf( stderr, "%x\n", GetBits( &p_vpar->bit_stream, 16 ) );
 856      exit(0);
 857 }
 858 }
 859
 860
 861 /*****************************************************************************
 862  * InitMacroblock : Initialize macroblock values
 863  *****************************************************************************/
 864 static __inline__ void InitMacroblock( vpar_thread_t * p_vpar,
 865                                        macroblock_t * p_mb )
 866 {
 867     static f_chroma_motion_t pf_chroma_motion[4] =
 868     { NULL, vdec_Motion420, vdec_Motion422, vdec_Motion444 };
 869
 870     p_mb->p_picture = p_vpar->picture.p_picture;
 871     p_mb->i_structure = p_vpar->picture.i_structure;
 872     p_mb->i_current_structure = p_vpar->picture.i_current_structure;
 873     p_mb->i_l_x = p_vpar->mb.i_l_x;
 874     p_mb->i_l_y = p_vpar->mb.i_l_y;
 875     p_mb->i_c_x = p_vpar->mb.i_c_x;
 876     p_mb->i_c_y = p_vpar->mb.i_c_y;
 877     p_mb->i_chroma_nb_blocks = p_vpar->sequence.i_chroma_nb_blocks;
 878     p_mb->pf_chroma_motion = pf_chroma_motion[p_vpar->sequence.i_chroma_format];
 879     p_mb->b_P_coding_type = ( p_vpar->picture.i_coding_type == P_CODING_TYPE );
 880
 881     p_mb->p_forward = p_vpar->sequence.p_forward;
 882     p_mb->p_backward = p_vpar->sequence.p_backward;
 883
 884     p_mb->i_addb_l_stride = p_mb->i_l_stride = p_vpar->picture.i_l_stride;
 885     p_mb->i_addb_c_stride = p_mb->i_c_stride = p_vpar->picture.i_c_stride;
 886
 887     /* Update macroblock real position. */
 888     p_vpar->mb.i_l_x += 16;
 889     p_vpar->mb.i_l_y += (p_vpar->mb.i_l_x / p_vpar->sequence.i_width)
 890                         * (2 - p_vpar->picture.b_frame_structure) * 16;
 891     p_vpar->mb.i_l_x %= p_vpar->sequence.i_width;
 892
 893     p_vpar->mb.i_c_x += p_vpar->sequence.i_chroma_mb_width;
 894     p_vpar->mb.i_c_y += (p_vpar->mb.i_c_x / p_vpar->sequence.i_chroma_width)
 895                         * (2 - p_vpar->picture.b_frame_structure)
 896                         * p_vpar->sequence.i_chroma_mb_height;
 897     p_vpar->mb.i_c_x %= p_vpar->sequence.i_chroma_width;
 898 }
 899
 900 /*****************************************************************************
 901  * MacroblockAddressIncrement : Get the macroblock_address_increment field
 902  *****************************************************************************/
 903 static __inline__ int MacroblockAddressIncrement( vpar_thread_t * p_vpar )
 904 {
 905     int i_addr_inc = 0;
 906     /* Index in the lookup table mb_addr_inc */
 907     int    i_index = ShowBits( &p_vpar->bit_stream, 11 );
 908
 909     /* Test the presence of the escape character */
 910     while( i_index == 8 )
 911     {
 912         RemoveBits( &p_vpar->bit_stream, 11 );
 913         i_addr_inc += 33;
 914         i_index = ShowBits( &p_vpar->bit_stream, 11 );
 915     }
 916
 917     /* Affect the value from the lookup table */
 918     i_addr_inc += p_vpar->pl_mb_addr_inc[i_index].i_value;
 919
 920     /* Dump the good number of bits */
 921     RemoveBits( &p_vpar->bit_stream, p_vpar->pl_mb_addr_inc[i_index].i_length );
 922
 923     return i_addr_inc;
 924 }
 925
 926 /*****************************************************************************
 927  * MacroblockModes : Get the macroblock_modes structure
 928  *****************************************************************************/
 929 static __inline__ void MacroblockModes( vpar_thread_t * p_vpar,
 930                                         macroblock_t * p_mb )
 931 {
 932     static f_motion_t   pf_motion[2][4] =
 933         { {NULL, vdec_FieldFieldRecon, vdec_Field16x8Recon, vdec_FieldDMVRecon},
 934           {NULL, vdec_FrameFieldRecon, vdec_FrameFrameRecon, vdec_FrameDMVRecon} };
 935     static int          ppi_mv_count[2][4] = { {0, 1, 2, 1}, {0, 2, 1, 1} };
 936     static int          ppi_mv_format[2][4] = { {0, 1, 1, 1}, {0, 1, 2, 1} };
 937
 938     /* Get macroblock_type. */
 939     p_vpar->mb.i_mb_type = (p_vpar->picture.pf_macroblock_type)( p_vpar );
 940     p_mb->i_mb_type = p_vpar->mb.i_mb_type;
 941
 942 //fprintf( stderr, "MB type : %d\n", p_mb->i_mb_type );
 943
 944     /* SCALABILITY : warning, we don't know if spatial_temporal_weight_code
 945      * has to be dropped, take care if you use scalable streams. */
 946     /* RemoveBits( &p_vpar->bit_stream, 2 ); */
 947
 948     if( !(p_vpar->mb.i_mb_type & (MB_MOTION_FORWARD | MB_MOTION_BACKWARD)) )
 949     {
 950         /* If mb_type has neither MOTION_FORWARD nor MOTION_BACKWARD, this
 951          * is useless, but also harmless. */
 952         p_vpar->mb.i_motion_type = MOTION_FRAME;
 953     }
 954     else
 955     {
 956         if( p_vpar->picture.i_structure == FRAME_STRUCTURE
 957             && p_vpar->picture.b_frame_pred_frame_dct )
 958         {
 959             p_vpar->mb.i_motion_type = MOTION_FRAME;
 960         }
 961         else
 962         {
 963             p_vpar->mb.i_motion_type = GetBits( &p_vpar->bit_stream, 2 );
 964         }
 965     }
 966
 967     if( p_mb->b_P_coding_type && !(p_vpar->mb.i_mb_type & (MB_MOTION_FORWARD|MB_INTRA)) )
 968     {
 969         /* Special No-MC macroblock in P pictures (7.6.3.5). */
 970         memset( p_vpar->slice.pppi_pmv, 0, 8*sizeof(int) );
 971         memset( p_mb->pppi_motion_vectors, 0, 8*sizeof(int) );
 972
 973         p_vpar->mb.i_motion_type = MOTION_FRAME;
 974         p_mb->ppi_field_select[0][0] = ( p_vpar->picture.i_current_structure == BOTTOM_FIELD );
 975     }
 976
 977      if( p_vpar->mb.i_mb_type & MB_INTRA )
 978     {
 979         /* For the intra macroblocks, we use an empty motion
 980          * compensation function */
 981         p_mb->pf_motion = vdec_DummyRecon;
 982     }
 983     else
 984     {
 985         p_mb->pf_motion = pf_motion[p_vpar->picture.b_frame_structure]
 986                                    [p_vpar->mb.i_motion_type];
 987     }
 988     p_vpar->mb.i_mv_count = ppi_mv_count[p_vpar->picture.b_frame_structure]
 989                                         [p_vpar->mb.i_motion_type];
 990     p_vpar->mb.i_mv_format = ppi_mv_format[p_vpar->picture.b_frame_structure]
 991                                           [p_vpar->mb.i_motion_type];
 992
 993     p_vpar->mb.b_dct_type = 0;
 994     if( (p_vpar->picture.i_structure == FRAME_STRUCTURE) &&
 995         (!p_vpar->picture.b_frame_pred_frame_dct) &&
 996         (p_vpar->mb.i_mb_type & (MB_PATTERN|MB_INTRA)) )
 997     {
 998         if( (p_vpar->mb.b_dct_type = GetBits( &p_vpar->bit_stream, 1 )) )
 999         {
1000             /* The DCT is coded on fields. Jump one line between each
1001              * sample. */
1002             p_mb->i_addb_l_stride <<= 1;
1003                 p_mb->i_addb_l_stride += 8;
1004             /* With CHROMA_420, the DCT is necessarily frame-coded. */
1005             if( p_vpar->sequence.i_chroma_format != CHROMA_420 )
1006             {
1007                     p_mb->i_addb_c_stride <<= 1;
1008                     p_mb->i_addb_c_stride += 8;
1009             }
1010         }
1011     }
1012     p_vpar->mb.b_dmv = p_vpar->mb.i_motion_type == MOTION_DMV;
1013 }
1014
1015 /*****************************************************************************
1016  * vpar_IMBType : macroblock_type in I pictures
1017  *****************************************************************************/
1018 int vpar_IMBType( vpar_thread_t * p_vpar )
1019 {
1020     /* Take two bits for testing */
1021     int                 i_type = ShowBits( &p_vpar->bit_stream, 2 );
1022
1023     /* Lookup table for macroblock_type */
1024     static lookup_t     pl_mb_Itype[4] = { {MB_ERROR, 0},
1025                                            {MB_QUANT|MB_INTRA, 2},
1026                                            {MB_INTRA, 1},
1027                                            {MB_INTRA, 1} };
1028     /* Dump the good number of bits */
1029     RemoveBits( &p_vpar->bit_stream, pl_mb_Itype[i_type].i_length );
1030     return pl_mb_Itype[i_type].i_value;
1031 }
1032
1033 /*****************************************************************************
1034  * vpar_PMBType : macroblock_type in P pictures
1035  *****************************************************************************/
1036 int vpar_PMBType( vpar_thread_t * p_vpar )
1037 {
1038     /* Testing on 6 bits */
1039     int                i_type = ShowBits( &p_vpar->bit_stream, 6 );
1040
1041     /* Dump the good number of bits */
1042     RemoveBits( &p_vpar->bit_stream, p_vpar->ppl_mb_type[0][i_type].i_length );
1043
1044     /* return the value from the lookup table for P type */
1045     return p_vpar->ppl_mb_type[0][i_type].i_value;
1046 }
1047
1048 /*****************************************************************************
1049  * vpar_BMBType : macroblock_type in B pictures
1050  *****************************************************************************/
1051 int vpar_BMBType( vpar_thread_t * p_vpar )
1052 {
1053      /* Testing on 6 bits */
1054     int                i_type = ShowBits( &p_vpar->bit_stream, 6 );
1055
1056     /* Dump the good number of bits */
1057     RemoveBits( &p_vpar->bit_stream, p_vpar->ppl_mb_type[1][i_type].i_length );
1058
1059     /* return the value from the lookup table for B type */
1060     return p_vpar->ppl_mb_type[1][i_type].i_value;
1061 }
1062
1063 /*****************************************************************************
1064  * vpar_DMBType : macroblock_type in D pictures
1065  *****************************************************************************/
1066 int vpar_DMBType( vpar_thread_t * p_vpar )
1067 {
1068     /* Taking 1 bit */
1069     int               i_type = GetBits( &p_vpar->bit_stream, 1 );
1070
1071     /* Lookup table */
1072     static int        pi_mb_Dtype[2] = { MB_ERROR, 1 };
1073
1074     return pi_mb_Dtype[i_type];
1075 }
1076
1077 /*****************************************************************************
1078  * vpar_CodedPattern420 : coded_block_pattern with 420 chroma
1079  *****************************************************************************/
1080 int vpar_CodedPattern420( vpar_thread_t * p_vpar )
1081 {
1082     /* Take the max 9 bits length vlc code for testing */
1083     int      i_vlc = ShowBits( &p_vpar->bit_stream, 9 );
1084
1085     /* Trash the good number of bits read in  the lookup table */
1086     RemoveBits( &p_vpar->bit_stream, p_vpar->pl_coded_pattern[i_vlc].i_length );
1087
1088     /* return the value from the vlc table */
1089     return p_vpar->pl_coded_pattern[i_vlc].i_value;
1090 }
1091
1092 /*****************************************************************************
1093  * vpar_CodedPattern422 : coded_block_pattern with 422 chroma
1094  *****************************************************************************/
1095 int vpar_CodedPattern422( vpar_thread_t * p_vpar )
1096 {
1097     int      i_vlc = ShowBits( &p_vpar->bit_stream, 9 );
1098
1099     /* Supplementary 2 bits long code for 422 format */
1100     int      i_coded_block_pattern_1;
1101
1102     RemoveBits( &p_vpar->bit_stream, p_vpar->pl_coded_pattern[i_vlc].i_length );
1103     i_coded_block_pattern_1 = GetBits( &p_vpar->bit_stream, 2 );
1104
1105     /* the code is just to be added to the value found in the table */
1106     return p_vpar->pl_coded_pattern[i_vlc].i_value |
1107            (i_coded_block_pattern_1 << 6);
1108 }
1109
1110 /*****************************************************************************
1111  * vpar_CodedPattern444 : coded_block_pattern with 444 chroma
1112  *****************************************************************************/
1113 int vpar_CodedPattern444( vpar_thread_t * p_vpar )
1114 {
1115     int      i_vlc = ShowBits( &p_vpar->bit_stream, 9 );
1116     int      i_coded_block_pattern_2;
1117
1118     RemoveBits( &p_vpar->bit_stream, p_vpar->pl_coded_pattern[i_vlc].i_length );
1119     i_coded_block_pattern_2 = GetBits( &p_vpar->bit_stream, 6 );
1120
1121     return p_vpar->pl_coded_pattern[i_vlc].i_value |
1122            ( i_coded_block_pattern_2 << 6 );
1123 }
1124
1125 /*****************************************************************************
1126  * vpar_DecodeMPEG1Non : decode MPEG-1 non-intra blocks
1127  *****************************************************************************/
1128 static void vpar_DecodeMPEG1Non( vpar_thread_t * p_vpar, macroblock_t * p_mb, int i_b )
1129 {
1130
1131     if( p_vpar->picture.i_coding_type == D_CODING_TYPE )
1132     {
1133         /* Remove end_of_macroblock (always 1, prevents startcode emulation)
1134          * ISO/IEC 11172-2 section 2.4.2.7 and 2.4.3.6 */
1135         RemoveBits( &p_vpar->bit_stream, 1 );
1136     }
1137 }
1138
1139 /*****************************************************************************
1140  * vpar_DecodeMPEG1Intra : decode MPEG-1 intra blocks
1141  *****************************************************************************/
1142 static void vpar_DecodeMPEG1Intra( vpar_thread_t * p_vpar, macroblock_t * p_mb, int i_b )
1143 {
1144
1145     if( p_vpar->picture.i_coding_type == D_CODING_TYPE )
1146     {
1147         /* Remove end_of_macroblock (always 1, prevents startcode emulation)
1148          * ISO/IEC 11172-2 section 2.4.2.7 and 2.4.3.6 */
1149         RemoveBits( &p_vpar->bit_stream, 1 );
1150     }
1151 }
1152
1153 /*****************************************************************************
1154  * vpar_DecodeMPEG2Non : decode MPEG-2 non-intra blocks
1155  *****************************************************************************/
1156 static void vpar_DecodeMPEG2Non( vpar_thread_t * p_vpar, macroblock_t * p_mb, int i_b )
1157 {
1158     int         i_parse;
1159     int         i_nc;
1160     int         i_cc;
1161     int         i_coef;
1162     int         i_type;
1163     int         i_code;
1164     int         i_length;
1165     int         i_pos;
1166     int         i_run;
1167     int         i_level;
1168     boolean_t   b_sign;
1169     int *       ppi_quant[2];
1170
1171     /* Lookup Table for the chromatic component */
1172     static int pi_cc_index[12] = { 0, 0, 0, 0, 1, 2, 1, 2, 1, 2 };
1173
1174     i_cc = pi_cc_index[i_b];
1175
1176     /* Determine whether it is luminance or not (chrominance) */
1177     i_type = ( i_cc + 1 ) >> 1;
1178
1179     /* Give a pointer to the quantization matrices for intra blocks */
1180     ppi_quant[0] = p_vpar->sequence.nonintra_quant.pi_matrix;
1181     ppi_quant[1] = p_vpar->sequence.chroma_nonintra_quant.pi_matrix;
1182
1183     /* Decoding of the AC coefficients */
1184
1185     i_nc = 0;
1186     i_coef = 0;
1187     for( i_parse = 0; ; i_parse++ )
1188     {
1189         i_code = ShowBits( &p_vpar->bit_stream, 16 );
1190         if( i_code >= 16384 )
1191         {
1192             if( i_parse == 0 )
1193             {
1194                 i_run =     pl_DCT_tab_dc[(i_code>>12)-4].i_run;
1195                 i_level =   pl_DCT_tab_dc[(i_code>>12)-4].i_level;
1196                 i_length =  pl_DCT_tab_dc[(i_code>>12)-4].i_length;
1197             }
1198             else
1199             {
1200                 i_run =     pl_DCT_tab_ac[(i_code>>12)-4].i_run;
1201                 i_level =   pl_DCT_tab_ac[(i_code>>12)-4].i_level;
1202                 i_length =  pl_DCT_tab_ac[(i_code>>12)-4].i_length;
1203              }
1204         }
1205         else
1206         {
1207             i_run =     ppl_dct_coef[0][i_code].i_run;
1208             i_length =  ppl_dct_coef[0][i_code].i_length;
1209             i_level =   ppl_dct_coef[0][i_code].i_level;
1210         }
1211
1212
1213         RemoveBits( &p_vpar->bit_stream, i_length );
1214
1215         switch( i_run )
1216         {
1217             case DCT_ESCAPE:
1218                 i_run = GetBits( &p_vpar->bit_stream, 6 );
1219                 i_level = GetBits( &p_vpar->bit_stream, 12 );
1220                 i_level = (b_sign = ( i_level > 2047 )) ? 4096 - i_level
1221                                                         : i_level;
1222                 break;
1223             case DCT_EOB:
1224                 if( i_nc <= 1 )
1225                 {
1226                     p_mb->pf_idct[i_b] = vdec_SparseIDCT;
1227                     p_mb->pi_sparse_pos[i_b] = i_coef;
1228                 }
1229                 else
1230                 {
1231                     p_mb->pf_idct[i_b] = vdec_IDCT;
1232                 }
1233                 return;
1234
1235                 break;
1236             default:
1237                 b_sign = GetBits( &p_vpar->bit_stream, 1 );
1238         }
1239         i_coef = i_parse;
1240         i_parse += i_run;
1241         i_nc ++;
1242
1243         i_pos = pi_scan[p_vpar->picture.b_alternate_scan][i_parse];
1244         i_level = ( i_level *
1245                     p_vpar->slice.i_quantizer_scale *
1246                     ppi_quant[i_type][i_pos] ) >> 5;
1247         p_mb->ppi_blocks[i_b][i_parse] = b_sign ? -i_level : i_level;
1248     }
1249 fprintf( stderr, "MPEG2 end (%d)\n", i_b );
1250 exit(0);
1251
1252 }
1253
1254 /*****************************************************************************
1255  * vpar_DecodeMPEG2Intra : decode MPEG-2 intra blocks
1256  *****************************************************************************/
1257 static void vpar_DecodeMPEG2Intra( vpar_thread_t * p_vpar, macroblock_t * p_mb, int i_b )
1258 {
1259     int         i_parse;
1260     int         i_nc;
1261     int         i_cc;
1262     int         i_coef;
1263     int         i_type;
1264     int         i_code;
1265     int         i_length;
1266     int         i_pos;
1267     int         i_dct_dc_size;
1268     int         i_dct_dc_diff;
1269     int         i_run;
1270     int         i_level;
1271     boolean_t   b_vlc_intra;
1272     boolean_t   b_sign;
1273     int *       ppi_quant[2];
1274
1275     /* Lookup Table for the chromatic component */
1276     static int pi_cc_index[12] = { 0, 0, 0, 0, 1, 2, 1, 2, 1, 2 };
1277
1278     i_cc = pi_cc_index[i_b];
1279
1280     /* Determine whether it is luminance or not (chrominance) */
1281     i_type = ( i_cc + 1 ) >> 1;
1282
1283     /* Give a pointer to the quantization matrices for intra blocks */
1284     ppi_quant[0] = p_vpar->sequence.intra_quant.pi_matrix;
1285     ppi_quant[1] = p_vpar->sequence.chroma_intra_quant.pi_matrix;
1286
1287 #if 0
1288     /* Decoding of the DC intra coefficient */
1289     /* The nb of bits to parse depends on i_type */
1290     i_code = ShowBits( &p_vpar->bit_stream, 9 + i_type );
1291
1292     /* To reduce memory occupation, there are two lookup tables
1293      * See InitDCT above */
1294     i_code5 = i_code >> (4+i_type);
1295
1296     /* Shall we lookup in the first or in the second table ? */
1297     i_select = ( i_code5 == 31 );
1298     /* Offset value for looking in the second table */
1299     i_offset = 0x1f0 + ( i_type * 0x1f0 );
1300     i_pos = ( i_code5 * ( ! i_select ) ) + ( ( i_code - i_offset ) * i_select );
1301     i_dct_dc_size = p_vpar->pppl_dct_dc_size[i_type][i_select][i_pos].i_value;
1302 #endif
1303
1304     if( !i_type/*i_b < 4*/ )
1305     {
1306         /* decode length */
1307         i_code = ShowBits(&p_vpar->bit_stream, 5);
1308         if (i_code<31)
1309         {
1310             i_dct_dc_size = pl_dct_dc_lum_init_table_1[i_code].i_value;
1311             i_length = pl_dct_dc_lum_init_table_1[i_code].i_length;
1312             RemoveBits( &p_vpar->bit_stream, i_length);
1313         }
1314         else
1315         {
1316             i_code = ShowBits(&p_vpar->bit_stream, 9) - 0x1f0;
1317             i_dct_dc_size = pl_dct_dc_lum_init_table_2[i_code].i_value;
1318             i_length = pl_dct_dc_lum_init_table_2[i_code].i_length;
1319             RemoveBits( &p_vpar->bit_stream, i_length);
1320         }
1321     }
1322     else
1323     {
1324         /* decode length */
1325         i_code = ShowBits(&p_vpar->bit_stream, 5);
1326
1327         if (i_code<31)
1328         {
1329             i_dct_dc_size = pl_dct_dc_chrom_init_table_1[i_code].i_value;
1330             i_length = pl_dct_dc_chrom_init_table_1[i_code].i_length;
1331             RemoveBits(&p_vpar->bit_stream, i_length);
1332         }
1333         else
1334         {
1335             i_code = ShowBits(&p_vpar->bit_stream, 10) - 0x3e0;
1336             i_dct_dc_size = pl_dct_dc_chrom_init_table_2[i_code].i_value;
1337             i_length = pl_dct_dc_chrom_init_table_2[i_code].i_length;
1338             RemoveBits( &p_vpar->bit_stream, i_length);
1339         }
1340     }
1341     if (i_dct_dc_size==0)
1342         i_dct_dc_diff = 0;
1343     else
1344     {
1345         i_dct_dc_diff = GetBits( &p_vpar->bit_stream, i_dct_dc_size);
1346         if ((i_dct_dc_diff & (1<<(i_dct_dc_size-1)))==0)
1347             i_dct_dc_diff-= (1<<i_dct_dc_size) - 1;
1348     }
1349
1350     /* Dump the variable length code */
1351     //RemoveBits( &p_vpar->bit_stream,
1352     //          p_vpar->pppl_dct_dc_size[i_type][i_select][i_pos].i_length );
1353
1354     /* Read the actual code with the good length */
1355     p_vpar->slice.pi_dc_dct_pred[i_cc] += i_dct_dc_diff;
1356
1357     p_mb->ppi_blocks[i_b][0] = ( p_vpar->slice.pi_dc_dct_pred[i_cc] <<
1358                                ( 3 - p_vpar->picture.i_intra_dc_precision ) );
1359     i_nc = ( p_vpar->slice.pi_dc_dct_pred[i_cc] != 0 );
1360
1361 //fprintf( stderr, "coucou\n" );
1362     /* Decoding of the AC coefficients */
1363
1364     i_coef = 0;
1365     b_vlc_intra = p_vpar->picture.b_intra_vlc_format;
1366     for( i_parse = 1; /*i_parse < 64*/; i_parse++ )
1367     {
1368         i_code = ShowBits( &p_vpar->bit_stream, 16 );
1369         if( i_code >= 16384 )
1370         {
1371             if( b_vlc_intra )
1372             {
1373                 i_run =     pl_DCT_tab0a[(i_code>>8)-4].i_run;
1374                 i_level =   pl_DCT_tab0a[(i_code>>8)-4].i_level;
1375                 i_length =  pl_DCT_tab0a[(i_code>>8)-4].i_length;
1376 //fprintf( stderr, "**********> %d, %d, %d *******\n", i_run, i_level, (i_code>>8)-4 );
1377             }
1378             else
1379             {
1380                 i_run =     pl_DCT_tab_ac[(i_code>>12)-4].i_run;
1381                 i_level =   pl_DCT_tab_ac[(i_code>>12)-4].i_level;
1382                 i_length =  pl_DCT_tab_ac[(i_code>>12)-4].i_length;
1383              }
1384         }
1385         else
1386         {
1387             i_run =     ppl_dct_coef[b_vlc_intra][i_code].i_run;
1388             i_length =  ppl_dct_coef[b_vlc_intra][i_code].i_length;
1389             i_level =   ppl_dct_coef[b_vlc_intra][i_code].i_level;
1390         }
1391
1392 #if 0
1393         {
1394             int code = i_code;
1395             int intra_vlc_format = b_vlc_intra;
1396             dct_lookup_t tab;
1397
1398     if (code>=16384 && !intra_vlc_format)
1399       tab = pl_DCT_tab_ac[(code>>12)-4];
1400     else if (code>=1024)
1401     {
1402       if (intra_vlc_format)
1403         tab = pl_DCT_tab0a[(code>>8)-4];
1404       else
1405         tab = pl_DCT_tab0[(code>>8)-4];
1406     }
1407     else if (code>=512)
1408     {
1409       if (intra_vlc_format)
1410         tab = pl_DCT_tab1a[(code>>6)-8];
1411       else
1412         tab = pl_DCT_tab1[(code>>6)-8];
1413     }
1414     else if (code>=256)
1415       tab = pl_DCT_tab2[(code>>4)-16];
1416     else if (code>=128)
1417       tab = pl_DCT_tab3[(code>>3)-16];
1418     else if (code>=64)
1419       tab = pl_DCT_tab4[(code>>2)-16];
1420     else if (code>=32)
1421       tab = pl_DCT_tab5[(code>>1)-16];
1422     else if (code>=16)
1423       tab = pl_DCT_tab6[code-16];
1424     else
1425     {
1426        fprintf( stderr, "invalid Huffman code in Decode_MPEG2_Intra_Block()\n");
1427     }
1428
1429     if( (i_run != tab.i_run) || (i_length != tab.i_length) || (i_level != tab.i_level) )
1430     {
1431         fprintf( stderr, "ET M....... !!!\n" );
1432         exit(0);
1433     }
1434         }
1435 #endif
1436
1437
1438
1439
1440         RemoveBits( &p_vpar->bit_stream, i_length );
1441
1442         switch( i_run )
1443         {
1444             case DCT_ESCAPE:
1445                 i_run = GetBits( &p_vpar->bit_stream, 6 );
1446                 i_level = GetBits( &p_vpar->bit_stream, 12 );
1447                 /*p_mb->ppi_blocks[i_b][i_parse] = ( b_sign = ( i_level > 2047 ) )
1448                                                           ? ( -4096 + i_level )
1449                                                           : i_level;*/
1450                 i_level = (b_sign = ( i_level > 2047 )) ? 4096 - i_level
1451                                                         : i_level;
1452                 break;
1453             case DCT_EOB:
1454                 if( i_nc <= 1 )
1455                 {
1456                     p_mb->pf_idct[i_b] = vdec_SparseIDCT;
1457                     p_mb->pi_sparse_pos[i_b] = i_coef;
1458                 }
1459                 else
1460                 {
1461                     p_mb->pf_idct[i_b] = vdec_IDCT;
1462                 }
1463                 return;
1464
1465                 break;
1466             default:
1467                 b_sign = GetBits( &p_vpar->bit_stream, 1 );
1468                 //p_mb->ppi_blocks[i_b][i_parse] = b_sign ? -i_level : i_level;
1469         }
1470 // fprintf( stderr, "i_code : %d (%d), run : %d, %d, %d (%4x) ", i_code ,  b_vlc_intra,
1471 //                  i_run, i_level, i_parse, ShowBits( &p_vpar->bit_stream, 16 ) );
1472
1473 //fprintf( stderr, "- %4x\n",ShowBits( &p_vpar->bit_stream, 16 ) );
1474
1475         i_coef = i_parse;
1476         i_parse += i_run;
1477         i_nc ++;
1478
1479         i_pos = pi_scan[p_vpar->picture.b_alternate_scan][i_parse];
1480         i_level = ( i_level *
1481                     p_vpar->slice.i_quantizer_scale *
1482                     ppi_quant[i_type][i_pos] ) >> 4;
1483         p_mb->ppi_blocks[i_b][i_parse] = b_sign ? -i_level : i_level;
1484     }
1485
1486 fprintf( stderr, "MPEG2 end (%d)\n", i_b );
1487 exit(0);
1488 }