]> git.sesse.net Git - vlc/blobdiff - src/video_parser/vpar_synchro.c
The motion compensation routines are now modules as well ; choose your
[vlc] / src / video_parser / vpar_synchro.c
index 81e495154421a35844a3f44d9ae591d346e65037..d96433e64db3f6e492a8dd620d2268df07da40db 100644 (file)
 /*****************************************************************************
- * vpar_motion.c : motion vectors parsing
- * (c)1999 VideoLAN
+ * vpar_synchro.c : frame dropping routines
+ *****************************************************************************
+ * Copyright (C) 1999, 2000 VideoLAN
+ * $Id: vpar_synchro.c,v 1.78 2001/01/18 05:13:23 sam Exp $
+ *
+ * Authors: Christophe Massiot <massiot@via.ecp.fr>
+ *          Samuel Hocevar <sam@via.ecp.fr>
+ *          Jean-Marc Dressler <polux@via.ecp.fr>
+ *
+ * This program is free software; you can redistribute it and/or modify
+ * it under the terms of the GNU General Public License as published by
+ * the Free Software Foundation; either version 2 of the License, or
+ * (at your option) any later version.
+ * 
+ * This program is distributed in the hope that it will be useful,
+ * but WITHOUT ANY WARRANTY; without even the implied warranty of
+ * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+ * GNU General Public License for more details.
+ *
+ * You should have received a copy of the GNU General Public License
+ * along with this program; if not, write to the Free Software
+ * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA  02111, USA.
  *****************************************************************************/
 
+/*
+ * DISCUSSION : How to Write an efficient Frame-Dropping Algorithm
+ * ==========
+ *
+ * This implementation is based on mathematical and statistical
+ * developments. Older implementations used an enslavement, considering
+ * that if we're late when reading an I picture, we will decode one frame
+ * less. It had a tendancy to derive, and wasn't responsive enough, which
+ * would have caused trouble with the stream control stuff.
+ *
+ * 1. Structure of a picture stream
+ *    =============================
+ * Between 2 I's, we have for instance :
+ *    I   B   P   B   P   B   P   B   P   B   P   B   I
+ *    t0  t1  t2  t3  t4  t5  t6  t7  t8  t9  t10 t11 t12
+ * Please bear in mind that B's and IP's will be inverted when displaying
+ * (decoding order != presentation order). Thus, t1 < t0.
+ *
+ * FIXME: write a few words about stream structure changes.
+ *
+ * 2. Definitions
+ *    ===========
+ * t[0..12]     : Presentation timestamps of pictures 0..12.
+ * t            : Current timestamp, at the moment of the decoding.
+ * T            : Picture period, T = 1/frame_rate.
+ * tau[I,P,B]   : Mean time to decode an [I,P,B] picture.
+ * tauYUV       : Mean time to render a picture (given by the video_output).
+ * tau´[I,P,B] = 2 * tau[I,P,B] + tauYUV
+ *              : Mean time + typical difference (estimated to tau/2, that
+ *                needs to be confirmed) + render time.
+ * DELTA        : A given error margin.
+ *
+ * 3. General considerations
+ *    ======================
+ * We define three types of machines :
+ *      14T > tauI : machines capable of decoding all I pictures
+ *      2T > tauP  : machines capable of decoding all P pictures
+ *      T > tauB   : machines capable of decoding all B pictures
+ *
+ * 4. Decoding of an I picture
+ *    ========================
+ * On fast machines, we decode all I's.
+ * Otherwise :
+ * We can decode an I picture if we simply have enough time to decode it 
+ * before displaying :
+ *      t0 - t > tau´I + DELTA
+ *
+ * 5. Decoding of a P picture
+ *    =======================
+ * On fast machines, we decode all P's.
+ * Otherwise :
+ * First criterion : have time to decode it.
+ *      t2 - t > tau´P + DELTA
+ *
+ * Second criterion : it shouldn't prevent us from displaying the forthcoming
+ * I picture, which is more important.
+ *      t12 - t > tau´P + tau´I + DELTA
+ *
+ * 6. Decoding of a B picture
+ *    =======================
+ * On fast machines, we decode all B's. Otherwise :
+ *      t1 - t > tau´B + DELTA
+ * Since the next displayed I or P is already decoded, we don't have to
+ * worry about it.
+ *
+ * I hope you will have a pleasant flight and do not forget your life
+ * jacket.
+ *                                                  --Meuuh (2000-12-29)
+ */
+
 /*****************************************************************************
  * Preamble
  *****************************************************************************/
-#include <errno.h>
-#include <stdlib.h>
-#include <stdio.h>
-#include <unistd.h>
-#include <string.h>
-#include <sys/uio.h>
+#include "defs.h"
 
 #include "config.h"
 #include "common.h"
+#include "threads.h"
 #include "mtime.h"
-#include "vlc_thread.h"
+#include "plugins.h"
 
 #include "intf_msg.h"
-#include "debug.h"                    /* ?? temporaire, requis par netlist.h */
 
-#include "input.h"
-#include "input_netlist.h"
-#include "decoder_fifo.h"
+#include "stream_control.h"
+#include "input_ext-dec.h"
+
 #include "video.h"
 #include "video_output.h"
 
-#include "vdec_idct.h"
 #include "video_decoder.h"
 #include "vdec_motion.h"
+#include "../video_decoder/vdec_idct.h"
 
 #include "vpar_blocks.h"
-#include "vpar_headers.h"
-#include "vpar_synchro.h"
-#include "video_parser.h"
-#include "video_fifo.h"
+#include "../video_decoder/vpar_headers.h"
+#include "../video_decoder/vpar_synchro.h"
+#include "../video_decoder/video_parser.h"
 
-#define MAX_COUNT 3
+#include "main.h"
 
 /*
  * Local prototypes
  */
+static int  SynchroType( void );
+
+/* Error margins */
+#define DELTA                   (int)(0.040*CLOCK_FREQ)
+
+#define DEFAULT_NB_P            5
+#define DEFAULT_NB_B            1
 
 /*****************************************************************************
- * vpar_SynchroUpdateTab : Update a mean table in the synchro structure
+ * vpar_SynchroInit : You know what ?
  *****************************************************************************/
-float vpar_SynchroUpdateTab( video_synchro_tab_t * tab, int count )
+void vpar_SynchroInit( vpar_thread_t * p_vpar )
 {
-       
-    tab->mean = ( tab->mean + MAX_COUNT * count ) / ( MAX_COUNT + 1 );
-    tab->deviation = ( tab->deviation + MAX_COUNT * abs (tab->mean - count) )
-                        / ( MAX_COUNT + 1 );
-
-    return tab->deviation;
+    p_vpar->synchro.i_type = SynchroType();
+    p_vpar->synchro.i_start = p_vpar->synchro.i_end = 0;
+    vlc_mutex_init( &p_vpar->synchro.fifo_lock );
+
+    /* We use a fake stream pattern, which is often right. */
+    p_vpar->synchro.i_n_p = p_vpar->synchro.i_eta_p = DEFAULT_NB_P;
+    p_vpar->synchro.i_n_b = p_vpar->synchro.i_eta_b = DEFAULT_NB_B;
+    memset( p_vpar->synchro.p_tau, 0, 4 * sizeof(mtime_t) );
+    memset( p_vpar->synchro.pi_meaningful, 0, 4 * sizeof(unsigned int) );
+    p_vpar->synchro.b_dropped_last = 0;
+    p_vpar->synchro.current_pts = mdate() + DEFAULT_PTS_DELAY;
+    p_vpar->synchro.backward_pts = 0;
+    p_vpar->synchro.i_current_period = p_vpar->synchro.i_backward_period = 0;
+#ifdef STATS
+    p_vpar->synchro.i_trashed_pic = p_vpar->synchro.i_not_chosen_pic = 
+        p_vpar->synchro.i_pic = 0;
+#endif
 }
 
 /*****************************************************************************
- * vpar_SynchroUpdateStructures : Update the synchro structures
+ * vpar_SynchroChoose : Decide whether we will decode a picture or not
  *****************************************************************************/
-void vpar_SynchroUpdateStructures( vpar_thread_t * p_vpar,
-                                   int i_coding_type, int dropped )
+boolean_t vpar_SynchroChoose( vpar_thread_t * p_vpar, int i_coding_type,
+                              int i_structure )
 {
-    float candidate_deviation;
-    float optimal_deviation;
-    float predict;
-    mtime_t i_current_pts;
-    mtime_t i_delay;
-    mtime_t i_displaydate;
-    decoder_fifo_t * decoder_fifo = p_vpar->bit_stream.p_decoder_fifo;
-
-    /* interpolate the current _decode_ PTS */
-    i_current_pts = decoder_fifo->buffer[decoder_fifo->i_start]->i_pts;
-    if( !i_current_pts )
-    {
-        i_current_pts = p_vpar->synchro.i_last_decode_pts
-                       + 1000000.0 / (1 + p_vpar->synchro.actual_fps);
-    }
-    p_vpar->synchro.i_last_decode_pts = i_current_pts;
-    /* see if the current image has a pts - if not, set to 0 */
-    p_vpar->synchro.fifo[p_vpar->synchro.i_fifo_stop].i_pts
-           = i_current_pts;
-
-    /* update display time */
-    i_displaydate = decoder_fifo->buffer[decoder_fifo->i_start]->i_pts;
-    if( !i_displaydate || i_coding_type != I_CODING_TYPE )
-    {
-        if (!p_vpar->synchro.i_images_since_pts )
-            p_vpar->synchro.i_images_since_pts = 10;
-
-        i_displaydate = p_vpar->synchro.i_last_display_pts
-                       + 1000000.0 / (p_vpar->synchro.theorical_fps);
-        //fprintf (stderr, "  ");
-    }
-    
-    /* else fprintf (stderr, "R ");
-    if (dropped) fprintf (stderr, "  "); else fprintf (stderr, "* ");
-    fprintf (stderr, "%i ", i_coding_type);
-    fprintf (stderr, "pts %lli delta %lli\n", i_displaydate, i_displaydate - p_vpar->synchro.i_last_display_pts); */
-
-    p_vpar->synchro.i_images_since_pts--;
-    p_vpar->synchro.i_last_display_pts = i_displaydate;
-
+    /* For clarity reasons, we separated the special synchros code from the
+     * mathematical synchro */
 
-
-    /* update structures */
-    switch(i_coding_type)
+    if( p_vpar->synchro.i_type != VPAR_SYNCHRO_DEFAULT )
     {
-        case P_CODING_TYPE:
-
-            p_vpar->synchro.current_p_count++;
-            if( !dropped ) p_vpar->synchro.nondropped_p_count++;
-            break;
-
-        case B_CODING_TYPE:
-            p_vpar->synchro.current_b_count++;
-            if( !dropped ) p_vpar->synchro.nondropped_b_count++;
-            break;
-
+        switch( i_coding_type )
+        {
         case I_CODING_TYPE:
+            /* I, IP, IP+, IPB */
+            if( p_vpar->synchro.i_type == VPAR_SYNCHRO_Iplus )
+            {
+                p_vpar->synchro.b_dropped_last = 1;
+            }
+            return( 1 );
 
-            /* update information about images we can decode */
-           if (i_current_pts != p_vpar->synchro.i_last_i_pts)
+        case P_CODING_TYPE:
+            if( p_vpar->synchro.i_type == VPAR_SYNCHRO_I ) /* I */
             {
-                if ( p_vpar->synchro.i_last_i_pts && i_current_pts != p_vpar->synchro.i_last_i_pts)
-                {
-                    p_vpar->synchro.theorical_fps = (p_vpar->synchro.theorical_fps + 1000000.0 * (1 + p_vpar->synchro.current_b_count + p_vpar->synchro.current_p_count) / (i_current_pts - p_vpar->synchro.i_last_i_pts)) / 2;
-                }
-                p_vpar->synchro.i_last_i_pts = i_current_pts;
+                return( 0 );
             }
 
-            if( !dropped )
+            if( p_vpar->synchro.i_type == VPAR_SYNCHRO_Iplus ) /* I+ */
             {
-                if ( p_vpar->synchro.i_last_nondropped_i_pts && i_current_pts != p_vpar->synchro.i_last_nondropped_i_pts)
+                if( p_vpar->synchro.b_dropped_last )
+                {
+                    p_vpar->synchro.b_dropped_last = 0;
+                    return( 1 );
+                }
+                else
                 {
-                    p_vpar->synchro.actual_fps = (p_vpar->synchro.actual_fps + 1000000.0 * (1 + p_vpar->synchro.nondropped_b_count + p_vpar->synchro.nondropped_p_count) / (i_current_pts - p_vpar->synchro.i_last_nondropped_i_pts)) / 2;
+                    return( 0 );
                 }
-    
             }
 
+            return( 1 ); /* IP, IP+, IPB */
 
-            /* update all the structures for P images */
-
-            /* period == 1 */
-            optimal_deviation = vpar_SynchroUpdateTab(
-                            &p_vpar->synchro.tab_p[0],
-                            p_vpar->synchro.current_p_count);
-            predict = p_vpar->synchro.tab_p[0].mean;
-
-            /* period == 2 */
-            candidate_deviation = vpar_SynchroUpdateTab(
-                            &p_vpar->synchro.tab_p[1 + (p_vpar->synchro.modulo & 0x1)],
-                            p_vpar->synchro.current_p_count);
-            if (candidate_deviation < optimal_deviation)
-           {
-                optimal_deviation = candidate_deviation;
-                predict = p_vpar->synchro.tab_p[1 + (p_vpar->synchro.modulo & 0x1)].mean;
+        case B_CODING_TYPE:
+            if( p_vpar->synchro.i_type <= VPAR_SYNCHRO_IP ) /* I, IP */
+            {
+                return( 0 );
             }
-
-            /* period == 3 */
-            candidate_deviation = vpar_SynchroUpdateTab(
-                            &p_vpar->synchro.tab_p[3 + (p_vpar->synchro.modulo % 3)],
-                            p_vpar->synchro.current_p_count);
-            if (candidate_deviation < optimal_deviation)
-           {
-                optimal_deviation = candidate_deviation;
-                predict = p_vpar->synchro.tab_p[1 + (p_vpar->synchro.modulo % 3)].mean;
+            else if( p_vpar->synchro.i_type == VPAR_SYNCHRO_IPB ) /* IPB */
+            {
+                return( 1 );
             }
 
-           p_vpar->synchro.p_count_predict = predict;
-            p_vpar->synchro.current_p_count = 0;
+            p_vpar->synchro.b_dropped_last ^= 1; /* IP+ */
+            return( !p_vpar->synchro.b_dropped_last );
+        }
+        return( 0 ); /* never reached but gcc yells at me */
+    }
+    else
+    {
+#define TAU_PRIME( coding_type )    (p_vpar->synchro.p_tau[(coding_type)] \
+                                 + (p_vpar->synchro.p_tau[(coding_type)] >> 1) \
+                                            + tau_yuv)
+#define S                           p_vpar->synchro
+        /* VPAR_SYNCHRO_DEFAULT */
+        mtime_t         now, pts, period, tau_yuv;
+        boolean_t       b_decode = 0;
+#ifdef DEBUG_VPAR
+        char            p_date[MSTRTIME_MAX_SIZE];
+#endif
 
+        now = mdate();
+        period = 1000000 * 1001 / p_vpar->sequence.i_frame_rate;
 
-            /* update all the structures for B images */
+        vlc_mutex_lock( &p_vpar->p_vout->change_lock );
+        tau_yuv = p_vpar->p_vout->render_time;
+        vlc_mutex_unlock( &p_vpar->p_vout->change_lock );
 
-            /* period == 1 */
-            optimal_deviation = vpar_SynchroUpdateTab(
-                            &p_vpar->synchro.tab_b[0],
-                            p_vpar->synchro.current_b_count);
-            predict = p_vpar->synchro.tab_b[0].mean;
+        vlc_mutex_lock( &p_vpar->synchro.fifo_lock );
 
-            /* period == 2 */
-            candidate_deviation = vpar_SynchroUpdateTab(
-                            &p_vpar->synchro.tab_b[1 + (p_vpar->synchro.modulo & 0x1)],
-                            p_vpar->synchro.current_b_count);
-            if (candidate_deviation < optimal_deviation)
-           {
-                optimal_deviation = candidate_deviation;
-                predict = p_vpar->synchro.tab_b[1 + (p_vpar->synchro.modulo & 0x1)].mean;
+        switch( i_coding_type )
+        {
+        case I_CODING_TYPE:
+            if( S.backward_pts )
+            {
+                pts = S.backward_pts;
             }
-
-            /* period == 3 */
-            candidate_deviation = vpar_SynchroUpdateTab(
-                            &p_vpar->synchro.tab_b[3 + (p_vpar->synchro.modulo % 3)],
-                            p_vpar->synchro.current_b_count);
-            if (candidate_deviation < optimal_deviation)
-           {
-                optimal_deviation = candidate_deviation;
-                predict = p_vpar->synchro.tab_b[1 + (p_vpar->synchro.modulo % 3)].mean;
+            else
+            {
+                /* displaying order : B B P B B I
+                 *                      ^       ^
+                 *                      |       +- current picture
+                 *                      +- current PTS
+                 */
+                pts = S.current_pts + period * (S.i_n_b + 2);
             }
 
-           p_vpar->synchro.b_count_predict = predict;
-            p_vpar->synchro.current_b_count = 0;
-           
-            /* now we calculated all statistics, it's time to
-             * decide what we have the time to display
-             */
-            i_delay = i_current_pts - p_vpar->synchro.i_last_nondropped_i_pts;
-
-           p_vpar->synchro.can_display_i
-                = ( p_vpar->synchro.i_mean_decode_time < i_delay );
-
-            p_vpar->synchro.can_display_p
-                    = ( p_vpar->synchro.i_mean_decode_time
-                    * (1 + p_vpar->synchro.p_count_predict) < i_delay );
-
-            if( !p_vpar->synchro.can_display_p )
+            if( (1 + S.i_n_p * (S.i_n_b + 1)) * period >
+                    S.p_tau[I_CODING_TYPE] )
             {
-                p_vpar->synchro.displayable_p
-                    = -1 + i_delay / p_vpar->synchro.i_mean_decode_time;
-                if( p_vpar->synchro.displayable_p < 0 )
-                    p_vpar->synchro.displayable_p = 0;
-            }
-           else
-                p_vpar->synchro.displayable_p = 0;
-
-           if( p_vpar->synchro.can_display_p
-                && !(p_vpar->synchro.can_display_b 
-                    = ( p_vpar->synchro.i_mean_decode_time
-                    * (1 + p_vpar->synchro.b_count_predict
-                        + p_vpar->synchro.p_count_predict)) < i_delay) )
-           {
-                p_vpar->synchro.displayable_b
-                    = -2.0 + i_delay / p_vpar->synchro.i_mean_decode_time
-                        - p_vpar->synchro.can_display_p;
+                b_decode = 1;
             }
             else
-                p_vpar->synchro.displayable_b = 0;
-
-#if 0
-            fprintf( stderr,
-                "I %i  P %i (%f)  B %i (%f)\n",
-                p_vpar->synchro.can_display_i,
-                p_vpar->synchro.can_display_p,
-                p_vpar->synchro.displayable_p,
-                p_vpar->synchro.can_display_b,
-                p_vpar->synchro.displayable_b );
-#endif
-
-            /* update some values */
-            if( !dropped )
             {
-                p_vpar->synchro.i_last_nondropped_i_pts = i_current_pts;
-                p_vpar->synchro.nondropped_p_count = 0;
-                p_vpar->synchro.nondropped_b_count = 0;
+                b_decode = (pts - now) > (TAU_PRIME(I_CODING_TYPE) + DELTA);
             }
-
+            if( !b_decode )
+                intf_WarnMsg( 3, "vpar synchro warning: trashing I" );
             break;
 
-    }
-
-    p_vpar->synchro.modulo++;
-
-}
-
-/*****************************************************************************
- * vpar_SynchroChoose : Decide whether we will decode a picture or not
- *****************************************************************************/
-boolean_t vpar_SynchroChoose( vpar_thread_t * p_vpar, int i_coding_type,
-                              int i_structure )
-{
-    mtime_t i_delay = p_vpar->synchro.i_last_decode_pts - mdate();
-
-    switch( i_coding_type )
-    {
-        case I_CODING_TYPE:
-
-            return( p_vpar->synchro.can_display_i );
-
         case P_CODING_TYPE:
+            if( S.backward_pts )
+            {
+                pts = S.backward_pts;
+            }
+            else
+            {
+                pts = S.current_pts + period * (S.i_n_b + 1);
+            }
 
-            if( p_vpar->synchro.can_display_p )
-                return( 1 );
-
-            if( p_vpar->synchro.displayable_p * i_delay
-                < p_vpar->synchro.i_mean_decode_time )
+            if( (1 + S.i_n_p * (S.i_n_b + 1)) * period >
+                    S.p_tau[I_CODING_TYPE] )
             {
-                //fprintf( stderr, "trashed a P\n");
-                return( 0 );
+                if( (S.i_n_b + 1) * period > S.p_tau[P_CODING_TYPE] )
+                {
+                    /* Security in case we're _really_ late */
+                    b_decode = (pts - now > 0);
+                }
+                else
+                {
+                    b_decode = (pts - now) > (TAU_PRIME(P_CODING_TYPE) + DELTA);
+                    /* next I */
+                    b_decode &= (pts - now
+                                  + period
+                              * ( (S.i_n_p - S.i_eta_p) * (1 + S.i_n_b) - 1 ))
+                                > (TAU_PRIME(P_CODING_TYPE)
+                                    + TAU_PRIME(I_CODING_TYPE) + DELTA);
+                }
+            }
+            else
+            {
+                b_decode = 0;
             }
+            break;
 
-            p_vpar->synchro.displayable_p--;
-            return( 1 );
-   
         case B_CODING_TYPE:
+            pts = S.current_pts;
 
-            if( p_vpar->synchro.can_display_b )
-                return( 1 );
-
-            /* modulo & 0x3 is here to add some randomness */
-            if( i_delay < (1 + (p_vpar->synchro.modulo & 0x3))
-                * p_vpar->synchro.i_mean_decode_time )
+            if( (S.i_n_b + 1) * period > S.p_tau[P_CODING_TYPE] )
             {
-                //fprintf( stderr, "trashed a B\n");
-                return( 0 );
+                b_decode = (pts - now) > (TAU_PRIME(B_CODING_TYPE) + DELTA);
             }
-            if( p_vpar->synchro.displayable_b <= 0 )
-                return( 0 );
-
-            p_vpar->synchro.displayable_b--;
-            return( 1 );
+            else
+            {
+                b_decode = 0;
+            }
+        }
+
+        vlc_mutex_unlock( &p_vpar->synchro.fifo_lock );
+#ifdef DEBUG_VPAR
+        intf_DbgMsg("vpar synchro debug: %s picture scheduled for %s, %s (%lld)",
+                    i_coding_type == B_CODING_TYPE ? "B" :
+                    (i_coding_type == P_CODING_TYPE ? "P" : "I"),
+                    mstrtime(p_date, pts), b_decode ? "decoding" : "trashed",
+                    S.p_tau[i_coding_type]);
+#endif
+#ifdef STATS
+        if( !b_decode )
+        {
+            S.i_not_chosen_pic++;
+        }
+#endif
+        return( b_decode );
+#undef S
+#undef TAU_PRIME
     }
-
-    return( 0 );
-
 }
 
 /*****************************************************************************
- * vpar_SynchroTrash : Update timers when we trash a picture
+ * vpar_SynchroTrash : Update counters when we trash a picture
  *****************************************************************************/
 void vpar_SynchroTrash( vpar_thread_t * p_vpar, int i_coding_type,
                         int i_structure )
 {
-    vpar_SynchroUpdateStructures (p_vpar, i_coding_type, 1);
-
+#ifdef STATS
+    p_vpar->synchro.i_trashed_pic++;
+#endif
 }
 
 /*****************************************************************************
  * vpar_SynchroDecode : Update timers when we decide to decode a picture
  *****************************************************************************/
 void vpar_SynchroDecode( vpar_thread_t * p_vpar, int i_coding_type,
-                            int i_structure )
+                         int i_structure )
 {
-    vpar_SynchroUpdateStructures (p_vpar, i_coding_type, 0);
+    vlc_mutex_lock( &p_vpar->synchro.fifo_lock );
 
-    p_vpar->synchro.fifo[p_vpar->synchro.i_fifo_stop].i_decode_date = mdate();
-    p_vpar->synchro.fifo[p_vpar->synchro.i_fifo_stop].i_image_type
-        = i_coding_type;
-
-    p_vpar->synchro.i_fifo_stop = (p_vpar->synchro.i_fifo_stop + 1) & 0xf;
+    if( ((p_vpar->synchro.i_end + 1 - p_vpar->synchro.i_start)
+            % MAX_DECODING_PIC) )
+    {
+        p_vpar->synchro.p_date_fifo[p_vpar->synchro.i_end] = mdate();
+        p_vpar->synchro.pi_coding_types[p_vpar->synchro.i_end] = i_coding_type;
 
+        FIFO_INCREMENT( i_end );
+    }
+    else
+    {
+        /* FIFO full, panic() */
+        intf_ErrMsg("vpar error: synchro fifo full, estimations will be biased");
+    }
+    vlc_mutex_unlock( &p_vpar->synchro.fifo_lock );
 }
 
 /*****************************************************************************
  * vpar_SynchroEnd : Called when the image is totally decoded
  *****************************************************************************/
-void vpar_SynchroEnd( vpar_thread_t * p_vpar )
+void vpar_SynchroEnd( vpar_thread_t * p_vpar, int i_garbage )
 {
-    mtime_t i_decode_time;
+    mtime_t     tau;
+    int         i_coding_type;
 
-    i_decode_time = (mdate() -
-            p_vpar->synchro.fifo[p_vpar->synchro.i_fifo_start].i_decode_date)
-        / (p_vpar->synchro.i_fifo_stop - p_vpar->synchro.i_fifo_start & 0x0f);
+    vlc_mutex_lock( &p_vpar->synchro.fifo_lock );
 
-    p_vpar->synchro.i_mean_decode_time =
-        ( 7 * p_vpar->synchro.i_mean_decode_time + i_decode_time ) / 8;
-
-    /* fprintf (stderr,
-        "decoding time was %lli\n",
-        p_vpar->synchro.i_mean_decode_time); */
+    if (!i_garbage)
+    {
+        tau = mdate() - p_vpar->synchro.p_date_fifo[p_vpar->synchro.i_start];
+        i_coding_type = p_vpar->synchro.pi_coding_types[p_vpar->synchro.i_start];
+
+        /* Mean with average tau, to ensure stability. */
+        p_vpar->synchro.p_tau[i_coding_type] =
+            (p_vpar->synchro.pi_meaningful[i_coding_type]
+             * p_vpar->synchro.p_tau[i_coding_type] + tau)
+            / (p_vpar->synchro.pi_meaningful[i_coding_type] + 1);
+        if( p_vpar->synchro.pi_meaningful[i_coding_type] < MAX_PIC_AVERAGE )
+        {
+            p_vpar->synchro.pi_meaningful[i_coding_type]++;
+        }
+#ifdef DEBUG_VPAR
+        intf_DbgMsg("vpar synchro debug: finished decoding %s (%lld)",
+                    i_coding_type == B_CODING_TYPE ? "B" :
+                    (i_coding_type == P_CODING_TYPE ? "P" : "I"), tau);
+#endif
+    }
 
-    p_vpar->synchro.i_fifo_start = (p_vpar->synchro.i_fifo_start + 1) & 0xf;
+    FIFO_INCREMENT( i_start );
 
+    vlc_mutex_unlock( &p_vpar->synchro.fifo_lock );
 }
 
 /*****************************************************************************
@@ -371,24 +411,205 @@ void vpar_SynchroEnd( vpar_thread_t * p_vpar )
  *****************************************************************************/
 mtime_t vpar_SynchroDate( vpar_thread_t * p_vpar )
 {
-    mtime_t i_displaydate = p_vpar->synchro.i_last_display_pts;
-    
-#if 0
-    static mtime_t i_delta = 0;
-
-    fprintf( stderr,
-        "displaying type %i with delay %lli and delta %lli\n",
-        p_vpar->synchro.fifo[p_vpar->synchro.i_fifo_start].i_image_type,
-        i_displaydate - mdate(),
-        i_displaydate - i_delta );
-
-    fprintf (stderr,
-        "theorical fps: %f - actual fps: %f \n",
-        p_vpar->synchro.theorical_fps, p_vpar->synchro.actual_fps );
-
-    i_delta = i_displaydate;
+    /* No need to lock, since PTS are only used by the video parser. */
+    return( p_vpar->synchro.current_pts );
+}
+
+/*****************************************************************************
+ * vpar_SynchroNewPicture: Update stream structure and PTS
+ *****************************************************************************/
+void vpar_SynchroNewPicture( vpar_thread_t * p_vpar, int i_coding_type,
+                             int i_repeat_field )
+{
+    mtime_t         period = 1000000 * 1001 / p_vpar->sequence.i_frame_rate;
+
+    switch( i_coding_type )
+    {
+    case I_CODING_TYPE:
+        if( p_vpar->synchro.i_eta_p
+                && p_vpar->synchro.i_eta_p != p_vpar->synchro.i_n_p )
+        {
+            intf_WarnMsg( 1, "Stream periodicity changed from P[%d] to P[%d]",
+                          p_vpar->synchro.i_n_p, p_vpar->synchro.i_eta_p );
+            p_vpar->synchro.i_n_p = p_vpar->synchro.i_eta_p;
+        }
+        p_vpar->synchro.i_eta_p = p_vpar->synchro.i_eta_b = 0;
+#ifdef STATS
+        if( p_vpar->synchro.i_type == VPAR_SYNCHRO_DEFAULT )
+        {
+            intf_Msg( "vpar synchro stats: I(%lld) P(%lld)[%d] B(%lld)[%d] YUV(%lld) : trashed %d:%d/%d",
+                  p_vpar->synchro.p_tau[I_CODING_TYPE],
+                  p_vpar->synchro.p_tau[P_CODING_TYPE],
+                  p_vpar->synchro.i_n_p,
+                  p_vpar->synchro.p_tau[B_CODING_TYPE],
+                  p_vpar->synchro.i_n_b,
+                  p_vpar->p_vout->render_time,
+                  p_vpar->synchro.i_not_chosen_pic,
+                  p_vpar->synchro.i_trashed_pic -
+                  p_vpar->synchro.i_not_chosen_pic,
+                  p_vpar->synchro.i_pic );
+            p_vpar->synchro.i_trashed_pic = p_vpar->synchro.i_not_chosen_pic
+                = p_vpar->synchro.i_pic = 0;
+        }
 #endif
+        break;
+    case P_CODING_TYPE:
+        p_vpar->synchro.i_eta_p++;
+        if( p_vpar->synchro.i_eta_b
+                && p_vpar->synchro.i_eta_b != p_vpar->synchro.i_n_b )
+        {
+            intf_WarnMsg( 1, "Stream periodicity changed from B[%d] to B[%d]",
+                          p_vpar->synchro.i_n_b, p_vpar->synchro.i_eta_b );
+            p_vpar->synchro.i_n_b = p_vpar->synchro.i_eta_b;
+        }
+        p_vpar->synchro.i_eta_b = 0;
+        break;
+    case B_CODING_TYPE:
+        p_vpar->synchro.i_eta_b++;
+        break;
+    }
+
+    p_vpar->synchro.current_pts += p_vpar->synchro.i_current_period
+                                        * (period >> 1);
+
+#define PTS_THRESHOLD   (period >> 2)
+    if( i_coding_type == B_CODING_TYPE )
+    {
+        /* A video frame can be displayed 1, 2 or 3 times, according to
+         * repeat_first_field, top_field_first, progressive_sequence and
+         * progressive_frame. */
+        p_vpar->synchro.i_current_period = i_repeat_field;
+
+        if( p_vpar->sequence.next_pts )
+        {
+            if( p_vpar->sequence.next_pts - p_vpar->synchro.current_pts
+                    > PTS_THRESHOLD
+                 || p_vpar->synchro.current_pts - p_vpar->sequence.next_pts
+                    > PTS_THRESHOLD )
+            {
+                intf_WarnMsg( 2,
+                        "vpar synchro warning: pts != current_date (%lld)",
+                        p_vpar->synchro.current_pts
+                            - p_vpar->sequence.next_pts );
+            }
+            p_vpar->synchro.current_pts = p_vpar->sequence.next_pts;
+            p_vpar->sequence.next_pts = 0;
+        }
+    }
+    else
+    {
+        p_vpar->synchro.i_current_period = p_vpar->synchro.i_backward_period;
+        p_vpar->synchro.i_backward_period = i_repeat_field;
+
+        if( p_vpar->synchro.backward_pts )
+        {
+            if( p_vpar->sequence.next_dts && 
+                (p_vpar->sequence.next_dts - p_vpar->synchro.backward_pts
+                    > PTS_THRESHOLD
+              || p_vpar->synchro.backward_pts - p_vpar->sequence.next_dts
+                    > PTS_THRESHOLD) )
+            {
+                intf_WarnMsg( 2,
+                        "vpar synchro warning: backward_pts != dts (%lld)",
+                        p_vpar->sequence.next_dts
+                            - p_vpar->synchro.backward_pts );
+            }
+            if( p_vpar->synchro.backward_pts - p_vpar->synchro.current_pts
+                    > PTS_THRESHOLD
+                 || p_vpar->synchro.current_pts - p_vpar->synchro.backward_pts
+                    > PTS_THRESHOLD )
+            {
+                intf_WarnMsg( 2,
+                   "vpar synchro warning: backward_pts != current_pts (%lld)",
+                   p_vpar->synchro.current_pts - p_vpar->synchro.backward_pts );
+            }
+            p_vpar->synchro.current_pts = p_vpar->synchro.backward_pts;
+            p_vpar->synchro.backward_pts = 0;
+        }
+        else if( p_vpar->sequence.next_dts )
+        {
+            if( p_vpar->sequence.next_dts - p_vpar->synchro.current_pts
+                    > PTS_THRESHOLD
+                 || p_vpar->synchro.current_pts - p_vpar->sequence.next_dts
+                    > PTS_THRESHOLD )
+            {
+                intf_WarnMsg( 2,
+                        "vpar synchro warning: dts != current_pts (%lld)",
+                        p_vpar->synchro.current_pts
+                            - p_vpar->sequence.next_dts );
+            }
+            /* By definition of a DTS. */
+            p_vpar->synchro.current_pts = p_vpar->sequence.next_dts;
+            p_vpar->sequence.next_dts = 0;
+        }
+
+        if( p_vpar->sequence.next_pts )
+        {
+            /* Store the PTS for the next time we have to date an I picture. */
+            p_vpar->synchro.backward_pts = p_vpar->sequence.next_pts;
+            p_vpar->sequence.next_pts = 0;
+        }
+    }
+#undef PTS_THRESHOLD
+
+#ifdef STATS
+    p_vpar->synchro.i_pic++;
+#endif
+}
+
+/*****************************************************************************
+ * SynchroType: Get the user's synchro type
+ *****************************************************************************
+ * This function is called at initialization.
+ *****************************************************************************/
+static int SynchroType( void )
+{
+    char * psz_synchro = main_GetPszVariable( VPAR_SYNCHRO_VAR, NULL );
+
+    if( psz_synchro == NULL )
+    {
+        return VPAR_SYNCHRO_DEFAULT;
+    }
+
+    switch( *psz_synchro++ )
+    {
+      case 'i':
+      case 'I':
+        switch( *psz_synchro++ )
+        {
+          case '\0':
+            return VPAR_SYNCHRO_I;
+
+          case '+':
+            if( *psz_synchro ) return 0;
+            return VPAR_SYNCHRO_Iplus;
+
+          case 'p':
+          case 'P':
+            switch( *psz_synchro++ )
+            {
+              case '\0':
+                return VPAR_SYNCHRO_IP;
+
+              case '+':
+                if( *psz_synchro ) return 0;
+                return VPAR_SYNCHRO_IPplus;
+
+              case 'b':
+              case 'B':
+                if( *psz_synchro ) return 0;
+                return VPAR_SYNCHRO_IPB;
+
+              default:
+                return VPAR_SYNCHRO_DEFAULT;
+                
+            }
+
+          default:
+            return VPAR_SYNCHRO_DEFAULT;
+        }
+    }
 
-    return i_displaydate;
+    return VPAR_SYNCHRO_DEFAULT;
 }