[cvs] / xvidcore / src / decoder.c Repository:
ViewVC logotype

Annotation of /xvidcore/src/decoder.c

Parent Directory Parent Directory | Revision Log Revision Log

Revision 1.37.2.30 - (view) (download)

1 : Isibaar 1.1 /**************************************************************************
2 :     *
3 : edgomez 1.12 * XVID MPEG-4 VIDEO CODEC
4 :     * - Decoder main module -
5 : Isibaar 1.1 *
6 : edgomez 1.12 * This program is an implementation of a part of one or more MPEG-4
7 :     * Video tools as specified in ISO/IEC 14496-2 standard. Those intending
8 :     * to use this software module in hardware or software products are
9 :     * advised that its use may infringe existing patents or copyrights, and
10 :     * any such use would be at such party's own risk. The original
11 :     * developer of this software module and his/her company, and subsequent
12 :     * editors and their companies, will have no liability for use of this
13 :     * software or modifications or derivatives thereof.
14 :     *
15 :     * This program is free software; you can redistribute it and/or modify
16 :     * it under the terms of the GNU General Public License as published by
17 :     * the Free Software Foundation; either version 2 of the License, or
18 :     * (at your option) any later version.
19 :     *
20 :     * This program is distributed in the hope that it will be useful,
21 :     * but WITHOUT ANY WARRANTY; without even the implied warranty of
22 :     * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
23 :     * GNU General Public License for more details.
24 :     *
25 :     * You should have received a copy of the GNU General Public License
26 :     * along with this program; if not, write to the Free Software
27 :     * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
28 : Isibaar 1.1 *
29 :     *************************************************************************/
30 :    
31 :     /**************************************************************************
32 :     *
33 : edgomez 1.12 * History:
34 : Isibaar 1.1 *
35 : chenm001 1.30 * 15.07.2002 fix a bug in B-frame decode at DIRECT mode
36 :     * MinChen <chenm001@163.com>
37 : chenm001 1.27 * 10.07.2002 added BFRAMES_DEC_DEBUG support
38 :     * Fix a little bug for low_delay flage
39 :     * MinChen <chenm001@163.com>
40 : chenm001 1.26 * 28.06.2002 added basic resync support to iframe/pframe_decode()
41 : chenm001 1.27 * 22.06.2002 added primative N_VOP support
42 : chenm001 1.30 * #define BFRAMES_DEC now enables Minchen's bframe decoder
43 : chenm001 1.17 * 08.05.2002 add low_delay support for B_VOP decode
44 :     * MinChen <chenm001@163.com>
45 :     * 05.05.2002 fix some B-frame decode problem
46 : chenm001 1.14 * 02.05.2002 add B-frame decode support(have some problem);
47 :     * MinChen <chenm001@163.com>
48 : chenm001 1.11 * 22.04.2002 add some B-frame decode support; chenm001 <chenm001@163.com>
49 : h 1.8 * 29.03.2002 interlacing fix - compensated block wasn't being used when
50 :     * reconstructing blocks, thus artifacts
51 :     * interlacing speedup - used transfers to re-interlace
52 :     * interlaced decoding should be as fast as progressive now
53 : h 1.5 * 26.03.2002 interlacing support - moved transfers outside decode loop
54 : edgomez 1.12 * 26.12.2001 decoder_mbinter: dequant/idct moved within if(coded) block
55 :     * 22.12.2001 lock based interpolation
56 :     * 01.12.2001 inital version; (c)2001 peter ross <pross@cs.rmit.edu.au>
57 :     *
58 : chl 1.37.2.30 * $Id: decoder.c,v 1.37.2.29 2003/01/11 21:22:24 chl Exp $
59 : Isibaar 1.1 *
60 :     *************************************************************************/
61 :    
62 : chl 1.37.2.27 #include <stdio.h>
63 : Isibaar 1.1 #include <stdlib.h>
64 : edgomez 1.12 #include <string.h>
65 : Isibaar 1.1
66 : chenm001 1.27 #ifdef BFRAMES_DEC_DEBUG
67 :     #define BFRAMES_DEC
68 :     #endif
69 :    
70 : Isibaar 1.1 #include "xvid.h"
71 :     #include "portab.h"
72 : suxen_drol 1.37.2.26 #include "global.h"
73 : Isibaar 1.1
74 :     #include "decoder.h"
75 :     #include "bitstream/bitstream.h"
76 :     #include "bitstream/mbcoding.h"
77 :    
78 :     #include "quant/quant_h263.h"
79 :     #include "quant/quant_mpeg4.h"
80 :     #include "dct/idct.h"
81 :     #include "dct/fdct.h"
82 :     #include "utils/mem_transfer.h"
83 :     #include "image/interpolate8x8.h"
84 : suxen_drol 1.37.2.15 #include "image/reduced.h"
85 : suxen_drol 1.37.2.19 #include "image/font.h"
86 : Isibaar 1.1
87 :     #include "bitstream/mbcoding.h"
88 :     #include "prediction/mbprediction.h"
89 :     #include "utils/timer.h"
90 :     #include "utils/emms.h"
91 : Isibaar 1.37.2.6 #include "motion/motion.h"
92 : Isibaar 1.1
93 :     #include "image/image.h"
94 :     #include "image/colorspace.h"
95 : Isibaar 1.3 #include "utils/mem_align.h"
96 : Isibaar 1.1
97 : edgomez 1.19 int
98 : suxen_drol 1.37.2.8 decoder_resize(DECODER * dec)
99 : Isibaar 1.1 {
100 : suxen_drol 1.37.2.8 /* free existing */
101 : Isibaar 1.1
102 : suxen_drol 1.37.2.21 image_destroy(&dec->cur, dec->edged_width, dec->edged_height);
103 : suxen_drol 1.37.2.8 image_destroy(&dec->refn[0], dec->edged_width, dec->edged_height);
104 :     image_destroy(&dec->refn[1], dec->edged_width, dec->edged_height);
105 : suxen_drol 1.37.2.21 image_destroy(&dec->tmp, dec->edged_width, dec->edged_height);
106 :     image_destroy(&dec->qtmp, dec->edged_width, dec->edged_height);
107 : Isibaar 1.1
108 : chl 1.37.2.27 image_destroy(&dec->gmc, dec->edged_width, dec->edged_height);
109 :    
110 : suxen_drol 1.37.2.8 if (dec->last_mbs)
111 :     xvid_free(dec->last_mbs);
112 :     if (dec->mbs)
113 :     xvid_free(dec->mbs);
114 :    
115 :     /* realloc */
116 : Isibaar 1.1
117 :     dec->mb_width = (dec->width + 15) / 16;
118 :     dec->mb_height = (dec->height + 15) / 16;
119 :    
120 :     dec->edged_width = 16 * dec->mb_width + 2 * EDGE_SIZE;
121 :     dec->edged_height = 16 * dec->mb_height + 2 * EDGE_SIZE;
122 : edgomez 1.19
123 :     if (image_create(&dec->cur, dec->edged_width, dec->edged_height)) {
124 : Isibaar 1.3 xvid_free(dec);
125 : Isibaar 1.1 return XVID_ERR_MEMORY;
126 :     }
127 :    
128 : edgomez 1.19 if (image_create(&dec->refn[0], dec->edged_width, dec->edged_height)) {
129 : Isibaar 1.1 image_destroy(&dec->cur, dec->edged_width, dec->edged_height);
130 : Isibaar 1.3 xvid_free(dec);
131 : Isibaar 1.1 return XVID_ERR_MEMORY;
132 :     }
133 : suxen_drol 1.37.2.8
134 : chenm001 1.11 // add by chenm001 <chenm001@163.com>
135 :     // for support B-frame to reference last 2 frame
136 : edgomez 1.19 if (image_create(&dec->refn[1], dec->edged_width, dec->edged_height)) {
137 : chenm001 1.11 image_destroy(&dec->cur, dec->edged_width, dec->edged_height);
138 :     image_destroy(&dec->refn[0], dec->edged_width, dec->edged_height);
139 :     xvid_free(dec);
140 :     return XVID_ERR_MEMORY;
141 :     }
142 : suxen_drol 1.37.2.21 if (image_create(&dec->tmp, dec->edged_width, dec->edged_height)) {
143 : chenm001 1.14 image_destroy(&dec->cur, dec->edged_width, dec->edged_height);
144 :     image_destroy(&dec->refn[0], dec->edged_width, dec->edged_height);
145 :     image_destroy(&dec->refn[1], dec->edged_width, dec->edged_height);
146 :     xvid_free(dec);
147 :     return XVID_ERR_MEMORY;
148 :     }
149 : Isibaar 1.1
150 : suxen_drol 1.37.2.21 if (image_create(&dec->qtmp, dec->edged_width, dec->edged_height)) {
151 : Isibaar 1.37.2.2 image_destroy(&dec->cur, dec->edged_width, dec->edged_height);
152 :     image_destroy(&dec->refn[0], dec->edged_width, dec->edged_height);
153 :     image_destroy(&dec->refn[1], dec->edged_width, dec->edged_height);
154 : suxen_drol 1.37.2.21 image_destroy(&dec->tmp, dec->edged_width, dec->edged_height);
155 : Isibaar 1.37.2.2 xvid_free(dec);
156 :     return XVID_ERR_MEMORY;
157 :     }
158 :    
159 : chl 1.37.2.27 if (image_create(&dec->gmc, dec->edged_width, dec->edged_height)) {
160 :     image_destroy(&dec->qtmp, dec->edged_width, dec->edged_height);
161 :     image_destroy(&dec->cur, dec->edged_width, dec->edged_height);
162 :     image_destroy(&dec->refn[0], dec->edged_width, dec->edged_height);
163 :     image_destroy(&dec->refn[1], dec->edged_width, dec->edged_height);
164 :     image_destroy(&dec->tmp, dec->edged_width, dec->edged_height);
165 :     xvid_free(dec);
166 :     return XVID_ERR_MEMORY;
167 :     }
168 :    
169 : edgomez 1.19 dec->mbs =
170 :     xvid_malloc(sizeof(MACROBLOCK) * dec->mb_width * dec->mb_height,
171 :     CACHE_LINE);
172 :     if (dec->mbs == NULL) {
173 : Isibaar 1.1 image_destroy(&dec->cur, dec->edged_width, dec->edged_height);
174 : chenm001 1.14 image_destroy(&dec->refn[0], dec->edged_width, dec->edged_height);
175 :     image_destroy(&dec->refn[1], dec->edged_width, dec->edged_height);
176 : suxen_drol 1.37.2.21 image_destroy(&dec->tmp, dec->edged_width, dec->edged_height);
177 :     image_destroy(&dec->qtmp, dec->edged_width, dec->edged_height);
178 : chenm001 1.14 xvid_free(dec);
179 :     return XVID_ERR_MEMORY;
180 :     }
181 : Isibaar 1.20 memset(dec->mbs, 0, sizeof(MACROBLOCK) * dec->mb_width * dec->mb_height);
182 :    
183 : chenm001 1.14 // add by chenm001 <chenm001@163.com>
184 :     // for skip MB flag
185 : edgomez 1.19 dec->last_mbs =
186 :     xvid_malloc(sizeof(MACROBLOCK) * dec->mb_width * dec->mb_height,
187 :     CACHE_LINE);
188 :     if (dec->last_mbs == NULL) {
189 : chenm001 1.14 xvid_free(dec->mbs);
190 :     image_destroy(&dec->cur, dec->edged_width, dec->edged_height);
191 :     image_destroy(&dec->refn[0], dec->edged_width, dec->edged_height);
192 :     image_destroy(&dec->refn[1], dec->edged_width, dec->edged_height);
193 : suxen_drol 1.37.2.21 image_destroy(&dec->tmp, dec->edged_width, dec->edged_height);
194 :     image_destroy(&dec->qtmp, dec->edged_width, dec->edged_height);
195 : Isibaar 1.3 xvid_free(dec);
196 : Isibaar 1.1 return XVID_ERR_MEMORY;
197 :     }
198 : Isibaar 1.20
199 :     memset(dec->last_mbs, 0, sizeof(MACROBLOCK) * dec->mb_width * dec->mb_height);
200 : Isibaar 1.1
201 : suxen_drol 1.37.2.8 return XVID_ERR_OK;
202 :     }
203 :    
204 :    
205 :     int
206 :     decoder_create(XVID_DEC_PARAM * param)
207 :     {
208 :     DECODER *dec;
209 :    
210 :     dec = xvid_malloc(sizeof(DECODER), CACHE_LINE);
211 :     if (dec == NULL) {
212 :     return XVID_ERR_MEMORY;
213 :     }
214 :     memset(dec, 0, sizeof(DECODER));
215 :    
216 :     param->handle = dec;
217 :    
218 :     dec->width = param->width;
219 :     dec->height = param->height;
220 :    
221 :     image_null(&dec->cur);
222 :     image_null(&dec->refn[0]);
223 :     image_null(&dec->refn[1]);
224 : suxen_drol 1.37.2.21 image_null(&dec->tmp);
225 :     image_null(&dec->qtmp);
226 : suxen_drol 1.37.2.8
227 : chl 1.37.2.27 /* image based GMC */
228 :     image_null(&dec->gmc);
229 :    
230 :    
231 : suxen_drol 1.37.2.8 dec->mbs = NULL;
232 :     dec->last_mbs = NULL;
233 :    
234 : Isibaar 1.1 init_timer();
235 :    
236 : chenm001 1.14 // add by chenm001 <chenm001@163.com>
237 :     // for support B-frame to save reference frame's time
238 : suxen_drol 1.37.2.19 dec->frames = 0;
239 : chenm001 1.14 dec->time = dec->time_base = dec->last_time_base = 0;
240 : suxen_drol 1.37.2.8 dec->low_delay = 0;
241 : suxen_drol 1.37.2.15 dec->packed_mode = 0;
242 : edgomez 1.19
243 : suxen_drol 1.37.2.8 dec->fixed_dimensions = (dec->width > 0 && dec->height > 0);
244 :    
245 :     if (dec->fixed_dimensions)
246 :     return decoder_resize(dec);
247 :     else
248 :     return XVID_ERR_OK;
249 : Isibaar 1.1 }
250 :    
251 :    
252 : edgomez 1.19 int
253 :     decoder_destroy(DECODER * dec)
254 : Isibaar 1.1 {
255 : chenm001 1.14 xvid_free(dec->last_mbs);
256 : Isibaar 1.3 xvid_free(dec->mbs);
257 : chl 1.37.2.27
258 :     image_destroy(&dec->gmc, dec->edged_width, dec->edged_height); /* image based GMC */
259 :    
260 : chenm001 1.11 image_destroy(&dec->refn[0], dec->edged_width, dec->edged_height);
261 : chenm001 1.14 image_destroy(&dec->refn[1], dec->edged_width, dec->edged_height);
262 : suxen_drol 1.37.2.21 image_destroy(&dec->tmp, dec->edged_width, dec->edged_height);
263 :     image_destroy(&dec->qtmp, dec->edged_width, dec->edged_height);
264 : Isibaar 1.1 image_destroy(&dec->cur, dec->edged_width, dec->edged_height);
265 : Isibaar 1.3 xvid_free(dec);
266 : Isibaar 1.1
267 :     write_timer();
268 :     return XVID_ERR_OK;
269 :     }
270 :    
271 :    
272 :    
273 : edgomez 1.19 static const int32_t dquant_table[4] = {
274 : Isibaar 1.1 -1, -2, 1, 2
275 :     };
276 :    
277 :    
278 : chenm001 1.14
279 :    
280 : Isibaar 1.1 // decode an intra macroblock
281 :    
282 : edgomez 1.19 void
283 :     decoder_mbintra(DECODER * dec,
284 :     MACROBLOCK * pMB,
285 :     const uint32_t x_pos,
286 :     const uint32_t y_pos,
287 :     const uint32_t acpred_flag,
288 :     const uint32_t cbp,
289 :     Bitstream * bs,
290 :     const uint32_t quant,
291 : chenm001 1.26 const uint32_t intra_dc_threshold,
292 : suxen_drol 1.37.2.15 const unsigned int bound,
293 :     const int reduced_resolution)
294 : Isibaar 1.1 {
295 : edgomez 1.7
296 :     DECLARE_ALIGNED_MATRIX(block, 6, 64, int16_t, CACHE_LINE);
297 : edgomez 1.19 DECLARE_ALIGNED_MATRIX(data, 6, 64, int16_t, CACHE_LINE);
298 : edgomez 1.7
299 : h 1.8 uint32_t stride = dec->edged_width;
300 :     uint32_t stride2 = stride / 2;
301 :     uint32_t next_block = stride * 8;
302 : h 1.5 uint32_t i;
303 :     uint32_t iQuant = pMB->quant;
304 :     uint8_t *pY_Cur, *pU_Cur, *pV_Cur;
305 :    
306 : suxen_drol 1.37.2.15 if (reduced_resolution) {
307 :     pY_Cur = dec->cur.y + (y_pos << 5) * stride + (x_pos << 5);
308 :     pU_Cur = dec->cur.u + (y_pos << 4) * stride2 + (x_pos << 4);
309 :     pV_Cur = dec->cur.v + (y_pos << 4) * stride2 + (x_pos << 4);
310 :     }else{
311 :     pY_Cur = dec->cur.y + (y_pos << 4) * stride + (x_pos << 4);
312 :     pU_Cur = dec->cur.u + (y_pos << 3) * stride2 + (x_pos << 3);
313 :     pV_Cur = dec->cur.v + (y_pos << 3) * stride2 + (x_pos << 3);
314 :     }
315 : edgomez 1.7
316 : edgomez 1.19 memset(block, 0, 6 * 64 * sizeof(int16_t)); // clear
317 : h 1.5
318 : edgomez 1.19 for (i = 0; i < 6; i++) {
319 : h 1.5 uint32_t iDcScaler = get_dc_scaler(iQuant, i < 4);
320 : Isibaar 1.1 int16_t predictors[8];
321 :     int start_coeff;
322 :    
323 :     start_timer();
324 : edgomez 1.19 predict_acdc(dec->mbs, x_pos, y_pos, dec->mb_width, i, &block[i * 64],
325 : chenm001 1.26 iQuant, iDcScaler, predictors, bound);
326 : edgomez 1.19 if (!acpred_flag) {
327 : h 1.5 pMB->acpred_directions[i] = 0;
328 : Isibaar 1.1 }
329 :     stop_prediction_timer();
330 :    
331 : edgomez 1.19 if (quant < intra_dc_threshold) {
332 : Isibaar 1.1 int dc_size;
333 :     int dc_dif;
334 :    
335 : edgomez 1.19 dc_size = i < 4 ? get_dc_size_lum(bs) : get_dc_size_chrom(bs);
336 :     dc_dif = dc_size ? get_dc_dif(bs, dc_size) : 0;
337 : Isibaar 1.1
338 : edgomez 1.19 if (dc_size > 8) {
339 :     BitstreamSkip(bs, 1); // marker
340 : Isibaar 1.1 }
341 : edgomez 1.19
342 :     block[i * 64 + 0] = dc_dif;
343 : Isibaar 1.1 start_coeff = 1;
344 : chenm001 1.26
345 :     DPRINTF(DPRINTF_COEFF,"block[0] %i", dc_dif);
346 : edgomez 1.19 } else {
347 : Isibaar 1.1 start_coeff = 0;
348 :     }
349 :    
350 :     start_timer();
351 : edgomez 1.19 if (cbp & (1 << (5 - i))) // coded
352 : Isibaar 1.1 {
353 : h 1.37.2.1 int direction = dec->alternate_vertical_scan ?
354 :     2 : pMB->acpred_directions[i];
355 :    
356 :     get_intra_block(bs, &block[i * 64], direction, start_coeff);
357 : Isibaar 1.1 }
358 :     stop_coding_timer();
359 :    
360 :     start_timer();
361 : edgomez 1.19 add_acdc(pMB, i, &block[i * 64], iDcScaler, predictors);
362 : Isibaar 1.1 stop_prediction_timer();
363 :    
364 :     start_timer();
365 : edgomez 1.19 if (dec->quant_type == 0) {
366 :     dequant_intra(&data[i * 64], &block[i * 64], iQuant, iDcScaler);
367 :     } else {
368 :     dequant4_intra(&data[i * 64], &block[i * 64], iQuant, iDcScaler);
369 : Isibaar 1.1 }
370 :     stop_iquant_timer();
371 :    
372 :     start_timer();
373 : edgomez 1.19 idct(&data[i * 64]);
374 : Isibaar 1.1 stop_idct_timer();
375 : suxen_drol 1.37.2.15
376 : h 1.5 }
377 : Isibaar 1.1
378 : edgomez 1.19 if (dec->interlacing && pMB->field_dct) {
379 : h 1.8 next_block = stride;
380 :     stride *= 2;
381 : Isibaar 1.1 }
382 : h 1.5
383 :     start_timer();
384 : suxen_drol 1.37.2.15
385 :     if (reduced_resolution)
386 :     {
387 :     next_block*=2;
388 :     copy_upsampled_8x8_16to8(pY_Cur, &data[0 * 64], stride);
389 :     copy_upsampled_8x8_16to8(pY_Cur + 16, &data[1 * 64], stride);
390 :     copy_upsampled_8x8_16to8(pY_Cur + next_block, &data[2 * 64], stride);
391 :     copy_upsampled_8x8_16to8(pY_Cur + 16 + next_block, &data[3 * 64], stride);
392 :     copy_upsampled_8x8_16to8(pU_Cur, &data[4 * 64], stride2);
393 :     copy_upsampled_8x8_16to8(pV_Cur, &data[5 * 64], stride2);
394 :     }else{
395 :     transfer_16to8copy(pY_Cur, &data[0 * 64], stride);
396 :     transfer_16to8copy(pY_Cur + 8, &data[1 * 64], stride);
397 :     transfer_16to8copy(pY_Cur + next_block, &data[2 * 64], stride);
398 :     transfer_16to8copy(pY_Cur + 8 + next_block, &data[3 * 64], stride);
399 :     transfer_16to8copy(pU_Cur, &data[4 * 64], stride2);
400 :     transfer_16to8copy(pV_Cur, &data[5 * 64], stride2);
401 :     }
402 : h 1.5 stop_transfer_timer();
403 : Isibaar 1.1 }
404 :    
405 :    
406 :    
407 :    
408 :     // decode an inter macroblock
409 : suxen_drol 1.37.2.15
410 : edgomez 1.19 void
411 :     decoder_mbinter(DECODER * dec,
412 :     const MACROBLOCK * pMB,
413 :     const uint32_t x_pos,
414 :     const uint32_t y_pos,
415 : chl 1.37.2.29 const uint32_t fcode,
416 : edgomez 1.19 const uint32_t cbp,
417 :     Bitstream * bs,
418 :     const uint32_t quant,
419 : suxen_drol 1.37.2.15 const uint32_t rounding,
420 : chl 1.37.2.28 const int reduced_resolution)
421 : Isibaar 1.1 {
422 : edgomez 1.7
423 : edgomez 1.19 DECLARE_ALIGNED_MATRIX(block, 6, 64, int16_t, CACHE_LINE);
424 : edgomez 1.7 DECLARE_ALIGNED_MATRIX(data, 6, 64, int16_t, CACHE_LINE);
425 : h 1.5
426 : h 1.8 uint32_t stride = dec->edged_width;
427 :     uint32_t stride2 = stride / 2;
428 : suxen_drol 1.37.2.15 uint32_t next_block = stride * (reduced_resolution ? 16 : 8);
429 : edgomez 1.7 uint32_t i;
430 :     uint32_t iQuant = pMB->quant;
431 : h 1.5 uint8_t *pY_Cur, *pU_Cur, *pV_Cur;
432 : suxen_drol 1.37.2.15
433 : Isibaar 1.1 int uv_dx, uv_dy;
434 : suxen_drol 1.37.2.18 VECTOR mv[4]; /* local copy of mvs */
435 : suxen_drol 1.37.2.15
436 :     if (reduced_resolution) {
437 :     pY_Cur = dec->cur.y + (y_pos << 5) * stride + (x_pos << 5);
438 :     pU_Cur = dec->cur.u + (y_pos << 4) * stride2 + (x_pos << 4);
439 :     pV_Cur = dec->cur.v + (y_pos << 4) * stride2 + (x_pos << 4);
440 : suxen_drol 1.37.2.18 for (i = 0; i < 4; i++) {
441 :     mv[i].x = RRV_MV_SCALEUP(pMB->mvs[i].x);
442 :     mv[i].y = RRV_MV_SCALEUP(pMB->mvs[i].y);
443 :     }
444 : chl 1.37.2.27 } else {
445 : suxen_drol 1.37.2.15 pY_Cur = dec->cur.y + (y_pos << 4) * stride + (x_pos << 4);
446 :     pU_Cur = dec->cur.u + (y_pos << 3) * stride2 + (x_pos << 3);
447 :     pV_Cur = dec->cur.v + (y_pos << 3) * stride2 + (x_pos << 3);
448 : suxen_drol 1.37.2.18 for (i = 0; i < 4; i++)
449 :     mv[i] = pMB->mvs[i];
450 : suxen_drol 1.37.2.15 }
451 : chl 1.37.2.27
452 : edgomez 1.19 if (pMB->mode == MODE_INTER || pMB->mode == MODE_INTER_Q) {
453 : Isibaar 1.37.2.7
454 : chl 1.37.2.27 uv_dx = mv[0].x / (1 + dec->quarterpel);
455 :     uv_dy = mv[0].y / (1 + dec->quarterpel);
456 :    
457 : Isibaar 1.37.2.7 uv_dx = (uv_dx >> 1) + roundtab_79[uv_dx & 0x3];
458 :     uv_dy = (uv_dy >> 1) + roundtab_79[uv_dy & 0x3];
459 : Isibaar 1.37.2.3
460 :     start_timer();
461 : suxen_drol 1.37.2.15 if (reduced_resolution)
462 :     {
463 : chl 1.37.2.28 interpolate32x32_switch(dec->cur.y, dec->refn[0].y, 32*x_pos, 32*y_pos,
464 : suxen_drol 1.37.2.15 mv[0].x, mv[0].y, stride, rounding);
465 : chl 1.37.2.28 interpolate16x16_switch(dec->cur.u, dec->refn[0].u, 16 * x_pos, 16 * y_pos,
466 : suxen_drol 1.37.2.15 uv_dx, uv_dy, stride2, rounding);
467 : chl 1.37.2.28 interpolate16x16_switch(dec->cur.v, dec->refn[0].v, 16 * x_pos, 16 * y_pos,
468 : suxen_drol 1.37.2.15 uv_dx, uv_dy, stride2, rounding);
469 :    
470 : Isibaar 1.37.2.3 }
471 : suxen_drol 1.37.2.15 else
472 :     {
473 :     if(dec->quarterpel) {
474 : chl 1.37.2.28 interpolate16x16_quarterpel(dec->cur.y, dec->refn[0].y, dec->qtmp.y, dec->qtmp.y + 64,
475 : suxen_drol 1.37.2.21 dec->qtmp.y + 128, 16*x_pos, 16*y_pos,
476 : chl 1.37.2.27 mv[0].x, mv[0].y, stride, rounding);
477 : suxen_drol 1.37.2.15 }
478 :     else {
479 : chl 1.37.2.28 interpolate16x16_switch(dec->cur.y, dec->refn[0].y, 16*x_pos, 16*y_pos,
480 : chl 1.37.2.27 mv[0].x, mv[0].y, stride, rounding);
481 : suxen_drol 1.37.2.15 }
482 : Isibaar 1.37.2.3
483 : chl 1.37.2.28 interpolate8x8_switch(dec->cur.u, dec->refn[0].u, 8 * x_pos, 8 * y_pos,
484 : suxen_drol 1.37.2.15 uv_dx, uv_dy, stride2, rounding);
485 : chl 1.37.2.28 interpolate8x8_switch(dec->cur.v, dec->refn[0].v, 8 * x_pos, 8 * y_pos,
486 : suxen_drol 1.37.2.15 uv_dx, uv_dy, stride2, rounding);
487 :     }
488 : Isibaar 1.37.2.3 stop_comp_timer();
489 :    
490 : suxen_drol 1.37.2.15 } else { /* MODE_INTER4V */
491 : Isibaar 1.1 int sum;
492 : Isibaar 1.37.2.6
493 :     if(dec->quarterpel)
494 : suxen_drol 1.37.2.15 sum = (mv[0].x / 2) + (mv[1].x / 2) + (mv[2].x / 2) + (mv[3].x / 2);
495 : Isibaar 1.37.2.6 else
496 : suxen_drol 1.37.2.15 sum = mv[0].x + mv[1].x + mv[2].x + mv[3].x;
497 : Isibaar 1.37.2.6
498 :     uv_dx = (sum >> 3) + roundtab_76[sum & 0xf];
499 :    
500 :     if(dec->quarterpel)
501 : suxen_drol 1.37.2.15 sum = (mv[0].y / 2) + (mv[1].y / 2) + (mv[2].y / 2) + (mv[3].y / 2);
502 : Isibaar 1.37.2.6 else
503 : suxen_drol 1.37.2.15 sum = mv[0].y + mv[1].y + mv[2].y + mv[3].y;
504 : Isibaar 1.34
505 : Isibaar 1.37.2.6 uv_dy = (sum >> 3) + roundtab_76[sum & 0xf];
506 : Isibaar 1.1
507 : Isibaar 1.37.2.3 start_timer();
508 : suxen_drol 1.37.2.15 if (reduced_resolution)
509 :     {
510 : chl 1.37.2.28 interpolate16x16_switch(dec->cur.y, dec->refn[0].y, 32*x_pos, 32*y_pos,
511 : suxen_drol 1.37.2.15 mv[0].x, mv[0].y, stride, rounding);
512 : chl 1.37.2.28 interpolate16x16_switch(dec->cur.y, dec->refn[0].y , 32*x_pos + 16, 32*y_pos,
513 : suxen_drol 1.37.2.15 mv[1].x, mv[1].y, stride, rounding);
514 : chl 1.37.2.28 interpolate16x16_switch(dec->cur.y, dec->refn[0].y , 32*x_pos, 32*y_pos + 16,
515 : suxen_drol 1.37.2.15 mv[2].x, mv[2].y, stride, rounding);
516 : chl 1.37.2.28 interpolate16x16_switch(dec->cur.y, dec->refn[0].y , 32*x_pos + 16, 32*y_pos + 16,
517 : suxen_drol 1.37.2.15 mv[3].x, mv[3].y, stride, rounding);
518 : chl 1.37.2.28 interpolate16x16_switch(dec->cur.u, dec->refn[0].u , 16 * x_pos, 16 * y_pos,
519 : suxen_drol 1.37.2.15 uv_dx, uv_dy, stride2, rounding);
520 : chl 1.37.2.28 interpolate16x16_switch(dec->cur.v, dec->refn[0].v , 16 * x_pos, 16 * y_pos,
521 : suxen_drol 1.37.2.15 uv_dx, uv_dy, stride2, rounding);
522 : Isibaar 1.37.2.3
523 : suxen_drol 1.37.2.15 // set_block(pY_Cur, stride, 32, 32, 127);
524 :     }
525 :     else
526 :     {
527 :     if(dec->quarterpel) {
528 : chl 1.37.2.28 interpolate8x8_quarterpel(dec->cur.y, dec->refn[0].y , dec->qtmp.y, dec->qtmp.y + 64,
529 : suxen_drol 1.37.2.21 dec->qtmp.y + 128, 16*x_pos, 16*y_pos,
530 : suxen_drol 1.37.2.15 mv[0].x, mv[0].y, stride, rounding);
531 : chl 1.37.2.28 interpolate8x8_quarterpel(dec->cur.y, dec->refn[0].y , dec->qtmp.y, dec->qtmp.y + 64,
532 : suxen_drol 1.37.2.21 dec->qtmp.y + 128, 16*x_pos + 8, 16*y_pos,
533 : suxen_drol 1.37.2.15 mv[1].x, mv[1].y, stride, rounding);
534 : chl 1.37.2.28 interpolate8x8_quarterpel(dec->cur.y, dec->refn[0].y , dec->qtmp.y, dec->qtmp.y + 64,
535 : suxen_drol 1.37.2.21 dec->qtmp.y + 128, 16*x_pos, 16*y_pos + 8,
536 : suxen_drol 1.37.2.15 mv[2].x, mv[2].y, stride, rounding);
537 : chl 1.37.2.28 interpolate8x8_quarterpel(dec->cur.y, dec->refn[0].y , dec->qtmp.y, dec->qtmp.y + 64,
538 : suxen_drol 1.37.2.21 dec->qtmp.y + 128, 16*x_pos + 8, 16*y_pos + 8,
539 : suxen_drol 1.37.2.15 mv[3].x, mv[3].y, stride, rounding);
540 :     }
541 :     else {
542 : chl 1.37.2.28 interpolate8x8_switch(dec->cur.y, dec->refn[0].y , 16*x_pos, 16*y_pos,
543 : suxen_drol 1.37.2.15 mv[0].x, mv[0].y, stride, rounding);
544 : chl 1.37.2.28 interpolate8x8_switch(dec->cur.y, dec->refn[0].y , 16*x_pos + 8, 16*y_pos,
545 : suxen_drol 1.37.2.15 mv[1].x, mv[1].y, stride, rounding);
546 : chl 1.37.2.28 interpolate8x8_switch(dec->cur.y, dec->refn[0].y , 16*x_pos, 16*y_pos + 8,
547 : suxen_drol 1.37.2.15 mv[2].x, mv[2].y, stride, rounding);
548 : chl 1.37.2.28 interpolate8x8_switch(dec->cur.y, dec->refn[0].y , 16*x_pos + 8, 16*y_pos + 8,
549 : suxen_drol 1.37.2.15 mv[3].x, mv[3].y, stride, rounding);
550 :     }
551 :    
552 : chl 1.37.2.28 interpolate8x8_switch(dec->cur.u, dec->refn[0].u , 8 * x_pos, 8 * y_pos,
553 : suxen_drol 1.37.2.15 uv_dx, uv_dy, stride2, rounding);
554 : chl 1.37.2.28 interpolate8x8_switch(dec->cur.v, dec->refn[0].v , 8 * x_pos, 8 * y_pos,
555 : suxen_drol 1.37.2.15 uv_dx, uv_dy, stride2, rounding);
556 :     }
557 : Isibaar 1.37.2.3 stop_comp_timer();
558 : Isibaar 1.34 }
559 : Isibaar 1.1
560 : edgomez 1.19 for (i = 0; i < 6; i++) {
561 : h 1.37.2.1 int direction = dec->alternate_vertical_scan ? 2 : 0;
562 :    
563 : edgomez 1.19 if (cbp & (1 << (5 - i))) // coded
564 : Isibaar 1.1 {
565 : edgomez 1.19 memset(&block[i * 64], 0, 64 * sizeof(int16_t)); // clear
566 : Isibaar 1.1
567 :     start_timer();
568 : h 1.37.2.1 get_inter_block(bs, &block[i * 64], direction);
569 : Isibaar 1.1 stop_coding_timer();
570 :    
571 :     start_timer();
572 : edgomez 1.19 if (dec->quant_type == 0) {
573 :     dequant_inter(&data[i * 64], &block[i * 64], iQuant);
574 :     } else {
575 :     dequant4_inter(&data[i * 64], &block[i * 64], iQuant);
576 : Isibaar 1.1 }
577 :     stop_iquant_timer();
578 :    
579 :     start_timer();
580 : edgomez 1.19 idct(&data[i * 64]);
581 : Isibaar 1.1 stop_idct_timer();
582 : h 1.5 }
583 :     }
584 : Isibaar 1.1
585 : edgomez 1.19 if (dec->interlacing && pMB->field_dct) {
586 : h 1.8 next_block = stride;
587 :     stride *= 2;
588 : Isibaar 1.1 }
589 : h 1.5
590 :     start_timer();
591 : suxen_drol 1.37.2.15 if (reduced_resolution)
592 :     {
593 :     if (cbp & 32)
594 :     add_upsampled_8x8_16to8(pY_Cur, &data[0 * 64], stride);
595 :     if (cbp & 16)
596 :     add_upsampled_8x8_16to8(pY_Cur + 16, &data[1 * 64], stride);
597 :     if (cbp & 8)
598 :     add_upsampled_8x8_16to8(pY_Cur + next_block, &data[2 * 64], stride);
599 :     if (cbp & 4)
600 :     add_upsampled_8x8_16to8(pY_Cur + 16 + next_block, &data[3 * 64], stride);
601 :     if (cbp & 2)
602 :     add_upsampled_8x8_16to8(pU_Cur, &data[4 * 64], stride2);
603 :     if (cbp & 1)
604 :     add_upsampled_8x8_16to8(pV_Cur, &data[5 * 64], stride2);
605 :     }
606 :     else
607 :     {
608 :     if (cbp & 32)
609 :     transfer_16to8add(pY_Cur, &data[0 * 64], stride);
610 :     if (cbp & 16)
611 :     transfer_16to8add(pY_Cur + 8, &data[1 * 64], stride);
612 :     if (cbp & 8)
613 :     transfer_16to8add(pY_Cur + next_block, &data[2 * 64], stride);
614 :     if (cbp & 4)
615 :     transfer_16to8add(pY_Cur + 8 + next_block, &data[3 * 64], stride);
616 :     if (cbp & 2)
617 :     transfer_16to8add(pU_Cur, &data[4 * 64], stride2);
618 :     if (cbp & 1)
619 :     transfer_16to8add(pV_Cur, &data[5 * 64], stride2);
620 :     }
621 : h 1.5 stop_transfer_timer();
622 : Isibaar 1.1 }
623 :    
624 : chl 1.37.2.29 static __inline int gmc_sanitize(int value, int quarterpel, int fcode)
625 :     {
626 :     int length = 1 << (fcode+4);
627 :    
628 : chl 1.37.2.30 // if (quarterpel) value *= 2;
629 : chl 1.37.2.29
630 :     if (value < -length)
631 :     return -length;
632 :     else if (value >= length)
633 :     return length-1;
634 :     else return value;
635 :     }
636 :    
637 : Isibaar 1.1
638 : chl 1.37.2.28 static void
639 :     decoder_mbgmc(DECODER * dec,
640 :     MACROBLOCK * const pMB,
641 :     const uint32_t x_pos,
642 :     const uint32_t y_pos,
643 : chl 1.37.2.29 const uint32_t fcode,
644 : chl 1.37.2.28 const uint32_t cbp,
645 :     Bitstream * bs,
646 :     const uint32_t quant,
647 :     const uint32_t rounding,
648 :     const int reduced_resolution) /* no reduced res support */
649 :     {
650 :    
651 :     DECLARE_ALIGNED_MATRIX(block, 6, 64, int16_t, CACHE_LINE);
652 :     DECLARE_ALIGNED_MATRIX(data, 6, 64, int16_t, CACHE_LINE);
653 :    
654 :     const uint32_t stride = dec->edged_width;
655 :     const uint32_t stride2 = stride / 2;
656 :     const uint32_t next_block = stride * (reduced_resolution ? 16 : 8);
657 :     uint32_t i;
658 :     const uint32_t iQuant = pMB->quant;
659 :     uint8_t *const pY_Cur=dec->cur.y + (y_pos << 4) * stride + (x_pos << 4);
660 :     uint8_t *const pU_Cur=dec->cur.u + (y_pos << 3) * stride2 + (x_pos << 3);
661 :     uint8_t *const pV_Cur=dec->cur.v + (y_pos << 3) * stride2 + (x_pos << 3);
662 :    
663 :     pMB->mvs[0] = pMB->mvs[1] = pMB->mvs[2] = pMB->mvs[3] = pMB->amv;
664 :    
665 :     start_timer();
666 : chl 1.37.2.29
667 :     /* this is where the calculations are done */
668 :    
669 :     {
670 :     pMB->amv = generate_GMCimageMB(&dec->gmc_data, &dec->refn[0], x_pos, y_pos,
671 :     stride, stride2, dec->quarterpel, rounding, &dec->cur);
672 :    
673 :     pMB->amv.x = gmc_sanitize(pMB->amv.x, dec->quarterpel, fcode);
674 :     pMB->amv.y = gmc_sanitize(pMB->amv.y, dec->quarterpel, fcode);
675 :     }
676 :     pMB->mvs[0] = pMB->mvs[1] = pMB->mvs[2] = pMB->mvs[3] = pMB->amv;
677 :    
678 :     /* transfer16x16_copy(pY_Cur, dec->gmc.y + (y_pos << 4)*stride + (x_pos << 4), stride);
679 : chl 1.37.2.28 transfer8x8_copy(pU_Cur, dec->gmc.u + (y_pos << 3)*stride2 + (x_pos << 3), stride2);
680 :     transfer8x8_copy(pV_Cur, dec->gmc.v + (y_pos << 3)*stride2 + (x_pos << 3), stride2);
681 : chl 1.37.2.29 */
682 :    
683 :    
684 : chl 1.37.2.28 stop_transfer_timer();
685 :    
686 :     if (!cbp) return;
687 :    
688 :     for (i = 0; i < 6; i++) {
689 :     int direction = dec->alternate_vertical_scan ? 2 : 0;
690 :    
691 :     if (cbp & (1 << (5 - i))) // coded
692 :     {
693 :     memset(&block[i * 64], 0, 64 * sizeof(int16_t)); // clear
694 :    
695 :     start_timer();
696 :     get_inter_block(bs, &block[i * 64], direction);
697 :     stop_coding_timer();
698 :    
699 :     start_timer();
700 :     if (dec->quant_type == 0) {
701 :     dequant_inter(&data[i * 64], &block[i * 64], iQuant);
702 :     } else {
703 :     dequant4_inter(&data[i * 64], &block[i * 64], iQuant);
704 :     }
705 :     stop_iquant_timer();
706 :    
707 :     start_timer();
708 :     idct(&data[i * 64]);
709 :     stop_idct_timer();
710 :     }
711 :     }
712 :    
713 :     /* interlace + GMC is this possible ??? */
714 :     /* if (dec->interlacing && pMB->field_dct) {
715 :     next_block = stride;
716 :     stride *= 2;
717 :     }
718 :     */
719 :     start_timer();
720 :     if (cbp & 32)
721 :     transfer_16to8add(pY_Cur, &data[0 * 64], stride);
722 :     if (cbp & 16)
723 :     transfer_16to8add(pY_Cur + 8, &data[1 * 64], stride);
724 :     if (cbp & 8)
725 :     transfer_16to8add(pY_Cur + next_block, &data[2 * 64], stride);
726 :     if (cbp & 4)
727 :     transfer_16to8add(pY_Cur + 8 + next_block, &data[3 * 64], stride);
728 :     if (cbp & 2)
729 :     transfer_16to8add(pU_Cur, &data[4 * 64], stride2);
730 :     if (cbp & 1)
731 :     transfer_16to8add(pV_Cur, &data[5 * 64], stride2);
732 :     stop_transfer_timer();
733 :     }
734 :    
735 :    
736 : edgomez 1.19 void
737 :     decoder_iframe(DECODER * dec,
738 :     Bitstream * bs,
739 : suxen_drol 1.37.2.8 int reduced_resolution,
740 : edgomez 1.19 int quant,
741 :     int intra_dc_threshold)
742 : Isibaar 1.1 {
743 : chenm001 1.26 uint32_t bound;
744 :     uint32_t x, y;
745 : suxen_drol 1.37.2.21 uint32_t mb_width = dec->mb_width;
746 :     uint32_t mb_height = dec->mb_height;
747 : suxen_drol 1.37.2.8
748 :     if (reduced_resolution)
749 :     {
750 : suxen_drol 1.37.2.17 mb_width = (dec->width + 31) / 32;
751 :     mb_height = (dec->height + 31) / 32;
752 : suxen_drol 1.37.2.8 }
753 : chenm001 1.25
754 : chenm001 1.26 bound = 0;
755 : edgomez 1.7
756 : suxen_drol 1.37.2.8 for (y = 0; y < mb_height; y++) {
757 :     for (x = 0; x < mb_width; x++) {
758 : chenm001 1.26 MACROBLOCK *mb;
759 : Isibaar 1.1 uint32_t mcbpc;
760 :     uint32_t cbpc;
761 :     uint32_t acpred_flag;
762 :     uint32_t cbpy;
763 :     uint32_t cbp;
764 :    
765 : chenm001 1.26 while (BitstreamShowBits(bs, 9) == 1)
766 :     BitstreamSkip(bs, 9);
767 :    
768 :     if (check_resync_marker(bs, 0))
769 :     {
770 : suxen_drol 1.37.2.8 bound = read_video_packet_header(bs, dec, 0,
771 :     &quant, NULL, NULL, &intra_dc_threshold);
772 : suxen_drol 1.37.2.15 x = bound % mb_width;
773 :     y = bound / mb_width;
774 : chenm001 1.26 }
775 :     mb = &dec->mbs[y * dec->mb_width + x];
776 :    
777 :     DPRINTF(DPRINTF_MB, "macroblock (%i,%i) %08x", x, y, BitstreamShowBits(bs, 32));
778 :    
779 : Isibaar 1.1 mcbpc = get_mcbpc_intra(bs);
780 :     mb->mode = mcbpc & 7;
781 :     cbpc = (mcbpc >> 4);
782 :    
783 :     acpred_flag = BitstreamGetBit(bs);
784 :    
785 :     cbpy = get_cbpy(bs, 1);
786 :     cbp = (cbpy << 2) | cbpc;
787 :    
788 : edgomez 1.19 if (mb->mode == MODE_INTRA_Q) {
789 :     quant += dquant_table[BitstreamGetBits(bs, 2)];
790 :     if (quant > 31) {
791 : Isibaar 1.1 quant = 31;
792 : edgomez 1.19 } else if (quant < 1) {
793 : Isibaar 1.1 quant = 1;
794 :     }
795 :     }
796 :     mb->quant = quant;
797 : chenm001 1.30 mb->mvs[0].x = mb->mvs[0].y =
798 :     mb->mvs[1].x = mb->mvs[1].y =
799 :     mb->mvs[2].x = mb->mvs[2].y =
800 :     mb->mvs[3].x = mb->mvs[3].y =0;
801 : h 1.5
802 : edgomez 1.19 if (dec->interlacing) {
803 : h 1.5 mb->field_dct = BitstreamGetBit(bs);
804 : suxen_drol 1.37.2.22 DPRINTF(DPRINTF_MB,"deci: field_dct: %i", mb->field_dct);
805 : h 1.5 }
806 : Isibaar 1.1
807 : edgomez 1.19 decoder_mbintra(dec, mb, x, y, acpred_flag, cbp, bs, quant,
808 : suxen_drol 1.37.2.15 intra_dc_threshold, bound, reduced_resolution);
809 : suxen_drol 1.37.2.8
810 : Isibaar 1.1 }
811 : albeu 1.32 if(dec->out_frm)
812 : suxen_drol 1.37.2.15 output_slice(&dec->cur, dec->edged_width,dec->width,dec->out_frm,0,y,mb_width);
813 : Isibaar 1.1 }
814 : edgomez 1.7
815 : Isibaar 1.1 }
816 :    
817 :    
818 : edgomez 1.19 void
819 :     get_motion_vector(DECODER * dec,
820 :     Bitstream * bs,
821 :     int x,
822 :     int y,
823 :     int k,
824 : suxen_drol 1.37.2.15 VECTOR * ret_mv,
825 : chenm001 1.26 int fcode,
826 :     const int bound)
827 : Isibaar 1.1 {
828 : edgomez 1.7
829 : Isibaar 1.1 int scale_fac = 1 << (fcode - 1);
830 :     int high = (32 * scale_fac) - 1;
831 :     int low = ((-32) * scale_fac);
832 :     int range = (64 * scale_fac);
833 :    
834 : chenm001 1.26 VECTOR pmv;
835 : suxen_drol 1.37.2.15 VECTOR mv;
836 : Isibaar 1.1
837 : chenm001 1.26 pmv = get_pmv2(dec->mbs, dec->mb_width, bound, x, y, k);
838 : Isibaar 1.1
839 : suxen_drol 1.37.2.15 mv.x = get_mv(bs, fcode);
840 :     mv.y = get_mv(bs, fcode);
841 : edgomez 1.19
842 : chl 1.37.2.27 DPRINTF(DPRINTF_MV,"mv_diff (%i,%i) pred (%i,%i) result (%i,%i)", mv.x, mv.y, pmv.x, pmv.y, mv.x+pmv.x, mv.y+pmv.y);
843 : chenm001 1.26
844 : suxen_drol 1.37.2.15 mv.x += pmv.x;
845 :     mv.y += pmv.y;
846 : Isibaar 1.1
847 : suxen_drol 1.37.2.15 if (mv.x < low) {
848 :     mv.x += range;
849 :     } else if (mv.x > high) {
850 :     mv.x -= range;
851 : Isibaar 1.1 }
852 :    
853 : suxen_drol 1.37.2.15 if (mv.y < low) {
854 :     mv.y += range;
855 :     } else if (mv.y > high) {
856 :     mv.y -= range;
857 : Isibaar 1.1 }
858 :    
859 : suxen_drol 1.37.2.15 ret_mv->x = mv.x;
860 :     ret_mv->y = mv.y;
861 : suxen_drol 1.37.2.8 }
862 : Isibaar 1.1
863 : suxen_drol 1.37.2.8
864 :    
865 : Isibaar 1.1
866 :    
867 : chl 1.37.2.27 /* for P_VOP set gmc_warp to NULL */
868 : edgomez 1.19 void
869 :     decoder_pframe(DECODER * dec,
870 :     Bitstream * bs,
871 :     int rounding,
872 : suxen_drol 1.37.2.8 int reduced_resolution,
873 : edgomez 1.19 int quant,
874 :     int fcode,
875 : suxen_drol 1.37.2.8 int intra_dc_threshold,
876 : chl 1.37.2.27 const WARPPOINTS *const gmc_warp)
877 : Isibaar 1.1 {
878 : edgomez 1.7
879 : Isibaar 1.1 uint32_t x, y;
880 : chenm001 1.26 uint32_t bound;
881 : albeu 1.32 int cp_mb, st_mb;
882 : suxen_drol 1.37.2.21 uint32_t mb_width = dec->mb_width;
883 :     uint32_t mb_height = dec->mb_height;
884 : suxen_drol 1.37.2.15
885 :     if (reduced_resolution)
886 :     {
887 : suxen_drol 1.37.2.17 mb_width = (dec->width + 31) / 32;
888 :     mb_height = (dec->height + 31) / 32;
889 : suxen_drol 1.37.2.15 }
890 : Isibaar 1.1
891 :     start_timer();
892 : edgomez 1.19 image_setedges(&dec->refn[0], dec->edged_width, dec->edged_height,
893 : h 1.37.2.1 dec->width, dec->height);
894 : Isibaar 1.1 stop_edges_timer();
895 :    
896 : chl 1.37.2.27 if (gmc_warp)
897 :     {
898 : chl 1.37.2.28
899 : chl 1.37.2.27 // accuracy: 0==1/2, 1=1/4, 2=1/8, 3=1/16
900 :     if ( (dec->sprite_warping_accuracy != 3) || (dec->sprite_warping_points != 2) )
901 :     {
902 :     fprintf(stderr,"Wrong GMC parameters acc=%d(-> 1/%d), %d!!!\n",
903 :     dec->sprite_warping_accuracy,(2<<dec->sprite_warping_accuracy),
904 :     dec->sprite_warping_points);
905 :     }
906 :    
907 :     generate_GMCparameters( dec->sprite_warping_points,
908 :     (2 << dec->sprite_warping_accuracy), gmc_warp,
909 :     dec->width, dec->height, &dec->gmc_data);
910 :    
911 : chl 1.37.2.29 /* image warping is done block-based in decoder_mbgmc(), now */
912 :     /*
913 :     generate_GMCimage(&dec->gmc_data, &dec->refn[0],
914 : chl 1.37.2.27 mb_width, mb_height,
915 :     dec->edged_width, dec->edged_width/2,
916 : chl 1.37.2.28 fcode, dec->quarterpel, 0,
917 : chl 1.37.2.27 rounding, dec->mbs, &dec->gmc);
918 : chl 1.37.2.29 */
919 : chl 1.37.2.27 }
920 :    
921 : chenm001 1.26 bound = 0;
922 :    
923 : suxen_drol 1.37.2.15 for (y = 0; y < mb_height; y++) {
924 : albeu 1.32 cp_mb = st_mb = 0;
925 : suxen_drol 1.37.2.15 for (x = 0; x < mb_width; x++) {
926 : chenm001 1.26 MACROBLOCK *mb;
927 :    
928 :     // skip stuffing
929 :     while (BitstreamShowBits(bs, 10) == 1)
930 :     BitstreamSkip(bs, 10);
931 :    
932 :     if (check_resync_marker(bs, fcode - 1))
933 :     {
934 : suxen_drol 1.37.2.8 bound = read_video_packet_header(bs, dec, fcode - 1,
935 :     &quant, &fcode, NULL, &intra_dc_threshold);
936 : suxen_drol 1.37.2.15 x = bound % mb_width;
937 :     y = bound / mb_width;
938 : chenm001 1.26 }
939 :     mb = &dec->mbs[y * dec->mb_width + x];
940 :    
941 :     DPRINTF(DPRINTF_MB, "macroblock (%i,%i) %08x", x, y, BitstreamShowBits(bs, 32));
942 : suxen_drol 1.23
943 : edgomez 1.19 //if (!(dec->mb_skip[y*dec->mb_width + x]=BitstreamGetBit(bs))) // not_coded
944 : chl 1.37.2.27 if (!(BitstreamGetBit(bs))) // block _is_ coded
945 : Isibaar 1.1 {
946 :     uint32_t mcbpc;
947 :     uint32_t cbpc;
948 :     uint32_t acpred_flag;
949 :     uint32_t cbpy;
950 :     uint32_t cbp;
951 :     uint32_t intra;
952 : suxen_drol 1.37.2.8 int mcsel = 0; // mcsel: '0'=local motion, '1'=GMC
953 : Isibaar 1.1
954 : albeu 1.32 cp_mb++;
955 : Isibaar 1.1 mcbpc = get_mcbpc_inter(bs);
956 :     mb->mode = mcbpc & 7;
957 :     cbpc = (mcbpc >> 4);
958 : chenm001 1.26
959 :     DPRINTF(DPRINTF_MB, "mode %i", mb->mode);
960 :     DPRINTF(DPRINTF_MB, "cbpc %i", cbpc);
961 : edgomez 1.2 acpred_flag = 0;
962 : Isibaar 1.1
963 :     intra = (mb->mode == MODE_INTRA || mb->mode == MODE_INTRA_Q);
964 : edgomez 1.19
965 :     if (intra) {
966 : Isibaar 1.1 acpred_flag = BitstreamGetBit(bs);
967 :     }
968 :    
969 : chl 1.37.2.27 if (gmc_warp && (mb->mode == MODE_INTER || mb->mode == MODE_INTER_Q))
970 : suxen_drol 1.37.2.8 {
971 :     mcsel = BitstreamGetBit(bs);
972 :     }
973 :    
974 : chenm001 1.26 cbpy = get_cbpy(bs, intra);
975 : chl 1.37.2.27 DPRINTF(DPRINTF_MB, "cbpy %i mcsel %i ", cbpy,mcsel);
976 : chenm001 1.25
977 : Isibaar 1.1 cbp = (cbpy << 2) | cbpc;
978 :    
979 : edgomez 1.19 if (mb->mode == MODE_INTER_Q || mb->mode == MODE_INTRA_Q) {
980 : chenm001 1.26 int dquant = dquant_table[BitstreamGetBits(bs, 2)];
981 :     DPRINTF(DPRINTF_MB, "dquant %i", dquant);
982 :     quant += dquant;
983 : edgomez 1.19 if (quant > 31) {
984 : Isibaar 1.1 quant = 31;
985 : chenm001 1.26 } else if (quant < 1) {
986 : Isibaar 1.1 quant = 1;
987 :     }
988 : chenm001 1.26 DPRINTF(DPRINTF_MB, "quant %i", quant);
989 : Isibaar 1.1 }
990 :     mb->quant = quant;
991 : h 1.5
992 : edgomez 1.19 if (dec->interlacing) {
993 : h 1.37 if (cbp || intra) {
994 :     mb->field_dct = BitstreamGetBit(bs);
995 : suxen_drol 1.37.2.22 DPRINTF(DPRINTF_MB,"decp: field_dct: %i", mb->field_dct);
996 : h 1.37 }
997 : h 1.5
998 : edgomez 1.19 if (mb->mode == MODE_INTER || mb->mode == MODE_INTER_Q) {
999 : h 1.5 mb->field_pred = BitstreamGetBit(bs);
1000 : suxen_drol 1.37.2.22 DPRINTF(DPRINTF_MB, "decp: field_pred: %i", mb->field_pred);
1001 : h 1.5
1002 : edgomez 1.19 if (mb->field_pred) {
1003 : h 1.5 mb->field_for_top = BitstreamGetBit(bs);
1004 : suxen_drol 1.37.2.22 DPRINTF(DPRINTF_MB,"decp: field_for_top: %i", mb->field_for_top);
1005 : h 1.5 mb->field_for_bot = BitstreamGetBit(bs);
1006 : suxen_drol 1.37.2.22 DPRINTF(DPRINTF_MB,"decp: field_for_bot: %i", mb->field_for_bot);
1007 : h 1.5 }
1008 :     }
1009 :     }
1010 : chl 1.37.2.28
1011 :     if (mcsel) {
1012 : chl 1.37.2.29 decoder_mbgmc(dec, mb, x, y, fcode, cbp, bs, quant,
1013 : chl 1.37.2.28 rounding, reduced_resolution);
1014 :     continue;
1015 : h 1.5
1016 : chl 1.37.2.28 } else if (mb->mode == MODE_INTER || mb->mode == MODE_INTER_Q) {
1017 : suxen_drol 1.37.2.8
1018 : chl 1.37.2.28 if (dec->interlacing && mb->field_pred) {
1019 : edgomez 1.19 get_motion_vector(dec, bs, x, y, 0, &mb->mvs[0],
1020 : chenm001 1.26 fcode, bound);
1021 : edgomez 1.19 get_motion_vector(dec, bs, x, y, 0, &mb->mvs[1],
1022 : chenm001 1.26 fcode, bound);
1023 : edgomez 1.19 } else {
1024 :     get_motion_vector(dec, bs, x, y, 0, &mb->mvs[0],
1025 : chenm001 1.26 fcode, bound);
1026 : chl 1.37.2.27 mb->mvs[1] = mb->mvs[2] = mb->mvs[3] = mb->mvs[0];
1027 : h 1.5 }
1028 : chl 1.35 } else if (mb->mode == MODE_INTER4V ) {
1029 :    
1030 : chenm001 1.26 get_motion_vector(dec, bs, x, y, 0, &mb->mvs[0], fcode, bound);
1031 :     get_motion_vector(dec, bs, x, y, 1, &mb->mvs[1], fcode, bound);
1032 :     get_motion_vector(dec, bs, x, y, 2, &mb->mvs[2], fcode, bound);
1033 :     get_motion_vector(dec, bs, x, y, 3, &mb->mvs[3], fcode, bound);
1034 : edgomez 1.19 } else // MODE_INTRA, MODE_INTRA_Q
1035 : Isibaar 1.1 {
1036 : edgomez 1.19 mb->mvs[0].x = mb->mvs[1].x = mb->mvs[2].x = mb->mvs[3].x =
1037 :     0;
1038 :     mb->mvs[0].y = mb->mvs[1].y = mb->mvs[2].y = mb->mvs[3].y =
1039 :     0;
1040 :     decoder_mbintra(dec, mb, x, y, acpred_flag, cbp, bs, quant,
1041 : suxen_drol 1.37.2.15 intra_dc_threshold, bound, reduced_resolution);
1042 : Isibaar 1.1 continue;
1043 :     }
1044 :    
1045 : chl 1.37.2.29 decoder_mbinter(dec, mb, x, y, fcode, cbp, bs, quant,
1046 : chl 1.37.2.28 rounding, reduced_resolution);
1047 : suxen_drol 1.37.2.8
1048 : chl 1.37.2.27 }
1049 :     else if (gmc_warp) /* a not coded S(GMC)-VOP macroblock */
1050 : suxen_drol 1.37.2.8 {
1051 : suxen_drol 1.37.2.24 mb->mode = MODE_NOT_CODED_GMC;
1052 : chl 1.37.2.27
1053 :     start_timer();
1054 :    
1055 : chl 1.37.2.29 decoder_mbgmc(dec, mb, x, y, fcode, 0x00, bs, quant,
1056 : chl 1.37.2.28 rounding, reduced_resolution);
1057 : chl 1.37.2.27
1058 :     stop_transfer_timer();
1059 :    
1060 :     if(dec->out_frm && cp_mb > 0) {
1061 :     output_slice(&dec->cur, dec->edged_width,dec->width,dec->out_frm,st_mb,y,cp_mb);
1062 :     cp_mb = 0;
1063 :     }
1064 :     st_mb = x+1;
1065 : suxen_drol 1.37.2.8 }
1066 :     else /* not coded P_VOP macroblock */
1067 : Isibaar 1.1 {
1068 : chenm001 1.14 mb->mode = MODE_NOT_CODED;
1069 : suxen_drol 1.37.2.15
1070 : Isibaar 1.1 mb->mvs[0].x = mb->mvs[1].x = mb->mvs[2].x = mb->mvs[3].x = 0;
1071 :     mb->mvs[0].y = mb->mvs[1].y = mb->mvs[2].y = mb->mvs[3].y = 0;
1072 :     // copy macroblock directly from ref to cur
1073 :    
1074 :     start_timer();
1075 :    
1076 : suxen_drol 1.37.2.15 if (reduced_resolution)
1077 :     {
1078 :     transfer32x32_copy(dec->cur.y + (32*y)*dec->edged_width + (32*x),
1079 :     dec->refn[0].y + (32*y)*dec->edged_width + (32*x),
1080 :     dec->edged_width);
1081 :    
1082 :     transfer16x16_copy(dec->cur.u + (16*y)*dec->edged_width/2 + (16*x),
1083 :     dec->refn[0].u + (16*y)*dec->edged_width/2 + (16*x),
1084 :     dec->edged_width/2);
1085 :    
1086 :     transfer16x16_copy(dec->cur.v + (16*y)*dec->edged_width/2 + (16*x),
1087 :     dec->refn[0].v + (16*y)*dec->edged_width/2 + (16*x),
1088 :     dec->edged_width/2);
1089 :     }
1090 :     else
1091 :     {
1092 :     transfer16x16_copy(dec->cur.y + (16*y)*dec->edged_width + (16*x),
1093 :     dec->refn[0].y + (16*y)*dec->edged_width + (16*x),
1094 :     dec->edged_width);
1095 :    
1096 :     transfer8x8_copy(dec->cur.u + (8*y)*dec->edged_width/2 + (8*x),
1097 :     dec->refn[0].u + (8*y)*dec->edged_width/2 + (8*x),
1098 :     dec->edged_width/2);
1099 :    
1100 :     transfer8x8_copy(dec->cur.v + (8*y)*dec->edged_width/2 + (8*x),
1101 :     dec->refn[0].v + (8*y)*dec->edged_width/2 + (8*x),
1102 :     dec->edged_width/2);
1103 :     }
1104 : suxen_drol 1.37.2.8
1105 : Isibaar 1.1 stop_transfer_timer();
1106 : suxen_drol 1.37.2.8
1107 : albeu 1.32 if(dec->out_frm && cp_mb > 0) {
1108 :     output_slice(&dec->cur, dec->edged_width,dec->width,dec->out_frm,st_mb,y,cp_mb);
1109 :     cp_mb = 0;
1110 :     }
1111 :     st_mb = x+1;
1112 : Isibaar 1.1 }
1113 :     }
1114 : albeu 1.32 if(dec->out_frm && cp_mb > 0)
1115 :     output_slice(&dec->cur, dec->edged_width,dec->width,dec->out_frm,st_mb,y,cp_mb);
1116 : Isibaar 1.1 }
1117 :     }
1118 :    
1119 : chenm001 1.14
1120 :     // add by MinChen <chenm001@163.com>
1121 :     // decode B-frame motion vector
1122 : edgomez 1.19 void
1123 :     get_b_motion_vector(DECODER * dec,
1124 :     Bitstream * bs,
1125 :     int x,
1126 :     int y,
1127 :     VECTOR * mv,
1128 :     int fcode,
1129 :     const VECTOR pmv)
1130 : chenm001 1.14 {
1131 :     int scale_fac = 1 << (fcode - 1);
1132 :     int high = (32 * scale_fac) - 1;
1133 :     int low = ((-32) * scale_fac);
1134 :     int range = (64 * scale_fac);
1135 :    
1136 :     int mv_x, mv_y;
1137 :     int pmv_x, pmv_y;
1138 :    
1139 :     pmv_x = pmv.x;
1140 :     pmv_y = pmv.y;
1141 :    
1142 :     mv_x = get_mv(bs, fcode);
1143 :     mv_y = get_mv(bs, fcode);
1144 : edgomez 1.19
1145 : chenm001 1.14 mv_x += pmv_x;
1146 :     mv_y += pmv_y;
1147 :    
1148 : edgomez 1.19 if (mv_x < low) {
1149 : chenm001 1.14 mv_x += range;
1150 : edgomez 1.19 } else if (mv_x > high) {
1151 : chenm001 1.14 mv_x -= range;
1152 :     }
1153 :    
1154 : edgomez 1.19 if (mv_y < low) {
1155 : chenm001 1.14 mv_y += range;
1156 : edgomez 1.19 } else if (mv_y > high) {
1157 : chenm001 1.14 mv_y -= range;
1158 :     }
1159 :    
1160 :     mv->x = mv_x;
1161 :     mv->y = mv_y;
1162 :     }
1163 :    
1164 :    
1165 :     // add by MinChen <chenm001@163.com>
1166 :     // decode an B-frame forward & backward inter macroblock
1167 : edgomez 1.19 void
1168 :     decoder_bf_mbinter(DECODER * dec,
1169 :     const MACROBLOCK * pMB,
1170 :     const uint32_t x_pos,
1171 :     const uint32_t y_pos,
1172 :     const uint32_t cbp,
1173 :     Bitstream * bs,
1174 :     const uint32_t quant,
1175 :     const uint8_t ref)
1176 : chenm001 1.14 {
1177 :    
1178 : edgomez 1.19 DECLARE_ALIGNED_MATRIX(block, 6, 64, int16_t, CACHE_LINE);
1179 : chenm001 1.14 DECLARE_ALIGNED_MATRIX(data, 6, 64, int16_t, CACHE_LINE);
1180 :    
1181 :     uint32_t stride = dec->edged_width;
1182 :     uint32_t stride2 = stride / 2;
1183 :     uint32_t next_block = stride * 8;
1184 :     uint32_t i;
1185 :     uint32_t iQuant = pMB->quant;
1186 :     uint8_t *pY_Cur, *pU_Cur, *pV_Cur;
1187 :     int uv_dx, uv_dy;
1188 :    
1189 :     pY_Cur = dec->cur.y + (y_pos << 4) * stride + (x_pos << 4);
1190 :     pU_Cur = dec->cur.u + (y_pos << 3) * stride2 + (x_pos << 3);
1191 :     pV_Cur = dec->cur.v + (y_pos << 3) * stride2 + (x_pos << 3);
1192 :    
1193 : chenm001 1.30
1194 : edgomez 1.19 if (!(pMB->mode == MODE_INTER || pMB->mode == MODE_INTER_Q)) {
1195 : chenm001 1.14 uv_dx = pMB->mvs[0].x;
1196 :     uv_dy = pMB->mvs[0].y;
1197 :    
1198 : Isibaar 1.37.2.13 if (dec->quarterpel)
1199 :     {
1200 :     uv_dx /= 2;
1201 :     uv_dy /= 2;
1202 :     }
1203 :    
1204 :     uv_dx = (uv_dx >> 1) + roundtab_79[uv_dx & 0x3];
1205 :     uv_dy = (uv_dy >> 1) + roundtab_79[uv_dy & 0x3];
1206 : edgomez 1.19 } else {
1207 : chenm001 1.14 int sum;
1208 : edgomez 1.19
1209 : Isibaar 1.37.2.13 if(dec->quarterpel)
1210 :     sum = (pMB->mvs[0].x / 2) + (pMB->mvs[1].x / 2) + (pMB->mvs[2].x / 2) + (pMB->mvs[3].x / 2);
1211 :     else
1212 :     sum = pMB->mvs[0].x + pMB->mvs[1].x + pMB->mvs[2].x + pMB->mvs[3].x;
1213 :    
1214 :     uv_dx = (sum >> 3) + roundtab_76[sum & 0xf];
1215 :    
1216 :     if(dec->quarterpel)
1217 :     sum = (pMB->mvs[0].y / 2) + (pMB->mvs[1].y / 2) + (pMB->mvs[2].y / 2) + (pMB->mvs[3].y / 2);
1218 :     else
1219 :     sum = pMB->mvs[0].y + pMB->mvs[1].y + pMB->mvs[2].y + pMB->mvs[3].y;
1220 :    
1221 :     uv_dy = (sum >> 3) + roundtab_76[sum & 0xf];
1222 : chenm001 1.14 }
1223 :    
1224 :     start_timer();
1225 : Isibaar 1.37.2.13 if(dec->quarterpel) {
1226 : suxen_drol 1.37.2.21 interpolate16x16_quarterpel(dec->cur.y, dec->refn[ref].y, dec->qtmp.y, dec->qtmp.y + 64,
1227 :     dec->qtmp.y + 128, 16*x_pos, 16*y_pos,
1228 : Isibaar 1.37.2.13 pMB->mvs[0].x, pMB->mvs[0].y, stride, 0);
1229 :     }
1230 :     else {
1231 :     interpolate8x8_switch(dec->cur.y, dec->refn[ref].y, 16*x_pos, 16*y_pos,
1232 :     pMB->mvs[0].x, pMB->mvs[0].y, stride, 0);
1233 :     interpolate8x8_switch(dec->cur.y, dec->refn[ref].y, 16*x_pos + 8, 16*y_pos,
1234 :     pMB->mvs[1].x, pMB->mvs[1].y, stride, 0);
1235 :     interpolate8x8_switch(dec->cur.y, dec->refn[ref].y, 16*x_pos, 16*y_pos + 8,
1236 :     pMB->mvs[2].x, pMB->mvs[2].y, stride, 0);
1237 :     interpolate8x8_switch(dec->cur.y, dec->refn[ref].y, 16*x_pos + 8, 16*y_pos + 8,
1238 :     pMB->mvs[3].x, pMB->mvs[3].y, stride, 0);
1239 :     }
1240 :    
1241 : edgomez 1.19 interpolate8x8_switch(dec->cur.u, dec->refn[ref].u, 8 * x_pos, 8 * y_pos,
1242 :     uv_dx, uv_dy, stride2, 0);
1243 :     interpolate8x8_switch(dec->cur.v, dec->refn[ref].v, 8 * x_pos, 8 * y_pos,
1244 :     uv_dx, uv_dy, stride2, 0);
1245 : chenm001 1.14 stop_comp_timer();
1246 :    
1247 : edgomez 1.19 for (i = 0; i < 6; i++) {
1248 : h 1.37.2.1 int direction = dec->alternate_vertical_scan ? 2 : 0;
1249 :    
1250 : edgomez 1.19 if (cbp & (1 << (5 - i))) // coded
1251 : chenm001 1.14 {
1252 : edgomez 1.19 memset(&block[i * 64], 0, 64 * sizeof(int16_t)); // clear
1253 : chenm001 1.14
1254 :     start_timer();
1255 : h 1.37.2.1 get_inter_block(bs, &block[i * 64], direction);
1256 : chenm001 1.14 stop_coding_timer();
1257 :    
1258 :     start_timer();
1259 : edgomez 1.19 if (dec->quant_type == 0) {
1260 :     dequant_inter(&data[i * 64], &block[i * 64], iQuant);
1261 :     } else {
1262 :     dequant4_inter(&data[i * 64], &block[i * 64], iQuant);
1263 : chenm001 1.14 }
1264 :     stop_iquant_timer();
1265 :    
1266 :     start_timer();
1267 : edgomez 1.19 idct(&data[i * 64]);
1268 : chenm001 1.14 stop_idct_timer();
1269 :     }
1270 :     }
1271 :    
1272 : edgomez 1.19 if (dec->interlacing && pMB->field_dct) {
1273 : chenm001 1.14 next_block = stride;
1274 :     stride *= 2;
1275 :     }
1276 :    
1277 :     start_timer();
1278 :     if (cbp & 32)
1279 : edgomez 1.19 transfer_16to8add(pY_Cur, &data[0 * 64], stride);
1280 : chenm001 1.14 if (cbp & 16)
1281 : edgomez 1.19 transfer_16to8add(pY_Cur + 8, &data[1 * 64], stride);
1282 : chenm001 1.14 if (cbp & 8)
1283 : edgomez 1.19 transfer_16to8add(pY_Cur + next_block, &data[2 * 64], stride);
1284 : chenm001 1.14 if (cbp & 4)
1285 : edgomez 1.19 transfer_16to8add(pY_Cur + 8 + next_block, &data[3 * 64], stride);
1286 : chenm001 1.14 if (cbp & 2)
1287 : edgomez 1.19 transfer_16to8add(pU_Cur, &data[4 * 64], stride2);
1288 : chenm001 1.14 if (cbp & 1)
1289 : edgomez 1.19 transfer_16to8add(pV_Cur, &data[5 * 64], stride2);
1290 : chenm001 1.14 stop_transfer_timer();
1291 :     }
1292 :    
1293 :     // add by MinChen <chenm001@163.com>
1294 :     // decode an B-frame direct & inter macroblock
1295 : edgomez 1.19 void
1296 :     decoder_bf_interpolate_mbinter(DECODER * dec,
1297 :     IMAGE forward,
1298 :     IMAGE backward,
1299 :     const MACROBLOCK * pMB,
1300 :     const uint32_t x_pos,
1301 :     const uint32_t y_pos,
1302 :     Bitstream * bs)
1303 : chenm001 1.14 {
1304 :    
1305 : edgomez 1.19 DECLARE_ALIGNED_MATRIX(block, 6, 64, int16_t, CACHE_LINE);
1306 : chenm001 1.14 DECLARE_ALIGNED_MATRIX(data, 6, 64, int16_t, CACHE_LINE);
1307 :    
1308 : edgomez 1.19 uint32_t stride = dec->edged_width;
1309 :     uint32_t stride2 = stride / 2;
1310 :     uint32_t next_block = stride * 8;
1311 :     uint32_t iQuant = pMB->quant;
1312 :     int uv_dx, uv_dy;
1313 :     int b_uv_dx, b_uv_dy;
1314 :     uint32_t i;
1315 :     uint8_t *pY_Cur, *pU_Cur, *pV_Cur;
1316 : chenm001 1.30 const uint32_t cbp = pMB->cbp;
1317 : chenm001 1.15
1318 :     pY_Cur = dec->cur.y + (y_pos << 4) * stride + (x_pos << 4);
1319 :     pU_Cur = dec->cur.u + (y_pos << 3) * stride2 + (x_pos << 3);
1320 :     pV_Cur = dec->cur.v + (y_pos << 3) * stride2 + (x_pos << 3);
1321 : chenm001 1.14
1322 : chenm001 1.30
1323 : edgomez 1.19 if ((pMB->mode == MODE_INTER || pMB->mode == MODE_INTER_Q)) {
1324 : chenm001 1.14 uv_dx = pMB->mvs[0].x;
1325 :     uv_dy = pMB->mvs[0].y;
1326 :    
1327 :     b_uv_dx = pMB->b_mvs[0].x;
1328 :     b_uv_dy = pMB->b_mvs[0].y;
1329 :    
1330 : Isibaar 1.37.2.13 if (dec->quarterpel)
1331 :     {
1332 :     uv_dx /= 2;
1333 :     uv_dy /= 2;
1334 :    
1335 :     b_uv_dx /= 2;
1336 :     b_uv_dy /= 2;
1337 :     }
1338 :    
1339 :     uv_dx = (uv_dx >> 1) + roundtab_79[uv_dx & 0x3];
1340 :     uv_dy = (uv_dy >> 1) + roundtab_79[uv_dy & 0x3];
1341 :    
1342 :     b_uv_dx = (b_uv_dx >> 1) + roundtab_79[b_uv_dx & 0x3];
1343 :     b_uv_dy = (b_uv_dy >> 1) + roundtab_79[b_uv_dy & 0x3];
1344 : edgomez 1.19 } else {
1345 : chenm001 1.14 int sum;
1346 : edgomez 1.19
1347 : Isibaar 1.37.2.13 if(dec->quarterpel)
1348 :     sum = (pMB->mvs[0].x / 2) + (pMB->mvs[1].x / 2) + (pMB->mvs[2].x / 2) + (pMB->mvs[3].x / 2);
1349 :     else
1350 :     sum = pMB->mvs[0].x + pMB->mvs[1].x + pMB->mvs[2].x + pMB->mvs[3].x;
1351 :    
1352 :     uv_dx = (sum >> 3) + roundtab_76[sum & 0xf];
1353 :    
1354 :     if(dec->quarterpel)
1355 :     sum = (pMB->mvs[0].y / 2) + (pMB->mvs[1].y / 2) + (pMB->mvs[2].y / 2) + (pMB->mvs[3].y / 2);
1356 :     else
1357 :     sum = pMB->mvs[0].y + pMB->mvs[1].y + pMB->mvs[2].y + pMB->mvs[3].y;
1358 :    
1359 :     uv_dy = (sum >> 3) + roundtab_76[sum & 0xf];
1360 :    
1361 :    
1362 :     if(dec->quarterpel)
1363 :     sum = (pMB->b_mvs[0].x / 2) + (pMB->b_mvs[1].x / 2) + (pMB->b_mvs[2].x / 2) + (pMB->b_mvs[3].x / 2);
1364 :     else
1365 :     sum = pMB->b_mvs[0].x + pMB->b_mvs[1].x + pMB->b_mvs[2].x + pMB->b_mvs[3].x;
1366 :    
1367 :     b_uv_dx = (sum >> 3) + roundtab_76[sum & 0xf];
1368 :    
1369 :     if(dec->quarterpel)
1370 :     sum = (pMB->b_mvs[0].y / 2) + (pMB->b_mvs[1].y / 2) + (pMB->b_mvs[2].y / 2) + (pMB->b_mvs[3].y / 2);
1371 :     else
1372 :     sum = pMB->b_mvs[0].y + pMB->b_mvs[1].y + pMB->b_mvs[2].y + pMB->b_mvs[3].y;
1373 :    
1374 :     b_uv_dy = (sum >> 3) + roundtab_76[sum & 0xf];
1375 : edgomez 1.19 }
1376 : chenm001 1.14
1377 :    
1378 :     start_timer();
1379 : Isibaar 1.37.2.13 if(dec->quarterpel) {
1380 :     if((pMB->mode == MODE_INTER || pMB->mode == MODE_INTER_Q))
1381 : suxen_drol 1.37.2.21 interpolate16x16_quarterpel(dec->cur.y, forward.y, dec->qtmp.y, dec->qtmp.y + 64,
1382 :     dec->qtmp.y + 128, 16*x_pos, 16*y_pos,
1383 : Isibaar 1.37.2.13 pMB->mvs[0].x, pMB->mvs[0].y, stride, 0);
1384 :     else {
1385 : suxen_drol 1.37.2.21 interpolate8x8_quarterpel(dec->cur.y, forward.y, dec->qtmp.y, dec->qtmp.y + 64,
1386 :     dec->qtmp.y + 128, 16*x_pos, 16*y_pos,
1387 : Isibaar 1.37.2.13 pMB->mvs[0].x, pMB->mvs[0].y, stride, 0);
1388 : suxen_drol 1.37.2.21 interpolate8x8_quarterpel(dec->cur.y, forward.y, dec->qtmp.y, dec->qtmp.y + 64,
1389 :     dec->qtmp.y + 128, 16*x_pos + 8, 16*y_pos,
1390 : Isibaar 1.37.2.13 pMB->mvs[1].x, pMB->mvs[1].y, stride, 0);
1391 : suxen_drol 1.37.2.21 interpolate8x8_quarterpel(dec->cur.y, forward.y, dec->qtmp.y, dec->qtmp.y + 64,
1392 :     dec->qtmp.y + 128, 16*x_pos, 16*y_pos + 8,
1393 : Isibaar 1.37.2.13 pMB->mvs[2].x, pMB->mvs[2].y, stride, 0);
1394 : suxen_drol 1.37.2.21 interpolate8x8_quarterpel(dec->cur.y, forward.y, dec->qtmp.y, dec->qtmp.y + 64,
1395 :     dec->qtmp.y + 128, 16*x_pos + 8, 16*y_pos + 8,
1396 : Isibaar 1.37.2.13 pMB->mvs[3].x, pMB->mvs[3].y, stride, 0);
1397 :     }
1398 :     }
1399 :     else {
1400 :     interpolate8x8_switch(dec->cur.y, forward.y, 16 * x_pos, 16 * y_pos,
1401 :     pMB->mvs[0].x, pMB->mvs[0].y, stride, 0);
1402 :     interpolate8x8_switch(dec->cur.y, forward.y, 16 * x_pos + 8, 16 * y_pos,
1403 :     pMB->mvs[1].x, pMB->mvs[1].y, stride, 0);
1404 :     interpolate8x8_switch(dec->cur.y, forward.y, 16 * x_pos, 16 * y_pos + 8,
1405 :     pMB->mvs[2].x, pMB->mvs[2].y, stride, 0);
1406 :     interpolate8x8_switch(dec->cur.y, forward.y, 16 * x_pos + 8,
1407 :     16 * y_pos + 8, pMB->mvs[3].x, pMB->mvs[3].y, stride,
1408 :     0);
1409 :     }
1410 :    
1411 : edgomez 1.19 interpolate8x8_switch(dec->cur.u, forward.u, 8 * x_pos, 8 * y_pos, uv_dx,
1412 :     uv_dy, stride2, 0);
1413 :     interpolate8x8_switch(dec->cur.v, forward.v, 8 * x_pos, 8 * y_pos, uv_dx,
1414 :     uv_dy, stride2, 0);
1415 :    
1416 :    
1417 : Isibaar 1.37.2.13 if(dec->quarterpel) {
1418 :     if((pMB->mode == MODE_INTER || pMB->mode == MODE_INTER_Q))
1419 : suxen_drol 1.37.2.21 interpolate16x16_quarterpel(dec->tmp.y, backward.y, dec->qtmp.y, dec->qtmp.y + 64,
1420 :     dec->qtmp.y + 128, 16*x_pos, 16*y_pos,
1421 : Isibaar 1.37.2.13 pMB->b_mvs[0].x, pMB->b_mvs[0].y, stride, 0);
1422 :     else {
1423 : suxen_drol 1.37.2.21 interpolate8x8_quarterpel(dec->tmp.y, backward.y, dec->qtmp.y, dec->qtmp.y + 64,
1424 :     dec->qtmp.y + 128, 16*x_pos, 16*y_pos,
1425 : Isibaar 1.37.2.13 pMB->b_mvs[0].x, pMB->b_mvs[0].y, stride, 0);
1426 : suxen_drol 1.37.2.21 interpolate8x8_quarterpel(dec->tmp.y, backward.y, dec->qtmp.y, dec->qtmp.y + 64,
1427 :     dec->qtmp.y + 128, 16*x_pos + 8, 16*y_pos,
1428 : Isibaar 1.37.2.13 pMB->b_mvs[1].x, pMB->b_mvs[1].y, stride, 0);
1429 : suxen_drol 1.37.2.21 interpolate8x8_quarterpel(dec->tmp.y, backward.y, dec->qtmp.y, dec->qtmp.y + 64,
1430 :     dec->qtmp.y + 128, 16*x_pos, 16*y_pos + 8,
1431 : Isibaar 1.37.2.13 pMB->b_mvs[2].x, pMB->b_mvs[2].y, stride, 0);
1432 : suxen_drol 1.37.2.21 interpolate8x8_quarterpel(dec->tmp.y, backward.y, dec->qtmp.y, dec->qtmp.y + 64,
1433 :     dec->qtmp.y + 128, 16*x_pos + 8, 16*y_pos + 8,
1434 : Isibaar 1.37.2.13 pMB->b_mvs[3].x, pMB->b_mvs[3].y, stride, 0);
1435 :     }
1436 :     }
1437 :     else {
1438 : suxen_drol 1.37.2.21 interpolate8x8_switch(dec->tmp.y, backward.y, 16 * x_pos, 16 * y_pos,
1439 : Isibaar 1.37.2.13 pMB->b_mvs[0].x, pMB->b_mvs[0].y, stride, 0);
1440 : suxen_drol 1.37.2.21 interpolate8x8_switch(dec->tmp.y, backward.y, 16 * x_pos + 8,
1441 : Isibaar 1.37.2.13 16 * y_pos, pMB->b_mvs[1].x, pMB->b_mvs[1].y, stride,
1442 :     0);
1443 : suxen_drol 1.37.2.21 interpolate8x8_switch(dec->tmp.y, backward.y, 16 * x_pos,
1444 : Isibaar 1.37.2.13 16 * y_pos + 8, pMB->b_mvs[2].x, pMB->b_mvs[2].y,
1445 :     stride, 0);
1446 : suxen_drol 1.37.2.21 interpolate8x8_switch(dec->tmp.y, backward.y, 16 * x_pos + 8,
1447 : Isibaar 1.37.2.13 16 * y_pos + 8, pMB->b_mvs[3].x, pMB->b_mvs[3].y,
1448 :     stride, 0);
1449 :     }
1450 :    
1451 : suxen_drol 1.37.2.21 interpolate8x8_switch(dec->tmp.u, backward.u, 8 * x_pos, 8 * y_pos,
1452 : edgomez 1.19 b_uv_dx, b_uv_dy, stride2, 0);
1453 : suxen_drol 1.37.2.21 interpolate8x8_switch(dec->tmp.v, backward.v, 8 * x_pos, 8 * y_pos,
1454 : edgomez 1.19 b_uv_dx, b_uv_dy, stride2, 0);
1455 :    
1456 : Isibaar 1.37.2.2 interpolate8x8_avg2(dec->cur.y + (16 * y_pos * stride) + 16 * x_pos,
1457 :     dec->cur.y + (16 * y_pos * stride) + 16 * x_pos,
1458 : suxen_drol 1.37.2.21 dec->tmp.y + (16 * y_pos * stride) + 16 * x_pos,
1459 : Isibaar 1.37.2.14 stride, 1, 8);
1460 : Isibaar 1.37.2.2
1461 :     interpolate8x8_avg2(dec->cur.y + (16 * y_pos * stride) + 16 * x_pos + 8,
1462 :     dec->cur.y + (16 * y_pos * stride) + 16 * x_pos + 8,
1463 : suxen_drol 1.37.2.21 dec->tmp.y + (16 * y_pos * stride) + 16 * x_pos + 8,
1464 : Isibaar 1.37.2.14 stride, 1, 8);
1465 : Isibaar 1.37.2.2
1466 : Isibaar 1.37.2.4 interpolate8x8_avg2(dec->cur.y + ((16 * y_pos + 8) * stride) + 16 * x_pos,
1467 :     dec->cur.y + ((16 * y_pos + 8) * stride) + 16 * x_pos,
1468 : suxen_drol 1.37.2.21 dec->tmp.y + ((16 * y_pos + 8) * stride) + 16 * x_pos,
1469 : Isibaar 1.37.2.14 stride, 1, 8);
1470 : Isibaar 1.37.2.2
1471 : Isibaar 1.37.2.4 interpolate8x8_avg2(dec->cur.y + ((16 * y_pos + 8) * stride) + 16 * x_pos + 8,
1472 :     dec->cur.y + ((16 * y_pos + 8) * stride) + 16 * x_pos + 8,
1473 : suxen_drol 1.37.2.21 dec->tmp.y + ((16 * y_pos + 8) * stride) + 16 * x_pos + 8,
1474 : Isibaar 1.37.2.14 stride, 1, 8);
1475 : Isibaar 1.37.2.2
1476 : Isibaar 1.37.2.4 interpolate8x8_avg2(dec->cur.u + (8 * y_pos * stride2) + 8 * x_pos,
1477 :     dec->cur.u + (8 * y_pos * stride2) + 8 * x_pos,
1478 : suxen_drol 1.37.2.21 dec->tmp.u + (8 * y_pos * stride2) + 8 * x_pos,
1479 : Isibaar 1.37.2.14 stride2, 1, 8);
1480 : Isibaar 1.37.2.2
1481 : Isibaar 1.37.2.4 interpolate8x8_avg2(dec->cur.v + (8 * y_pos * stride2) + 8 * x_pos,
1482 :     dec->cur.v + (8 * y_pos * stride2) + 8 * x_pos,
1483 : suxen_drol 1.37.2.21 dec->tmp.v + (8 * y_pos * stride2) + 8 * x_pos,
1484 : Isibaar 1.37.2.14 stride2, 1, 8);
1485 : Isibaar 1.37.2.2
1486 : chenm001 1.14 stop_comp_timer();
1487 : edgomez 1.19
1488 :     for (i = 0; i < 6; i++) {
1489 : h 1.37.2.1 int direction = dec->alternate_vertical_scan ? 2 : 0;
1490 :    
1491 : edgomez 1.19 if (cbp & (1 << (5 - i))) // coded
1492 : chenm001 1.15 {
1493 : edgomez 1.19 memset(&block[i * 64], 0, 64 * sizeof(int16_t)); // clear
1494 : chenm001 1.14
1495 : chenm001 1.15 start_timer();
1496 : h 1.37.2.1 get_inter_block(bs, &block[i * 64], direction);
1497 : chenm001 1.15 stop_coding_timer();
1498 :    
1499 :     start_timer();
1500 : edgomez 1.19 if (dec->quant_type == 0) {
1501 :     dequant_inter(&data[i * 64], &block[i * 64], iQuant);
1502 :     } else {
1503 :     dequant4_inter(&data[i * 64], &block[i * 64], iQuant);
1504 : chenm001 1.15 }
1505 :     stop_iquant_timer();
1506 :    
1507 :     start_timer();
1508 : edgomez 1.19 idct(&data[i * 64]);
1509 : chenm001 1.15 stop_idct_timer();
1510 :     }
1511 :     }
1512 :    
1513 : edgomez 1.19 if (dec->interlacing && pMB->field_dct) {
1514 : chenm001 1.15 next_block = stride;
1515 :     stride *= 2;
1516 :     }
1517 :    
1518 :     start_timer();
1519 :     if (cbp & 32)
1520 : edgomez 1.19 transfer_16to8add(pY_Cur, &data[0 * 64], stride);
1521 : chenm001 1.15 if (cbp & 16)
1522 : edgomez 1.19 transfer_16to8add(pY_Cur + 8, &data[1 * 64], stride);
1523 : chenm001 1.15 if (cbp & 8)
1524 : edgomez 1.19 transfer_16to8add(pY_Cur + next_block, &data[2 * 64], stride);
1525 : chenm001 1.15 if (cbp & 4)
1526 : edgomez 1.19 transfer_16to8add(pY_Cur + 8 + next_block, &data[3 * 64], stride);
1527 : chenm001 1.15 if (cbp & 2)
1528 : edgomez 1.19 transfer_16to8add(pU_Cur, &data[4 * 64], stride2);
1529 : chenm001 1.15 if (cbp & 1)
1530 : edgomez 1.19 transfer_16to8add(pV_Cur, &data[5 * 64], stride2);
1531 : chenm001 1.15 stop_transfer_timer();
1532 : chenm001 1.14 }
1533 :    
1534 :    
1535 :     // add by MinChen <chenm001@163.com>
1536 :     // for decode B-frame dbquant
1537 : edgomez 1.19 int32_t __inline
1538 :     get_dbquant(Bitstream * bs)
1539 : chenm001 1.14 {
1540 : edgomez 1.19 if (!BitstreamGetBit(bs)) // '0'
1541 :     return (0);
1542 : chenm001 1.14 else if (!BitstreamGetBit(bs)) // '10'
1543 : edgomez 1.19 return (-2);
1544 : chenm001 1.14 else
1545 : edgomez 1.19 return (2); // '11'
1546 : chenm001 1.14 }
1547 :    
1548 :     // add by MinChen <chenm001@163.com>
1549 :     // for decode B-frame mb_type
1550 :     // bit ret_value
1551 :     // 1 0
1552 :     // 01 1
1553 :     // 001 2
1554 :     // 0001 3
1555 : edgomez 1.19 int32_t __inline
1556 :     get_mbtype(Bitstream * bs)
1557 : chenm001 1.14 {
1558 : chenm001 1.15 int32_t mb_type;
1559 : edgomez 1.19
1560 :     for (mb_type = 0; mb_type <= 3; mb_type++) {
1561 :     if (BitstreamGetBit(bs))
1562 : chenm001 1.15 break;
1563 :     }
1564 :    
1565 : edgomez 1.19 if (mb_type <= 3)
1566 :     return (mb_type);
1567 : chenm001 1.14 else
1568 : edgomez 1.19 return (-1);
1569 : chenm001 1.14 }
1570 :    
1571 : edgomez 1.19 void
1572 :     decoder_bframe(DECODER * dec,
1573 :     Bitstream * bs,
1574 :     int quant,
1575 :     int fcode_forward,
1576 :     int fcode_backward)
1577 : chenm001 1.14 {
1578 : edgomez 1.19 uint32_t x, y;
1579 : chenm001 1.30 VECTOR mv;
1580 :     const VECTOR zeromv = {0,0};
1581 : chenm001 1.27 #ifdef BFRAMES_DEC_DEBUG
1582 :     FILE *fp;
1583 :     static char first=0;
1584 :     #define BFRAME_DEBUG if (!first && fp){ \
1585 :     fprintf(fp,"Y=%3d X=%3d MB=%2d CBP=%02X\n",y,x,mb->mb_type,mb->cbp); \
1586 :     }
1587 :     #endif
1588 : chenm001 1.14
1589 :     start_timer();
1590 : edgomez 1.19 image_setedges(&dec->refn[0], dec->edged_width, dec->edged_height,
1591 : h 1.37.2.1 dec->width, dec->height);
1592 : chenm001 1.30 image_setedges(&dec->refn[1], dec->edged_width, dec->edged_height,
1593 : h 1.37.2.1 dec->width, dec->height);
1594 : chenm001 1.14 stop_edges_timer();
1595 :    
1596 : chenm001 1.27 #ifdef BFRAMES_DEC_DEBUG
1597 :     if (!first){
1598 :     fp=fopen("C:\\XVIDDBG.TXT","w");
1599 :     }
1600 :     #endif
1601 : chenm001 1.15
1602 : edgomez 1.19 for (y = 0; y < dec->mb_height; y++) {
1603 : chenm001 1.15 // Initialize Pred Motion Vector
1604 : chenm001 1.30 dec->p_fmv = dec->p_bmv = zeromv;
1605 : edgomez 1.19 for (x = 0; x < dec->mb_width; x++) {
1606 :     MACROBLOCK *mb = &dec->mbs[y * dec->mb_width + x];
1607 :     MACROBLOCK *last_mb = &dec->last_mbs[y * dec->mb_width + x];
1608 : chenm001 1.14
1609 : chenm001 1.30 mv =
1610 :     mb->b_mvs[0] = mb->b_mvs[1] = mb->b_mvs[2] = mb->b_mvs[3] =
1611 :     mb->mvs[0] = mb->mvs[1] = mb->mvs[2] = mb->mvs[3] = zeromv;
1612 : chenm001 1.15
1613 : suxen_drol 1.37.2.24 // skip if the co-located P_VOP macroblock is not coded
1614 : chl 1.37.2.28 // if not codec in co-located S_VOP macroblock is _not_ automatically skipped
1615 : suxen_drol 1.37.2.24
1616 : edgomez 1.19 if (last_mb->mode == MODE_NOT_CODED) {
1617 : chenm001 1.15 //DEBUG2("Skip MB in B-frame at (X,Y)=!",x,y);
1618 : chenm001 1.27 mb->cbp = 0;
1619 :     #ifdef BFRAMES_DEC_DEBUG
1620 : chenm001 1.29 mb->mb_type = MODE_NOT_CODED;
1621 : chenm001 1.27 BFRAME_DEBUG
1622 :     #endif
1623 : chenm001 1.15 mb->mb_type = MODE_FORWARD;
1624 : chenm001 1.30 mb->quant = last_mb->quant;
1625 :     //mb->mvs[1].x = mb->mvs[2].x = mb->mvs[3].x = mb->mvs[0].x;
1626 :     //mb->mvs[1].y = mb->mvs[2].y = mb->mvs[3].y = mb->mvs[0].y;
1627 : chenm001 1.14
1628 : chenm001 1.30 decoder_bf_mbinter(dec, mb, x, y, mb->cbp, bs, mb->quant, 1);
1629 : chenm001 1.14 continue;
1630 :     }
1631 : chenm001 1.15
1632 : edgomez 1.19 if (!BitstreamGetBit(bs)) { // modb=='0'
1633 :     const uint8_t modb2 = BitstreamGetBit(bs);
1634 : chenm001 1.14
1635 :     mb->mb_type = get_mbtype(bs);
1636 : edgomez 1.19
1637 :     if (!modb2) { // modb=='00'
1638 :     mb->cbp = BitstreamGetBits(bs, 6);
1639 : chenm001 1.14 } else {
1640 :     mb->cbp = 0;
1641 :     }
1642 : edgomez 1.19 if (mb->mb_type && mb->cbp) {
1643 : chenm001 1.14 quant += get_dbquant(bs);
1644 :    
1645 : edgomez 1.19 if (quant > 31) {
1646 : chenm001 1.14 quant = 31;
1647 : chenm001 1.26 } else if (quant < 1) {
1648 : chenm001 1.14 quant = 1;
1649 :     }
1650 :     }
1651 :     } else {
1652 :     mb->mb_type = MODE_DIRECT_NONE_MV;
1653 : edgomez 1.19 mb->cbp = 0;
1654 : chenm001 1.14 }
1655 :    
1656 : chenm001 1.30 mb->quant = quant;
1657 :     mb->mode = MODE_INTER4V;
1658 : chenm001 1.14 //DEBUG1("Switch bm_type=",mb->mb_type);
1659 :    
1660 : chenm001 1.27 #ifdef BFRAMES_DEC_DEBUG
1661 :     BFRAME_DEBUG
1662 :     #endif
1663 : chl 1.33
1664 : edgomez 1.19 switch (mb->mb_type) {
1665 : chenm001 1.14 case MODE_DIRECT:
1666 : chenm001 1.30 get_b_motion_vector(dec, bs, x, y, &mv, 1, zeromv);
1667 : chenm001 1.14
1668 :     case MODE_DIRECT_NONE_MV:
1669 : chl 1.33 {
1670 : chenm001 1.30 const int64_t TRB = dec->time_pp - dec->time_bp, TRD = dec->time_pp;
1671 : chenm001 1.14 int i;
1672 : edgomez 1.19
1673 :     for (i = 0; i < 4; i++) {
1674 : chenm001 1.30 mb->mvs[i].x = (int32_t) ((TRB * last_mb->mvs[i].x)
1675 :     / TRD + mv.x);
1676 :     mb->b_mvs[i].x = (int32_t) ((mv.x == 0)
1677 :     ? ((TRB - TRD) * last_mb->mvs[i].x)
1678 :     / TRD
1679 :     : mb->mvs[i].x - last_mb->mvs[i].x);
1680 :     mb->mvs[i].y = (int32_t) ((TRB * last_mb->mvs[i].y)
1681 :     / TRD + mv.y);
1682 :     mb->b_mvs[i].y = (int32_t) ((mv.y == 0)
1683 :     ? ((TRB - TRD) * last_mb->mvs[i].y)
1684 :     / TRD
1685 :     : mb->mvs[i].y - last_mb->mvs[i].y);
1686 : chenm001 1.14 }
1687 :     //DEBUG("B-frame Direct!\n");
1688 :     }
1689 : edgomez 1.19 decoder_bf_interpolate_mbinter(dec, dec->refn[1], dec->refn[0],
1690 : chenm001 1.30 mb, x, y, bs);
1691 : chenm001 1.14 break;
1692 :    
1693 :     case MODE_INTERPOLATE:
1694 : edgomez 1.19 get_b_motion_vector(dec, bs, x, y, &mb->mvs[0], fcode_forward,
1695 :     dec->p_fmv);
1696 : chenm001 1.30 dec->p_fmv = mb->mvs[1] = mb->mvs[2] = mb->mvs[3] = mb->mvs[0];
1697 : edgomez 1.19
1698 :     get_b_motion_vector(dec, bs, x, y, &mb->b_mvs[0],
1699 :     fcode_backward, dec->p_bmv);
1700 : chenm001 1.30 dec->p_bmv = mb->b_mvs[1] = mb->b_mvs[2] =
1701 :     mb->b_mvs[3] = mb->b_mvs[0];
1702 : chenm001 1.14
1703 : edgomez 1.19 decoder_bf_interpolate_mbinter(dec, dec->refn[1], dec->refn[0],
1704 : chenm001 1.30 mb, x, y, bs);
1705 : chenm001 1.14 //DEBUG("B-frame Bidir!\n");
1706 :     break;
1707 :    
1708 :     case MODE_BACKWARD:
1709 : edgomez 1.19 get_b_motion_vector(dec, bs, x, y, &mb->mvs[0], fcode_backward,
1710 :     dec->p_bmv);
1711 : chenm001 1.30 dec->p_bmv = mb->mvs[1] = mb->mvs[2] = mb->mvs[3] = mb->mvs[0];
1712 : chenm001 1.15
1713 : chenm001 1.30 mb->mode = MODE_INTER;
1714 : chenm001 1.14 decoder_bf_mbinter(dec, mb, x, y, mb->cbp, bs, quant, 0);
1715 :     //DEBUG("B-frame Backward!\n");
1716 :     break;
1717 :    
1718 :     case MODE_FORWARD:
1719 : edgomez 1.19 get_b_motion_vector(dec, bs, x, y, &mb->mvs[0], fcode_forward,
1720 :     dec->p_fmv);
1721 : chenm001 1.30 dec->p_fmv = mb->mvs[1] = mb->mvs[2] = mb->mvs[3] = mb->mvs[0];
1722 : chenm001 1.15
1723 : chenm001 1.30 mb->mode = MODE_INTER;
1724 : chenm001 1.14 decoder_bf_mbinter(dec, mb, x, y, mb->cbp, bs, quant, 1);
1725 : chenm001 1.15 //DEBUG("B-frame Forward!\n");
1726 : chenm001 1.14 break;
1727 :    
1728 :     default:
1729 : suxen_drol 1.37.2.22 DPRINTF(DPRINTF_ERROR,"Not support B-frame mb_type = %i", mb->mb_type);
1730 : chl 1.33 }
1731 : chenm001 1.14
1732 : edgomez 1.19 } // end of FOR
1733 : chenm001 1.14 }
1734 : chenm001 1.27 #ifdef BFRAMES_DEC_DEBUG
1735 :     if (!first){
1736 :     first=1;
1737 :     if (fp)
1738 :     fclose(fp);
1739 :     }
1740 :     #endif
1741 : chenm001 1.14 }
1742 :    
1743 :     // swap two MACROBLOCK array
1744 : edgomez 1.19 void
1745 :     mb_swap(MACROBLOCK ** mb1,
1746 :     MACROBLOCK ** mb2)
1747 : chenm001 1.14 {
1748 : edgomez 1.19 MACROBLOCK *temp = *mb1;
1749 :    
1750 :     *mb1 = *mb2;
1751 :     *mb2 = temp;
1752 : chenm001 1.14 }
1753 :    
1754 : suxen_drol 1.37.2.19
1755 : suxen_drol 1.37.2.21 /* perform post processing if necessary, and output the image */
1756 :     void decoder_output(DECODER * dec, IMAGE * img, MACROBLOCK * mbs,
1757 :     const XVID_DEC_FRAME * frame, int pp_disable)
1758 :     {
1759 :    
1760 :     if ((frame->general & (XVID_DEC_DEBLOCKY|XVID_DEC_DEBLOCKUV)) && !pp_disable) /* post process */
1761 :     {
1762 :     /* note: image is stored to tmp */
1763 :     image_copy(&dec->tmp, img, dec->edged_width, dec->height);
1764 :     image_deblock_rrv(&dec->tmp, dec->edged_width,
1765 :     mbs, dec->mb_width, dec->mb_height, dec->mb_width,
1766 :     8, frame->general);
1767 :     img = &dec->tmp;
1768 :     }
1769 :    
1770 :     image_output(img, dec->width, dec->height,
1771 :     dec->edged_width, frame->image, frame->stride,
1772 :     frame->colorspace, dec->interlacing);
1773 :     }
1774 :    
1775 :    
1776 : edgomez 1.19 int
1777 :     decoder_decode(DECODER * dec,
1778 : suxen_drol 1.37.2.8 XVID_DEC_FRAME * frame, XVID_DEC_STATS * stats)
1779 : Isibaar 1.1 {
1780 : edgomez 1.7
1781 : Isibaar 1.1 Bitstream bs;
1782 :     uint32_t rounding;
1783 : suxen_drol 1.37.2.8 uint32_t reduced_resolution;
1784 : Isibaar 1.1 uint32_t quant;
1785 : chenm001 1.14 uint32_t fcode_forward;
1786 :     uint32_t fcode_backward;
1787 : Isibaar 1.1 uint32_t intra_dc_threshold;
1788 : chl 1.37.2.27 WARPPOINTS gmc_warp;
1789 : chenm001 1.11 uint32_t vop_type;
1790 : suxen_drol 1.37.2.12 int success = 0;
1791 : suxen_drol 1.37.2.19 int output = 0;
1792 :     int seen_something = 0;
1793 : Isibaar 1.1
1794 :     start_global_timer();
1795 : albeu 1.32
1796 : suxen_drol 1.37.2.19 dec->low_delay_default = (frame->general & XVID_DEC_LOWDELAY);
1797 : albeu 1.32 dec->out_frm = (frame->colorspace == XVID_CSP_EXTERN) ? frame->image : NULL;
1798 : edgomez 1.19
1799 : suxen_drol 1.37.2.23 if ((frame->general & XVID_DEC_DISCONTINUITY))
1800 :     dec->frames = 0;
1801 :    
1802 : suxen_drol 1.37.2.21 if (frame->length < 0) /* decoder flush */
1803 :     {
1804 :     /* if not decoding "low_delay/packed", and this isn't low_delay and
1805 :     we have a reference frame, then outout the reference frame */
1806 :     if (!(dec->low_delay_default && dec->packed_mode) && !dec->low_delay && dec->frames>0)
1807 :     {
1808 :     decoder_output(dec, &dec->refn[0], dec->mbs, frame, dec->last_reduced_resolution);
1809 :     output = 1;
1810 :     }
1811 :    
1812 :     frame->length = 0;
1813 :     if (stats)
1814 :     {
1815 :     stats->notify = output ? XVID_DEC_VOP : XVID_DEC_NOTHING;
1816 :     stats->data.vop.time_base = (int)dec->time_base;
1817 :     stats->data.vop.time_increment = 0; //XXX: todo
1818 :     }
1819 :    
1820 :     emms();
1821 :    
1822 :     stop_global_timer();
1823 :     return XVID_ERR_OK;
1824 :     }
1825 :    
1826 : Isibaar 1.1 BitstreamInit(&bs, frame->bitstream, frame->length);
1827 :    
1828 : suxen_drol 1.37.2.8 // XXX: 0x7f is only valid whilst decoding vfw xvid/divx5 avi's
1829 : suxen_drol 1.37.2.21 if(dec->low_delay_default && frame->length == 1 && BitstreamShowBits(&bs, 8) == 0x7f)
1830 : suxen_drol 1.37.2.8 {
1831 :     if (stats)
1832 :     stats->notify = XVID_DEC_VOP;
1833 :     frame->length = 1;
1834 : Isibaar 1.37.2.10 image_output(&dec->refn[0], dec->width, dec->height, dec->edged_width,
1835 : suxen_drol 1.37.2.8 frame->image, frame->stride, frame->colorspace, dec->interlacing);
1836 : Isibaar 1.37.2.10 emms();
1837 : Isibaar 1.37.2.5 return XVID_ERR_OK;
1838 : suxen_drol 1.37.2.8 }
1839 : Isibaar 1.37.2.4
1840 : suxen_drol 1.37.2.21 repeat:
1841 : suxen_drol 1.37.2.8
1842 : suxen_drol 1.37.2.21 vop_type = BitstreamReadHeaders(&bs, dec, &rounding, &reduced_resolution,
1843 : chl 1.37.2.27 &quant, &fcode_forward, &fcode_backward, &intra_dc_threshold, &gmc_warp);
1844 : suxen_drol 1.37.2.8
1845 : chl 1.37.2.27 DPRINTF(DPRINTF_HEADER, "vop_type=%i, packed=%i, time=%lli, time_pp=%i, time_bp=%i",
1846 : suxen_drol 1.37.2.19 vop_type, dec->packed_mode, dec->time, dec->time_pp, dec->time_bp);
1847 : suxen_drol 1.37.2.8
1848 : suxen_drol 1.37.2.19 if (vop_type == - 1)
1849 :     {
1850 :     if (success) goto done;
1851 : suxen_drol 1.37.2.25 emms();
1852 : suxen_drol 1.37.2.19 return XVID_ERR_FAIL;
1853 :     }
1854 : suxen_drol 1.37.2.12
1855 : suxen_drol 1.37.2.8 if (vop_type == -2 || vop_type == -3)
1856 :     {
1857 :     if (vop_type == -3)
1858 :     decoder_resize(dec);
1859 :    
1860 :     if (stats)
1861 :     {
1862 :     stats->notify = XVID_DEC_VOL;
1863 :     stats->data.vol.general = 0;
1864 :     if (dec->interlacing)
1865 :     stats->data.vol.general |= XVID_INTERLACING;
1866 :     stats->data.vol.width = dec->width;
1867 :     stats->data.vol.height = dec->height;
1868 :     stats->data.vol.aspect_ratio = dec->aspect_ratio;
1869 :     stats->data.vol.par_width = dec->par_width;
1870 :     stats->data.vol.par_height = dec->par_height;
1871 :     frame->length = BitstreamPos(&bs) / 8;
1872 : suxen_drol 1.37.2.25 emms();
1873 : suxen_drol 1.37.2.8 return XVID_ERR_OK;
1874 :     }
1875 : suxen_drol 1.37.2.21 goto repeat;
1876 : suxen_drol 1.37.2.8 }
1877 : edgomez 1.19
1878 :     dec->p_bmv.x = dec->p_bmv.y = dec->p_fmv.y = dec->p_fmv.y = 0; // init pred vector to 0
1879 :    
1880 : Isibaar 1.1
1881 : suxen_drol 1.37.2.19 /* packed_mode: special-N_VOP treament */
1882 :     if (dec->packed_mode && vop_type == N_VOP)
1883 :     {
1884 : suxen_drol 1.37.2.21 if (dec->low_delay_default && dec->frames > 0)
1885 : suxen_drol 1.37.2.19 {
1886 : suxen_drol 1.37.2.21 decoder_output(dec, &dec->refn[0], dec->last_mbs, frame, dec->last_reduced_resolution);
1887 : suxen_drol 1.37.2.19 output = 1;
1888 : chenm001 1.14 }
1889 : suxen_drol 1.37.2.19 /* ignore otherwise */
1890 : Isibaar 1.1 }
1891 : suxen_drol 1.37.2.19 else if (vop_type != B_VOP)
1892 : suxen_drol 1.37.2.15 {
1893 : suxen_drol 1.37.2.19 switch(vop_type)
1894 :     {
1895 :     case I_VOP :
1896 :     decoder_iframe(dec, &bs, reduced_resolution, quant, intra_dc_threshold);
1897 :     break;
1898 :     case P_VOP :
1899 :     decoder_pframe(dec, &bs, rounding, reduced_resolution, quant,
1900 :     fcode_forward, intra_dc_threshold, NULL);
1901 :     break;
1902 :     case S_VOP :
1903 :     decoder_pframe(dec, &bs, rounding, reduced_resolution, quant,
1904 : chl 1.37.2.27 fcode_forward, intra_dc_threshold, &gmc_warp);
1905 : suxen_drol 1.37.2.19 break;
1906 :     case N_VOP :
1907 :     image_copy(&dec->cur, &dec->refn[0], dec->edged_width, dec->height);
1908 :     break;
1909 :     }
1910 : chenm001 1.27
1911 : suxen_drol 1.37.2.19 if (reduced_resolution)
1912 :     {
1913 : suxen_drol 1.37.2.20 image_deblock_rrv(&dec->cur, dec->edged_width, dec->mbs,
1914 : suxen_drol 1.37.2.21 (dec->width + 31) / 32, (dec->height + 31) / 32, dec->mb_width,
1915 :     16, XVID_DEC_DEBLOCKY|XVID_DEC_DEBLOCKUV);
1916 : suxen_drol 1.37.2.19 }
1917 : chenm001 1.27
1918 : suxen_drol 1.37.2.19 /* note: for packed_mode, output is performed when the special-N_VOP is decoded */
1919 :     if (!(dec->low_delay_default && dec->packed_mode))
1920 :     {
1921 :     if (dec->low_delay)
1922 :     {
1923 : suxen_drol 1.37.2.21 decoder_output(dec, &dec->cur, dec->mbs, frame, reduced_resolution);
1924 : suxen_drol 1.37.2.19 output = 1;
1925 :     }
1926 :     else if (dec->frames > 0) /* is the reference frame valid? */
1927 :     {
1928 : suxen_drol 1.37.2.21 /* output the reference frame */
1929 :     decoder_output(dec, &dec->refn[0], dec->last_mbs, frame, dec->last_reduced_resolution);
1930 : suxen_drol 1.37.2.19 output = 1;
1931 : chenm001 1.17 }
1932 : chenm001 1.14 }
1933 : Isibaar 1.18
1934 : chenm001 1.14 image_swap(&dec->refn[0], &dec->refn[1]);
1935 :     image_swap(&dec->cur, &dec->refn[0]);
1936 : suxen_drol 1.37.2.19 mb_swap(&dec->mbs, &dec->last_mbs);
1937 : suxen_drol 1.37.2.21 dec->last_reduced_resolution = reduced_resolution;
1938 : suxen_drol 1.37.2.19
1939 :     dec->frames++;
1940 :     seen_something = 1;
1941 :    
1942 :     }else{ /* B_VOP */
1943 :    
1944 :     if (dec->low_delay)
1945 :     {
1946 :     DPRINTF(DPRINTF_ERROR, "warning: bvop found in low_delay==1 stream");
1947 :     dec->low_delay = 1;
1948 :     }
1949 : chenm001 1.31
1950 : suxen_drol 1.37.2.21 if (dec->frames < 2)
1951 :     {
1952 :     /* attemping to decode a bvop without atleast 2 reference frames */
1953 :     image_printf(&dec->cur, dec->edged_width, dec->height, 16, 16,
1954 :     "broken b-frame, mising ref frames");
1955 :     }else if (dec->time_pp <= dec->time_bp) {
1956 : suxen_drol 1.37.2.19 /* this occurs when dx50_bvop_compatibility==0 sequences are
1957 :     decoded in vfw. */
1958 :     image_printf(&dec->cur, dec->edged_width, dec->height, 16, 16,
1959 :     "broken b-frame, tpp=%i tbp=%i", dec->time_pp, dec->time_bp);
1960 :     }else{
1961 :     decoder_bframe(dec, &bs, quant, fcode_forward, fcode_backward);
1962 :     }
1963 : suxen_drol 1.37.2.21
1964 :     decoder_output(dec, &dec->cur, dec->mbs, frame, reduced_resolution);
1965 : suxen_drol 1.37.2.19 output = 1;
1966 :     dec->frames++;
1967 : suxen_drol 1.37.2.8 }
1968 : Isibaar 1.37.2.11
1969 : suxen_drol 1.37.2.19 BitstreamByteAlign(&bs);
1970 : suxen_drol 1.37.2.12
1971 : suxen_drol 1.37.2.19 /* low_delay_default mode: repeat in packed_mode */
1972 :     if (dec->low_delay_default && dec->packed_mode && output == 0 && success == 0)
1973 : suxen_drol 1.37.2.15 {
1974 :     success = 1;
1975 : suxen_drol 1.37.2.21 goto repeat;
1976 : suxen_drol 1.37.2.15 }
1977 : suxen_drol 1.37.2.12
1978 :     done :
1979 : Isibaar 1.37.2.11
1980 : suxen_drol 1.37.2.19 /* low_delay_default mode: if we've gotten here without outputting anything,
1981 :     then output the recently decoded frame, or print an error message */
1982 :     if (dec->low_delay_default && output == 0)
1983 :     {
1984 :     if (dec->packed_mode && seen_something)
1985 :     {
1986 :     /* output the recently decoded frame */
1987 : suxen_drol 1.37.2.21 decoder_output(dec, &dec->refn[0], dec->last_mbs, frame, dec->last_reduced_resolution);
1988 : suxen_drol 1.37.2.23 output = 1;
1989 : suxen_drol 1.37.2.19 }
1990 :     else
1991 :     {
1992 : suxen_drol 1.37.2.21 image_clear(&dec->cur, dec->width, dec->height, dec->edged_width, 0, 128, 128);
1993 : suxen_drol 1.37.2.19 image_printf(&dec->cur, dec->edged_width, dec->height, 16, 16,
1994 :     "warning: nothing to output");
1995 :     image_printf(&dec->cur, dec->edged_width, dec->height, 16, 64,
1996 :     "bframe decoder lag");
1997 :    
1998 : suxen_drol 1.37.2.21 decoder_output(dec, &dec->cur, NULL, frame, 1 /*disable pp*/);
1999 : suxen_drol 1.37.2.19 }
2000 :     }
2001 :    
2002 : Isibaar 1.37.2.11 frame->length = BitstreamPos(&bs) / 8;
2003 : suxen_drol 1.37.2.8
2004 :     if (stats)
2005 :     {
2006 : suxen_drol 1.37.2.19 stats->notify = output ? XVID_DEC_VOP : XVID_DEC_NOTHING;
2007 : suxen_drol 1.37.2.8 stats->data.vop.time_base = (int)dec->time_base;
2008 :     stats->data.vop.time_increment = 0; //XXX: todo
2009 : chenm001 1.14 }
2010 : suxen_drol 1.37.2.21
2011 : Isibaar 1.1 emms();
2012 :    
2013 :     stop_global_timer();
2014 :    
2015 :     return XVID_ERR_OK;
2016 :     }
No admin address has been configured
 ViewVC Help
Powered by ViewVC 1.0.4