[svn] / branches / dev-api-4 / xvidcore / src / utils / mbtransquant.c Repository:
ViewVC logotype

Diff of /branches/dev-api-4/xvidcore/src/utils/mbtransquant.c

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

trunk/xvidcore/src/utils/mbtransquant.c revision 3, Fri Mar 8 02:46:11 2002 UTC branches/dev-api-4/xvidcore/src/utils/mbtransquant.c revision 965, Sun Mar 30 13:16:42 2003 UTC
# Line 1  Line 1 
1   /******************************************************************************  /*****************************************************************************
2    *                                                                            *   *
3    *  This file is part of XviD, a free MPEG-4 video encoder/decoder            *   *  XVID MPEG-4 VIDEO CODEC
4    *                                                                            *   *  - MB Transfert/Quantization functions -
5    *  XviD is an implementation of a part of one or more MPEG-4 Video tools     *   *
6    *  as specified in ISO/IEC 14496-2 standard.  Those intending to use this    *   *  Copyright(C) 2001-2003  Peter Ross <pross@xvid.org>
7    *  software module in hardware or software products are advised that its     *   *               2001-2003  Michael Militzer <isibaar@xvid.org>
8    *  use may infringe existing patents or copyrights, and any such use         *   *               2003       Edouard Gomez <ed.gomez@free.fr>
9    *  would be at such party's own risk.  The original developer of this        *   *
10    *  software module and his/her company, and subsequent editors and their     *   *  This program is free software ; you can redistribute it and/or modify
11    *  companies, will have no liability for use of this software or             *   *  it under the terms of the GNU General Public License as published by
12    *  modifications or derivatives thereof.                                     *   *  the Free Software Foundation ; either version 2 of the License, or
13    *                                                                            *   *  (at your option) any later version.
14    *  XviD is free software; you can redistribute it and/or modify it           *   *
15    *  under the terms of the GNU General Public License as published by         *   *  This program is distributed in the hope that it will be useful,
16    *  the Free Software Foundation; either version 2 of the License, or         *   *  but WITHOUT ANY WARRANTY ; without even the implied warranty of
17    *  (at your option) any later version.                                       *   *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
18    *                                                                            *   *  GNU General Public License for more details.
19    *  XviD is distributed in the hope that it will be useful, but               *   *
20    *  WITHOUT ANY WARRANTY; without even the implied warranty of                *   *  You should have received a copy of the GNU General Public License
21    *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the             *   *  along with this program ; if not, write to the Free Software
22    *  GNU General Public License for more details.                              *   *  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
23    *                                                                            *   *
24    *  You should have received a copy of the GNU General Public License         *   * $Id: mbtransquant.c,v 1.21.2.6 2003-03-30 13:16:41 edgomez Exp $
25    *  along with this program; if not, write to the Free Software               *   *
26    *  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA  *   ****************************************************************************/
27    *                                                                            *  
28    ******************************************************************************/  #include <string.h>
   
  /******************************************************************************  
   *                                                                            *  
   *  mbtransquant.c                                                            *  
   *                                                                            *  
   *  Copyright (C) 2001 - Peter Ross <pross@cs.rmit.edu.au>                    *  
   *  Copyright (C) 2001 - Michael Militzer <isibaar@xvid.org>                  *  
   *                                                                            *  
   *  For more information visit the XviD homepage: http://www.xvid.org         *  
   *                                                                            *  
   ******************************************************************************/  
   
  /******************************************************************************  
   *                                                                            *  
   *  Revision history:                                                         *  
   *                                                                            *  
   *  22.12.2001 get_dc_scaler() moved to common.h  
   *  19.11.2001 introduced coefficient thresholding (Isibaar)                  *  
   *  17.11.2001 initial version                                                *  
   *                                                                            *  
   ******************************************************************************/  
29    
30  #include "../portab.h"  #include "../portab.h"
31  #include "mbfunctions.h"  #include "mbfunctions.h"
# Line 60  Line 39 
39  #include "../quant/quant_h263.h"  #include "../quant/quant_h263.h"
40  #include "../encoder.h"  #include "../encoder.h"
41    
42  #define MIN(X, Y) ((X)<(Y)?(X):(Y))  #include "../image/reduced.h"
 #define MAX(X, Y) ((X)>(Y)?(X):(Y))  
   
 #define TOOSMALL_LIMIT 1 /* skip blocks having a coefficient sum below this value */  
43    
44  /* this isnt pretty, but its better than 20 ifdefs */  MBFIELDTEST_PTR MBFieldTest;
45    
46  void MBTransQuantIntra(const MBParam *pParam,  /*
47                         const uint32_t x_pos,   * Skip blocks having a coefficient sum below this value. This value will be
48                         const uint32_t y_pos,   * corrected according to the MB quantizer to avoid artifacts for quant==1
49                         int16_t data[][64],   */
50                             int16_t qcoeff[][64],  #define PVOP_TOOSMALL_LIMIT 1
51                             IMAGE * const pCurrent)  #define BVOP_TOOSMALL_LIMIT 3
52    
53    /*****************************************************************************
54     * Local functions
55     ****************************************************************************/
56    
57    /* permute block and return field dct choice */
58    static __inline uint32_t
59    MBDecideFieldDCT(int16_t data[6 * 64])
60  {  {
61          const uint32_t stride = pParam->edged_width;          uint32_t field = MBFieldTest(data);
         uint32_t i;  
         uint32_t iQuant = pParam->quant;  
         uint8_t *pY_Cur, *pU_Cur, *pV_Cur;  
62    
63      pY_Cur = pCurrent->y + (y_pos << 4) * stride + (x_pos << 4);          if (field)
64      pU_Cur = pCurrent->u + (y_pos << 3) * (stride >> 1) + (x_pos << 3);                  MBFrameToField(data);
     pV_Cur = pCurrent->v + (y_pos << 3) * (stride >> 1) + (x_pos << 3);  
65    
66          for(i = 0; i < 6; i++) {          return field;
67                  uint32_t iDcScaler = get_dc_scaler(iQuant, i < 4);  }
68    
69    /* Performs Forward DCT on all blocks */
70    static __inline void
71    MBfDCT(const MBParam * pParam,
72               FRAMEINFO * frame,
73               MACROBLOCK * pMB,
74               uint32_t x_pos,
75               uint32_t y_pos,
76               int16_t data[6 * 64])
77    {
78            /* Handles interlacing */
79                  start_timer();                  start_timer();
80            pMB->field_dct = 0;
81                  switch(i) {          if ((frame->vol_flags & XVID_VOL_INTERLACING) &&
82                  case 0 :                  (x_pos>0) && (x_pos<pParam->mb_width-1) &&
83                          transfer_8to16copy(data[0], pY_Cur, stride);                  (y_pos>0) && (y_pos<pParam->mb_height-1)) {
84                          break;                  pMB->field_dct = MBDecideFieldDCT(data);
85                  case 1 :          }
86                          transfer_8to16copy(data[1], pY_Cur + 8, stride);          stop_interlacing_timer();
87                          break;  
88                  case 2 :          /* Perform DCT */
89                      transfer_8to16copy(data[2], pY_Cur + 8 * stride, stride);          start_timer();
90                          break;          fdct(&data[0 * 64]);
91                  case 3 :          fdct(&data[1 * 64]);
92                          transfer_8to16copy(data[3], pY_Cur + 8 * stride + 8, stride);          fdct(&data[2 * 64]);
93                          break;          fdct(&data[3 * 64]);
94                  case 4 :          fdct(&data[4 * 64]);
95                          transfer_8to16copy(data[4], pU_Cur, stride / 2);          fdct(&data[5 * 64]);
96                          break;          stop_dct_timer();
                 case 5 :  
                         transfer_8to16copy(data[5], pV_Cur, stride / 2);  
                         break;  
97                  }                  }
                 stop_transfer_timer();  
98    
99    /* Performs Inverse DCT on all blocks */
100    static __inline void
101    MBiDCT(int16_t data[6 * 64],
102               const uint8_t cbp)
103    {
104                  start_timer();                  start_timer();
105                  fdct(data[i]);          if(cbp & (1 << (5 - 0))) idct(&data[0 * 64]);
106                  stop_dct_timer();          if(cbp & (1 << (5 - 1))) idct(&data[1 * 64]);
107            if(cbp & (1 << (5 - 2))) idct(&data[2 * 64]);
108            if(cbp & (1 << (5 - 3))) idct(&data[3 * 64]);
109            if(cbp & (1 << (5 - 4))) idct(&data[4 * 64]);
110            if(cbp & (1 << (5 - 5))) idct(&data[5 * 64]);
111            stop_idct_timer();
112    }
113    
114                  if (pParam->quant_type == H263_QUANT)  /* Quantize all blocks -- Intra mode */
115    static __inline void
116    MBQuantIntra(const MBParam * pParam,
117                             const MACROBLOCK * pMB,
118                         int16_t qcoeff[6 * 64],
119                             int16_t data[6*64])
120                  {                  {
121            int i;
122    
123            for (i = 0; i < 6; i++) {
124                    uint32_t iDcScaler = get_dc_scaler(pMB->quant, i < 4);
125    
126                    /* Quantize the block */
127                          start_timer();                          start_timer();
128                          quant_intra(qcoeff[i], data[i], iQuant, iDcScaler);                  if (!(pParam->vol_flags & XVID_VOL_MPEGQUANT))
129                            quant_intra(&data[i * 64], &qcoeff[i * 64], pMB->quant, iDcScaler);
130                    else
131                            quant4_intra(&data[i * 64], &qcoeff[i * 64], pMB->quant, iDcScaler);
132                          stop_quant_timer();                          stop_quant_timer();
133            }
134    }
135    
136    /* DeQuantize all blocks -- Intra mode */
137    static __inline void
138    MBDeQuantIntra(const MBParam * pParam,
139                               const int iQuant,
140                               int16_t qcoeff[6 * 64],
141                               int16_t data[6*64])
142    {
143            int i;
144    
145            for (i = 0; i < 6; i++) {
146                    uint32_t iDcScaler = get_dc_scaler(iQuant, i < 4);
147    
148                          start_timer();                          start_timer();
149                          dequant_intra(data[i], qcoeff[i], iQuant, iDcScaler);                  if (!(pParam->vol_flags & XVID_VOL_MPEGQUANT))
150                            dequant_intra(&qcoeff[i * 64], &data[i * 64], iQuant, iDcScaler);
151                    else
152                            dequant4_intra(&qcoeff[i * 64], &data[i * 64], iQuant, iDcScaler);
153                          stop_iquant_timer();                          stop_iquant_timer();
154                  }                  }
155                  else  }
156    
157    /* Quantize all blocks -- Inter mode */
158    static __inline uint8_t
159    MBQuantInter(const MBParam * pParam,
160                             const MACROBLOCK * pMB,
161                             int16_t data[6 * 64],
162                             int16_t qcoeff[6 * 64],
163                             int bvop,
164                             int limit)
165                  {                  {
166    
167            int i;
168            uint8_t cbp = 0;
169            int sum;
170            int code_block;
171    
172            for (i = 0; i < 6; i++) {
173    
174                    /* Quantize the block */
175                          start_timer();                          start_timer();
176                          quant4_intra(qcoeff[i], data[i], iQuant, iDcScaler);                  if (!(pParam->vol_flags & XVID_VOL_MPEGQUANT))
177                            sum = quant_inter(&qcoeff[i * 64], &data[i * 64], pMB->quant);
178                    else
179                            sum = quant4_inter(&qcoeff[i * 64], &data[i * 64], pMB->quant);
180                          stop_quant_timer();                          stop_quant_timer();
181    
182                          start_timer();                  /*
183                          dequant4_intra(data[i], qcoeff[i], iQuant, iDcScaler);                   * We code the block if the sum is higher than the limit and if the first
184                          stop_iquant_timer();                   * two AC coefficients in zig zag order are not zero.
185                     */
186                    code_block = 0;
187                    if ((sum >= limit) || (qcoeff[i*64+1] != 0) || (qcoeff[i*64+8] != 0)) {
188                            code_block = 1;
189                    } else {
190    
191                            if (bvop && (pMB->mode == MODE_DIRECT || pMB->mode == MODE_DIRECT_NO4V)) {
192                                    /* dark blocks prevention for direct mode */
193                                    if ((qcoeff[i*64] < -1) || (qcoeff[i*64] > 0))
194                                            code_block = 1;
195                            } else {
196                                    /* not direct mode */
197                                    if (qcoeff[i*64] != 0)
198                                            code_block = 1;
199                            }
200                  }                  }
201    
202                  start_timer();                  /* Set the corresponding cbp bit */
203                  idct(data[i]);                  cbp |= code_block << (5 - i);
                 stop_idct_timer();  
204    
205                  start_timer();          }
206    
207                  switch(i) {          return(cbp);
208                  case 0:  }
209                          transfer_16to8copy(pY_Cur, data[0], stride);  
210                          break;  /* DeQuantize all blocks -- Inter mode */
211                  case 1:  static __inline void
212                          transfer_16to8copy(pY_Cur + 8, data[1], stride);  MBDeQuantInter(const MBParam * pParam,
213                          break;                             const int iQuant,
214                  case 2:                             int16_t data[6 * 64],
215                          transfer_16to8copy(pY_Cur + 8 * stride, data[2], stride);                             int16_t qcoeff[6 * 64],
216                          break;                             const uint8_t cbp)
217                  case 3:  {
218                          transfer_16to8copy(pY_Cur + 8 + 8 * stride, data[3], stride);          int i;
219                          break;  
220                  case 4:          for (i = 0; i < 6; i++) {
221                          transfer_16to8copy(pU_Cur, data[4], stride / 2);                  if (cbp & (1 << (5 - i))) {
222                          break;                          start_timer();
223                  case 5:                          if (!(pParam->vol_flags & XVID_VOL_MPEGQUANT))
224                          transfer_16to8copy(pV_Cur, data[5], stride / 2);                                  dequant_inter(&data[i * 64], &qcoeff[i * 64], iQuant);
225                          break;                          else
226                                    dequant4_inter(&data[i * 64], &qcoeff[i * 64], iQuant);
227                            stop_iquant_timer();
228                  }                  }
                 stop_transfer_timer();  
229      }      }
230  }  }
231    
232    typedef void (transfer_operation_8to16_t) (int16_t *Dst, const uint8_t *Src, int BpS);
233    typedef void (transfer_operation_16to8_t) (uint8_t *Dst, const int16_t *Src, int BpS);
234    
 uint8_t MBTransQuantInter(const MBParam *pParam,  
                                         const uint32_t x_pos, const uint32_t y_pos,  
                                         int16_t data[][64],  
                                         int16_t qcoeff[][64],  
                                         IMAGE * const pCurrent)  
235    
236    static __inline void
237    MBTrans8to16(const MBParam * pParam,
238                             FRAMEINFO * frame,
239                             MACROBLOCK * pMB,
240                             const uint32_t x_pos,
241                             const uint32_t y_pos,
242                             int16_t data[6 * 64])
243  {  {
244          const uint32_t stride = pParam->edged_width;          uint32_t stride = pParam->edged_width;
245      uint8_t i;          uint32_t stride2 = stride / 2;
246      uint8_t iQuant = pParam->quant;          uint32_t next_block = stride * 8;
247            int32_t cst;
248          uint8_t *pY_Cur, *pU_Cur, *pV_Cur;          uint8_t *pY_Cur, *pU_Cur, *pV_Cur;
249      uint8_t cbp = 0;          IMAGE *pCurrent = &frame->image;
250          uint32_t sum;          transfer_operation_8to16_t *transfer_op = NULL;
251    
252            if ((frame->vop_flags & XVID_VOP_REDUCED)) {
253    
254                    /* Image pointers */
255                    pY_Cur = pCurrent->y + (y_pos << 5) * stride  + (x_pos << 5);
256                    pU_Cur = pCurrent->u + (y_pos << 4) * stride2 + (x_pos << 4);
257                    pV_Cur = pCurrent->v + (y_pos << 4) * stride2 + (x_pos << 4);
258    
259                    /* Block size */
260                    cst = 16;
261    
262                    /* Operation function */
263                    transfer_op = (transfer_operation_8to16_t*)filter_18x18_to_8x8;
264            } else {
265    
266                    /* Image pointers */
267      pY_Cur = pCurrent->y + (y_pos << 4) * stride + (x_pos << 4);      pY_Cur = pCurrent->y + (y_pos << 4) * stride + (x_pos << 4);
268      pU_Cur = pCurrent->u + (y_pos << 3) * (stride >> 1) + (x_pos << 3);                  pU_Cur = pCurrent->u + (y_pos << 3) * stride2 + (x_pos << 3);
269      pV_Cur = pCurrent->v + (y_pos << 3) * (stride >> 1) + (x_pos << 3);                  pV_Cur = pCurrent->v + (y_pos << 3) * stride2 + (x_pos << 3);
270    
271      for(i = 0; i < 6; i++) {                  /* Block size */
272                  /*                  cst = 8;
273                  no need to transfer 8->16-bit  
274                  (this is performed already in motion compensation)                  /* Operation function */
275                  */                  transfer_op = (transfer_operation_8to16_t*)transfer_8to16copy;
276            }
277    
278            /* Do the transfer */
279                  start_timer();                  start_timer();
280                  fdct(data[i]);          transfer_op(&data[0 * 64], pY_Cur, stride);
281                  stop_dct_timer();          transfer_op(&data[1 * 64], pY_Cur + cst, stride);
282            transfer_op(&data[2 * 64], pY_Cur + next_block, stride);
283            transfer_op(&data[3 * 64], pY_Cur + next_block + cst, stride);
284            transfer_op(&data[4 * 64], pU_Cur, stride2);
285            transfer_op(&data[5 * 64], pV_Cur, stride2);
286            stop_transfer_timer();
287    }
288    
289                  if (pParam->quant_type == 0)  static __inline void
290    MBTrans16to8(const MBParam * pParam,
291                             FRAMEINFO * frame,
292                             MACROBLOCK * pMB,
293                             const uint32_t x_pos,
294                             const uint32_t y_pos,
295                             int16_t data[6 * 64],
296                             const uint32_t add,
297                             const uint8_t cbp)
298                  {                  {
299            uint8_t *pY_Cur, *pU_Cur, *pV_Cur;
300            uint32_t stride = pParam->edged_width;
301            uint32_t stride2 = stride / 2;
302            uint32_t next_block = stride * 8;
303            uint32_t cst;
304            IMAGE *pCurrent = &frame->image;
305            transfer_operation_16to8_t *transfer_op = NULL;
306    
307            if (pMB->field_dct) {
308                    next_block = stride;
309                    stride *= 2;
310            }
311    
312            if ((frame->vop_flags & XVID_VOP_REDUCED)) {
313    
314                    /* Image pointers */
315                    pY_Cur = pCurrent->y + (y_pos << 5) * stride  + (x_pos << 5);
316                    pU_Cur = pCurrent->u + (y_pos << 4) * stride2 + (x_pos << 4);
317                    pV_Cur = pCurrent->v + (y_pos << 4) * stride2 + (x_pos << 4);
318    
319                    /* Block size */
320                    cst = 16;
321    
322                    /* Operation function */
323                    if(add)
324                            transfer_op = (transfer_operation_16to8_t*)add_upsampled_8x8_16to8;
325                    else
326                            transfer_op = (transfer_operation_16to8_t*)copy_upsampled_8x8_16to8;
327            } else {
328    
329                    /* Image pointers */
330                    pY_Cur = pCurrent->y + (y_pos << 4) * stride  + (x_pos << 4);
331                    pU_Cur = pCurrent->u + (y_pos << 3) * stride2 + (x_pos << 3);
332                    pV_Cur = pCurrent->v + (y_pos << 3) * stride2 + (x_pos << 3);
333    
334                    /* Block size */
335                    cst = 8;
336    
337                    /* Operation function */
338                    if(add)
339                            transfer_op = (transfer_operation_16to8_t*)transfer_16to8add;
340                    else
341                            transfer_op = (transfer_operation_16to8_t*)transfer_16to8copy;
342            }
343    
344            /* Do the operation */
345                          start_timer();                          start_timer();
346                          sum = quant_inter(qcoeff[i], data[i], iQuant);          if (cbp&32) transfer_op(pY_Cur, &data[0 * 64], stride);
347                          stop_quant_timer();          if (cbp&16) transfer_op(pY_Cur + cst, &data[1 * 64], stride);
348            if (cbp& 8) transfer_op(pY_Cur + next_block, &data[2 * 64], stride);
349            if (cbp& 4) transfer_op(pY_Cur + next_block + cst, &data[3 * 64], stride);
350            if (cbp& 2) transfer_op(pU_Cur, &data[4 * 64], stride2);
351            if (cbp& 1) transfer_op(pV_Cur, &data[5 * 64], stride2);
352            stop_transfer_timer();
353                  }                  }
354                  else  
355    /*****************************************************************************
356     * Module functions
357     ****************************************************************************/
358    
359    void
360    MBTransQuantIntra(const MBParam * pParam,
361                                      FRAMEINFO * frame,
362                                      MACROBLOCK * pMB,
363                                      const uint32_t x_pos,
364                                      const uint32_t y_pos,
365                                      int16_t data[6 * 64],
366                                      int16_t qcoeff[6 * 64])
367                  {                  {
368                          start_timer();  
369                          sum = quant4_inter(qcoeff[i], data[i], iQuant);          /* Transfer data */
370                          stop_quant_timer();          MBTrans8to16(pParam, frame, pMB, x_pos, y_pos, data);
371    
372            /* Perform DCT (and field decision) */
373            MBfDCT(pParam, frame, pMB, x_pos, y_pos, data);
374    
375            /* Quantize the block */
376            MBQuantIntra(pParam, pMB, data, qcoeff);
377    
378            /* DeQuantize the block */
379            MBDeQuantIntra(pParam, pMB->quant, data, qcoeff);
380    
381            /* Perform inverse DCT*/
382            MBiDCT(data, 0x3F);
383    
384            /* Transfer back the data -- Don't add data */
385            MBTrans16to8(pParam, frame, pMB, x_pos, y_pos, data, 0, 0x3F);
386                  }                  }
387    
                 if(sum >= TOOSMALL_LIMIT) { // skip block ?  
388    
389                          if (pParam->quant_type == H263_QUANT)  uint8_t
390    MBTransQuantInter(const MBParam * pParam,
391                                      FRAMEINFO * frame,
392                                      MACROBLOCK * pMB,
393                                      const uint32_t x_pos,
394                                      const uint32_t y_pos,
395                                      int16_t data[6 * 64],
396                                      int16_t qcoeff[6 * 64])
397                          {                          {
398                                  start_timer();          uint8_t cbp;
399                                  dequant_inter(data[i], qcoeff[i], iQuant);          uint32_t limit;
400                                  stop_iquant_timer();  
401            /*
402             * There is no MBTrans8to16 for Inter block, that's done in motion compensation
403             * already
404             */
405    
406            /* Perform DCT (and field decision) */
407            MBfDCT(pParam, frame, pMB, x_pos, y_pos, data);
408    
409            /* Set the limit threshold */
410            limit = PVOP_TOOSMALL_LIMIT + ((pMB->quant == 1)? 1 : 0);
411    
412            /* Quantize the block */
413            cbp = MBQuantInter(pParam, pMB, data, qcoeff, 0, limit);
414    
415            /* DeQuantize the block */
416            MBDeQuantInter(pParam, pMB->quant, data, qcoeff, cbp);
417    
418            /* Perform inverse DCT*/
419            MBiDCT(data, cbp);
420    
421            /* Transfer back the data -- Add the data */
422            MBTrans16to8(pParam, frame, pMB, x_pos, y_pos, data, 1, cbp);
423    
424            return(cbp);
425                          }                          }
426                          else  
427    uint8_t
428    MBTransQuantInterBVOP(const MBParam * pParam,
429                                      FRAMEINFO * frame,
430                                      MACROBLOCK * pMB,
431                                      const uint32_t x_pos,
432                                      const uint32_t y_pos,
433                                      int16_t data[6 * 64],
434                                      int16_t qcoeff[6 * 64])
435                          {                          {
436                                  start_timer();          uint8_t cbp;
437                                  dequant4_inter(data[i], qcoeff[i], iQuant);          uint32_t limit;
                                 stop_iquant_timer();  
                         }  
438    
439                          cbp |= 1 << (5 - i);          /*
440             * There is no MBTrans8to16 for Inter block, that's done in motion compensation
441             * already
442             */
443    
444                          start_timer();          /* Perform DCT (and field decision) */
445                          idct(data[i]);          MBfDCT(pParam, frame, pMB, x_pos, y_pos, data);
                         stop_idct_timer();  
446    
447                          start_timer();          /* Set the limit threshold */
448            limit = BVOP_TOOSMALL_LIMIT;
449    
450            /* Quantize the block */
451            cbp = MBQuantInter(pParam, pMB, data, qcoeff, 1, limit);
452    
453            /*
454             * History comment:
455             * We don't have to DeQuant, iDCT and Transfer back data for B-frames.
456             *
457             * BUT some plugins require the original frame to be passed so we have
458             * to take care of that here
459             */
460            if((pParam->plugin_flags & XVID_REQORIGINAL)) {
461    
462                    /* DeQuantize the block */
463                    MBDeQuantInter(pParam, pMB->quant, data, qcoeff, cbp);
464    
465                          switch(i) {                  /* Perform inverse DCT*/
466                          case 0:                  MBiDCT(data, cbp);
467                                  transfer_16to8add(pY_Cur, data[0], stride);  
468                                  break;                  /* Transfer back the data -- Add the data */
469                          case 1:                  MBTrans16to8(pParam, frame, pMB, x_pos, y_pos, data, 1, cbp);
                                 transfer_16to8add(pY_Cur + 8, data[1], stride);  
                                 break;  
                         case 2:  
                                 transfer_16to8add(pY_Cur + 8 * stride, data[2], stride);  
                                 break;  
                         case 3:  
                                 transfer_16to8add(pY_Cur + 8 + 8 * stride, data[3], stride);  
                                 break;  
                         case 4:  
                                 transfer_16to8add(pU_Cur, data[4], stride / 2);  
                                 break;  
                         case 5:  
                                 transfer_16to8add(pV_Cur, data[5], stride / 2);  
                                 break;  
470                          }                          }
471                          stop_transfer_timer();  
472            return(cbp);
473    }
474    
475    /* if sum(diff between field lines) < sum(diff between frame lines), use field dct */
476    uint32_t
477    MBFieldTest_c(int16_t data[6 * 64])
478    {
479            const uint8_t blocks[] =
480                    { 0 * 64, 0 * 64, 0 * 64, 0 * 64, 2 * 64, 2 * 64, 2 * 64, 2 * 64 };
481            const uint8_t lines[] = { 0, 16, 32, 48, 0, 16, 32, 48 };
482    
483            int frame = 0, field = 0;
484            int i, j;
485    
486            for (i = 0; i < 7; ++i) {
487                    for (j = 0; j < 8; ++j) {
488                            frame +=
489                                    ABS(data[0 * 64 + (i + 1) * 8 + j] - data[0 * 64 + i * 8 + j]);
490                            frame +=
491                                    ABS(data[1 * 64 + (i + 1) * 8 + j] - data[1 * 64 + i * 8 + j]);
492                            frame +=
493                                    ABS(data[2 * 64 + (i + 1) * 8 + j] - data[2 * 64 + i * 8 + j]);
494                            frame +=
495                                    ABS(data[3 * 64 + (i + 1) * 8 + j] - data[3 * 64 + i * 8 + j]);
496    
497                            field +=
498                                    ABS(data[blocks[i + 1] + lines[i + 1] + j] -
499                                            data[blocks[i] + lines[i] + j]);
500                            field +=
501                                    ABS(data[blocks[i + 1] + lines[i + 1] + 8 + j] -
502                                            data[blocks[i] + lines[i] + 8 + j]);
503                            field +=
504                                    ABS(data[blocks[i + 1] + 64 + lines[i + 1] + j] -
505                                            data[blocks[i] + 64 + lines[i] + j]);
506                            field +=
507                                    ABS(data[blocks[i + 1] + 64 + lines[i + 1] + 8 + j] -
508                                            data[blocks[i] + 64 + lines[i] + 8 + j]);
509                  }                  }
510          }          }
511      return cbp;  
512            return (frame >= (field + 350));
513    }
514    
515    
516    /* deinterlace Y blocks vertically */
517    
518    #define MOVLINE(X,Y) memcpy(X, Y, sizeof(tmp))
519    #define LINE(X,Y)    &data[X*64 + Y*8]
520    
521    void
522    MBFrameToField(int16_t data[6 * 64])
523    {
524            int16_t tmp[8];
525    
526            /* left blocks */
527    
528            // 1=2, 2=4, 4=8, 8=1
529            MOVLINE(tmp, LINE(0, 1));
530            MOVLINE(LINE(0, 1), LINE(0, 2));
531            MOVLINE(LINE(0, 2), LINE(0, 4));
532            MOVLINE(LINE(0, 4), LINE(2, 0));
533            MOVLINE(LINE(2, 0), tmp);
534    
535            // 3=6, 6=12, 12=9, 9=3
536            MOVLINE(tmp, LINE(0, 3));
537            MOVLINE(LINE(0, 3), LINE(0, 6));
538            MOVLINE(LINE(0, 6), LINE(2, 4));
539            MOVLINE(LINE(2, 4), LINE(2, 1));
540            MOVLINE(LINE(2, 1), tmp);
541    
542            // 5=10, 10=5
543            MOVLINE(tmp, LINE(0, 5));
544            MOVLINE(LINE(0, 5), LINE(2, 2));
545            MOVLINE(LINE(2, 2), tmp);
546    
547            // 7=14, 14=13, 13=11, 11=7
548            MOVLINE(tmp, LINE(0, 7));
549            MOVLINE(LINE(0, 7), LINE(2, 6));
550            MOVLINE(LINE(2, 6), LINE(2, 5));
551            MOVLINE(LINE(2, 5), LINE(2, 3));
552            MOVLINE(LINE(2, 3), tmp);
553    
554            /* right blocks */
555    
556            // 1=2, 2=4, 4=8, 8=1
557            MOVLINE(tmp, LINE(1, 1));
558            MOVLINE(LINE(1, 1), LINE(1, 2));
559            MOVLINE(LINE(1, 2), LINE(1, 4));
560            MOVLINE(LINE(1, 4), LINE(3, 0));
561            MOVLINE(LINE(3, 0), tmp);
562    
563            // 3=6, 6=12, 12=9, 9=3
564            MOVLINE(tmp, LINE(1, 3));
565            MOVLINE(LINE(1, 3), LINE(1, 6));
566            MOVLINE(LINE(1, 6), LINE(3, 4));
567            MOVLINE(LINE(3, 4), LINE(3, 1));
568            MOVLINE(LINE(3, 1), tmp);
569    
570            // 5=10, 10=5
571            MOVLINE(tmp, LINE(1, 5));
572            MOVLINE(LINE(1, 5), LINE(3, 2));
573            MOVLINE(LINE(3, 2), tmp);
574    
575            // 7=14, 14=13, 13=11, 11=7
576            MOVLINE(tmp, LINE(1, 7));
577            MOVLINE(LINE(1, 7), LINE(3, 6));
578            MOVLINE(LINE(3, 6), LINE(3, 5));
579            MOVLINE(LINE(3, 5), LINE(3, 3));
580            MOVLINE(LINE(3, 3), tmp);
581  }  }
Legend:
Removed from v.3  
changed lines
  Added in v.965
No admin address has been configured
 ViewVC Help
Powered by ViewVC 1.0.4