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