]> git.sesse.net Git - ffmpeg/blobdiff - libavcodec/dv.c
* preAlpha DV encoding support -- there's still a truckload
[ffmpeg] / libavcodec / dv.c
index 327149068222a09d21d74da0e3de9a3f076e2872..b9a39965c4c5f2e389943840cc4c16d3ffcfa00c 100644 (file)
@@ -2,6 +2,12 @@
  * DV decoder
  * Copyright (c) 2002 Fabrice Bellard.
  *
+ * DV encoder 
+ * Copyright (c) 2003 Roman Shaposhnik.
+ *
+ * Many thanks to Dan Dennedy <dan@dennedy.org> for providing wealth
+ * of DV technical info.
+ *
  * This library is free software; you can redistribute it and/or
  * modify it under the terms of the GNU Lesser General Public
  * License as published by the Free Software Foundation; either
  * License along with this library; if not, write to the Free Software
  * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
  */
+
+/**
+ * @file dv.c
+ * DV codec.
+ */
 #include "avcodec.h"
 #include "dsputil.h"
 #include "mpegvideo.h"
 #include "simple_idct.h"
-
-#define NTSC_FRAME_SIZE 120000
-#define PAL_FRAME_SIZE  144000
-
-#define TEX_VLC_BITS 9
+#include "dvdata.h"
 
 typedef struct DVVideoDecodeContext {
-    AVCodecContext *avctx;
+    const DVprofile* sys;
     GetBitContext gb;
-    VLC *vlc;
-    int sampling_411; /* 0 = 420, 1 = 411 */
-    int width, height;
-    UINT8 *current_picture[3]; /* picture structure */
-    int linesize[3];
+    AVFrame picture;
     DCTELEM block[5*6][64] __align8;
-    UINT8 dv_zigzag[2][64];
-    UINT8 idct_permutation[64];
+    
+    /* FIXME: the following is extracted from DSP */
+    uint8_t dv_zigzag[2][64];
+    uint8_t idct_permutation[64];
+    void (*get_pixels)(DCTELEM *block, const uint8_t *pixels, int line_size);
+    void (*fdct)(DCTELEM *block);
+    
     /* XXX: move it to static storage ? */
-    UINT8 dv_shift[2][22][64];
-    void (*idct_put[2])(UINT8 *dest, int line_size, DCTELEM *block);
+    uint8_t dv_shift[2][22][64];
+    void (*idct_put[2])(uint8_t *dest, int line_size, DCTELEM *block);
 } DVVideoDecodeContext;
 
-#include "dvdata.h"
-
-static VLC dv_vlc;
+#define TEX_VLC_BITS 9
 /* XXX: also include quantization */
 static RL_VLC_ELEM *dv_rl_vlc[1];
+static VLC_TYPE dv_vlc_codes[15][23];
 
 static void dv_build_unquantize_tables(DVVideoDecodeContext *s)
 {
@@ -71,14 +78,15 @@ static void dv_build_unquantize_tables(DVVideoDecodeContext *s)
     }
 }
 
-static int dvvideo_decode_init(AVCodecContext *avctx)
+static int dvvideo_init(AVCodecContext *avctx)
 {
     DVVideoDecodeContext *s = avctx->priv_data;
     MpegEncContext s2;
-    static int done;
+    static int done=0;
 
     if (!done) {
         int i;
+        VLC dv_vlc;
 
         done = 1;
 
@@ -108,17 +116,27 @@ static int dvvideo_decode_init(AVCodecContext *avctx)
             dv_rl_vlc[0][i].level = level;
             dv_rl_vlc[0][i].run = run;
         }
+
+       memset(dv_vlc_codes, 0xff, sizeof(dv_vlc_codes));
+       for (i = 0; i < NB_DV_VLC - 1; i++) {
+          if (dv_vlc_run[i] < 15 && dv_vlc_level[i] < 23 && dv_vlc_len[i] < 15)
+              dv_vlc_codes[dv_vlc_run[i]][dv_vlc_level[i]] = i;
+       }
     }
 
     /* ugly way to get the idct & scantable */
     /* XXX: fix it */
     memset(&s2, 0, sizeof(MpegEncContext));
     s2.avctx = avctx;
+    dsputil_init(&s2.dsp, avctx);
     if (DCT_common_init(&s2) < 0)
        return -1;
 
-    s->idct_put[0] = s2.idct_put;
-    memcpy(s->idct_permutation, s2.idct_permutation, 64);
+    s->get_pixels = s2.dsp.get_pixels;
+    s->fdct = s2.dsp.fdct;
+    
+    s->idct_put[0] = s2.dsp.idct_put;
+    memcpy(s->idct_permutation, s2.dsp.idct_permutation, 64);
     memcpy(s->dv_zigzag[0], s2.intra_scantable.permutated, 64);
 
     /* XXX: use MMX also for idct248 */
@@ -128,43 +146,50 @@ static int dvvideo_decode_init(AVCodecContext *avctx)
     /* XXX: do it only for constant case */
     dv_build_unquantize_tables(s);
 
+    /* FIXME: I really don't think this should be here */
+    if (dv_codec_profile(avctx))
+       avctx->pix_fmt = dv_codec_profile(avctx)->pix_fmt; 
+    
     return 0;
 }
 
-//#define VLC_DEBUG
+// #define VLC_DEBUG
 
 typedef struct BlockInfo {
-    const UINT8 *shift_table;
-    const UINT8 *scan_table;
-    UINT8 pos; /* position in block */
-    UINT8 eob_reached; /* true if EOB has been reached */
-    UINT8 dct_mode;
-    UINT8 partial_bit_count;
-    UINT16 partial_bit_buffer;
+    const uint8_t *shift_table;
+    const uint8_t *scan_table;
+    uint8_t pos; /* position in block */
+    uint8_t eob_reached; /* true if EOB has been reached */
+    uint8_t dct_mode;
+    uint8_t partial_bit_count;
+    uint16_t partial_bit_buffer;
     int shift_offset;
 } BlockInfo;
 
 /* block size in bits */
-static const UINT16 block_sizes[6] = {
+static const uint16_t block_sizes[6] = {
     112, 112, 112, 112, 80, 80
 };
 
 #ifndef ALT_BITSTREAM_READER
-#error only works with ALT_BITSTREAM_READER
+#warning only works with ALT_BITSTREAM_READER
 #endif
 
 /* decode ac coefs */
 static void dv_decode_ac(DVVideoDecodeContext *s, 
-                         BlockInfo *mb, INT16 *block, int last_index)
+                         BlockInfo *mb, DCTELEM *block, int last_index)
 {
     int last_re_index;
     int shift_offset = mb->shift_offset;
-    const UINT8 *scan_table = mb->scan_table;
-    const UINT8 *shift_table = mb->shift_table;
+    const uint8_t *scan_table = mb->scan_table;
+    const uint8_t *shift_table = mb->shift_table;
     int pos = mb->pos;
     int level, pos1, sign, run;
     int partial_bit_count;
-
+#ifndef ALT_BITSTREAM_READER //FIXME
+    int re_index=0; 
+    int re1_index=0;
+#endif
     OPEN_READER(re, &s->gb);
     
 #ifdef VLC_DEBUG
@@ -174,8 +199,8 @@ static void dv_decode_ac(DVVideoDecodeContext *s,
     /* if we must parse a partial vlc, we do it here */
     partial_bit_count = mb->partial_bit_count;
     if (partial_bit_count > 0) {
-        UINT8 buf[4];
-        UINT32 v;
+        uint8_t buf[4];
+        uint32_t v;
         int l, l1;
         GetBitContext gb1;
 
@@ -193,7 +218,7 @@ static void dv_decode_ac(DVVideoDecodeContext *s,
                v, partial_bit_count, (mb->partial_bit_buffer << l));
 #endif
         /* try to read the codeword */
-        init_get_bits(&gb1, buf, 4);
+        init_get_bits(&gb1, buf, 4*8);
         {
             OPEN_READER(re1, &gb1);
             UPDATE_CACHE(re1, &gb1);
@@ -271,7 +296,7 @@ static void dv_decode_ac(DVVideoDecodeContext *s,
             if (pos >= 64) {
             read_error:
 #if defined(VLC_DEBUG) || 1
-                printf("error pos=%d\n", pos);
+                fprintf(stderr, "error pos=%d\n", pos);
 #endif
                 /* for errors, we consider the eob is reached */
                 mb->eob_reached = 1;
@@ -296,21 +321,21 @@ static inline void bit_copy(PutBitContext *pb, GetBitContext *gb, int bits_left)
 
 /* mb_x and mb_y are in units of 8 pixels */
 static inline void dv_decode_video_segment(DVVideoDecodeContext *s, 
-                                           UINT8 *buf_ptr1, 
-                                           const UINT16 *mb_pos_ptr)
+                                           uint8_t *buf_ptr1, 
+                                           const uint16_t *mb_pos_ptr)
 {
     int quant, dc, dct_mode, class1, j;
     int mb_index, mb_x, mb_y, v, last_index;
     DCTELEM *block, *block1;
     int c_offset, bits_left;
-    UINT8 *y_ptr;
+    uint8_t *y_ptr;
     BlockInfo mb_data[5 * 6], *mb, *mb1;
-    void (*idct_put)(UINT8 *dest, int line_size, DCTELEM *block);
-    UINT8 *buf_ptr;
+    void (*idct_put)(uint8_t *dest, int line_size, DCTELEM *block);
+    uint8_t *buf_ptr;
     PutBitContext pb, vs_pb;
-    UINT8 mb_bit_buffer[80 + 4]; /* allow some slack */
+    uint8_t mb_bit_buffer[80 + 4]; /* allow some slack */
     int mb_bit_count;
-    UINT8 vs_bit_buffer[5 * 80 + 4]; /* allow some slack */
+    uint8_t vs_bit_buffer[5 * 80 + 4]; /* allow some slack */
     int vs_bit_count;
     
     memset(s->block, 0, sizeof(s->block));
@@ -331,7 +356,7 @@ static inline void dv_decode_video_segment(DVVideoDecodeContext *s,
         block = block1;
         for(j = 0;j < 6; j++) {
             /* NOTE: size is not important here */
-            init_get_bits(&s->gb, buf_ptr, 14);
+            init_get_bits(&s->gb, buf_ptr, 14*8);
             
             /* get the dc */
             dc = get_bits(&s->gb, 9);
@@ -353,11 +378,14 @@ static inline void dv_decode_video_segment(DVVideoDecodeContext *s,
             mb->pos = 0;
             mb->partial_bit_count = 0;
 
+#ifdef VLC_DEBUG
+            printf("MB block: %d, %d ", mb_index, j);
+#endif
             dv_decode_ac(s, mb, block, last_index);
 
             /* write the remaining bits  in a new buffer only if the
                block is finished */
-            bits_left = last_index - s->gb.index;
+            bits_left = last_index - get_bits_count(&s->gb);
             if (mb->eob_reached) {
                 mb->partial_bit_count = 0;
                 mb_bit_count += bits_left;
@@ -376,19 +404,19 @@ static inline void dv_decode_video_segment(DVVideoDecodeContext *s,
 
         /* pass 2 : we can do it just after */
 #ifdef VLC_DEBUG
-        printf("***pass 2 size=%d\n", mb_bit_count);
+        printf("***pass 2 size=%d MB#=%d\n", mb_bit_count, mb_index);
 #endif
         block = block1;
         mb = mb1;
-        init_get_bits(&s->gb, mb_bit_buffer, 80);
+        init_get_bits(&s->gb, mb_bit_buffer, 80*8);
         for(j = 0;j < 6; j++) {
-            if (!mb->eob_reached && s->gb.index < mb_bit_count) {
+            if (!mb->eob_reached && get_bits_count(&s->gb) < mb_bit_count) {
                 dv_decode_ac(s, mb, block, mb_bit_count);
                 /* if still not finished, no need to parse other blocks */
                 if (!mb->eob_reached) {
                     /* we could not parse the current AC coefficient,
                        so we add the remaining bytes */
-                    bits_left = mb_bit_count - s->gb.index;
+                    bits_left = mb_bit_count - get_bits_count(&s->gb);
                     if (bits_left > 0) {
                         mb->partial_bit_count += bits_left;
                         mb->partial_bit_buffer = 
@@ -403,7 +431,7 @@ static inline void dv_decode_video_segment(DVVideoDecodeContext *s,
         }
         /* all blocks are finished, so the extra bytes can be used at
            the video segment level */
-        bits_left = mb_bit_count - s->gb.index;
+        bits_left = mb_bit_count - get_bits_count(&s->gb);
         vs_bit_count += bits_left;
         bit_copy(&vs_pb, &s->gb, bits_left);
     next_mb:
@@ -419,7 +447,7 @@ static inline void dv_decode_video_segment(DVVideoDecodeContext *s,
 #endif
     block = &s->block[0][0];
     mb = mb_data;
-    init_get_bits(&s->gb, vs_bit_buffer, 5 * 80);
+    init_get_bits(&s->gb, vs_bit_buffer, 5 * 80*8);
     for(mb_index = 0; mb_index < 5; mb_index++) {
         for(j = 0;j < 6; j++) {
             if (!mb->eob_reached) {
@@ -440,44 +468,44 @@ static inline void dv_decode_video_segment(DVVideoDecodeContext *s,
         v = *mb_pos_ptr++;
         mb_x = v & 0xff;
         mb_y = v >> 8;
-        y_ptr = s->current_picture[0] + (mb_y * s->linesize[0] * 8) + (mb_x * 8);
-        if (s->sampling_411)
-            c_offset = (mb_y * s->linesize[1] * 8) + ((mb_x >> 2) * 8);
+        y_ptr = s->picture.data[0] + (mb_y * s->picture.linesize[0] * 8) + (mb_x * 8);
+        if (s->sys->pix_fmt == PIX_FMT_YUV411P)
+            c_offset = (mb_y * s->picture.linesize[1] * 8) + ((mb_x >> 2) * 8);
         else
-            c_offset = ((mb_y >> 1) * s->linesize[1] * 8) + ((mb_x >> 1) * 8);
+            c_offset = ((mb_y >> 1) * s->picture.linesize[1] * 8) + ((mb_x >> 1) * 8);
         for(j = 0;j < 6; j++) {
             idct_put = s->idct_put[mb->dct_mode];
             if (j < 4) {
-                if (s->sampling_411 && mb_x < (704 / 8)) {
+                if (s->sys->pix_fmt == PIX_FMT_YUV411P && mb_x < (704 / 8)) {
                     /* NOTE: at end of line, the macroblock is handled as 420 */
-                    idct_put(y_ptr + (j * 8), s->linesize[0], block);
+                    idct_put(y_ptr + (j * 8), s->picture.linesize[0], block);
                 } else {
-                    idct_put(y_ptr + ((j & 1) * 8) + ((j >> 1) * 8 * s->linesize[0]),
-                             s->linesize[0], block);
+                    idct_put(y_ptr + ((j & 1) * 8) + ((j >> 1) * 8 * s->picture.linesize[0]),
+                             s->picture.linesize[0], block);
                 }
             } else {
-                if (s->sampling_411 && mb_x >= (704 / 8)) {
+                if (s->sys->pix_fmt == PIX_FMT_YUV411P && mb_x >= (704 / 8)) {
                     uint8_t pixels[64], *c_ptr, *c_ptr1, *ptr;
                     int y, linesize;
                     /* NOTE: at end of line, the macroblock is handled as 420 */
                     idct_put(pixels, 8, block);
-                    linesize = s->linesize[6 - j];
-                    c_ptr = s->current_picture[6 - j] + c_offset;
+                    linesize = s->picture.linesize[6 - j];
+                    c_ptr = s->picture.data[6 - j] + c_offset;
                     ptr = pixels;
                     for(y = 0;y < 8; y++) {
                         /* convert to 411P */
-                        c_ptr1 = c_ptr + linesize;
-                        c_ptr1[0] = c_ptr[0] = (ptr[0] + ptr[1]) >> 1;
-                        c_ptr1[1] = c_ptr[1] = (ptr[2] + ptr[3]) >> 1;
-                        c_ptr1[2] = c_ptr[2] = (ptr[4] + ptr[5]) >> 1;
-                        c_ptr1[3] = c_ptr[3] = (ptr[6] + ptr[7]) >> 1;
-                        c_ptr += linesize * 2;
+                        c_ptr1 = c_ptr + 8*linesize;
+                        c_ptr[0]= ptr[0]; c_ptr1[0]= ptr[4];
+                        c_ptr[1]= ptr[1]; c_ptr1[1]= ptr[5];
+                        c_ptr[2]= ptr[2]; c_ptr1[2]= ptr[6];
+                        c_ptr[3]= ptr[3]; c_ptr1[3]= ptr[7];
+                        c_ptr += linesize;
                         ptr += 8;
                     }
                 } else {
                     /* don't ask me why they inverted Cb and Cr ! */
-                    idct_put(s->current_picture[6 - j] + c_offset, 
-                             s->linesize[6 - j], block);
+                    idct_put(s->picture.data[6 - j] + c_offset, 
+                             s->picture.linesize[6 - j], block);
                 }
             }
             block += 64;
@@ -486,102 +514,356 @@ static inline void dv_decode_video_segment(DVVideoDecodeContext *s,
     }
 }
 
+/* Converts run and level (where level != 0) pair into vlc, returning bit size */
+static inline int dv_rl2vlc(int run, int l, uint32_t* vlc)
+{
+    int sign = l >> 8;
+    int level = (l ^ sign) - sign;
+    int size;
+    
+    sign = (sign & 1);
+
+    if (run < 15 && level < 23 && dv_vlc_codes[run][level] != -1) {
+        *vlc = (dv_vlc_bits[dv_vlc_codes[run][level]] << 1) | sign; 
+       size = dv_vlc_len[dv_vlc_codes[run][level]] + 1;
+    }
+    else { 
+       if (level < 23) {
+           *vlc = (dv_vlc_bits[dv_vlc_codes[0][level]] << 1) | sign; 
+           size = dv_vlc_len[dv_vlc_codes[0][level]] + 1;
+       } else {
+           *vlc = 0xfe00 | (level << 1) | sign;
+           size = 16;
+       }
+
+       switch(run) {
+       case 0:
+           break;
+       case 1:
+       case 2:
+           *vlc |= ((0x7ce | (run - 1)) << size);
+           size += 11;
+           break;
+       case 3:
+       case 4:
+       case 5:
+       case 6:
+           *vlc |= ((0xfac | (run - 3)) << size);
+           size += 12;
+           break;
+       default:
+           *vlc |= ((0x1f80 | (run - 1)) << size);
+           size += 13;
+           break;
+       }
+    }
+    
+    return size;
+}
+
+typedef struct EncBlockInfo {
+    int qno;
+    int cno;
+    int dct_mode;
+    int block_size;
+    DCTELEM *mb;
+    PutBitContext pb;
+} EncBlockInfo;
+
+static inline int dv_bits_left(EncBlockInfo* bi)
+{
+    return (bi->block_size - get_bit_count(&bi->pb));
+}
+
+static inline void dv_encode_ac(EncBlockInfo* bi, PutBitContext* heap)
+{
+    int i, level, size, run = 0;
+    uint32_t vlc;
+    PutBitContext* cpb = &bi->pb;
+    
+    for (i=1; i<64; i++) {
+       level = bi->mb[ff_zigzag_direct[i]] / 
+               (1<<(dv_quant_shifts[bi->qno + dv_quant_offset[bi->cno]]
+                              [dv_88_areas[ff_zigzag_direct[i]]] + 4 + (bi->cno == 3)));
+       if (level != 0) {
+          size = dv_rl2vlc(run, level, &vlc);
+put_vlc:
+
+#ifdef VLC_DEBUG
+           printf(" %3d:%3d", run, level);
+#endif
+          if (cpb == &bi->pb && size > dv_bits_left(bi)) {
+              size -= dv_bits_left(bi);
+              put_bits(cpb, dv_bits_left(bi), vlc >> size);
+              vlc = vlc & ((1<<size)-1);
+              cpb = heap;
+          }
+          put_bits(cpb, size, vlc);
+          run = 0;
+       } else
+          run++;
+    }
+   
+    if (i == 64) {
+        size = 4; vlc = 6; /* End Of Block stamp */
+       goto put_vlc;
+    }
+}
+
+static inline void dv_redistr_bits(EncBlockInfo* bi, int count, uint8_t* extra_data, int extra_bits, PutBitContext* heap)
+{
+    int i;
+    GetBitContext gb;
+    
+    init_get_bits(&gb, extra_data, extra_bits);
+    
+    for (i=0; i<count; i++) {
+       int bits_left = dv_bits_left(bi);
+#ifdef VLC_DEBUG
+       if (bits_left)
+           printf("------------> inserting %d bytes in %d:%d\n", bits_left, i/6, i%6);
+#endif
+       if (bits_left > extra_bits) {
+           bit_copy(&bi->pb, &gb, extra_bits); 
+          extra_bits = 0;
+          break;
+       } else
+           bit_copy(&bi->pb, &gb, bits_left);
+          
+       extra_bits -= bits_left;
+       bi++;
+    }
+    
+    if (extra_bits > 0 && heap)
+       bit_copy(heap, &gb, extra_bits);
+}
+
+static inline void dv_set_class_number(EncBlockInfo* bi, int j)
+{
+    int i, max_ac = 0;
+
+    for (i=1; i<64; i++) {
+       int ac = abs(bi->mb[ff_zigzag_direct[i]]) / 4;
+       if (max_ac < ac)
+           max_ac = ac;
+    }
+    if (max_ac < 12)
+        bi->cno = j;
+    else if (max_ac < 24)
+        bi->cno = j + 1;
+    else if (max_ac < 36)
+        bi->cno = j + 2;
+    else
+        bi->cno = j + 3;
+    
+    if (bi->cno > 3)
+        bi->cno = 3;
+}
+
+/*
+ * This is a very rough initial implementaion. The performance is
+ * horrible and some features are missing, mainly 2-4-8 DCT encoding.
+ * The weighting is missing as well, but it's missing from the decoding
+ * step also -- so at least we're on the same page with decoder ;-)
+ */
+static inline void dv_encode_video_segment(DVVideoDecodeContext *s, 
+                                           uint8_t *dif, 
+                                           const uint16_t *mb_pos_ptr)
+{
+    int mb_index, i, j, v;
+    int mb_x, mb_y, c_offset, linesize; 
+    uint8_t*  y_ptr;
+    uint8_t*  data;
+    int       do_edge_wrap;
+    DCTELEM  *block;
+    EncBlockInfo  enc_blks[5*6];
+    EncBlockInfo* enc_blk;
+    int       free_vs_bits;
+    int extra_bits;
+    PutBitContext extra_vs;
+    uint8_t   extra_vs_data[5*6*128];
+    uint8_t   extra_mb_data[6*128];
+
+    int       QNO = 15;
+   
+    /* Stage 1 -- doing DCT on 5 MBs */
+    block = &s->block[0][0];
+    for(mb_index = 0; mb_index < 5; mb_index++) {
+        v = *mb_pos_ptr++;
+        mb_x = v & 0xff;
+        mb_y = v >> 8;
+        y_ptr = s->picture.data[0] + (mb_y * s->picture.linesize[0] * 8) + (mb_x * 8);
+       c_offset = (s->sys->pix_fmt == PIX_FMT_YUV411P) ?
+                  ((mb_y * s->picture.linesize[1] * 8) + ((mb_x >> 2) * 8)) :
+                  (((mb_y >> 1) * s->picture.linesize[1] * 8) + ((mb_x >> 1) * 8));
+       do_edge_wrap = 0;
+        for(j = 0;j < 6; j++) {
+            if (j < 4) {  /* Four Y blocks */
+               /* NOTE: at end of line, the macroblock is handled as 420 */
+               if (s->sys->pix_fmt == PIX_FMT_YUV411P && mb_x < (704 / 8)) {
+                    data = y_ptr + (j * 8);
+                } else {
+                    data = y_ptr + ((j & 1) * 8) + ((j >> 1) * 8 * s->picture.linesize[0]);
+                }
+               linesize = s->picture.linesize[0];
+            } else {      /* Cr and Cb blocks */
+               /* don't ask Fabrice why they inverted Cb and Cr ! */
+               data = s->picture.data[6 - j] + c_offset;
+               linesize = s->picture.linesize[6 - j];
+               if (s->sys->pix_fmt == PIX_FMT_YUV411P && mb_x >= (704 / 8))
+                   do_edge_wrap = 1;
+           }   
+            
+           /* Everything is set up -- now just copy data -> DCT block */
+           if (do_edge_wrap) {  /* Edge wrap copy: 4x16 -> 8x8 */
+               uint8_t* d;
+               DCTELEM *b = block;
+               for (i=0;i<8;i++) {
+                  d = data + 8 * linesize;
+                  b[0] = data[0]; b[1] = data[1]; b[2] = data[2]; b[3] = data[3];
+                   b[4] =    d[0]; b[5] =    d[1]; b[6] =    d[2]; b[7] =    d[3];
+                  data += linesize;
+                  b += 8;
+               }
+           } else {             /* Simple copy: 8x8 -> 8x8 */
+               s->get_pixels(block, data, linesize);
+           }
+            
+           s->fdct(block);
+           
+           block += 64;
+        }
+    }
+
+    /* Stage 2 -- setup for encoding phase */
+    enc_blk = &enc_blks[0];
+    block = &s->block[0][0];
+    for (i=0; i<5; i++) {
+       for (j=0; j<6; j++) {
+         enc_blk->mb = block;
+         enc_blk->dct_mode = 0;
+         enc_blk->block_size = block_sizes[j];
+         
+         dv_set_class_number(enc_blk, j/4*(j%2));
+         
+         block += 64;
+         enc_blk++;
+       }
+    }
+   
+    /* Stage 3 -- encoding by trial-and-error */
+encode_vs:
+    enc_blk = &enc_blks[0];
+    for (i=0; i<5; i++) {
+       uint8_t* p = dif + i*80 + 4;
+       for (j=0; j<6; j++) {
+          enc_blk->qno = QNO;
+         init_put_bits(&enc_blk->pb, p, block_sizes[j]/8, NULL, NULL);
+         enc_blk++;
+         p += block_sizes[j]/8;
+       }
+    }
+
+    init_put_bits(&extra_vs, extra_vs_data, sizeof(extra_vs_data), NULL, NULL);
+    free_vs_bits = 0;
+    enc_blk = &enc_blks[0];
+    for (i=0; i<5; i++) {
+       PutBitContext extra_mb;
+       EncBlockInfo* enc_blk2 = enc_blk;
+       int free_mb_bits = 0;
+
+       init_put_bits(&extra_mb, extra_mb_data, sizeof(extra_mb_data), NULL, NULL);
+       dif[i*80 + 3] = enc_blk->qno;
+       
+       for (j=0; j<6; j++) {
+         uint16_t dc = ((enc_blk->mb[0] >> 3) - 1024) >> 2;
+
+         put_bits(&enc_blk->pb, 9, dc);
+         put_bits(&enc_blk->pb, 1, enc_blk->dct_mode);
+         put_bits(&enc_blk->pb, 2, enc_blk->cno);
+
+#ifdef VLC_DEBUG
+          printf("[%d, %d]: ", i, j);
+#endif
+         dv_encode_ac(enc_blk, &extra_mb);
+#ifdef VLC_DEBUG
+          printf("\n");
+#endif
+         
+         free_mb_bits += dv_bits_left(enc_blk);
+         enc_blk++;
+       }
+       
+       /* We can't flush extra_mb just yet -- since it'll round up bit number */
+       extra_bits = get_bit_count(&extra_mb);
+       if (free_mb_bits > extra_bits)
+           free_vs_bits += free_mb_bits - extra_bits;
+    
+       if (extra_bits) {  /* FIXME: speed up things when free_mb_bits == 0 */
+           flush_put_bits(&extra_mb);
+           dv_redistr_bits(enc_blk2, 6, extra_mb_data, extra_bits, &extra_vs);
+       }
+    }
+    
+    /* We can't flush extra_mb just yet -- since it'll round up bit number */
+    extra_bits = get_bit_count(&extra_vs);
+    if (extra_bits > free_vs_bits && QNO) { /* FIXME: very crude trial-and-error */
+        QNO--;
+       goto encode_vs;
+    }
+    
+    if (extra_bits) {
+        flush_put_bits(&extra_vs);
+        dv_redistr_bits(&enc_blks[0], 5*6, extra_vs_data, extra_bits, NULL);
+    }
+
+    for (i=0; i<6*5; i++) {
+       flush_put_bits(&enc_blks[i].pb);
+#ifdef VLC_DEBUG
+       printf("[%d:%d] qno=%d cno=%d\n", i/6, i%6, enc_blks[i].qno, enc_blks[i].cno);
+#endif
+    }
+}
 
 /* NOTE: exactly one frame must be given (120000 bytes for NTSC,
    144000 bytes for PAL) */
 static int dvvideo_decode_frame(AVCodecContext *avctx, 
                                  void *data, int *data_size,
-                                 UINT8 *buf, int buf_size)
+                                 uint8_t *buf, int buf_size)
 {
     DVVideoDecodeContext *s = avctx->priv_data;
-    int sct, dsf, apt, ds, nb_dif_segs, vs, size, width, height, i, packet_size;
-    UINT8 *buf_ptr;
-    const UINT16 *mb_pos_ptr;
-    AVPicture *picture;
+    int ds, vs;
+    const uint16_t *mb_pos_ptr;
     
-    /* parse id */
-    init_get_bits(&s->gb, buf, buf_size);
-    sct = get_bits(&s->gb, 3);
-    if (sct != 0)
-        return -1;
-    skip_bits(&s->gb, 5);
-    get_bits(&s->gb, 4); /* dsn (sequence number */
-    get_bits(&s->gb, 1); /* fsc (channel number) */
-    skip_bits(&s->gb, 3);
-    get_bits(&s->gb, 8); /* dbn (diff block number 0-134) */
-
-    dsf = get_bits(&s->gb, 1); /* 0 = NTSC 1 = PAL */
-    if (get_bits(&s->gb, 1) != 0)
-        return -1;
-    skip_bits(&s->gb, 11);
-    apt = get_bits(&s->gb, 3); /* apt */
-
-    get_bits(&s->gb, 1); /* tf1 */
-    skip_bits(&s->gb, 4);
-    get_bits(&s->gb, 3); /* ap1 */
+    s->sys = dv_frame_profile(buf);
+    if (!s->sys || buf_size < s->sys->frame_size)
+        return -1; /* NOTE: we only accept several full frames */
 
-    get_bits(&s->gb, 1); /* tf2 */
-    skip_bits(&s->gb, 4);
-    get_bits(&s->gb, 3); /* ap2 */
-
-    get_bits(&s->gb, 1); /* tf3 */
-    skip_bits(&s->gb, 4);
-    get_bits(&s->gb, 3); /* ap3 */
+       
+    if(s->picture.data[0])
+        avctx->release_buffer(avctx, &s->picture);
     
-    /* init size */
-    width = 720;
-    if (dsf) {
-        packet_size = PAL_FRAME_SIZE;
-        height = 576;
-        nb_dif_segs = 12;
-    } else {
-        packet_size = NTSC_FRAME_SIZE;
-        height = 480;
-        nb_dif_segs = 10;
-    }
-    /* NOTE: we only accept several full frames */
-    if (buf_size < packet_size)
+    s->picture.reference = 0;
+    avctx->pix_fmt = s->sys->pix_fmt;
+    if(avctx->get_buffer(avctx, &s->picture) < 0) {
+        fprintf(stderr, "get_buffer() failed\n");
         return -1;
-    
-    /* XXX: is it correct to assume that 420 is always used in PAL
-       mode ? */
-    s->sampling_411 = !dsf;
-    if (s->sampling_411)
-        mb_pos_ptr = dv_place_411;
-    else
-        mb_pos_ptr = dv_place_420;
-
-    /* (re)alloc picture if needed */
-    if (s->width != width || s->height != height) {
-        for(i=0;i<3;i++)
-            av_freep(&s->current_picture[i]);
-        for(i=0;i<3;i++) {
-            size = width * height;
-            s->linesize[i] = width;
-            if (i >= 1) {
-                size >>= 2;
-                s->linesize[i] >>= s->sampling_411 ? 2 : 1;
-            }
-            s->current_picture[i] = av_malloc(size);
-            if (!s->current_picture[i])
-                return -1;
-        }
-        s->width = width;
-        s->height = height;
     }
 
     /* for each DIF segment */
-    buf_ptr = buf;
-    for (ds = 0; ds < nb_dif_segs; ds++) {
-        buf_ptr += 6 * 80; /* skip DIF segment header */
+    mb_pos_ptr = s->sys->video_place;
+    for (ds = 0; ds < s->sys->difseg_size; ds++) {
+        buf += 6 * 80; /* skip DIF segment header */
         
         for(vs = 0; vs < 27; vs++) {
-            if ((vs % 3) == 0) {
-                /* skip audio block */
-                buf_ptr += 80;
-            }
-            dv_decode_video_segment(s, buf_ptr, mb_pos_ptr);
-            buf_ptr += 5 * 80;
+            if ((vs % 3) == 0)
+               buf += 80; /* skip audio block */
+            
+#ifdef VLC_DEBUG
+            printf("********************* %d, %d **********************\n", ds, vs);
+#endif
+           dv_decode_video_segment(s, buf, mb_pos_ptr);
+            buf += 5 * 80;
             mb_pos_ptr += 5;
         }
     }
@@ -589,85 +871,63 @@ static int dvvideo_decode_frame(AVCodecContext *avctx,
     emms_c();
 
     /* return image */
-    avctx->width = width;
-    avctx->height = height;
-    if (s->sampling_411)
-        avctx->pix_fmt = PIX_FMT_YUV411P;
-    else
-        avctx->pix_fmt = PIX_FMT_YUV420P;
-    if (dsf)
-        avctx->frame_rate = 25 * FRAME_RATE_BASE;
-    else
-        avctx->frame_rate = 30 * FRAME_RATE_BASE;
-    *data_size = sizeof(AVPicture);
-    picture = data;
-    for(i=0;i<3;i++) {
-        picture->data[i] = s->current_picture[i];
-        picture->linesize[i] = s->linesize[i];
-    }
-    return packet_size;
+    *data_size = sizeof(AVFrame);
+    *(AVFrame*)data= s->picture;
+    
+    return s->sys->frame_size;
 }
 
-static int dvvideo_decode_end(AVCodecContext *avctx)
+static int dvvideo_encode_frame(AVCodecContext *c, uint8_t *buf, int buf_size, 
+                                void *data)
 {
-    DVVideoDecodeContext *s = avctx->priv_data;
-    int i;
-
-    for(i=0;i<3;i++)
-        av_freep(&s->current_picture[i]);
-    return 0;
-}
+    DVVideoDecodeContext *s = c->priv_data;
+    const uint16_t *mb_pos_ptr;
+    int ds, vs;
 
-AVCodec dvvideo_decoder = {
-    "dvvideo",
-    CODEC_TYPE_VIDEO,
-    CODEC_ID_DVVIDEO,
-    sizeof(DVVideoDecodeContext),
-    dvvideo_decode_init,
-    NULL,
-    dvvideo_decode_end,
-    dvvideo_decode_frame,
-    0,
-    NULL
-};
-
-typedef struct DVAudioDecodeContext {
-    AVCodecContext *avctx;
-    GetBitContext gb;
+    s->sys = dv_codec_profile(c);
+    if (!s->sys)
+       return -1;
+    
+    c->pix_fmt = s->sys->pix_fmt;
+    s->picture = *((AVFrame *)data);
 
-} DVAudioDecodeContext;
+    /* for each DIF segment */
+    mb_pos_ptr = s->sys->video_place;
+    for (ds = 0; ds < s->sys->difseg_size; ds++) {
+        buf += 6 * 80; /* skip DIF segment header */
+        
+        for(vs = 0; vs < 27; vs++) {
+            if ((vs % 3) == 0)
+               buf += 80; /* skip audio block */
 
-static int dvaudio_decode_init(AVCodecContext *avctx)
-{
-    //    DVAudioDecodeContext *s = avctx->priv_data;
-    return 0;
-}
+#ifdef VLC_DEBUG
+            printf("********************* %d, %d **********************\n", ds, vs);
+#endif
+           dv_encode_video_segment(s, buf, mb_pos_ptr);
+            buf += 5 * 80;
+            mb_pos_ptr += 5;
+        }
+    }
 
-/* NOTE: exactly one frame must be given (120000 bytes for NTSC,
-   144000 bytes for PAL) */
-static int dvaudio_decode_frame(AVCodecContext *avctx, 
-                                 void *data, int *data_size,
-                                 UINT8 *buf, int buf_size)
-{
-    //    DVAudioDecodeContext *s = avctx->priv_data;
-    return buf_size;
+    emms_c();
+    return s->sys->frame_size;
 }
 
-static int dvaudio_decode_end(AVCodecContext *avctx)
+static int dvvideo_end(AVCodecContext *avctx)
 {
-    //    DVAudioDecodeContext *s = avctx->priv_data;
+    avcodec_default_free_buffers(avctx);    
     return 0;
 }
 
-AVCodec dvaudio_decoder = {
-    "dvaudio",
-    CODEC_TYPE_AUDIO,
-    CODEC_ID_DVAUDIO,
-    sizeof(DVAudioDecodeContext),
-    dvaudio_decode_init,
-    NULL,
-    dvaudio_decode_end,
-    dvaudio_decode_frame,
-    0,
+AVCodec dvvideo_decoder = {
+    "dvvideo",
+    CODEC_TYPE_VIDEO,
+    CODEC_ID_DVVIDEO,
+    sizeof(DVVideoDecodeContext),
+    dvvideo_init,
+    dvvideo_encode_frame,
+    dvvideo_end,
+    dvvideo_decode_frame,
+    CODEC_CAP_DR1,
     NULL
 };