Subversion Repositories Kolibri OS

Rev

Blame | Last modification | View Log | RSS feed

  1. /*
  2.  * The simplest mpeg encoder (well, it was the simplest!)
  3.  * Copyright (c) 2000,2001 Fabrice Bellard
  4.  * Copyright (c) 2002-2004 Michael Niedermayer <michaelni@gmx.at>
  5.  *
  6.  * 4MV & hq & B-frame encoding stuff by Michael Niedermayer <michaelni@gmx.at>
  7.  *
  8.  * This file is part of FFmpeg.
  9.  *
  10.  * FFmpeg is free software; you can redistribute it and/or
  11.  * modify it under the terms of the GNU Lesser General Public
  12.  * License as published by the Free Software Foundation; either
  13.  * version 2.1 of the License, or (at your option) any later version.
  14.  *
  15.  * FFmpeg is distributed in the hope that it will be useful,
  16.  * but WITHOUT ANY WARRANTY; without even the implied warranty of
  17.  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  18.  * Lesser General Public License for more details.
  19.  *
  20.  * You should have received a copy of the GNU Lesser General Public
  21.  * License along with FFmpeg; if not, write to the Free Software
  22.  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
  23.  */
  24.  
  25. /**
  26.  * @file
  27.  * The simplest mpeg encoder (well, it was the simplest!).
  28.  */
  29.  
  30. #include "libavutil/attributes.h"
  31. #include "libavutil/avassert.h"
  32. #include "libavutil/imgutils.h"
  33. #include "libavutil/internal.h"
  34. #include "libavutil/motion_vector.h"
  35. #include "libavutil/timer.h"
  36. #include "avcodec.h"
  37. #include "blockdsp.h"
  38. #include "h264chroma.h"
  39. #include "idctdsp.h"
  40. #include "internal.h"
  41. #include "mathops.h"
  42. #include "mpeg_er.h"
  43. #include "mpegutils.h"
  44. #include "mpegvideo.h"
  45. #include "mpegvideodata.h"
  46. #include "mjpegenc.h"
  47. #include "msmpeg4.h"
  48. #include "qpeldsp.h"
  49. #include "thread.h"
  50. #include "wmv2.h"
  51. #include <limits.h>
  52.  
  53. static void dct_unquantize_mpeg1_intra_c(MpegEncContext *s,
  54.                                    int16_t *block, int n, int qscale)
  55. {
  56.     int i, level, nCoeffs;
  57.     const uint16_t *quant_matrix;
  58.  
  59.     nCoeffs= s->block_last_index[n];
  60.  
  61.     block[0] *= n < 4 ? s->y_dc_scale : s->c_dc_scale;
  62.     /* XXX: only mpeg1 */
  63.     quant_matrix = s->intra_matrix;
  64.     for(i=1;i<=nCoeffs;i++) {
  65.         int j= s->intra_scantable.permutated[i];
  66.         level = block[j];
  67.         if (level) {
  68.             if (level < 0) {
  69.                 level = -level;
  70.                 level = (int)(level * qscale * quant_matrix[j]) >> 3;
  71.                 level = (level - 1) | 1;
  72.                 level = -level;
  73.             } else {
  74.                 level = (int)(level * qscale * quant_matrix[j]) >> 3;
  75.                 level = (level - 1) | 1;
  76.             }
  77.             block[j] = level;
  78.         }
  79.     }
  80. }
  81.  
  82. static void dct_unquantize_mpeg1_inter_c(MpegEncContext *s,
  83.                                    int16_t *block, int n, int qscale)
  84. {
  85.     int i, level, nCoeffs;
  86.     const uint16_t *quant_matrix;
  87.  
  88.     nCoeffs= s->block_last_index[n];
  89.  
  90.     quant_matrix = s->inter_matrix;
  91.     for(i=0; i<=nCoeffs; i++) {
  92.         int j= s->intra_scantable.permutated[i];
  93.         level = block[j];
  94.         if (level) {
  95.             if (level < 0) {
  96.                 level = -level;
  97.                 level = (((level << 1) + 1) * qscale *
  98.                          ((int) (quant_matrix[j]))) >> 4;
  99.                 level = (level - 1) | 1;
  100.                 level = -level;
  101.             } else {
  102.                 level = (((level << 1) + 1) * qscale *
  103.                          ((int) (quant_matrix[j]))) >> 4;
  104.                 level = (level - 1) | 1;
  105.             }
  106.             block[j] = level;
  107.         }
  108.     }
  109. }
  110.  
  111. static void dct_unquantize_mpeg2_intra_c(MpegEncContext *s,
  112.                                    int16_t *block, int n, int qscale)
  113. {
  114.     int i, level, nCoeffs;
  115.     const uint16_t *quant_matrix;
  116.  
  117.     if(s->alternate_scan) nCoeffs= 63;
  118.     else nCoeffs= s->block_last_index[n];
  119.  
  120.     block[0] *= n < 4 ? s->y_dc_scale : s->c_dc_scale;
  121.     quant_matrix = s->intra_matrix;
  122.     for(i=1;i<=nCoeffs;i++) {
  123.         int j= s->intra_scantable.permutated[i];
  124.         level = block[j];
  125.         if (level) {
  126.             if (level < 0) {
  127.                 level = -level;
  128.                 level = (int)(level * qscale * quant_matrix[j]) >> 3;
  129.                 level = -level;
  130.             } else {
  131.                 level = (int)(level * qscale * quant_matrix[j]) >> 3;
  132.             }
  133.             block[j] = level;
  134.         }
  135.     }
  136. }
  137.  
  138. static void dct_unquantize_mpeg2_intra_bitexact(MpegEncContext *s,
  139.                                    int16_t *block, int n, int qscale)
  140. {
  141.     int i, level, nCoeffs;
  142.     const uint16_t *quant_matrix;
  143.     int sum=-1;
  144.  
  145.     if(s->alternate_scan) nCoeffs= 63;
  146.     else nCoeffs= s->block_last_index[n];
  147.  
  148.     block[0] *= n < 4 ? s->y_dc_scale : s->c_dc_scale;
  149.     sum += block[0];
  150.     quant_matrix = s->intra_matrix;
  151.     for(i=1;i<=nCoeffs;i++) {
  152.         int j= s->intra_scantable.permutated[i];
  153.         level = block[j];
  154.         if (level) {
  155.             if (level < 0) {
  156.                 level = -level;
  157.                 level = (int)(level * qscale * quant_matrix[j]) >> 3;
  158.                 level = -level;
  159.             } else {
  160.                 level = (int)(level * qscale * quant_matrix[j]) >> 3;
  161.             }
  162.             block[j] = level;
  163.             sum+=level;
  164.         }
  165.     }
  166.     block[63]^=sum&1;
  167. }
  168.  
  169. static void dct_unquantize_mpeg2_inter_c(MpegEncContext *s,
  170.                                    int16_t *block, int n, int qscale)
  171. {
  172.     int i, level, nCoeffs;
  173.     const uint16_t *quant_matrix;
  174.     int sum=-1;
  175.  
  176.     if(s->alternate_scan) nCoeffs= 63;
  177.     else nCoeffs= s->block_last_index[n];
  178.  
  179.     quant_matrix = s->inter_matrix;
  180.     for(i=0; i<=nCoeffs; i++) {
  181.         int j= s->intra_scantable.permutated[i];
  182.         level = block[j];
  183.         if (level) {
  184.             if (level < 0) {
  185.                 level = -level;
  186.                 level = (((level << 1) + 1) * qscale *
  187.                          ((int) (quant_matrix[j]))) >> 4;
  188.                 level = -level;
  189.             } else {
  190.                 level = (((level << 1) + 1) * qscale *
  191.                          ((int) (quant_matrix[j]))) >> 4;
  192.             }
  193.             block[j] = level;
  194.             sum+=level;
  195.         }
  196.     }
  197.     block[63]^=sum&1;
  198. }
  199.  
  200. static void dct_unquantize_h263_intra_c(MpegEncContext *s,
  201.                                   int16_t *block, int n, int qscale)
  202. {
  203.     int i, level, qmul, qadd;
  204.     int nCoeffs;
  205.  
  206.     av_assert2(s->block_last_index[n]>=0 || s->h263_aic);
  207.  
  208.     qmul = qscale << 1;
  209.  
  210.     if (!s->h263_aic) {
  211.         block[0] *= n < 4 ? s->y_dc_scale : s->c_dc_scale;
  212.         qadd = (qscale - 1) | 1;
  213.     }else{
  214.         qadd = 0;
  215.     }
  216.     if(s->ac_pred)
  217.         nCoeffs=63;
  218.     else
  219.         nCoeffs= s->inter_scantable.raster_end[ s->block_last_index[n] ];
  220.  
  221.     for(i=1; i<=nCoeffs; i++) {
  222.         level = block[i];
  223.         if (level) {
  224.             if (level < 0) {
  225.                 level = level * qmul - qadd;
  226.             } else {
  227.                 level = level * qmul + qadd;
  228.             }
  229.             block[i] = level;
  230.         }
  231.     }
  232. }
  233.  
  234. static void dct_unquantize_h263_inter_c(MpegEncContext *s,
  235.                                   int16_t *block, int n, int qscale)
  236. {
  237.     int i, level, qmul, qadd;
  238.     int nCoeffs;
  239.  
  240.     av_assert2(s->block_last_index[n]>=0);
  241.  
  242.     qadd = (qscale - 1) | 1;
  243.     qmul = qscale << 1;
  244.  
  245.     nCoeffs= s->inter_scantable.raster_end[ s->block_last_index[n] ];
  246.  
  247.     for(i=0; i<=nCoeffs; i++) {
  248.         level = block[i];
  249.         if (level) {
  250.             if (level < 0) {
  251.                 level = level * qmul - qadd;
  252.             } else {
  253.                 level = level * qmul + qadd;
  254.             }
  255.             block[i] = level;
  256.         }
  257.     }
  258. }
  259.  
  260.  
  261. static void gray16(uint8_t *dst, const uint8_t *src, ptrdiff_t linesize, int h)
  262. {
  263.     while(h--)
  264.         memset(dst + h*linesize, 128, 16);
  265. }
  266.  
  267. static void gray8(uint8_t *dst, const uint8_t *src, ptrdiff_t linesize, int h)
  268. {
  269.     while(h--)
  270.         memset(dst + h*linesize, 128, 8);
  271. }
  272.  
  273. /* init common dct for both encoder and decoder */
  274. static av_cold int dct_init(MpegEncContext *s)
  275. {
  276.     ff_blockdsp_init(&s->bdsp, s->avctx);
  277.     ff_h264chroma_init(&s->h264chroma, 8); //for lowres
  278.     ff_hpeldsp_init(&s->hdsp, s->avctx->flags);
  279.     ff_mpegvideodsp_init(&s->mdsp);
  280.     ff_videodsp_init(&s->vdsp, s->avctx->bits_per_raw_sample);
  281.  
  282.     if (s->avctx->debug & FF_DEBUG_NOMC) {
  283.         int i;
  284.         for (i=0; i<4; i++) {
  285.             s->hdsp.avg_pixels_tab[0][i] = gray16;
  286.             s->hdsp.put_pixels_tab[0][i] = gray16;
  287.             s->hdsp.put_no_rnd_pixels_tab[0][i] = gray16;
  288.  
  289.             s->hdsp.avg_pixels_tab[1][i] = gray8;
  290.             s->hdsp.put_pixels_tab[1][i] = gray8;
  291.             s->hdsp.put_no_rnd_pixels_tab[1][i] = gray8;
  292.         }
  293.     }
  294.  
  295.     s->dct_unquantize_h263_intra = dct_unquantize_h263_intra_c;
  296.     s->dct_unquantize_h263_inter = dct_unquantize_h263_inter_c;
  297.     s->dct_unquantize_mpeg1_intra = dct_unquantize_mpeg1_intra_c;
  298.     s->dct_unquantize_mpeg1_inter = dct_unquantize_mpeg1_inter_c;
  299.     s->dct_unquantize_mpeg2_intra = dct_unquantize_mpeg2_intra_c;
  300.     if (s->avctx->flags & AV_CODEC_FLAG_BITEXACT)
  301.         s->dct_unquantize_mpeg2_intra = dct_unquantize_mpeg2_intra_bitexact;
  302.     s->dct_unquantize_mpeg2_inter = dct_unquantize_mpeg2_inter_c;
  303.  
  304.     if (HAVE_INTRINSICS_NEON)
  305.         ff_mpv_common_init_neon(s);
  306.  
  307.     if (ARCH_ALPHA)
  308.         ff_mpv_common_init_axp(s);
  309.     if (ARCH_ARM)
  310.         ff_mpv_common_init_arm(s);
  311.     if (ARCH_PPC)
  312.         ff_mpv_common_init_ppc(s);
  313.     if (ARCH_X86)
  314.         ff_mpv_common_init_x86(s);
  315.     if (ARCH_MIPS)
  316.         ff_mpv_common_init_mips(s);
  317.  
  318.     return 0;
  319. }
  320.  
  321. av_cold void ff_mpv_idct_init(MpegEncContext *s)
  322. {
  323.     ff_idctdsp_init(&s->idsp, s->avctx);
  324.  
  325.     /* load & permutate scantables
  326.      * note: only wmv uses different ones
  327.      */
  328.     if (s->alternate_scan) {
  329.         ff_init_scantable(s->idsp.idct_permutation, &s->inter_scantable, ff_alternate_vertical_scan);
  330.         ff_init_scantable(s->idsp.idct_permutation, &s->intra_scantable, ff_alternate_vertical_scan);
  331.     } else {
  332.         ff_init_scantable(s->idsp.idct_permutation, &s->inter_scantable, ff_zigzag_direct);
  333.         ff_init_scantable(s->idsp.idct_permutation, &s->intra_scantable, ff_zigzag_direct);
  334.     }
  335.     ff_init_scantable(s->idsp.idct_permutation, &s->intra_h_scantable, ff_alternate_horizontal_scan);
  336.     ff_init_scantable(s->idsp.idct_permutation, &s->intra_v_scantable, ff_alternate_vertical_scan);
  337. }
  338.  
  339. static int alloc_picture(MpegEncContext *s, Picture *pic, int shared)
  340. {
  341.     return ff_alloc_picture(s->avctx, pic, &s->me, &s->sc, shared, 0,
  342.                             s->chroma_x_shift, s->chroma_y_shift, s->out_format,
  343.                             s->mb_stride, s->mb_width, s->mb_height, s->b8_stride,
  344.                             &s->linesize, &s->uvlinesize);
  345. }
  346.  
  347. static int init_duplicate_context(MpegEncContext *s)
  348. {
  349.     int y_size = s->b8_stride * (2 * s->mb_height + 1);
  350.     int c_size = s->mb_stride * (s->mb_height + 1);
  351.     int yc_size = y_size + 2 * c_size;
  352.     int i;
  353.  
  354.     if (s->mb_height & 1)
  355.         yc_size += 2*s->b8_stride + 2*s->mb_stride;
  356.  
  357.     s->sc.edge_emu_buffer =
  358.     s->me.scratchpad   =
  359.     s->me.temp         =
  360.     s->sc.rd_scratchpad   =
  361.     s->sc.b_scratchpad    =
  362.     s->sc.obmc_scratchpad = NULL;
  363.  
  364.     if (s->encoding) {
  365.         FF_ALLOCZ_OR_GOTO(s->avctx, s->me.map,
  366.                           ME_MAP_SIZE * sizeof(uint32_t), fail)
  367.         FF_ALLOCZ_OR_GOTO(s->avctx, s->me.score_map,
  368.                           ME_MAP_SIZE * sizeof(uint32_t), fail)
  369.         if (s->avctx->noise_reduction) {
  370.             FF_ALLOCZ_OR_GOTO(s->avctx, s->dct_error_sum,
  371.                               2 * 64 * sizeof(int), fail)
  372.         }
  373.     }
  374.     FF_ALLOCZ_OR_GOTO(s->avctx, s->blocks, 64 * 12 * 2 * sizeof(int16_t), fail)
  375.     s->block = s->blocks[0];
  376.  
  377.     for (i = 0; i < 12; i++) {
  378.         s->pblocks[i] = &s->block[i];
  379.     }
  380.     if (s->avctx->codec_tag == AV_RL32("VCR2")) {
  381.         // exchange uv
  382.         FFSWAP(void *, s->pblocks[4], s->pblocks[5]);
  383.     }
  384.  
  385.     if (s->out_format == FMT_H263) {
  386.         /* ac values */
  387.         FF_ALLOCZ_OR_GOTO(s->avctx, s->ac_val_base,
  388.                           yc_size * sizeof(int16_t) * 16, fail);
  389.         s->ac_val[0] = s->ac_val_base + s->b8_stride + 1;
  390.         s->ac_val[1] = s->ac_val_base + y_size + s->mb_stride + 1;
  391.         s->ac_val[2] = s->ac_val[1] + c_size;
  392.     }
  393.  
  394.     return 0;
  395. fail:
  396.     return -1; // free() through ff_mpv_common_end()
  397. }
  398.  
  399. static void free_duplicate_context(MpegEncContext *s)
  400. {
  401.     if (!s)
  402.         return;
  403.  
  404.     av_freep(&s->sc.edge_emu_buffer);
  405.     av_freep(&s->me.scratchpad);
  406.     s->me.temp =
  407.     s->sc.rd_scratchpad =
  408.     s->sc.b_scratchpad =
  409.     s->sc.obmc_scratchpad = NULL;
  410.  
  411.     av_freep(&s->dct_error_sum);
  412.     av_freep(&s->me.map);
  413.     av_freep(&s->me.score_map);
  414.     av_freep(&s->blocks);
  415.     av_freep(&s->ac_val_base);
  416.     s->block = NULL;
  417. }
  418.  
  419. static void backup_duplicate_context(MpegEncContext *bak, MpegEncContext *src)
  420. {
  421. #define COPY(a) bak->a = src->a
  422.     COPY(sc.edge_emu_buffer);
  423.     COPY(me.scratchpad);
  424.     COPY(me.temp);
  425.     COPY(sc.rd_scratchpad);
  426.     COPY(sc.b_scratchpad);
  427.     COPY(sc.obmc_scratchpad);
  428.     COPY(me.map);
  429.     COPY(me.score_map);
  430.     COPY(blocks);
  431.     COPY(block);
  432.     COPY(start_mb_y);
  433.     COPY(end_mb_y);
  434.     COPY(me.map_generation);
  435.     COPY(pb);
  436.     COPY(dct_error_sum);
  437.     COPY(dct_count[0]);
  438.     COPY(dct_count[1]);
  439.     COPY(ac_val_base);
  440.     COPY(ac_val[0]);
  441.     COPY(ac_val[1]);
  442.     COPY(ac_val[2]);
  443. #undef COPY
  444. }
  445.  
  446. int ff_update_duplicate_context(MpegEncContext *dst, MpegEncContext *src)
  447. {
  448.     MpegEncContext bak;
  449.     int i, ret;
  450.     // FIXME copy only needed parts
  451.     // START_TIMER
  452.     backup_duplicate_context(&bak, dst);
  453.     memcpy(dst, src, sizeof(MpegEncContext));
  454.     backup_duplicate_context(dst, &bak);
  455.     for (i = 0; i < 12; i++) {
  456.         dst->pblocks[i] = &dst->block[i];
  457.     }
  458.     if (dst->avctx->codec_tag == AV_RL32("VCR2")) {
  459.         // exchange uv
  460.         FFSWAP(void *, dst->pblocks[4], dst->pblocks[5]);
  461.     }
  462.     if (!dst->sc.edge_emu_buffer &&
  463.         (ret = ff_mpeg_framesize_alloc(dst->avctx, &dst->me,
  464.                                        &dst->sc, dst->linesize)) < 0) {
  465.         av_log(dst->avctx, AV_LOG_ERROR, "failed to allocate context "
  466.                "scratch buffers.\n");
  467.         return ret;
  468.     }
  469.     // STOP_TIMER("update_duplicate_context")
  470.     // about 10k cycles / 0.01 sec for  1000frames on 1ghz with 2 threads
  471.     return 0;
  472. }
  473.  
  474. int ff_mpeg_update_thread_context(AVCodecContext *dst,
  475.                                   const AVCodecContext *src)
  476. {
  477.     int i, ret;
  478.     MpegEncContext *s = dst->priv_data, *s1 = src->priv_data;
  479.  
  480.     if (dst == src)
  481.         return 0;
  482.  
  483.     av_assert0(s != s1);
  484.  
  485.     // FIXME can parameters change on I-frames?
  486.     // in that case dst may need a reinit
  487.     if (!s->context_initialized) {
  488.         int err;
  489.         memcpy(s, s1, sizeof(MpegEncContext));
  490.  
  491.         s->avctx                 = dst;
  492.         s->bitstream_buffer      = NULL;
  493.         s->bitstream_buffer_size = s->allocated_bitstream_buffer_size = 0;
  494.  
  495.         if (s1->context_initialized){
  496. //             s->picture_range_start  += MAX_PICTURE_COUNT;
  497. //             s->picture_range_end    += MAX_PICTURE_COUNT;
  498.             ff_mpv_idct_init(s);
  499.             if((err = ff_mpv_common_init(s)) < 0){
  500.                 memset(s, 0, sizeof(MpegEncContext));
  501.                 s->avctx = dst;
  502.                 return err;
  503.             }
  504.         }
  505.     }
  506.  
  507.     if (s->height != s1->height || s->width != s1->width || s->context_reinit) {
  508.         s->context_reinit = 0;
  509.         s->height = s1->height;
  510.         s->width  = s1->width;
  511.         if ((ret = ff_mpv_common_frame_size_change(s)) < 0)
  512.             return ret;
  513.     }
  514.  
  515.     s->avctx->coded_height  = s1->avctx->coded_height;
  516.     s->avctx->coded_width   = s1->avctx->coded_width;
  517.     s->avctx->width         = s1->avctx->width;
  518.     s->avctx->height        = s1->avctx->height;
  519.  
  520.     s->coded_picture_number = s1->coded_picture_number;
  521.     s->picture_number       = s1->picture_number;
  522.  
  523.     av_assert0(!s->picture || s->picture != s1->picture);
  524.     if(s->picture)
  525.     for (i = 0; i < MAX_PICTURE_COUNT; i++) {
  526.         ff_mpeg_unref_picture(s->avctx, &s->picture[i]);
  527.         if (s1->picture && s1->picture[i].f->buf[0] &&
  528.             (ret = ff_mpeg_ref_picture(s->avctx, &s->picture[i], &s1->picture[i])) < 0)
  529.             return ret;
  530.     }
  531.  
  532. #define UPDATE_PICTURE(pic)\
  533. do {\
  534.     ff_mpeg_unref_picture(s->avctx, &s->pic);\
  535.     if (s1->pic.f && s1->pic.f->buf[0])\
  536.         ret = ff_mpeg_ref_picture(s->avctx, &s->pic, &s1->pic);\
  537.     else\
  538.         ret = ff_update_picture_tables(&s->pic, &s1->pic);\
  539.     if (ret < 0)\
  540.         return ret;\
  541. } while (0)
  542.  
  543.     UPDATE_PICTURE(current_picture);
  544.     UPDATE_PICTURE(last_picture);
  545.     UPDATE_PICTURE(next_picture);
  546.  
  547. #define REBASE_PICTURE(pic, new_ctx, old_ctx)                                 \
  548.     ((pic && pic >= old_ctx->picture &&                                       \
  549.       pic < old_ctx->picture + MAX_PICTURE_COUNT) ?                           \
  550.         &new_ctx->picture[pic - old_ctx->picture] : NULL)
  551.  
  552.     s->last_picture_ptr    = REBASE_PICTURE(s1->last_picture_ptr,    s, s1);
  553.     s->current_picture_ptr = REBASE_PICTURE(s1->current_picture_ptr, s, s1);
  554.     s->next_picture_ptr    = REBASE_PICTURE(s1->next_picture_ptr,    s, s1);
  555.  
  556.     // Error/bug resilience
  557.     s->next_p_frame_damaged = s1->next_p_frame_damaged;
  558.     s->workaround_bugs      = s1->workaround_bugs;
  559.     s->padding_bug_score    = s1->padding_bug_score;
  560.  
  561.     // MPEG4 timing info
  562.     memcpy(&s->last_time_base, &s1->last_time_base,
  563.            (char *) &s1->pb_field_time + sizeof(s1->pb_field_time) -
  564.            (char *) &s1->last_time_base);
  565.  
  566.     // B-frame info
  567.     s->max_b_frames = s1->max_b_frames;
  568.     s->low_delay    = s1->low_delay;
  569.     s->droppable    = s1->droppable;
  570.  
  571.     // DivX handling (doesn't work)
  572.     s->divx_packed  = s1->divx_packed;
  573.  
  574.     if (s1->bitstream_buffer) {
  575.         if (s1->bitstream_buffer_size +
  576.             AV_INPUT_BUFFER_PADDING_SIZE > s->allocated_bitstream_buffer_size) {
  577.             av_fast_malloc(&s->bitstream_buffer,
  578.                            &s->allocated_bitstream_buffer_size,
  579.                            s1->allocated_bitstream_buffer_size);
  580.             if (!s->bitstream_buffer) {
  581.                 s->bitstream_buffer_size = 0;
  582.                 return AVERROR(ENOMEM);
  583.             }
  584.         }
  585.         s->bitstream_buffer_size = s1->bitstream_buffer_size;
  586.         memcpy(s->bitstream_buffer, s1->bitstream_buffer,
  587.                s1->bitstream_buffer_size);
  588.         memset(s->bitstream_buffer + s->bitstream_buffer_size, 0,
  589.                AV_INPUT_BUFFER_PADDING_SIZE);
  590.     }
  591.  
  592.     // linesize dependend scratch buffer allocation
  593.     if (!s->sc.edge_emu_buffer)
  594.         if (s1->linesize) {
  595.             if (ff_mpeg_framesize_alloc(s->avctx, &s->me,
  596.                                         &s->sc, s1->linesize) < 0) {
  597.                 av_log(s->avctx, AV_LOG_ERROR, "Failed to allocate context "
  598.                        "scratch buffers.\n");
  599.                 return AVERROR(ENOMEM);
  600.             }
  601.         } else {
  602.             av_log(s->avctx, AV_LOG_ERROR, "Context scratch buffers could not "
  603.                    "be allocated due to unknown size.\n");
  604.         }
  605.  
  606.     // MPEG2/interlacing info
  607.     memcpy(&s->progressive_sequence, &s1->progressive_sequence,
  608.            (char *) &s1->rtp_mode - (char *) &s1->progressive_sequence);
  609.  
  610.     if (!s1->first_field) {
  611.         s->last_pict_type = s1->pict_type;
  612.         if (s1->current_picture_ptr)
  613.             s->last_lambda_for[s1->pict_type] = s1->current_picture_ptr->f->quality;
  614.     }
  615.  
  616.     return 0;
  617. }
  618.  
  619. /**
  620.  * Set the given MpegEncContext to common defaults
  621.  * (same for encoding and decoding).
  622.  * The changed fields will not depend upon the
  623.  * prior state of the MpegEncContext.
  624.  */
  625. void ff_mpv_common_defaults(MpegEncContext *s)
  626. {
  627.     s->y_dc_scale_table      =
  628.     s->c_dc_scale_table      = ff_mpeg1_dc_scale_table;
  629.     s->chroma_qscale_table   = ff_default_chroma_qscale_table;
  630.     s->progressive_frame     = 1;
  631.     s->progressive_sequence  = 1;
  632.     s->picture_structure     = PICT_FRAME;
  633.  
  634.     s->coded_picture_number  = 0;
  635.     s->picture_number        = 0;
  636.  
  637.     s->f_code                = 1;
  638.     s->b_code                = 1;
  639.  
  640.     s->slice_context_count   = 1;
  641. }
  642.  
  643. /**
  644.  * Set the given MpegEncContext to defaults for decoding.
  645.  * the changed fields will not depend upon
  646.  * the prior state of the MpegEncContext.
  647.  */
  648. void ff_mpv_decode_defaults(MpegEncContext *s)
  649. {
  650.     ff_mpv_common_defaults(s);
  651. }
  652.  
  653. void ff_mpv_decode_init(MpegEncContext *s, AVCodecContext *avctx)
  654. {
  655.     s->avctx           = avctx;
  656.     s->width           = avctx->coded_width;
  657.     s->height          = avctx->coded_height;
  658.     s->codec_id        = avctx->codec->id;
  659.     s->workaround_bugs = avctx->workaround_bugs;
  660.  
  661.     /* convert fourcc to upper case */
  662.     s->codec_tag          = avpriv_toupper4(avctx->codec_tag);
  663. }
  664.  
  665. /**
  666.  * Initialize and allocates MpegEncContext fields dependent on the resolution.
  667.  */
  668. static int init_context_frame(MpegEncContext *s)
  669. {
  670.     int y_size, c_size, yc_size, i, mb_array_size, mv_table_size, x, y;
  671.  
  672.     s->mb_width   = (s->width + 15) / 16;
  673.     s->mb_stride  = s->mb_width + 1;
  674.     s->b8_stride  = s->mb_width * 2 + 1;
  675.     mb_array_size = s->mb_height * s->mb_stride;
  676.     mv_table_size = (s->mb_height + 2) * s->mb_stride + 1;
  677.  
  678.     /* set default edge pos, will be overridden
  679.      * in decode_header if needed */
  680.     s->h_edge_pos = s->mb_width * 16;
  681.     s->v_edge_pos = s->mb_height * 16;
  682.  
  683.     s->mb_num     = s->mb_width * s->mb_height;
  684.  
  685.     s->block_wrap[0] =
  686.     s->block_wrap[1] =
  687.     s->block_wrap[2] =
  688.     s->block_wrap[3] = s->b8_stride;
  689.     s->block_wrap[4] =
  690.     s->block_wrap[5] = s->mb_stride;
  691.  
  692.     y_size  = s->b8_stride * (2 * s->mb_height + 1);
  693.     c_size  = s->mb_stride * (s->mb_height + 1);
  694.     yc_size = y_size + 2   * c_size;
  695.  
  696.     if (s->mb_height & 1)
  697.         yc_size += 2*s->b8_stride + 2*s->mb_stride;
  698.  
  699.     FF_ALLOCZ_OR_GOTO(s->avctx, s->mb_index2xy, (s->mb_num + 1) * sizeof(int), fail); // error ressilience code looks cleaner with this
  700.     for (y = 0; y < s->mb_height; y++)
  701.         for (x = 0; x < s->mb_width; x++)
  702.             s->mb_index2xy[x + y * s->mb_width] = x + y * s->mb_stride;
  703.  
  704.     s->mb_index2xy[s->mb_height * s->mb_width] = (s->mb_height - 1) * s->mb_stride + s->mb_width; // FIXME really needed?
  705.  
  706.     if (s->encoding) {
  707.         /* Allocate MV tables */
  708.         FF_ALLOCZ_OR_GOTO(s->avctx, s->p_mv_table_base,                 mv_table_size * 2 * sizeof(int16_t), fail)
  709.         FF_ALLOCZ_OR_GOTO(s->avctx, s->b_forw_mv_table_base,            mv_table_size * 2 * sizeof(int16_t), fail)
  710.         FF_ALLOCZ_OR_GOTO(s->avctx, s->b_back_mv_table_base,            mv_table_size * 2 * sizeof(int16_t), fail)
  711.         FF_ALLOCZ_OR_GOTO(s->avctx, s->b_bidir_forw_mv_table_base,      mv_table_size * 2 * sizeof(int16_t), fail)
  712.         FF_ALLOCZ_OR_GOTO(s->avctx, s->b_bidir_back_mv_table_base,      mv_table_size * 2 * sizeof(int16_t), fail)
  713.         FF_ALLOCZ_OR_GOTO(s->avctx, s->b_direct_mv_table_base,          mv_table_size * 2 * sizeof(int16_t), fail)
  714.         s->p_mv_table            = s->p_mv_table_base + s->mb_stride + 1;
  715.         s->b_forw_mv_table       = s->b_forw_mv_table_base + s->mb_stride + 1;
  716.         s->b_back_mv_table       = s->b_back_mv_table_base + s->mb_stride + 1;
  717.         s->b_bidir_forw_mv_table = s->b_bidir_forw_mv_table_base + s->mb_stride + 1;
  718.         s->b_bidir_back_mv_table = s->b_bidir_back_mv_table_base + s->mb_stride + 1;
  719.         s->b_direct_mv_table     = s->b_direct_mv_table_base + s->mb_stride + 1;
  720.  
  721.         /* Allocate MB type table */
  722.         FF_ALLOCZ_OR_GOTO(s->avctx, s->mb_type, mb_array_size * sizeof(uint16_t), fail) // needed for encoding
  723.  
  724.         FF_ALLOCZ_OR_GOTO(s->avctx, s->lambda_table, mb_array_size * sizeof(int), fail)
  725.  
  726.         FF_ALLOC_OR_GOTO(s->avctx, s->cplx_tab,
  727.                          mb_array_size * sizeof(float), fail);
  728.         FF_ALLOC_OR_GOTO(s->avctx, s->bits_tab,
  729.                          mb_array_size * sizeof(float), fail);
  730.  
  731.     }
  732.  
  733.     if (s->codec_id == AV_CODEC_ID_MPEG4 ||
  734.         (s->avctx->flags & AV_CODEC_FLAG_INTERLACED_ME)) {
  735.         /* interlaced direct mode decoding tables */
  736.         for (i = 0; i < 2; i++) {
  737.             int j, k;
  738.             for (j = 0; j < 2; j++) {
  739.                 for (k = 0; k < 2; k++) {
  740.                     FF_ALLOCZ_OR_GOTO(s->avctx,
  741.                                       s->b_field_mv_table_base[i][j][k],
  742.                                       mv_table_size * 2 * sizeof(int16_t),
  743.                                       fail);
  744.                     s->b_field_mv_table[i][j][k] = s->b_field_mv_table_base[i][j][k] +
  745.                                                    s->mb_stride + 1;
  746.                 }
  747.                 FF_ALLOCZ_OR_GOTO(s->avctx, s->b_field_select_table [i][j], mb_array_size * 2 * sizeof(uint8_t), fail)
  748.                 FF_ALLOCZ_OR_GOTO(s->avctx, s->p_field_mv_table_base[i][j], mv_table_size * 2 * sizeof(int16_t), fail)
  749.                 s->p_field_mv_table[i][j] = s->p_field_mv_table_base[i][j] + s->mb_stride + 1;
  750.             }
  751.             FF_ALLOCZ_OR_GOTO(s->avctx, s->p_field_select_table[i], mb_array_size * 2 * sizeof(uint8_t), fail)
  752.         }
  753.     }
  754.     if (s->out_format == FMT_H263) {
  755.         /* cbp values */
  756.         FF_ALLOCZ_OR_GOTO(s->avctx, s->coded_block_base, y_size + (s->mb_height&1)*2*s->b8_stride, fail);
  757.         s->coded_block = s->coded_block_base + s->b8_stride + 1;
  758.  
  759.         /* cbp, ac_pred, pred_dir */
  760.         FF_ALLOCZ_OR_GOTO(s->avctx, s->cbp_table     , mb_array_size * sizeof(uint8_t), fail);
  761.         FF_ALLOCZ_OR_GOTO(s->avctx, s->pred_dir_table, mb_array_size * sizeof(uint8_t), fail);
  762.     }
  763.  
  764.     if (s->h263_pred || s->h263_plus || !s->encoding) {
  765.         /* dc values */
  766.         // MN: we need these for  error resilience of intra-frames
  767.         FF_ALLOCZ_OR_GOTO(s->avctx, s->dc_val_base, yc_size * sizeof(int16_t), fail);
  768.         s->dc_val[0] = s->dc_val_base + s->b8_stride + 1;
  769.         s->dc_val[1] = s->dc_val_base + y_size + s->mb_stride + 1;
  770.         s->dc_val[2] = s->dc_val[1] + c_size;
  771.         for (i = 0; i < yc_size; i++)
  772.             s->dc_val_base[i] = 1024;
  773.     }
  774.  
  775.     /* which mb is a intra block */
  776.     FF_ALLOCZ_OR_GOTO(s->avctx, s->mbintra_table, mb_array_size, fail);
  777.     memset(s->mbintra_table, 1, mb_array_size);
  778.  
  779.     /* init macroblock skip table */
  780.     FF_ALLOCZ_OR_GOTO(s->avctx, s->mbskip_table, mb_array_size + 2, fail);
  781.     // Note the + 1 is for  a quicker mpeg4 slice_end detection
  782.  
  783.     return ff_mpeg_er_init(s);
  784. fail:
  785.     return AVERROR(ENOMEM);
  786. }
  787.  
  788. static void clear_context(MpegEncContext *s)
  789. {
  790.     int i, j, k;
  791.  
  792.     memset(&s->next_picture, 0, sizeof(s->next_picture));
  793.     memset(&s->last_picture, 0, sizeof(s->last_picture));
  794.     memset(&s->current_picture, 0, sizeof(s->current_picture));
  795.     memset(&s->new_picture, 0, sizeof(s->new_picture));
  796.  
  797.     memset(s->thread_context, 0, sizeof(s->thread_context));
  798.  
  799.     s->me.map = NULL;
  800.     s->me.score_map = NULL;
  801.     s->dct_error_sum = NULL;
  802.     s->block = NULL;
  803.     s->blocks = NULL;
  804.     memset(s->pblocks, 0, sizeof(s->pblocks));
  805.     s->ac_val_base = NULL;
  806.     s->ac_val[0] =
  807.     s->ac_val[1] =
  808.     s->ac_val[2] =NULL;
  809.     s->sc.edge_emu_buffer = NULL;
  810.     s->me.scratchpad = NULL;
  811.     s->me.temp =
  812.     s->sc.rd_scratchpad =
  813.     s->sc.b_scratchpad =
  814.     s->sc.obmc_scratchpad = NULL;
  815.  
  816.     s->parse_context.buffer = NULL;
  817.     s->parse_context.buffer_size = 0;
  818.     s->parse_context.overread = 0;
  819.     s->bitstream_buffer = NULL;
  820.     s->allocated_bitstream_buffer_size = 0;
  821.     s->picture          = NULL;
  822.     s->mb_type          = NULL;
  823.     s->p_mv_table_base  = NULL;
  824.     s->b_forw_mv_table_base = NULL;
  825.     s->b_back_mv_table_base = NULL;
  826.     s->b_bidir_forw_mv_table_base = NULL;
  827.     s->b_bidir_back_mv_table_base = NULL;
  828.     s->b_direct_mv_table_base = NULL;
  829.     s->p_mv_table            = NULL;
  830.     s->b_forw_mv_table       = NULL;
  831.     s->b_back_mv_table       = NULL;
  832.     s->b_bidir_forw_mv_table = NULL;
  833.     s->b_bidir_back_mv_table = NULL;
  834.     s->b_direct_mv_table     = NULL;
  835.     for (i = 0; i < 2; i++) {
  836.         for (j = 0; j < 2; j++) {
  837.             for (k = 0; k < 2; k++) {
  838.                 s->b_field_mv_table_base[i][j][k] = NULL;
  839.                 s->b_field_mv_table[i][j][k] = NULL;
  840.             }
  841.             s->b_field_select_table[i][j] = NULL;
  842.             s->p_field_mv_table_base[i][j] = NULL;
  843.             s->p_field_mv_table[i][j] = NULL;
  844.         }
  845.         s->p_field_select_table[i] = NULL;
  846.     }
  847.  
  848.     s->dc_val_base = NULL;
  849.     s->coded_block_base = NULL;
  850.     s->mbintra_table = NULL;
  851.     s->cbp_table = NULL;
  852.     s->pred_dir_table = NULL;
  853.  
  854.     s->mbskip_table = NULL;
  855.  
  856.     s->er.error_status_table = NULL;
  857.     s->er.er_temp_buffer = NULL;
  858.     s->mb_index2xy = NULL;
  859.     s->lambda_table = NULL;
  860.  
  861.     s->cplx_tab = NULL;
  862.     s->bits_tab = NULL;
  863. }
  864.  
  865. /**
  866.  * init common structure for both encoder and decoder.
  867.  * this assumes that some variables like width/height are already set
  868.  */
  869. av_cold int ff_mpv_common_init(MpegEncContext *s)
  870. {
  871.     int i;
  872.     int nb_slices = (HAVE_THREADS &&
  873.                      s->avctx->active_thread_type & FF_THREAD_SLICE) ?
  874.                     s->avctx->thread_count : 1;
  875.  
  876.     clear_context(s);
  877.  
  878.     if (s->encoding && s->avctx->slices)
  879.         nb_slices = s->avctx->slices;
  880.  
  881.     if (s->codec_id == AV_CODEC_ID_MPEG2VIDEO && !s->progressive_sequence)
  882.         s->mb_height = (s->height + 31) / 32 * 2;
  883.     else
  884.         s->mb_height = (s->height + 15) / 16;
  885.  
  886.     if (s->avctx->pix_fmt == AV_PIX_FMT_NONE) {
  887.         av_log(s->avctx, AV_LOG_ERROR,
  888.                "decoding to AV_PIX_FMT_NONE is not supported.\n");
  889.         return -1;
  890.     }
  891.  
  892.     if (nb_slices > MAX_THREADS || (nb_slices > s->mb_height && s->mb_height)) {
  893.         int max_slices;
  894.         if (s->mb_height)
  895.             max_slices = FFMIN(MAX_THREADS, s->mb_height);
  896.         else
  897.             max_slices = MAX_THREADS;
  898.         av_log(s->avctx, AV_LOG_WARNING, "too many threads/slices (%d),"
  899.                " reducing to %d\n", nb_slices, max_slices);
  900.         nb_slices = max_slices;
  901.     }
  902.  
  903.     if ((s->width || s->height) &&
  904.         av_image_check_size(s->width, s->height, 0, s->avctx))
  905.         return -1;
  906.  
  907.     dct_init(s);
  908.  
  909.     /* set chroma shifts */
  910.     avcodec_get_chroma_sub_sample(s->avctx->pix_fmt,
  911.                                   &s->chroma_x_shift,
  912.                                   &s->chroma_y_shift);
  913.  
  914.  
  915.     FF_ALLOCZ_OR_GOTO(s->avctx, s->picture,
  916.                       MAX_PICTURE_COUNT * sizeof(Picture), fail);
  917.     for (i = 0; i < MAX_PICTURE_COUNT; i++) {
  918.         s->picture[i].f = av_frame_alloc();
  919.         if (!s->picture[i].f)
  920.             goto fail;
  921.     }
  922.     s->next_picture.f = av_frame_alloc();
  923.     if (!s->next_picture.f)
  924.         goto fail;
  925.     s->last_picture.f = av_frame_alloc();
  926.     if (!s->last_picture.f)
  927.         goto fail;
  928.     s->current_picture.f = av_frame_alloc();
  929.     if (!s->current_picture.f)
  930.         goto fail;
  931.     s->new_picture.f = av_frame_alloc();
  932.     if (!s->new_picture.f)
  933.         goto fail;
  934.  
  935.         if (init_context_frame(s))
  936.             goto fail;
  937.  
  938.         s->parse_context.state = -1;
  939.  
  940.         s->context_initialized = 1;
  941.         memset(s->thread_context, 0, sizeof(s->thread_context));
  942.         s->thread_context[0]   = s;
  943.  
  944. //     if (s->width && s->height) {
  945.         if (nb_slices > 1) {
  946.             for (i = 0; i < nb_slices; i++) {
  947.                 if (i) {
  948.                     s->thread_context[i] = av_memdup(s, sizeof(MpegEncContext));
  949.                     if (!s->thread_context[i])
  950.                         goto fail;
  951.                 }
  952.                 if (init_duplicate_context(s->thread_context[i]) < 0)
  953.                     goto fail;
  954.                     s->thread_context[i]->start_mb_y =
  955.                         (s->mb_height * (i) + nb_slices / 2) / nb_slices;
  956.                     s->thread_context[i]->end_mb_y   =
  957.                         (s->mb_height * (i + 1) + nb_slices / 2) / nb_slices;
  958.             }
  959.         } else {
  960.             if (init_duplicate_context(s) < 0)
  961.                 goto fail;
  962.             s->start_mb_y = 0;
  963.             s->end_mb_y   = s->mb_height;
  964.         }
  965.         s->slice_context_count = nb_slices;
  966. //     }
  967.  
  968.     return 0;
  969.  fail:
  970.     ff_mpv_common_end(s);
  971.     return -1;
  972. }
  973.  
  974. /**
  975.  * Frees and resets MpegEncContext fields depending on the resolution.
  976.  * Is used during resolution changes to avoid a full reinitialization of the
  977.  * codec.
  978.  */
  979. static void free_context_frame(MpegEncContext *s)
  980. {
  981.     int i, j, k;
  982.  
  983.     av_freep(&s->mb_type);
  984.     av_freep(&s->p_mv_table_base);
  985.     av_freep(&s->b_forw_mv_table_base);
  986.     av_freep(&s->b_back_mv_table_base);
  987.     av_freep(&s->b_bidir_forw_mv_table_base);
  988.     av_freep(&s->b_bidir_back_mv_table_base);
  989.     av_freep(&s->b_direct_mv_table_base);
  990.     s->p_mv_table            = NULL;
  991.     s->b_forw_mv_table       = NULL;
  992.     s->b_back_mv_table       = NULL;
  993.     s->b_bidir_forw_mv_table = NULL;
  994.     s->b_bidir_back_mv_table = NULL;
  995.     s->b_direct_mv_table     = NULL;
  996.     for (i = 0; i < 2; i++) {
  997.         for (j = 0; j < 2; j++) {
  998.             for (k = 0; k < 2; k++) {
  999.                 av_freep(&s->b_field_mv_table_base[i][j][k]);
  1000.                 s->b_field_mv_table[i][j][k] = NULL;
  1001.             }
  1002.             av_freep(&s->b_field_select_table[i][j]);
  1003.             av_freep(&s->p_field_mv_table_base[i][j]);
  1004.             s->p_field_mv_table[i][j] = NULL;
  1005.         }
  1006.         av_freep(&s->p_field_select_table[i]);
  1007.     }
  1008.  
  1009.     av_freep(&s->dc_val_base);
  1010.     av_freep(&s->coded_block_base);
  1011.     av_freep(&s->mbintra_table);
  1012.     av_freep(&s->cbp_table);
  1013.     av_freep(&s->pred_dir_table);
  1014.  
  1015.     av_freep(&s->mbskip_table);
  1016.  
  1017.     av_freep(&s->er.error_status_table);
  1018.     av_freep(&s->er.er_temp_buffer);
  1019.     av_freep(&s->mb_index2xy);
  1020.     av_freep(&s->lambda_table);
  1021.  
  1022.     av_freep(&s->cplx_tab);
  1023.     av_freep(&s->bits_tab);
  1024.  
  1025.     s->linesize = s->uvlinesize = 0;
  1026. }
  1027.  
  1028. int ff_mpv_common_frame_size_change(MpegEncContext *s)
  1029. {
  1030.     int i, err = 0;
  1031.  
  1032.     if (!s->context_initialized)
  1033.         return AVERROR(EINVAL);
  1034.  
  1035.     if (s->slice_context_count > 1) {
  1036.         for (i = 0; i < s->slice_context_count; i++) {
  1037.             free_duplicate_context(s->thread_context[i]);
  1038.         }
  1039.         for (i = 1; i < s->slice_context_count; i++) {
  1040.             av_freep(&s->thread_context[i]);
  1041.         }
  1042.     } else
  1043.         free_duplicate_context(s);
  1044.  
  1045.     free_context_frame(s);
  1046.  
  1047.     if (s->picture)
  1048.         for (i = 0; i < MAX_PICTURE_COUNT; i++) {
  1049.                 s->picture[i].needs_realloc = 1;
  1050.         }
  1051.  
  1052.     s->last_picture_ptr         =
  1053.     s->next_picture_ptr         =
  1054.     s->current_picture_ptr      = NULL;
  1055.  
  1056.     // init
  1057.     if (s->codec_id == AV_CODEC_ID_MPEG2VIDEO && !s->progressive_sequence)
  1058.         s->mb_height = (s->height + 31) / 32 * 2;
  1059.     else
  1060.         s->mb_height = (s->height + 15) / 16;
  1061.  
  1062.     if ((s->width || s->height) &&
  1063.         (err = av_image_check_size(s->width, s->height, 0, s->avctx)) < 0)
  1064.         goto fail;
  1065.  
  1066.     if ((err = init_context_frame(s)))
  1067.         goto fail;
  1068.  
  1069.     memset(s->thread_context, 0, sizeof(s->thread_context));
  1070.     s->thread_context[0]   = s;
  1071.  
  1072.     if (s->width && s->height) {
  1073.         int nb_slices = s->slice_context_count;
  1074.         if (nb_slices > 1) {
  1075.             for (i = 0; i < nb_slices; i++) {
  1076.                 if (i) {
  1077.                     s->thread_context[i] = av_memdup(s, sizeof(MpegEncContext));
  1078.                     if (!s->thread_context[i]) {
  1079.                         err = AVERROR(ENOMEM);
  1080.                         goto fail;
  1081.                     }
  1082.                 }
  1083.                 if ((err = init_duplicate_context(s->thread_context[i])) < 0)
  1084.                     goto fail;
  1085.                     s->thread_context[i]->start_mb_y =
  1086.                         (s->mb_height * (i) + nb_slices / 2) / nb_slices;
  1087.                     s->thread_context[i]->end_mb_y   =
  1088.                         (s->mb_height * (i + 1) + nb_slices / 2) / nb_slices;
  1089.             }
  1090.         } else {
  1091.             err = init_duplicate_context(s);
  1092.             if (err < 0)
  1093.                 goto fail;
  1094.             s->start_mb_y = 0;
  1095.             s->end_mb_y   = s->mb_height;
  1096.         }
  1097.         s->slice_context_count = nb_slices;
  1098.     }
  1099.  
  1100.     return 0;
  1101.  fail:
  1102.     ff_mpv_common_end(s);
  1103.     return err;
  1104. }
  1105.  
  1106. /* init common structure for both encoder and decoder */
  1107. void ff_mpv_common_end(MpegEncContext *s)
  1108. {
  1109.     int i;
  1110.  
  1111.     if (!s)
  1112.         return ;
  1113.  
  1114.     if (s->slice_context_count > 1) {
  1115.         for (i = 0; i < s->slice_context_count; i++) {
  1116.             free_duplicate_context(s->thread_context[i]);
  1117.         }
  1118.         for (i = 1; i < s->slice_context_count; i++) {
  1119.             av_freep(&s->thread_context[i]);
  1120.         }
  1121.         s->slice_context_count = 1;
  1122.     } else free_duplicate_context(s);
  1123.  
  1124.     av_freep(&s->parse_context.buffer);
  1125.     s->parse_context.buffer_size = 0;
  1126.  
  1127.     av_freep(&s->bitstream_buffer);
  1128.     s->allocated_bitstream_buffer_size = 0;
  1129.  
  1130.     if (s->picture) {
  1131.         for (i = 0; i < MAX_PICTURE_COUNT; i++) {
  1132.             ff_free_picture_tables(&s->picture[i]);
  1133.             ff_mpeg_unref_picture(s->avctx, &s->picture[i]);
  1134.             av_frame_free(&s->picture[i].f);
  1135.         }
  1136.     }
  1137.     av_freep(&s->picture);
  1138.     ff_free_picture_tables(&s->last_picture);
  1139.     ff_mpeg_unref_picture(s->avctx, &s->last_picture);
  1140.     av_frame_free(&s->last_picture.f);
  1141.     ff_free_picture_tables(&s->current_picture);
  1142.     ff_mpeg_unref_picture(s->avctx, &s->current_picture);
  1143.     av_frame_free(&s->current_picture.f);
  1144.     ff_free_picture_tables(&s->next_picture);
  1145.     ff_mpeg_unref_picture(s->avctx, &s->next_picture);
  1146.     av_frame_free(&s->next_picture.f);
  1147.     ff_free_picture_tables(&s->new_picture);
  1148.     ff_mpeg_unref_picture(s->avctx, &s->new_picture);
  1149.     av_frame_free(&s->new_picture.f);
  1150.  
  1151.     free_context_frame(s);
  1152.  
  1153.     s->context_initialized      = 0;
  1154.     s->last_picture_ptr         =
  1155.     s->next_picture_ptr         =
  1156.     s->current_picture_ptr      = NULL;
  1157.     s->linesize = s->uvlinesize = 0;
  1158. }
  1159.  
  1160.  
  1161. static void gray_frame(AVFrame *frame)
  1162. {
  1163.     int i, h_chroma_shift, v_chroma_shift;
  1164.  
  1165.     av_pix_fmt_get_chroma_sub_sample(frame->format, &h_chroma_shift, &v_chroma_shift);
  1166.  
  1167.     for(i=0; i<frame->height; i++)
  1168.         memset(frame->data[0] + frame->linesize[0]*i, 0x80, frame->width);
  1169.     for(i=0; i<FF_CEIL_RSHIFT(frame->height, v_chroma_shift); i++) {
  1170.         memset(frame->data[1] + frame->linesize[1]*i,
  1171.                0x80, FF_CEIL_RSHIFT(frame->width, h_chroma_shift));
  1172.         memset(frame->data[2] + frame->linesize[2]*i,
  1173.                0x80, FF_CEIL_RSHIFT(frame->width, h_chroma_shift));
  1174.     }
  1175. }
  1176.  
  1177. /**
  1178.  * generic function called after decoding
  1179.  * the header and before a frame is decoded.
  1180.  */
  1181. int ff_mpv_frame_start(MpegEncContext *s, AVCodecContext *avctx)
  1182. {
  1183.     int i, ret;
  1184.     Picture *pic;
  1185.     s->mb_skipped = 0;
  1186.  
  1187.     if (!ff_thread_can_start_frame(avctx)) {
  1188.         av_log(avctx, AV_LOG_ERROR, "Attempt to start a frame outside SETUP state\n");
  1189.         return -1;
  1190.     }
  1191.  
  1192.     /* mark & release old frames */
  1193.     if (s->pict_type != AV_PICTURE_TYPE_B && s->last_picture_ptr &&
  1194.         s->last_picture_ptr != s->next_picture_ptr &&
  1195.         s->last_picture_ptr->f->buf[0]) {
  1196.         ff_mpeg_unref_picture(s->avctx, s->last_picture_ptr);
  1197.     }
  1198.  
  1199.     /* release forgotten pictures */
  1200.     /* if (mpeg124/h263) */
  1201.     for (i = 0; i < MAX_PICTURE_COUNT; i++) {
  1202.         if (&s->picture[i] != s->last_picture_ptr &&
  1203.             &s->picture[i] != s->next_picture_ptr &&
  1204.             s->picture[i].reference && !s->picture[i].needs_realloc) {
  1205.             if (!(avctx->active_thread_type & FF_THREAD_FRAME))
  1206.                 av_log(avctx, AV_LOG_ERROR,
  1207.                        "releasing zombie picture\n");
  1208.             ff_mpeg_unref_picture(s->avctx, &s->picture[i]);
  1209.         }
  1210.     }
  1211.  
  1212.     ff_mpeg_unref_picture(s->avctx, &s->current_picture);
  1213.  
  1214.     /* release non reference frames */
  1215.     for (i = 0; i < MAX_PICTURE_COUNT; i++) {
  1216.         if (!s->picture[i].reference)
  1217.             ff_mpeg_unref_picture(s->avctx, &s->picture[i]);
  1218.     }
  1219.  
  1220.     if (s->current_picture_ptr && !s->current_picture_ptr->f->buf[0]) {
  1221.         // we already have a unused image
  1222.         // (maybe it was set before reading the header)
  1223.         pic = s->current_picture_ptr;
  1224.     } else {
  1225.         i   = ff_find_unused_picture(s->avctx, s->picture, 0);
  1226.         if (i < 0) {
  1227.             av_log(s->avctx, AV_LOG_ERROR, "no frame buffer available\n");
  1228.             return i;
  1229.         }
  1230.         pic = &s->picture[i];
  1231.     }
  1232.  
  1233.     pic->reference = 0;
  1234.     if (!s->droppable) {
  1235.         if (s->pict_type != AV_PICTURE_TYPE_B)
  1236.             pic->reference = 3;
  1237.     }
  1238.  
  1239.     pic->f->coded_picture_number = s->coded_picture_number++;
  1240.  
  1241.     if (alloc_picture(s, pic, 0) < 0)
  1242.         return -1;
  1243.  
  1244.     s->current_picture_ptr = pic;
  1245.     // FIXME use only the vars from current_pic
  1246.     s->current_picture_ptr->f->top_field_first = s->top_field_first;
  1247.     if (s->codec_id == AV_CODEC_ID_MPEG1VIDEO ||
  1248.         s->codec_id == AV_CODEC_ID_MPEG2VIDEO) {
  1249.         if (s->picture_structure != PICT_FRAME)
  1250.             s->current_picture_ptr->f->top_field_first =
  1251.                 (s->picture_structure == PICT_TOP_FIELD) == s->first_field;
  1252.     }
  1253.     s->current_picture_ptr->f->interlaced_frame = !s->progressive_frame &&
  1254.                                                  !s->progressive_sequence;
  1255.     s->current_picture_ptr->field_picture      =  s->picture_structure != PICT_FRAME;
  1256.  
  1257.     s->current_picture_ptr->f->pict_type = s->pict_type;
  1258.     // if (s->avctx->flags && AV_CODEC_FLAG_QSCALE)
  1259.     //     s->current_picture_ptr->quality = s->new_picture_ptr->quality;
  1260.     s->current_picture_ptr->f->key_frame = s->pict_type == AV_PICTURE_TYPE_I;
  1261.  
  1262.     if ((ret = ff_mpeg_ref_picture(s->avctx, &s->current_picture,
  1263.                                    s->current_picture_ptr)) < 0)
  1264.         return ret;
  1265.  
  1266.     if (s->pict_type != AV_PICTURE_TYPE_B) {
  1267.         s->last_picture_ptr = s->next_picture_ptr;
  1268.         if (!s->droppable)
  1269.             s->next_picture_ptr = s->current_picture_ptr;
  1270.     }
  1271.     ff_dlog(s->avctx, "L%p N%p C%p L%p N%p C%p type:%d drop:%d\n",
  1272.             s->last_picture_ptr, s->next_picture_ptr,s->current_picture_ptr,
  1273.             s->last_picture_ptr    ? s->last_picture_ptr->f->data[0]    : NULL,
  1274.             s->next_picture_ptr    ? s->next_picture_ptr->f->data[0]    : NULL,
  1275.             s->current_picture_ptr ? s->current_picture_ptr->f->data[0] : NULL,
  1276.             s->pict_type, s->droppable);
  1277.  
  1278.     if ((!s->last_picture_ptr || !s->last_picture_ptr->f->buf[0]) &&
  1279.         (s->pict_type != AV_PICTURE_TYPE_I ||
  1280.          s->picture_structure != PICT_FRAME)) {
  1281.         int h_chroma_shift, v_chroma_shift;
  1282.         av_pix_fmt_get_chroma_sub_sample(s->avctx->pix_fmt,
  1283.                                          &h_chroma_shift, &v_chroma_shift);
  1284.         if (s->pict_type == AV_PICTURE_TYPE_B && s->next_picture_ptr && s->next_picture_ptr->f->buf[0])
  1285.             av_log(avctx, AV_LOG_DEBUG,
  1286.                    "allocating dummy last picture for B frame\n");
  1287.         else if (s->pict_type != AV_PICTURE_TYPE_I)
  1288.             av_log(avctx, AV_LOG_ERROR,
  1289.                    "warning: first frame is no keyframe\n");
  1290.         else if (s->picture_structure != PICT_FRAME)
  1291.             av_log(avctx, AV_LOG_DEBUG,
  1292.                    "allocate dummy last picture for field based first keyframe\n");
  1293.  
  1294.         /* Allocate a dummy frame */
  1295.         i = ff_find_unused_picture(s->avctx, s->picture, 0);
  1296.         if (i < 0) {
  1297.             av_log(s->avctx, AV_LOG_ERROR, "no frame buffer available\n");
  1298.             return i;
  1299.         }
  1300.         s->last_picture_ptr = &s->picture[i];
  1301.  
  1302.         s->last_picture_ptr->reference   = 3;
  1303.         s->last_picture_ptr->f->key_frame = 0;
  1304.         s->last_picture_ptr->f->pict_type = AV_PICTURE_TYPE_P;
  1305.  
  1306.         if (alloc_picture(s, s->last_picture_ptr, 0) < 0) {
  1307.             s->last_picture_ptr = NULL;
  1308.             return -1;
  1309.         }
  1310.  
  1311.         if (!avctx->hwaccel
  1312. #if FF_API_CAP_VDPAU
  1313.             && !(avctx->codec->capabilities&AV_CODEC_CAP_HWACCEL_VDPAU)
  1314. #endif
  1315.             ) {
  1316.             for(i=0; i<avctx->height; i++)
  1317.                 memset(s->last_picture_ptr->f->data[0] + s->last_picture_ptr->f->linesize[0]*i,
  1318.                        0x80, avctx->width);
  1319.             if (s->last_picture_ptr->f->data[2]) {
  1320.                 for(i=0; i<FF_CEIL_RSHIFT(avctx->height, v_chroma_shift); i++) {
  1321.                     memset(s->last_picture_ptr->f->data[1] + s->last_picture_ptr->f->linesize[1]*i,
  1322.                         0x80, FF_CEIL_RSHIFT(avctx->width, h_chroma_shift));
  1323.                     memset(s->last_picture_ptr->f->data[2] + s->last_picture_ptr->f->linesize[2]*i,
  1324.                         0x80, FF_CEIL_RSHIFT(avctx->width, h_chroma_shift));
  1325.                 }
  1326.             }
  1327.  
  1328.             if(s->codec_id == AV_CODEC_ID_FLV1 || s->codec_id == AV_CODEC_ID_H263){
  1329.                 for(i=0; i<avctx->height; i++)
  1330.                 memset(s->last_picture_ptr->f->data[0] + s->last_picture_ptr->f->linesize[0]*i, 16, avctx->width);
  1331.             }
  1332.         }
  1333.  
  1334.         ff_thread_report_progress(&s->last_picture_ptr->tf, INT_MAX, 0);
  1335.         ff_thread_report_progress(&s->last_picture_ptr->tf, INT_MAX, 1);
  1336.     }
  1337.     if ((!s->next_picture_ptr || !s->next_picture_ptr->f->buf[0]) &&
  1338.         s->pict_type == AV_PICTURE_TYPE_B) {
  1339.         /* Allocate a dummy frame */
  1340.         i = ff_find_unused_picture(s->avctx, s->picture, 0);
  1341.         if (i < 0) {
  1342.             av_log(s->avctx, AV_LOG_ERROR, "no frame buffer available\n");
  1343.             return i;
  1344.         }
  1345.         s->next_picture_ptr = &s->picture[i];
  1346.  
  1347.         s->next_picture_ptr->reference   = 3;
  1348.         s->next_picture_ptr->f->key_frame = 0;
  1349.         s->next_picture_ptr->f->pict_type = AV_PICTURE_TYPE_P;
  1350.  
  1351.         if (alloc_picture(s, s->next_picture_ptr, 0) < 0) {
  1352.             s->next_picture_ptr = NULL;
  1353.             return -1;
  1354.         }
  1355.         ff_thread_report_progress(&s->next_picture_ptr->tf, INT_MAX, 0);
  1356.         ff_thread_report_progress(&s->next_picture_ptr->tf, INT_MAX, 1);
  1357.     }
  1358.  
  1359. #if 0 // BUFREF-FIXME
  1360.     memset(s->last_picture.f->data, 0, sizeof(s->last_picture.f->data));
  1361.     memset(s->next_picture.f->data, 0, sizeof(s->next_picture.f->data));
  1362. #endif
  1363.     if (s->last_picture_ptr) {
  1364.         ff_mpeg_unref_picture(s->avctx, &s->last_picture);
  1365.         if (s->last_picture_ptr->f->buf[0] &&
  1366.             (ret = ff_mpeg_ref_picture(s->avctx, &s->last_picture,
  1367.                                        s->last_picture_ptr)) < 0)
  1368.             return ret;
  1369.     }
  1370.     if (s->next_picture_ptr) {
  1371.         ff_mpeg_unref_picture(s->avctx, &s->next_picture);
  1372.         if (s->next_picture_ptr->f->buf[0] &&
  1373.             (ret = ff_mpeg_ref_picture(s->avctx, &s->next_picture,
  1374.                                        s->next_picture_ptr)) < 0)
  1375.             return ret;
  1376.     }
  1377.  
  1378.     av_assert0(s->pict_type == AV_PICTURE_TYPE_I || (s->last_picture_ptr &&
  1379.                                                  s->last_picture_ptr->f->buf[0]));
  1380.  
  1381.     if (s->picture_structure!= PICT_FRAME) {
  1382.         int i;
  1383.         for (i = 0; i < 4; i++) {
  1384.             if (s->picture_structure == PICT_BOTTOM_FIELD) {
  1385.                 s->current_picture.f->data[i] +=
  1386.                     s->current_picture.f->linesize[i];
  1387.             }
  1388.             s->current_picture.f->linesize[i] *= 2;
  1389.             s->last_picture.f->linesize[i]    *= 2;
  1390.             s->next_picture.f->linesize[i]    *= 2;
  1391.         }
  1392.     }
  1393.  
  1394.     /* set dequantizer, we can't do it during init as
  1395.      * it might change for mpeg4 and we can't do it in the header
  1396.      * decode as init is not called for mpeg4 there yet */
  1397.     if (s->mpeg_quant || s->codec_id == AV_CODEC_ID_MPEG2VIDEO) {
  1398.         s->dct_unquantize_intra = s->dct_unquantize_mpeg2_intra;
  1399.         s->dct_unquantize_inter = s->dct_unquantize_mpeg2_inter;
  1400.     } else if (s->out_format == FMT_H263 || s->out_format == FMT_H261) {
  1401.         s->dct_unquantize_intra = s->dct_unquantize_h263_intra;
  1402.         s->dct_unquantize_inter = s->dct_unquantize_h263_inter;
  1403.     } else {
  1404.         s->dct_unquantize_intra = s->dct_unquantize_mpeg1_intra;
  1405.         s->dct_unquantize_inter = s->dct_unquantize_mpeg1_inter;
  1406.     }
  1407.  
  1408.     if (s->avctx->debug & FF_DEBUG_NOMC) {
  1409.         gray_frame(s->current_picture_ptr->f);
  1410.     }
  1411.  
  1412.     return 0;
  1413. }
  1414.  
  1415. /* called after a frame has been decoded. */
  1416. void ff_mpv_frame_end(MpegEncContext *s)
  1417. {
  1418.     emms_c();
  1419.  
  1420.     if (s->current_picture.reference)
  1421.         ff_thread_report_progress(&s->current_picture_ptr->tf, INT_MAX, 0);
  1422. }
  1423.  
  1424.  
  1425. #if FF_API_VISMV
  1426. static int clip_line(int *sx, int *sy, int *ex, int *ey, int maxx)
  1427. {
  1428.     if(*sx > *ex)
  1429.         return clip_line(ex, ey, sx, sy, maxx);
  1430.  
  1431.     if (*sx < 0) {
  1432.         if (*ex < 0)
  1433.             return 1;
  1434.         *sy = *ey + (*sy - *ey) * (int64_t)*ex / (*ex - *sx);
  1435.         *sx = 0;
  1436.     }
  1437.  
  1438.     if (*ex > maxx) {
  1439.         if (*sx > maxx)
  1440.             return 1;
  1441.         *ey = *sy + (*ey - *sy) * (int64_t)(maxx - *sx) / (*ex - *sx);
  1442.         *ex = maxx;
  1443.     }
  1444.     return 0;
  1445. }
  1446.  
  1447.  
  1448. /**
  1449.  * Draw a line from (ex, ey) -> (sx, sy).
  1450.  * @param w width of the image
  1451.  * @param h height of the image
  1452.  * @param stride stride/linesize of the image
  1453.  * @param color color of the arrow
  1454.  */
  1455. static void draw_line(uint8_t *buf, int sx, int sy, int ex, int ey,
  1456.                       int w, int h, int stride, int color)
  1457. {
  1458.     int x, y, fr, f;
  1459.  
  1460.     if (clip_line(&sx, &sy, &ex, &ey, w - 1))
  1461.         return;
  1462.     if (clip_line(&sy, &sx, &ey, &ex, h - 1))
  1463.         return;
  1464.  
  1465.     sx = av_clip(sx, 0, w - 1);
  1466.     sy = av_clip(sy, 0, h - 1);
  1467.     ex = av_clip(ex, 0, w - 1);
  1468.     ey = av_clip(ey, 0, h - 1);
  1469.  
  1470.     buf[sy * stride + sx] += color;
  1471.  
  1472.     if (FFABS(ex - sx) > FFABS(ey - sy)) {
  1473.         if (sx > ex) {
  1474.             FFSWAP(int, sx, ex);
  1475.             FFSWAP(int, sy, ey);
  1476.         }
  1477.         buf += sx + sy * stride;
  1478.         ex  -= sx;
  1479.         f    = ((ey - sy) << 16) / ex;
  1480.         for (x = 0; x <= ex; x++) {
  1481.             y  = (x * f) >> 16;
  1482.             fr = (x * f) & 0xFFFF;
  1483.             buf[y * stride + x]       += (color * (0x10000 - fr)) >> 16;
  1484.             if(fr) buf[(y + 1) * stride + x] += (color *            fr ) >> 16;
  1485.         }
  1486.     } else {
  1487.         if (sy > ey) {
  1488.             FFSWAP(int, sx, ex);
  1489.             FFSWAP(int, sy, ey);
  1490.         }
  1491.         buf += sx + sy * stride;
  1492.         ey  -= sy;
  1493.         if (ey)
  1494.             f = ((ex - sx) << 16) / ey;
  1495.         else
  1496.             f = 0;
  1497.         for(y= 0; y <= ey; y++){
  1498.             x  = (y*f) >> 16;
  1499.             fr = (y*f) & 0xFFFF;
  1500.             buf[y * stride + x]     += (color * (0x10000 - fr)) >> 16;
  1501.             if(fr) buf[y * stride + x + 1] += (color *            fr ) >> 16;
  1502.         }
  1503.     }
  1504. }
  1505.  
  1506. /**
  1507.  * Draw an arrow from (ex, ey) -> (sx, sy).
  1508.  * @param w width of the image
  1509.  * @param h height of the image
  1510.  * @param stride stride/linesize of the image
  1511.  * @param color color of the arrow
  1512.  */
  1513. static void draw_arrow(uint8_t *buf, int sx, int sy, int ex,
  1514.                        int ey, int w, int h, int stride, int color, int tail, int direction)
  1515. {
  1516.     int dx,dy;
  1517.  
  1518.     if (direction) {
  1519.         FFSWAP(int, sx, ex);
  1520.         FFSWAP(int, sy, ey);
  1521.     }
  1522.  
  1523.     sx = av_clip(sx, -100, w + 100);
  1524.     sy = av_clip(sy, -100, h + 100);
  1525.     ex = av_clip(ex, -100, w + 100);
  1526.     ey = av_clip(ey, -100, h + 100);
  1527.  
  1528.     dx = ex - sx;
  1529.     dy = ey - sy;
  1530.  
  1531.     if (dx * dx + dy * dy > 3 * 3) {
  1532.         int rx =  dx + dy;
  1533.         int ry = -dx + dy;
  1534.         int length = ff_sqrt((rx * rx + ry * ry) << 8);
  1535.  
  1536.         // FIXME subpixel accuracy
  1537.         rx = ROUNDED_DIV(rx * 3 << 4, length);
  1538.         ry = ROUNDED_DIV(ry * 3 << 4, length);
  1539.  
  1540.         if (tail) {
  1541.             rx = -rx;
  1542.             ry = -ry;
  1543.         }
  1544.  
  1545.         draw_line(buf, sx, sy, sx + rx, sy + ry, w, h, stride, color);
  1546.         draw_line(buf, sx, sy, sx - ry, sy + rx, w, h, stride, color);
  1547.     }
  1548.     draw_line(buf, sx, sy, ex, ey, w, h, stride, color);
  1549. }
  1550. #endif
  1551.  
  1552. static int add_mb(AVMotionVector *mb, uint32_t mb_type,
  1553.                   int dst_x, int dst_y,
  1554.                   int src_x, int src_y,
  1555.                   int direction)
  1556. {
  1557.     mb->w = IS_8X8(mb_type) || IS_8X16(mb_type) ? 8 : 16;
  1558.     mb->h = IS_8X8(mb_type) || IS_16X8(mb_type) ? 8 : 16;
  1559.     mb->src_x = src_x;
  1560.     mb->src_y = src_y;
  1561.     mb->dst_x = dst_x;
  1562.     mb->dst_y = dst_y;
  1563.     mb->source = direction ? 1 : -1;
  1564.     mb->flags = 0; // XXX: does mb_type contain extra information that could be exported here?
  1565.     return 1;
  1566. }
  1567.  
  1568. /**
  1569.  * Print debugging info for the given picture.
  1570.  */
  1571. void ff_print_debug_info2(AVCodecContext *avctx, AVFrame *pict, uint8_t *mbskip_table,
  1572.                          uint32_t *mbtype_table, int8_t *qscale_table, int16_t (*motion_val[2])[2],
  1573.                          int *low_delay,
  1574.                          int mb_width, int mb_height, int mb_stride, int quarter_sample)
  1575. {
  1576.     if ((avctx->flags2 & AV_CODEC_FLAG2_EXPORT_MVS) && mbtype_table && motion_val[0]) {
  1577.         const int shift = 1 + quarter_sample;
  1578.         const int mv_sample_log2 = avctx->codec_id == AV_CODEC_ID_H264 || avctx->codec_id == AV_CODEC_ID_SVQ3 ? 2 : 1;
  1579.         const int mv_stride      = (mb_width << mv_sample_log2) +
  1580.                                    (avctx->codec->id == AV_CODEC_ID_H264 ? 0 : 1);
  1581.         int mb_x, mb_y, mbcount = 0;
  1582.  
  1583.         /* size is width * height * 2 * 4 where 2 is for directions and 4 is
  1584.          * for the maximum number of MB (4 MB in case of IS_8x8) */
  1585.         AVMotionVector *mvs = av_malloc_array(mb_width * mb_height, 2 * 4 * sizeof(AVMotionVector));
  1586.         if (!mvs)
  1587.             return;
  1588.  
  1589.         for (mb_y = 0; mb_y < mb_height; mb_y++) {
  1590.             for (mb_x = 0; mb_x < mb_width; mb_x++) {
  1591.                 int i, direction, mb_type = mbtype_table[mb_x + mb_y * mb_stride];
  1592.                 for (direction = 0; direction < 2; direction++) {
  1593.                     if (!USES_LIST(mb_type, direction))
  1594.                         continue;
  1595.                     if (IS_8X8(mb_type)) {
  1596.                         for (i = 0; i < 4; i++) {
  1597.                             int sx = mb_x * 16 + 4 + 8 * (i & 1);
  1598.                             int sy = mb_y * 16 + 4 + 8 * (i >> 1);
  1599.                             int xy = (mb_x * 2 + (i & 1) +
  1600.                                       (mb_y * 2 + (i >> 1)) * mv_stride) << (mv_sample_log2 - 1);
  1601.                             int mx = (motion_val[direction][xy][0] >> shift) + sx;
  1602.                             int my = (motion_val[direction][xy][1] >> shift) + sy;
  1603.                             mbcount += add_mb(mvs + mbcount, mb_type, sx, sy, mx, my, direction);
  1604.                         }
  1605.                     } else if (IS_16X8(mb_type)) {
  1606.                         for (i = 0; i < 2; i++) {
  1607.                             int sx = mb_x * 16 + 8;
  1608.                             int sy = mb_y * 16 + 4 + 8 * i;
  1609.                             int xy = (mb_x * 2 + (mb_y * 2 + i) * mv_stride) << (mv_sample_log2 - 1);
  1610.                             int mx = (motion_val[direction][xy][0] >> shift);
  1611.                             int my = (motion_val[direction][xy][1] >> shift);
  1612.  
  1613.                             if (IS_INTERLACED(mb_type))
  1614.                                 my *= 2;
  1615.  
  1616.                             mbcount += add_mb(mvs + mbcount, mb_type, sx, sy, mx + sx, my + sy, direction);
  1617.                         }
  1618.                     } else if (IS_8X16(mb_type)) {
  1619.                         for (i = 0; i < 2; i++) {
  1620.                             int sx = mb_x * 16 + 4 + 8 * i;
  1621.                             int sy = mb_y * 16 + 8;
  1622.                             int xy = (mb_x * 2 + i + mb_y * 2 * mv_stride) << (mv_sample_log2 - 1);
  1623.                             int mx = motion_val[direction][xy][0] >> shift;
  1624.                             int my = motion_val[direction][xy][1] >> shift;
  1625.  
  1626.                             if (IS_INTERLACED(mb_type))
  1627.                                 my *= 2;
  1628.  
  1629.                             mbcount += add_mb(mvs + mbcount, mb_type, sx, sy, mx + sx, my + sy, direction);
  1630.                         }
  1631.                     } else {
  1632.                           int sx = mb_x * 16 + 8;
  1633.                           int sy = mb_y * 16 + 8;
  1634.                           int xy = (mb_x + mb_y * mv_stride) << mv_sample_log2;
  1635.                           int mx = (motion_val[direction][xy][0]>>shift) + sx;
  1636.                           int my = (motion_val[direction][xy][1]>>shift) + sy;
  1637.                           mbcount += add_mb(mvs + mbcount, mb_type, sx, sy, mx, my, direction);
  1638.                     }
  1639.                 }
  1640.             }
  1641.         }
  1642.  
  1643.         if (mbcount) {
  1644.             AVFrameSideData *sd;
  1645.  
  1646.             av_log(avctx, AV_LOG_DEBUG, "Adding %d MVs info to frame %d\n", mbcount, avctx->frame_number);
  1647.             sd = av_frame_new_side_data(pict, AV_FRAME_DATA_MOTION_VECTORS, mbcount * sizeof(AVMotionVector));
  1648.             if (!sd) {
  1649.                 av_freep(&mvs);
  1650.                 return;
  1651.             }
  1652.             memcpy(sd->data, mvs, mbcount * sizeof(AVMotionVector));
  1653.         }
  1654.  
  1655.         av_freep(&mvs);
  1656.     }
  1657.  
  1658.     /* TODO: export all the following to make them accessible for users (and filters) */
  1659.     if (avctx->hwaccel || !mbtype_table
  1660. #if FF_API_CAP_VDPAU
  1661.         || (avctx->codec->capabilities&AV_CODEC_CAP_HWACCEL_VDPAU)
  1662. #endif
  1663.         )
  1664.         return;
  1665.  
  1666.  
  1667.     if (avctx->debug & (FF_DEBUG_SKIP | FF_DEBUG_QP | FF_DEBUG_MB_TYPE)) {
  1668.         int x,y;
  1669.  
  1670.         av_log(avctx, AV_LOG_DEBUG, "New frame, type: %c\n",
  1671.                av_get_picture_type_char(pict->pict_type));
  1672.         for (y = 0; y < mb_height; y++) {
  1673.             for (x = 0; x < mb_width; x++) {
  1674.                 if (avctx->debug & FF_DEBUG_SKIP) {
  1675.                     int count = mbskip_table ? mbskip_table[x + y * mb_stride] : 0;
  1676.                     if (count > 9)
  1677.                         count = 9;
  1678.                     av_log(avctx, AV_LOG_DEBUG, "%1d", count);
  1679.                 }
  1680.                 if (avctx->debug & FF_DEBUG_QP) {
  1681.                     av_log(avctx, AV_LOG_DEBUG, "%2d",
  1682.                            qscale_table[x + y * mb_stride]);
  1683.                 }
  1684.                 if (avctx->debug & FF_DEBUG_MB_TYPE) {
  1685.                     int mb_type = mbtype_table[x + y * mb_stride];
  1686.                     // Type & MV direction
  1687.                     if (IS_PCM(mb_type))
  1688.                         av_log(avctx, AV_LOG_DEBUG, "P");
  1689.                     else if (IS_INTRA(mb_type) && IS_ACPRED(mb_type))
  1690.                         av_log(avctx, AV_LOG_DEBUG, "A");
  1691.                     else if (IS_INTRA4x4(mb_type))
  1692.                         av_log(avctx, AV_LOG_DEBUG, "i");
  1693.                     else if (IS_INTRA16x16(mb_type))
  1694.                         av_log(avctx, AV_LOG_DEBUG, "I");
  1695.                     else if (IS_DIRECT(mb_type) && IS_SKIP(mb_type))
  1696.                         av_log(avctx, AV_LOG_DEBUG, "d");
  1697.                     else if (IS_DIRECT(mb_type))
  1698.                         av_log(avctx, AV_LOG_DEBUG, "D");
  1699.                     else if (IS_GMC(mb_type) && IS_SKIP(mb_type))
  1700.                         av_log(avctx, AV_LOG_DEBUG, "g");
  1701.                     else if (IS_GMC(mb_type))
  1702.                         av_log(avctx, AV_LOG_DEBUG, "G");
  1703.                     else if (IS_SKIP(mb_type))
  1704.                         av_log(avctx, AV_LOG_DEBUG, "S");
  1705.                     else if (!USES_LIST(mb_type, 1))
  1706.                         av_log(avctx, AV_LOG_DEBUG, ">");
  1707.                     else if (!USES_LIST(mb_type, 0))
  1708.                         av_log(avctx, AV_LOG_DEBUG, "<");
  1709.                     else {
  1710.                         av_assert2(USES_LIST(mb_type, 0) && USES_LIST(mb_type, 1));
  1711.                         av_log(avctx, AV_LOG_DEBUG, "X");
  1712.                     }
  1713.  
  1714.                     // segmentation
  1715.                     if (IS_8X8(mb_type))
  1716.                         av_log(avctx, AV_LOG_DEBUG, "+");
  1717.                     else if (IS_16X8(mb_type))
  1718.                         av_log(avctx, AV_LOG_DEBUG, "-");
  1719.                     else if (IS_8X16(mb_type))
  1720.                         av_log(avctx, AV_LOG_DEBUG, "|");
  1721.                     else if (IS_INTRA(mb_type) || IS_16X16(mb_type))
  1722.                         av_log(avctx, AV_LOG_DEBUG, " ");
  1723.                     else
  1724.                         av_log(avctx, AV_LOG_DEBUG, "?");
  1725.  
  1726.  
  1727.                     if (IS_INTERLACED(mb_type))
  1728.                         av_log(avctx, AV_LOG_DEBUG, "=");
  1729.                     else
  1730.                         av_log(avctx, AV_LOG_DEBUG, " ");
  1731.                 }
  1732.             }
  1733.             av_log(avctx, AV_LOG_DEBUG, "\n");
  1734.         }
  1735.     }
  1736.  
  1737.     if ((avctx->debug & (FF_DEBUG_VIS_QP | FF_DEBUG_VIS_MB_TYPE)) ||
  1738.         (avctx->debug_mv)) {
  1739.         int mb_y;
  1740.         int i;
  1741.         int h_chroma_shift, v_chroma_shift, block_height;
  1742. #if FF_API_VISMV
  1743.         const int shift = 1 + quarter_sample;
  1744.         uint8_t *ptr;
  1745.         const int width          = avctx->width;
  1746.         const int height         = avctx->height;
  1747. #endif
  1748.         const int mv_sample_log2 = avctx->codec_id == AV_CODEC_ID_H264 || avctx->codec_id == AV_CODEC_ID_SVQ3 ? 2 : 1;
  1749.         const int mv_stride      = (mb_width << mv_sample_log2) +
  1750.                                    (avctx->codec->id == AV_CODEC_ID_H264 ? 0 : 1);
  1751.  
  1752.         *low_delay = 0; // needed to see the vectors without trashing the buffers
  1753.  
  1754.         avcodec_get_chroma_sub_sample(avctx->pix_fmt, &h_chroma_shift, &v_chroma_shift);
  1755.  
  1756.         av_frame_make_writable(pict);
  1757.  
  1758.         pict->opaque = NULL;
  1759. #if FF_API_VISMV
  1760.         ptr          = pict->data[0];
  1761. #endif
  1762.         block_height = 16 >> v_chroma_shift;
  1763.  
  1764.         for (mb_y = 0; mb_y < mb_height; mb_y++) {
  1765.             int mb_x;
  1766.             for (mb_x = 0; mb_x < mb_width; mb_x++) {
  1767.                 const int mb_index = mb_x + mb_y * mb_stride;
  1768. #if FF_API_VISMV
  1769.                 if ((avctx->debug_mv) && motion_val[0]) {
  1770.                     int type;
  1771.                     for (type = 0; type < 3; type++) {
  1772.                         int direction = 0;
  1773.                         switch (type) {
  1774.                         case 0:
  1775.                             if ((!(avctx->debug_mv & FF_DEBUG_VIS_MV_P_FOR)) ||
  1776.                                 (pict->pict_type!= AV_PICTURE_TYPE_P))
  1777.                                 continue;
  1778.                             direction = 0;
  1779.                             break;
  1780.                         case 1:
  1781.                             if ((!(avctx->debug_mv & FF_DEBUG_VIS_MV_B_FOR)) ||
  1782.                                 (pict->pict_type!= AV_PICTURE_TYPE_B))
  1783.                                 continue;
  1784.                             direction = 0;
  1785.                             break;
  1786.                         case 2:
  1787.                             if ((!(avctx->debug_mv & FF_DEBUG_VIS_MV_B_BACK)) ||
  1788.                                 (pict->pict_type!= AV_PICTURE_TYPE_B))
  1789.                                 continue;
  1790.                             direction = 1;
  1791.                             break;
  1792.                         }
  1793.                         if (!USES_LIST(mbtype_table[mb_index], direction))
  1794.                             continue;
  1795.  
  1796.                         if (IS_8X8(mbtype_table[mb_index])) {
  1797.                             int i;
  1798.                             for (i = 0; i < 4; i++) {
  1799.                                 int sx = mb_x * 16 + 4 + 8 * (i & 1);
  1800.                                 int sy = mb_y * 16 + 4 + 8 * (i >> 1);
  1801.                                 int xy = (mb_x * 2 + (i & 1) +
  1802.                                           (mb_y * 2 + (i >> 1)) * mv_stride) << (mv_sample_log2 - 1);
  1803.                                 int mx = (motion_val[direction][xy][0] >> shift) + sx;
  1804.                                 int my = (motion_val[direction][xy][1] >> shift) + sy;
  1805.                                 draw_arrow(ptr, sx, sy, mx, my, width,
  1806.                                            height, pict->linesize[0], 100, 0, direction);
  1807.                             }
  1808.                         } else if (IS_16X8(mbtype_table[mb_index])) {
  1809.                             int i;
  1810.                             for (i = 0; i < 2; i++) {
  1811.                                 int sx = mb_x * 16 + 8;
  1812.                                 int sy = mb_y * 16 + 4 + 8 * i;
  1813.                                 int xy = (mb_x * 2 + (mb_y * 2 + i) * mv_stride) << (mv_sample_log2 - 1);
  1814.                                 int mx = (motion_val[direction][xy][0] >> shift);
  1815.                                 int my = (motion_val[direction][xy][1] >> shift);
  1816.  
  1817.                                 if (IS_INTERLACED(mbtype_table[mb_index]))
  1818.                                     my *= 2;
  1819.  
  1820.                                 draw_arrow(ptr, sx, sy, mx + sx, my + sy, width,
  1821.                                            height, pict->linesize[0], 100, 0, direction);
  1822.                             }
  1823.                         } else if (IS_8X16(mbtype_table[mb_index])) {
  1824.                             int i;
  1825.                             for (i = 0; i < 2; i++) {
  1826.                                 int sx = mb_x * 16 + 4 + 8 * i;
  1827.                                 int sy = mb_y * 16 + 8;
  1828.                                 int xy = (mb_x * 2 + i + mb_y * 2 * mv_stride) << (mv_sample_log2 - 1);
  1829.                                 int mx = motion_val[direction][xy][0] >> shift;
  1830.                                 int my = motion_val[direction][xy][1] >> shift;
  1831.  
  1832.                                 if (IS_INTERLACED(mbtype_table[mb_index]))
  1833.                                     my *= 2;
  1834.  
  1835.                                 draw_arrow(ptr, sx, sy, mx + sx, my + sy, width,
  1836.                                            height, pict->linesize[0], 100, 0, direction);
  1837.                             }
  1838.                         } else {
  1839.                               int sx= mb_x * 16 + 8;
  1840.                               int sy= mb_y * 16 + 8;
  1841.                               int xy= (mb_x + mb_y * mv_stride) << mv_sample_log2;
  1842.                               int mx= (motion_val[direction][xy][0]>>shift) + sx;
  1843.                               int my= (motion_val[direction][xy][1]>>shift) + sy;
  1844.                               draw_arrow(ptr, sx, sy, mx, my, width, height, pict->linesize[0], 100, 0, direction);
  1845.                         }
  1846.                     }
  1847.                 }
  1848. #endif
  1849.                 if ((avctx->debug & FF_DEBUG_VIS_QP)) {
  1850.                     uint64_t c = (qscale_table[mb_index] * 128 / 31) *
  1851.                                  0x0101010101010101ULL;
  1852.                     int y;
  1853.                     for (y = 0; y < block_height; y++) {
  1854.                         *(uint64_t *)(pict->data[1] + 8 * mb_x +
  1855.                                       (block_height * mb_y + y) *
  1856.                                       pict->linesize[1]) = c;
  1857.                         *(uint64_t *)(pict->data[2] + 8 * mb_x +
  1858.                                       (block_height * mb_y + y) *
  1859.                                       pict->linesize[2]) = c;
  1860.                     }
  1861.                 }
  1862.                 if ((avctx->debug & FF_DEBUG_VIS_MB_TYPE) &&
  1863.                     motion_val[0]) {
  1864.                     int mb_type = mbtype_table[mb_index];
  1865.                     uint64_t u,v;
  1866.                     int y;
  1867. #define COLOR(theta, r) \
  1868.     u = (int)(128 + r * cos(theta * 3.141592 / 180)); \
  1869.     v = (int)(128 + r * sin(theta * 3.141592 / 180));
  1870.  
  1871.  
  1872.                     u = v = 128;
  1873.                     if (IS_PCM(mb_type)) {
  1874.                         COLOR(120, 48)
  1875.                     } else if ((IS_INTRA(mb_type) && IS_ACPRED(mb_type)) ||
  1876.                                IS_INTRA16x16(mb_type)) {
  1877.                         COLOR(30, 48)
  1878.                     } else if (IS_INTRA4x4(mb_type)) {
  1879.                         COLOR(90, 48)
  1880.                     } else if (IS_DIRECT(mb_type) && IS_SKIP(mb_type)) {
  1881.                         // COLOR(120, 48)
  1882.                     } else if (IS_DIRECT(mb_type)) {
  1883.                         COLOR(150, 48)
  1884.                     } else if (IS_GMC(mb_type) && IS_SKIP(mb_type)) {
  1885.                         COLOR(170, 48)
  1886.                     } else if (IS_GMC(mb_type)) {
  1887.                         COLOR(190, 48)
  1888.                     } else if (IS_SKIP(mb_type)) {
  1889.                         // COLOR(180, 48)
  1890.                     } else if (!USES_LIST(mb_type, 1)) {
  1891.                         COLOR(240, 48)
  1892.                     } else if (!USES_LIST(mb_type, 0)) {
  1893.                         COLOR(0, 48)
  1894.                     } else {
  1895.                         av_assert2(USES_LIST(mb_type, 0) && USES_LIST(mb_type, 1));
  1896.                         COLOR(300,48)
  1897.                     }
  1898.  
  1899.                     u *= 0x0101010101010101ULL;
  1900.                     v *= 0x0101010101010101ULL;
  1901.                     for (y = 0; y < block_height; y++) {
  1902.                         *(uint64_t *)(pict->data[1] + 8 * mb_x +
  1903.                                       (block_height * mb_y + y) * pict->linesize[1]) = u;
  1904.                         *(uint64_t *)(pict->data[2] + 8 * mb_x +
  1905.                                       (block_height * mb_y + y) * pict->linesize[2]) = v;
  1906.                     }
  1907.  
  1908.                     // segmentation
  1909.                     if (IS_8X8(mb_type) || IS_16X8(mb_type)) {
  1910.                         *(uint64_t *)(pict->data[0] + 16 * mb_x + 0 +
  1911.                                       (16 * mb_y + 8) * pict->linesize[0]) ^= 0x8080808080808080ULL;
  1912.                         *(uint64_t *)(pict->data[0] + 16 * mb_x + 8 +
  1913.                                       (16 * mb_y + 8) * pict->linesize[0]) ^= 0x8080808080808080ULL;
  1914.                     }
  1915.                     if (IS_8X8(mb_type) || IS_8X16(mb_type)) {
  1916.                         for (y = 0; y < 16; y++)
  1917.                             pict->data[0][16 * mb_x + 8 + (16 * mb_y + y) *
  1918.                                           pict->linesize[0]] ^= 0x80;
  1919.                     }
  1920.                     if (IS_8X8(mb_type) && mv_sample_log2 >= 2) {
  1921.                         int dm = 1 << (mv_sample_log2 - 2);
  1922.                         for (i = 0; i < 4; i++) {
  1923.                             int sx = mb_x * 16 + 8 * (i & 1);
  1924.                             int sy = mb_y * 16 + 8 * (i >> 1);
  1925.                             int xy = (mb_x * 2 + (i & 1) +
  1926.                                      (mb_y * 2 + (i >> 1)) * mv_stride) << (mv_sample_log2 - 1);
  1927.                             // FIXME bidir
  1928.                             int32_t *mv = (int32_t *) &motion_val[0][xy];
  1929.                             if (mv[0] != mv[dm] ||
  1930.                                 mv[dm * mv_stride] != mv[dm * (mv_stride + 1)])
  1931.                                 for (y = 0; y < 8; y++)
  1932.                                     pict->data[0][sx + 4 + (sy + y) * pict->linesize[0]] ^= 0x80;
  1933.                             if (mv[0] != mv[dm * mv_stride] || mv[dm] != mv[dm * (mv_stride + 1)])
  1934.                                 *(uint64_t *)(pict->data[0] + sx + (sy + 4) *
  1935.                                               pict->linesize[0]) ^= 0x8080808080808080ULL;
  1936.                         }
  1937.                     }
  1938.  
  1939.                     if (IS_INTERLACED(mb_type) &&
  1940.                         avctx->codec->id == AV_CODEC_ID_H264) {
  1941.                         // hmm
  1942.                     }
  1943.                 }
  1944.                 if (mbskip_table)
  1945.                     mbskip_table[mb_index] = 0;
  1946.             }
  1947.         }
  1948.     }
  1949. }
  1950.  
  1951. void ff_print_debug_info(MpegEncContext *s, Picture *p, AVFrame *pict)
  1952. {
  1953.     ff_print_debug_info2(s->avctx, pict, s->mbskip_table, p->mb_type,
  1954.                          p->qscale_table, p->motion_val, &s->low_delay,
  1955.                          s->mb_width, s->mb_height, s->mb_stride, s->quarter_sample);
  1956. }
  1957.  
  1958. int ff_mpv_export_qp_table(MpegEncContext *s, AVFrame *f, Picture *p, int qp_type)
  1959. {
  1960.     AVBufferRef *ref = av_buffer_ref(p->qscale_table_buf);
  1961.     int offset = 2*s->mb_stride + 1;
  1962.     if(!ref)
  1963.         return AVERROR(ENOMEM);
  1964.     av_assert0(ref->size >= offset + s->mb_stride * ((f->height+15)/16));
  1965.     ref->size -= offset;
  1966.     ref->data += offset;
  1967.     return av_frame_set_qp_table(f, ref, s->mb_stride, qp_type);
  1968. }
  1969.  
  1970. static inline int hpel_motion_lowres(MpegEncContext *s,
  1971.                                      uint8_t *dest, uint8_t *src,
  1972.                                      int field_based, int field_select,
  1973.                                      int src_x, int src_y,
  1974.                                      int width, int height, ptrdiff_t stride,
  1975.                                      int h_edge_pos, int v_edge_pos,
  1976.                                      int w, int h, h264_chroma_mc_func *pix_op,
  1977.                                      int motion_x, int motion_y)
  1978. {
  1979.     const int lowres   = s->avctx->lowres;
  1980.     const int op_index = FFMIN(lowres, 3);
  1981.     const int s_mask   = (2 << lowres) - 1;
  1982.     int emu = 0;
  1983.     int sx, sy;
  1984.  
  1985.     if (s->quarter_sample) {
  1986.         motion_x /= 2;
  1987.         motion_y /= 2;
  1988.     }
  1989.  
  1990.     sx = motion_x & s_mask;
  1991.     sy = motion_y & s_mask;
  1992.     src_x += motion_x >> lowres + 1;
  1993.     src_y += motion_y >> lowres + 1;
  1994.  
  1995.     src   += src_y * stride + src_x;
  1996.  
  1997.     if ((unsigned)src_x > FFMAX( h_edge_pos - (!!sx) - w,                 0) ||
  1998.         (unsigned)src_y > FFMAX((v_edge_pos >> field_based) - (!!sy) - h, 0)) {
  1999.         s->vdsp.emulated_edge_mc(s->sc.edge_emu_buffer, src,
  2000.                                  s->linesize, s->linesize,
  2001.                                  w + 1, (h + 1) << field_based,
  2002.                                  src_x, src_y   << field_based,
  2003.                                  h_edge_pos, v_edge_pos);
  2004.         src = s->sc.edge_emu_buffer;
  2005.         emu = 1;
  2006.     }
  2007.  
  2008.     sx = (sx << 2) >> lowres;
  2009.     sy = (sy << 2) >> lowres;
  2010.     if (field_select)
  2011.         src += s->linesize;
  2012.     pix_op[op_index](dest, src, stride, h, sx, sy);
  2013.     return emu;
  2014. }
  2015.  
  2016. /* apply one mpeg motion vector to the three components */
  2017. static av_always_inline void mpeg_motion_lowres(MpegEncContext *s,
  2018.                                                 uint8_t *dest_y,
  2019.                                                 uint8_t *dest_cb,
  2020.                                                 uint8_t *dest_cr,
  2021.                                                 int field_based,
  2022.                                                 int bottom_field,
  2023.                                                 int field_select,
  2024.                                                 uint8_t **ref_picture,
  2025.                                                 h264_chroma_mc_func *pix_op,
  2026.                                                 int motion_x, int motion_y,
  2027.                                                 int h, int mb_y)
  2028. {
  2029.     uint8_t *ptr_y, *ptr_cb, *ptr_cr;
  2030.     int mx, my, src_x, src_y, uvsrc_x, uvsrc_y, sx, sy, uvsx, uvsy;
  2031.     ptrdiff_t uvlinesize, linesize;
  2032.     const int lowres     = s->avctx->lowres;
  2033.     const int op_index   = FFMIN(lowres-1+s->chroma_x_shift, 3);
  2034.     const int block_s    = 8>>lowres;
  2035.     const int s_mask     = (2 << lowres) - 1;
  2036.     const int h_edge_pos = s->h_edge_pos >> lowres;
  2037.     const int v_edge_pos = s->v_edge_pos >> lowres;
  2038.     linesize   = s->current_picture.f->linesize[0] << field_based;
  2039.     uvlinesize = s->current_picture.f->linesize[1] << field_based;
  2040.  
  2041.     // FIXME obviously not perfect but qpel will not work in lowres anyway
  2042.     if (s->quarter_sample) {
  2043.         motion_x /= 2;
  2044.         motion_y /= 2;
  2045.     }
  2046.  
  2047.     if(field_based){
  2048.         motion_y += (bottom_field - field_select)*((1 << lowres)-1);
  2049.     }
  2050.  
  2051.     sx = motion_x & s_mask;
  2052.     sy = motion_y & s_mask;
  2053.     src_x = s->mb_x * 2 * block_s + (motion_x >> lowres + 1);
  2054.     src_y = (mb_y * 2 * block_s >> field_based) + (motion_y >> lowres + 1);
  2055.  
  2056.     if (s->out_format == FMT_H263) {
  2057.         uvsx    = ((motion_x >> 1) & s_mask) | (sx & 1);
  2058.         uvsy    = ((motion_y >> 1) & s_mask) | (sy & 1);
  2059.         uvsrc_x = src_x >> 1;
  2060.         uvsrc_y = src_y >> 1;
  2061.     } else if (s->out_format == FMT_H261) {
  2062.         // even chroma mv's are full pel in H261
  2063.         mx      = motion_x / 4;
  2064.         my      = motion_y / 4;
  2065.         uvsx    = (2 * mx) & s_mask;
  2066.         uvsy    = (2 * my) & s_mask;
  2067.         uvsrc_x = s->mb_x * block_s + (mx >> lowres);
  2068.         uvsrc_y =    mb_y * block_s + (my >> lowres);
  2069.     } else {
  2070.         if(s->chroma_y_shift){
  2071.             mx      = motion_x / 2;
  2072.             my      = motion_y / 2;
  2073.             uvsx    = mx & s_mask;
  2074.             uvsy    = my & s_mask;
  2075.             uvsrc_x = s->mb_x * block_s                 + (mx >> lowres + 1);
  2076.             uvsrc_y =   (mb_y * block_s >> field_based) + (my >> lowres + 1);
  2077.         } else {
  2078.             if(s->chroma_x_shift){
  2079.             //Chroma422
  2080.                 mx = motion_x / 2;
  2081.                 uvsx = mx & s_mask;
  2082.                 uvsy = motion_y & s_mask;
  2083.                 uvsrc_y = src_y;
  2084.                 uvsrc_x = s->mb_x*block_s               + (mx >> (lowres+1));
  2085.             } else {
  2086.             //Chroma444
  2087.                 uvsx = motion_x & s_mask;
  2088.                 uvsy = motion_y & s_mask;
  2089.                 uvsrc_x = src_x;
  2090.                 uvsrc_y = src_y;
  2091.             }
  2092.         }
  2093.     }
  2094.  
  2095.     ptr_y  = ref_picture[0] + src_y   * linesize   + src_x;
  2096.     ptr_cb = ref_picture[1] + uvsrc_y * uvlinesize + uvsrc_x;
  2097.     ptr_cr = ref_picture[2] + uvsrc_y * uvlinesize + uvsrc_x;
  2098.  
  2099.     if ((unsigned) src_x > FFMAX( h_edge_pos - (!!sx) - 2 * block_s,       0) || uvsrc_y<0 ||
  2100.         (unsigned) src_y > FFMAX((v_edge_pos >> field_based) - (!!sy) - h, 0)) {
  2101.         s->vdsp.emulated_edge_mc(s->sc.edge_emu_buffer, ptr_y,
  2102.                                  linesize >> field_based, linesize >> field_based,
  2103.                                  17, 17 + field_based,
  2104.                                 src_x, src_y << field_based, h_edge_pos,
  2105.                                 v_edge_pos);
  2106.         ptr_y = s->sc.edge_emu_buffer;
  2107.         if (!CONFIG_GRAY || !(s->avctx->flags & AV_CODEC_FLAG_GRAY)) {
  2108.             uint8_t *ubuf = s->sc.edge_emu_buffer + 18 * s->linesize;
  2109.             uint8_t *vbuf =ubuf + 9 * s->uvlinesize;
  2110.             s->vdsp.emulated_edge_mc(ubuf,  ptr_cb,
  2111.                                      uvlinesize >> field_based, uvlinesize >> field_based,
  2112.                                      9, 9 + field_based,
  2113.                                     uvsrc_x, uvsrc_y << field_based,
  2114.                                     h_edge_pos >> 1, v_edge_pos >> 1);
  2115.             s->vdsp.emulated_edge_mc(vbuf,  ptr_cr,
  2116.                                      uvlinesize >> field_based,uvlinesize >> field_based,
  2117.                                      9, 9 + field_based,
  2118.                                     uvsrc_x, uvsrc_y << field_based,
  2119.                                     h_edge_pos >> 1, v_edge_pos >> 1);
  2120.             ptr_cb = ubuf;
  2121.             ptr_cr = vbuf;
  2122.         }
  2123.     }
  2124.  
  2125.     // FIXME use this for field pix too instead of the obnoxious hack which changes picture.f->data
  2126.     if (bottom_field) {
  2127.         dest_y  += s->linesize;
  2128.         dest_cb += s->uvlinesize;
  2129.         dest_cr += s->uvlinesize;
  2130.     }
  2131.  
  2132.     if (field_select) {
  2133.         ptr_y   += s->linesize;
  2134.         ptr_cb  += s->uvlinesize;
  2135.         ptr_cr  += s->uvlinesize;
  2136.     }
  2137.  
  2138.     sx = (sx << 2) >> lowres;
  2139.     sy = (sy << 2) >> lowres;
  2140.     pix_op[lowres - 1](dest_y, ptr_y, linesize, h, sx, sy);
  2141.  
  2142.     if (!CONFIG_GRAY || !(s->avctx->flags & AV_CODEC_FLAG_GRAY)) {
  2143.         int hc = s->chroma_y_shift ? (h+1-bottom_field)>>1 : h;
  2144.         uvsx = (uvsx << 2) >> lowres;
  2145.         uvsy = (uvsy << 2) >> lowres;
  2146.         if (hc) {
  2147.             pix_op[op_index](dest_cb, ptr_cb, uvlinesize, hc, uvsx, uvsy);
  2148.             pix_op[op_index](dest_cr, ptr_cr, uvlinesize, hc, uvsx, uvsy);
  2149.         }
  2150.     }
  2151.     // FIXME h261 lowres loop filter
  2152. }
  2153.  
  2154. static inline void chroma_4mv_motion_lowres(MpegEncContext *s,
  2155.                                             uint8_t *dest_cb, uint8_t *dest_cr,
  2156.                                             uint8_t **ref_picture,
  2157.                                             h264_chroma_mc_func * pix_op,
  2158.                                             int mx, int my)
  2159. {
  2160.     const int lowres     = s->avctx->lowres;
  2161.     const int op_index   = FFMIN(lowres, 3);
  2162.     const int block_s    = 8 >> lowres;
  2163.     const int s_mask     = (2 << lowres) - 1;
  2164.     const int h_edge_pos = s->h_edge_pos >> lowres + 1;
  2165.     const int v_edge_pos = s->v_edge_pos >> lowres + 1;
  2166.     int emu = 0, src_x, src_y, sx, sy;
  2167.     ptrdiff_t offset;
  2168.     uint8_t *ptr;
  2169.  
  2170.     if (s->quarter_sample) {
  2171.         mx /= 2;
  2172.         my /= 2;
  2173.     }
  2174.  
  2175.     /* In case of 8X8, we construct a single chroma motion vector
  2176.        with a special rounding */
  2177.     mx = ff_h263_round_chroma(mx);
  2178.     my = ff_h263_round_chroma(my);
  2179.  
  2180.     sx = mx & s_mask;
  2181.     sy = my & s_mask;
  2182.     src_x = s->mb_x * block_s + (mx >> lowres + 1);
  2183.     src_y = s->mb_y * block_s + (my >> lowres + 1);
  2184.  
  2185.     offset = src_y * s->uvlinesize + src_x;
  2186.     ptr = ref_picture[1] + offset;
  2187.     if ((unsigned) src_x > FFMAX(h_edge_pos - (!!sx) - block_s, 0) ||
  2188.         (unsigned) src_y > FFMAX(v_edge_pos - (!!sy) - block_s, 0)) {
  2189.         s->vdsp.emulated_edge_mc(s->sc.edge_emu_buffer, ptr,
  2190.                                  s->uvlinesize, s->uvlinesize,
  2191.                                  9, 9,
  2192.                                  src_x, src_y, h_edge_pos, v_edge_pos);
  2193.         ptr = s->sc.edge_emu_buffer;
  2194.         emu = 1;
  2195.     }
  2196.     sx = (sx << 2) >> lowres;
  2197.     sy = (sy << 2) >> lowres;
  2198.     pix_op[op_index](dest_cb, ptr, s->uvlinesize, block_s, sx, sy);
  2199.  
  2200.     ptr = ref_picture[2] + offset;
  2201.     if (emu) {
  2202.         s->vdsp.emulated_edge_mc(s->sc.edge_emu_buffer, ptr,
  2203.                                  s->uvlinesize, s->uvlinesize,
  2204.                                  9, 9,
  2205.                                  src_x, src_y, h_edge_pos, v_edge_pos);
  2206.         ptr = s->sc.edge_emu_buffer;
  2207.     }
  2208.     pix_op[op_index](dest_cr, ptr, s->uvlinesize, block_s, sx, sy);
  2209. }
  2210.  
  2211. /**
  2212.  * motion compensation of a single macroblock
  2213.  * @param s context
  2214.  * @param dest_y luma destination pointer
  2215.  * @param dest_cb chroma cb/u destination pointer
  2216.  * @param dest_cr chroma cr/v destination pointer
  2217.  * @param dir direction (0->forward, 1->backward)
  2218.  * @param ref_picture array[3] of pointers to the 3 planes of the reference picture
  2219.  * @param pix_op halfpel motion compensation function (average or put normally)
  2220.  * the motion vectors are taken from s->mv and the MV type from s->mv_type
  2221.  */
  2222. static inline void MPV_motion_lowres(MpegEncContext *s,
  2223.                                      uint8_t *dest_y, uint8_t *dest_cb,
  2224.                                      uint8_t *dest_cr,
  2225.                                      int dir, uint8_t **ref_picture,
  2226.                                      h264_chroma_mc_func *pix_op)
  2227. {
  2228.     int mx, my;
  2229.     int mb_x, mb_y, i;
  2230.     const int lowres  = s->avctx->lowres;
  2231.     const int block_s = 8 >>lowres;
  2232.  
  2233.     mb_x = s->mb_x;
  2234.     mb_y = s->mb_y;
  2235.  
  2236.     switch (s->mv_type) {
  2237.     case MV_TYPE_16X16:
  2238.         mpeg_motion_lowres(s, dest_y, dest_cb, dest_cr,
  2239.                            0, 0, 0,
  2240.                            ref_picture, pix_op,
  2241.                            s->mv[dir][0][0], s->mv[dir][0][1],
  2242.                            2 * block_s, mb_y);
  2243.         break;
  2244.     case MV_TYPE_8X8:
  2245.         mx = 0;
  2246.         my = 0;
  2247.         for (i = 0; i < 4; i++) {
  2248.             hpel_motion_lowres(s, dest_y + ((i & 1) + (i >> 1) *
  2249.                                s->linesize) * block_s,
  2250.                                ref_picture[0], 0, 0,
  2251.                                (2 * mb_x + (i & 1)) * block_s,
  2252.                                (2 * mb_y + (i >> 1)) * block_s,
  2253.                                s->width, s->height, s->linesize,
  2254.                                s->h_edge_pos >> lowres, s->v_edge_pos >> lowres,
  2255.                                block_s, block_s, pix_op,
  2256.                                s->mv[dir][i][0], s->mv[dir][i][1]);
  2257.  
  2258.             mx += s->mv[dir][i][0];
  2259.             my += s->mv[dir][i][1];
  2260.         }
  2261.  
  2262.         if (!CONFIG_GRAY || !(s->avctx->flags & AV_CODEC_FLAG_GRAY))
  2263.             chroma_4mv_motion_lowres(s, dest_cb, dest_cr, ref_picture,
  2264.                                      pix_op, mx, my);
  2265.         break;
  2266.     case MV_TYPE_FIELD:
  2267.         if (s->picture_structure == PICT_FRAME) {
  2268.             /* top field */
  2269.             mpeg_motion_lowres(s, dest_y, dest_cb, dest_cr,
  2270.                                1, 0, s->field_select[dir][0],
  2271.                                ref_picture, pix_op,
  2272.                                s->mv[dir][0][0], s->mv[dir][0][1],
  2273.                                block_s, mb_y);
  2274.             /* bottom field */
  2275.             mpeg_motion_lowres(s, dest_y, dest_cb, dest_cr,
  2276.                                1, 1, s->field_select[dir][1],
  2277.                                ref_picture, pix_op,
  2278.                                s->mv[dir][1][0], s->mv[dir][1][1],
  2279.                                block_s, mb_y);
  2280.         } else {
  2281.             if (s->picture_structure != s->field_select[dir][0] + 1 &&
  2282.                 s->pict_type != AV_PICTURE_TYPE_B && !s->first_field) {
  2283.                 ref_picture = s->current_picture_ptr->f->data;
  2284.  
  2285.             }
  2286.             mpeg_motion_lowres(s, dest_y, dest_cb, dest_cr,
  2287.                                0, 0, s->field_select[dir][0],
  2288.                                ref_picture, pix_op,
  2289.                                s->mv[dir][0][0],
  2290.                                s->mv[dir][0][1], 2 * block_s, mb_y >> 1);
  2291.             }
  2292.         break;
  2293.     case MV_TYPE_16X8:
  2294.         for (i = 0; i < 2; i++) {
  2295.             uint8_t **ref2picture;
  2296.  
  2297.             if (s->picture_structure == s->field_select[dir][i] + 1 ||
  2298.                 s->pict_type == AV_PICTURE_TYPE_B || s->first_field) {
  2299.                 ref2picture = ref_picture;
  2300.             } else {
  2301.                 ref2picture = s->current_picture_ptr->f->data;
  2302.             }
  2303.  
  2304.             mpeg_motion_lowres(s, dest_y, dest_cb, dest_cr,
  2305.                                0, 0, s->field_select[dir][i],
  2306.                                ref2picture, pix_op,
  2307.                                s->mv[dir][i][0], s->mv[dir][i][1] +
  2308.                                2 * block_s * i, block_s, mb_y >> 1);
  2309.  
  2310.             dest_y  +=  2 * block_s *  s->linesize;
  2311.             dest_cb += (2 * block_s >> s->chroma_y_shift) * s->uvlinesize;
  2312.             dest_cr += (2 * block_s >> s->chroma_y_shift) * s->uvlinesize;
  2313.         }
  2314.         break;
  2315.     case MV_TYPE_DMV:
  2316.         if (s->picture_structure == PICT_FRAME) {
  2317.             for (i = 0; i < 2; i++) {
  2318.                 int j;
  2319.                 for (j = 0; j < 2; j++) {
  2320.                     mpeg_motion_lowres(s, dest_y, dest_cb, dest_cr,
  2321.                                        1, j, j ^ i,
  2322.                                        ref_picture, pix_op,
  2323.                                        s->mv[dir][2 * i + j][0],
  2324.                                        s->mv[dir][2 * i + j][1],
  2325.                                        block_s, mb_y);
  2326.                 }
  2327.                 pix_op = s->h264chroma.avg_h264_chroma_pixels_tab;
  2328.             }
  2329.         } else {
  2330.             for (i = 0; i < 2; i++) {
  2331.                 mpeg_motion_lowres(s, dest_y, dest_cb, dest_cr,
  2332.                                    0, 0, s->picture_structure != i + 1,
  2333.                                    ref_picture, pix_op,
  2334.                                    s->mv[dir][2 * i][0],s->mv[dir][2 * i][1],
  2335.                                    2 * block_s, mb_y >> 1);
  2336.  
  2337.                 // after put we make avg of the same block
  2338.                 pix_op = s->h264chroma.avg_h264_chroma_pixels_tab;
  2339.  
  2340.                 // opposite parity is always in the same
  2341.                 // frame if this is second field
  2342.                 if (!s->first_field) {
  2343.                     ref_picture = s->current_picture_ptr->f->data;
  2344.                 }
  2345.             }
  2346.         }
  2347.         break;
  2348.     default:
  2349.         av_assert2(0);
  2350.     }
  2351. }
  2352.  
  2353. /**
  2354.  * find the lowest MB row referenced in the MVs
  2355.  */
  2356. static int lowest_referenced_row(MpegEncContext *s, int dir)
  2357. {
  2358.     int my_max = INT_MIN, my_min = INT_MAX, qpel_shift = !s->quarter_sample;
  2359.     int my, off, i, mvs;
  2360.  
  2361.     if (s->picture_structure != PICT_FRAME || s->mcsel)
  2362.         goto unhandled;
  2363.  
  2364.     switch (s->mv_type) {
  2365.         case MV_TYPE_16X16:
  2366.             mvs = 1;
  2367.             break;
  2368.         case MV_TYPE_16X8:
  2369.             mvs = 2;
  2370.             break;
  2371.         case MV_TYPE_8X8:
  2372.             mvs = 4;
  2373.             break;
  2374.         default:
  2375.             goto unhandled;
  2376.     }
  2377.  
  2378.     for (i = 0; i < mvs; i++) {
  2379.         my = s->mv[dir][i][1];
  2380.         my_max = FFMAX(my_max, my);
  2381.         my_min = FFMIN(my_min, my);
  2382.     }
  2383.  
  2384.     off = ((FFMAX(-my_min, my_max)<<qpel_shift) + 63) >> 6;
  2385.  
  2386.     return av_clip(s->mb_y + off, 0, s->mb_height - 1);
  2387. unhandled:
  2388.     return s->mb_height-1;
  2389. }
  2390.  
  2391. /* put block[] to dest[] */
  2392. static inline void put_dct(MpegEncContext *s,
  2393.                            int16_t *block, int i, uint8_t *dest, int line_size, int qscale)
  2394. {
  2395.     s->dct_unquantize_intra(s, block, i, qscale);
  2396.     s->idsp.idct_put(dest, line_size, block);
  2397. }
  2398.  
  2399. /* add block[] to dest[] */
  2400. static inline void add_dct(MpegEncContext *s,
  2401.                            int16_t *block, int i, uint8_t *dest, int line_size)
  2402. {
  2403.     if (s->block_last_index[i] >= 0) {
  2404.         s->idsp.idct_add(dest, line_size, block);
  2405.     }
  2406. }
  2407.  
  2408. static inline void add_dequant_dct(MpegEncContext *s,
  2409.                            int16_t *block, int i, uint8_t *dest, int line_size, int qscale)
  2410. {
  2411.     if (s->block_last_index[i] >= 0) {
  2412.         s->dct_unquantize_inter(s, block, i, qscale);
  2413.  
  2414.         s->idsp.idct_add(dest, line_size, block);
  2415.     }
  2416. }
  2417.  
  2418. /**
  2419.  * Clean dc, ac, coded_block for the current non-intra MB.
  2420.  */
  2421. void ff_clean_intra_table_entries(MpegEncContext *s)
  2422. {
  2423.     int wrap = s->b8_stride;
  2424.     int xy = s->block_index[0];
  2425.  
  2426.     s->dc_val[0][xy           ] =
  2427.     s->dc_val[0][xy + 1       ] =
  2428.     s->dc_val[0][xy     + wrap] =
  2429.     s->dc_val[0][xy + 1 + wrap] = 1024;
  2430.     /* ac pred */
  2431.     memset(s->ac_val[0][xy       ], 0, 32 * sizeof(int16_t));
  2432.     memset(s->ac_val[0][xy + wrap], 0, 32 * sizeof(int16_t));
  2433.     if (s->msmpeg4_version>=3) {
  2434.         s->coded_block[xy           ] =
  2435.         s->coded_block[xy + 1       ] =
  2436.         s->coded_block[xy     + wrap] =
  2437.         s->coded_block[xy + 1 + wrap] = 0;
  2438.     }
  2439.     /* chroma */
  2440.     wrap = s->mb_stride;
  2441.     xy = s->mb_x + s->mb_y * wrap;
  2442.     s->dc_val[1][xy] =
  2443.     s->dc_val[2][xy] = 1024;
  2444.     /* ac pred */
  2445.     memset(s->ac_val[1][xy], 0, 16 * sizeof(int16_t));
  2446.     memset(s->ac_val[2][xy], 0, 16 * sizeof(int16_t));
  2447.  
  2448.     s->mbintra_table[xy]= 0;
  2449. }
  2450.  
  2451. /* generic function called after a macroblock has been parsed by the
  2452.    decoder or after it has been encoded by the encoder.
  2453.  
  2454.    Important variables used:
  2455.    s->mb_intra : true if intra macroblock
  2456.    s->mv_dir   : motion vector direction
  2457.    s->mv_type  : motion vector type
  2458.    s->mv       : motion vector
  2459.    s->interlaced_dct : true if interlaced dct used (mpeg2)
  2460.  */
  2461. static av_always_inline
  2462. void mpv_decode_mb_internal(MpegEncContext *s, int16_t block[12][64],
  2463.                             int lowres_flag, int is_mpeg12)
  2464. {
  2465.     const int mb_xy = s->mb_y * s->mb_stride + s->mb_x;
  2466.  
  2467.     if (CONFIG_XVMC &&
  2468.         s->avctx->hwaccel && s->avctx->hwaccel->decode_mb) {
  2469.         s->avctx->hwaccel->decode_mb(s);//xvmc uses pblocks
  2470.         return;
  2471.     }
  2472.  
  2473.     if(s->avctx->debug&FF_DEBUG_DCT_COEFF) {
  2474.        /* print DCT coefficients */
  2475.        int i,j;
  2476.        av_log(s->avctx, AV_LOG_DEBUG, "DCT coeffs of MB at %dx%d:\n", s->mb_x, s->mb_y);
  2477.        for(i=0; i<6; i++){
  2478.            for(j=0; j<64; j++){
  2479.                av_log(s->avctx, AV_LOG_DEBUG, "%5d",
  2480.                       block[i][s->idsp.idct_permutation[j]]);
  2481.            }
  2482.            av_log(s->avctx, AV_LOG_DEBUG, "\n");
  2483.        }
  2484.     }
  2485.  
  2486.     s->current_picture.qscale_table[mb_xy] = s->qscale;
  2487.  
  2488.     /* update DC predictors for P macroblocks */
  2489.     if (!s->mb_intra) {
  2490.         if (!is_mpeg12 && (s->h263_pred || s->h263_aic)) {
  2491.             if(s->mbintra_table[mb_xy])
  2492.                 ff_clean_intra_table_entries(s);
  2493.         } else {
  2494.             s->last_dc[0] =
  2495.             s->last_dc[1] =
  2496.             s->last_dc[2] = 128 << s->intra_dc_precision;
  2497.         }
  2498.     }
  2499.     else if (!is_mpeg12 && (s->h263_pred || s->h263_aic))
  2500.         s->mbintra_table[mb_xy]=1;
  2501.  
  2502.     if ((s->avctx->flags & AV_CODEC_FLAG_PSNR) || s->avctx->frame_skip_threshold || s->avctx->frame_skip_factor ||
  2503.         !(s->encoding && (s->intra_only || s->pict_type == AV_PICTURE_TYPE_B) &&
  2504.           s->avctx->mb_decision != FF_MB_DECISION_RD)) { // FIXME precalc
  2505.         uint8_t *dest_y, *dest_cb, *dest_cr;
  2506.         int dct_linesize, dct_offset;
  2507.         op_pixels_func (*op_pix)[4];
  2508.         qpel_mc_func (*op_qpix)[16];
  2509.         const int linesize   = s->current_picture.f->linesize[0]; //not s->linesize as this would be wrong for field pics
  2510.         const int uvlinesize = s->current_picture.f->linesize[1];
  2511.         const int readable= s->pict_type != AV_PICTURE_TYPE_B || s->encoding || s->avctx->draw_horiz_band || lowres_flag;
  2512.         const int block_size= lowres_flag ? 8>>s->avctx->lowres : 8;
  2513.  
  2514.         /* avoid copy if macroblock skipped in last frame too */
  2515.         /* skip only during decoding as we might trash the buffers during encoding a bit */
  2516.         if(!s->encoding){
  2517.             uint8_t *mbskip_ptr = &s->mbskip_table[mb_xy];
  2518.  
  2519.             if (s->mb_skipped) {
  2520.                 s->mb_skipped= 0;
  2521.                 av_assert2(s->pict_type!=AV_PICTURE_TYPE_I);
  2522.                 *mbskip_ptr = 1;
  2523.             } else if(!s->current_picture.reference) {
  2524.                 *mbskip_ptr = 1;
  2525.             } else{
  2526.                 *mbskip_ptr = 0; /* not skipped */
  2527.             }
  2528.         }
  2529.  
  2530.         dct_linesize = linesize << s->interlaced_dct;
  2531.         dct_offset   = s->interlaced_dct ? linesize : linesize * block_size;
  2532.  
  2533.         if(readable){
  2534.             dest_y=  s->dest[0];
  2535.             dest_cb= s->dest[1];
  2536.             dest_cr= s->dest[2];
  2537.         }else{
  2538.             dest_y = s->sc.b_scratchpad;
  2539.             dest_cb= s->sc.b_scratchpad+16*linesize;
  2540.             dest_cr= s->sc.b_scratchpad+32*linesize;
  2541.         }
  2542.  
  2543.         if (!s->mb_intra) {
  2544.             /* motion handling */
  2545.             /* decoding or more than one mb_type (MC was already done otherwise) */
  2546.             if(!s->encoding){
  2547.  
  2548.                 if(HAVE_THREADS && s->avctx->active_thread_type&FF_THREAD_FRAME) {
  2549.                     if (s->mv_dir & MV_DIR_FORWARD) {
  2550.                         ff_thread_await_progress(&s->last_picture_ptr->tf,
  2551.                                                  lowest_referenced_row(s, 0),
  2552.                                                  0);
  2553.                     }
  2554.                     if (s->mv_dir & MV_DIR_BACKWARD) {
  2555.                         ff_thread_await_progress(&s->next_picture_ptr->tf,
  2556.                                                  lowest_referenced_row(s, 1),
  2557.                                                  0);
  2558.                     }
  2559.                 }
  2560.  
  2561.                 if(lowres_flag){
  2562.                     h264_chroma_mc_func *op_pix = s->h264chroma.put_h264_chroma_pixels_tab;
  2563.  
  2564.                     if (s->mv_dir & MV_DIR_FORWARD) {
  2565.                         MPV_motion_lowres(s, dest_y, dest_cb, dest_cr, 0, s->last_picture.f->data, op_pix);
  2566.                         op_pix = s->h264chroma.avg_h264_chroma_pixels_tab;
  2567.                     }
  2568.                     if (s->mv_dir & MV_DIR_BACKWARD) {
  2569.                         MPV_motion_lowres(s, dest_y, dest_cb, dest_cr, 1, s->next_picture.f->data, op_pix);
  2570.                     }
  2571.                 }else{
  2572.                     op_qpix = s->me.qpel_put;
  2573.                     if ((!s->no_rounding) || s->pict_type==AV_PICTURE_TYPE_B){
  2574.                         op_pix = s->hdsp.put_pixels_tab;
  2575.                     }else{
  2576.                         op_pix = s->hdsp.put_no_rnd_pixels_tab;
  2577.                     }
  2578.                     if (s->mv_dir & MV_DIR_FORWARD) {
  2579.                         ff_mpv_motion(s, dest_y, dest_cb, dest_cr, 0, s->last_picture.f->data, op_pix, op_qpix);
  2580.                         op_pix = s->hdsp.avg_pixels_tab;
  2581.                         op_qpix= s->me.qpel_avg;
  2582.                     }
  2583.                     if (s->mv_dir & MV_DIR_BACKWARD) {
  2584.                         ff_mpv_motion(s, dest_y, dest_cb, dest_cr, 1, s->next_picture.f->data, op_pix, op_qpix);
  2585.                     }
  2586.                 }
  2587.             }
  2588.  
  2589.             /* skip dequant / idct if we are really late ;) */
  2590.             if(s->avctx->skip_idct){
  2591.                 if(  (s->avctx->skip_idct >= AVDISCARD_NONREF && s->pict_type == AV_PICTURE_TYPE_B)
  2592.                    ||(s->avctx->skip_idct >= AVDISCARD_NONKEY && s->pict_type != AV_PICTURE_TYPE_I)
  2593.                    || s->avctx->skip_idct >= AVDISCARD_ALL)
  2594.                     goto skip_idct;
  2595.             }
  2596.  
  2597.             /* add dct residue */
  2598.             if(s->encoding || !(   s->msmpeg4_version || s->codec_id==AV_CODEC_ID_MPEG1VIDEO || s->codec_id==AV_CODEC_ID_MPEG2VIDEO
  2599.                                 || (s->codec_id==AV_CODEC_ID_MPEG4 && !s->mpeg_quant))){
  2600.                 add_dequant_dct(s, block[0], 0, dest_y                          , dct_linesize, s->qscale);
  2601.                 add_dequant_dct(s, block[1], 1, dest_y              + block_size, dct_linesize, s->qscale);
  2602.                 add_dequant_dct(s, block[2], 2, dest_y + dct_offset             , dct_linesize, s->qscale);
  2603.                 add_dequant_dct(s, block[3], 3, dest_y + dct_offset + block_size, dct_linesize, s->qscale);
  2604.  
  2605.                 if (!CONFIG_GRAY || !(s->avctx->flags & AV_CODEC_FLAG_GRAY)) {
  2606.                     if (s->chroma_y_shift){
  2607.                         add_dequant_dct(s, block[4], 4, dest_cb, uvlinesize, s->chroma_qscale);
  2608.                         add_dequant_dct(s, block[5], 5, dest_cr, uvlinesize, s->chroma_qscale);
  2609.                     }else{
  2610.                         dct_linesize >>= 1;
  2611.                         dct_offset >>=1;
  2612.                         add_dequant_dct(s, block[4], 4, dest_cb,              dct_linesize, s->chroma_qscale);
  2613.                         add_dequant_dct(s, block[5], 5, dest_cr,              dct_linesize, s->chroma_qscale);
  2614.                         add_dequant_dct(s, block[6], 6, dest_cb + dct_offset, dct_linesize, s->chroma_qscale);
  2615.                         add_dequant_dct(s, block[7], 7, dest_cr + dct_offset, dct_linesize, s->chroma_qscale);
  2616.                     }
  2617.                 }
  2618.             } else if(is_mpeg12 || (s->codec_id != AV_CODEC_ID_WMV2)){
  2619.                 add_dct(s, block[0], 0, dest_y                          , dct_linesize);
  2620.                 add_dct(s, block[1], 1, dest_y              + block_size, dct_linesize);
  2621.                 add_dct(s, block[2], 2, dest_y + dct_offset             , dct_linesize);
  2622.                 add_dct(s, block[3], 3, dest_y + dct_offset + block_size, dct_linesize);
  2623.  
  2624.                 if (!CONFIG_GRAY || !(s->avctx->flags & AV_CODEC_FLAG_GRAY)) {
  2625.                     if(s->chroma_y_shift){//Chroma420
  2626.                         add_dct(s, block[4], 4, dest_cb, uvlinesize);
  2627.                         add_dct(s, block[5], 5, dest_cr, uvlinesize);
  2628.                     }else{
  2629.                         //chroma422
  2630.                         dct_linesize = uvlinesize << s->interlaced_dct;
  2631.                         dct_offset   = s->interlaced_dct ? uvlinesize : uvlinesize*block_size;
  2632.  
  2633.                         add_dct(s, block[4], 4, dest_cb, dct_linesize);
  2634.                         add_dct(s, block[5], 5, dest_cr, dct_linesize);
  2635.                         add_dct(s, block[6], 6, dest_cb+dct_offset, dct_linesize);
  2636.                         add_dct(s, block[7], 7, dest_cr+dct_offset, dct_linesize);
  2637.                         if(!s->chroma_x_shift){//Chroma444
  2638.                             add_dct(s, block[8], 8, dest_cb+block_size, dct_linesize);
  2639.                             add_dct(s, block[9], 9, dest_cr+block_size, dct_linesize);
  2640.                             add_dct(s, block[10], 10, dest_cb+block_size+dct_offset, dct_linesize);
  2641.                             add_dct(s, block[11], 11, dest_cr+block_size+dct_offset, dct_linesize);
  2642.                         }
  2643.                     }
  2644.                 }//fi gray
  2645.             }
  2646.             else if (CONFIG_WMV2_DECODER || CONFIG_WMV2_ENCODER) {
  2647.                 ff_wmv2_add_mb(s, block, dest_y, dest_cb, dest_cr);
  2648.             }
  2649.         } else {
  2650.             /* dct only in intra block */
  2651.             if(s->encoding || !(s->codec_id==AV_CODEC_ID_MPEG1VIDEO || s->codec_id==AV_CODEC_ID_MPEG2VIDEO)){
  2652.                 put_dct(s, block[0], 0, dest_y                          , dct_linesize, s->qscale);
  2653.                 put_dct(s, block[1], 1, dest_y              + block_size, dct_linesize, s->qscale);
  2654.                 put_dct(s, block[2], 2, dest_y + dct_offset             , dct_linesize, s->qscale);
  2655.                 put_dct(s, block[3], 3, dest_y + dct_offset + block_size, dct_linesize, s->qscale);
  2656.  
  2657.                 if (!CONFIG_GRAY || !(s->avctx->flags & AV_CODEC_FLAG_GRAY)) {
  2658.                     if(s->chroma_y_shift){
  2659.                         put_dct(s, block[4], 4, dest_cb, uvlinesize, s->chroma_qscale);
  2660.                         put_dct(s, block[5], 5, dest_cr, uvlinesize, s->chroma_qscale);
  2661.                     }else{
  2662.                         dct_offset >>=1;
  2663.                         dct_linesize >>=1;
  2664.                         put_dct(s, block[4], 4, dest_cb,              dct_linesize, s->chroma_qscale);
  2665.                         put_dct(s, block[5], 5, dest_cr,              dct_linesize, s->chroma_qscale);
  2666.                         put_dct(s, block[6], 6, dest_cb + dct_offset, dct_linesize, s->chroma_qscale);
  2667.                         put_dct(s, block[7], 7, dest_cr + dct_offset, dct_linesize, s->chroma_qscale);
  2668.                     }
  2669.                 }
  2670.             }else{
  2671.                 s->idsp.idct_put(dest_y,                           dct_linesize, block[0]);
  2672.                 s->idsp.idct_put(dest_y              + block_size, dct_linesize, block[1]);
  2673.                 s->idsp.idct_put(dest_y + dct_offset,              dct_linesize, block[2]);
  2674.                 s->idsp.idct_put(dest_y + dct_offset + block_size, dct_linesize, block[3]);
  2675.  
  2676.                 if (!CONFIG_GRAY || !(s->avctx->flags & AV_CODEC_FLAG_GRAY)) {
  2677.                     if(s->chroma_y_shift){
  2678.                         s->idsp.idct_put(dest_cb, uvlinesize, block[4]);
  2679.                         s->idsp.idct_put(dest_cr, uvlinesize, block[5]);
  2680.                     }else{
  2681.  
  2682.                         dct_linesize = uvlinesize << s->interlaced_dct;
  2683.                         dct_offset   = s->interlaced_dct ? uvlinesize : uvlinesize*block_size;
  2684.  
  2685.                         s->idsp.idct_put(dest_cb,              dct_linesize, block[4]);
  2686.                         s->idsp.idct_put(dest_cr,              dct_linesize, block[5]);
  2687.                         s->idsp.idct_put(dest_cb + dct_offset, dct_linesize, block[6]);
  2688.                         s->idsp.idct_put(dest_cr + dct_offset, dct_linesize, block[7]);
  2689.                         if(!s->chroma_x_shift){//Chroma444
  2690.                             s->idsp.idct_put(dest_cb + block_size,              dct_linesize, block[8]);
  2691.                             s->idsp.idct_put(dest_cr + block_size,              dct_linesize, block[9]);
  2692.                             s->idsp.idct_put(dest_cb + block_size + dct_offset, dct_linesize, block[10]);
  2693.                             s->idsp.idct_put(dest_cr + block_size + dct_offset, dct_linesize, block[11]);
  2694.                         }
  2695.                     }
  2696.                 }//gray
  2697.             }
  2698.         }
  2699. skip_idct:
  2700.         if(!readable){
  2701.             s->hdsp.put_pixels_tab[0][0](s->dest[0], dest_y ,   linesize,16);
  2702.             if (!CONFIG_GRAY || !(s->avctx->flags & AV_CODEC_FLAG_GRAY)) {
  2703.                 s->hdsp.put_pixels_tab[s->chroma_x_shift][0](s->dest[1], dest_cb, uvlinesize,16 >> s->chroma_y_shift);
  2704.                 s->hdsp.put_pixels_tab[s->chroma_x_shift][0](s->dest[2], dest_cr, uvlinesize,16 >> s->chroma_y_shift);
  2705.             }
  2706.         }
  2707.     }
  2708. }
  2709.  
  2710. void ff_mpv_decode_mb(MpegEncContext *s, int16_t block[12][64])
  2711. {
  2712. #if !CONFIG_SMALL
  2713.     if(s->out_format == FMT_MPEG1) {
  2714.         if(s->avctx->lowres) mpv_decode_mb_internal(s, block, 1, 1);
  2715.         else                 mpv_decode_mb_internal(s, block, 0, 1);
  2716.     } else
  2717. #endif
  2718.     if(s->avctx->lowres) mpv_decode_mb_internal(s, block, 1, 0);
  2719.     else                  mpv_decode_mb_internal(s, block, 0, 0);
  2720. }
  2721.  
  2722. void ff_mpeg_draw_horiz_band(MpegEncContext *s, int y, int h)
  2723. {
  2724.     ff_draw_horiz_band(s->avctx, s->current_picture_ptr->f,
  2725.                        s->last_picture_ptr ? s->last_picture_ptr->f : NULL, y, h, s->picture_structure,
  2726.                        s->first_field, s->low_delay);
  2727. }
  2728.  
  2729. void ff_init_block_index(MpegEncContext *s){ //FIXME maybe rename
  2730.     const int linesize   = s->current_picture.f->linesize[0]; //not s->linesize as this would be wrong for field pics
  2731.     const int uvlinesize = s->current_picture.f->linesize[1];
  2732.     const int mb_size= 4 - s->avctx->lowres;
  2733.  
  2734.     s->block_index[0]= s->b8_stride*(s->mb_y*2    ) - 2 + s->mb_x*2;
  2735.     s->block_index[1]= s->b8_stride*(s->mb_y*2    ) - 1 + s->mb_x*2;
  2736.     s->block_index[2]= s->b8_stride*(s->mb_y*2 + 1) - 2 + s->mb_x*2;
  2737.     s->block_index[3]= s->b8_stride*(s->mb_y*2 + 1) - 1 + s->mb_x*2;
  2738.     s->block_index[4]= s->mb_stride*(s->mb_y + 1)                + s->b8_stride*s->mb_height*2 + s->mb_x - 1;
  2739.     s->block_index[5]= s->mb_stride*(s->mb_y + s->mb_height + 2) + s->b8_stride*s->mb_height*2 + s->mb_x - 1;
  2740.     //block_index is not used by mpeg2, so it is not affected by chroma_format
  2741.  
  2742.     s->dest[0] = s->current_picture.f->data[0] + (int)((s->mb_x - 1U) <<  mb_size);
  2743.     s->dest[1] = s->current_picture.f->data[1] + (int)((s->mb_x - 1U) << (mb_size - s->chroma_x_shift));
  2744.     s->dest[2] = s->current_picture.f->data[2] + (int)((s->mb_x - 1U) << (mb_size - s->chroma_x_shift));
  2745.  
  2746.     if(!(s->pict_type==AV_PICTURE_TYPE_B && s->avctx->draw_horiz_band && s->picture_structure==PICT_FRAME))
  2747.     {
  2748.         if(s->picture_structure==PICT_FRAME){
  2749.         s->dest[0] += s->mb_y *   linesize << mb_size;
  2750.         s->dest[1] += s->mb_y * uvlinesize << (mb_size - s->chroma_y_shift);
  2751.         s->dest[2] += s->mb_y * uvlinesize << (mb_size - s->chroma_y_shift);
  2752.         }else{
  2753.             s->dest[0] += (s->mb_y>>1) *   linesize << mb_size;
  2754.             s->dest[1] += (s->mb_y>>1) * uvlinesize << (mb_size - s->chroma_y_shift);
  2755.             s->dest[2] += (s->mb_y>>1) * uvlinesize << (mb_size - s->chroma_y_shift);
  2756.             av_assert1((s->mb_y&1) == (s->picture_structure == PICT_BOTTOM_FIELD));
  2757.         }
  2758.     }
  2759. }
  2760.  
  2761. void ff_mpeg_flush(AVCodecContext *avctx){
  2762.     int i;
  2763.     MpegEncContext *s = avctx->priv_data;
  2764.  
  2765.     if (!s || !s->picture)
  2766.         return;
  2767.  
  2768.     for (i = 0; i < MAX_PICTURE_COUNT; i++)
  2769.         ff_mpeg_unref_picture(s->avctx, &s->picture[i]);
  2770.     s->current_picture_ptr = s->last_picture_ptr = s->next_picture_ptr = NULL;
  2771.  
  2772.     ff_mpeg_unref_picture(s->avctx, &s->current_picture);
  2773.     ff_mpeg_unref_picture(s->avctx, &s->last_picture);
  2774.     ff_mpeg_unref_picture(s->avctx, &s->next_picture);
  2775.  
  2776.     s->mb_x= s->mb_y= 0;
  2777.     s->closed_gop= 0;
  2778.  
  2779.     s->parse_context.state= -1;
  2780.     s->parse_context.frame_start_found= 0;
  2781.     s->parse_context.overread= 0;
  2782.     s->parse_context.overread_index= 0;
  2783.     s->parse_context.index= 0;
  2784.     s->parse_context.last_index= 0;
  2785.     s->bitstream_buffer_size=0;
  2786.     s->pp_time=0;
  2787. }
  2788.  
  2789. /**
  2790.  * set qscale and update qscale dependent variables.
  2791.  */
  2792. void ff_set_qscale(MpegEncContext * s, int qscale)
  2793. {
  2794.     if (qscale < 1)
  2795.         qscale = 1;
  2796.     else if (qscale > 31)
  2797.         qscale = 31;
  2798.  
  2799.     s->qscale = qscale;
  2800.     s->chroma_qscale= s->chroma_qscale_table[qscale];
  2801.  
  2802.     s->y_dc_scale= s->y_dc_scale_table[ qscale ];
  2803.     s->c_dc_scale= s->c_dc_scale_table[ s->chroma_qscale ];
  2804. }
  2805.  
  2806. void ff_mpv_report_decode_progress(MpegEncContext *s)
  2807. {
  2808.     if (s->pict_type != AV_PICTURE_TYPE_B && !s->partitioned_frame && !s->er.error_occurred)
  2809.         ff_thread_report_progress(&s->current_picture_ptr->tf, s->mb_y, 0);
  2810. }
  2811.