]> git.sesse.net Git - vlc/blob - src/video_parser/vpar_blocks.c
* Fonction de d�codage des Blocs MPEG2intra.
[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  * Local prototypes
44  */
45 static __inline__ void InitMacroblock( vpar_thread_t * p_vpar,
46                                        macroblock_t * p_mb );
47 static __inline__ int MacroblockAddressIncrement( vpar_thread_t * p_vpar );
48 static __inline__ void MacroblockModes( vpar_thread_t * p_vpar,
49                                         macroblock_t * p_mb );
50 typedef void (*f_decode_block_t)( vpar_thread_t *, macroblock_t *, int );
51 static void vpar_DecodeMPEG1Non( vpar_thread_t * p_vpar, macroblock_t * p_mb, int i_b );
52 static void vpar_DecodeMPEG1Intra( vpar_thread_t * p_vpar, macroblock_t * p_mb, int i_b );
53 static void vpar_DecodeMPEG2Non( vpar_thread_t * p_vpar, macroblock_t * p_mb, int i_b );
54 static void vpar_DecodeMPEG2Intra( vpar_thread_t * p_vpar, macroblock_t * p_mb, int i_b );
55
56 /*
57  * Initialisation tables
58  */
59     /* Table for coded_block_pattern resolution */
60 lookup_t     pl_coded_pattern_init_table[512] = 
61     { {MB_ERROR, 0}, {0, 9}, {39, 9}, {27, 9}, {59, 9}, {55, 9}, {47, 9}, {31, 9},
62     {58, 8}, {58, 8}, {54, 8}, {54, 8}, {46, 8}, {46, 8}, {30, 8}, {30, 8},
63     {57, 8}, {57, 8}, {53, 8}, {53, 8}, {45, 8}, {45, 8}, {29, 8}, {29, 8},
64     {38, 8}, {38, 8}, {26, 8}, {26, 8}, {37, 8}, {37, 8}, {25, 8}, {25, 8},
65     {43, 8}, {43, 8}, {23, 8}, {23, 8}, {51, 8}, {51, 8}, {15, 8}, {15, 8},
66     {42, 8}, {42, 8}, {22, 8}, {22, 8}, {50, 8}, {50, 8}, {14, 8}, {14, 8},
67     {41, 8}, {41, 8}, {21, 8}, {21, 8}, {49, 8}, {49, 8}, {13, 8}, {13, 8},
68     {35, 8}, {35, 8}, {19, 8}, {19, 8}, {11, 8}, {11, 8}, {7, 8}, {7, 8},
69     {34, 7}, {34, 7}, {34, 7}, {34, 7}, {18, 7}, {18, 7}, {18, 7}, {18, 7},
70     {10, 7}, {10, 7}, {10, 7}, {10, 7}, {6, 7}, {6, 7}, {6, 7}, {6, 7}, 
71     {33, 7}, {33, 7}, {33, 7}, {33, 7}, {17, 7}, {17, 7}, {17, 7}, {17, 7}, 
72     {9, 7}, {9, 7}, {9, 7}, {9, 7}, {5, 7}, {5, 7}, {5, 7}, {5, 7}, 
73     {63, 6}, {63, 6}, {63, 6}, {63, 6}, {63, 6}, {63, 6}, {63, 6}, {63, 6}, 
74     {3, 6}, {3, 6}, {3, 6}, {3, 6}, {3, 6}, {3, 6}, {3, 6}, {3, 6}, 
75     {36, 6}, {36, 6}, {36, 6}, {36, 6}, {36, 6}, {36, 6}, {36, 6}, {36, 6}, 
76     {24, 6}, {24, 6}, {24, 6}, {24, 6}, {24, 6}, {24, 6}, {24, 6}, {24, 6}, 
77     {62, 5}, {62, 5}, {62, 5}, {62, 5}, {62, 5}, {62, 5}, {62, 5}, {62, 5},
78     {62, 5}, {62, 5}, {62, 5}, {62, 5}, {62, 5}, {62, 5}, {62, 5}, {62, 5},
79     {2, 5}, {2, 5}, {2, 5}, {2, 5}, {2, 5}, {2, 5}, {2, 5}, {2, 5}, 
80     {2, 5}, {2, 5}, {2, 5}, {2, 5}, {2, 5}, {2, 5}, {2, 5}, {2, 5}, 
81     {61, 5}, {61, 5}, {61, 5}, {61, 5}, {61, 5}, {61, 5}, {61, 5}, {61, 5}, 
82     {61, 5}, {61, 5}, {61, 5}, {61, 5}, {61, 5}, {61, 5}, {61, 5}, {61, 5}, 
83     {1, 5}, {1, 5}, {1, 5}, {1, 5}, {1, 5}, {1, 5}, {1, 5}, {1, 5}, 
84     {1, 5}, {1, 5}, {1, 5}, {1, 5}, {1, 5}, {1, 5}, {1, 5}, {1, 5}, 
85     {56, 5}, {56, 5}, {56, 5}, {56, 5}, {56, 5}, {56, 5}, {56, 5}, {56, 5}, 
86     {56, 5}, {56, 5}, {56, 5}, {56, 5}, {56, 5}, {56, 5}, {56, 5}, {56, 5}, 
87     {52, 5}, {52, 5}, {52, 5}, {52, 5}, {52, 5}, {52, 5}, {52, 5}, {52, 5}, 
88     {52, 5}, {52, 5}, {52, 5}, {52, 5}, {52, 5}, {52, 5}, {52, 5}, {52, 5}, 
89     {44, 5}, {44, 5}, {44, 5}, {44, 5}, {44, 5}, {44, 5}, {44, 5}, {44, 5}, 
90     {44, 5}, {44, 5}, {44, 5}, {44, 5}, {44, 5}, {44, 5}, {44, 5}, {44, 5}, 
91     {28, 5}, {28, 5}, {28, 5}, {28, 5}, {28, 5}, {28, 5}, {28, 5}, {28, 5}, 
92     {28, 5}, {28, 5}, {28, 5}, {28, 5}, {28, 5}, {28, 5}, {28, 5}, {28, 5}, 
93     {40, 5}, {40, 5}, {40, 5}, {40, 5}, {40, 5}, {40, 5}, {40, 5}, {40, 5}, 
94     {40, 5}, {40, 5}, {40, 5}, {40, 5}, {40, 5}, {40, 5}, {40, 5}, {40, 5}, 
95     {20, 5}, {20, 5}, {20, 5}, {20, 5}, {20, 5}, {20, 5}, {20, 5}, {20, 5}, 
96     {20, 5}, {20, 5}, {20, 5}, {20, 5}, {20, 5}, {20, 5}, {20, 5}, {20, 5}, 
97     {48, 5}, {48, 5}, {48, 5}, {48, 5}, {48, 5}, {48, 5}, {48, 5}, {48, 5}, 
98     {48, 5}, {48, 5}, {48, 5}, {48, 5}, {48, 5}, {48, 5}, {48, 5}, {48, 5}, 
99     {12, 5}, {12, 5}, {12, 5}, {12, 5}, {12, 5}, {12, 5}, {12, 5}, {12, 5}, 
100     {12, 5}, {12, 5}, {12, 5}, {12, 5}, {12, 5}, {12, 5}, {12, 5}, {12, 5}, 
101     {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4}, 
102     {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4}, 
103     {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4}, 
104     {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4}, {32, 4}, 
105     {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4}, 
106     {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4}, 
107     {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4}, 
108     {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4}, {16, 4}, 
109     {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4}, 
110     {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4}, 
111     {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4}, 
112     {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4}, {8, 4},
113     {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4},
114     {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4}, 
115     {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4}, 
116     {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4}, {4, 4},
117     {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, 
118     {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, 
119     {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, 
120     {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, 
121     {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, 
122     {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, 
123     {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, 
124     {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3}, {60, 3} };
125
126     /* Tables for dc DCT coefficients 
127      * Tables are cut in two parts to reduce memory occupation
128      */
129     
130     /* Table B-12, dct_dc_size_luminance, codes 00xxx ... 11110 */
131 lookup_t     pl_dct_dc_lum_init_table_1[32] =
132     { {1, 2}, {1, 2}, {1, 2}, {1, 2}, {1, 2}, {1, 2}, {1, 2}, {1, 2},
133       {2, 2}, {2, 2}, {2, 2}, {2, 2}, {2, 2}, {2, 2}, {2, 2}, {2, 2},
134       {0, 3}, {0, 3}, {0, 3}, {0, 3}, {3, 3}, {3, 3}, {3, 3}, {3, 3},
135       {4, 3}, {4, 3}, {4, 3}, {4, 3}, {5, 4}, {5, 4}, {6, 5}, {MB_ERROR, 0}
136     };
137
138     /* Table B-12, dct_dc_size_luminance, codes 111110xxx ... 111111111 */
139 lookup_t     pl_dct_dc_lum_init_table_2[32] =
140     { {7, 6}, {7, 6}, {7, 6}, {7, 6}, {7, 6}, {7, 6}, {7, 6}, {7, 6},
141       {8, 7}, {8, 7}, {8, 7}, {8, 7}, {9, 8}, {9, 8}, {10,9}, {11,9},
142       {MB_ERROR, 0}, {MB_ERROR, 0}, {MB_ERROR, 0}, {MB_ERROR, 0},
143       {MB_ERROR, 0}, {MB_ERROR, 0}, {MB_ERROR, 0}, {MB_ERROR, 0},
144       {MB_ERROR, 0}, {MB_ERROR, 0}, {MB_ERROR, 0}, {MB_ERROR, 0},
145       {MB_ERROR, 0}, {MB_ERROR, 0}, {MB_ERROR, 0}, {MB_ERROR, 0}
146     };
147
148     /* Table B-13, dct_dc_size_chrominance, codes 00xxx ... 11110 */
149 lookup_t     pl_dct_dc_chrom_init_table_1[32] =
150     { {0, 2}, {0, 2}, {0, 2}, {0, 2}, {0, 2}, {0, 2}, {0, 2}, {0, 2},
151       {1, 2}, {1, 2}, {1, 2}, {1, 2}, {1, 2}, {1, 2}, {1, 2}, {1, 2},
152       {2, 2}, {2, 2}, {2, 2}, {2, 2}, {2, 2}, {2, 2}, {2, 2}, {2, 2},
153       {3, 3}, {3, 3}, {3, 3}, {3, 3}, {4, 4}, {4, 4}, {5, 5}, {MB_ERROR, 0}
154     };
155
156     /* Table B-13, dct_dc_size_chrominance, codes 111110xxxx ... 1111111111 */
157 lookup_t     pl_dct_dc_chrom_init_table_2[32] =
158     { {6, 6}, {6, 6}, {6, 6}, {6, 6}, {6, 6}, {6, 6}, {6, 6}, {6, 6},
159       {6, 6}, {6, 6}, {6, 6}, {6, 6}, {6, 6}, {6, 6}, {6, 6}, {6, 6},
160       {7, 7}, {7, 7}, {7, 7}, {7, 7}, {7, 7}, {7, 7}, {7, 7}, {7, 7},
161       {8, 8}, {8, 8}, {8, 8}, {8, 8}, {9, 9}, {9, 9}, {10,10}, {11,10}
162     };
163
164 /* 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 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 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 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 dct_lookup_t pl_DCT_tab0a[252] =
212     {
213         {DCT_ESCAPE,0,6}, {DCT_ESCAPE,0,6}, {DCT_ESCAPE,0,6}, {DCT_ESCAPE,0,6},
214         /* Escape */
215         {7,1,7}, {7,1,7}, {8,1,7}, {8,1,7},
216         {6,1,7}, {6,1,7}, {2,2,7}, {2,2,7},
217         {0,7,6}, {0,7,6}, {0,7,6}, {0,7,6},
218         {0,6,6}, {0,6,6}, {0,6,6}, {0,6,6},
219         {4,1,6}, {4,1,6}, {4,1,6}, {4,1,6},
220         {5,1,6}, {5,1,6}, {5,1,6}, {5,1,6},
221         {1,5,8}, {11,1,8}, {0,11,8}, {0,10,8},
222         {13,1,8}, {12,1,8}, {3,2,8}, {1,4,8},
223         {2,1,5}, {2,1,5}, {2,1,5}, {2,1,5},
224         {2,1,5}, {2,1,5}, {2,1,5}, {2,1,5},
225         {1,2,5}, {1,2,5}, {1,2,5}, {1,2,5},
226         {1,2,5}, {1,2,5}, {1,2,5}, {1,2,5},
227         {3,1,5}, {3,1,5}, {3,1,5}, {3,1,5},
228         {3,1,5}, {3,1,5}, {3,1,5}, {3,1,5},
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         {1,1,3}, {1,1,3}, {1,1,3}, {1,1,3},
237         {DCT_EOB,0,4}, {DCT_EOB,0,4}, {DCT_EOB,0,4}, {DCT_EOB,0,4}, /* EOB */
238         {DCT_EOB,0,4}, {DCT_EOB,0,4}, {DCT_EOB,0,4}, {DCT_EOB,0,4},
239         {DCT_EOB,0,4}, {DCT_EOB,0,4}, {DCT_EOB,0,4}, {DCT_EOB,0,4},
240         {DCT_EOB,0,4}, {DCT_EOB,0,4}, {DCT_EOB,0,4}, {DCT_EOB,0,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,3,4}, {0,3,4}, {0,3,4}, {0,3,4},
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 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 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 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 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 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 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 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 + 32768;
364
365     for( i_dummy = -32768; 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 < 32767; 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 /* Fonction for filling up the lookup table for mb_addr_inc */
385 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_dummy, i_dummy2;
390     for( i_dummy = i_start ; i_dummy < i_end ; i_dummy += i_step )
391         for( i_dummy2 = 0 ; i_dummy2 < i_step ; i_dummy2 ++ )
392         {
393             p_vpar->pl_mb_addr_inc[i_dummy + i_dummy2].i_value = * pi_value;
394             p_vpar->pl_mb_addr_inc[i_dummy + i_dummy2].i_length = i_length;
395         }
396     (*pi_value)--;
397 }
398     
399 /* Fonction that initialize the table using the last one */
400 void InitMbAddrInc( vpar_thread_t * p_vpar )
401 {
402     int i_dummy;
403     int *  pi_value;
404     for( i_dummy = 0 ; i_dummy < 8 ; i_dummy++ )
405     {
406         p_vpar->pl_mb_addr_inc[i_dummy].i_value = MB_ERROR;
407         p_vpar->pl_mb_addr_inc[i_dummy].i_length = 0;
408     }
409
410     p_vpar->pl_mb_addr_inc[8].i_value = MB_ADDRINC_ESCAPE;
411     p_vpar->pl_mb_addr_inc[8].i_length = 11;
412
413     for( i_dummy = 9 ; i_dummy < 15 ; i_dummy ++ ) 
414     {
415         p_vpar->pl_mb_addr_inc[i_dummy].i_value =  MB_ERROR;
416         p_vpar->pl_mb_addr_inc[i_dummy].i_length = 0;
417     }
418
419     p_vpar->pl_mb_addr_inc[15].i_value = MB_ADDRINC_STUFFING;
420     p_vpar->pl_mb_addr_inc[15].i_length = 11;
421
422     for( i_dummy = 16; i_dummy < 24; i_dummy++ ) 
423     {
424         p_vpar->pl_mb_addr_inc[i_dummy].i_value =  MB_ERROR;
425         p_vpar->pl_mb_addr_inc[i_dummy].i_length = 0;
426     }
427     pi_value = (int *) malloc( sizeof( int ) );
428     * pi_value = 33;
429     FillMbAddrIncTable( p_vpar, 1024, 2048, 1024, pi_value, 1 );
430     FillMbAddrIncTable( p_vpar, 512, 1024, 256, pi_value, 3 );
431     FillMbAddrIncTable( p_vpar, 256, 512, 128, pi_value, 4 );
432     FillMbAddrIncTable( p_vpar, 128, 256, 64, pi_value, 5 );
433     FillMbAddrIncTable( p_vpar, 96, 128, 16, pi_value, 7 );
434     FillMbAddrIncTable( p_vpar, 48, 96, 8, pi_value, 8 );
435     FillMbAddrIncTable( p_vpar, 36, 48, 2, pi_value, 10 );
436     FillMbAddrIncTable( p_vpar, 24, 36, 1, pi_value, 11 );
437 }
438
439 /*****************************************************************************
440  * Init*MBType : Initialize lookup table for the Macroblock type
441  *****************************************************************************/
442
443 /* Fonction for filling up the tables */
444 static void __inline__ FillMBType( vpar_thread_t * p_vpar,
445                                    int           i_mb_type,
446                                    int           i_start,
447                                    int           i_end,
448                                    int           i_value,
449                                    int           i_length )
450 {
451     int i_dummy;
452     
453     for( i_dummy = i_start ; i_dummy < i_end ; i_dummy++ )
454     {
455         p_vpar->ppl_mb_type[i_mb_type][i_dummy].i_value = i_value;
456         p_vpar->ppl_mb_type[i_mb_type][i_dummy].i_length = i_length;
457     }
458 }
459
460 /* Fonction that fills the table for P MB_Type */
461 void InitPMBType( vpar_thread_t * p_vpar )
462 {
463     FillMBType( p_vpar, 0, 32, 64, MB_MOTION_FORWARD|MB_PATTERN, 1 );
464     FillMBType( p_vpar, 0, 16, 32, MB_PATTERN, 2 );
465     FillMBType( p_vpar, 0, 8, 16, MB_MOTION_FORWARD, 3 );
466     FillMBType( p_vpar, 0, 6, 8, MB_INTRA, 5 );
467     FillMBType( p_vpar, 0, 4, 6, MB_QUANT|MB_MOTION_FORWARD|MB_PATTERN, 5 );
468     FillMBType( p_vpar, 0, 2, 4, MB_QUANT|MB_PATTERN, 5 );
469     p_vpar->ppl_mb_type[0][1].i_value = MB_QUANT|MB_INTRA;
470     p_vpar->ppl_mb_type[0][1].i_length = 6;
471     p_vpar->ppl_mb_type[0][0].i_value = MB_ERROR;
472     p_vpar->ppl_mb_type[0][0].i_length = 0;
473 }
474
475 /* Fonction that fills the table for B MB_Type */
476 void InitBMBType( vpar_thread_t * p_vpar )
477 {
478     FillMBType( p_vpar, 1, 48, 64, MB_MOTION_FORWARD
479                                   |MB_MOTION_BACKWARD|MB_PATTERN, 2 );
480     FillMBType( p_vpar, 1, 32, 48, MB_MOTION_FORWARD|MB_MOTION_BACKWARD, 2 );
481     FillMBType( p_vpar, 1, 24, 32, MB_MOTION_BACKWARD|MB_PATTERN, 3 );
482     FillMBType( p_vpar, 1, 16, 24, MB_MOTION_BACKWARD, 3 );
483     FillMBType( p_vpar, 1, 12, 16, MB_MOTION_FORWARD|MB_PATTERN, 4 );
484     FillMBType( p_vpar, 1, 8, 12, MB_MOTION_FORWARD, 4 );
485     FillMBType( p_vpar, 1, 6, 8, MB_INTRA, 5 );
486     FillMBType( p_vpar, 1, 4, 6, MB_QUANT|MB_MOTION_FORWARD
487                                 |MB_MOTION_BACKWARD|MB_PATTERN, 5 );
488     p_vpar->ppl_mb_type[1][3].i_value = MB_QUANT|MB_MOTION_FORWARD|MB_PATTERN;
489     p_vpar->ppl_mb_type[1][3].i_length = 6;
490     p_vpar->ppl_mb_type[1][2].i_value = MB_QUANT|MB_MOTION_BACKWARD|MB_PATTERN;
491     p_vpar->ppl_mb_type[1][2].i_length = 6;
492     p_vpar->ppl_mb_type[1][1].i_value = MB_QUANT|MB_INTRA;
493     p_vpar->ppl_mb_type[1][1].i_length = 6;
494     p_vpar->ppl_mb_type[1][0].i_value =MB_ERROR;
495     p_vpar->ppl_mb_type[1][0].i_length = 0;
496 }
497
498 /*****************************************************************************
499  * InitCodedPattern : Initialize the lookup table for decoding 
500  *                    coded block pattern
501  *****************************************************************************/
502 void InitCodedPattern( vpar_thread_t * p_vpar )
503 {
504    memcpy( p_vpar->pl_coded_pattern, pl_coded_pattern_init_table ,
505             sizeof(pl_coded_pattern_init_table) );
506 }
507
508 /*****************************************************************************
509  * InitDCT : Initialize tables giving the length of the dct coefficient
510  *           from the vlc code
511  *****************************************************************************/
512
513 void InitDCTTables( vpar_thread_t * p_vpar )
514 {
515     /* Tables are cut in two parts to reduce memory occupation */
516     memcpy( p_vpar->pppl_dct_dc_size[0][0], pl_dct_dc_lum_init_table_1,
517             sizeof( pl_dct_dc_lum_init_table_1 ) );
518     memcpy( p_vpar->pppl_dct_dc_size[0][1], pl_dct_dc_lum_init_table_2,
519             sizeof( pl_dct_dc_lum_init_table_2 ) );
520     memcpy( p_vpar->pppl_dct_dc_size[1][0], pl_dct_dc_chrom_init_table_1,
521             sizeof( pl_dct_dc_chrom_init_table_1 ) );
522     memcpy( p_vpar->pppl_dct_dc_size[1][1], pl_dct_dc_chrom_init_table_2,
523             sizeof( pl_dct_dc_chrom_init_table_2 ) );
524     /* For table B14 & B15, we have a pointer to tables */
525     memcpy( p_vpar->ppl_dct_coef[0], pl_DCT_tab_dc, sizeof( pl_DCT_tab_dc ) );
526     memcpy( p_vpar->ppl_dct_coef[1], pl_DCT_tab_ac, sizeof( pl_DCT_tab_ac ) );
527     memcpy( p_vpar->ppl_dct_coef[2], pl_DCT_tab0, sizeof( pl_DCT_tab0 ) );
528     memcpy( p_vpar->ppl_dct_coef[3], pl_DCT_tab1, sizeof( pl_DCT_tab1 ) );
529     memcpy( p_vpar->ppl_dct_coef[4], pl_DCT_tab2, sizeof( pl_DCT_tab2 ) );
530     memcpy( p_vpar->ppl_dct_coef[5], pl_DCT_tab3, sizeof( pl_DCT_tab3 ) );
531     memcpy( p_vpar->ppl_dct_coef[6], pl_DCT_tab4, sizeof( pl_DCT_tab4 ) );
532     memcpy( p_vpar->ppl_dct_coef[7], pl_DCT_tab5, sizeof( pl_DCT_tab5 ) );
533     memcpy( p_vpar->ppl_dct_coef[8], pl_DCT_tab6, sizeof( pl_DCT_tab6 ) );
534     memcpy( p_vpar->ppl_dct_coef[9], pl_DCT_tab0a, sizeof( pl_DCT_tab0a ) );
535     memcpy( p_vpar->ppl_dct_coef[10], pl_DCT_tab1a, sizeof( pl_DCT_tab1a ) );
536
537
538 }
539
540 /*
541  * Macroblock parsing functions
542  */
543
544 /*****************************************************************************
545  * vpar_ParseMacroblock : Parse the next macroblock
546  *****************************************************************************/
547 void vpar_ParseMacroblock( vpar_thread_t * p_vpar, int * pi_mb_address,
548                            int i_mb_previous, int i_mb_base )
549 {
550     static f_addb_t ppf_addb_intra[2] = {vdec_AddBlock, vdec_CopyBlock};
551     static f_decode_block_t pppf_decode_block[2][2] =
552                 { {vpar_DecodeMPEG1Non, vpar_DecodeMPEG1Intra},
553                   {vpar_DecodeMPEG2Non, vpar_DecodeMPEG2Intra} };
554     static int      pi_x[12] = {0,8,0,8,0,0,0,0,8,8,8,8};
555     static int      pi_y[2][12] = { {0,0,8,8,0,0,8,8,0,0,8,8},
556                                     {0,0,1,1,0,0,1,1,0,0,1,1} };
557     static int      pi_chroma_hor[4] = { 0, 1, 1, 0 };
558     static int      pi_chroma_ver[4] = { 0, 1, 0, 0 };
559
560     int             i_mb, i_b, i_mask;
561     macroblock_t *  p_mb;
562     f_addb_t        pf_addb;
563     elem_t *        p_data1;
564     elem_t *        p_data2;
565
566     *pi_mb_address += MacroblockAddressIncrement( p_vpar );
567
568     for( i_mb = i_mb_previous; i_mb < *pi_mb_address; i_mb++ )
569     {
570         /* Skipped macroblock (ISO/IEC 13818-2 7.6.6). */
571         static int          pi_dc_dct_reinit[4] = {128,256,512,1024};
572         static f_motion_t   pf_motion_skipped[4] = {NULL, vdec_MotionField,
573                                 vdec_MotionField, vdec_MotionFrame};
574
575         /* Reset DC predictors (7.2.1). */
576         p_vpar->slice.pi_dc_dct_pred[0] = p_vpar->slice.pi_dc_dct_pred[1]
577             = p_vpar->slice.pi_dc_dct_pred[2]
578             = pi_dc_dct_reinit[p_vpar->picture.i_intra_dc_precision];
579
580         if( p_vpar->picture.i_coding_type == P_CODING_TYPE )
581         {
582             /* Reset motion vector predictors (ISO/IEC 13818-2 7.6.3.4). */
583             bzero( p_vpar->slice.pppi_pmv, 8*sizeof(int) );
584         }
585
586         if( (p_mb = p_vpar->picture.pp_mb[i_mb_base + i_mb] =
587              vpar_NewMacroblock( &p_vpar->vfifo )) == NULL )
588         {
589             p_vpar->picture.b_error = 1;
590             intf_ErrMsg("vpar error: macroblock list is empty !");
591             return;
592         }
593
594         InitMacroblock( p_vpar, p_mb );
595        
596         /* No IDCT nor AddBlock. */
597         for( i_b = 0; i_b < 12; i_b++ )
598         {
599             p_mb->pf_idct[i_b] = vdec_DummyIDCT;
600             p_mb->pf_addb[i_b] = vdec_DummyBlock;
601         }
602
603         /* Motion type is picture structure. */
604         p_mb->pf_motion = pf_motion_skipped[p_vpar->picture.i_structure];
605
606         /* Set the field we use for motion compensation */
607         p_mb->ppi_field_select[0][0] = p_mb->ppi_field_select[0][1]
608                                      = ( p_vpar->picture.i_current_structure == BOTTOM_FIELD );
609         
610         /* Predict from field of same parity. */
611         /* ??? */
612     }
613
614     /* Get a macroblock structure. */
615     if( (p_mb = p_vpar->picture.pp_mb[i_mb_base + *pi_mb_address] =
616          vpar_NewMacroblock( &p_vpar->vfifo )) == NULL )
617     {
618         p_vpar->picture.b_error = 1;
619         intf_ErrMsg("vpar error: macroblock list is empty !");
620         return;
621     }
622
623     InitMacroblock( p_vpar, p_mb );
624
625     /* Parse off macroblock_modes structure. */
626     MacroblockModes( p_vpar, p_mb );
627
628     if( p_vpar->mb.i_mb_type & MB_QUANT )
629     {
630         LoadQuantizerScale( p_vpar );
631     }
632
633     if( p_vpar->mb.i_mb_type & MB_MOTION_FORWARD )
634     {
635         (*p_vpar->sequence.pf_decode_mv)( p_vpar, p_mb, 0 );
636     }
637
638     if( p_vpar->mb.i_mb_type & MB_MOTION_BACKWARD )
639     {
640         (*p_vpar->sequence.pf_decode_mv)( p_vpar, p_mb, 1 );
641     }
642
643     if( p_vpar->picture.b_concealment_mv && (p_vpar->mb.i_mb_type & MB_INTRA) )
644     {
645         DumpBits( &p_vpar->bit_stream, 1 );
646     }
647
648     if( p_vpar->mb.i_mb_type & MB_PATTERN )
649     {
650         (*p_vpar->sequence.pf_decode_pattern)( p_vpar );
651     }
652     else
653     {
654         int     pi_coded_block_pattern[2] = {0,
655                     (1 << 4+2*p_vpar->sequence.i_chroma_nb_blocks) - 1};
656         p_vpar->mb.i_coded_block_pattern = pi_coded_block_pattern
657                                     [p_vpar->mb.i_mb_type & MB_INTRA];
658     }
659
660     pf_addb = ppf_addb_intra[p_vpar->mb.i_mb_type & MB_INTRA];
661
662     /*
663      * Effectively decode blocks.
664      */
665
666     i_mask = 1 << (3 + 2*p_vpar->sequence.i_chroma_nb_blocks);
667
668     /* luminance */
669     p_data1 = (elem_t*) p_mb->p_picture->p_y;
670              // + p_mb->i_l_x + p_mb->i_l_y*(p_vpar->sequence.i_width);
671
672     for( i_b = 0 ; i_b < 4 ; i_b++, i_mask >>= 1 )
673     {
674         if( p_vpar->mb.i_coded_block_pattern & i_mask )
675         {
676             memset( p_mb->ppi_blocks[i_b], 0, 64*sizeof(elem_t) );
677             (*pppf_decode_block[p_vpar->sequence.b_mpeg2]
678                                [p_vpar->mb.i_mb_type & MB_INTRA])
679                 ( p_vpar, p_mb, i_b );
680         
681             /* decode_block has already set pf_idct and pi_sparse_pos. */
682             p_mb->pf_addb[i_b] = pf_addb;
683      
684             /* Calculate block coordinates. */
685             p_mb->p_data[i_b] = p_data1
686                                  + pi_y[p_vpar->mb.b_dct_type][i_b]
687                                    * p_vpar->sequence.i_chroma_width;
688         }
689         else
690         {
691             /* Block not coded, so no IDCT, nor AddBlock */
692             p_mb->pf_addb[i_b] = vdec_DummyBlock;
693             p_mb->pf_idct[i_b] = vdec_DummyIDCT;
694         }
695     }
696
697     /* chrominance U */
698     p_data1 = (elem_t*) p_mb->p_picture->p_u
699               + (p_mb->i_c_x >> pi_chroma_hor[p_vpar->sequence.i_chroma_format])
700               + (p_mb->i_c_y >> pi_chroma_ver[p_vpar->sequence.i_chroma_format])
701                 * (p_vpar->sequence.i_chroma_width);
702     p_data2 = (elem_t*) p_mb->p_picture->p_v
703               + (p_mb->i_c_x >> pi_chroma_hor[p_vpar->sequence.i_chroma_format])
704               + (p_mb->i_c_y >> pi_chroma_ver[p_vpar->sequence.i_chroma_format])
705                 * (p_vpar->sequence.i_chroma_width);
706
707     for( i_b = 4; i_b < 4 + 2*p_vpar->sequence.i_chroma_nb_blocks;
708          i_b++, i_mask >>= 1 )
709     {
710         elem_t *    pp_data[2] = {p_data1, p_data2};
711
712         if( p_vpar->mb.i_coded_block_pattern & i_mask )
713         {
714             memset( p_mb->ppi_blocks[i_b], 0, 64*sizeof(elem_t) );
715             (*pppf_decode_block[p_vpar->sequence.b_mpeg2]
716                                [p_vpar->mb.i_mb_type & MB_INTRA])
717                 ( p_vpar, p_mb, i_b );
718
719             /* decode_block has already set pf_idct and pi_sparse_pos. */
720             p_mb->pf_addb[i_b] = pf_addb;
721
722             /* Calculate block coordinates. */
723             p_mb->p_data[i_b] = pp_data[i_b & 1]
724                                  + pi_y[p_vpar->mb.b_dct_type][i_b]
725                                    * p_vpar->sequence.i_chroma_width;
726         }
727         else
728         {
729             /* Block not coded, so no IDCT, nor AddBlock */
730             p_mb->pf_addb[i_b] = vdec_DummyBlock;
731             p_mb->pf_idct[i_b] = vdec_DummyIDCT;
732         }
733     }
734
735     if( !( p_vpar->mb.i_mb_type & MB_INTRA ) )
736     {
737         static int          pi_dc_dct_reinit[4] = {128,256,512,1024};
738
739         /* Reset DC predictors (7.2.1). */
740         p_vpar->slice.pi_dc_dct_pred[0] = p_vpar->slice.pi_dc_dct_pred[1]
741             = p_vpar->slice.pi_dc_dct_pred[2]
742             = pi_dc_dct_reinit[p_vpar->picture.i_intra_dc_precision];
743     }
744     else if( !p_vpar->picture.b_concealment_mv )
745     {
746         /* Reset MV predictors. */
747         bzero( p_vpar->slice.pppi_pmv, 8*sizeof(int) );
748     }
749 }
750
751 /*****************************************************************************
752  * InitMacroblock : Initialize macroblock values
753  *****************************************************************************/
754 static __inline__ void InitMacroblock( vpar_thread_t * p_vpar,
755                                        macroblock_t * p_mb )
756 {
757     static f_chroma_motion_t pf_chroma_motion[4] =
758     { NULL, vdec_Motion420, vdec_Motion422, vdec_Motion444 };
759     
760     p_mb->p_picture = p_vpar->picture.p_picture;
761     p_mb->i_structure = p_vpar->picture.i_structure;
762     p_mb->i_l_x = p_vpar->mb.i_l_x;
763     p_mb->i_l_y = p_vpar->mb.i_l_y;
764     p_mb->i_c_x = p_vpar->mb.i_c_x;
765     p_mb->i_c_y = p_vpar->mb.i_c_y;
766     p_mb->i_chroma_nb_blocks = p_vpar->sequence.i_chroma_nb_blocks;
767     p_mb->pf_chroma_motion = pf_chroma_motion[p_vpar->sequence.i_chroma_format];
768
769     p_mb->p_forward = p_vpar->sequence.p_forward;
770     p_mb->p_backward = p_vpar->sequence.p_backward;
771     
772     p_mb->i_addb_l_stride = p_mb->i_l_stride = p_vpar->picture.i_l_stride;
773     p_mb->i_addb_c_stride = p_mb->i_c_stride = p_vpar->picture.i_c_stride;
774
775     /* Update macroblock real position. */
776     p_vpar->mb.i_l_x += 16;
777     p_vpar->mb.i_l_y += (p_vpar->mb.i_l_x / p_vpar->sequence.i_width)
778                         * (2 - p_vpar->picture.b_frame_structure) * 16;
779     p_vpar->mb.i_l_x %= p_vpar->sequence.i_width;
780
781     p_vpar->mb.i_c_x += p_vpar->sequence.i_chroma_mb_width;
782     p_vpar->mb.i_c_y += (p_vpar->mb.i_c_x / p_vpar->sequence.i_chroma_width)
783                         * (2 - p_vpar->picture.b_frame_structure)
784                         * p_vpar->sequence.i_chroma_mb_height;
785     p_vpar->mb.i_c_x %= p_vpar->sequence.i_chroma_width;
786 }
787
788 /*****************************************************************************
789  * MacroblockAddressIncrement : Get the macroblock_address_increment field
790  *****************************************************************************/
791 static __inline__ int MacroblockAddressIncrement( vpar_thread_t * p_vpar )
792 {
793     /* Index in the lookup table mb_addr_inc */
794     int    i_index = ShowBits( &p_vpar->bit_stream, 11 );
795     p_vpar->mb.i_addr_inc = 0;
796     /* Test the presence of the escape character */
797     while( i_index == 8 )
798     {
799         DumpBits( &p_vpar->bit_stream, 11 );
800         p_vpar->mb.i_addr_inc += 33;
801         i_index = ShowBits( &p_vpar->bit_stream, 11 );
802     }
803     /* Affect the value from the lookup table */
804     p_vpar->mb.i_addr_inc += p_vpar->pl_mb_addr_inc[i_index].i_value;
805     /* Dump the good number of bits */
806     DumpBits( &p_vpar->bit_stream, p_vpar->pl_mb_addr_inc[i_index].i_length );
807 }
808
809 /*****************************************************************************
810  * MacroblockModes : Get the macroblock_modes structure
811  *****************************************************************************/
812 static __inline__ void MacroblockModes( vpar_thread_t * p_vpar,
813                                         macroblock_t * p_mb )
814 {
815     static f_motion_t   pf_motion[2][4] =
816         { {NULL, vdec_FieldRecon, vdec_16x8Recon, vdec_DMVRecon},
817           {NULL, vdec_FieldRecon, vdec_FrameRecon, vdec_DMVRecon} };
818     static int          ppi_mv_count[2][4] = { {0, 1, 2, 1}, {0, 2, 1, 1} };
819     static int          ppi_mv_format[2][4] = { {0, 1, 1, 1}, {0, 1, 2, 1} };
820
821     /* Get macroblock_type. */
822     p_vpar->mb.i_mb_type = (p_vpar->picture.pf_macroblock_type)( p_vpar );
823
824     /* SCALABILITY : warning, we don't know if spatial_temporal_weight_code
825      * has to be dropped, take care if you use scalable streams. */
826     /* DumpBits( &p_vpar->bit_stream, 2 ); */
827     
828     if( (p_vpar->picture.i_coding_type == P_CODING_TYPE) &&
829         !(p_vpar->mb.i_mb_type & (MB_MOTION_FORWARD|MB_INTRA)) )
830     {
831         /* Special No-MC macroblock in P pictures (7.6.3.5). */
832         memset( p_vpar->slice.pppi_pmv, 0, 8*sizeof(int) );
833         memset( p_mb->pppi_motion_vectors, 0, 8*sizeof(int) );
834         
835         p_vpar->mb.i_motion_type = MOTION_FRAME;
836         p_mb->ppi_field_select[0][0] = ( p_vpar->picture.i_current_structure == BOTTOM_FIELD );
837     }
838     else if( !(p_vpar->mb.i_mb_type & (MB_MOTION_FORWARD | MB_MOTION_BACKWARD))
839              || p_vpar->picture.b_frame_pred_frame_dct )
840     {
841         /* If mb_type has neither MOTION_FORWARD nor MOTION_BACKWARD, this
842          * is useless, but also harmless. */
843         p_vpar->mb.i_motion_type = MOTION_FRAME;
844     }
845     else
846     {
847         p_vpar->mb.i_motion_type = GetBits( &p_vpar->bit_stream, 2 );
848     }
849
850      if( p_vpar->mb.i_mb_type & MB_INTRA )
851     {
852         /* For the intra macroblocks, we use an empty motion
853          * compensation function */
854         p_mb->pf_motion = vdec_DummyRecon;
855     }
856     else
857     {
858         p_mb->pf_motion = pf_motion[p_vpar->picture.b_frame_structure]
859                                    [p_vpar->mb.i_motion_type];
860     }
861     p_vpar->mb.i_mv_count = ppi_mv_count[p_vpar->picture.b_frame_structure]
862                                         [p_vpar->mb.i_motion_type];
863     p_vpar->mb.i_mv_format = ppi_mv_format[p_vpar->picture.b_frame_structure]
864                                           [p_vpar->mb.i_motion_type];
865
866     p_vpar->mb.b_dct_type = 0;
867     if( (p_vpar->picture.i_structure == FRAME_STRUCTURE) &&
868         (!p_vpar->picture.b_frame_pred_frame_dct) &&
869         (p_vpar->mb.i_mb_type & (MB_PATTERN|MB_INTRA)) )
870     {
871         if( p_vpar->mb.b_dct_type = GetBits( &p_vpar->bit_stream, 1 ) )
872         {
873             /* The DCT is coded on fields. Jump one line between each
874              * sample. */
875             p_mb->i_addb_l_stride <<= 1;
876             p_mb->i_addb_l_stride += 8;
877             /* With CHROMA_420, the DCT is necessarily frame-coded. */
878             if( p_vpar->sequence.i_chroma_format != CHROMA_420 )
879             {
880                 p_mb->i_addb_c_stride <<= 1;
881                 p_mb->i_addb_c_stride += 8;
882             }
883         }
884     }
885
886     p_vpar->mb.b_dmv = p_vpar->mb.i_motion_type == MOTION_DMV;
887 }
888
889 /*****************************************************************************
890  * vpar_IMBType : macroblock_type in I pictures
891  *****************************************************************************/
892 int vpar_IMBType( vpar_thread_t * p_vpar )
893 {
894     /* Take two bits for testing */
895     int                 i_type = ShowBits( &p_vpar->bit_stream, 2 );
896     /* Lookup table for macroblock_type */
897     static lookup_t     pl_mb_Itype[4] = { {MB_ERROR, 0},
898                                            {MB_QUANT|MB_INTRA, 2},
899                                            {MB_INTRA, 1},
900                                            {MB_INTRA, 2} };
901     /* Dump the good number of bits */
902     DumpBits( &p_vpar->bit_stream, pl_mb_Itype[i_type].i_length );
903     return pl_mb_Itype[i_type].i_value;
904 }
905
906 /*****************************************************************************
907  * vpar_PMBType : macroblock_type in P pictures
908  *****************************************************************************/
909 int vpar_PMBType( vpar_thread_t * p_vpar )
910 {
911     /* Testing on 6 bits */
912     int                i_type = ShowBits( &p_vpar->bit_stream, 6 );
913     /* Dump the good number of bits */
914     DumpBits( &p_vpar->bit_stream, p_vpar->ppl_mb_type[0][i_type].i_length );
915     /* return the value from the lookup table for P type */
916     return p_vpar->ppl_mb_type[0][i_type].i_value;
917 }
918
919 /*****************************************************************************
920  * vpar_BMBType : macroblock_type in B pictures
921  *****************************************************************************/
922 int vpar_BMBType( vpar_thread_t * p_vpar )
923 {
924      /* Testing on 6 bits */
925     int                i_type = ShowBits( &p_vpar->bit_stream, 6 );
926     /* Dump the good number of bits */
927     DumpBits( &p_vpar->bit_stream, p_vpar->ppl_mb_type[1][i_type].i_length );
928     /* return the value from the lookup table for B type */
929     return p_vpar->ppl_mb_type[1][i_type].i_value;
930 }
931
932 /*****************************************************************************
933  * vpar_DMBType : macroblock_type in D pictures
934  *****************************************************************************/
935 int vpar_DMBType( vpar_thread_t * p_vpar )
936 {
937     /* Taking 1 bit */
938     int               i_type = GetBits( &p_vpar->bit_stream, 1 );
939     /* Lookup table */
940     static int        pi_mb_Dtype[2] = { MB_ERROR, 1 };
941     return pi_mb_Dtype[i_type];
942 }
943
944 /*****************************************************************************
945  * vpar_CodedPattern420 : coded_block_pattern with 420 chroma
946  *****************************************************************************/
947 int vpar_CodedPattern420( vpar_thread_t * p_vpar )
948 {
949     /* Take the max 9 bits length vlc code for testing */
950     int      i_vlc = ShowBits( &p_vpar->bit_stream, 9 );
951     /* Trash the good number of bits read in  the lookup table */
952     DumpBits( &p_vpar->bit_stream, p_vpar->pl_coded_pattern[i_vlc].i_length );
953     /* return the value from the vlc table */
954     return p_vpar->pl_coded_pattern[i_vlc].i_value;
955 }
956
957 /*****************************************************************************
958  * vpar_CodedPattern422 : coded_block_pattern with 422 chroma
959  *****************************************************************************/
960 int vpar_CodedPattern422( vpar_thread_t * p_vpar )
961 {
962     int      i_vlc = ShowBits( &p_vpar->bit_stream, 9 );
963     /* Supplementary 2 bits long code for 422 format */
964     int      i_coded_block_pattern_1;
965     DumpBits( &p_vpar->bit_stream, p_vpar->pl_coded_pattern[i_vlc].i_length );
966     i_coded_block_pattern_1 = GetBits( &p_vpar->bit_stream, 2 );
967     /* the code is just to be added to the value found in the table */
968     return p_vpar->pl_coded_pattern[i_vlc].i_value |
969            (i_coded_block_pattern_1 << 6);
970 }
971
972 /*****************************************************************************
973  * vpar_CodedPattern444 : coded_block_pattern with 444 chroma
974  *****************************************************************************/
975 int vpar_CodedPattern444( vpar_thread_t * p_vpar )
976 {
977     int      i_vlc = ShowBits( &p_vpar->bit_stream, 9 );
978     int      i_coded_block_pattern_2;
979     DumpBits( &p_vpar->bit_stream, p_vpar->pl_coded_pattern[i_vlc].i_length );
980     i_coded_block_pattern_2 = GetBits( &p_vpar->bit_stream, 6 );
981     return p_vpar->pl_coded_pattern[i_vlc].i_value |
982            ( i_coded_block_pattern_2 << 6 );
983 }
984
985 /*****************************************************************************
986  * vpar_DecodeMPEG1Non : decode MPEG-1 non-intra blocks
987  *****************************************************************************/
988 static void vpar_DecodeMPEG1Non( vpar_thread_t * p_vpar, macroblock_t * p_mb, int i_b )
989 {
990     /* À pomper dans Berkeley. Pour toutes ces fonctions, il faut mettre
991        p_mb->pf_idct[i_b] à :
992         - vdec_IDCT ou
993         - vdec_SparseIDCT si la matrice n'a qu'un coefficient non nul.
994        Dans le deuxième cas, p_mb->pi_sparse_pos[i_b] contient le numéro
995        de ce coefficient. */
996
997     if( p_vpar->picture.i_coding_type == D_CODING_TYPE )
998     {
999         /* Remove end_of_macroblock (always 1, prevents startcode emulation)
1000          * ISO/IEC 11172-2 section 2.4.2.7 and 2.4.3.6 */
1001         DumpBits( &p_vpar->bit_stream, 1 );
1002     }
1003 }
1004
1005 /*****************************************************************************
1006  * vpar_DecodeMPEG1Intra : decode MPEG-1 intra blocks
1007  *****************************************************************************/
1008 static void vpar_DecodeMPEG1Intra( vpar_thread_t * p_vpar, macroblock_t * p_mb, int i_b )
1009 {
1010     /* À pomper dans Berkeley. */
1011
1012     if( p_vpar->picture.i_coding_type == D_CODING_TYPE )
1013     {
1014         /* Remove end_of_macroblock (always 1, prevents startcode emulation)
1015          * ISO/IEC 11172-2 section 2.4.2.7 and 2.4.3.6 */
1016         DumpBits( &p_vpar->bit_stream, 1 );
1017     }
1018 }
1019
1020 /*****************************************************************************
1021  * vpar_DecodeMPEG2Non : decode MPEG-2 non-intra blocks
1022  *****************************************************************************/
1023 static void vpar_DecodeMPEG2Non( vpar_thread_t * p_vpar, macroblock_t * p_mb, int i_b )
1024 {
1025     /* À pomper dans Berkeley. Bien sûr les matrices seront différentes... */
1026 }
1027
1028 /*****************************************************************************
1029  * vpar_DecodeMPEG2Intra : decode MPEG-2 intra blocks
1030  *****************************************************************************/
1031 static void vpar_DecodeMPEG2Intra( vpar_thread_t * p_vpar, macroblock_t * p_mb, int i_b )
1032 {
1033     int         i_dummy;
1034     int         i_nc;
1035     int         i_cc;
1036     int         i_coef;
1037     int         i_type;
1038     int         i_code;
1039     int         i_code5;
1040     int         i_select;
1041     int         i_offset;
1042     int         i_pos;
1043     int         i_dct_dc_size;
1044     int         i_dct_dc_diff;
1045     int         i_run;
1046     int         i_level;
1047     int         i_2b;
1048     int         i_4b;
1049     boolean_t   b_sign;
1050     /* Lookup Table for the chromatic component */
1051     static int pi_cc_index[12] = { 0, 0, 0, 0, 1, 2, 1, 2, 1, 2 };
1052     /* Lookup table for the offset in the tables for ac coef decoding */
1053     static lookup_t pl_offset_table[8] = { { 12, 4 }, { 8, 4 }, { 6, 8 },
1054                                            { 4, 16 }, { 3, 16 }, { 2, 16 },
1055                                            { 1, 16 }, { 0, 16 } };
1056     
1057     i_cc = pi_cc_index[i_b];
1058     /* Determine whether it is luminance or not (chrominance) */
1059     i_type = ( i_cc + 1 ) >> 1;
1060
1061     /* Decoding of the DC intra coefficient */
1062     /* The nb of bits to parse depends on i_type */
1063     i_code = ShowBits( &p_vpar->bit_stream, 9 + i_type );
1064     /* To reduce memory occupation, there are two lookup tables
1065      * See InitDCT above */
1066     i_code5 = i_code >> 4;
1067     /* Shall we lookup in the first or in the second table ? */
1068     i_select = ( i_code5 - 1 ) / 31;
1069     /* Offset value for looking in the second table */
1070     i_offset = 0x1f0 + ( i_type * 0x1f0 );
1071     i_pos = ( i_code5 * ( ! i_select ) ) + ( ( i_code - i_offset ) * i_select );
1072     i_dct_dc_size = p_vpar->pppl_dct_dc_size[i_type][i_select][i_pos].i_value;
1073     /* Dump the variable length code */
1074     DumpBits( &p_vpar->bit_stream, 
1075               p_vpar->pppl_dct_dc_size[i_type][i_select][i_pos].i_length );
1076     /* Read the actual code with the good length */
1077     i_dct_dc_diff = GetBits( &p_vpar->bit_stream, i_dct_dc_size );
1078     p_vpar->slice.pi_dc_dct_pred[i_cc] += i_dct_dc_diff;
1079 #ifndef VDEC_DFT
1080     p_mb->ppi_blocks[i_b][0] = ( p_vpar->slice.pi_dc_dct_pred[i_cc] <<
1081                                ( 3 - p_vpar->picture.i_intra_dc_precision ) );
1082     i_nc = 1;
1083 #else
1084     p_mb->ppi_blocks[i_b][0] = ( p_vpar->slice.pi_dc_dct_pred[i_cc] <<
1085                                ( 11 - p_vpar->picture.i_intra_dc_precision ) );
1086     i_nc = ( p_vpar->slice.pi_dc_dct_pred[i_cc] != 0 );
1087 #endif
1088
1089     /* Decoding of the AC coefficients */
1090     i_coef = 0;
1091     for( i_dummy = 1; i_dummy < 64; i_dummy++ )
1092     {
1093         i_code = ShowBits( &p_vpar->bit_stream, 16 );
1094         i_2b = ! ( i_code >> 14 );
1095         i_4b = ! ( ( i_code & 0x3C00 ) >> 10 );
1096         /* Will contain the number of table to use in ppl_dct_coef */
1097         i_select = i_2b + i_4b + ( i_code & 0x0200 ) + ( i_code & 0x0100 ) + 
1098                                  ( i_code & 0x0080 ) + ( i_code & 0x0040 ) +
1099                                  ( i_code & 0x0020 ) + ( i_code & 0x0010 );
1100         /* Shall we use B15 table or not */
1101         i_type = p_vpar->picture.b_intra_vlc_format & ( ( i_select & 3 ) >> 1 );
1102         i_select += (i_type * 7 );
1103         i_offset = ( i_code >> pl_offset_table[i_select].i_value ) - 
1104                                pl_offset_table[i_select].i_length;
1105         i_run = ( * ( p_vpar->ppl_dct_coef[i_select]
1106                     + ( i_offset * sizeof(dct_lookup_t) ) ) ).i_run;
1107         DumpBits( &p_vpar->bit_stream, 
1108                   ( * ( p_vpar->ppl_dct_coef[i_select]
1109                     + ( i_offset * sizeof(dct_lookup_t) ) ) ).i_length );
1110         switch( i_run )
1111         {
1112             case DCT_ESCAPE:
1113                 i_run = GetBits( &p_vpar->bit_stream, 6 );
1114                 i_level = GetBits( &p_vpar->bit_stream, 12 );
1115                 p_mb->ppi_blocks[i_b][i_dummy] = ( b_sign = ( i_level > 2047 ) ) 
1116                                                           ? ( -4096 + i_level )
1117                                                           : i_level;
1118                 i_coef = i_dummy;
1119                 i_dummy += i_run;
1120                 i_nc ++;
1121                 break;
1122             case DCT_EOB:
1123                 i_dummy = 64;
1124                 break;
1125             default:
1126                 i_level = ( * ( p_vpar->ppl_dct_coef[i_select]
1127                               + ( i_offset * sizeof(dct_lookup_t) ) ) ).i_level;
1128                 b_sign = Getbits( &p_vpar->bit_stream, 1 );
1129                 p_mb->ppi_blocks[i_b][i_dummy] = b_sign ? -i_level : i_level;
1130                 i_coef = i_dummy;
1131                 i_dummy += i_run;
1132                 i_nc ++;
1133         }
1134     }
1135     if( i_nc == 1 )
1136     {
1137         p_mb->pf_idct[i_b] = vdec_SparseIDCT;
1138         p_mb->pi_sparse_pos[i_b] = i_coef;
1139     }
1140     else
1141     {
1142         p_mb->pf_idct[i_b] = vdec_IDCT;
1143     }
1144 }