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

Diff of /xvidcore/src/encoder.c

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

revision 1.19, Thu Apr 4 13:58:02 2002 UTC revision 1.96, Tue Mar 4 10:55:21 2003 UTC
# Line 1  Line 1 
1    /*****************************************************************************
2     *
3     *  XVID MPEG-4 VIDEO CODEC
4     *  -  Encoder main module  -
5     *
6     *  This program is an implementation of a part of one or more MPEG-4
7     *  Video tools as specified in ISO/IEC 14496-2 standard.  Those intending
8     *  to use this software module in hardware or software products are
9     *  advised that its use may infringe existing patents or copyrights, and
10     *  any such use would be at such party's own risk.  The original
11     *  developer of this software module and his/her company, and subsequent
12     *  editors and their companies, will have no liability for use of this
13     *  software or modifications or derivatives thereof.
14     *
15     *  This program is free software; you can redistribute it and/or modify
16     *  it under the terms of the GNU General Public License as published by
17     *  the Free Software Foundation; either version 2 of the License, or
18     *  (at your option) any later version.
19     *
20     *  This program is distributed in the hope that it will be useful,
21     *  but WITHOUT ANY WARRANTY; without even the implied warranty of
22     *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
23     *  GNU General Public License for more details.
24     *
25     *  You should have received a copy of the GNU General Public License
26     *  along with this program; if not, write to the Free Software
27     *  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
28     *
29     *  $Id$
30     *
31     ****************************************************************************/
32    
33  #include <stdlib.h>  #include <stdlib.h>
34  #include <stdio.h>  #include <stdio.h>
35  #include <math.h>  #include <math.h>
36    #include <string.h>
37    
38  #include "encoder.h"  #include "encoder.h"
39  #include "prediction/mbprediction.h"  #include "prediction/mbprediction.h"
40  #include "global.h"  #include "global.h"
41  #include "utils/timer.h"  #include "utils/timer.h"
42  #include "image/image.h"  #include "image/image.h"
43    #include "image/font.h"
44    #include "motion/sad.h"
45    #include "motion/motion.h"
46  #include "bitstream/cbp.h"  #include "bitstream/cbp.h"
47  #include "utils/mbfunctions.h"  #include "utils/mbfunctions.h"
48  #include "bitstream/bitstream.h"  #include "bitstream/bitstream.h"
# Line 18  Line 54 
54  #include "quant/quant_matrix.h"  #include "quant/quant_matrix.h"
55  #include "utils/mem_align.h"  #include "utils/mem_align.h"
56    
57  #define ENC_CHECK(X) if(!(X)) return XVID_ERR_FORMAT  /*****************************************************************************
58     * Local macros
59     ****************************************************************************/
60    
61    #define ENC_CHECK(X) if(!(X)) return XVID_ERR_FORMAT
62    #define SWAP(_T_,A,B)   { _T_ tmp = A; A = B; B = tmp; }
63    
64  static int FrameCodeI(Encoder * pEnc, Bitstream * bs, uint32_t *pBits);  /*****************************************************************************
65  static int FrameCodeP(Encoder * pEnc, Bitstream * bs, uint32_t *pBits, bool force_inter, bool vol_header);   * Local function prototypes
66     ****************************************************************************/
67    
68    static int FrameCodeI(Encoder * pEnc,
69                                              Bitstream * bs,
70                                              uint32_t * pBits);
71    
72    static int FrameCodeP(Encoder * pEnc,
73                                              Bitstream * bs,
74                                              uint32_t * pBits,
75                                              bool force_inter,
76                                              bool vol_header);
77    
78    static void FrameCodeB(Encoder * pEnc,
79                                               FRAMEINFO * frame,
80                                               Bitstream * bs,
81                                               uint32_t * pBits);
82    
83    /*****************************************************************************
84     * Local data
85     ****************************************************************************/
86    
87  static int DQtab[4] =  static int DQtab[4] = {
 {  
88          -1, -2, 1, 2          -1, -2, 1, 2
89  };  };
90    
91  static int iDQtab[5] =  static int iDQtab[5] = {
 {  
92          1, 0, NO_CHANGE, 2, 3          1, 0, NO_CHANGE, 2, 3
93  };  };
94    
95    
96  int encoder_create(XVID_ENC_PARAM * pParam)  /*****************************************************************************
97     * Encoder creation
98     *
99     * This function creates an Encoder instance, it allocates all necessary
100     * image buffers (reference, current and bframes) and initialize the internal
101     * xvid encoder paremeters according to the XVID_ENC_PARAM input parameter.
102     *
103     * The code seems to be very long but is very basic, mainly memory allocation
104     * and cleaning code.
105     *
106     * Returned values :
107     *      - XVID_ERR_OK    - no errors
108     *      - XVID_ERR_MEMORY - the libc could not allocate memory, the function
109     *                                              cleans the structure before exiting.
110     *                                              pParam->handle is also set to NULL.
111     *
112     ****************************************************************************/
113    
114    int
115    encoder_create(XVID_ENC_PARAM * pParam)
116  {  {
117          Encoder *pEnc;          Encoder *pEnc;
118          uint32_t i;          int i;
   
119          pParam->handle = NULL;          pParam->handle = NULL;
120    
121          ENC_CHECK(pParam);          ENC_CHECK(pParam);
# Line 49  Line 125 
125          ENC_CHECK(!(pParam->width % 2));          ENC_CHECK(!(pParam->width % 2));
126          ENC_CHECK(!(pParam->height % 2));          ENC_CHECK(!(pParam->height % 2));
127    
128          if (pParam->fincr <= 0 || pParam->fbase <= 0)          /* Fps */
129          {  
130            if (pParam->fincr <= 0 || pParam->fbase <= 0) {
131                  pParam->fincr = 1;                  pParam->fincr = 1;
132                  pParam->fbase = 25;                  pParam->fbase = 25;
133          }          }
134    
135          // simplify the "fincr/fbase" fraction          /*
136          // (neccessary, since windows supplies us with huge numbers)           * Simplify the "fincr/fbase" fraction
137             * (neccessary, since windows supplies us with huge numbers)
138             */
139    
140          i = pParam->fincr;          i = pParam->fincr;
141          while (i > 1)          while (i > 1) {
142          {                  if (pParam->fincr % i == 0 && pParam->fbase % i == 0) {
                 if (pParam->fincr % i == 0 && pParam->fbase % i == 0)  
                 {  
143                          pParam->fincr /= i;                          pParam->fincr /= i;
144                          pParam->fbase /= i;                          pParam->fbase /= i;
145                          i = pParam->fincr;                          i = pParam->fincr;
# Line 71  Line 148 
148                  i--;                  i--;
149          }          }
150    
151          if (pParam->fbase > 65535)          if (pParam->fbase > 65535) {
         {  
152                  float div = (float)pParam->fbase / 65535;                  float div = (float)pParam->fbase / 65535;
153    
154                  pParam->fbase = (int)(pParam->fbase / div);                  pParam->fbase = (int)(pParam->fbase / div);
155                  pParam->fincr = (int)(pParam->fincr / div);                  pParam->fincr = (int)(pParam->fincr / div);
156          }          }
157    
158          if (pParam->bitrate <= 0)          /* Bitrate allocator defaults */
159                  pParam->bitrate = 900000;  
160            if (pParam->rc_bitrate <= 0)
161                    pParam->rc_bitrate = 900000;
162    
163            if (pParam->rc_reaction_delay_factor <= 0)
164                    pParam->rc_reaction_delay_factor = 16;
165    
166            if (pParam->rc_averaging_period <= 0)
167                    pParam->rc_averaging_period = 100;
168    
169          if (pParam->rc_buffersize <= 0)          if (pParam->rc_buffer <= 0)
170                  pParam->rc_buffersize = 16;                  pParam->rc_buffer = 100;
171    
172            /* Max and min quantizers */
173    
174          if ((pParam->min_quantizer <= 0) || (pParam->min_quantizer > 31))          if ((pParam->min_quantizer <= 0) || (pParam->min_quantizer > 31))
175                  pParam->min_quantizer = 1;                  pParam->min_quantizer = 1;
# Line 90  Line 177 
177          if ((pParam->max_quantizer <= 0) || (pParam->max_quantizer > 31))          if ((pParam->max_quantizer <= 0) || (pParam->max_quantizer > 31))
178                  pParam->max_quantizer = 31;                  pParam->max_quantizer = 31;
179    
         if (pParam->max_key_interval == 0)              /* 1 keyframe each 10 seconds */  
                 pParam->max_key_interval = 10 * pParam->fincr / pParam->fbase;  
   
180          if (pParam->max_quantizer < pParam->min_quantizer)          if (pParam->max_quantizer < pParam->min_quantizer)
181                  pParam->max_quantizer = pParam->min_quantizer;                  pParam->max_quantizer = pParam->min_quantizer;
182    
183          if ((pEnc = (Encoder *) xvid_malloc(sizeof(Encoder), CACHE_LINE)) == NULL)          /* 1 keyframe each 10 seconds */
184    
185            if (pParam->max_key_interval <= 0)
186                    pParam->max_key_interval = 10 * pParam->fincr / pParam->fbase;
187    
188            pEnc = (Encoder *) xvid_malloc(sizeof(Encoder), CACHE_LINE);
189            if (pEnc == NULL)
190                  return XVID_ERR_MEMORY;                  return XVID_ERR_MEMORY;
191    
192            /* Zero the Encoder Structure */
193    
194            memset(pEnc, 0, sizeof(Encoder));
195    
196          /* Fill members of Encoder structure */          /* Fill members of Encoder structure */
197    
198          pEnc->mbParam.width = pParam->width;          pEnc->mbParam.width = pParam->width;
# Line 110  Line 204 
204          pEnc->mbParam.edged_width = 16 * pEnc->mbParam.mb_width + 2 * EDGE_SIZE;          pEnc->mbParam.edged_width = 16 * pEnc->mbParam.mb_width + 2 * EDGE_SIZE;
205          pEnc->mbParam.edged_height = 16 * pEnc->mbParam.mb_height + 2 * EDGE_SIZE;          pEnc->mbParam.edged_height = 16 * pEnc->mbParam.mb_height + 2 * EDGE_SIZE;
206    
207          pEnc->sStat.fMvPrevSigma = -1;          pEnc->mbParam.fbase = pParam->fbase;
208            pEnc->mbParam.fincr = pParam->fincr;
209    
210            pEnc->mbParam.m_quant_type = H263_QUANT;
211    
212            pEnc->fMvPrevSigma = -1;
213    
214          /* Fill rate control parameters */          /* Fill rate control parameters */
215    
216          pEnc->mbParam.quant = 4;          pEnc->bitrate = pParam->rc_bitrate;
217    
218          pEnc->bitrate = pParam->bitrate;          pEnc->iFrameNum = -1;
219            pEnc->mbParam.iMaxKeyInterval = pParam->max_key_interval;
220    
221          pEnc->iFrameNum = 0;          /* try to allocate frame memory */
         pEnc->iMaxKeyInterval = pParam->max_key_interval;  
222    
223          /* try to allocate memory */          pEnc->current = xvid_malloc(sizeof(FRAMEINFO), CACHE_LINE);
224            pEnc->reference = xvid_malloc(sizeof(FRAMEINFO), CACHE_LINE);
225    
226            if (pEnc->current == NULL || pEnc->reference == NULL)
227                    goto xvid_err_memory1;
228    
229            /* try to allocate mb memory */
230    
231            pEnc->current->mbs =
232                    xvid_malloc(sizeof(MACROBLOCK) * pEnc->mbParam.mb_width *
233                                            pEnc->mbParam.mb_height, CACHE_LINE);
234            pEnc->reference->mbs =
235                    xvid_malloc(sizeof(MACROBLOCK) * pEnc->mbParam.mb_width *
236                                            pEnc->mbParam.mb_height, CACHE_LINE);
237    
238            if (pEnc->current->mbs == NULL || pEnc->reference->mbs == NULL)
239                    goto xvid_err_memory2;
240    
241            /* try to allocate image memory */
242    
243            if (pParam->global & XVID_GLOBAL_EXTRASTATS)
244                    image_null(&pEnc->sOriginal);
245    
246            image_null(&pEnc->f_refh);
247            image_null(&pEnc->f_refv);
248            image_null(&pEnc->f_refhv);
249    
250            image_null(&pEnc->current->image);
251            image_null(&pEnc->reference->image);
252            image_null(&pEnc->vInterH);
253            image_null(&pEnc->vInterV);
254            image_null(&pEnc->vInterVf);
255            image_null(&pEnc->vInterHV);
256            image_null(&pEnc->vInterHVf);
257    
258            if (pParam->global & XVID_GLOBAL_EXTRASTATS)
259            {       if (image_create
260                            (&pEnc->sOriginal, pEnc->mbParam.edged_width,
261                             pEnc->mbParam.edged_height) < 0)
262                            goto xvid_err_memory3;
263            }
264    
265            if (image_create
266                    (&pEnc->f_refh, pEnc->mbParam.edged_width,
267                     pEnc->mbParam.edged_height) < 0)
268                    goto xvid_err_memory3;
269            if (image_create
270                    (&pEnc->f_refv, pEnc->mbParam.edged_width,
271                     pEnc->mbParam.edged_height) < 0)
272                    goto xvid_err_memory3;
273            if (image_create
274                    (&pEnc->f_refhv, pEnc->mbParam.edged_width,
275                     pEnc->mbParam.edged_height) < 0)
276                    goto xvid_err_memory3;
277    
278            if (image_create
279                    (&pEnc->current->image, pEnc->mbParam.edged_width,
280                     pEnc->mbParam.edged_height) < 0)
281                    goto xvid_err_memory3;
282            if (image_create
283                    (&pEnc->reference->image, pEnc->mbParam.edged_width,
284                     pEnc->mbParam.edged_height) < 0)
285                    goto xvid_err_memory3;
286            if (image_create
287                    (&pEnc->vInterH, pEnc->mbParam.edged_width,
288                     pEnc->mbParam.edged_height) < 0)
289                    goto xvid_err_memory3;
290            if (image_create
291                    (&pEnc->vInterV, pEnc->mbParam.edged_width,
292                     pEnc->mbParam.edged_height) < 0)
293                    goto xvid_err_memory3;
294            if (image_create
295                    (&pEnc->vInterVf, pEnc->mbParam.edged_width,
296                     pEnc->mbParam.edged_height) < 0)
297                    goto xvid_err_memory3;
298            if (image_create
299                    (&pEnc->vInterHV, pEnc->mbParam.edged_width,
300                     pEnc->mbParam.edged_height) < 0)
301                    goto xvid_err_memory3;
302            if (image_create
303                    (&pEnc->vInterHVf, pEnc->mbParam.edged_width,
304                     pEnc->mbParam.edged_height) < 0)
305                    goto xvid_err_memory3;
306    
307    /* Create full bitplane for GMC, this might be wasteful */
308            if (image_create
309                    (&pEnc->vGMC, pEnc->mbParam.edged_width,
310                     pEnc->mbParam.edged_height) < 0)
311                    goto xvid_err_memory3;
312    
313    
314    
315            pEnc->mbParam.global = pParam->global;
316    
317            /* B Frames specific init */
318            pEnc->mbParam.max_bframes = pParam->max_bframes;
319            pEnc->mbParam.bquant_ratio = pParam->bquant_ratio;
320            pEnc->mbParam.bquant_offset = pParam->bquant_offset;
321            pEnc->mbParam.frame_drop_ratio = pParam->frame_drop_ratio;
322            pEnc->bframes = NULL;
323    
324            if (pEnc->mbParam.max_bframes > 0) {
325                    int n;
326    
327                    pEnc->bframes =
328                            xvid_malloc(pEnc->mbParam.max_bframes * sizeof(FRAMEINFO *),
329                                                    CACHE_LINE);
330    
331                    if (pEnc->bframes == NULL)
332                            goto xvid_err_memory3;
333    
334                    for (n = 0; n < pEnc->mbParam.max_bframes; n++)
335                            pEnc->bframes[n] = NULL;
336    
337    
338                    for (n = 0; n < pEnc->mbParam.max_bframes; n++) {
339                            pEnc->bframes[n] = xvid_malloc(sizeof(FRAMEINFO), CACHE_LINE);
340    
341                            if (pEnc->bframes[n] == NULL)
342                                    goto xvid_err_memory4;
343    
344                            pEnc->bframes[n]->mbs =
345                                    xvid_malloc(sizeof(MACROBLOCK) * pEnc->mbParam.mb_width *
346                                                            pEnc->mbParam.mb_height, CACHE_LINE);
347    
348                            if (pEnc->bframes[n]->mbs == NULL)
349                                    goto xvid_err_memory4;
350    
351                            image_null(&pEnc->bframes[n]->image);
352    
353                            if (image_create
354                                    (&pEnc->bframes[n]->image, pEnc->mbParam.edged_width,
355                                     pEnc->mbParam.edged_height) < 0)
356                                    goto xvid_err_memory4;
357    
         pEnc->sCurrent.y        =       pEnc->sCurrent.u        =       pEnc->sCurrent.v        = NULL;  
         pEnc->sReference.y      =       pEnc->sReference.u      =       pEnc->sReference.v      = NULL;  
         pEnc->vInterH.y         =       pEnc->vInterH.u         =       pEnc->vInterH.v         = NULL;  
         pEnc->vInterV.y         =       pEnc->vInterV.u         =       pEnc->vInterV.v         = NULL;  
         pEnc->vInterVf.y        =       pEnc->vInterVf.u        =       pEnc->vInterVf.v        = NULL;  
         pEnc->vInterHV.y        =       pEnc->vInterHV.u        =       pEnc->vInterHV.v        = NULL;  
         pEnc->vInterHVf.y       =       pEnc->vInterHVf.u       =       pEnc->vInterHVf.v       = NULL;  
   
         pEnc->pMBs = NULL;  
   
         if (image_create(&pEnc->sCurrent, pEnc->mbParam.edged_width, pEnc->mbParam.edged_height) < 0 ||  
                 image_create(&pEnc->sReference, pEnc->mbParam.edged_width, pEnc->mbParam.edged_height) < 0 ||  
                 image_create(&pEnc->vInterH, pEnc->mbParam.edged_width, pEnc->mbParam.edged_height) < 0 ||  
                 image_create(&pEnc->vInterV, pEnc->mbParam.edged_width, pEnc->mbParam.edged_height) < 0 ||  
                 image_create(&pEnc->vInterVf, pEnc->mbParam.edged_width, pEnc->mbParam.edged_height) < 0 ||  
                 image_create(&pEnc->vInterHV, pEnc->mbParam.edged_width, pEnc->mbParam.edged_height) < 0 ||  
                 image_create(&pEnc->vInterHVf, pEnc->mbParam.edged_width, pEnc->mbParam.edged_height) < 0 ||  
                 (pEnc->pMBs = xvid_malloc(sizeof(MACROBLOCK) * pEnc->mbParam.mb_width * pEnc->mbParam.mb_height, CACHE_LINE)) == NULL)  
         {  
                 image_destroy(&pEnc->sCurrent, pEnc->mbParam.edged_width, pEnc->mbParam.edged_height);  
                 image_destroy(&pEnc->sReference, pEnc->mbParam.edged_width, pEnc->mbParam.edged_height);  
                 image_destroy(&pEnc->vInterH, pEnc->mbParam.edged_width, pEnc->mbParam.edged_height);  
                 image_destroy(&pEnc->vInterV, pEnc->mbParam.edged_width, pEnc->mbParam.edged_height);  
                 image_destroy(&pEnc->vInterVf, pEnc->mbParam.edged_width, pEnc->mbParam.edged_height);  
                 image_destroy(&pEnc->vInterHV, pEnc->mbParam.edged_width, pEnc->mbParam.edged_height);  
                 image_destroy(&pEnc->vInterHVf, pEnc->mbParam.edged_width, pEnc->mbParam.edged_height);  
                 if (pEnc)  
                 {  
                         xvid_free(pEnc);  
358                  }                  }
                 return XVID_ERR_MEMORY;  
359          }          }
360    
361          // init macroblock array          pEnc->bframenum_head = 0;
362          for (i = 0; i < pEnc->mbParam.mb_width * pEnc->mbParam.mb_height; i++)          pEnc->bframenum_tail = 0;
363          {          pEnc->flush_bframes = 0;
364                  pEnc->pMBs[i].dquant = NO_CHANGE;          pEnc->bframenum_dx50bvop = -1;
365    
366            pEnc->queue = NULL;
367    
368            if (pEnc->mbParam.max_bframes > 0) {
369                    int n;
370    
371                    pEnc->queue =
372                            xvid_malloc(pEnc->mbParam.max_bframes * sizeof(IMAGE),
373                                                    CACHE_LINE);
374    
375                    if (pEnc->queue == NULL)
376                            goto xvid_err_memory4;
377    
378                    for (n = 0; n < pEnc->mbParam.max_bframes; n++)
379                            image_null(&pEnc->queue[n]);
380    
381                    for (n = 0; n < pEnc->mbParam.max_bframes; n++) {
382                            if (image_create
383                                    (&pEnc->queue[n], pEnc->mbParam.edged_width,
384                                     pEnc->mbParam.edged_height) < 0)
385                                    goto xvid_err_memory5;
386    
387          }          }
388            }
389    
390            pEnc->queue_head = 0;
391            pEnc->queue_tail = 0;
392            pEnc->queue_size = 0;
393    
394            pEnc->mbParam.m_stamp = 0;
395    
396            pEnc->m_framenum = 0;
397            pEnc->current->stamp = 0;
398            pEnc->reference->stamp = 0;
399    
400          pParam->handle = (void *)pEnc;          pParam->handle = (void *)pEnc;
401    
402          if (pParam->bitrate)          if (pParam->rc_bitrate) {
403          {                  RateControlInit(&pEnc->rate_control, pParam->rc_bitrate,
404                  RateControlInit(pParam->bitrate, pParam->rc_buffersize, pParam->fbase * 100 / pParam->fincr,                                                  pParam->rc_reaction_delay_factor,
405                                                    pParam->rc_averaging_period, pParam->rc_buffer,
406                                                    pParam->fbase * 1000 / pParam->fincr,
407                                  pParam->max_quantizer, pParam->min_quantizer);                                  pParam->max_quantizer, pParam->min_quantizer);
408          }          }
409    
410          init_timer();          init_timer();
411    
412          return XVID_ERR_OK;          return XVID_ERR_OK;
 }  
413    
414            /*
415             * We handle all XVID_ERR_MEMORY here, this makes the code lighter
416             */
417    
418  int encoder_destroy(Encoder * pEnc)    xvid_err_memory5:
 {  
         ENC_CHECK(pEnc);  
         ENC_CHECK(pEnc->sCurrent.y);  
         ENC_CHECK(pEnc->sReference.y);  
419    
420          xvid_free(pEnc->pMBs);          if (pEnc->mbParam.max_bframes > 0) {
         image_destroy(&pEnc->sCurrent, pEnc->mbParam.edged_width, pEnc->mbParam.edged_height);  
         image_destroy(&pEnc->sReference, pEnc->mbParam.edged_width, pEnc->mbParam.edged_height);  
         image_destroy(&pEnc->vInterH, pEnc->mbParam.edged_width, pEnc->mbParam.edged_height);  
         image_destroy(&pEnc->vInterV, pEnc->mbParam.edged_width, pEnc->mbParam.edged_height);  
         image_destroy(&pEnc->vInterVf, pEnc->mbParam.edged_width, pEnc->mbParam.edged_height);  
         image_destroy(&pEnc->vInterHV, pEnc->mbParam.edged_width, pEnc->mbParam.edged_height);  
         image_destroy(&pEnc->vInterHVf, pEnc->mbParam.edged_width, pEnc->mbParam.edged_height);  
         xvid_free(pEnc);  
421    
422          return XVID_ERR_OK;                  for (i = 0; i < pEnc->mbParam.max_bframes; i++) {
423                            image_destroy(&pEnc->queue[i], pEnc->mbParam.edged_width,
424                                                      pEnc->mbParam.edged_height);
425                    }
426                    xvid_free(pEnc->queue);
427  }  }
428    
429  int encoder_encode(Encoder * pEnc, XVID_ENC_FRAME * pFrame, XVID_ENC_STATS * pResult)    xvid_err_memory4:
 {  
         uint16_t x, y;  
         Bitstream bs;  
         uint32_t bits;  
         uint16_t write_vol_header = 0;  
430    
431          start_global_timer();          if (pEnc->mbParam.max_bframes > 0) {
432    
433          ENC_CHECK(pEnc);                  for (i = 0; i < pEnc->mbParam.max_bframes; i++) {
434          ENC_CHECK(pFrame);  
435          ENC_CHECK(pFrame->bitstream);                          if (pEnc->bframes[i] == NULL)
436          ENC_CHECK(pFrame->image);                                  continue;
437    
438          pEnc->mbParam.global_flags = pFrame->general;                          image_destroy(&pEnc->bframes[i]->image, pEnc->mbParam.edged_width,
439          pEnc->mbParam.motion_flags = pFrame->motion;                                                    pEnc->mbParam.edged_height);
440    
441                            xvid_free(pEnc->bframes[i]->mbs);
442    
443                            xvid_free(pEnc->bframes[i]);
444    
         start_timer();  
         if (image_input(&pEnc->sCurrent, pEnc->mbParam.width, pEnc->mbParam.height, pEnc->mbParam.edged_width,  
                         pFrame->image, pFrame->colorspace))  
         {  
                 return XVID_ERR_FORMAT;  
445          }          }
         stop_conv_timer();  
446    
447          EMMS();                  xvid_free(pEnc->bframes);
448            }
449    
450          BitstreamInit(&bs, pFrame->bitstream, 0);    xvid_err_memory3:
451    
452          if (pFrame->quant == 0)          if (pEnc->mbParam.global & XVID_GLOBAL_EXTRASTATS)
453          {          {       image_destroy(&pEnc->sOriginal, pEnc->mbParam.edged_width,
454                  pEnc->mbParam.quant = RateControlGetQ(0);                                            pEnc->mbParam.edged_height);
455          }          }
456          else  
457          {          image_destroy(&pEnc->f_refh, pEnc->mbParam.edged_width,
458                  pEnc->mbParam.quant = pFrame->quant;                                    pEnc->mbParam.edged_height);
459            image_destroy(&pEnc->f_refv, pEnc->mbParam.edged_width,
460                                      pEnc->mbParam.edged_height);
461            image_destroy(&pEnc->f_refhv, pEnc->mbParam.edged_width,
462                                      pEnc->mbParam.edged_height);
463    
464            image_destroy(&pEnc->current->image, pEnc->mbParam.edged_width,
465                                      pEnc->mbParam.edged_height);
466            image_destroy(&pEnc->reference->image, pEnc->mbParam.edged_width,
467                                      pEnc->mbParam.edged_height);
468            image_destroy(&pEnc->vInterH, pEnc->mbParam.edged_width,
469                                      pEnc->mbParam.edged_height);
470            image_destroy(&pEnc->vInterV, pEnc->mbParam.edged_width,
471                                      pEnc->mbParam.edged_height);
472            image_destroy(&pEnc->vInterVf, pEnc->mbParam.edged_width,
473                                      pEnc->mbParam.edged_height);
474            image_destroy(&pEnc->vInterHV, pEnc->mbParam.edged_width,
475                                      pEnc->mbParam.edged_height);
476            image_destroy(&pEnc->vInterHVf, pEnc->mbParam.edged_width,
477                                      pEnc->mbParam.edged_height);
478    
479    /* destroy GMC image */
480            image_destroy(&pEnc->vGMC, pEnc->mbParam.edged_width,
481                                      pEnc->mbParam.edged_height);
482    
483    
484      xvid_err_memory2:
485            xvid_free(pEnc->current->mbs);
486            xvid_free(pEnc->reference->mbs);
487    
488      xvid_err_memory1:
489            xvid_free(pEnc->current);
490            xvid_free(pEnc->reference);
491            xvid_free(pEnc);
492    
493            pParam->handle = NULL;
494    
495            return XVID_ERR_MEMORY;
496          }          }
497    
498          if ((pEnc->mbParam.global_flags & XVID_LUMIMASKING) > 0)  /*****************************************************************************
499     * Encoder destruction
500     *
501     * This function destroy the entire encoder structure created by a previous
502     * successful encoder_create call.
503     *
504     * Returned values (for now only one returned value) :
505     *      - XVID_ERR_OK    - no errors
506     *
507     ****************************************************************************/
508    
509    int
510    encoder_destroy(Encoder * pEnc)
511          {          {
512                  int * temp_dquants = (int *) xvid_malloc(pEnc->mbParam.mb_width * pEnc->mbParam.mb_height * sizeof(int), CACHE_LINE);          int i;
513    
514                  pEnc->mbParam.quant = adaptive_quantization(pEnc->sCurrent.y,          ENC_CHECK(pEnc);
                                                             pEnc->mbParam.width,  
                                                             temp_dquants,  
                                                             pEnc->mbParam.quant,  
                                                             pEnc->mbParam.quant,  
                                                             2*pEnc->mbParam.quant,  
                                                             pEnc->mbParam.mb_width,  
                                                             pEnc->mbParam.mb_height);  
515    
516                  for (y = 0; y < pEnc->mbParam.mb_height; y++)          /* B Frames specific */
517                          for (x = 0; x < pEnc->mbParam.mb_width; x++)          if (pEnc->mbParam.max_bframes > 0) {
518                          {  
519                                  MACROBLOCK *pMB = &pEnc->pMBs[x + y * pEnc->mbParam.mb_width];                  for (i = 0; i < pEnc->mbParam.max_bframes; i++) {
520                                  pMB->dquant = iDQtab[(temp_dquants[y * pEnc->mbParam.mb_width + x] + 2)];  
521                            image_destroy(&pEnc->queue[i], pEnc->mbParam.edged_width,
522                                              pEnc->mbParam.edged_height);
523                          }                          }
524                  xvid_free(temp_dquants);                  xvid_free(pEnc->queue);
525          }          }
526    
527          if(pEnc->mbParam.global_flags & XVID_H263QUANT) {          if (pEnc->mbParam.max_bframes > 0) {
                 if(pEnc->mbParam.quant_type != H263_QUANT)  
                         write_vol_header = 1;  
                 pEnc->mbParam.quant_type = H263_QUANT;  
         }  
         else if(pEnc->mbParam.global_flags & XVID_MPEGQUANT) {  
                 int ret1, ret2;  
528    
529                  ret1 = ret2 = 0;                  for (i = 0; i < pEnc->mbParam.max_bframes; i++) {
530    
531                  if(pEnc->mbParam.quant_type != MPEG4_QUANT)                          if (pEnc->bframes[i] == NULL)
532                          write_vol_header = 1;                                  continue;
533    
534                  pEnc->mbParam.quant_type = MPEG4_QUANT;                          image_destroy(&pEnc->bframes[i]->image, pEnc->mbParam.edged_width,
535                                              pEnc->mbParam.edged_height);
536    
537                  if ((pEnc->mbParam.global_flags & XVID_CUSTOM_QMATRIX) > 0) {                          xvid_free(pEnc->bframes[i]->mbs);
538                          if(pFrame->quant_intra_matrix != NULL)  
539                                  ret1 = set_intra_matrix(pFrame->quant_intra_matrix);                          xvid_free(pEnc->bframes[i]);
                         if(pFrame->quant_inter_matrix != NULL)  
                                 ret2 = set_inter_matrix(pFrame->quant_inter_matrix);  
                 }  
                 else {  
                         ret1 = set_intra_matrix(get_default_intra_matrix());  
                         ret2 = set_inter_matrix(get_default_inter_matrix());  
540                  }                  }
541                  if(write_vol_header == 0)  
542                          write_vol_header = ret1 | ret2;                  xvid_free(pEnc->bframes);
543    
544          }          }
545    
546          if (pFrame->intra < 0)          /* All images, reference, current etc ... */
         {  
                 if ((pEnc->iFrameNum == 0) || ((pEnc->iMaxKeyInterval > 0)  
                                                && (pEnc->iFrameNum >= pEnc->iMaxKeyInterval)))  
547    
548                          pFrame->intra = FrameCodeI(pEnc, &bs, &bits);          image_destroy(&pEnc->current->image, pEnc->mbParam.edged_width,
549                  else                                    pEnc->mbParam.edged_height);
550                          pFrame->intra = FrameCodeP(pEnc, &bs, &bits, 0, write_vol_header);          image_destroy(&pEnc->reference->image, pEnc->mbParam.edged_width,
551                                      pEnc->mbParam.edged_height);
552            image_destroy(&pEnc->vInterH, pEnc->mbParam.edged_width,
553                                      pEnc->mbParam.edged_height);
554            image_destroy(&pEnc->vInterV, pEnc->mbParam.edged_width,
555                                      pEnc->mbParam.edged_height);
556            image_destroy(&pEnc->vInterVf, pEnc->mbParam.edged_width,
557                                      pEnc->mbParam.edged_height);
558            image_destroy(&pEnc->vInterHV, pEnc->mbParam.edged_width,
559                                      pEnc->mbParam.edged_height);
560            image_destroy(&pEnc->vInterHVf, pEnc->mbParam.edged_width,
561                                      pEnc->mbParam.edged_height);
562    
563            image_destroy(&pEnc->f_refh, pEnc->mbParam.edged_width,
564                                      pEnc->mbParam.edged_height);
565            image_destroy(&pEnc->f_refv, pEnc->mbParam.edged_width,
566                                      pEnc->mbParam.edged_height);
567            image_destroy(&pEnc->f_refhv, pEnc->mbParam.edged_width,
568                                      pEnc->mbParam.edged_height);
569    
570            if (pEnc->mbParam.global & XVID_GLOBAL_EXTRASTATS)
571            {       image_destroy(&pEnc->sOriginal, pEnc->mbParam.edged_width,
572                                              pEnc->mbParam.edged_height);
573          }          }
574          else  
575          {          /* Encoder structure */
576                  if (pFrame->intra == 1)  
577                          pFrame->intra = FrameCodeI(pEnc, &bs, &bits);          xvid_free(pEnc->current->mbs);
578                  else          xvid_free(pEnc->current);
579                          pFrame->intra = FrameCodeP(pEnc, &bs, &bits, 1, write_vol_header);  
580            xvid_free(pEnc->reference->mbs);
581            xvid_free(pEnc->reference);
582    
583            xvid_free(pEnc);
584    
585            return XVID_ERR_OK;
586          }          }
587    
         BitstreamPutBits(&bs, 0xFFFF, 16);  
         BitstreamPutBits(&bs, 0xFFFF, 16);  
         BitstreamPad(&bs);  
         pFrame->length = BitstreamLength(&bs);  
588    
589          if (pResult)  static __inline void inc_frame_num(Encoder * pEnc)
590          {          {
591                  pResult->quant = pEnc->mbParam.quant;          pEnc->current->stamp = pEnc->mbParam.m_stamp;   /* first frame is zero */
592                  pResult->hlength = pFrame->length - (pEnc->sStat.iTextBits / 8);          pEnc->mbParam.m_stamp += pEnc->mbParam.fincr;
                 pResult->kblks = pEnc->sStat.kblks;  
                 pResult->mblks = pEnc->sStat.mblks;  
                 pResult->ublks = pEnc->sStat.ublks;  
593          }          }
594    
         EMMS();  
595    
596          if (pFrame->quant == 0)  static __inline void
597    queue_image(Encoder * pEnc, XVID_ENC_FRAME * pFrame)
598    {
599            if (pEnc->queue_size >= pEnc->mbParam.max_bframes)
600          {          {
601                  RateControlUpdate(pEnc->mbParam.quant, pFrame->length, pFrame->intra);                  DPRINTF(DPRINTF_DEBUG,"FATAL: QUEUE FULL");
602                    return;
603          }          }
604    
605          pEnc->iFrameNum++;          DPRINTF(DPRINTF_DEBUG,"*** QUEUE bf: head=%i tail=%i   queue: head=%i tail=%i size=%i",
606          image_swap(&pEnc->sCurrent, &pEnc->sReference);                                  pEnc->bframenum_head, pEnc->bframenum_tail,
607                                    pEnc->queue_head, pEnc->queue_tail, pEnc->queue_size);
608    
         stop_global_timer();  
         write_timer();  
609    
610          return XVID_ERR_OK;          start_timer();
611            if (image_input
612                    (&pEnc->queue[pEnc->queue_tail], pEnc->mbParam.width, pEnc->mbParam.height,
613                     pEnc->mbParam.edged_width, pFrame->image, pFrame->stride, pFrame->colorspace, pFrame->general & XVID_INTERLACING))
614                    return;
615            stop_conv_timer();
616    
617            if ((pFrame->general & XVID_CHROMAOPT)) {
618                    image_chroma_optimize(&pEnc->queue[pEnc->queue_tail],
619                            pEnc->mbParam.width, pEnc->mbParam.height, pEnc->mbParam.edged_width);
620  }  }
621    
622            pEnc->queue_size++;
623            pEnc->queue_tail =  (pEnc->queue_tail + 1) % pEnc->mbParam.max_bframes;
624    }
625    
626  static __inline void CodeIntraMB(Encoder *pEnc, MACROBLOCK *pMB) {  static __inline void
627    set_timecodes(FRAMEINFO* pCur,FRAMEINFO *pRef, int32_t time_base)
628    {
629    
630          pMB->mode = MODE_INTRA;                  pCur->ticks = (int32_t)pCur->stamp % time_base;
631                    pCur->seconds =  ((int32_t)pCur->stamp / time_base)     - ((int32_t)pRef->stamp / time_base) ;
632    
633          if ((pEnc->mbParam.global_flags & XVID_LUMIMASKING) > 0) {                  /* HEAVY DEBUG OUTPUT remove when timecodes prove to be stable */
                 if(pMB->dquant != NO_CHANGE)  
                 {  
                         pMB->mode = MODE_INTRA_Q;  
                         pEnc->mbParam.quant += DQtab[pMB->dquant];  
634    
635                          if (pEnc->mbParam.quant > 31) pEnc->mbParam.quant = 31;  /*              fprintf(stderr,"WriteVop:   %d - %d \n",
636                          if (pEnc->mbParam.quant < 1) pEnc->mbParam.quant = 1;                          ((int32_t)pCur->stamp / time_base), ((int32_t)pRef->stamp / time_base));
637                  }                  fprintf(stderr,"set_timecodes: VOP %1d   stamp=%lld ref_stamp=%lld  base=%d\n",
638          }                          pCur->coding_type, pCur->stamp, pRef->stamp, time_base);
639                    fprintf(stderr,"set_timecodes: VOP %1d   seconds=%d   ticks=%d   (ref-sec=%d  ref-tick=%d)\n",
640                            pCur->coding_type, pCur->seconds, pCur->ticks, pRef->seconds, pRef->ticks);
641    
642          pMB->quant = pEnc->mbParam.quant;  */
643  }  }
644    
645    
646  static int FrameCodeI(Encoder * pEnc, Bitstream * bs, uint32_t *pBits)  
647    /* convert pFrame->intra to coding_type */
648    static int intra2coding_type(int intra)
649  {  {
650            if (intra < 0)  return -1;
651            if (intra == 1) return I_VOP;
652            if (intra == 2) return B_VOP;
653    
654            return P_VOP;
655    }
656    
657    
         DECLARE_ALIGNED_MATRIX(dct_codes, 6, 64, int16_t, CACHE_LINE);  
         DECLARE_ALIGNED_MATRIX(qcoeff,    6, 64, int16_t, CACHE_LINE);  
658    
659    /*****************************************************************************
660     * IPB frame encoder entry point
661     *
662     * Returned values :
663     *      - XVID_ERR_OK    - no errors
664     *      - XVID_ERR_FORMAT - the image subsystem reported the image had a wrong
665     *                                              format
666     ****************************************************************************/
667    
668    int
669    encoder_encode_bframes(Encoder * pEnc,
670                                    XVID_ENC_FRAME * pFrame,
671                                    XVID_ENC_STATS * pResult)
672    {
673          uint16_t x, y;          uint16_t x, y;
674            Bitstream bs;
675            uint32_t bits;
676            int mode = -1; /* Just to shut up compiler warning */
677    
678          pEnc->iFrameNum = 0;          int input_valid = 1;
679          pEnc->mbParam.rounding_type = 1;          int bframes_count = 0;
         pEnc->mbParam.coding_type = I_VOP;  
680    
681          BitstreamWriteVolHeader(bs, &pEnc->mbParam);          ENC_CHECK(pEnc);
682          BitstreamWriteVopHeader(bs, &pEnc->mbParam);          ENC_CHECK(pFrame);
683            ENC_CHECK(pFrame->image);
684    
685          *pBits = BitstreamPos(bs);          start_global_timer();
686    
687          pEnc->sStat.iTextBits = 0;          BitstreamInit(&bs, pFrame->bitstream, 0);
         pEnc->sStat.kblks = pEnc->mbParam.mb_width * pEnc->mbParam.mb_height;  
         pEnc->sStat.mblks = pEnc->sStat.ublks = 0;  
688    
689          for (y = 0; y < pEnc->mbParam.mb_height; y++)  ipvop_loop:
                 for (x = 0; x < pEnc->mbParam.mb_width; x++)  
                 {  
                         MACROBLOCK *pMB = &pEnc->pMBs[x + y * pEnc->mbParam.mb_width];  
690    
691                          CodeIntraMB(pEnc, pMB);          /*
692             * bframe "flush" code
693             */
694    
695                          MBTransQuantIntra(&pEnc->mbParam, pMB, x, y, dct_codes, qcoeff, &pEnc->sCurrent);          if ((pFrame->image == NULL || pEnc->flush_bframes)
696                    && (pEnc->bframenum_head < pEnc->bframenum_tail)) {
697    
698                          start_timer();                  if (pEnc->flush_bframes == 0) {
699                          MBPrediction(&pEnc->mbParam, x, y, pEnc->mbParam.mb_width, qcoeff, pEnc->pMBs);                          /*
700                          stop_prediction_timer();                           * we have reached the end of stream without getting
701                             * a future reference frame... so encode last final
702                             * frame as a pframe
703                             */
704    
705                          start_timer();                          DPRINTF(DPRINTF_DEBUG,"*** BFRAME (final frame) bf: head=%i tail=%i   queue: head=%i tail=%i size=%i",
706                          MBCoding(&pEnc->mbParam, pMB, qcoeff, bs, &pEnc->sStat);                                  pEnc->bframenum_head, pEnc->bframenum_tail,
707                          stop_coding_timer();                                  pEnc->queue_head, pEnc->queue_tail, pEnc->queue_size);
                 }  
708    
709          emms();                          pEnc->bframenum_tail--;
710                            SWAP(FRAMEINFO *, pEnc->current, pEnc->reference);
711    
712          *pBits = BitstreamPos(bs) - *pBits;                          SWAP(FRAMEINFO *, pEnc->current, pEnc->bframes[pEnc->bframenum_tail]);
         pEnc->sStat.fMvPrevSigma = -1;  
         pEnc->sStat.iMvSum = 0;  
         pEnc->sStat.iMvCount = 0;  
         pEnc->mbParam.fixed_code = 2;  
713    
714          return 1;                                        // intra                          FrameCodeP(pEnc, &bs, &bits, 1, 0);
715  }                          bframes_count = 0;
716    
717                            BitstreamPadAlways(&bs);
718                            pFrame->length = BitstreamLength(&bs);
719                            pFrame->intra = 0;
720    
 #define INTRA_THRESHOLD 0.5  
721    
722  static int FrameCodeP(Encoder * pEnc, Bitstream * bs, uint32_t *pBits, bool force_inter, bool vol_header)                          emms();
 {  
         float fSigma;  
723    
724          DECLARE_ALIGNED_MATRIX(dct_codes, 6, 64, int16_t, CACHE_LINE);                          if (pResult) {
725          DECLARE_ALIGNED_MATRIX(qcoeff,    6, 64, int16_t, CACHE_LINE);                                  pResult->quant = pEnc->current->quant;
726                                    pResult->hlength = pFrame->length - (pEnc->current->sStat.iTextBits / 8);
727                                    pResult->kblks = pEnc->current->sStat.kblks;
728                                    pResult->mblks = pEnc->current->sStat.mblks;
729                                    pResult->ublks = pEnc->current->sStat.ublks;
730                            }
731    
732          int iLimit;                          return XVID_ERR_OK;
733          uint32_t x, y;                  }
         int iSearchRange;  
         bool bIntra;  
734    
         IMAGE *pCurrent = &pEnc->sCurrent;  
         IMAGE *pRef = &pEnc->sReference;  
735    
736          start_timer();                  DPRINTF(DPRINTF_DEBUG,"*** BFRAME (flush) bf: head=%i tail=%i   queue: head=%i tail=%i size=%i",
737          image_setedges(pRef,                                  pEnc->bframenum_head, pEnc->bframenum_tail,
738                         pEnc->mbParam.edged_width,                                  pEnc->queue_head, pEnc->queue_tail, pEnc->queue_size);
                        pEnc->mbParam.edged_height,  
                        pEnc->mbParam.width,  
                        pEnc->mbParam.height,  
                        pEnc->mbParam.global_flags & XVID_INTERLACING);  
         stop_edges_timer();  
739    
740          pEnc->mbParam.rounding_type = 1 - pEnc->mbParam.rounding_type;                  FrameCodeB(pEnc, pEnc->bframes[pEnc->bframenum_head], &bs, &bits);
741                    pEnc->bframenum_head++;
742    
743          if (!force_inter)                  BitstreamPadAlways(&bs);
744                  iLimit = (int)(pEnc->mbParam.mb_width * pEnc->mbParam.mb_height * INTRA_THRESHOLD);                  pFrame->length = BitstreamLength(&bs);
745          else                  pFrame->intra = 2;
                 iLimit = pEnc->mbParam.mb_width * pEnc->mbParam.mb_height + 1;  
746    
747          if ((pEnc->mbParam.global_flags & XVID_HALFPEL) > 0) {                  if (pResult) {
748                  start_timer();                          pResult->quant = pEnc->current->quant;
749                  image_interpolate(pRef, &pEnc->vInterH, &pEnc->vInterV, &pEnc->vInterHV,                          pResult->hlength = pFrame->length - (pEnc->current->sStat.iTextBits / 8);
750                                    pEnc->mbParam.edged_width, pEnc->mbParam.edged_height,                          pResult->kblks = pEnc->current->sStat.kblks;
751                                    pEnc->mbParam.rounding_type);                          pResult->mblks = pEnc->current->sStat.mblks;
752                  stop_inter_timer();                          pResult->ublks = pEnc->current->sStat.ublks;
753          }          }
754    
755          start_timer();                  if (input_valid)
756          bIntra = MotionEstimation(pEnc->pMBs, &pEnc->mbParam, &pEnc->sReference,                          queue_image(pEnc, pFrame);
                                   &pEnc->vInterH, &pEnc->vInterV,  
                                   &pEnc->vInterHV, &pEnc->sCurrent, iLimit);  
         stop_motion_timer();  
757    
758          if (bIntra == 1)                  emms();
                 return FrameCodeI(pEnc, bs, pBits);  
759    
760          pEnc->mbParam.coding_type = P_VOP;                  return XVID_ERR_OK;
761            }
762    
763          if(vol_header)          if (pEnc->bframenum_head > 0) {
764                  BitstreamWriteVolHeader(bs, &pEnc->mbParam);                  pEnc->bframenum_head = pEnc->bframenum_tail = 0;
765    
766          BitstreamWriteVopHeader(bs, &pEnc->mbParam);                  /* write an empty marker to the bitstream.
767                       for divx5 decoder compatibility, this marker must consist
768                       of a not-coded p-vop, with a time_base of zero, and time_increment
769                       indentical to the future-referece frame.
770                    */
771    
772          *pBits = BitstreamPos(bs);                  if ((pEnc->mbParam.global & XVID_GLOBAL_PACKED)) {
773                            int tmp;
774    
775          pEnc->sStat.iTextBits = 0;                          DPRINTF(DPRINTF_DEBUG,"*** EMPTY bf: head=%i tail=%i   queue: head=%i tail=%i size=%i",
776          pEnc->sStat.iMvSum = 0;                                  pEnc->bframenum_head, pEnc->bframenum_tail,
777          pEnc->sStat.iMvCount = 0;                                  pEnc->queue_head, pEnc->queue_tail, pEnc->queue_size);
         pEnc->sStat.kblks = pEnc->sStat.mblks = pEnc->sStat.ublks = 0;  
778    
779          for(y = 0; y < pEnc->mbParam.mb_height; y++)                          tmp = pEnc->current->seconds;
780          {                          pEnc->current->seconds = 0; /* force time_base = 0 */
                 for(x = 0; x < pEnc->mbParam.mb_width; x++)  
                 {  
                         MACROBLOCK * pMB = &pEnc->pMBs[x + y * pEnc->mbParam.mb_width];  
781    
782                          bIntra = (pMB->mode == MODE_INTRA) || (pMB->mode == MODE_INTRA_Q);                          BitstreamWriteVopHeader(&bs, &pEnc->mbParam, pEnc->current, 0);
783                            pEnc->current->seconds = tmp;
784    
785                          if (!bIntra)                          BitstreamPadAlways(&bs);
786                          {                          pFrame->length = BitstreamLength(&bs);
787                                  start_timer();                          pFrame->intra = 4;
                                 MBMotionCompensation(pMB,  
                                                      x, y,  
                                                      &pEnc->sReference,  
                                                      &pEnc->vInterH,  
                                                      &pEnc->vInterV,  
                                                      &pEnc->vInterHV,  
                                                      &pEnc->sCurrent,  
                                                      dct_codes,  
                                                      pEnc->mbParam.width,  
                                                      pEnc->mbParam.height,  
                                                      pEnc->mbParam.edged_width,  
                                                      pEnc->mbParam.rounding_type);  
                                 stop_comp_timer();  
788    
789                                  if ((pEnc->mbParam.global_flags & XVID_LUMIMASKING) > 0) {                          if (pResult) {
790                                          if(pMB->dquant != NO_CHANGE) {                                  pResult->quant = pEnc->current->quant;
791                                                  pMB->mode = MODE_INTER_Q;                                  pResult->hlength = pFrame->length - (pEnc->current->sStat.iTextBits / 8);
792                                                  pEnc->mbParam.quant += DQtab[pMB->dquant];                                  pResult->kblks = pEnc->current->sStat.kblks;
793                                                  if (pEnc->mbParam.quant > 31) pEnc->mbParam.quant = 31;                                  pResult->mblks = pEnc->current->sStat.mblks;
794                                                  else if(pEnc->mbParam.quant < 1) pEnc->mbParam.quant = 1;                                  pResult->ublks = pEnc->current->sStat.ublks;
795                                          }                                          }
                                 }  
                                 pMB->quant = pEnc->mbParam.quant;  
796    
797                                  pMB->field_pred = 0;                          if (input_valid)
798                                    queue_image(pEnc, pFrame);
799    
800                            emms();
801    
802                                  pMB->cbp = MBTransQuantInter(&pEnc->mbParam, pMB, x, y, dct_codes, qcoeff, pCurrent);                          return XVID_ERR_OK;
803                          }                          }
                         else  
                         {  
                                 CodeIntraMB(pEnc, pMB);  
                                 MBTransQuantIntra(&pEnc->mbParam, pMB, x, y, dct_codes, qcoeff, pCurrent);  
804                          }                          }
805    
                         start_timer();  
                         MBPrediction(&pEnc->mbParam, x, y, pEnc->mbParam.mb_width, qcoeff, pEnc->pMBs);  
                         stop_prediction_timer();  
806    
807                          if (pMB->mode == MODE_INTRA || pMB->mode == MODE_INTRA_Q)  bvop_loop:
808    
809            if (pEnc->bframenum_dx50bvop != -1)
810                          {                          {
811                                  pEnc->sStat.kblks++;  
812                    SWAP(FRAMEINFO *, pEnc->current, pEnc->reference);
813                    SWAP(FRAMEINFO *, pEnc->current, pEnc->bframes[pEnc->bframenum_dx50bvop]);
814    
815                    if ((pEnc->mbParam.global & XVID_GLOBAL_DEBUG)) {
816                            image_printf(&pEnc->current->image, pEnc->mbParam.edged_width, pEnc->mbParam.height, 5, 100, "DX50 IVOP");
817                          }                          }
818                          else if (pMB->cbp ||  
819                                   pMB->mvs[0].x || pMB->mvs[0].y ||                  if (input_valid)
                                  pMB->mvs[1].x || pMB->mvs[1].y ||  
                                  pMB->mvs[2].x || pMB->mvs[2].y ||  
                                  pMB->mvs[3].x || pMB->mvs[3].y)  
820                          {                          {
821                                  pEnc->sStat.mblks++;                          queue_image(pEnc, pFrame);
822                            input_valid = 0;
823                          }                          }
824                          else  
825                          {          } else if (input_valid) {
826                                  pEnc->sStat.ublks++;  
827                    SWAP(FRAMEINFO *, pEnc->current, pEnc->reference);
828    
829                    start_timer();
830                    if (image_input
831                            (&pEnc->current->image, pEnc->mbParam.width, pEnc->mbParam.height,
832                            pEnc->mbParam.edged_width, pFrame->image, pFrame->stride, pFrame->colorspace, pFrame->general & XVID_INTERLACING))
833                    {
834                            emms();
835                            return XVID_ERR_FORMAT;
836                    }
837                    stop_conv_timer();
838    
839                    if ((pFrame->general & XVID_CHROMAOPT)) {
840                            image_chroma_optimize(&pEnc->current->image,
841                                    pEnc->mbParam.width, pEnc->mbParam.height, pEnc->mbParam.edged_width);
842                    }
843    
844                    /* queue input frame, and dequue next image */
845                    if (pEnc->queue_size > 0)
846                    {
847                            image_swap(&pEnc->current->image, &pEnc->queue[pEnc->queue_tail]);
848                            if (pEnc->queue_head != pEnc->queue_tail)
849                            {
850                                    image_swap(&pEnc->current->image, &pEnc->queue[pEnc->queue_head]);
851                            }
852                            pEnc->queue_head =  (pEnc->queue_head + 1) % pEnc->mbParam.max_bframes;
853                            pEnc->queue_tail =  (pEnc->queue_tail + 1) % pEnc->mbParam.max_bframes;
854                    }
855    
856            } else if (pEnc->queue_size > 0) {
857    
858                    SWAP(FRAMEINFO *, pEnc->current, pEnc->reference);
859    
860                    image_swap(&pEnc->current->image, &pEnc->queue[pEnc->queue_head]);
861                    pEnc->queue_head =  (pEnc->queue_head + 1) % pEnc->mbParam.max_bframes;
862                    pEnc->queue_size--;
863    
864            } else {
865    
866                    /* if nothing was encoded, write an 'ignore this frame' flag
867                       to the bitstream */
868    
869                    if (BitstreamPos(&bs) == 0) {
870    
871                            DPRINTF(DPRINTF_DEBUG,"*** SKIP bf: head=%i tail=%i   queue: head=%i tail=%i size=%i",
872                                    pEnc->bframenum_head, pEnc->bframenum_tail,
873                                    pEnc->queue_head, pEnc->queue_tail, pEnc->queue_size);
874    
875                            /* That disabled line of code was supposed to inform VirtualDub
876                             * that the frame was a dummy delay frame - now disabled (thx god :-)
877                             */
878                            //BitstreamPutBits(&bs, 0x7f, 8);
879                            pFrame->intra = 5;
880    
881                            if (pResult) {
882                                    /*
883                                     * We must decide what to put there because i know some apps
884                                     * are storing statistics about quantizers and just do
885                                     * stats[quant]++ or stats[quant-1]++
886                                     * transcode is one of these app with its 2pass module
887                                     */
888    
889                                    /*
890                                     * For now i prefer 31 than 0 that could lead to a segfault
891                                     * in transcode
892                                     */
893                                    pResult->quant = 31;
894    
895                                    pResult->hlength = 0;
896                                    pResult->kblks = 0;
897                                    pResult->mblks = 0;
898                                    pResult->ublks = 0;
899                            }
900                    } else {
901    
902                            if (pResult) {
903                                    pResult->quant = pEnc->current->quant;
904                                    pResult->hlength = pFrame->length - (pEnc->current->sStat.iTextBits / 8);
905                                    pResult->kblks = pEnc->current->sStat.kblks;
906                                    pResult->mblks = pEnc->current->sStat.mblks;
907                                    pResult->ublks = pEnc->current->sStat.ublks;
908                            }
909    
910                    }
911    
912                    pFrame->length = BitstreamLength(&bs);
913    
914                    emms();
915    
916                    return XVID_ERR_OK;
917            }
918    
919            pEnc->flush_bframes = 0;
920    
921            emms();
922    
923            /* only inc frame num, adapt quant, etc. if we havent seen it before */
924            if (pEnc->bframenum_dx50bvop < 0 )
925            {
926                    mode = intra2coding_type(pFrame->intra);
927                    if (pFrame->quant == 0)
928                            pEnc->current->quant = RateControlGetQ(&pEnc->rate_control, 0);
929                    else
930                            pEnc->current->quant = pFrame->quant;
931    
932    /*              if (pEnc->current->quant < 1)
933                            pEnc->current->quant = 1;
934    
935                    if (pEnc->current->quant > 31)
936                            pEnc->current->quant = 31;
937    */
938                    pEnc->current->global_flags = pFrame->general;
939                    pEnc->current->motion_flags = pFrame->motion;
940    
941                    /* ToDo : dynamic fcode (in both directions) */
942                    pEnc->current->fcode = pEnc->mbParam.m_fcode;
943                    pEnc->current->bcode = pEnc->mbParam.m_fcode;
944    
945                    inc_frame_num(pEnc);
946    
947                    if (pFrame->general & XVID_EXTRASTATS)
948                    {       image_copy(&pEnc->sOriginal, &pEnc->current->image,
949                                       pEnc->mbParam.edged_width, pEnc->mbParam.height);
950                    }
951    
952                    emms();
953    
954                    if ((pEnc->mbParam.global & XVID_GLOBAL_DEBUG)) {
955                            image_printf(&pEnc->current->image, pEnc->mbParam.edged_width, pEnc->mbParam.height, 5, 5,
956                                    "%i  if:%i  st:%i", pEnc->m_framenum++, pEnc->iFrameNum, pEnc->current->stamp);
957                    }
958    
959            /* %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
960             * Luminance masking
961             * %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% */
962    
963                    if ((pEnc->current->global_flags & XVID_LUMIMASKING)) {
964                            int *temp_dquants =
965                                    (int *) xvid_malloc(pEnc->mbParam.mb_width *
966                                                                    pEnc->mbParam.mb_height * sizeof(int),
967                                                                    CACHE_LINE);
968    
969                            pEnc->current->quant =
970                                    adaptive_quantization(pEnc->current->image.y,
971                                                                      pEnc->mbParam.edged_width, temp_dquants,
972                                                                      pEnc->current->quant, pEnc->current->quant,
973                                                                      2 * pEnc->current->quant,
974                                                                      pEnc->mbParam.mb_width,
975                                                                      pEnc->mbParam.mb_height);
976    
977                            for (y = 0; y < pEnc->mbParam.mb_height; y++) {
978    
979    #define OFFSET(x,y) ((x) + (y)*pEnc->mbParam.mb_width)
980    
981                                    for (x = 0; x < pEnc->mbParam.mb_width; x++) {
982                                            MACROBLOCK *pMB = &pEnc->current->mbs[OFFSET(x, y)];
983    
984                                            pMB->dquant = iDQtab[temp_dquants[OFFSET(x, y)] + 2];
985                                    }
986    #undef OFFSET
987                            }
988                            xvid_free(temp_dquants);
989                    }
990            }
991    
992            /* %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
993             * ivop/pvop/bvop selection
994             * %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% */
995            pEnc->iFrameNum++;
996    
997            if (pEnc->iFrameNum == 0 || pEnc->bframenum_dx50bvop >= 0 ||
998                    (mode < 0 && pEnc->mbParam.iMaxKeyInterval > 0 &&
999                            pEnc->iFrameNum >= pEnc->mbParam.iMaxKeyInterval)) {
1000                    mode = I_VOP;
1001            } else {
1002                    mode = MEanalysis(&pEnc->reference->image, pEnc->current,
1003                                            &pEnc->mbParam, pEnc->mbParam.iMaxKeyInterval,
1004                                            (/*mode < 0*/1/*hack*/) ? pEnc->iFrameNum : 0,
1005                                            bframes_count++, pFrame->bframe_threshold);
1006            }
1007    
1008            if (mode == I_VOP) {
1009                    /*
1010                     * This will be coded as an Intra Frame
1011                     */
1012                    if ((pEnc->current->global_flags & XVID_QUARTERPEL))
1013                            pEnc->mbParam.m_quarterpel = 1;
1014                    else
1015                            pEnc->mbParam.m_quarterpel = 0;
1016    
1017                    if (pEnc->current->global_flags & XVID_MPEGQUANT) pEnc->mbParam.m_quant_type = MPEG4_QUANT;
1018    
1019                    if ((pEnc->current->global_flags & XVID_CUSTOM_QMATRIX) > 0) {
1020                            if (pFrame->quant_intra_matrix != NULL)
1021                                    set_intra_matrix(pFrame->quant_intra_matrix);
1022                            if (pFrame->quant_inter_matrix != NULL)
1023                                    set_inter_matrix(pFrame->quant_inter_matrix);
1024                    }
1025    
1026    
1027                    DPRINTF(DPRINTF_DEBUG,"*** IFRAME bf: head=%i tail=%i   queue: head=%i tail=%i size=%i",
1028                                    pEnc->bframenum_head, pEnc->bframenum_tail,
1029                                    pEnc->queue_head, pEnc->queue_tail, pEnc->queue_size);
1030    
1031                    if ((pEnc->mbParam.global & XVID_GLOBAL_DEBUG)) {
1032                            image_printf(&pEnc->current->image, pEnc->mbParam.edged_width, pEnc->mbParam.height, 5, 200, "IVOP");
1033                    }
1034    
1035                    /* when we reach an iframe in DX50BVOP mode, encode the last bframe as a pframe */
1036    
1037                    if ((pEnc->mbParam.global & XVID_GLOBAL_DX50BVOP) && pEnc->bframenum_tail > 0) {
1038    
1039                            pEnc->bframenum_tail--;
1040                            pEnc->bframenum_dx50bvop = pEnc->bframenum_tail;
1041    
1042                            SWAP(FRAMEINFO *, pEnc->current, pEnc->bframes[pEnc->bframenum_dx50bvop]);
1043                            if ((pEnc->mbParam.global & XVID_GLOBAL_DEBUG)) {
1044                                    image_printf(&pEnc->current->image, pEnc->mbParam.edged_width, pEnc->mbParam.height, 5, 100, "DX50 BVOP->PVOP");
1045                            }
1046                            FrameCodeP(pEnc, &bs, &bits, 1, 0);
1047                            bframes_count = 0;
1048                            pFrame->intra = 0;
1049    
1050                    } else {
1051    
1052                            FrameCodeI(pEnc, &bs, &bits);
1053                            bframes_count = 0;
1054                            pFrame->intra = 1;
1055    
1056                            pEnc->bframenum_dx50bvop = -1;
1057                    }
1058    
1059                    pEnc->flush_bframes = 1;
1060    
1061                    if ((pEnc->mbParam.global & XVID_GLOBAL_PACKED) && pEnc->bframenum_tail > 0) {
1062                            BitstreamPadAlways(&bs);
1063                            input_valid = 0;
1064                            goto ipvop_loop;
1065                    }
1066    
1067                    /*
1068                     * NB : sequences like "IIBB" decode fine with msfdam but,
1069                     *        go screwy with divx 5.00
1070                     */
1071            } else if (mode == P_VOP || mode == S_VOP || pEnc->bframenum_tail >= pEnc->mbParam.max_bframes) {
1072                    /*
1073                     * This will be coded as a Predicted Frame
1074                     */
1075    
1076                    DPRINTF(DPRINTF_DEBUG,"*** PFRAME bf: head=%i tail=%i   queue: head=%i tail=%i size=%i",
1077                                    pEnc->bframenum_head, pEnc->bframenum_tail,
1078                                    pEnc->queue_head, pEnc->queue_tail, pEnc->queue_size);
1079    
1080                    if ((pEnc->mbParam.global & XVID_GLOBAL_DEBUG)) {
1081                            image_printf(&pEnc->current->image, pEnc->mbParam.edged_width, pEnc->mbParam.height, 5, 200, "PVOP");
1082                    }
1083    
1084                    FrameCodeP(pEnc, &bs, &bits, 1, 0);
1085                    bframes_count = 0;
1086                    pFrame->intra = 0;
1087                    pEnc->flush_bframes = 1;
1088    
1089                    if ((pEnc->mbParam.global & XVID_GLOBAL_PACKED) && (pEnc->bframenum_tail > 0)) {
1090                            BitstreamPadAlways(&bs);
1091                            input_valid = 0;
1092                            goto ipvop_loop;
1093                    }
1094    
1095            } else {        /* mode == B_VOP */
1096                    /*
1097                     * This will be coded as a Bidirectional Frame
1098                     */
1099    
1100                    if ((pEnc->mbParam.global & XVID_GLOBAL_DEBUG)) {
1101                            image_printf(&pEnc->current->image, pEnc->mbParam.edged_width, pEnc->mbParam.height, 5, 200, "BVOP");
1102                    }
1103    
1104                    if (pFrame->bquant < 1) {
1105                            pEnc->current->quant = ((((pEnc->reference->quant + pEnc->current->quant) *
1106                                    pEnc->mbParam.bquant_ratio) / 2) + pEnc->mbParam.bquant_offset)/100;
1107    
1108                    } else {
1109                            pEnc->current->quant = pFrame->bquant;
1110                    }
1111    
1112                    if (pEnc->current->quant < 1)
1113                            pEnc->current->quant = 1;
1114                    else if (pEnc->current->quant > 31)
1115                            pEnc->current->quant = 31;
1116    
1117                    DPRINTF(DPRINTF_DEBUG,"*** BFRAME (store) bf: head=%i tail=%i   queue: head=%i tail=%i size=%i  quant=%i\n",
1118                                    pEnc->bframenum_head, pEnc->bframenum_tail,
1119                                    pEnc->queue_head, pEnc->queue_tail, pEnc->queue_size,pEnc->current->quant);
1120    
1121                    /* store frame into bframe buffer & swap ref back to current */
1122                    SWAP(FRAMEINFO *, pEnc->current, pEnc->bframes[pEnc->bframenum_tail]);
1123                    SWAP(FRAMEINFO *, pEnc->current, pEnc->reference);
1124    
1125                    pEnc->bframenum_tail++;
1126    
1127                    /* bframe report by koepi */
1128                    pFrame->intra = 2;
1129                    pFrame->length = 0;
1130    
1131                    input_valid = 0;
1132                    goto bvop_loop;
1133            }
1134    
1135            BitstreamPadAlways(&bs);
1136            pFrame->length = BitstreamLength(&bs);
1137    
1138            if (pResult) {
1139                    pResult->quant = pEnc->current->quant;
1140                    pResult->hlength = pFrame->length - (pEnc->current->sStat.iTextBits / 8);
1141                    pResult->kblks = pEnc->current->sStat.kblks;
1142                    pResult->mblks = pEnc->current->sStat.mblks;
1143                    pResult->ublks = pEnc->current->sStat.ublks;
1144    
1145                    if (pFrame->general & XVID_EXTRASTATS)
1146                    {       pResult->sse_y =
1147                                    plane_sse( pEnc->sOriginal.y, pEnc->current->image.y,
1148                                                       pEnc->mbParam.edged_width, pEnc->mbParam.width,
1149                                                       pEnc->mbParam.height);
1150    
1151                            pResult->sse_u =
1152                                    plane_sse( pEnc->sOriginal.u, pEnc->current->image.u,
1153                                                       pEnc->mbParam.edged_width/2, pEnc->mbParam.width/2,
1154                                                       pEnc->mbParam.height/2);
1155    
1156                            pResult->sse_v =
1157                                    plane_sse( pEnc->sOriginal.v, pEnc->current->image.v,
1158                                                       pEnc->mbParam.edged_width/2, pEnc->mbParam.width/2,
1159                                                       pEnc->mbParam.height/2);
1160                    }
1161            }
1162    
1163            emms();
1164    
1165            if (pFrame->quant == 0) {
1166                    RateControlUpdate(&pEnc->rate_control, pEnc->current->quant,
1167                                                      pFrame->length, pFrame->intra);
1168            }
1169    
1170            stop_global_timer();
1171            write_timer();
1172    
1173            emms();
1174            return XVID_ERR_OK;
1175    }
1176    
1177    
1178    
1179    /*****************************************************************************
1180     * "original" IP frame encoder entry point
1181     *
1182     * Returned values :
1183     *      - XVID_ERR_OK    - no errors
1184     *      - XVID_ERR_FORMAT - the image subsystem reported the image had a wrong
1185     *                                              format
1186     ****************************************************************************/
1187    
1188    int
1189    encoder_encode(Encoder * pEnc,
1190                               XVID_ENC_FRAME * pFrame,
1191                               XVID_ENC_STATS * pResult)
1192    {
1193            uint16_t x, y;
1194            Bitstream bs;
1195            uint32_t bits;
1196            uint16_t write_vol_header = 0;
1197    
1198            float psnr;
1199            char temp[128];
1200    
1201            start_global_timer();
1202    
1203            ENC_CHECK(pEnc);
1204            ENC_CHECK(pFrame);
1205            ENC_CHECK(pFrame->bitstream);
1206            ENC_CHECK(pFrame->image);
1207    
1208            SWAP(FRAMEINFO *, pEnc->current, pEnc->reference);
1209    
1210            pEnc->current->global_flags = pFrame->general;
1211            pEnc->current->motion_flags = pFrame->motion;
1212            pEnc->mbParam.hint = &pFrame->hint;
1213    
1214            inc_frame_num(pEnc);
1215    
1216            /* disable alternate scan flag if interlacing is not enabled */
1217            if ((pEnc->current->global_flags & XVID_ALTERNATESCAN) &&
1218                    !(pEnc->current->global_flags & XVID_INTERLACING))
1219            {
1220                    pEnc->current->global_flags -= XVID_ALTERNATESCAN;
1221            }
1222    
1223            start_timer();
1224            if (image_input
1225                    (&pEnc->current->image, pEnc->mbParam.width, pEnc->mbParam.height,
1226                     pEnc->mbParam.edged_width, pFrame->image, pFrame->stride, pFrame->colorspace, pFrame->general & XVID_INTERLACING) < 0)
1227                    return XVID_ERR_FORMAT;
1228            stop_conv_timer();
1229    
1230            if ((pFrame->general & XVID_CHROMAOPT)) {
1231                    image_chroma_optimize(&pEnc->current->image,
1232                            pEnc->mbParam.width, pEnc->mbParam.height, pEnc->mbParam.edged_width);
1233            }
1234    
1235            if (pFrame->general & XVID_EXTRASTATS)
1236            {       image_copy(&pEnc->sOriginal, &pEnc->current->image,
1237                                       pEnc->mbParam.edged_width, pEnc->mbParam.height);
1238            }
1239    
1240            emms();
1241    
1242            BitstreamInit(&bs, pFrame->bitstream, 0);
1243    
1244            if (pFrame->quant == 0) {
1245                    pEnc->current->quant = RateControlGetQ(&pEnc->rate_control, 0);
1246            } else {
1247                    pEnc->current->quant = pFrame->quant;
1248            }
1249    
1250            if ((pEnc->current->global_flags & XVID_QUARTERPEL))
1251                    pEnc->mbParam.m_quarterpel = 1;
1252            else
1253                    pEnc->mbParam.m_quarterpel = 0;
1254    
1255            if ((pEnc->current->global_flags & XVID_LUMIMASKING)) {
1256                    int *temp_dquants =
1257                            (int *) xvid_malloc(pEnc->mbParam.mb_width *
1258                                                                    pEnc->mbParam.mb_height * sizeof(int),
1259                                                                    CACHE_LINE);
1260    
1261                    pEnc->current->quant =
1262                            adaptive_quantization(pEnc->current->image.y,
1263                                                                      pEnc->mbParam.edged_width, temp_dquants,
1264                                                                      pEnc->current->quant, pEnc->current->quant,
1265                                                                      2 * pEnc->current->quant,
1266                                                                      pEnc->mbParam.mb_width,
1267                                                                      pEnc->mbParam.mb_height);
1268    
1269                    for (y = 0; y < pEnc->mbParam.mb_height; y++) {
1270    
1271    #define OFFSET(x,y) ((x) + (y)*pEnc->mbParam.mb_width)
1272    
1273                            for (x = 0; x < pEnc->mbParam.mb_width; x++) {
1274    
1275    
1276                                    MACROBLOCK *pMB = &pEnc->current->mbs[OFFSET(x, y)];
1277    
1278                                    pMB->dquant = iDQtab[temp_dquants[OFFSET(x, y)] + 2];
1279                            }
1280    
1281    #undef OFFSET
1282                    }
1283    
1284                    xvid_free(temp_dquants);
1285            }
1286    
1287            if (pEnc->current->global_flags & XVID_H263QUANT) {
1288                    if (pEnc->mbParam.m_quant_type != H263_QUANT)
1289                            write_vol_header = 1;
1290                    pEnc->mbParam.m_quant_type = H263_QUANT;
1291            } else if (pEnc->current->global_flags & XVID_MPEGQUANT) {
1292                    int matrix1_changed, matrix2_changed;
1293    
1294                    matrix1_changed = matrix2_changed = 0;
1295    
1296                    if (pEnc->mbParam.m_quant_type != MPEG4_QUANT)
1297                            write_vol_header = 1;
1298    
1299                    pEnc->mbParam.m_quant_type = MPEG4_QUANT;
1300    
1301                    if ((pEnc->current->global_flags & XVID_CUSTOM_QMATRIX) > 0) {
1302                            if (pFrame->quant_intra_matrix != NULL)
1303                                    matrix1_changed = set_intra_matrix(pFrame->quant_intra_matrix);
1304                            if (pFrame->quant_inter_matrix != NULL)
1305                                    matrix2_changed = set_inter_matrix(pFrame->quant_inter_matrix);
1306                    } else {
1307                            matrix1_changed = set_intra_matrix(get_default_intra_matrix());
1308                            matrix2_changed = set_inter_matrix(get_default_inter_matrix());
1309                    }
1310                    if (write_vol_header == 0)
1311                            write_vol_header = matrix1_changed | matrix2_changed;
1312            }
1313    
1314            if (pFrame->intra < 0) {
1315                    if ((pEnc->iFrameNum == -1)
1316                            || ((pEnc->mbParam.iMaxKeyInterval > 0)
1317                                    && (pEnc->iFrameNum >= pEnc->mbParam.iMaxKeyInterval))) {
1318                            pFrame->intra = FrameCodeI(pEnc, &bs, &bits);
1319                    } else {
1320                            pFrame->intra = FrameCodeP(pEnc, &bs, &bits, 0, write_vol_header);
1321                    }
1322            } else {
1323                    if (pFrame->intra == 1) {
1324                            pFrame->intra = FrameCodeI(pEnc, &bs, &bits);
1325                    } else {
1326                            pFrame->intra = FrameCodeP(pEnc, &bs, &bits, 1, write_vol_header);
1327                    }
1328    
1329            }
1330    
1331            /* Relic from OpenDivX - now disabled
1332            BitstreamPutBits(&bs, 0xFFFF, 16);
1333            BitstreamPutBits(&bs, 0xFFFF, 16);
1334            */
1335    
1336            BitstreamPadAlways(&bs);
1337            pFrame->length = BitstreamLength(&bs);
1338    
1339            if (pResult) {
1340                    pResult->quant = pEnc->current->quant;
1341                    pResult->hlength = pFrame->length - (pEnc->current->sStat.iTextBits / 8);
1342                    pResult->kblks = pEnc->current->sStat.kblks;
1343                    pResult->mblks = pEnc->current->sStat.mblks;
1344                    pResult->ublks = pEnc->current->sStat.ublks;
1345            }
1346    
1347            emms();
1348    
1349            if (pFrame->quant == 0) {
1350                    RateControlUpdate(&pEnc->rate_control, pEnc->current->quant,
1351                                                      pFrame->length, pFrame->intra);
1352            }
1353            if (pFrame->general & XVID_EXTRASTATS)
1354            {
1355                    psnr =
1356                            image_psnr(&pEnc->sOriginal, &pEnc->current->image,
1357                                               pEnc->mbParam.edged_width, pEnc->mbParam.width,
1358                                               pEnc->mbParam.height);
1359    
1360                    snprintf(temp, 127, "PSNR: %f\n", psnr);
1361            }
1362    
1363            pEnc->iFrameNum++;
1364    
1365            stop_global_timer();
1366            write_timer();
1367    
1368            return XVID_ERR_OK;
1369    }
1370    
1371    
1372    static __inline void
1373    CodeIntraMB(Encoder * pEnc,
1374                            MACROBLOCK * pMB)
1375    {
1376    
1377            pMB->mode = MODE_INTRA;
1378    
1379            /* zero mv statistics */
1380            pMB->mvs[0].x = pMB->mvs[1].x = pMB->mvs[2].x = pMB->mvs[3].x = 0;
1381            pMB->mvs[0].y = pMB->mvs[1].y = pMB->mvs[2].y = pMB->mvs[3].y = 0;
1382            pMB->sad8[0] = pMB->sad8[1] = pMB->sad8[2] = pMB->sad8[3] = 0;
1383            pMB->sad16 = 0;
1384    
1385            if ((pEnc->current->global_flags & XVID_LUMIMASKING)) {
1386                    if (pMB->dquant != NO_CHANGE) {
1387                            pMB->mode = MODE_INTRA_Q;
1388                            pEnc->current->quant += DQtab[pMB->dquant];
1389    
1390                            if (pEnc->current->quant > 31)
1391                                    pEnc->current->quant = 31;
1392                            if (pEnc->current->quant < 1)
1393                                    pEnc->current->quant = 1;
1394                    }
1395            }
1396    
1397            pMB->quant = pEnc->current->quant;
1398    }
1399    
1400    
1401    #define FCODEBITS       3
1402    #define MODEBITS        5
1403    
1404    void
1405    HintedMESet(Encoder * pEnc,
1406                            int *intra)
1407    {
1408            HINTINFO *hint;
1409            Bitstream bs;
1410            int length, high;
1411            uint32_t x, y;
1412    
1413            hint = pEnc->mbParam.hint;
1414    
1415            if (hint->rawhints) {
1416                    *intra = hint->mvhint.intra;
1417            } else {
1418                    BitstreamInit(&bs, hint->hintstream, hint->hintlength);
1419                    *intra = BitstreamGetBit(&bs);
1420            }
1421    
1422            if (*intra) {
1423                    return;
1424            }
1425    
1426            pEnc->current->fcode = (hint->rawhints) ?
1427                    (uint32_t)hint->mvhint.fcode : BitstreamGetBits(&bs, FCODEBITS);
1428    
1429            length = pEnc->current->fcode + 5;
1430            high = 1 << (length - 1);
1431    
1432            for (y = 0; y < pEnc->mbParam.mb_height; ++y) {
1433                    for (x = 0; x < pEnc->mbParam.mb_width; ++x) {
1434                            MACROBLOCK *pMB =
1435                                    &pEnc->current->mbs[x + y * pEnc->mbParam.mb_width];
1436                            MVBLOCKHINT *bhint =
1437                                    &hint->mvhint.block[x + y * pEnc->mbParam.mb_width];
1438                            VECTOR pred;
1439                            VECTOR tmp;
1440                            int vec;
1441    
1442                            pMB->mode =     (hint->rawhints) ?
1443                                    (uint32_t)bhint->mode : BitstreamGetBits(&bs, MODEBITS);
1444    
1445                            pMB->mode = (pMB->mode == MODE_INTER_Q) ? MODE_INTER : pMB->mode;
1446                            pMB->mode = (pMB->mode == MODE_INTRA_Q) ? MODE_INTRA : pMB->mode;
1447    
1448                            if (pMB->mode == MODE_INTER) {
1449                                    tmp.x = (hint->rawhints) ?
1450                                            bhint->mvs[0].x : (int)BitstreamGetBits(&bs, length);
1451                                    tmp.y = (hint->rawhints) ?
1452                                            bhint->mvs[0].y : (int)BitstreamGetBits(&bs, length);
1453                                    tmp.x -= (tmp.x >= high) ? high * 2 : 0;
1454                                    tmp.y -= (tmp.y >= high) ? high * 2 : 0;
1455    
1456                                    pred = get_pmv2(pEnc->current->mbs,pEnc->mbParam.mb_width,0,x,y,0);
1457    
1458                                    for (vec = 0; vec < 4; ++vec) {
1459                                            pMB->mvs[vec].x = tmp.x;
1460                                            pMB->mvs[vec].y = tmp.y;
1461                                            pMB->pmvs[vec].x = pMB->mvs[0].x - pred.x;
1462                                            pMB->pmvs[vec].y = pMB->mvs[0].y - pred.y;
1463                                    }
1464                            } else if (pMB->mode == MODE_INTER4V) {
1465                                    for (vec = 0; vec < 4; ++vec) {
1466                                            tmp.x = (hint->rawhints) ?
1467                                                    bhint->mvs[vec].x : (int)BitstreamGetBits(&bs, length);
1468                                            tmp.y = (hint->rawhints) ?
1469                                                    bhint->mvs[vec].y : (int)BitstreamGetBits(&bs, length);
1470                                            tmp.x -= (tmp.x >= high) ? high * 2 : 0;
1471                                            tmp.y -= (tmp.y >= high) ? high * 2 : 0;
1472    
1473                                            pred = get_pmv2(pEnc->current->mbs,pEnc->mbParam.mb_width,0,x,y,vec);
1474    
1475                                            pMB->mvs[vec].x = tmp.x;
1476                                            pMB->mvs[vec].y = tmp.y;
1477                                            pMB->pmvs[vec].x = pMB->mvs[vec].x - pred.x;
1478                                            pMB->pmvs[vec].y = pMB->mvs[vec].y - pred.y;
1479                                    }
1480                            } else                          /* intra / stuffing / not_coded */
1481                            {
1482                                    for (vec = 0; vec < 4; ++vec) {
1483                                            pMB->mvs[vec].x = pMB->mvs[vec].y = 0;
1484                                    }
1485                            }
1486    
1487                            if (pMB->mode == MODE_INTER4V &&
1488                                    (pEnc->current->global_flags & XVID_LUMIMASKING)
1489                                    && pMB->dquant != NO_CHANGE) {
1490                                    pMB->mode = MODE_INTRA;
1491    
1492                                    for (vec = 0; vec < 4; ++vec) {
1493                                            pMB->mvs[vec].x = pMB->mvs[vec].y = 0;
1494                                    }
1495                            }
1496                    }
1497            }
1498    }
1499    
1500    
1501    void
1502    HintedMEGet(Encoder * pEnc,
1503                            int intra)
1504    {
1505            HINTINFO *hint;
1506            Bitstream bs;
1507            uint32_t x, y;
1508            int length, high;
1509    
1510            hint = pEnc->mbParam.hint;
1511    
1512            if (hint->rawhints) {
1513                    hint->mvhint.intra = intra;
1514            } else {
1515                    BitstreamInit(&bs, hint->hintstream, 0);
1516                    BitstreamPutBit(&bs, intra);
1517            }
1518    
1519            if (intra) {
1520                    if (!hint->rawhints) {
1521                            BitstreamPadAlways(&bs);
1522                            hint->hintlength = BitstreamLength(&bs);
1523                    }
1524                    return;
1525            }
1526    
1527            length = pEnc->current->fcode + 5;
1528            high = 1 << (length - 1);
1529    
1530            if (hint->rawhints) {
1531                    hint->mvhint.fcode = pEnc->current->fcode;
1532            } else {
1533                    BitstreamPutBits(&bs, pEnc->current->fcode, FCODEBITS);
1534            }
1535    
1536            for (y = 0; y < pEnc->mbParam.mb_height; ++y) {
1537                    for (x = 0; x < pEnc->mbParam.mb_width; ++x) {
1538                            MACROBLOCK *pMB =
1539                                    &pEnc->current->mbs[x + y * pEnc->mbParam.mb_width];
1540                            MVBLOCKHINT *bhint =
1541                                    &hint->mvhint.block[x + y * pEnc->mbParam.mb_width];
1542                            VECTOR tmp;
1543    
1544                            if (hint->rawhints) {
1545                                    bhint->mode = pMB->mode;
1546                            } else {
1547                                    BitstreamPutBits(&bs, pMB->mode, MODEBITS);
1548                            }
1549    
1550                            if (pMB->mode == MODE_INTER || pMB->mode == MODE_INTER_Q) {
1551                                    tmp.x = pMB->mvs[0].x;
1552                                    tmp.y = pMB->mvs[0].y;
1553                                    tmp.x += (tmp.x < 0) ? high * 2 : 0;
1554                                    tmp.y += (tmp.y < 0) ? high * 2 : 0;
1555    
1556                                    if (hint->rawhints) {
1557                                            bhint->mvs[0].x = tmp.x;
1558                                            bhint->mvs[0].y = tmp.y;
1559                                    } else {
1560                                            BitstreamPutBits(&bs, tmp.x, length);
1561                                            BitstreamPutBits(&bs, tmp.y, length);
1562                                    }
1563                            } else if (pMB->mode == MODE_INTER4V) {
1564                                    int vec;
1565    
1566                                    for (vec = 0; vec < 4; ++vec) {
1567                                            tmp.x = pMB->mvs[vec].x;
1568                                            tmp.y = pMB->mvs[vec].y;
1569                                            tmp.x += (tmp.x < 0) ? high * 2 : 0;
1570                                            tmp.y += (tmp.y < 0) ? high * 2 : 0;
1571    
1572                                            if (hint->rawhints) {
1573                                                    bhint->mvs[vec].x = tmp.x;
1574                                                    bhint->mvs[vec].y = tmp.y;
1575                                            } else {
1576                                                    BitstreamPutBits(&bs, tmp.x, length);
1577                                                    BitstreamPutBits(&bs, tmp.y, length);
1578                                            }
1579                                    }
1580                            }
1581                    }
1582            }
1583    
1584            if (!hint->rawhints) {
1585                    BitstreamPad(&bs);
1586                    hint->hintlength = BitstreamLength(&bs);
1587            }
1588    }
1589    
1590    
1591    static int
1592    FrameCodeI(Encoder * pEnc,
1593                       Bitstream * bs,
1594                       uint32_t * pBits)
1595    {
1596            int mb_width = pEnc->mbParam.mb_width;
1597            int mb_height = pEnc->mbParam.mb_height;
1598    
1599            DECLARE_ALIGNED_MATRIX(dct_codes, 6, 64, int16_t, CACHE_LINE);
1600            DECLARE_ALIGNED_MATRIX(qcoeff, 6, 64, int16_t, CACHE_LINE);
1601    
1602            uint16_t x, y;
1603    
1604            if ((pEnc->current->global_flags & XVID_REDUCED))
1605            {
1606                    mb_width = (pEnc->mbParam.width + 31) / 32;
1607                    mb_height = (pEnc->mbParam.height + 31) / 32;
1608    
1609                    /* 16x16->8x8 downsample requires 1 additional edge pixel*/
1610                    /* XXX: setedges is overkill */
1611                    start_timer();
1612                    image_setedges(&pEnc->current->image,
1613                            pEnc->mbParam.edged_width, pEnc->mbParam.edged_height,
1614                            pEnc->mbParam.width, pEnc->mbParam.height);
1615                    stop_edges_timer();
1616            }
1617            pEnc->iFrameNum = 0;
1618            pEnc->mbParam.m_rounding_type = 1;
1619            pEnc->current->rounding_type = pEnc->mbParam.m_rounding_type;
1620            pEnc->current->quarterpel =  pEnc->mbParam.m_quarterpel;
1621            pEnc->current->coding_type = I_VOP;
1622    
1623            BitstreamWriteVolHeader(bs, &pEnc->mbParam, pEnc->current);
1624    
1625            set_timecodes(pEnc->current,pEnc->reference,pEnc->mbParam.fbase);
1626    
1627            BitstreamPadAlways(bs);
1628            BitstreamWriteVopHeader(bs, &pEnc->mbParam, pEnc->current, 1);
1629    
1630            *pBits = BitstreamPos(bs);
1631    
1632            pEnc->current->sStat.iTextBits = 0;
1633            pEnc->current->sStat.kblks = mb_width * mb_height;
1634            pEnc->current->sStat.mblks = pEnc->current->sStat.ublks = 0;
1635    
1636            for (y = 0; y < mb_height; y++)
1637                    for (x = 0; x < mb_width; x++) {
1638                            MACROBLOCK *pMB =
1639                                    &pEnc->current->mbs[x + y * pEnc->mbParam.mb_width];
1640    
1641                            CodeIntraMB(pEnc, pMB);
1642    
1643                            MBTransQuantIntra(&pEnc->mbParam, pEnc->current, pMB, x, y,
1644                                                              dct_codes, qcoeff);
1645    
1646                            start_timer();
1647                            MBPrediction(pEnc->current, x, y, pEnc->mbParam.mb_width, qcoeff);
1648                            stop_prediction_timer();
1649    
1650                            start_timer();
1651                            if (pEnc->current->global_flags & XVID_GREYSCALE)
1652                            {       pMB->cbp &= 0x3C;               /* keep only bits 5-2 */
1653                                    qcoeff[4*64+0]=0;               /* zero, because for INTRA MBs DC value is saved */
1654                                    qcoeff[5*64+0]=0;
1655                            }
1656                            MBCoding(pEnc->current, pMB, qcoeff, bs, &pEnc->current->sStat);
1657                            stop_coding_timer();
1658                    }
1659    
1660            if ((pEnc->current->global_flags & XVID_REDUCED))
1661            {
1662                    image_deblock_rrv(&pEnc->current->image, pEnc->mbParam.edged_width,
1663                            pEnc->current->mbs, mb_width, mb_height, pEnc->mbParam.mb_width,
1664                            16, XVID_DEC_DEBLOCKY|XVID_DEC_DEBLOCKUV);
1665            }
1666            emms();
1667    
1668            *pBits = BitstreamPos(bs) - *pBits;
1669            pEnc->fMvPrevSigma = -1;
1670            pEnc->mbParam.m_fcode = 2;
1671    
1672            if (pEnc->current->global_flags & XVID_HINTEDME_GET) {
1673                    HintedMEGet(pEnc, 1);
1674            }
1675    
1676            return 1;                                       /* intra */
1677    }
1678    
1679    
1680    #define INTRA_THRESHOLD 0.5
1681    #define BFRAME_SKIP_THRESHHOLD 30
1682    
1683    
1684    /* FrameCodeP also handles S(GMC)-VOPs */
1685    static int
1686    FrameCodeP(Encoder * pEnc,
1687                       Bitstream * bs,
1688                       uint32_t * pBits,
1689                       bool force_inter,
1690                       bool vol_header)
1691    {
1692            float fSigma;
1693    
1694            DECLARE_ALIGNED_MATRIX(dct_codes, 6, 64, int16_t, CACHE_LINE);
1695            DECLARE_ALIGNED_MATRIX(qcoeff, 6, 64, int16_t, CACHE_LINE);
1696    
1697            int mb_width = pEnc->mbParam.mb_width;
1698            int mb_height = pEnc->mbParam.mb_height;
1699    
1700            int iLimit;
1701            int x, y, k;
1702            int iSearchRange;
1703            int bIntra, skip_possible;
1704    
1705            /* IMAGE *pCurrent = &pEnc->current->image; */
1706            IMAGE *pRef = &pEnc->reference->image;
1707    
1708            if ((pEnc->current->global_flags & XVID_REDUCED))
1709            {
1710                    mb_width = (pEnc->mbParam.width + 31) / 32;
1711                    mb_height = (pEnc->mbParam.height + 31) / 32;
1712            }
1713    
1714            start_timer();
1715            image_setedges(pRef, pEnc->mbParam.edged_width, pEnc->mbParam.edged_height,
1716                                       pEnc->mbParam.width, pEnc->mbParam.height);
1717            stop_edges_timer();
1718    
1719            pEnc->mbParam.m_rounding_type = 1 - pEnc->mbParam.m_rounding_type;
1720            pEnc->current->rounding_type = pEnc->mbParam.m_rounding_type;
1721            pEnc->current->quarterpel =  pEnc->mbParam.m_quarterpel;
1722            pEnc->current->fcode = pEnc->mbParam.m_fcode;
1723    
1724            if (!force_inter)
1725                    iLimit = (int)(mb_width * mb_height *  INTRA_THRESHOLD);
1726            else
1727                    iLimit = mb_width * mb_height + 1;
1728    
1729            if ((pEnc->current->global_flags & XVID_HALFPEL)) {
1730                    start_timer();
1731                    image_interpolate(pRef, &pEnc->vInterH, &pEnc->vInterV,
1732                                                      &pEnc->vInterHV, pEnc->mbParam.edged_width,
1733                                                      pEnc->mbParam.edged_height,
1734                                                      pEnc->mbParam.m_quarterpel,
1735                                                      pEnc->current->rounding_type);
1736                    stop_inter_timer();
1737            }
1738    
1739            pEnc->current->coding_type = P_VOP;
1740    
1741            start_timer();
1742            if (pEnc->current->global_flags & XVID_HINTEDME_SET)
1743                    HintedMESet(pEnc, &bIntra);
1744            else
1745                    bIntra =
1746                            MotionEstimation(&pEnc->mbParam, pEnc->current, pEnc->reference,
1747                                                     &pEnc->vInterH, &pEnc->vInterV, &pEnc->vInterHV,
1748                                                     iLimit);
1749    
1750            stop_motion_timer();
1751    
1752            if (bIntra == 1) return FrameCodeI(pEnc, bs, pBits);
1753    
1754            if ( ( pEnc->current->global_flags & XVID_GMC )
1755                    && ( (pEnc->current->warp.duv[1].x != 0) || (pEnc->current->warp.duv[1].y != 0) ) )
1756            {
1757                    pEnc->current->coding_type = S_VOP;
1758    
1759                    generate_GMCparameters( 2, 16, &pEnc->current->warp,
1760                                            pEnc->mbParam.width, pEnc->mbParam.height,
1761                                            &pEnc->current->gmc_data);
1762    
1763                    generate_GMCimage(&pEnc->current->gmc_data, &pEnc->reference->image,
1764                                    pEnc->mbParam.mb_width, pEnc->mbParam.mb_height,
1765                                    pEnc->mbParam.edged_width, pEnc->mbParam.edged_width/2,
1766                                    pEnc->mbParam.m_fcode, pEnc->mbParam.m_quarterpel, 0,
1767                                    pEnc->current->rounding_type, pEnc->current->mbs, &pEnc->vGMC);
1768    
1769            }
1770    
1771            set_timecodes(pEnc->current,pEnc->reference,pEnc->mbParam.fbase);
1772            if (vol_header)
1773            {       BitstreamWriteVolHeader(bs, &pEnc->mbParam, pEnc->current);
1774                    BitstreamPadAlways(bs);
1775            }
1776    
1777            BitstreamWriteVopHeader(bs, &pEnc->mbParam, pEnc->current, 1);
1778    
1779            *pBits = BitstreamPos(bs);
1780    
1781            pEnc->current->sStat.iTextBits = pEnc->current->sStat.iMvSum = pEnc->current->sStat.iMvCount =
1782                    pEnc->current->sStat.kblks = pEnc->current->sStat.mblks = pEnc->current->sStat.ublks = 0;
1783    
1784    
1785            for (y = 0; y < mb_height; y++) {
1786                    for (x = 0; x < mb_width; x++) {
1787                            MACROBLOCK *pMB =
1788                                    &pEnc->current->mbs[x + y * pEnc->mbParam.mb_width];
1789    
1790    /* Mode decision: Check, if the block should be INTRA / INTER or GMC-coded */
1791    /* For a start, leave INTRA decision as is, only choose only between INTER/GMC  - gruel, 9.1.2002 */
1792    
1793                            bIntra = (pMB->mode == MODE_INTRA) || (pMB->mode == MODE_INTRA_Q);
1794    
1795                            if (bIntra) {
1796                                    CodeIntraMB(pEnc, pMB);
1797                                    MBTransQuantIntra(&pEnc->mbParam, pEnc->current, pMB, x, y,
1798                                                                      dct_codes, qcoeff);
1799    
1800                                    start_timer();
1801                                    MBPrediction(pEnc->current, x, y, pEnc->mbParam.mb_width, qcoeff);
1802                                    stop_prediction_timer();
1803    
1804                                    pEnc->current->sStat.kblks++;
1805    
1806                                    MBCoding(pEnc->current, pMB, qcoeff, bs, &pEnc->current->sStat);
1807                                    stop_coding_timer();
1808                                    continue;
1809                            }
1810    
1811                            if (pEnc->current->coding_type == S_VOP) {
1812    
1813                                    int32_t iSAD = sad16(pEnc->current->image.y + 16*y*pEnc->mbParam.edged_width + 16*x,
1814                                            pEnc->vGMC.y + 16*y*pEnc->mbParam.edged_width + 16*x,
1815                                            pEnc->mbParam.edged_width, 65536);
1816    
1817                                    if (pEnc->current->motion_flags & PMV_CHROMA16) {
1818                                            iSAD += sad8(pEnc->current->image.u + 8*y*(pEnc->mbParam.edged_width/2) + 8*x,
1819                                            pEnc->vGMC.u + 8*y*(pEnc->mbParam.edged_width/2) + 8*x, pEnc->mbParam.edged_width/2);
1820    
1821                                            iSAD += sad8(pEnc->current->image.v + 8*y*(pEnc->mbParam.edged_width/2) + 8*x,
1822                                            pEnc->vGMC.v + 8*y*(pEnc->mbParam.edged_width/2) + 8*x, pEnc->mbParam.edged_width/2);
1823                                    }
1824    
1825                                    if (iSAD <= pMB->sad16) {               /* mode decision GMC */
1826    
1827                                            if (pEnc->mbParam.m_quarterpel)
1828                                                    pMB->qmvs[0] = pMB->qmvs[1] = pMB->qmvs[2] = pMB->qmvs[3] = pMB->amv;
1829                                            else
1830                                                    pMB->mvs[0] = pMB->mvs[1] = pMB->mvs[2] = pMB->mvs[3] = pMB->amv;
1831    
1832                                            pMB->mode = MODE_INTER;
1833                                            pMB->mcsel = 1;
1834                                            pMB->sad16 = iSAD;
1835                                    } else {
1836                                            pMB->mcsel = 0;
1837                                    }
1838                            } else {
1839                                    pMB->mcsel = 0; /* just a precaution */
1840                            }
1841    
1842                            start_timer();
1843                            MBMotionCompensation(pMB, x, y, &pEnc->reference->image,
1844                                                                     &pEnc->vInterH, &pEnc->vInterV,
1845                                                                     &pEnc->vInterHV, &pEnc->vGMC,
1846                                                                     &pEnc->current->image,
1847                                                                     dct_codes, pEnc->mbParam.width,
1848                                                                     pEnc->mbParam.height,
1849                                                                     pEnc->mbParam.edged_width,
1850                                                                     pEnc->mbParam.m_quarterpel,
1851                                                                     (pEnc->current->global_flags & XVID_REDUCED),
1852                                                                     pEnc->current->rounding_type);
1853    
1854                            stop_comp_timer();
1855    
1856                            if ((pEnc->current->global_flags & XVID_LUMIMASKING)) {
1857                                    if (pMB->dquant != NO_CHANGE) {
1858                                            pMB->mode = MODE_INTER_Q;
1859                                            pEnc->current->quant += DQtab[pMB->dquant];
1860                                            if (pEnc->current->quant > 31)
1861                                                    pEnc->current->quant = 31;
1862                                            else if (pEnc->current->quant < 1)
1863                                                    pEnc->current->quant = 1;
1864                                    }
1865                            }
1866                            pMB->quant = pEnc->current->quant;
1867    
1868                            pMB->field_pred = 0;
1869    
1870                            if (pMB->mode != MODE_NOT_CODED)
1871                            {       pMB->cbp =
1872                                            MBTransQuantInter(&pEnc->mbParam, pEnc->current, pMB, x, y,
1873                                                                              dct_codes, qcoeff);
1874                            }
1875    
1876                            if (pMB->cbp || pMB->mvs[0].x || pMB->mvs[0].y ||
1877                                       pMB->mvs[1].x || pMB->mvs[1].y || pMB->mvs[2].x ||
1878                                       pMB->mvs[2].y || pMB->mvs[3].x || pMB->mvs[3].y) {
1879                                    pEnc->current->sStat.mblks++;
1880                            }  else {
1881                                    pEnc->current->sStat.ublks++;
1882                          }                          }
1883    
1884                          start_timer();                          start_timer();
1885                          MBCoding(&pEnc->mbParam, pMB, qcoeff, bs, &pEnc->sStat);  
1886                            /* Finished processing the MB, now check if to CODE or SKIP */
1887    
1888                            skip_possible = (pMB->cbp == 0) && (pMB->mode == MODE_INTER) &&
1889                                                            (pMB->dquant == NO_CHANGE);
1890    
1891                            if (pEnc->current->coding_type == S_VOP)
1892                                    skip_possible &= (pMB->mcsel == 1);
1893                            else if (pEnc->current->coding_type == P_VOP) {
1894                                    if (pEnc->mbParam.m_quarterpel)
1895                                            skip_possible &= ( (pMB->qmvs[0].x == 0) && (pMB->qmvs[0].y == 0) );
1896                                    else
1897                                            skip_possible &= ( (pMB->mvs[0].x == 0) && (pMB->mvs[0].y == 0) );
1898                            }
1899    
1900                            if ( (pMB->mode == MODE_NOT_CODED) || (skip_possible)) {
1901    
1902    /* This is a candidate for SKIPping, but for P-VOPs check intermediate B-frames first */
1903    
1904                                    if (pEnc->current->coding_type == P_VOP)        /* special rule for P-VOP's SKIP */
1905                                    {
1906                                            int bSkip = 1;
1907    
1908                                            for (k=pEnc->bframenum_head; k< pEnc->bframenum_tail; k++)
1909                                            {
1910                                                    int iSAD;
1911                                                    iSAD = sad16(pEnc->reference->image.y + 16*y*pEnc->mbParam.edged_width + 16*x,
1912                                                                            pEnc->bframes[k]->image.y + 16*y*pEnc->mbParam.edged_width + 16*x,
1913                                                                    pEnc->mbParam.edged_width,BFRAME_SKIP_THRESHHOLD);
1914                                                    if (iSAD >= BFRAME_SKIP_THRESHHOLD * pMB->quant)
1915                                                    {       bSkip = 0;
1916                                                            break;
1917                                                    }
1918                                            }
1919    
1920                                            if (!bSkip) {   /* no SKIP, but trivial block */
1921                                                    if(pEnc->mbParam.m_quarterpel) {
1922                                                            VECTOR predMV = get_qpmv2(pEnc->current->mbs, pEnc->mbParam.mb_width, 0, x, y, 0);
1923                                                            pMB->pmvs[0].x = - predMV.x;
1924                                                            pMB->pmvs[0].y = - predMV.y;
1925                                                    }
1926                                                    else {
1927                                                            VECTOR predMV = get_pmv2(pEnc->current->mbs, pEnc->mbParam.mb_width, 0, x, y, 0);
1928                                                            pMB->pmvs[0].x = - predMV.x;
1929                                                            pMB->pmvs[0].y = - predMV.y;
1930                                                    }
1931                                                    pMB->mode = MODE_INTER;
1932                                                    pMB->cbp = 0;
1933                                                    MBCoding(pEnc->current, pMB, qcoeff, bs, &pEnc->current->sStat);
1934                                                    stop_coding_timer();
1935    
1936                                                    continue;       /* next MB */
1937                                            }
1938                                    }
1939                                    /* do SKIP */
1940    
1941                                    pMB->mode = MODE_NOT_CODED;
1942                                    MBSkip(bs);
1943                                    stop_coding_timer();
1944                                    continue;       /* next MB */
1945                            }
1946                            /* ordinary case: normal coded INTER/INTER4V block */
1947    
1948                            if (pEnc->current->global_flags & XVID_GREYSCALE)
1949                            {       pMB->cbp &= 0x3C;               /* keep only bits 5-2 */
1950                                    qcoeff[4*64+0]=0;               /* zero, because DC for INTRA MBs DC value is saved */
1951                                    qcoeff[5*64+0]=0;
1952                            }
1953    
1954                            if(pEnc->mbParam.m_quarterpel) {
1955                                    VECTOR predMV = get_qpmv2(pEnc->current->mbs, pEnc->mbParam.mb_width, 0, x, y, 0);
1956                                    pMB->pmvs[0].x = pMB->qmvs[0].x - predMV.x;
1957                                    pMB->pmvs[0].y = pMB->qmvs[0].y - predMV.y;
1958                                    DPRINTF(DPRINTF_MV,"mv_diff (%i,%i) pred (%i,%i) result (%i,%i)", pMB->pmvs[0].x, pMB->pmvs[0].y, predMV.x, predMV.y, pMB->mvs[0].x, pMB->mvs[0].y);
1959                            } else {
1960                                    VECTOR predMV = get_pmv2(pEnc->current->mbs, pEnc->mbParam.mb_width, 0, x, y, 0);
1961                                    pMB->pmvs[0].x = pMB->mvs[0].x - predMV.x;
1962                                    pMB->pmvs[0].y = pMB->mvs[0].y - predMV.y;
1963                                    DPRINTF(DPRINTF_MV,"mv_diff (%i,%i) pred (%i,%i) result (%i,%i)", pMB->pmvs[0].x, pMB->pmvs[0].y, predMV.x, predMV.y, pMB->mvs[0].x, pMB->mvs[0].y);
1964                            }
1965    
1966    
1967                            if (pMB->mode == MODE_INTER4V)
1968                            {       int k;
1969                                    for (k=1;k<4;k++)
1970                                    {
1971                                            if(pEnc->mbParam.m_quarterpel) {
1972                                                    VECTOR predMV = get_qpmv2(pEnc->current->mbs, pEnc->mbParam.mb_width, 0, x, y, k);
1973                                                    pMB->pmvs[k].x = pMB->qmvs[k].x - predMV.x;
1974                                                    pMB->pmvs[k].y = pMB->qmvs[k].y - predMV.y;
1975                                    DPRINTF(DPRINTF_MV,"mv_diff (%i,%i) pred (%i,%i) result (%i,%i)", pMB->pmvs[k].x, pMB->pmvs[k].y, predMV.x, predMV.y, pMB->mvs[k].x, pMB->mvs[k].y);
1976                                            } else {
1977                                                    VECTOR predMV = get_pmv2(pEnc->current->mbs, pEnc->mbParam.mb_width, 0, x, y, k);
1978                                                    pMB->pmvs[k].x = pMB->mvs[k].x - predMV.x;
1979                                                    pMB->pmvs[k].y = pMB->mvs[k].y - predMV.y;
1980                                    DPRINTF(DPRINTF_MV,"mv_diff (%i,%i) pred (%i,%i) result (%i,%i)", pMB->pmvs[k].x, pMB->pmvs[k].y, predMV.x, predMV.y, pMB->mvs[k].x, pMB->mvs[k].y);
1981                                            }
1982    
1983                                    }
1984                            }
1985    
1986                            MBCoding(pEnc->current, pMB, qcoeff, bs, &pEnc->current->sStat);
1987                          stop_coding_timer();                          stop_coding_timer();
1988    
1989                    }
1990                  }                  }
1991    
1992            if ((pEnc->current->global_flags & XVID_REDUCED))
1993            {
1994                    image_deblock_rrv(&pEnc->current->image, pEnc->mbParam.edged_width,
1995                            pEnc->current->mbs, mb_width, mb_height, pEnc->mbParam.mb_width,
1996                            16, XVID_DEC_DEBLOCKY|XVID_DEC_DEBLOCKUV);
1997          }          }
1998    
1999          emms();          emms();
2000    
2001          if (pEnc->sStat.iMvCount == 0)          if (pEnc->current->global_flags & XVID_HINTEDME_GET) {
2002                  pEnc->sStat.iMvCount = 1;                  HintedMEGet(pEnc, 0);
2003            }
2004    
2005            if (pEnc->current->sStat.iMvCount == 0)
2006                    pEnc->current->sStat.iMvCount = 1;
2007    
2008          fSigma = (float)sqrt((float) pEnc->sStat.iMvSum / pEnc->sStat.iMvCount);          fSigma = (float) sqrt((float) pEnc->current->sStat.iMvSum / pEnc->current->sStat.iMvCount);
2009    
2010          iSearchRange = 1 << (3 + pEnc->mbParam.fixed_code);          iSearchRange = 1 << (3 + pEnc->mbParam.m_fcode);
2011    
2012          if ((fSigma > iSearchRange / 3)          if ((fSigma > iSearchRange / 3)
2013              && (pEnc->mbParam.fixed_code <= 3)) // maximum search range 128                  && (pEnc->mbParam.m_fcode <= (3 + pEnc->mbParam.m_quarterpel))) /* maximum search range 128 */
2014          {          {
2015                  pEnc->mbParam.fixed_code++;                  pEnc->mbParam.m_fcode++;
2016                  iSearchRange *= 2;                  iSearchRange *= 2;
2017          }          } else if ((fSigma < iSearchRange / 6)
2018          else if ((fSigma < iSearchRange / 6)                             && (pEnc->fMvPrevSigma >= 0)
2019                   && (pEnc->sStat.fMvPrevSigma >= 0)                             && (pEnc->fMvPrevSigma < iSearchRange / 6)
2020                   && (pEnc->sStat.fMvPrevSigma < iSearchRange / 6)                             && (pEnc->mbParam.m_fcode >= (2 + pEnc->mbParam.m_quarterpel)))      /* minimum search range 16 */
                  && (pEnc->mbParam.fixed_code >= 2))    // minimum search range 16  
2021          {          {
2022                  pEnc->mbParam.fixed_code--;                  pEnc->mbParam.m_fcode--;
2023                  iSearchRange /= 2;                  iSearchRange /= 2;
2024          }          }
2025    
2026          pEnc->sStat.fMvPrevSigma = fSigma;          pEnc->fMvPrevSigma = fSigma;
2027    
2028            /* frame drop code */
2029            DPRINTF(DPRINTF_DEBUG, "kmu %i %i %i", pEnc->current->sStat.kblks, pEnc->current->sStat.mblks, pEnc->current->sStat.ublks);
2030            if (pEnc->current->sStat.kblks + pEnc->current->sStat.mblks <
2031                    (pEnc->mbParam.frame_drop_ratio * mb_width * mb_height) / 100)
2032            {
2033                    pEnc->current->sStat.kblks = pEnc->current->sStat.mblks = 0;
2034                    pEnc->current->sStat.ublks = mb_width * mb_height;
2035    
2036                    BitstreamReset(bs);
2037    
2038                    set_timecodes(pEnc->current,pEnc->reference,pEnc->mbParam.fbase);
2039                    BitstreamWriteVopHeader(bs, &pEnc->mbParam, pEnc->current, 0);
2040    
2041                    /* copy reference frame details into the current frame */
2042                    pEnc->current->quant = pEnc->reference->quant;
2043                    pEnc->current->motion_flags = pEnc->reference->motion_flags;
2044                    pEnc->current->rounding_type = pEnc->reference->rounding_type;
2045                    pEnc->current->quarterpel =  pEnc->reference->quarterpel;
2046                    pEnc->current->fcode = pEnc->reference->fcode;
2047                    pEnc->current->bcode = pEnc->reference->bcode;
2048                    image_copy(&pEnc->current->image, &pEnc->reference->image, pEnc->mbParam.edged_width, pEnc->mbParam.height);
2049                    memcpy(pEnc->current->mbs, pEnc->reference->mbs, sizeof(MACROBLOCK) * mb_width * mb_height);
2050            }
2051    
2052            /* XXX: debug
2053            {
2054                    char s[100];
2055                    sprintf(s, "\\%05i_cur.pgm", pEnc->m_framenum);
2056                    image_dump_yuvpgm(&pEnc->current->image,
2057                            pEnc->mbParam.edged_width,
2058                            pEnc->mbParam.width, pEnc->mbParam.height, s);
2059    
2060                    sprintf(s, "\\%05i_ref.pgm", pEnc->m_framenum);
2061                    image_dump_yuvpgm(&pEnc->reference->image,
2062                            pEnc->mbParam.edged_width,
2063                            pEnc->mbParam.width, pEnc->mbParam.height, s);
2064            }
2065            */
2066    
2067    
2068            *pBits = BitstreamPos(bs) - *pBits;
2069    
2070            return 0;                                       /* inter */
2071    }
2072    
2073    
2074    static void
2075    FrameCodeB(Encoder * pEnc,
2076                       FRAMEINFO * frame,
2077                       Bitstream * bs,
2078                       uint32_t * pBits)
2079    {
2080            DECLARE_ALIGNED_MATRIX(dct_codes, 6, 64, int16_t, CACHE_LINE);
2081            DECLARE_ALIGNED_MATRIX(qcoeff, 6, 64, int16_t, CACHE_LINE);
2082            uint32_t x, y;
2083    
2084            IMAGE *f_ref = &pEnc->reference->image;
2085            IMAGE *b_ref = &pEnc->current->image;
2086    
2087    #ifdef BFRAMES_DEC_DEBUG
2088            FILE *fp;
2089            static char first=0;
2090    #define BFRAME_DEBUG    if (!first && fp){ \
2091                    fprintf(fp,"Y=%3d   X=%3d   MB=%2d   CBP=%02X\n",y,x,mb->mode,mb->cbp); \
2092            }
2093    
2094            pEnc->current->global_flags &= ~XVID_REDUCED;   /* reduced resoltion not yet supported */
2095    
2096            if (!first){
2097                    fp=fopen("C:\\XVIDDBGE.TXT","w");
2098            }
2099    #endif
2100    
2101            frame->quarterpel =  pEnc->mbParam.m_quarterpel;
2102    
2103            /* forward  */
2104            image_setedges(f_ref, pEnc->mbParam.edged_width,
2105                                       pEnc->mbParam.edged_height, pEnc->mbParam.width,
2106                                       pEnc->mbParam.height);
2107            start_timer();
2108            image_interpolate(f_ref, &pEnc->f_refh, &pEnc->f_refv, &pEnc->f_refhv,
2109                                              pEnc->mbParam.edged_width, pEnc->mbParam.edged_height,
2110                                              pEnc->mbParam.m_quarterpel, 0);
2111            stop_inter_timer();
2112    
2113            /* backward */
2114            image_setedges(b_ref, pEnc->mbParam.edged_width,
2115                                       pEnc->mbParam.edged_height, pEnc->mbParam.width,
2116                                       pEnc->mbParam.height);
2117            start_timer();
2118            image_interpolate(b_ref, &pEnc->vInterH, &pEnc->vInterV, &pEnc->vInterHV,
2119                                              pEnc->mbParam.edged_width, pEnc->mbParam.edged_height,
2120                                              pEnc->mbParam.m_quarterpel, 0);
2121            stop_inter_timer();
2122    
2123            start_timer();
2124    
2125            MotionEstimationBVOP(&pEnc->mbParam, frame,
2126                                                     ((int32_t)(pEnc->current->stamp - frame->stamp)),                              /* time_bp */
2127                                                     ((int32_t)(pEnc->current->stamp - pEnc->reference->stamp)),    /* time_pp */
2128                                                     pEnc->reference->mbs, f_ref,
2129                                                     &pEnc->f_refh, &pEnc->f_refv, &pEnc->f_refhv,
2130                                                     pEnc->current, b_ref, &pEnc->vInterH,
2131                                                     &pEnc->vInterV, &pEnc->vInterHV);
2132    
2133    
2134            stop_motion_timer();
2135            /*
2136            if (test_quant_type(&pEnc->mbParam, pEnc->current)) {
2137                    BitstreamWriteVolHeader(bs, pEnc->mbParam.width, pEnc->mbParam.height, pEnc->mbParam.quant_type);
2138            }
2139            */
2140    
2141            frame->coding_type = B_VOP;
2142    
2143            set_timecodes(frame, pEnc->reference,pEnc->mbParam.fbase);
2144            BitstreamWriteVopHeader(bs, &pEnc->mbParam, frame, 1);
2145    
2146            *pBits = BitstreamPos(bs);
2147    
2148            frame->sStat.iTextBits = 0;
2149            frame->sStat.iMvSum = 0;
2150            frame->sStat.iMvCount = 0;
2151            frame->sStat.kblks = frame->sStat.mblks = frame->sStat.ublks = 0;
2152    
2153    
2154            for (y = 0; y < pEnc->mbParam.mb_height; y++) {
2155                    for (x = 0; x < pEnc->mbParam.mb_width; x++) {
2156                            MACROBLOCK * const mb = &frame->mbs[x + y * pEnc->mbParam.mb_width];
2157                            int direction = pEnc->mbParam.global & XVID_ALTERNATESCAN ? 2 : 0;
2158    
2159                            /* decoder ignores mb when refence block is INTER(0,0), CBP=0 */
2160                            if (mb->mode == MODE_NOT_CODED) {
2161                                    /* mb->mvs[0].x = mb->mvs[0].y = mb->cbp = 0; */
2162                                    continue;
2163                            }
2164    
2165                            if (mb->mode != MODE_DIRECT_NONE_MV) {
2166                                    MBMotionCompensationBVOP(&pEnc->mbParam, mb, x, y, &frame->image,
2167                                                                             f_ref, &pEnc->f_refh, &pEnc->f_refv,
2168                                                                             &pEnc->f_refhv, b_ref, &pEnc->vInterH,
2169                                                                             &pEnc->vInterV, &pEnc->vInterHV,
2170                                                                             dct_codes);
2171    
2172                                    if (mb->mode == MODE_DIRECT_NO4V) mb->mode = MODE_DIRECT;
2173                                    mb->quant = frame->quant;
2174    
2175                                    mb->cbp =
2176                                            MBTransQuantInterBVOP(&pEnc->mbParam, frame, mb, dct_codes, qcoeff);
2177    
2178                                    if ( (mb->mode == MODE_DIRECT) && (mb->cbp == 0)
2179                                            && (mb->pmvs[3].x == 0) && (mb->pmvs[3].y == 0) ) {
2180                                            mb->mode = MODE_DIRECT_NONE_MV; /* skipped */
2181                                    }
2182                            }
2183    
2184    #ifdef BFRAMES_DEC_DEBUG
2185            BFRAME_DEBUG
2186    #endif
2187                            start_timer();
2188                            MBCodingBVOP(mb, qcoeff, frame->fcode, frame->bcode, bs,
2189                                                     &frame->sStat, direction);
2190                            stop_coding_timer();
2191                    }
2192            }
2193    
2194            emms();
2195    
2196            /* TODO: dynamic fcode/bcode ??? */
2197    
2198          *pBits = BitstreamPos(bs) - *pBits;          *pBits = BitstreamPos(bs) - *pBits;
2199    
2200          return 0;                                        // inter  #ifdef BFRAMES_DEC_DEBUG
2201            if (!first){
2202                    first=1;
2203                    if (fp)
2204                            fclose(fp);
2205            }
2206    #endif
2207  }  }

Legend:
Removed from v.1.19  
changed lines
  Added in v.1.96

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