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

Diff of /xvidcore/src/decoder.c

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.37.2.5, Fri Oct 11 15:07:32 2002 UTC revision 1.37.2.28, Sat Jan 11 20:37:46 2003 UTC
# Line 59  Line 59 
59   *   *
60   *************************************************************************/   *************************************************************************/
61    
62    #include <stdio.h>
63  #include <stdlib.h>  #include <stdlib.h>
64  #include <string.h>  #include <string.h>
65    
# Line 68  Line 69 
69    
70  #include "xvid.h"  #include "xvid.h"
71  #include "portab.h"  #include "portab.h"
72    #include "global.h"
73    
74  #include "decoder.h"  #include "decoder.h"
75  #include "bitstream/bitstream.h"  #include "bitstream/bitstream.h"
# Line 79  Line 81 
81  #include "dct/fdct.h"  #include "dct/fdct.h"
82  #include "utils/mem_transfer.h"  #include "utils/mem_transfer.h"
83  #include "image/interpolate8x8.h"  #include "image/interpolate8x8.h"
84    #include "image/reduced.h"
85    #include "image/font.h"
86    
87  #include "bitstream/mbcoding.h"  #include "bitstream/mbcoding.h"
88  #include "prediction/mbprediction.h"  #include "prediction/mbprediction.h"
89  #include "utils/timer.h"  #include "utils/timer.h"
90  #include "utils/emms.h"  #include "utils/emms.h"
91    #include "motion/motion.h"
92    
93  #include "image/image.h"  #include "image/image.h"
94  #include "image/colorspace.h"  #include "image/colorspace.h"
95  #include "utils/mem_align.h"  #include "utils/mem_align.h"
96    
97  int  int
98  decoder_create(XVID_DEC_PARAM * param)  decoder_resize(DECODER * dec)
99  {  {
100          DECODER *dec;          /* free existing */
101    
102          dec = xvid_malloc(sizeof(DECODER), CACHE_LINE);          image_destroy(&dec->cur, dec->edged_width, dec->edged_height);
103          if (dec == NULL) {          image_destroy(&dec->refn[0], dec->edged_width, dec->edged_height);
104                  return XVID_ERR_MEMORY;          image_destroy(&dec->refn[1], dec->edged_width, dec->edged_height);
105          }          image_destroy(&dec->tmp, dec->edged_width, dec->edged_height);
106          param->handle = dec;          image_destroy(&dec->qtmp, dec->edged_width, dec->edged_height);
107    
108          dec->width = param->width;          image_destroy(&dec->gmc, dec->edged_width, dec->edged_height);
109          dec->height = param->height;  
110            if (dec->last_mbs)
111                    xvid_free(dec->last_mbs);
112            if (dec->mbs)
113                    xvid_free(dec->mbs);
114    
115            /* realloc */
116    
117          dec->mb_width = (dec->width + 15) / 16;          dec->mb_width = (dec->width + 15) / 16;
118          dec->mb_height = (dec->height + 15) / 16;          dec->mb_height = (dec->height + 15) / 16;
119    
120          dec->edged_width = 16 * dec->mb_width + 2 * EDGE_SIZE;          dec->edged_width = 16 * dec->mb_width + 2 * EDGE_SIZE;
121          dec->edged_height = 16 * dec->mb_height + 2 * EDGE_SIZE;          dec->edged_height = 16 * dec->mb_height + 2 * EDGE_SIZE;
         dec->low_delay = 0;  
122    
123          if (image_create(&dec->cur, dec->edged_width, dec->edged_height)) {          if (image_create(&dec->cur, dec->edged_width, dec->edged_height)) {
124                  xvid_free(dec);                  xvid_free(dec);
# Line 120  Line 130 
130                  xvid_free(dec);                  xvid_free(dec);
131                  return XVID_ERR_MEMORY;                  return XVID_ERR_MEMORY;
132          }          }
133    
134          // add by chenm001 <chenm001@163.com>          // add by chenm001 <chenm001@163.com>
135          // for support B-frame to reference last 2 frame          // for support B-frame to reference last 2 frame
136          if (image_create(&dec->refn[1], dec->edged_width, dec->edged_height)) {          if (image_create(&dec->refn[1], dec->edged_width, dec->edged_height)) {
# Line 128  Line 139 
139                  xvid_free(dec);                  xvid_free(dec);
140                  return XVID_ERR_MEMORY;                  return XVID_ERR_MEMORY;
141          }          }
142          if (image_create(&dec->refn[2], dec->edged_width, dec->edged_height)) {          if (image_create(&dec->tmp, dec->edged_width, dec->edged_height)) {
143                    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    
150            if (image_create(&dec->qtmp, dec->edged_width, dec->edged_height)) {
151                  image_destroy(&dec->cur, dec->edged_width, dec->edged_height);                  image_destroy(&dec->cur, dec->edged_width, dec->edged_height);
152                  image_destroy(&dec->refn[0], dec->edged_width, dec->edged_height);                  image_destroy(&dec->refn[0], dec->edged_width, dec->edged_height);
153                  image_destroy(&dec->refn[1], dec->edged_width, dec->edged_height);                  image_destroy(&dec->refn[1], dec->edged_width, dec->edged_height);
154                    image_destroy(&dec->tmp, dec->edged_width, dec->edged_height);
155                  xvid_free(dec);                  xvid_free(dec);
156                  return XVID_ERR_MEMORY;                  return XVID_ERR_MEMORY;
157          }          }
158    
159          if (image_create(&dec->refh, dec->edged_width, dec->edged_height)) {          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);                  image_destroy(&dec->cur, dec->edged_width, dec->edged_height);
162                  image_destroy(&dec->refn[0], dec->edged_width, dec->edged_height);                  image_destroy(&dec->refn[0], dec->edged_width, dec->edged_height);
163                  image_destroy(&dec->refn[1], dec->edged_width, dec->edged_height);                  image_destroy(&dec->refn[1], dec->edged_width, dec->edged_height);
164                  image_destroy(&dec->refn[2], dec->edged_width, dec->edged_height);                  image_destroy(&dec->tmp, dec->edged_width, dec->edged_height);
165                  xvid_free(dec);                  xvid_free(dec);
166                  return XVID_ERR_MEMORY;                  return XVID_ERR_MEMORY;
167          }          }
# Line 152  Line 173 
173                  image_destroy(&dec->cur, dec->edged_width, dec->edged_height);                  image_destroy(&dec->cur, dec->edged_width, dec->edged_height);
174                  image_destroy(&dec->refn[0], dec->edged_width, dec->edged_height);                  image_destroy(&dec->refn[0], dec->edged_width, dec->edged_height);
175                  image_destroy(&dec->refn[1], dec->edged_width, dec->edged_height);                  image_destroy(&dec->refn[1], dec->edged_width, dec->edged_height);
176                  image_destroy(&dec->refn[2], dec->edged_width, dec->edged_height);                  image_destroy(&dec->tmp, dec->edged_width, dec->edged_height);
177                  image_destroy(&dec->refh, dec->edged_width, dec->edged_height);                  image_destroy(&dec->qtmp, dec->edged_width, dec->edged_height);
178                  xvid_free(dec);                  xvid_free(dec);
179                  return XVID_ERR_MEMORY;                  return XVID_ERR_MEMORY;
180          }          }
   
181          memset(dec->mbs, 0, sizeof(MACROBLOCK) * dec->mb_width * dec->mb_height);          memset(dec->mbs, 0, sizeof(MACROBLOCK) * dec->mb_width * dec->mb_height);
182    
183          // add by chenm001 <chenm001@163.com>          // add by chenm001 <chenm001@163.com>
# Line 170  Line 190 
190                  image_destroy(&dec->cur, dec->edged_width, dec->edged_height);                  image_destroy(&dec->cur, dec->edged_width, dec->edged_height);
191                  image_destroy(&dec->refn[0], dec->edged_width, dec->edged_height);                  image_destroy(&dec->refn[0], dec->edged_width, dec->edged_height);
192                  image_destroy(&dec->refn[1], dec->edged_width, dec->edged_height);                  image_destroy(&dec->refn[1], dec->edged_width, dec->edged_height);
193                  image_destroy(&dec->refn[2], dec->edged_width, dec->edged_height);                  image_destroy(&dec->tmp, dec->edged_width, dec->edged_height);
194                  image_destroy(&dec->refh, dec->edged_width, dec->edged_height);                  image_destroy(&dec->qtmp, dec->edged_width, dec->edged_height);
195                  xvid_free(dec);                  xvid_free(dec);
196                  return XVID_ERR_MEMORY;                  return XVID_ERR_MEMORY;
197          }          }
198    
199          memset(dec->last_mbs, 0, sizeof(MACROBLOCK) * dec->mb_width * dec->mb_height);          memset(dec->last_mbs, 0, sizeof(MACROBLOCK) * dec->mb_width * dec->mb_height);
200    
201            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            image_null(&dec->tmp);
225            image_null(&dec->qtmp);
226    
227    /* image based GMC */
228            image_null(&dec->gmc);
229    
230    
231            dec->mbs = NULL;
232            dec->last_mbs = NULL;
233    
234          init_timer();          init_timer();
235    
236          // add by chenm001 <chenm001@163.com>          // add by chenm001 <chenm001@163.com>
237          // for support B-frame to save reference frame's time          // for support B-frame to save reference frame's time
238          dec->frames = -1;          dec->frames = 0;
239          dec->time = dec->time_base = dec->last_time_base = 0;          dec->time = dec->time_base = dec->last_time_base = 0;
240            dec->low_delay = 0;
241            dec->packed_mode = 0;
242    
243            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;          return XVID_ERR_OK;
249  }  }
250    
# Line 194  Line 254 
254  {  {
255          xvid_free(dec->last_mbs);          xvid_free(dec->last_mbs);
256          xvid_free(dec->mbs);          xvid_free(dec->mbs);
257    
258            image_destroy(&dec->gmc, dec->edged_width, dec->edged_height);          /* image based GMC */
259    
260          image_destroy(&dec->refn[0], dec->edged_width, dec->edged_height);          image_destroy(&dec->refn[0], dec->edged_width, dec->edged_height);
261          image_destroy(&dec->refn[1], dec->edged_width, dec->edged_height);          image_destroy(&dec->refn[1], dec->edged_width, dec->edged_height);
262          image_destroy(&dec->refn[2], dec->edged_width, dec->edged_height);          image_destroy(&dec->tmp, dec->edged_width, dec->edged_height);
263          image_destroy(&dec->refh, dec->edged_width, dec->edged_height);          image_destroy(&dec->qtmp, dec->edged_width, dec->edged_height);
264          image_destroy(&dec->cur, dec->edged_width, dec->edged_height);          image_destroy(&dec->cur, dec->edged_width, dec->edged_height);
265          xvid_free(dec);          xvid_free(dec);
266    
# Line 226  Line 289 
289                                  Bitstream * bs,                                  Bitstream * bs,
290                                  const uint32_t quant,                                  const uint32_t quant,
291                                  const uint32_t intra_dc_threshold,                                  const uint32_t intra_dc_threshold,
292                                  const unsigned int bound)                                  const unsigned int bound,
293                                    const int reduced_resolution)
294  {  {
295    
296          DECLARE_ALIGNED_MATRIX(block, 6, 64, int16_t, CACHE_LINE);          DECLARE_ALIGNED_MATRIX(block, 6, 64, int16_t, CACHE_LINE);
# Line 239  Line 303 
303          uint32_t iQuant = pMB->quant;          uint32_t iQuant = pMB->quant;
304          uint8_t *pY_Cur, *pU_Cur, *pV_Cur;          uint8_t *pY_Cur, *pU_Cur, *pV_Cur;
305    
306            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);          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);          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);          pV_Cur = dec->cur.v + (y_pos << 3) * stride2 + (x_pos << 3);
314            }
315    
316          memset(block, 0, 6 * 64 * sizeof(int16_t));     // clear          memset(block, 0, 6 * 64 * sizeof(int16_t));     // clear
317    
# Line 302  Line 372 
372                  start_timer();                  start_timer();
373                  idct(&data[i * 64]);                  idct(&data[i * 64]);
374                  stop_idct_timer();                  stop_idct_timer();
375    
376          }          }
377    
378          if (dec->interlacing && pMB->field_dct) {          if (dec->interlacing && pMB->field_dct) {
# Line 310  Line 381 
381          }          }
382    
383          start_timer();          start_timer();
384    
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);          transfer_16to8copy(pY_Cur, &data[0 * 64], stride);
396          transfer_16to8copy(pY_Cur + 8, &data[1 * 64], stride);          transfer_16to8copy(pY_Cur + 8, &data[1 * 64], stride);
397          transfer_16to8copy(pY_Cur + next_block, &data[2 * 64], stride);          transfer_16to8copy(pY_Cur + next_block, &data[2 * 64], stride);
398          transfer_16to8copy(pY_Cur + 8 + next_block, &data[3 * 64], stride);          transfer_16to8copy(pY_Cur + 8 + next_block, &data[3 * 64], stride);
399          transfer_16to8copy(pU_Cur, &data[4 * 64], stride2);          transfer_16to8copy(pU_Cur, &data[4 * 64], stride2);
400          transfer_16to8copy(pV_Cur, &data[5 * 64], stride2);          transfer_16to8copy(pV_Cur, &data[5 * 64], stride2);
401            }
402          stop_transfer_timer();          stop_transfer_timer();
403  }  }
404    
405    
406    
407    
   
 #define SIGN(X) (((X)>0)?1:-1)  
 #define ABS(X) (((X)>0)?(X):-(X))  
 static const uint32_t roundtab[16] =  
         { 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2 };  
   
   
408  // decode an inter macroblock  // decode an inter macroblock
409    
410  void  void
# Line 340  Line 416 
416                                  const uint32_t cbp,                                  const uint32_t cbp,
417                                  Bitstream * bs,                                  Bitstream * bs,
418                                  const uint32_t quant,                                  const uint32_t quant,
419                                  const uint32_t rounding)                                  const uint32_t rounding,
420                                    const int reduced_resolution)
421  {  {
422    
423          DECLARE_ALIGNED_MATRIX(block, 6, 64, int16_t, CACHE_LINE);          DECLARE_ALIGNED_MATRIX(block, 6, 64, int16_t, CACHE_LINE);
# Line 348  Line 425 
425    
426          uint32_t stride = dec->edged_width;          uint32_t stride = dec->edged_width;
427          uint32_t stride2 = stride / 2;          uint32_t stride2 = stride / 2;
428          uint32_t next_block = stride * 8;          uint32_t next_block = stride * (reduced_resolution ? 16 : 8);
429          uint32_t i;          uint32_t i;
430          uint32_t iQuant = pMB->quant;          uint32_t iQuant = pMB->quant;
431          uint8_t *pY_Cur, *pU_Cur, *pV_Cur;          uint8_t *pY_Cur, *pU_Cur, *pV_Cur;
432    
433          int uv_dx, uv_dy;          int uv_dx, uv_dy;
434            VECTOR mv[4];   /* local copy of mvs */
435    
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                    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            } else {
445          pY_Cur = dec->cur.y + (y_pos << 4) * stride + (x_pos << 4);          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);          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);          pV_Cur = dec->cur.v + (y_pos << 3) * stride2 + (x_pos << 3);
448                    for (i = 0; i < 4; i++)
449                            mv[i] = pMB->mvs[i];
450            }
451    
452          if (pMB->mode == MODE_INTER || pMB->mode == MODE_INTER_Q) {          if (pMB->mode == MODE_INTER || pMB->mode == MODE_INTER_Q) {
                 uv_dx = pMB->mvs[0].x;  
                 uv_dy = pMB->mvs[0].y;  
453    
454                  if (dec->quarterpel)                  uv_dx = mv[0].x / (1 + dec->quarterpel);
455                  {                  uv_dy = mv[0].y / (1 + dec->quarterpel);
                         uv_dx = (uv_dx >> 1) | (uv_dx & 1);  
                         uv_dy = (uv_dy >> 1) | (uv_dy & 1);  
                 }  
456    
457                  uv_dx = (uv_dx & 3) ? (uv_dx >> 1) | 1 : uv_dx / 2;                  uv_dx = (uv_dx >> 1) + roundtab_79[uv_dx & 0x3];
458                  uv_dy = (uv_dy & 3) ? (uv_dy >> 1) | 1 : uv_dy / 2;                  uv_dy = (uv_dy >> 1) + roundtab_79[uv_dy & 0x3];
459    
460                  start_timer();                  start_timer();
461                    if (reduced_resolution)
462                    {
463                            interpolate32x32_switch(dec->cur.y, dec->refn[0].y, 32*x_pos, 32*y_pos,
464                                                                      mv[0].x, mv[0].y, stride,  rounding);
465                            interpolate16x16_switch(dec->cur.u, dec->refn[0].u, 16 * x_pos, 16 * y_pos,
466                                                                      uv_dx, uv_dy, stride2, rounding);
467                            interpolate16x16_switch(dec->cur.v, dec->refn[0].v, 16 * x_pos, 16 * y_pos,
468                                                                      uv_dx, uv_dy, stride2, rounding);
469    
470                    }
471                    else
472                    {
473                  if(dec->quarterpel) {                  if(dec->quarterpel) {
474                          interpolate16x16_quarterpel(dec->cur.y, dec->refn[0].y, dec->refh.y, dec->refh.y + 64,                                  interpolate16x16_quarterpel(dec->cur.y, dec->refn[0].y, dec->qtmp.y, dec->qtmp.y + 64,
475                                                                              dec->refh.y + 128, 16*x_pos, 16*y_pos,                                                                                          dec->qtmp.y + 128, 16*x_pos, 16*y_pos,
476                                                                              pMB->mvs[0].x, pMB->mvs[0].y, stride,  rounding);                                                                                          mv[0].x, mv[0].y, stride, rounding);
477                  }                  }
478                  else {                  else {
479                          interpolate8x8_switch(dec->cur.y, dec->refn[0].y, 16*x_pos, 16*y_pos,                                  interpolate16x16_switch(dec->cur.y, dec->refn[0].y, 16*x_pos, 16*y_pos,
480                                                                    pMB->mvs[0].x, pMB->mvs[0].y, stride,  rounding);                                                                            mv[0].x, mv[0].y, stride, rounding);
                         interpolate8x8_switch(dec->cur.y, dec->refn[0].y, 16*x_pos + 8, 16*y_pos,  
                                                               pMB->mvs[1].x, pMB->mvs[1].y, stride,  rounding);  
                         interpolate8x8_switch(dec->cur.y, dec->refn[0].y, 16*x_pos, 16*y_pos + 8,  
                                                                   pMB->mvs[2].x, pMB->mvs[2].y, stride,  rounding);  
                         interpolate8x8_switch(dec->cur.y, dec->refn[0].y, 16*x_pos + 8, 16*y_pos + 8,  
                                                                   pMB->mvs[3].x, pMB->mvs[3].y, stride,  rounding);  
481                  }                  }
482    
483                  interpolate8x8_switch(dec->cur.u, dec->refn[0].u, 8 * x_pos, 8 * y_pos,                  interpolate8x8_switch(dec->cur.u, dec->refn[0].u, 8 * x_pos, 8 * y_pos,
484                                                            uv_dx, uv_dy, stride2, rounding);                                                            uv_dx, uv_dy, stride2, rounding);
485                  interpolate8x8_switch(dec->cur.v, dec->refn[0].v, 8 * x_pos, 8 * y_pos,                  interpolate8x8_switch(dec->cur.v, dec->refn[0].v, 8 * x_pos, 8 * y_pos,
486                                                            uv_dx, uv_dy, stride2, rounding);                                                            uv_dx, uv_dy, stride2, rounding);
487                    }
488                  stop_comp_timer();                  stop_comp_timer();
489    
490          } else {          } else {        /* MODE_INTER4V */
491                  int sum;                  int sum;
                 sum = pMB->mvs[0].x + pMB->mvs[1].x + pMB->mvs[2].x + pMB->mvs[3].x;  
492    
493                  if (dec->quarterpel)                  if (dec->quarterpel)
494                  {                          sum = (mv[0].x / 2) + (mv[1].x / 2) + (mv[2].x / 2) + (mv[3].x / 2);
495                          sum /= 2;                  else
496                  }                          sum = mv[0].x + mv[1].x + mv[2].x + mv[3].x;
   
                 uv_dx = (sum == 0 ? 0 : SIGN(sum) * (roundtab[ABS(sum) % 16] + (ABS(sum) / 16) * 2));  
497    
498                  sum = pMB->mvs[0].y + pMB->mvs[1].y + pMB->mvs[2].y + pMB->mvs[3].y;                  uv_dx = (sum >> 3) + roundtab_76[sum & 0xf];
499    
500                  if (dec->quarterpel)                  if (dec->quarterpel)
501                  {                          sum = (mv[0].y / 2) + (mv[1].y / 2) + (mv[2].y / 2) + (mv[3].y / 2);
502                          sum /= 2;                  else
503                  }                          sum = mv[0].y + mv[1].y + mv[2].y + mv[3].y;
504    
505                  uv_dy = (sum == 0 ? 0 : SIGN(sum) * (roundtab[ABS(sum) % 16] + (ABS(sum) / 16) * 2));                  uv_dy = (sum >> 3) + roundtab_76[sum & 0xf];
506    
507                  start_timer();                  start_timer();
508                    if (reduced_resolution)
509                    {
510                            interpolate16x16_switch(dec->cur.y, dec->refn[0].y, 32*x_pos, 32*y_pos,
511                                                                      mv[0].x, mv[0].y, stride,  rounding);
512                            interpolate16x16_switch(dec->cur.y, dec->refn[0].y , 32*x_pos + 16, 32*y_pos,
513                                                                      mv[1].x, mv[1].y, stride,  rounding);
514                            interpolate16x16_switch(dec->cur.y, dec->refn[0].y , 32*x_pos, 32*y_pos + 16,
515                                                                      mv[2].x, mv[2].y, stride,  rounding);
516                            interpolate16x16_switch(dec->cur.y, dec->refn[0].y , 32*x_pos + 16, 32*y_pos + 16,
517                                                                      mv[3].x, mv[3].y, stride,  rounding);
518                            interpolate16x16_switch(dec->cur.u, dec->refn[0].u , 16 * x_pos, 16 * y_pos,
519                                                                      uv_dx, uv_dy, stride2, rounding);
520                            interpolate16x16_switch(dec->cur.v, dec->refn[0].v , 16 * x_pos, 16 * y_pos,
521                                                                      uv_dx, uv_dy, stride2, rounding);
522    
523                            // set_block(pY_Cur, stride, 32, 32, 127);
524                    }
525                    else
526                    {
527                  if(dec->quarterpel) {                  if(dec->quarterpel) {
528                          interpolate8x8_quarterpel(dec->cur.y, dec->refn[0].y, dec->refh.y, dec->refh.y + 64,                                  interpolate8x8_quarterpel(dec->cur.y, dec->refn[0].y , dec->qtmp.y, dec->qtmp.y + 64,
529                                                                            dec->refh.y + 128, 16*x_pos, 16*y_pos,                                                                                    dec->qtmp.y + 128, 16*x_pos, 16*y_pos,
530                                                                            pMB->mvs[0].x, pMB->mvs[0].y, stride,  rounding);                                                                                    mv[0].x, mv[0].y, stride,  rounding);
531                          interpolate8x8_quarterpel(dec->cur.y, dec->refn[0].y, dec->refh.y, dec->refh.y + 64,                                  interpolate8x8_quarterpel(dec->cur.y, dec->refn[0].y , dec->qtmp.y, dec->qtmp.y + 64,
532                                                                            dec->refh.y + 128, 16*x_pos + 8, 16*y_pos,                                                                                    dec->qtmp.y + 128, 16*x_pos + 8, 16*y_pos,
533                                                                            pMB->mvs[1].x, pMB->mvs[1].y, stride,  rounding);                                                                                    mv[1].x, mv[1].y, stride,  rounding);
534                          interpolate8x8_quarterpel(dec->cur.y, dec->refn[0].y, dec->refh.y, dec->refh.y + 64,                                  interpolate8x8_quarterpel(dec->cur.y, dec->refn[0].y , dec->qtmp.y, dec->qtmp.y + 64,
535                                                                            dec->refh.y + 128, 16*x_pos, 16*y_pos + 8,                                                                                    dec->qtmp.y + 128, 16*x_pos, 16*y_pos + 8,
536                                                                            pMB->mvs[2].x, pMB->mvs[2].y, stride,  rounding);                                                                                    mv[2].x, mv[2].y, stride,  rounding);
537                          interpolate8x8_quarterpel(dec->cur.y, dec->refn[0].y, dec->refh.y, dec->refh.y + 64,                                  interpolate8x8_quarterpel(dec->cur.y, dec->refn[0].y , dec->qtmp.y, dec->qtmp.y + 64,
538                                                                            dec->refh.y + 128, 16*x_pos + 8, 16*y_pos + 8,                                                                                    dec->qtmp.y + 128, 16*x_pos + 8, 16*y_pos + 8,
539                                                                            pMB->mvs[3].x, pMB->mvs[3].y, stride,  rounding);                                                                                    mv[3].x, mv[3].y, stride,  rounding);
540                  }                  }
541                  else {                  else {
542                          interpolate8x8_switch(dec->cur.y, dec->refn[0].y, 16*x_pos, 16*y_pos,                          interpolate8x8_switch(dec->cur.y, dec->refn[0].y, 16*x_pos, 16*y_pos,
543                                                                    pMB->mvs[0].x, pMB->mvs[0].y, stride,  rounding);                                                                            mv[0].x, mv[0].y, stride,  rounding);
544                          interpolate8x8_switch(dec->cur.y, dec->refn[0].y, 16*x_pos + 8, 16*y_pos,                          interpolate8x8_switch(dec->cur.y, dec->refn[0].y, 16*x_pos + 8, 16*y_pos,
545                                                                    pMB->mvs[1].x, pMB->mvs[1].y, stride,  rounding);                                                                            mv[1].x, mv[1].y, stride,  rounding);
546                          interpolate8x8_switch(dec->cur.y, dec->refn[0].y, 16*x_pos, 16*y_pos + 8,                          interpolate8x8_switch(dec->cur.y, dec->refn[0].y, 16*x_pos, 16*y_pos + 8,
547                                                                    pMB->mvs[2].x, pMB->mvs[2].y, stride,  rounding);                                                                            mv[2].x, mv[2].y, stride,  rounding);
548                          interpolate8x8_switch(dec->cur.y, dec->refn[0].y, 16*x_pos + 8, 16*y_pos + 8,                          interpolate8x8_switch(dec->cur.y, dec->refn[0].y, 16*x_pos + 8, 16*y_pos + 8,
549                                                                    pMB->mvs[3].x, pMB->mvs[3].y, stride,  rounding);                                                                            mv[3].x, mv[3].y, stride,  rounding);
550                  }                  }
551    
552                  interpolate8x8_switch(dec->cur.u, dec->refn[0].u, 8 * x_pos, 8 * y_pos,                  interpolate8x8_switch(dec->cur.u, dec->refn[0].u, 8 * x_pos, 8 * y_pos,
553                                                            uv_dx, uv_dy, stride2, rounding);                                                            uv_dx, uv_dy, stride2, rounding);
554                  interpolate8x8_switch(dec->cur.v, dec->refn[0].v, 8 * x_pos, 8 * y_pos,                  interpolate8x8_switch(dec->cur.v, dec->refn[0].v, 8 * x_pos, 8 * y_pos,
555                                                            uv_dx, uv_dy, stride2, rounding);                                                            uv_dx, uv_dy, stride2, rounding);
556                    }
557                  stop_comp_timer();                  stop_comp_timer();
558          }          }
559    
# Line 478  Line 588 
588          }          }
589    
590          start_timer();          start_timer();
591            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            stop_transfer_timer();
622    }
623    
624    
625    static void
626    decoder_mbgmc(DECODER * dec,
627                                    MACROBLOCK * const pMB,
628                                    const uint32_t x_pos,
629                                    const uint32_t y_pos,
630                                    const uint32_t acpred_flag,
631                                    const uint32_t cbp,
632                                    Bitstream * bs,
633                                    const uint32_t quant,
634                                    const uint32_t rounding,
635                                    const int reduced_resolution)   /* no reduced res support */
636    {
637    
638            DECLARE_ALIGNED_MATRIX(block, 6, 64, int16_t, CACHE_LINE);
639            DECLARE_ALIGNED_MATRIX(data, 6, 64, int16_t, CACHE_LINE);
640    
641            const uint32_t stride = dec->edged_width;
642            const uint32_t stride2 = stride / 2;
643            const uint32_t next_block = stride * (reduced_resolution ? 16 : 8);
644            uint32_t i;
645            const uint32_t iQuant = pMB->quant;
646            uint8_t *const pY_Cur=dec->cur.y + (y_pos << 4) * stride + (x_pos << 4);
647            uint8_t *const pU_Cur=dec->cur.u + (y_pos << 3) * stride2 + (x_pos << 3);
648            uint8_t *const pV_Cur=dec->cur.v + (y_pos << 3) * stride2 + (x_pos << 3);
649    
650            pMB->mvs[0] = pMB->mvs[1] = pMB->mvs[2] = pMB->mvs[3] = pMB->amv;
651    
652            start_timer();
653            transfer16x16_copy(pY_Cur, dec->gmc.y + (y_pos << 4)*stride + (x_pos  << 4), stride);
654            transfer8x8_copy(pU_Cur, dec->gmc.u + (y_pos << 3)*stride2 + (x_pos  << 3), stride2);
655            transfer8x8_copy(pV_Cur, dec->gmc.v + (y_pos << 3)*stride2 + (x_pos << 3), stride2);
656            stop_transfer_timer();
657    
658            if (!cbp) return;
659    
660            for (i = 0; i < 6; i++) {
661                    int direction = dec->alternate_vertical_scan ? 2 : 0;
662    
663                    if (cbp & (1 << (5 - i)))       // coded
664                    {
665                            memset(&block[i * 64], 0, 64 * sizeof(int16_t));        // clear
666    
667                            start_timer();
668                            get_inter_block(bs, &block[i * 64], direction);
669                            stop_coding_timer();
670    
671                            start_timer();
672                            if (dec->quant_type == 0) {
673                                    dequant_inter(&data[i * 64], &block[i * 64], iQuant);
674                            } else {
675                                    dequant4_inter(&data[i * 64], &block[i * 64], iQuant);
676                            }
677                            stop_iquant_timer();
678    
679                            start_timer();
680                            idct(&data[i * 64]);
681                            stop_idct_timer();
682                    }
683            }
684    
685    /* interlace + GMC is this possible ??? */
686    /*      if (dec->interlacing && pMB->field_dct) {
687                    next_block = stride;
688                    stride *= 2;
689            }
690    */
691            start_timer();
692          if (cbp & 32)          if (cbp & 32)
693                  transfer_16to8add(pY_Cur, &data[0 * 64], stride);                  transfer_16to8add(pY_Cur, &data[0 * 64], stride);
694          if (cbp & 16)          if (cbp & 16)
# Line 497  Line 708 
708  void  void
709  decoder_iframe(DECODER * dec,  decoder_iframe(DECODER * dec,
710                             Bitstream * bs,                             Bitstream * bs,
711                               int reduced_resolution,
712                             int quant,                             int quant,
713                             int intra_dc_threshold)                             int intra_dc_threshold)
714  {  {
715          uint32_t bound;          uint32_t bound;
716          uint32_t x, y;          uint32_t x, y;
717            uint32_t mb_width = dec->mb_width;
718            uint32_t mb_height = dec->mb_height;
719    
720            if (reduced_resolution)
721            {
722                    mb_width = (dec->width + 31) / 32;
723                    mb_height = (dec->height + 31) / 32;
724            }
725    
726          bound = 0;          bound = 0;
727    
728          for (y = 0; y < dec->mb_height; y++) {          for (y = 0; y < mb_height; y++) {
729                  for (x = 0; x < dec->mb_width; x++) {                  for (x = 0; x < mb_width; x++) {
730                          MACROBLOCK *mb;                          MACROBLOCK *mb;
731                          uint32_t mcbpc;                          uint32_t mcbpc;
732                          uint32_t cbpc;                          uint32_t cbpc;
# Line 519  Line 739 
739    
740                          if (check_resync_marker(bs, 0))                          if (check_resync_marker(bs, 0))
741                          {                          {
742                                  bound = read_video_packet_header(bs, 0, &quant);                                  bound = read_video_packet_header(bs, dec, 0,
743                                  x = bound % dec->mb_width;                                                          &quant, NULL, NULL, &intra_dc_threshold);
744                                  y = bound / dec->mb_width;                                  x = bound % mb_width;
745                                    y = bound / mb_width;
746                          }                          }
747                          mb = &dec->mbs[y * dec->mb_width + x];                          mb = &dec->mbs[y * dec->mb_width + x];
748    
# Line 552  Line 773 
773    
774                          if (dec->interlacing) {                          if (dec->interlacing) {
775                                  mb->field_dct = BitstreamGetBit(bs);                                  mb->field_dct = BitstreamGetBit(bs);
776                                  DEBUG1("deci: field_dct: ", mb->field_dct);                                  DPRINTF(DPRINTF_MB,"deci: field_dct: %i", mb->field_dct);
777                          }                          }
778    
779                          decoder_mbintra(dec, mb, x, y, acpred_flag, cbp, bs, quant,                          decoder_mbintra(dec, mb, x, y, acpred_flag, cbp, bs, quant,
780                                                          intra_dc_threshold, bound);                                                          intra_dc_threshold, bound, reduced_resolution);
781    
782                  }                  }
783                  if(dec->out_frm)                  if(dec->out_frm)
784                    output_slice(&dec->cur, dec->edged_width,dec->width,dec->out_frm,0,y,dec->mb_width);                    output_slice(&dec->cur, dec->edged_width,dec->width,dec->out_frm,0,y,mb_width);
   
785          }          }
786    
787  }  }
# Line 572  Line 793 
793                                    int x,                                    int x,
794                                    int y,                                    int y,
795                                    int k,                                    int k,
796                                    VECTOR * mv,                                    VECTOR * ret_mv,
797                                    int fcode,                                    int fcode,
798                                    const int bound)                                    const int bound)
799  {  {
# Line 583  Line 804 
804          int range = (64 * scale_fac);          int range = (64 * scale_fac);
805    
806          VECTOR pmv;          VECTOR pmv;
807          int mv_x, mv_y;          VECTOR mv;
808    
809          pmv = get_pmv2(dec->mbs, dec->mb_width, bound, x, y, k);          pmv = get_pmv2(dec->mbs, dec->mb_width, bound, x, y, k);
810    
811          mv_x = get_mv(bs, fcode);          mv.x = get_mv(bs, fcode);
812          mv_y = get_mv(bs, fcode);          mv.y = get_mv(bs, fcode);
813    
814          DPRINTF(DPRINTF_MV,"mv_diff (%i,%i) pred (%i,%i)", mv_x, mv_y, pmv.x, pmv.y);          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);
815    
816          mv_x += pmv.x;          mv.x += pmv.x;
817          mv_y += pmv.y;          mv.y += pmv.y;
818    
819          if (mv_x < low) {          if (mv.x < low) {
820                  mv_x += range;                  mv.x += range;
821          } else if (mv_x > high) {          } else if (mv.x > high) {
822                  mv_x -= range;                  mv.x -= range;
823          }          }
824    
825          if (mv_y < low) {          if (mv.y < low) {
826                  mv_y += range;                  mv.y += range;
827          } else if (mv_y > high) {          } else if (mv.y > high) {
828                  mv_y -= range;                  mv.y -= range;
829          }          }
830    
831          mv->x = mv_x;          ret_mv->x = mv.x;
832          mv->y = mv_y;          ret_mv->y = mv.y;
833    }
834    
835    
836    
837    static __inline int gmc_sanitize(int value, int quarterpel, int fcode)
838    {
839            int length = 1 << (fcode+4);
840    
841            if (quarterpel) value *= 2;
842    
843            if (value < -length)
844                    return -length;
845            else if (value >= length)
846                    return length-1;
847            else return value;
848  }  }
849    
850    
851    /* for P_VOP set gmc_warp to NULL */
852  void  void
853  decoder_pframe(DECODER * dec,  decoder_pframe(DECODER * dec,
854                             Bitstream * bs,                             Bitstream * bs,
855                             int rounding,                             int rounding,
856                               int reduced_resolution,
857                             int quant,                             int quant,
858                             int fcode,                             int fcode,
859                             int intra_dc_threshold)                             int intra_dc_threshold,
860                               const WARPPOINTS *const gmc_warp)
861  {  {
862    
863          uint32_t x, y;          uint32_t x, y;
864          uint32_t bound;          uint32_t bound;
865          int cp_mb, st_mb;          int cp_mb, st_mb;
866            uint32_t mb_width = dec->mb_width;
867            uint32_t mb_height = dec->mb_height;
868    
869            static int framecount=0;
870            if (reduced_resolution)
871            {
872                    mb_width = (dec->width + 31) / 32;
873                    mb_height = (dec->height + 31) / 32;
874            }
875    
876          start_timer();          start_timer();
877          image_setedges(&dec->refn[0], dec->edged_width, dec->edged_height,          image_setedges(&dec->refn[0], dec->edged_width, dec->edged_height,
878                                     dec->width, dec->height);                                     dec->width, dec->height);
879          stop_edges_timer();          stop_edges_timer();
880    
881            if (gmc_warp)
882            {
883    
884                    // accuracy:  0==1/2, 1=1/4, 2=1/8, 3=1/16
885                    if ( (dec->sprite_warping_accuracy != 3) || (dec->sprite_warping_points != 2) )
886                    {
887                            fprintf(stderr,"Wrong GMC parameters acc=%d(-> 1/%d), %d!!!\n",
888                                    dec->sprite_warping_accuracy,(2<<dec->sprite_warping_accuracy),
889                                    dec->sprite_warping_points);
890                    }
891    
892                    generate_GMCparameters( dec->sprite_warping_points,
893                                    (2 << dec->sprite_warping_accuracy), gmc_warp,
894                                    dec->width, dec->height, &dec->gmc_data);
895    
896                    generate_GMCimage(&dec->gmc_data, &dec->refn[0],
897                                            mb_width, mb_height,
898                                            dec->edged_width, dec->edged_width/2,
899                                            fcode, dec->quarterpel, 0,
900                                            rounding, dec->mbs, &dec->gmc);
901    
902            }
903    
904          bound = 0;          bound = 0;
905    
906          for (y = 0; y < dec->mb_height; y++) {          for (y = 0; y < mb_height; y++) {
907                  cp_mb = st_mb = 0;                  cp_mb = st_mb = 0;
908                  for (x = 0; x < dec->mb_width; x++) {                  for (x = 0; x < mb_width; x++) {
909                          MACROBLOCK *mb;                          MACROBLOCK *mb;
910    
911                          // skip stuffing                          // skip stuffing
# Line 644  Line 914 
914    
915                          if (check_resync_marker(bs, fcode - 1))                          if (check_resync_marker(bs, fcode - 1))
916                          {                          {
917                                  bound = read_video_packet_header(bs, fcode - 1, &quant);                                  bound = read_video_packet_header(bs, dec, fcode - 1,
918                                  x = bound % dec->mb_width;                                          &quant, &fcode, NULL, &intra_dc_threshold);
919                                  y = bound / dec->mb_width;                                  x = bound % mb_width;
920                                    y = bound / mb_width;
921                          }                          }
922                          mb = &dec->mbs[y * dec->mb_width + x];                          mb = &dec->mbs[y * dec->mb_width + x];
923    
924                          DPRINTF(DPRINTF_MB, "macroblock (%i,%i) %08x", x, y, BitstreamShowBits(bs, 32));                          DPRINTF(DPRINTF_MB, "macroblock (%i,%i) %08x", x, y, BitstreamShowBits(bs, 32));
925    
926                          //if (!(dec->mb_skip[y*dec->mb_width + x]=BitstreamGetBit(bs)))         // not_coded                          //if (!(dec->mb_skip[y*dec->mb_width + x]=BitstreamGetBit(bs)))         // not_coded
927                          if (!(BitstreamGetBit(bs)))     // not_coded                          if (!(BitstreamGetBit(bs)))     // block _is_ coded
928                          {                          {
929                                  uint32_t mcbpc;                                  uint32_t mcbpc;
930                                  uint32_t cbpc;                                  uint32_t cbpc;
# Line 661  Line 932 
932                                  uint32_t cbpy;                                  uint32_t cbpy;
933                                  uint32_t cbp;                                  uint32_t cbp;
934                                  uint32_t intra;                                  uint32_t intra;
935                                    int mcsel = 0;          // mcsel: '0'=local motion, '1'=GMC
936    
937                                  cp_mb++;                                  cp_mb++;
938                                  mcbpc = get_mcbpc_inter(bs);                                  mcbpc = get_mcbpc_inter(bs);
# Line 677  Line 949 
949                                          acpred_flag = BitstreamGetBit(bs);                                          acpred_flag = BitstreamGetBit(bs);
950                                  }                                  }
951    
952                                    if (gmc_warp && (mb->mode == MODE_INTER || mb->mode == MODE_INTER_Q))
953                                    {
954                                            mcsel = BitstreamGetBit(bs);
955                                    }
956    
957                                  cbpy = get_cbpy(bs, intra);                                  cbpy = get_cbpy(bs, intra);
958                                  DPRINTF(DPRINTF_MB, "cbpy %i", cbpy);                                  DPRINTF(DPRINTF_MB, "cbpy %i  mcsel %i ", cbpy,mcsel);
959    
960                                  cbp = (cbpy << 2) | cbpc;                                  cbp = (cbpy << 2) | cbpc;
961    
# Line 698  Line 975 
975                                  if (dec->interlacing) {                                  if (dec->interlacing) {
976                                          if (cbp || intra) {                                          if (cbp || intra) {
977                                                  mb->field_dct = BitstreamGetBit(bs);                                                  mb->field_dct = BitstreamGetBit(bs);
978                                                  DEBUG1("decp: field_dct: ", mb->field_dct);                                                  DPRINTF(DPRINTF_MB,"decp: field_dct: %i", mb->field_dct);
979                                          }                                          }
980    
981                                          if (mb->mode == MODE_INTER || mb->mode == MODE_INTER_Q) {                                          if (mb->mode == MODE_INTER || mb->mode == MODE_INTER_Q) {
982                                                  mb->field_pred = BitstreamGetBit(bs);                                                  mb->field_pred = BitstreamGetBit(bs);
983                                                  DEBUG1("decp: field_pred: ", mb->field_pred);                                                  DPRINTF(DPRINTF_MB, "decp: field_pred: %i", mb->field_pred);
984    
985                                                  if (mb->field_pred) {                                                  if (mb->field_pred) {
986                                                          mb->field_for_top = BitstreamGetBit(bs);                                                          mb->field_for_top = BitstreamGetBit(bs);
987                                                          DEBUG1("decp: field_for_top: ", mb->field_for_top);                                                          DPRINTF(DPRINTF_MB,"decp: field_for_top: %i", mb->field_for_top);
988                                                          mb->field_for_bot = BitstreamGetBit(bs);                                                          mb->field_for_bot = BitstreamGetBit(bs);
989                                                          DEBUG1("decp: field_for_bot: ", mb->field_for_bot);                                                          DPRINTF(DPRINTF_MB,"decp: field_for_bot: %i", mb->field_for_bot);
990                                                  }                                                  }
991                                          }                                          }
992                                  }                                  }
993    
994                                  if (mb->mode == MODE_INTER || mb->mode == MODE_INTER_Q) {                                  if (mcsel) {
995                                            decoder_mbgmc(dec, mb, x, y, 0, cbp, bs, quant,
996                                                                    rounding, reduced_resolution);
997                                            continue;
998    
999                                    } else if (mb->mode == MODE_INTER || mb->mode == MODE_INTER_Q) {
1000    
1001                                          if (dec->interlacing && mb->field_pred) {                                          if (dec->interlacing && mb->field_pred) {
1002                                                  get_motion_vector(dec, bs, x, y, 0, &mb->mvs[0],                                                  get_motion_vector(dec, bs, x, y, 0, &mb->mvs[0],
1003                                                                                    fcode, bound);                                                                                    fcode, bound);
# Line 723  Line 1006 
1006                                          } else {                                          } else {
1007                                                  get_motion_vector(dec, bs, x, y, 0, &mb->mvs[0],                                                  get_motion_vector(dec, bs, x, y, 0, &mb->mvs[0],
1008                                                                                    fcode, bound);                                                                                    fcode, bound);
1009                                                  mb->mvs[1].x = mb->mvs[2].x = mb->mvs[3].x =                                                  mb->mvs[1] = mb->mvs[2] = mb->mvs[3] = mb->mvs[0];
                                                         mb->mvs[0].x;  
                                                 mb->mvs[1].y = mb->mvs[2].y = mb->mvs[3].y =  
                                                         mb->mvs[0].y;  
1010                                          }                                          }
1011                                  } else if (mb->mode == MODE_INTER4V ) {                                  } else if (mb->mode == MODE_INTER4V ) {
1012    
# Line 741  Line 1021 
1021                                          mb->mvs[0].y = mb->mvs[1].y = mb->mvs[2].y = mb->mvs[3].y =                                          mb->mvs[0].y = mb->mvs[1].y = mb->mvs[2].y = mb->mvs[3].y =
1022                                                  0;                                                  0;
1023                                          decoder_mbintra(dec, mb, x, y, acpred_flag, cbp, bs, quant,                                          decoder_mbintra(dec, mb, x, y, acpred_flag, cbp, bs, quant,
1024                                                                          intra_dc_threshold, bound);                                                                          intra_dc_threshold, bound, reduced_resolution);
1025                                          continue;                                          continue;
1026                                  }                                  }
1027    
1028                                  decoder_mbinter(dec, mb, x, y, acpred_flag, cbp, bs, quant,                                  decoder_mbinter(dec, mb, x, y, 0, cbp, bs, quant,
1029                                                                  rounding);                                                                  rounding, reduced_resolution);
1030                          } else                          // not coded  
1031                            }
1032                            else if (gmc_warp)      /* a not coded S(GMC)-VOP macroblock */
1033                            {
1034                                    mb->mode = MODE_NOT_CODED_GMC;
1035    
1036                                    start_timer();
1037    
1038                                    decoder_mbgmc(dec, mb, x, y, 0, 0x00, bs, quant,
1039                                                                    rounding, reduced_resolution);
1040    
1041                                    stop_transfer_timer();
1042    
1043                                    if(dec->out_frm && cp_mb > 0) {
1044                                      output_slice(&dec->cur, dec->edged_width,dec->width,dec->out_frm,st_mb,y,cp_mb);
1045                                      cp_mb = 0;
1046                                    }
1047                                    st_mb = x+1;
1048                            }
1049                            else    /* not coded P_VOP macroblock */
1050                          {                          {
1051                                  mb->mode = MODE_NOT_CODED;                                  mb->mode = MODE_NOT_CODED;
1052    
1053                                  mb->mvs[0].x = mb->mvs[1].x = mb->mvs[2].x = mb->mvs[3].x = 0;                                  mb->mvs[0].x = mb->mvs[1].x = mb->mvs[2].x = mb->mvs[3].x = 0;
1054                                  mb->mvs[0].y = mb->mvs[1].y = mb->mvs[2].y = mb->mvs[3].y = 0;                                  mb->mvs[0].y = mb->mvs[1].y = mb->mvs[2].y = mb->mvs[3].y = 0;
   
1055                                  // copy macroblock directly from ref to cur                                  // copy macroblock directly from ref to cur
1056    
1057                                  start_timer();                                  start_timer();
1058    
1059                                  transfer8x8_copy(dec->cur.y + (16 * y) * dec->edged_width +                                  if (reduced_resolution)
1060                                                                   (16 * x),                                  {
1061                                                                   dec->refn[0].y + (16 * y) * dec->edged_width +                                          transfer32x32_copy(dec->cur.y + (32*y)*dec->edged_width + (32*x),
1062                                                                   (16 * x), dec->edged_width);                                                                           dec->refn[0].y + (32*y)*dec->edged_width + (32*x),
1063                                                                             dec->edged_width);
1064                                  transfer8x8_copy(dec->cur.y + (16 * y) * dec->edged_width +  
1065                                                                   (16 * x + 8),                                          transfer16x16_copy(dec->cur.u + (16*y)*dec->edged_width/2 + (16*x),
1066                                                                   dec->refn[0].y + (16 * y) * dec->edged_width +                                                                          dec->refn[0].u + (16*y)*dec->edged_width/2 + (16*x),
                                                                  (16 * x + 8), dec->edged_width);  
   
                                 transfer8x8_copy(dec->cur.y + (16 * y + 8) * dec->edged_width +  
                                                                  (16 * x),  
                                                                  dec->refn[0].y + (16 * y +  
                                                                                                    8) * dec->edged_width +  
                                                                  (16 * x), dec->edged_width);  
   
                                 transfer8x8_copy(dec->cur.y + (16 * y + 8) * dec->edged_width +  
                                                                  (16 * x + 8),  
                                                                  dec->refn[0].y + (16 * y +  
                                                                                                    8) * dec->edged_width +  
                                                                  (16 * x + 8), dec->edged_width);  
   
                                 transfer8x8_copy(dec->cur.u + (8 * y) * dec->edged_width / 2 +  
                                                                  (8 * x),  
                                                                  dec->refn[0].u +  
                                                                  (8 * y) * dec->edged_width / 2 + (8 * x),  
1067                                                                   dec->edged_width / 2);                                                                   dec->edged_width / 2);
1068    
1069                                  transfer8x8_copy(dec->cur.v + (8 * y) * dec->edged_width / 2 +                                          transfer16x16_copy(dec->cur.v + (16*y)*dec->edged_width/2 + (16*x),
1070                                                                   (8 * x),                                                                           dec->refn[0].v + (16*y)*dec->edged_width/2 + (16*x),
                                                                  dec->refn[0].v +  
                                                                  (8 * y) * dec->edged_width / 2 + (8 * x),  
1071                                                                   dec->edged_width / 2);                                                                   dec->edged_width / 2);
1072                                    }
1073                                    else
1074                                    {
1075                                            transfer16x16_copy(dec->cur.y + (16*y)*dec->edged_width + (16*x),
1076                                                                             dec->refn[0].y + (16*y)*dec->edged_width + (16*x),
1077                                                                             dec->edged_width);
1078    
1079                                            transfer8x8_copy(dec->cur.u + (8*y)*dec->edged_width/2 + (8*x),
1080                                                                            dec->refn[0].u + (8*y)*dec->edged_width/2 + (8*x),
1081                                                                            dec->edged_width/2);
1082    
1083                                            transfer8x8_copy(dec->cur.v + (8*y)*dec->edged_width/2 + (8*x),
1084                                                                             dec->refn[0].v + (8*y)*dec->edged_width/2 + (8*x),
1085                                                                             dec->edged_width/2);
1086                                    }
1087    
1088                                  stop_transfer_timer();                                  stop_transfer_timer();
1089    
1090                                  if(dec->out_frm && cp_mb > 0) {                                  if(dec->out_frm && cp_mb > 0) {
1091                                    output_slice(&dec->cur, dec->edged_width,dec->width,dec->out_frm,st_mb,y,cp_mb);                                    output_slice(&dec->cur, dec->edged_width,dec->width,dec->out_frm,st_mb,y,cp_mb);
1092                                    cp_mb = 0;                                    cp_mb = 0;
# Line 882  Line 1178 
1178                  uv_dx = pMB->mvs[0].x;                  uv_dx = pMB->mvs[0].x;
1179                  uv_dy = pMB->mvs[0].y;                  uv_dy = pMB->mvs[0].y;
1180    
1181                  uv_dx = (uv_dx & 3) ? (uv_dx >> 1) | 1 : uv_dx / 2;                  if (dec->quarterpel)
1182                  uv_dy = (uv_dy & 3) ? (uv_dy >> 1) | 1 : uv_dy / 2;                  {
1183                            uv_dx /= 2;
1184                            uv_dy /= 2;
1185                    }
1186    
1187                    uv_dx = (uv_dx >> 1) + roundtab_79[uv_dx & 0x3];
1188                    uv_dy = (uv_dy >> 1) + roundtab_79[uv_dy & 0x3];
1189          } else {          } else {
1190                  int sum;                  int sum;
1191    
1192                    if(dec->quarterpel)
1193                            sum = (pMB->mvs[0].x / 2) + (pMB->mvs[1].x / 2) + (pMB->mvs[2].x / 2) + (pMB->mvs[3].x / 2);
1194                    else
1195                  sum = pMB->mvs[0].x + pMB->mvs[1].x + pMB->mvs[2].x + pMB->mvs[3].x;                  sum = pMB->mvs[0].x + pMB->mvs[1].x + pMB->mvs[2].x + pMB->mvs[3].x;
                 uv_dx =  
                         (sum ==  
                          0 ? 0 : SIGN(sum) * (roundtab[ABS(sum) % 16] +  
                                                                   (ABS(sum) / 16) * 2));  
1196    
1197                    uv_dx = (sum >> 3) + roundtab_76[sum & 0xf];
1198    
1199                    if(dec->quarterpel)
1200                            sum = (pMB->mvs[0].y / 2) + (pMB->mvs[1].y / 2) + (pMB->mvs[2].y / 2) + (pMB->mvs[3].y / 2);
1201                    else
1202                  sum = pMB->mvs[0].y + pMB->mvs[1].y + pMB->mvs[2].y + pMB->mvs[3].y;                  sum = pMB->mvs[0].y + pMB->mvs[1].y + pMB->mvs[2].y + pMB->mvs[3].y;
1203                  uv_dy =  
1204                          (sum ==                  uv_dy = (sum >> 3) + roundtab_76[sum & 0xf];
                          0 ? 0 : SIGN(sum) * (roundtab[ABS(sum) % 16] +  
                                                                   (ABS(sum) / 16) * 2));  
1205          }          }
1206    
1207          start_timer();          start_timer();
1208            if(dec->quarterpel) {
1209                    interpolate16x16_quarterpel(dec->cur.y, dec->refn[ref].y, dec->qtmp.y, dec->qtmp.y + 64,
1210                                                                        dec->qtmp.y + 128, 16*x_pos, 16*y_pos,
1211                                                                        pMB->mvs[0].x, pMB->mvs[0].y, stride, 0);
1212            }
1213            else {
1214          interpolate8x8_switch(dec->cur.y, dec->refn[ref].y, 16 * x_pos, 16 * y_pos,          interpolate8x8_switch(dec->cur.y, dec->refn[ref].y, 16 * x_pos, 16 * y_pos,
1215                                                    pMB->mvs[0].x, pMB->mvs[0].y, stride, 0);                                                    pMB->mvs[0].x, pMB->mvs[0].y, stride, 0);
1216          interpolate8x8_switch(dec->cur.y, dec->refn[ref].y, 16 * x_pos + 8,                  interpolate8x8_switch(dec->cur.y, dec->refn[ref].y, 16*x_pos + 8, 16*y_pos,
1217                                                    16 * y_pos, pMB->mvs[1].x, pMB->mvs[1].y, stride, 0);                                                        pMB->mvs[1].x, pMB->mvs[1].y, stride, 0);
1218          interpolate8x8_switch(dec->cur.y, dec->refn[ref].y, 16 * x_pos,                  interpolate8x8_switch(dec->cur.y, dec->refn[ref].y, 16*x_pos, 16*y_pos + 8,
1219                                                    16 * y_pos + 8, pMB->mvs[2].x, pMB->mvs[2].y, stride,                                                            pMB->mvs[2].x, pMB->mvs[2].y, stride, 0);
1220                                                    0);                  interpolate8x8_switch(dec->cur.y, dec->refn[ref].y, 16*x_pos + 8, 16*y_pos + 8,
1221          interpolate8x8_switch(dec->cur.y, dec->refn[ref].y, 16 * x_pos + 8,                                                            pMB->mvs[3].x, pMB->mvs[3].y, stride, 0);
1222                                                    16 * y_pos + 8, pMB->mvs[3].x, pMB->mvs[3].y, stride,          }
1223                                                    0);  
1224          interpolate8x8_switch(dec->cur.u, dec->refn[ref].u, 8 * x_pos, 8 * y_pos,          interpolate8x8_switch(dec->cur.u, dec->refn[ref].u, 8 * x_pos, 8 * y_pos,
1225                                                    uv_dx, uv_dy, stride2, 0);                                                    uv_dx, uv_dy, stride2, 0);
1226          interpolate8x8_switch(dec->cur.v, dec->refn[ref].v, 8 * x_pos, 8 * y_pos,          interpolate8x8_switch(dec->cur.v, dec->refn[ref].v, 8 * x_pos, 8 * y_pos,
# Line 997  Line 1307 
1307                  uv_dx = pMB->mvs[0].x;                  uv_dx = pMB->mvs[0].x;
1308                  uv_dy = pMB->mvs[0].y;                  uv_dy = pMB->mvs[0].y;
1309    
                 uv_dx = (uv_dx & 3) ? (uv_dx >> 1) | 1 : uv_dx / 2;  
                 uv_dy = (uv_dy & 3) ? (uv_dy >> 1) | 1 : uv_dy / 2;  
   
1310                  b_uv_dx = pMB->b_mvs[0].x;                  b_uv_dx = pMB->b_mvs[0].x;
1311                  b_uv_dy = pMB->b_mvs[0].y;                  b_uv_dy = pMB->b_mvs[0].y;
1312    
1313                  b_uv_dx = (uv_dx & 3) ? (uv_dx >> 1) | 1 : uv_dx / 2;                  if (dec->quarterpel)
1314                  b_uv_dy = (uv_dy & 3) ? (uv_dy >> 1) | 1 : uv_dy / 2;                  {
1315                            uv_dx /= 2;
1316                            uv_dy /= 2;
1317    
1318                            b_uv_dx /= 2;
1319                            b_uv_dy /= 2;
1320                    }
1321    
1322                    uv_dx = (uv_dx >> 1) + roundtab_79[uv_dx & 0x3];
1323                    uv_dy = (uv_dy >> 1) + roundtab_79[uv_dy & 0x3];
1324    
1325                    b_uv_dx = (b_uv_dx >> 1) + roundtab_79[b_uv_dx & 0x3];
1326                    b_uv_dy = (b_uv_dy >> 1) + roundtab_79[b_uv_dy & 0x3];
1327          } else {          } else {
1328                  int sum;                  int sum;
1329    
1330                    if(dec->quarterpel)
1331                            sum = (pMB->mvs[0].x / 2) + (pMB->mvs[1].x / 2) + (pMB->mvs[2].x / 2) + (pMB->mvs[3].x / 2);
1332                    else
1333                  sum = pMB->mvs[0].x + pMB->mvs[1].x + pMB->mvs[2].x + pMB->mvs[3].x;                  sum = pMB->mvs[0].x + pMB->mvs[1].x + pMB->mvs[2].x + pMB->mvs[3].x;
                 uv_dx =  
                         (sum ==  
                          0 ? 0 : SIGN(sum) * (roundtab[ABS(sum) % 16] +  
                                                                   (ABS(sum) / 16) * 2));  
1334    
1335                    uv_dx = (sum >> 3) + roundtab_76[sum & 0xf];
1336    
1337                    if(dec->quarterpel)
1338                            sum = (pMB->mvs[0].y / 2) + (pMB->mvs[1].y / 2) + (pMB->mvs[2].y / 2) + (pMB->mvs[3].y / 2);
1339                    else
1340                  sum = pMB->mvs[0].y + pMB->mvs[1].y + pMB->mvs[2].y + pMB->mvs[3].y;                  sum = pMB->mvs[0].y + pMB->mvs[1].y + pMB->mvs[2].y + pMB->mvs[3].y;
1341                  uv_dy =  
1342                          (sum ==                  uv_dy = (sum >> 3) + roundtab_76[sum & 0xf];
1343                           0 ? 0 : SIGN(sum) * (roundtab[ABS(sum) % 16] +  
1344                                                                    (ABS(sum) / 16) * 2));  
1345                    if(dec->quarterpel)
1346                  sum =                          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);
1347                          pMB->b_mvs[0].x + pMB->b_mvs[1].x + pMB->b_mvs[2].x +                  else
1348                          pMB->b_mvs[3].x;                          sum = pMB->b_mvs[0].x + pMB->b_mvs[1].x + pMB->b_mvs[2].x + pMB->b_mvs[3].x;
1349                  b_uv_dx =  
1350                          (sum ==                  b_uv_dx = (sum >> 3) + roundtab_76[sum & 0xf];
1351                           0 ? 0 : SIGN(sum) * (roundtab[ABS(sum) % 16] +  
1352                                                                    (ABS(sum) / 16) * 2));                  if(dec->quarterpel)
1353                            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);
1354                  sum =                  else
1355                          pMB->b_mvs[0].y + pMB->b_mvs[1].y + pMB->b_mvs[2].y +                          sum = pMB->b_mvs[0].y + pMB->b_mvs[1].y + pMB->b_mvs[2].y + pMB->b_mvs[3].y;
1356                          pMB->b_mvs[3].y;  
1357                  b_uv_dy =                  b_uv_dy = (sum >> 3) + roundtab_76[sum & 0xf];
                         (sum ==  
                          0 ? 0 : SIGN(sum) * (roundtab[ABS(sum) % 16] +  
                                                                   (ABS(sum) / 16) * 2));  
1358          }          }
1359    
1360    
1361          start_timer();          start_timer();
1362            if(dec->quarterpel) {
1363                    if((pMB->mode == MODE_INTER || pMB->mode == MODE_INTER_Q))
1364                            interpolate16x16_quarterpel(dec->cur.y, forward.y, dec->qtmp.y, dec->qtmp.y + 64,
1365                                                                                dec->qtmp.y + 128, 16*x_pos, 16*y_pos,
1366                                                                                pMB->mvs[0].x, pMB->mvs[0].y, stride, 0);
1367                    else {
1368                            interpolate8x8_quarterpel(dec->cur.y, forward.y, dec->qtmp.y, dec->qtmp.y + 64,
1369                                                                                dec->qtmp.y + 128, 16*x_pos, 16*y_pos,
1370                                                                                pMB->mvs[0].x, pMB->mvs[0].y, stride, 0);
1371                            interpolate8x8_quarterpel(dec->cur.y, forward.y, dec->qtmp.y, dec->qtmp.y + 64,
1372                                                                                dec->qtmp.y + 128, 16*x_pos + 8, 16*y_pos,
1373                                                                                pMB->mvs[1].x, pMB->mvs[1].y, stride, 0);
1374                            interpolate8x8_quarterpel(dec->cur.y, forward.y, dec->qtmp.y, dec->qtmp.y + 64,
1375                                                                                dec->qtmp.y + 128, 16*x_pos, 16*y_pos + 8,
1376                                                                                pMB->mvs[2].x, pMB->mvs[2].y, stride, 0);
1377                            interpolate8x8_quarterpel(dec->cur.y, forward.y, dec->qtmp.y, dec->qtmp.y + 64,
1378                                                                                dec->qtmp.y + 128, 16*x_pos + 8, 16*y_pos + 8,
1379                                                                                pMB->mvs[3].x, pMB->mvs[3].y, stride, 0);
1380                    }
1381            }
1382            else {
1383          interpolate8x8_switch(dec->cur.y, forward.y, 16 * x_pos, 16 * y_pos,          interpolate8x8_switch(dec->cur.y, forward.y, 16 * x_pos, 16 * y_pos,
1384                                                    pMB->mvs[0].x, pMB->mvs[0].y, stride, 0);                                                    pMB->mvs[0].x, pMB->mvs[0].y, stride, 0);
1385          interpolate8x8_switch(dec->cur.y, forward.y, 16 * x_pos + 8, 16 * y_pos,          interpolate8x8_switch(dec->cur.y, forward.y, 16 * x_pos + 8, 16 * y_pos,
# Line 1048  Line 1389 
1389          interpolate8x8_switch(dec->cur.y, forward.y, 16 * x_pos + 8,          interpolate8x8_switch(dec->cur.y, forward.y, 16 * x_pos + 8,
1390                                                    16 * y_pos + 8, pMB->mvs[3].x, pMB->mvs[3].y, stride,                                                    16 * y_pos + 8, pMB->mvs[3].x, pMB->mvs[3].y, stride,
1391                                                    0);                                                    0);
1392            }
1393    
1394          interpolate8x8_switch(dec->cur.u, forward.u, 8 * x_pos, 8 * y_pos, uv_dx,          interpolate8x8_switch(dec->cur.u, forward.u, 8 * x_pos, 8 * y_pos, uv_dx,
1395                                                    uv_dy, stride2, 0);                                                    uv_dy, stride2, 0);
1396          interpolate8x8_switch(dec->cur.v, forward.v, 8 * x_pos, 8 * y_pos, uv_dx,          interpolate8x8_switch(dec->cur.v, forward.v, 8 * x_pos, 8 * y_pos, uv_dx,
1397                                                    uv_dy, stride2, 0);                                                    uv_dy, stride2, 0);
1398    
1399    
1400          interpolate8x8_switch(dec->refn[2].y, backward.y, 16 * x_pos, 16 * y_pos,          if(dec->quarterpel) {
1401                    if((pMB->mode == MODE_INTER || pMB->mode == MODE_INTER_Q))
1402                            interpolate16x16_quarterpel(dec->tmp.y, backward.y, dec->qtmp.y, dec->qtmp.y + 64,
1403                                                                                dec->qtmp.y + 128, 16*x_pos, 16*y_pos,
1404                                                    pMB->b_mvs[0].x, pMB->b_mvs[0].y, stride, 0);                                                    pMB->b_mvs[0].x, pMB->b_mvs[0].y, stride, 0);
1405          interpolate8x8_switch(dec->refn[2].y, backward.y, 16 * x_pos + 8,                  else {
1406                            interpolate8x8_quarterpel(dec->tmp.y, backward.y, dec->qtmp.y, dec->qtmp.y + 64,
1407                                                                                dec->qtmp.y + 128, 16*x_pos, 16*y_pos,
1408                                                                                pMB->b_mvs[0].x, pMB->b_mvs[0].y, stride, 0);
1409                            interpolate8x8_quarterpel(dec->tmp.y, backward.y, dec->qtmp.y, dec->qtmp.y + 64,
1410                                                                                dec->qtmp.y + 128, 16*x_pos + 8, 16*y_pos,
1411                                                                                pMB->b_mvs[1].x, pMB->b_mvs[1].y, stride, 0);
1412                            interpolate8x8_quarterpel(dec->tmp.y, backward.y, dec->qtmp.y, dec->qtmp.y + 64,
1413                                                                                dec->qtmp.y + 128, 16*x_pos, 16*y_pos + 8,
1414                                                                                pMB->b_mvs[2].x, pMB->b_mvs[2].y, stride, 0);
1415                            interpolate8x8_quarterpel(dec->tmp.y, backward.y, dec->qtmp.y, dec->qtmp.y + 64,
1416                                                                                dec->qtmp.y + 128, 16*x_pos + 8, 16*y_pos + 8,
1417                                                                                pMB->b_mvs[3].x, pMB->b_mvs[3].y, stride, 0);
1418                    }
1419            }
1420            else {
1421                    interpolate8x8_switch(dec->tmp.y, backward.y, 16 * x_pos, 16 * y_pos,
1422                                                              pMB->b_mvs[0].x, pMB->b_mvs[0].y, stride, 0);
1423                    interpolate8x8_switch(dec->tmp.y, backward.y, 16 * x_pos + 8,
1424                                                    16 * y_pos, pMB->b_mvs[1].x, pMB->b_mvs[1].y, stride,                                                    16 * y_pos, pMB->b_mvs[1].x, pMB->b_mvs[1].y, stride,
1425                                                    0);                                                    0);
1426          interpolate8x8_switch(dec->refn[2].y, backward.y, 16 * x_pos,                  interpolate8x8_switch(dec->tmp.y, backward.y, 16 * x_pos,
1427                                                    16 * y_pos + 8, pMB->b_mvs[2].x, pMB->b_mvs[2].y,                                                    16 * y_pos + 8, pMB->b_mvs[2].x, pMB->b_mvs[2].y,
1428                                                    stride, 0);                                                    stride, 0);
1429          interpolate8x8_switch(dec->refn[2].y, backward.y, 16 * x_pos + 8,                  interpolate8x8_switch(dec->tmp.y, backward.y, 16 * x_pos + 8,
1430                                                    16 * y_pos + 8, pMB->b_mvs[3].x, pMB->b_mvs[3].y,                                                    16 * y_pos + 8, pMB->b_mvs[3].x, pMB->b_mvs[3].y,
1431                                                    stride, 0);                                                    stride, 0);
1432          interpolate8x8_switch(dec->refn[2].u, backward.u, 8 * x_pos, 8 * y_pos,          }
1433    
1434            interpolate8x8_switch(dec->tmp.u, backward.u, 8 * x_pos, 8 * y_pos,
1435                                                    b_uv_dx, b_uv_dy, stride2, 0);                                                    b_uv_dx, b_uv_dy, stride2, 0);
1436          interpolate8x8_switch(dec->refn[2].v, backward.v, 8 * x_pos, 8 * y_pos,          interpolate8x8_switch(dec->tmp.v, backward.v, 8 * x_pos, 8 * y_pos,
1437                                                    b_uv_dx, b_uv_dy, stride2, 0);                                                    b_uv_dx, b_uv_dy, stride2, 0);
1438    
1439          interpolate8x8_avg2(dec->cur.y + (16 * y_pos * stride) + 16 * x_pos,          interpolate8x8_avg2(dec->cur.y + (16 * y_pos * stride) + 16 * x_pos,
1440                                                  dec->cur.y + (16 * y_pos * stride) + 16 * x_pos,                                                  dec->cur.y + (16 * y_pos * stride) + 16 * x_pos,
1441                                                  dec->refn[2].y + (16 * y_pos * stride) + 16 * x_pos,                                                  dec->tmp.y + (16 * y_pos * stride) + 16 * x_pos,
1442                                                  stride, 0);                                                  stride, 1, 8);
1443    
1444          interpolate8x8_avg2(dec->cur.y + (16 * y_pos * stride) + 16 * x_pos + 8,          interpolate8x8_avg2(dec->cur.y + (16 * y_pos * stride) + 16 * x_pos + 8,
1445                                                  dec->cur.y + (16 * y_pos * stride) + 16 * x_pos + 8,                                                  dec->cur.y + (16 * y_pos * stride) + 16 * x_pos + 8,
1446                                                  dec->refn[2].y + (16 * y_pos * stride) + 16 * x_pos + 8,                                                  dec->tmp.y + (16 * y_pos * stride) + 16 * x_pos + 8,
1447                                                  stride, 0);                                                  stride, 1, 8);
1448    
1449          interpolate8x8_avg2(dec->cur.y + ((16 * y_pos + 8) * stride) + 16 * x_pos,          interpolate8x8_avg2(dec->cur.y + ((16 * y_pos + 8) * stride) + 16 * x_pos,
1450                                                  dec->cur.y + ((16 * y_pos + 8) * stride) + 16 * x_pos,                                                  dec->cur.y + ((16 * y_pos + 8) * stride) + 16 * x_pos,
1451                                                  dec->refn[2].y + ((16 * y_pos + 8) * stride) + 16 * x_pos,                                                  dec->tmp.y + ((16 * y_pos + 8) * stride) + 16 * x_pos,
1452                                                  stride, 0);                                                  stride, 1, 8);
1453    
1454          interpolate8x8_avg2(dec->cur.y + ((16 * y_pos + 8) * stride) + 16 * x_pos + 8,          interpolate8x8_avg2(dec->cur.y + ((16 * y_pos + 8) * stride) + 16 * x_pos + 8,
1455                                                  dec->cur.y + ((16 * y_pos + 8) * stride) + 16 * x_pos + 8,                                                  dec->cur.y + ((16 * y_pos + 8) * stride) + 16 * x_pos + 8,
1456                                                  dec->refn[2].y + ((16 * y_pos + 8) * stride) + 16 * x_pos + 8,                                                  dec->tmp.y + ((16 * y_pos + 8) * stride) + 16 * x_pos + 8,
1457                                                  stride, 0);                                                  stride, 1, 8);
1458    
1459          interpolate8x8_avg2(dec->cur.u + (8 * y_pos * stride2) + 8 * x_pos,          interpolate8x8_avg2(dec->cur.u + (8 * y_pos * stride2) + 8 * x_pos,
1460                                                  dec->cur.u + (8 * y_pos * stride2) + 8 * x_pos,                                                  dec->cur.u + (8 * y_pos * stride2) + 8 * x_pos,
1461                                                  dec->refn[2].u + (8 * y_pos * stride2) + 8 * x_pos,                                                  dec->tmp.u + (8 * y_pos * stride2) + 8 * x_pos,
1462                                                  stride2, 0);                                                  stride2, 1, 8);
1463    
1464          interpolate8x8_avg2(dec->cur.v + (8 * y_pos * stride2) + 8 * x_pos,          interpolate8x8_avg2(dec->cur.v + (8 * y_pos * stride2) + 8 * x_pos,
1465                                                  dec->cur.v + (8 * y_pos * stride2) + 8 * x_pos,                                                  dec->cur.v + (8 * y_pos * stride2) + 8 * x_pos,
1466                                                  dec->refn[2].v + (8 * y_pos * stride2) + 8 * x_pos,                                                  dec->tmp.v + (8 * y_pos * stride2) + 8 * x_pos,
1467                                                  stride2, 0);                                                  stride2, 1, 8);
1468    
1469          stop_comp_timer();          stop_comp_timer();
1470    
# Line 1227  Line 1593 
1593                          mb->b_mvs[0] = mb->b_mvs[1] = mb->b_mvs[2] = mb->b_mvs[3] =                          mb->b_mvs[0] = mb->b_mvs[1] = mb->b_mvs[2] = mb->b_mvs[3] =
1594                          mb->mvs[0] = mb->mvs[1] = mb->mvs[2] = mb->mvs[3] = zeromv;                          mb->mvs[0] = mb->mvs[1] = mb->mvs[2] = mb->mvs[3] = zeromv;
1595    
1596                          // the last P_VOP is skip macroblock ?                          // skip if the co-located P_VOP macroblock is not coded
1597                            // if not codec in co-located S_VOP macroblock is _not_ automatically skipped
1598    
1599                          if (last_mb->mode == MODE_NOT_CODED) {                          if (last_mb->mode == MODE_NOT_CODED) {
1600                                  //DEBUG2("Skip MB in B-frame at (X,Y)=!",x,y);                                  //DEBUG2("Skip MB in B-frame at (X,Y)=!",x,y);
1601                                  mb->cbp = 0;                                  mb->cbp = 0;
# Line 1341  Line 1709 
1709                                  break;                                  break;
1710    
1711                          default:                          default:
1712                                  DEBUG1("Not support B-frame mb_type =", mb->mb_type);                                  DPRINTF(DPRINTF_ERROR,"Not support B-frame mb_type = %i", mb->mb_type);
1713                          }                          }
1714    
1715                  }                                               // end of FOR                  }                                               // end of FOR
# Line 1366  Line 1734 
1734          *mb2 = temp;          *mb2 = temp;
1735  }  }
1736    
1737    
1738    /* perform post processing if necessary, and output the image */
1739    void decoder_output(DECODER * dec, IMAGE * img, MACROBLOCK * mbs,
1740                                            const XVID_DEC_FRAME * frame, int pp_disable)
1741    {
1742    
1743            if ((frame->general & (XVID_DEC_DEBLOCKY|XVID_DEC_DEBLOCKUV)) && !pp_disable)   /* post process */
1744            {
1745                    /* note: image is stored to tmp */
1746                    image_copy(&dec->tmp, img, dec->edged_width, dec->height);
1747                    image_deblock_rrv(&dec->tmp, dec->edged_width,
1748                                                    mbs, dec->mb_width, dec->mb_height, dec->mb_width,
1749                                                    8, frame->general);
1750                    img = &dec->tmp;
1751            }
1752    
1753            image_output(img, dec->width, dec->height,
1754                                     dec->edged_width, frame->image, frame->stride,
1755                                     frame->colorspace, dec->interlacing);
1756    }
1757    
1758    
1759  int  int
1760  decoder_decode(DECODER * dec,  decoder_decode(DECODER * dec,
1761                             XVID_DEC_FRAME * frame)                             XVID_DEC_FRAME * frame, XVID_DEC_STATS * stats)
1762  {  {
1763    
1764          Bitstream bs;          Bitstream bs;
1765          uint32_t rounding;          uint32_t rounding;
1766            uint32_t reduced_resolution;
1767          uint32_t quant;          uint32_t quant;
1768          uint32_t fcode_forward;          uint32_t fcode_forward;
1769          uint32_t fcode_backward;          uint32_t fcode_backward;
1770          uint32_t intra_dc_threshold;          uint32_t intra_dc_threshold;
1771            WARPPOINTS gmc_warp;
1772          uint32_t vop_type;          uint32_t vop_type;
1773            int success = 0;
1774            int output = 0;
1775            int seen_something = 0;
1776    
1777          start_global_timer();          start_global_timer();
1778    
1779            dec->low_delay_default = (frame->general & XVID_DEC_LOWDELAY);
1780          dec->out_frm = (frame->colorspace == XVID_CSP_EXTERN) ? frame->image : NULL;          dec->out_frm = (frame->colorspace == XVID_CSP_EXTERN) ? frame->image : NULL;
1781    
1782            if ((frame->general & XVID_DEC_DISCONTINUITY))
1783                    dec->frames = 0;
1784    
1785            if (frame->length < 0)  /* decoder flush */
1786            {
1787                    /* if  not decoding "low_delay/packed", and this isn't low_delay and
1788                        we have a reference frame, then outout the reference frame */
1789                    if (!(dec->low_delay_default && dec->packed_mode) && !dec->low_delay && dec->frames>0)
1790                    {
1791                            decoder_output(dec, &dec->refn[0], dec->mbs, frame, dec->last_reduced_resolution);
1792                            output = 1;
1793                    }
1794    
1795                    frame->length = 0;
1796                    if (stats)
1797                    {
1798                            stats->notify = output ? XVID_DEC_VOP : XVID_DEC_NOTHING;
1799                            stats->data.vop.time_base = (int)dec->time_base;
1800                            stats->data.vop.time_increment = 0;     //XXX: todo
1801                    }
1802    
1803                    emms();
1804    
1805                    stop_global_timer();
1806                    return XVID_ERR_OK;
1807            }
1808    
1809          BitstreamInit(&bs, frame->bitstream, frame->length);          BitstreamInit(&bs, frame->bitstream, frame->length);
1810    
1811          if(BitstreamShowBits(&bs, 8) == 0x7f)          // XXX: 0x7f is only valid whilst decoding vfw xvid/divx5 avi's
1812            if(dec->low_delay_default && frame->length == 1 && BitstreamShowBits(&bs, 8) == 0x7f)
1813            {
1814                    if (stats)
1815                            stats->notify = XVID_DEC_VOP;
1816                    frame->length = 1;
1817                    image_output(&dec->refn[0], dec->width, dec->height, dec->edged_width,
1818                                             frame->image, frame->stride, frame->colorspace, dec->interlacing);
1819                    emms();
1820                  return XVID_ERR_OK;                  return XVID_ERR_OK;
1821            }
1822    
1823          // add by chenm001 <chenm001@163.com>  repeat:
1824          // for support B-frame to reference last 2 frame  
1825          dec->frames++;          vop_type =      BitstreamReadHeaders(&bs, dec, &rounding, &reduced_resolution,
1826          vop_type =                          &quant, &fcode_forward, &fcode_backward, &intra_dc_threshold, &gmc_warp);
1827                  BitstreamReadHeaders(&bs, dec, &rounding, &quant, &fcode_forward,  
1828                                                           &fcode_backward, &intra_dc_threshold);          DPRINTF(DPRINTF_HEADER, "vop_type=%i,  packed=%i,  time=%lli,  time_pp=%i,  time_bp=%i",
1829                                                            vop_type,       dec->packed_mode, dec->time, dec->time_pp, dec->time_bp);
1830    
1831            if (vop_type == - 1)
1832            {
1833                    if (success) goto done;
1834                    emms();
1835                    return XVID_ERR_FAIL;
1836            }
1837    
1838            if (vop_type == -2 || vop_type == -3)
1839            {
1840                    if (vop_type == -3)
1841                            decoder_resize(dec);
1842    
1843                    if (stats)
1844                    {
1845                            stats->notify = XVID_DEC_VOL;
1846                            stats->data.vol.general = 0;
1847                            if (dec->interlacing)
1848                                    stats->data.vol.general |= XVID_INTERLACING;
1849                            stats->data.vol.width = dec->width;
1850                            stats->data.vol.height = dec->height;
1851                            stats->data.vol.aspect_ratio = dec->aspect_ratio;
1852                            stats->data.vol.par_width = dec->par_width;
1853                            stats->data.vol.par_height = dec->par_height;
1854                            frame->length = BitstreamPos(&bs) / 8;
1855                            emms();
1856                            return XVID_ERR_OK;
1857                    }
1858                    goto repeat;
1859            }
1860    
1861          dec->p_bmv.x = dec->p_bmv.y = dec->p_fmv.y = dec->p_fmv.y = 0;  // init pred vector to 0          dec->p_bmv.x = dec->p_bmv.y = dec->p_fmv.y = dec->p_fmv.y = 0;  // init pred vector to 0
1862    
         switch (vop_type) {  
         case P_VOP:  
                 decoder_pframe(dec, &bs, rounding, quant, fcode_forward,  
                                            intra_dc_threshold);  
 #ifdef BFRAMES_DEC  
                 DEBUG1("P_VOP  Time=", dec->time);  
 #endif  
                 break;  
1863    
1864            /* packed_mode: special-N_VOP treament */
1865            if (dec->packed_mode && vop_type == N_VOP)
1866            {
1867                    if (dec->low_delay_default && dec->frames > 0)
1868                    {
1869                            decoder_output(dec, &dec->refn[0], dec->last_mbs, frame, dec->last_reduced_resolution);
1870                            output = 1;
1871                    }
1872                    /* ignore otherwise */
1873            }
1874            else if (vop_type != B_VOP)
1875            {
1876                    switch(vop_type)
1877                    {
1878          case I_VOP:          case I_VOP:
1879                  decoder_iframe(dec, &bs, quant, intra_dc_threshold);                          decoder_iframe(dec, &bs, reduced_resolution, quant, intra_dc_threshold);
 #ifdef BFRAMES_DEC  
                 DEBUG1("I_VOP  Time=", dec->time);  
 #endif  
1880                  break;                  break;
1881                    case P_VOP :
1882          case B_VOP:                          decoder_pframe(dec, &bs, rounding, reduced_resolution, quant,
1883  #ifdef BFRAMES_DEC                                                  fcode_forward, intra_dc_threshold, NULL);
                 if (dec->time_pp > dec->time_bp) {  
                         DEBUG1("B_VOP  Time=", dec->time);  
                         decoder_bframe(dec, &bs, quant, fcode_forward, fcode_backward);  
                 } else {  
                         DEBUG("broken B-frame!");  
                 }  
 #else  
                 image_copy(&dec->cur, &dec->refn[0], dec->edged_width, dec->height);  
 #endif  
1884                  break;                  break;
1885                    case S_VOP :
1886          case N_VOP:                             // vop not coded                          decoder_pframe(dec, &bs, rounding, reduced_resolution, quant,
1887                  // when low_delay==0, N_VOP's should interpolate between the past and future frames                                                  fcode_forward, intra_dc_threshold, &gmc_warp);
1888                            break;
1889                    case N_VOP :
1890                  image_copy(&dec->cur, &dec->refn[0], dec->edged_width, dec->height);                  image_copy(&dec->cur, &dec->refn[0], dec->edged_width, dec->height);
 #ifdef BFRAMES_DEC  
                 DEBUG1("N_VOP  Time=", dec->time);  
 #endif  
1891                  break;                  break;
1892                    }
1893    
1894          default:                  if (reduced_resolution)
1895                  return XVID_ERR_FAIL;                  {
1896                            image_deblock_rrv(&dec->cur, dec->edged_width, dec->mbs,
1897                                    (dec->width + 31) / 32, (dec->height + 31) / 32, dec->mb_width,
1898                                    16, XVID_DEC_DEBLOCKY|XVID_DEC_DEBLOCKUV);
1899          }          }
1900    
1901  #ifdef BFRAMES_DEC_DEBUG                  /* note: for packed_mode, output is performed when the special-N_VOP is decoded */
1902          if (frame->length != BitstreamPos(&bs) / 8){                  if (!(dec->low_delay_default && dec->packed_mode))
1903                  DEBUG2("InLen/UseLen",frame->length, BitstreamPos(&bs) / 8);                  {
1904                            if (dec->low_delay)
1905                            {
1906                                    decoder_output(dec, &dec->cur, dec->mbs, frame, reduced_resolution);
1907                                    output = 1;
1908                            }
1909                            else if (dec->frames > 0)       /* is the reference frame valid? */
1910                            {
1911                                    /* output the reference frame */
1912                                    decoder_output(dec, &dec->refn[0], dec->last_mbs, frame, dec->last_reduced_resolution);
1913                                    output = 1;
1914                            }
1915          }          }
 #endif  
         frame->length = BitstreamPos(&bs) / 8;  
1916    
1917                    image_swap(&dec->refn[0], &dec->refn[1]);
1918                    image_swap(&dec->cur, &dec->refn[0]);
1919                    mb_swap(&dec->mbs, &dec->last_mbs);
1920                    dec->last_reduced_resolution = reduced_resolution;
1921    
1922  #ifdef BFRAMES_DEC                  dec->frames++;
1923          // test if no B_VOP                  seen_something = 1;
1924          if (dec->low_delay || dec->frames == 0) {  
1925  #endif          }else{  /* B_VOP */
1926          image_output(&dec->cur, dec->width, dec->height, dec->edged_width,  
1927                                           frame->image, frame->stride, frame->colorspace);                  if (dec->low_delay)
1928                    {
1929                            DPRINTF(DPRINTF_ERROR, "warning: bvop found in low_delay==1 stream");
1930                            dec->low_delay = 1;
1931                    }
1932    
1933  #ifdef BFRAMES_DEC                  if (dec->frames < 2)
1934                    {
1935                            /* attemping to decode a bvop without atleast 2 reference frames */
1936                            image_printf(&dec->cur, dec->edged_width, dec->height, 16, 16,
1937                                                    "broken b-frame, mising ref frames");
1938                    }else if (dec->time_pp <= dec->time_bp) {
1939                            /* this occurs when dx50_bvop_compatibility==0 sequences are
1940                            decoded in vfw. */
1941                            image_printf(&dec->cur, dec->edged_width, dec->height, 16, 16,
1942                                                    "broken b-frame, tpp=%i tbp=%i", dec->time_pp, dec->time_bp);
1943          } else {          } else {
1944                  if (dec->frames >= 1) {                          decoder_bframe(dec, &bs, quant, fcode_forward, fcode_backward);
                         start_timer();  
                         if ((vop_type == I_VOP || vop_type == P_VOP)) {  
                                 image_output(&dec->refn[0], dec->width, dec->height,  
                                                          dec->edged_width, frame->image, frame->stride,  
                                                          frame->colorspace);  
                         } else if (vop_type == B_VOP) {  
                                 image_output(&dec->cur, dec->width, dec->height,  
                                                          dec->edged_width, frame->image, frame->stride,  
                                                          frame->colorspace);  
1945                          }                          }
1946                          stop_conv_timer();  
1947                    decoder_output(dec, &dec->cur, dec->mbs, frame, reduced_resolution);
1948                    output = 1;
1949                    dec->frames++;
1950                  }                  }
1951    
1952            BitstreamByteAlign(&bs);
1953    
1954            /* low_delay_default mode: repeat in packed_mode */
1955            if (dec->low_delay_default && dec->packed_mode && output == 0 && success == 0)
1956            {
1957                    success = 1;
1958                    goto repeat;
1959          }          }
 #endif  
1960    
1961          if (vop_type == I_VOP || vop_type == P_VOP) {  done :
                 image_swap(&dec->refn[0], &dec->refn[1]);  
                 image_swap(&dec->cur, &dec->refn[0]);  
1962    
1963                  // swap MACROBLOCK          /* low_delay_default mode: if we've gotten here without outputting anything,
1964                  // the Divx will not set the low_delay flage some times             then output the recently decoded frame, or print an error message  */
1965                  // so follow code will wrong to not swap at that time          if (dec->low_delay_default && output == 0)
1966                  // this will broken bitstream! so I'm change it,          {
1967                  // But that is not the best way! can anyone tell me how                  if (dec->packed_mode && seen_something)
1968                  // to do another way?                  {
1969                  // 18-07-2002   MinChen<chenm001@163.com>                          /* output the recently decoded frame */
1970                  //if (!dec->low_delay && vop_type == P_VOP)                          decoder_output(dec, &dec->refn[0], dec->last_mbs, frame, dec->last_reduced_resolution);
1971                  if (vop_type == P_VOP)                          output = 1;
1972                          mb_swap(&dec->mbs, &dec->last_mbs);                  }
1973                    else
1974                    {
1975                            image_clear(&dec->cur, dec->width, dec->height, dec->edged_width, 0, 128, 128);
1976                            image_printf(&dec->cur, dec->edged_width, dec->height, 16, 16,
1977                                    "warning: nothing to output");
1978                            image_printf(&dec->cur, dec->edged_width, dec->height, 16, 64,
1979                                    "bframe decoder lag");
1980    
1981                            decoder_output(dec, &dec->cur, NULL, frame, 1 /*disable pp*/);
1982                    }
1983            }
1984    
1985            frame->length = BitstreamPos(&bs) / 8;
1986    
1987            if (stats)
1988            {
1989                    stats->notify = output ? XVID_DEC_VOP : XVID_DEC_NOTHING;
1990                    stats->data.vop.time_base = (int)dec->time_base;
1991                    stats->data.vop.time_increment = 0;     //XXX: todo
1992          }          }
1993    
1994          emms();          emms();

Legend:
Removed from v.1.37.2.5  
changed lines
  Added in v.1.37.2.28

No admin address has been configured
ViewVC Help
Powered by ViewVC 1.0.4