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