falkTX
/
FFmpeg
mirror of https://github.com/falkTX/FFmpeg.git
1
0
Fork 0
Code Issues 0 Releases 338 Wiki Activity
You can not select more than 25 topics Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
3195 Commits
32 Branches
484MB
Tree: cac0a56c55
Branches Tags
${ item.name }
Create branch ${ searchTerm }
from 'cac0a56c55'
${ noResults }
FFmpeg/libavcodec/motion_est.c

1889 lines
66KB
Raw Blame History 

  1. /*

  2.  * Motion estimation 

  3.  * Copyright (c) 2000,2001 Fabrice Bellard.

  4.  * Copyright (c) 2002-2004 Michael Niedermayer

  5.  * 

  6.  *

  7.  * This library is free software; you can redistribute it and/or

  8.  * modify it under the terms of the GNU Lesser General Public

  9.  * License as published by the Free Software Foundation; either

  10.  * version 2 of the License, or (at your option) any later version.

  11.  *

  12.  * This library is distributed in the hope that it will be useful,

  13.  * but WITHOUT ANY WARRANTY; without even the implied warranty of

  14.  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU

  15.  * Lesser General Public License for more details.

  16.  *

  17.  * You should have received a copy of the GNU Lesser General Public

  18.  * License along with this library; if not, write to the Free Software

  19.  * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA

  20.  *

  21.  * new Motion Estimation (X1/EPZS) by Michael Niedermayer <michaelni@gmx.at>

  22.  */

  23.  

  24. /**

  25.  * @file motion_est.c

  26.  * Motion estimation.

  27.  */

  28.  

  29. #include <stdlib.h>

  30. #include <stdio.h>

  31. #include <limits.h>

  32. #include "avcodec.h"

  33. #include "dsputil.h"

  34. #include "mpegvideo.h"

  35. 

  36. //#undef NDEBUG

  37. //#include <assert.h>

  38. 

  39. #define SQ(a) ((a)*(a))

  40. 

  41. #define P_LEFT P[1]

  42. #define P_TOP P[2]

  43. #define P_TOPRIGHT P[3]

  44. #define P_MEDIAN P[4]

  45. #define P_MV1 P[9]

  46. 

  47. static inline int sad_hpel_motion_search(MpegEncContext * s,

  48. 				  int *mx_ptr, int *my_ptr, int dmin,

  49.                                   int pred_x, int pred_y, uint8_t *src_data[3],

  50.                                   uint8_t *ref_data[6], int stride, int uvstride,

  51.                                   int size, int h, uint8_t * const mv_penalty);

  52. 

  53. static inline int update_map_generation(MpegEncContext * s)

  54. {

  55.     s->me.map_generation+= 1<<(ME_MAP_MV_BITS*2);

  56.     if(s->me.map_generation==0){

  57.         s->me.map_generation= 1<<(ME_MAP_MV_BITS*2);

  58.         memset(s->me.map, 0, sizeof(uint32_t)*ME_MAP_SIZE);

  59.     }

  60.     return s->me.map_generation;

  61. }

  62. 

  63. /* shape adaptive search stuff */

  64. typedef struct Minima{

  65.     int height;

  66.     int x, y;

  67.     int checked;

  68. }Minima;

  69. 

  70. static int minima_cmp(const void *a, const void *b){

  71.     const Minima *da = (const Minima *) a;

  72.     const Minima *db = (const Minima *) b;

  73.     

  74.     return da->height - db->height;

  75. }

  76.                                   

  77. /* SIMPLE */

  78. #define RENAME(a) simple_ ## a

  79. 

  80. #define CMP(d, x, y, size)\

  81. d = cmp(s, src_y, (ref_y) + (x) + (y)*(stride), stride, h);

  82. 

  83. #define CMP_HPEL(d, dx, dy, x, y, size)\

  84. {\

  85.     const int dxy= (dx) + 2*(dy);\

  86.     hpel_put[0][dxy](s->me.scratchpad, (ref_y) + (x) + (y)*(stride), stride, h);\

  87.     d = cmp_sub(s, s->me.scratchpad, src_y, stride, h);\

  88. }

  89. 

  90. 

  91. #define CMP_QPEL(d, dx, dy, x, y, size)\

  92. {\

  93.     const int dxy= (dx) + 4*(dy);\

  94.     qpel_put[0][dxy](s->me.scratchpad, (ref_y) + (x) + (y)*(stride), stride);\

  95.     d = cmp_sub(s, s->me.scratchpad, src_y, stride, h);\

  96. }

  97. 

  98. #include "motion_est_template.c"

  99. #undef RENAME

  100. #undef CMP

  101. #undef CMP_HPEL

  102. #undef CMP_QPEL

  103. #undef INIT

  104. 

  105. /* SIMPLE CHROMA */

  106. #define RENAME(a) simple_chroma_ ## a

  107. 

  108. #define CMP(d, x, y, size)\

  109. d = cmp(s, src_y, (ref_y) + (x) + (y)*(stride), stride, h);\

  110. if(chroma_cmp){\

  111.     int dxy= ((x)&1) + 2*((y)&1);\

  112.     int c= ((x)>>1) + ((y)>>1)*uvstride;\

  113. \

  114.     chroma_hpel_put[0][dxy](s->me.scratchpad, ref_u + c, uvstride, h>>1);\

  115.     d += chroma_cmp(s, s->me.scratchpad, src_u, uvstride, h>>1);\

  116.     chroma_hpel_put[0][dxy](s->me.scratchpad, ref_v + c, uvstride, h>>1);\

  117.     d += chroma_cmp(s, s->me.scratchpad, src_v, uvstride, h>>1);\

  118. }

  119. 

  120. #define CMP_HPEL(d, dx, dy, x, y, size)\

  121. {\

  122.     const int dxy= (dx) + 2*(dy);\

  123.     hpel_put[0][dxy](s->me.scratchpad, (ref_y) + (x) + (y)*(stride), stride, h);\

  124.     d = cmp_sub(s, s->me.scratchpad, src_y, stride, h);\

  125.     if(chroma_cmp_sub){\

  126.         int cxy= (dxy) | ((x)&1) | (2*((y)&1));\

  127.         int c= ((x)>>1) + ((y)>>1)*uvstride;\

  128.         chroma_hpel_put[0][cxy](s->me.scratchpad, ref_u + c, uvstride, h>>1);\

  129.         d += chroma_cmp_sub(s, s->me.scratchpad, src_u, uvstride, h>>1);\

  130.         chroma_hpel_put[0][cxy](s->me.scratchpad, ref_v + c, uvstride, h>>1);\

  131.         d += chroma_cmp_sub(s, s->me.scratchpad, src_v, uvstride, h>>1);\

  132.     }\

  133. }

  134. 

  135. #define CMP_QPEL(d, dx, dy, x, y, size)\

  136. {\

  137.     const int dxy= (dx) + 4*(dy);\

  138.     qpel_put[0][dxy](s->me.scratchpad, (ref_y) + (x) + (y)*(stride), stride);\

  139.     d = cmp_sub(s, s->me.scratchpad, src_y, stride, h);\

  140.     if(chroma_cmp_sub){\

  141.         int cxy, c;\

  142.         int cx= (4*(x) + (dx))/2;\

  143.         int cy= (4*(y) + (dy))/2;\

  144.         cx= (cx>>1)|(cx&1);\

  145.         cy= (cy>>1)|(cy&1);\

  146.         cxy= (cx&1) + 2*(cy&1);\

  147.         c= ((cx)>>1) + ((cy)>>1)*uvstride;\

  148.         chroma_hpel_put[0][cxy](s->me.scratchpad, ref_u + c, uvstride, h>>1);\

  149.         d += chroma_cmp_sub(s, s->me.scratchpad, src_u, uvstride, h>>1);\

  150.         chroma_hpel_put[0][cxy](s->me.scratchpad, ref_v + c, uvstride, h>>1);\

  151.         d += chroma_cmp_sub(s, s->me.scratchpad, src_v, uvstride, h>>1);\

  152.     }\

  153. }

  154. 

  155. #include "motion_est_template.c"

  156. #undef RENAME

  157. #undef CMP

  158. #undef CMP_HPEL

  159. #undef CMP_QPEL

  160. #undef INIT

  161. 

  162. /* SIMPLE DIRECT HPEL */

  163. #define RENAME(a) simple_direct_hpel_ ## a

  164. //FIXME precalc divisions stuff

  165. 

  166. #define CMP_DIRECT(d, dx, dy, x, y, size, cmp_func)\

  167. if((x) >= xmin && 2*(x) + (dx) <= 2*xmax && (y) >= ymin && 2*(y) + (dy) <= 2*ymax){\

  168.     const int hx= 2*(x) + (dx);\

  169.     const int hy= 2*(y) + (dy);\

  170.     if(s->mv_type==MV_TYPE_8X8){\

  171.         int i;\

  172.         for(i=0; i<4; i++){\

  173.             int fx = s->me.direct_basis_mv[i][0] + hx;\

  174.             int fy = s->me.direct_basis_mv[i][1] + hy;\

  175.             int bx = hx ? fx - s->me.co_located_mv[i][0] : s->me.co_located_mv[i][0]*(time_pb - time_pp)/time_pp + (i &1)*16;\

  176.             int by = hy ? fy - s->me.co_located_mv[i][1] : s->me.co_located_mv[i][1]*(time_pb - time_pp)/time_pp + (i>>1)*16;\

  177.             int fxy= (fx&1) + 2*(fy&1);\

  178.             int bxy= (bx&1) + 2*(by&1);\

  179. \

  180.             uint8_t *dst= s->me.scratchpad + 8*(i&1) + 8*stride*(i>>1);\

  181.             hpel_put[1][fxy](dst, (ref_y ) + (fx>>1) + (fy>>1)*(stride), stride, 8);\

  182.             hpel_avg[1][bxy](dst, (ref_data[3]) + (bx>>1) + (by>>1)*(stride), stride, 8);\

  183.         }\

  184.     }else{\

  185.         int fx = s->me.direct_basis_mv[0][0] + hx;\

  186.         int fy = s->me.direct_basis_mv[0][1] + hy;\

  187.         int bx = hx ? fx - s->me.co_located_mv[0][0] : (s->me.co_located_mv[0][0]*(time_pb - time_pp)/time_pp);\

  188.         int by = hy ? fy - s->me.co_located_mv[0][1] : (s->me.co_located_mv[0][1]*(time_pb - time_pp)/time_pp);\

  189.         int fxy= (fx&1) + 2*(fy&1);\

  190.         int bxy= (bx&1) + 2*(by&1);\

  191.         \

  192.         assert((fx>>1) + 16*s->mb_x >= -16);\

  193.         assert((fy>>1) + 16*s->mb_y >= -16);\

  194.         assert((fx>>1) + 16*s->mb_x <= s->width);\

  195.         assert((fy>>1) + 16*s->mb_y <= s->height);\

  196.         assert((bx>>1) + 16*s->mb_x >= -16);\

  197.         assert((by>>1) + 16*s->mb_y >= -16);\

  198.         assert((bx>>1) + 16*s->mb_x <= s->width);\

  199.         assert((by>>1) + 16*s->mb_y <= s->height);\

  200. \

  201.         hpel_put[0][fxy](s->me.scratchpad, (ref_y ) + (fx>>1) + (fy>>1)*(stride), stride, 16);\

  202.         hpel_avg[0][bxy](s->me.scratchpad, (ref_data[3]) + (bx>>1) + (by>>1)*(stride), stride, 16);\

  203.     }\

  204.     d = cmp_func(s, s->me.scratchpad, src_y, stride, 16);\

  205. }else\

  206.     d= 256*256*256*32;

  207. 

  208. 

  209. #define CMP_HPEL(d, dx, dy, x, y, size)\

  210.     CMP_DIRECT(d, dx, dy, x, y, size, cmp_sub)

  211. 

  212. #define CMP(d, x, y, size)\

  213.     CMP_DIRECT(d, 0, 0, x, y, size, cmp)

  214.     

  215. #include "motion_est_template.c"

  216. #undef RENAME

  217. #undef CMP

  218. #undef CMP_HPEL

  219. #undef CMP_QPEL

  220. #undef INIT

  221. #undef CMP_DIRECT

  222. 

  223. /* SIMPLE DIRECT QPEL */

  224. #define RENAME(a) simple_direct_qpel_ ## a

  225. 

  226. #define CMP_DIRECT(d, dx, dy, x, y, size, cmp_func)\

  227. if((x) >= xmin && 4*(x) + (dx) <= 4*xmax && (y) >= ymin && 4*(y) + (dy) <= 4*ymax){\

  228.     const int qx= 4*(x) + (dx);\

  229.     const int qy= 4*(y) + (dy);\

  230.     if(s->mv_type==MV_TYPE_8X8){\

  231.         int i;\

  232.         for(i=0; i<4; i++){\

  233.             int fx = s->me.direct_basis_mv[i][0] + qx;\

  234.             int fy = s->me.direct_basis_mv[i][1] + qy;\

  235.             int bx = qx ? fx - s->me.co_located_mv[i][0] : s->me.co_located_mv[i][0]*(time_pb - time_pp)/time_pp + (i &1)*16;\

  236.             int by = qy ? fy - s->me.co_located_mv[i][1] : s->me.co_located_mv[i][1]*(time_pb - time_pp)/time_pp + (i>>1)*16;\

  237.             int fxy= (fx&3) + 4*(fy&3);\

  238.             int bxy= (bx&3) + 4*(by&3);\

  239. \

  240.             uint8_t *dst= s->me.scratchpad + 8*(i&1) + 8*stride*(i>>1);\

  241.             qpel_put[1][fxy](dst, (ref_y ) + (fx>>2) + (fy>>2)*(stride), stride);\

  242.             qpel_avg[1][bxy](dst, (ref_data[3]) + (bx>>2) + (by>>2)*(stride), stride);\

  243.         }\

  244.     }else{\

  245.         int fx = s->me.direct_basis_mv[0][0] + qx;\

  246.         int fy = s->me.direct_basis_mv[0][1] + qy;\

  247.         int bx = qx ? fx - s->me.co_located_mv[0][0] : s->me.co_located_mv[0][0]*(time_pb - time_pp)/time_pp;\

  248.         int by = qy ? fy - s->me.co_located_mv[0][1] : s->me.co_located_mv[0][1]*(time_pb - time_pp)/time_pp;\

  249.         int fxy= (fx&3) + 4*(fy&3);\

  250.         int bxy= (bx&3) + 4*(by&3);\

  251. \

  252.         qpel_put[1][fxy](s->me.scratchpad               , (ref_y ) + (fx>>2) + (fy>>2)*(stride)               , stride);\

  253.         qpel_put[1][fxy](s->me.scratchpad + 8           , (ref_y ) + (fx>>2) + (fy>>2)*(stride) + 8           , stride);\

  254.         qpel_put[1][fxy](s->me.scratchpad     + 8*stride, (ref_y ) + (fx>>2) + (fy>>2)*(stride)     + 8*stride, stride);\

  255.         qpel_put[1][fxy](s->me.scratchpad + 8 + 8*stride, (ref_y ) + (fx>>2) + (fy>>2)*(stride) + 8 + 8*stride, stride);\

  256.         qpel_avg[1][bxy](s->me.scratchpad               , (ref_data[3]) + (bx>>2) + (by>>2)*(stride)               , stride);\

  257.         qpel_avg[1][bxy](s->me.scratchpad + 8           , (ref_data[3]) + (bx>>2) + (by>>2)*(stride) + 8           , stride);\

  258.         qpel_avg[1][bxy](s->me.scratchpad     + 8*stride, (ref_data[3]) + (bx>>2) + (by>>2)*(stride)     + 8*stride, stride);\

  259.         qpel_avg[1][bxy](s->me.scratchpad + 8 + 8*stride, (ref_data[3]) + (bx>>2) + (by>>2)*(stride) + 8 + 8*stride, stride);\

  260.     }\

  261.     d = cmp_func(s, s->me.scratchpad, src_y, stride, 16);\

  262. }else\

  263.     d= 256*256*256*32;

  264. 

  265. 

  266. #define CMP_QPEL(d, dx, dy, x, y, size)\

  267.     CMP_DIRECT(d, dx, dy, x, y, size, cmp_sub)

  268. 

  269. #define CMP(d, x, y, size)\

  270.     CMP_DIRECT(d, 0, 0, x, y, size, cmp)

  271. 

  272. #include "motion_est_template.c"

  273. #undef RENAME

  274. #undef CMP

  275. #undef CMP_HPEL

  276. #undef CMP_QPEL

  277. #undef INIT

  278. #undef CMP__DIRECT

  279. 

  280. static inline int get_penalty_factor(MpegEncContext *s, int type){

  281.     switch(type&0xFF){

  282.     default:

  283.     case FF_CMP_SAD:

  284.         return s->qscale*2;

  285.     case FF_CMP_DCT:

  286.         return s->qscale*3;

  287.     case FF_CMP_SATD:

  288.         return s->qscale*6;

  289.     case FF_CMP_SSE:

  290.         return s->qscale*s->qscale*2;

  291.     case FF_CMP_BIT:

  292.         return 1;

  293.     case FF_CMP_RD:

  294.     case FF_CMP_PSNR:

  295.         return (s->qscale*s->qscale*185 + 64)>>7;

  296.     }

  297. }

  298. 

  299. void ff_init_me(MpegEncContext *s){

  300.     ff_set_cmp(&s->dsp, s->dsp.me_pre_cmp, s->avctx->me_pre_cmp);

  301.     ff_set_cmp(&s->dsp, s->dsp.me_cmp, s->avctx->me_cmp);

  302.     ff_set_cmp(&s->dsp, s->dsp.me_sub_cmp, s->avctx->me_sub_cmp);

  303.     ff_set_cmp(&s->dsp, s->dsp.mb_cmp, s->avctx->mb_cmp);

  304. 

  305.     if(s->flags&CODEC_FLAG_QPEL){

  306.         if(s->avctx->me_sub_cmp&FF_CMP_CHROMA)

  307.             s->me.sub_motion_search= simple_chroma_qpel_motion_search;

  308.         else

  309.             s->me.sub_motion_search= simple_qpel_motion_search;

  310.     }else{

  311.         if(s->avctx->me_sub_cmp&FF_CMP_CHROMA)

  312.             s->me.sub_motion_search= simple_chroma_hpel_motion_search;

  313.         else if(   s->avctx->me_sub_cmp == FF_CMP_SAD 

  314.                 && s->avctx->    me_cmp == FF_CMP_SAD 

  315.                 && s->avctx->    mb_cmp == FF_CMP_SAD)

  316.             s->me.sub_motion_search= sad_hpel_motion_search; // 2050 vs. 2450 cycles

  317.         else

  318.             s->me.sub_motion_search= simple_hpel_motion_search;

  319.     }

  320. 

  321.     if(s->avctx->me_cmp&FF_CMP_CHROMA){

  322.         s->me.motion_search[0]= simple_chroma_epzs_motion_search;

  323.         s->me.motion_search[1]= simple_chroma_epzs_motion_search4;

  324.         s->me.motion_search[4]= simple_chroma_epzs_motion_search2;

  325.     }else{

  326.         s->me.motion_search[0]= simple_epzs_motion_search;

  327.         s->me.motion_search[1]= simple_epzs_motion_search4;

  328.         s->me.motion_search[4]= simple_epzs_motion_search2;

  329.     }

  330.     

  331.     if(s->avctx->me_pre_cmp&FF_CMP_CHROMA){

  332.         s->me.pre_motion_search= simple_chroma_epzs_motion_search;

  333.     }else{

  334.         s->me.pre_motion_search= simple_epzs_motion_search;

  335.     }

  336.     

  337.     if(s->flags&CODEC_FLAG_QPEL){

  338.         if(s->avctx->mb_cmp&FF_CMP_CHROMA)

  339.             s->me.get_mb_score= simple_chroma_qpel_get_mb_score;

  340.         else

  341.             s->me.get_mb_score= simple_qpel_get_mb_score;

  342.     }else{

  343.         if(s->avctx->mb_cmp&FF_CMP_CHROMA)

  344.             s->me.get_mb_score= simple_chroma_hpel_get_mb_score;

  345.         else

  346.             s->me.get_mb_score= simple_hpel_get_mb_score;

  347.     }

  348. }

  349.       

  350. #if 0

  351. static int pix_dev(uint8_t * pix, int line_size, int mean)

  352. {

  353.     int s, i, j;

  354. 

  355.     s = 0;

  356.     for (i = 0; i < 16; i++) {

  357. 	for (j = 0; j < 16; j += 8) {

  358. 	    s += ABS(pix[0]-mean);

  359. 	    s += ABS(pix[1]-mean);

  360. 	    s += ABS(pix[2]-mean);

  361. 	    s += ABS(pix[3]-mean);

  362. 	    s += ABS(pix[4]-mean);

  363. 	    s += ABS(pix[5]-mean);

  364. 	    s += ABS(pix[6]-mean);

  365. 	    s += ABS(pix[7]-mean);

  366. 	    pix += 8;

  367. 	}

  368. 	pix += line_size - 16;

  369.     }

  370.     return s;

  371. }

  372. #endif

  373. 

  374. static inline void no_motion_search(MpegEncContext * s,

  375. 				    int *mx_ptr, int *my_ptr)

  376. {

  377.     *mx_ptr = 16 * s->mb_x;

  378.     *my_ptr = 16 * s->mb_y;

  379. }

  380. 

  381. static int full_motion_search(MpegEncContext * s,

  382.                               int *mx_ptr, int *my_ptr, int range,

  383.                               int xmin, int ymin, int xmax, int ymax, uint8_t *ref_picture)

  384. {

  385.     int x1, y1, x2, y2, xx, yy, x, y;

  386.     int mx, my, dmin, d;

  387.     uint8_t *pix;

  388. 

  389.     xx = 16 * s->mb_x;

  390.     yy = 16 * s->mb_y;

  391.     x1 = xx - range + 1;	/* we loose one pixel to avoid boundary pb with half pixel pred */

  392.     if (x1 < xmin)

  393. 	x1 = xmin;

  394.     x2 = xx + range - 1;

  395.     if (x2 > xmax)

  396. 	x2 = xmax;

  397.     y1 = yy - range + 1;

  398.     if (y1 < ymin)

  399. 	y1 = ymin;

  400.     y2 = yy + range - 1;

  401.     if (y2 > ymax)

  402. 	y2 = ymax;

  403.     pix = s->new_picture.data[0] + (yy * s->linesize) + xx;

  404.     dmin = 0x7fffffff;

  405.     mx = 0;

  406.     my = 0;

  407.     for (y = y1; y <= y2; y++) {

  408. 	for (x = x1; x <= x2; x++) {

  409. 	    d = s->dsp.pix_abs[0][0](NULL, pix, ref_picture + (y * s->linesize) + x,

  410. 			     s->linesize, 16);

  411. 	    if (d < dmin ||

  412. 		(d == dmin &&

  413. 		 (abs(x - xx) + abs(y - yy)) <

  414. 		 (abs(mx - xx) + abs(my - yy)))) {

  415. 		dmin = d;

  416. 		mx = x;

  417. 		my = y;

  418. 	    }

  419. 	}

  420.     }

  421. 

  422.     *mx_ptr = mx;

  423.     *my_ptr = my;

  424. 

  425. #if 0

  426.     if (*mx_ptr < -(2 * range) || *mx_ptr >= (2 * range) ||

  427. 	*my_ptr < -(2 * range) || *my_ptr >= (2 * range)) {

  428. 	fprintf(stderr, "error %d %d\n", *mx_ptr, *my_ptr);

  429.     }

  430. #endif

  431.     return dmin;

  432. }

  433. 

  434. 

  435. static int log_motion_search(MpegEncContext * s,

  436.                              int *mx_ptr, int *my_ptr, int range,

  437.                              int xmin, int ymin, int xmax, int ymax, uint8_t *ref_picture)

  438. {

  439.     int x1, y1, x2, y2, xx, yy, x, y;

  440.     int mx, my, dmin, d;

  441.     uint8_t *pix;

  442. 

  443.     xx = s->mb_x << 4;

  444.     yy = s->mb_y << 4;

  445. 

  446.     /* Left limit */

  447.     x1 = xx - range;

  448.     if (x1 < xmin)

  449. 	x1 = xmin;

  450. 

  451.     /* Right limit */

  452.     x2 = xx + range;

  453.     if (x2 > xmax)

  454. 	x2 = xmax;

  455. 

  456.     /* Upper limit */

  457.     y1 = yy - range;

  458.     if (y1 < ymin)

  459. 	y1 = ymin;

  460. 

  461.     /* Lower limit */

  462.     y2 = yy + range;

  463.     if (y2 > ymax)

  464. 	y2 = ymax;

  465. 

  466.     pix = s->new_picture.data[0] + (yy * s->linesize) + xx;

  467.     dmin = 0x7fffffff;

  468.     mx = 0;

  469.     my = 0;

  470. 

  471.     do {

  472. 	for (y = y1; y <= y2; y += range) {

  473. 	    for (x = x1; x <= x2; x += range) {

  474. 		d = s->dsp.pix_abs[0][0](NULL, pix, ref_picture + (y * s->linesize) + x, s->linesize, 16);

  475. 		if (d < dmin || (d == dmin && (abs(x - xx) + abs(y - yy)) < (abs(mx - xx) + abs(my - yy)))) {

  476. 		    dmin = d;

  477. 		    mx = x;

  478. 		    my = y;

  479. 		}

  480. 	    }

  481. 	}

  482. 

  483. 	range = range >> 1;

  484. 

  485. 	x1 = mx - range;

  486. 	if (x1 < xmin)

  487. 	    x1 = xmin;

  488. 

  489. 	x2 = mx + range;

  490. 	if (x2 > xmax)

  491. 	    x2 = xmax;

  492. 

  493. 	y1 = my - range;

  494. 	if (y1 < ymin)

  495. 	    y1 = ymin;

  496. 

  497. 	y2 = my + range;

  498. 	if (y2 > ymax)

  499. 	    y2 = ymax;

  500. 

  501.     } while (range >= 1);

  502. 

  503. #ifdef DEBUG

  504.     fprintf(stderr, "log       - MX: %d\tMY: %d\n", mx, my);

  505. #endif

  506.     *mx_ptr = mx;

  507.     *my_ptr = my;

  508.     return dmin;

  509. }

  510. 

  511. static int phods_motion_search(MpegEncContext * s,

  512.                                int *mx_ptr, int *my_ptr, int range,

  513.                                int xmin, int ymin, int xmax, int ymax, uint8_t *ref_picture)

  514. {

  515.     int x1, y1, x2, y2, xx, yy, x, y, lastx, d;

  516.     int mx, my, dminx, dminy;

  517.     uint8_t *pix;

  518. 

  519.     xx = s->mb_x << 4;

  520.     yy = s->mb_y << 4;

  521. 

  522.     /* Left limit */

  523.     x1 = xx - range;

  524.     if (x1 < xmin)

  525. 	x1 = xmin;

  526. 

  527.     /* Right limit */

  528.     x2 = xx + range;

  529.     if (x2 > xmax)

  530. 	x2 = xmax;

  531. 

  532.     /* Upper limit */

  533.     y1 = yy - range;

  534.     if (y1 < ymin)

  535. 	y1 = ymin;

  536. 

  537.     /* Lower limit */

  538.     y2 = yy + range;

  539.     if (y2 > ymax)

  540. 	y2 = ymax;

  541. 

  542.     pix = s->new_picture.data[0] + (yy * s->linesize) + xx;

  543.     mx = 0;

  544.     my = 0;

  545. 

  546.     x = xx;

  547.     y = yy;

  548.     do {

  549.         dminx = 0x7fffffff;

  550.         dminy = 0x7fffffff;

  551. 

  552. 	lastx = x;

  553. 	for (x = x1; x <= x2; x += range) {

  554. 	    d = s->dsp.pix_abs[0][0](NULL, pix, ref_picture + (y * s->linesize) + x, s->linesize, 16);

  555. 	    if (d < dminx || (d == dminx && (abs(x - xx) + abs(y - yy)) < (abs(mx - xx) + abs(my - yy)))) {

  556. 		dminx = d;

  557. 		mx = x;

  558. 	    }

  559. 	}

  560. 

  561. 	x = lastx;

  562. 	for (y = y1; y <= y2; y += range) {

  563. 	    d = s->dsp.pix_abs[0][0](NULL, pix, ref_picture + (y * s->linesize) + x, s->linesize, 16);

  564. 	    if (d < dminy || (d == dminy && (abs(x - xx) + abs(y - yy)) < (abs(mx - xx) + abs(my - yy)))) {

  565. 		dminy = d;

  566. 		my = y;

  567. 	    }

  568. 	}

  569. 

  570. 	range = range >> 1;

  571. 

  572. 	x = mx;

  573. 	y = my;

  574. 	x1 = mx - range;

  575. 	if (x1 < xmin)

  576. 	    x1 = xmin;

  577. 

  578. 	x2 = mx + range;

  579. 	if (x2 > xmax)

  580. 	    x2 = xmax;

  581. 

  582. 	y1 = my - range;

  583. 	if (y1 < ymin)

  584. 	    y1 = ymin;

  585. 

  586. 	y2 = my + range;

  587. 	if (y2 > ymax)

  588. 	    y2 = ymax;

  589. 

  590.     } while (range >= 1);

  591. 

  592. #ifdef DEBUG

  593.     fprintf(stderr, "phods     - MX: %d\tMY: %d\n", mx, my);

  594. #endif

  595. 

  596.     /* half pixel search */

  597.     *mx_ptr = mx;

  598.     *my_ptr = my;

  599.     return dminy;

  600. }

  601. 

  602. 

  603. #define Z_THRESHOLD 256

  604. 

  605. #define CHECK_SAD_HALF_MV(suffix, x, y) \

  606. {\

  607.     d= s->dsp.pix_abs[size][(x?1:0)+(y?2:0)](NULL, pix, ptr+((x)>>1), stride, h);\

  608.     d += (mv_penalty[pen_x + x] + mv_penalty[pen_y + y])*penalty_factor;\

  609.     COPY3_IF_LT(dminh, d, dx, x, dy, y)\

  610. }

  611. 

  612. static inline int sad_hpel_motion_search(MpegEncContext * s,

  613. 				  int *mx_ptr, int *my_ptr, int dmin,

  614.                                   int pred_x, int pred_y, uint8_t *src_data[3],

  615.                                   uint8_t *ref_data[6], int stride, int uvstride,                                  

  616.                                   int size, int h, uint8_t * const mv_penalty)

  617. {

  618.     uint32_t *score_map= s->me.score_map;

  619.     const int penalty_factor= s->me.sub_penalty_factor;

  620.     int mx, my, dminh;

  621.     uint8_t *pix, *ptr;

  622.     const int xmin= s->me.xmin;

  623.     const int ymin= s->me.ymin;

  624.     const int xmax= s->me.xmax;

  625.     const int ymax= s->me.ymax;

  626. 

  627.     if(s->me.skip){

  628. //    printf("S");

  629.         *mx_ptr = 0;

  630.         *my_ptr = 0;

  631.         return dmin;

  632.     }

  633. //    printf("N");

  634.         

  635.     pix = src_data[0];

  636. 

  637.     mx = *mx_ptr;

  638.     my = *my_ptr;

  639.     ptr = ref_data[0] + (my * stride) + mx;

  640.     

  641.     dminh = dmin;

  642. 

  643.     if (mx > xmin && mx < xmax && 

  644.         my > ymin && my < ymax) {

  645.         int dx=0, dy=0;

  646.         int d, pen_x, pen_y; 

  647.         const int index= (my<<ME_MAP_SHIFT) + mx;

  648.         const int t= score_map[(index-(1<<ME_MAP_SHIFT))&(ME_MAP_SIZE-1)];

  649.         const int l= score_map[(index- 1               )&(ME_MAP_SIZE-1)];

  650.         const int r= score_map[(index+ 1               )&(ME_MAP_SIZE-1)];

  651.         const int b= score_map[(index+(1<<ME_MAP_SHIFT))&(ME_MAP_SIZE-1)];

  652.         mx<<=1;

  653.         my<<=1;

  654. 

  655.         

  656.         pen_x= pred_x + mx;

  657.         pen_y= pred_y + my;

  658. 

  659.         ptr-= stride;

  660.         if(t<=b){

  661.             CHECK_SAD_HALF_MV(y2 , 0, -1)

  662.             if(l<=r){

  663.                 CHECK_SAD_HALF_MV(xy2, -1, -1)

  664.                 if(t+r<=b+l){

  665.                     CHECK_SAD_HALF_MV(xy2, +1, -1)

  666.                     ptr+= stride;

  667.                 }else{

  668.                     ptr+= stride;

  669.                     CHECK_SAD_HALF_MV(xy2, -1, +1)

  670.                 }

  671.                 CHECK_SAD_HALF_MV(x2 , -1,  0)

  672.             }else{

  673.                 CHECK_SAD_HALF_MV(xy2, +1, -1)

  674.                 if(t+l<=b+r){

  675.                     CHECK_SAD_HALF_MV(xy2, -1, -1)

  676.                     ptr+= stride;

  677.                 }else{

  678.                     ptr+= stride;

  679.                     CHECK_SAD_HALF_MV(xy2, +1, +1)

  680.                 }

  681.                 CHECK_SAD_HALF_MV(x2 , +1,  0)

  682.             }

  683.         }else{

  684.             if(l<=r){

  685.                 if(t+l<=b+r){

  686.                     CHECK_SAD_HALF_MV(xy2, -1, -1)

  687.                     ptr+= stride;

  688.                 }else{

  689.                     ptr+= stride;

  690.                     CHECK_SAD_HALF_MV(xy2, +1, +1)

  691.                 }

  692.                 CHECK_SAD_HALF_MV(x2 , -1,  0)

  693.                 CHECK_SAD_HALF_MV(xy2, -1, +1)

  694.             }else{

  695.                 if(t+r<=b+l){

  696.                     CHECK_SAD_HALF_MV(xy2, +1, -1)

  697.                     ptr+= stride;

  698.                 }else{

  699.                     ptr+= stride;

  700.                     CHECK_SAD_HALF_MV(xy2, -1, +1)

  701.                 }

  702.                 CHECK_SAD_HALF_MV(x2 , +1,  0)

  703.                 CHECK_SAD_HALF_MV(xy2, +1, +1)

  704.             }

  705.             CHECK_SAD_HALF_MV(y2 ,  0, +1)

  706.         }

  707.         mx+=dx;

  708.         my+=dy;

  709. 

  710.     }else{

  711.         mx<<=1;

  712.         my<<=1;

  713.     }

  714. 

  715.     *mx_ptr = mx;

  716.     *my_ptr = my;

  717.     return dminh;

  718. }

  719. 

  720. static inline void set_p_mv_tables(MpegEncContext * s, int mx, int my, int mv4)

  721. {

  722.     const int xy= s->mb_x + s->mb_y*s->mb_stride;

  723.     

  724.     s->p_mv_table[xy][0] = mx;

  725.     s->p_mv_table[xy][1] = my;

  726. 

  727.     /* has allready been set to the 4 MV if 4MV is done */

  728.     if(mv4){

  729.         int mot_xy= s->block_index[0];

  730. 

  731.         s->current_picture.motion_val[0][mot_xy  ][0]= mx;

  732.         s->current_picture.motion_val[0][mot_xy  ][1]= my;

  733.         s->current_picture.motion_val[0][mot_xy+1][0]= mx;

  734.         s->current_picture.motion_val[0][mot_xy+1][1]= my;

  735. 

  736.         mot_xy += s->block_wrap[0];

  737.         s->current_picture.motion_val[0][mot_xy  ][0]= mx;

  738.         s->current_picture.motion_val[0][mot_xy  ][1]= my;

  739.         s->current_picture.motion_val[0][mot_xy+1][0]= mx;

  740.         s->current_picture.motion_val[0][mot_xy+1][1]= my;

  741.     }

  742. }

  743. 

  744. /**

  745.  * get fullpel ME search limits.

  746.  */

  747. static inline void get_limits(MpegEncContext *s, int x, int y)

  748. {

  749. /*

  750.     if(s->avctx->me_range) s->me.range= s->avctx->me_range >> 1;

  751.     else                   s->me.range= 16;

  752. */

  753.     if (s->unrestricted_mv) {

  754.         s->me.xmin = - x - 16;

  755.         s->me.ymin = - y - 16;

  756.         s->me.xmax = - x + s->mb_width *16;

  757.         s->me.ymax = - y + s->mb_height*16;

  758.     } else {

  759.         s->me.xmin = - x;

  760.         s->me.ymin = - y;

  761.         s->me.xmax = - x + s->mb_width *16 - 16;

  762.         s->me.ymax = - y + s->mb_height*16 - 16;

  763.     }

  764. }

  765. 

  766. static inline int h263_mv4_search(MpegEncContext *s, int mx, int my, int shift)

  767. {

  768.     const int size= 1;

  769.     const int h=8;

  770.     int block;

  771.     int P[10][2];

  772.     int dmin_sum=0, mx4_sum=0, my4_sum=0;

  773.     uint8_t * const mv_penalty= s->me.mv_penalty[s->f_code] + MAX_MV;

  774.     int same=1;

  775.     const int stride= s->linesize;

  776.     const int uvstride= s->uvlinesize;

  777. 

  778.     for(block=0; block<4; block++){

  779.         int mx4, my4;

  780.         int pred_x4, pred_y4;

  781.         int dmin4;

  782.         static const int off[4]= {2, 1, 1, -1};

  783.         const int mot_stride = s->block_wrap[0];

  784.         const int mot_xy = s->block_index[block];

  785.         const int block_x= (block&1);

  786.         const int block_y= (block>>1);

  787.         uint8_t *src_data[3]= {

  788.             s->new_picture.data[0] + 8*(2*s->mb_x + block_x) + stride  *8*(2*s->mb_y + block_y), //FIXME chroma?

  789.             s->new_picture.data[1] + 4*(2*s->mb_x + block_x) + uvstride*4*(2*s->mb_y + block_y),

  790.             s->new_picture.data[2] + 4*(2*s->mb_x + block_x) + uvstride*4*(2*s->mb_y + block_y)

  791.         };

  792.         uint8_t *ref_data[3]= {

  793.             s->last_picture.data[0] + 8*(2*s->mb_x + block_x) + stride  *8*(2*s->mb_y + block_y), //FIXME chroma?

  794.             s->last_picture.data[1] + 4*(2*s->mb_x + block_x) + uvstride*4*(2*s->mb_y + block_y),

  795.             s->last_picture.data[2] + 4*(2*s->mb_x + block_x) + uvstride*4*(2*s->mb_y + block_y)

  796.         };

  797. 

  798.         P_LEFT[0] = s->current_picture.motion_val[0][mot_xy - 1][0];

  799.         P_LEFT[1] = s->current_picture.motion_val[0][mot_xy - 1][1];

  800. 

  801.         if(P_LEFT[0]       > (s->me.xmax<<shift)) P_LEFT[0]       = (s->me.xmax<<shift);

  802. 

  803.         /* special case for first line */

  804.         if (s->first_slice_line && block<2) {

  805.             pred_x4= P_LEFT[0];

  806.             pred_y4= P_LEFT[1];

  807.         } else {

  808.             P_TOP[0]      = s->current_picture.motion_val[0][mot_xy - mot_stride             ][0];

  809.             P_TOP[1]      = s->current_picture.motion_val[0][mot_xy - mot_stride             ][1];

  810.             P_TOPRIGHT[0] = s->current_picture.motion_val[0][mot_xy - mot_stride + off[block]][0];

  811.             P_TOPRIGHT[1] = s->current_picture.motion_val[0][mot_xy - mot_stride + off[block]][1];

  812.             if(P_TOP[1]      > (s->me.ymax<<shift)) P_TOP[1]     = (s->me.ymax<<shift);

  813.             if(P_TOPRIGHT[0] < (s->me.xmin<<shift)) P_TOPRIGHT[0]= (s->me.xmin<<shift);

  814.             if(P_TOPRIGHT[0] > (s->me.xmax<<shift)) P_TOPRIGHT[0]= (s->me.xmax<<shift);

  815.             if(P_TOPRIGHT[1] > (s->me.ymax<<shift)) P_TOPRIGHT[1]= (s->me.ymax<<shift);

  816.     

  817.             P_MEDIAN[0]= mid_pred(P_LEFT[0], P_TOP[0], P_TOPRIGHT[0]);

  818.             P_MEDIAN[1]= mid_pred(P_LEFT[1], P_TOP[1], P_TOPRIGHT[1]);

  819. 

  820. //            if(s->out_format == FMT_H263){

  821.                 pred_x4 = P_MEDIAN[0];

  822.                 pred_y4 = P_MEDIAN[1];

  823. #if 0

  824.             }else { /* mpeg1 at least */

  825.                 pred_x4= P_LEFT[0];

  826.                 pred_y4= P_LEFT[1];

  827.             }

  828. #endif

  829.         }

  830.         P_MV1[0]= mx;

  831.         P_MV1[1]= my;

  832. 

  833.         dmin4 = s->me.motion_search[1](s, &mx4, &my4, P, pred_x4, pred_y4, 

  834.                                        src_data, ref_data, stride, uvstride, s->p_mv_table, (1<<16)>>shift, mv_penalty);

  835. 

  836.         dmin4= s->me.sub_motion_search(s, &mx4, &my4, dmin4, 

  837. 					  pred_x4, pred_y4, src_data, ref_data, stride, uvstride, size, h, mv_penalty);

  838.         

  839.         if(s->dsp.me_sub_cmp[0] != s->dsp.mb_cmp[0]

  840.            && s->avctx->mb_decision == FF_MB_DECISION_SIMPLE){

  841.             int dxy;

  842.             const int offset= ((block&1) + (block>>1)*stride)*8;

  843.             uint8_t *dest_y = s->me.scratchpad + offset;

  844.             if(s->quarter_sample){

  845.                 uint8_t *ref= ref_data[0] + (mx4>>2) + (my4>>2)*stride;

  846.                 dxy = ((my4 & 3) << 2) | (mx4 & 3);

  847. 

  848.                 if(s->no_rounding)

  849.                     s->dsp.put_no_rnd_qpel_pixels_tab[1][dxy](dest_y   , ref    , stride);

  850.                 else

  851.                     s->dsp.put_qpel_pixels_tab       [1][dxy](dest_y   , ref    , stride);

  852.             }else{

  853.                 uint8_t *ref= ref_data[0] + (mx4>>1) + (my4>>1)*stride;

  854.                 dxy = ((my4 & 1) << 1) | (mx4 & 1);

  855. 

  856.                 if(s->no_rounding)

  857.                     s->dsp.put_no_rnd_pixels_tab[1][dxy](dest_y    , ref    , stride, h);

  858.                 else

  859.                     s->dsp.put_pixels_tab       [1][dxy](dest_y    , ref    , stride, h);

  860.             }

  861.             dmin_sum+= (mv_penalty[mx4-pred_x4] + mv_penalty[my4-pred_y4])*s->me.mb_penalty_factor;

  862.         }else

  863.             dmin_sum+= dmin4;

  864. 

  865.         if(s->quarter_sample){

  866.             mx4_sum+= mx4/2;

  867.             my4_sum+= my4/2;

  868.         }else{

  869.             mx4_sum+= mx4;

  870.             my4_sum+= my4;

  871.         }

  872.             

  873.         s->current_picture.motion_val[0][ s->block_index[block] ][0]= mx4;

  874.         s->current_picture.motion_val[0][ s->block_index[block] ][1]= my4;

  875. 

  876.         if(mx4 != mx || my4 != my) same=0;

  877.     }

  878.     

  879.     if(same)

  880.         return INT_MAX;

  881.     

  882.     if(s->dsp.me_sub_cmp[0] != s->dsp.mb_cmp[0]){

  883.         dmin_sum += s->dsp.mb_cmp[0](s, s->new_picture.data[0] + s->mb_x*16 + s->mb_y*16*stride, s->me.scratchpad, stride, 16);

  884.     }

  885.     

  886.     if(s->avctx->mb_cmp&FF_CMP_CHROMA){

  887.         int dxy;

  888.         int mx, my;

  889.         int offset;

  890. 

  891.         mx= ff_h263_round_chroma(mx4_sum);

  892.         my= ff_h263_round_chroma(my4_sum);

  893.         dxy = ((my & 1) << 1) | (mx & 1);

  894.         

  895.         offset= (s->mb_x*8 + (mx>>1)) + (s->mb_y*8 + (my>>1))*s->uvlinesize;

  896.        

  897.         if(s->no_rounding){

  898.             s->dsp.put_no_rnd_pixels_tab[1][dxy](s->me.scratchpad    , s->last_picture.data[1] + offset, s->uvlinesize, 8);

  899.             s->dsp.put_no_rnd_pixels_tab[1][dxy](s->me.scratchpad+8  , s->last_picture.data[2] + offset, s->uvlinesize, 8);

  900.         }else{

  901.             s->dsp.put_pixels_tab       [1][dxy](s->me.scratchpad    , s->last_picture.data[1] + offset, s->uvlinesize, 8);

  902.             s->dsp.put_pixels_tab       [1][dxy](s->me.scratchpad+8  , s->last_picture.data[2] + offset, s->uvlinesize, 8);

  903.         }

  904. 

  905.         dmin_sum += s->dsp.mb_cmp[1](s, s->new_picture.data[1] + s->mb_x*8 + s->mb_y*8*s->uvlinesize, s->me.scratchpad  , s->uvlinesize, 8);

  906.         dmin_sum += s->dsp.mb_cmp[1](s, s->new_picture.data[2] + s->mb_x*8 + s->mb_y*8*s->uvlinesize, s->me.scratchpad+8, s->uvlinesize, 8);

  907.     }

  908. 

  909.     switch(s->avctx->mb_cmp&0xFF){

  910.     /*case FF_CMP_SSE:

  911.         return dmin_sum+ 32*s->qscale*s->qscale;*/

  912.     case FF_CMP_RD:

  913.         return dmin_sum;

  914.     default:

  915.         return dmin_sum+ 11*s->me.mb_penalty_factor;

  916.     }

  917. }

  918. 

  919. static int interlaced_search(MpegEncContext *s, uint8_t *frame_src_data[3], uint8_t *frame_ref_data[3], 

  920.                              int16_t (*mv_tables[2][2])[2], uint8_t *field_select_tables[2], int f_code, int mx, int my)

  921. {

  922.     const int size=0;

  923.     const int h=8;

  924.     int block;

  925.     int P[10][2];

  926.     uint8_t * const mv_penalty= s->me.mv_penalty[f_code] + MAX_MV;

  927.     int same=1;

  928.     const int stride= 2*s->linesize;

  929.     const int uvstride= 2*s->uvlinesize;

  930.     int dmin_sum= 0;

  931.     const int mot_stride= s->mb_stride;

  932.     const int xy= s->mb_x + s->mb_y*mot_stride;

  933.     

  934.     s->me.ymin>>=1;

  935.     s->me.ymax>>=1;

  936.     

  937.     for(block=0; block<2; block++){

  938.         int field_select;

  939.         int best_dmin= INT_MAX;

  940.         int best_field= -1;

  941. 

  942.         uint8_t *src_data[3]= {

  943.             frame_src_data[0] + s->  linesize*block,

  944.             frame_src_data[1] + s->uvlinesize*block,

  945.             frame_src_data[2] + s->uvlinesize*block

  946.         };

  947. 

  948.         for(field_select=0; field_select<2; field_select++){

  949.             int dmin, mx_i, my_i, pred_x, pred_y;

  950.             uint8_t *ref_data[3]= {

  951.                 frame_ref_data[0] + s->  linesize*field_select,

  952.                 frame_ref_data[1] + s->uvlinesize*field_select,

  953.                 frame_ref_data[2] + s->uvlinesize*field_select

  954.             };

  955.             int16_t (*mv_table)[2]= mv_tables[block][field_select];

  956.             

  957.             P_LEFT[0] = mv_table[xy - 1][0];

  958.             P_LEFT[1] = mv_table[xy - 1][1];

  959.             if(P_LEFT[0]       > (s->me.xmax<<1)) P_LEFT[0]       = (s->me.xmax<<1);

  960.             

  961.             pred_x= P_LEFT[0];

  962.             pred_y= P_LEFT[1];

  963.             

  964.             if(!s->first_slice_line){

  965.                 P_TOP[0]      = mv_table[xy - mot_stride][0];

  966.                 P_TOP[1]      = mv_table[xy - mot_stride][1];

  967.                 P_TOPRIGHT[0] = mv_table[xy - mot_stride + 1][0];

  968.                 P_TOPRIGHT[1] = mv_table[xy - mot_stride + 1][1];

  969.                 if(P_TOP[1]      > (s->me.ymax<<1)) P_TOP[1]     = (s->me.ymax<<1);

  970.                 if(P_TOPRIGHT[0] < (s->me.xmin<<1)) P_TOPRIGHT[0]= (s->me.xmin<<1);

  971.                 if(P_TOPRIGHT[0] > (s->me.xmax<<1)) P_TOPRIGHT[0]= (s->me.xmax<<1);

  972.                 if(P_TOPRIGHT[1] > (s->me.ymax<<1)) P_TOPRIGHT[1]= (s->me.ymax<<1);

  973.     

  974.                 P_MEDIAN[0]= mid_pred(P_LEFT[0], P_TOP[0], P_TOPRIGHT[0]);

  975.                 P_MEDIAN[1]= mid_pred(P_LEFT[1], P_TOP[1], P_TOPRIGHT[1]);

  976.             }

  977.             P_MV1[0]= mx; //FIXME not correct if block != field_select

  978.             P_MV1[1]= my / 2;

  979.             

  980.             dmin = s->me.motion_search[4](s, &mx_i, &my_i, P, pred_x, pred_y, 

  981.                                            src_data, ref_data, stride, uvstride, mv_table, (1<<16)>>1, mv_penalty);

  982. 

  983.             dmin= s->me.sub_motion_search(s, &mx_i, &my_i, dmin, 

  984.                                            pred_x, pred_y, src_data, ref_data, stride, uvstride, size, h, mv_penalty);

  985.             

  986.             mv_table[xy][0]= mx_i;

  987.             mv_table[xy][1]= my_i;

  988.             

  989.             if(s->dsp.me_sub_cmp[0] != s->dsp.mb_cmp[0]

  990.                && s->avctx->mb_decision == FF_MB_DECISION_SIMPLE){

  991.                 int dxy;

  992. 

  993.                 //FIXME chroma ME

  994.                 uint8_t *ref= ref_data[0] + (mx_i>>1) + (my_i>>1)*stride;

  995.                 dxy = ((my_i & 1) << 1) | (mx_i & 1);

  996. 

  997.                 if(s->no_rounding){

  998.                     s->dsp.put_no_rnd_pixels_tab[size][dxy](s->me.scratchpad, ref    , stride, h);

  999.                 }else{

  1000.                     s->dsp.put_pixels_tab       [size][dxy](s->me.scratchpad, ref    , stride, h);

  1001.                 }

  1002.                 dmin= s->dsp.mb_cmp[size](s, src_data[0], s->me.scratchpad, stride, h);

  1003.                 dmin+= (mv_penalty[mx_i-pred_x] + mv_penalty[my_i-pred_y] + 1)*s->me.mb_penalty_factor;

  1004.             }else

  1005.                 dmin+= s->me.mb_penalty_factor; //field_select bits

  1006.                 

  1007.             dmin += field_select != block; //slightly prefer same field

  1008.             

  1009.             if(dmin < best_dmin){

  1010.                 best_dmin= dmin;

  1011.                 best_field= field_select;

  1012.             }

  1013.         }

  1014.         {

  1015.             int16_t (*mv_table)[2]= mv_tables[block][best_field];

  1016. 

  1017.             if(mv_table[xy][0] != mx) same=0; //FIXME check if these checks work and are any good at all

  1018.             if(mv_table[xy][1]&1) same=0;

  1019.             if(mv_table[xy][1]*2 != my) same=0; 

  1020.             if(best_field != block) same=0;

  1021.         }

  1022. 

  1023.         field_select_tables[block][xy]= best_field;

  1024.         dmin_sum += best_dmin;

  1025.     }

  1026.     

  1027.     s->me.ymin<<=1;

  1028.     s->me.ymax<<=1;

  1029. 

  1030.     if(same)

  1031.         return INT_MAX;

  1032.     

  1033.     switch(s->avctx->mb_cmp&0xFF){

  1034.     /*case FF_CMP_SSE:

  1035.         return dmin_sum+ 32*s->qscale*s->qscale;*/

  1036.     case FF_CMP_RD:

  1037.         return dmin_sum;

  1038.     default:

  1039.         return dmin_sum+ 11*s->me.mb_penalty_factor;

  1040.     }

  1041. }

  1042. 

  1043. void ff_estimate_p_frame_motion(MpegEncContext * s,

  1044.                                 int mb_x, int mb_y)

  1045. {

  1046.     uint8_t *pix, *ppix;

  1047.     int sum, varc, vard, mx, my, dmin, xx, yy;

  1048.     int pred_x=0, pred_y=0;

  1049.     int P[10][2];

  1050.     const int shift= 1+s->quarter_sample;

  1051.     int mb_type=0;

  1052.     uint8_t *ref_picture= s->last_picture.data[0];

  1053.     Picture * const pic= &s->current_picture;

  1054.     uint8_t * const mv_penalty= s->me.mv_penalty[s->f_code] + MAX_MV;

  1055.     const int stride= s->linesize;

  1056.     const int uvstride= s->uvlinesize;

  1057.     uint8_t *src_data[3]= {

  1058.         s->new_picture.data[0] + 16*(mb_x + stride*mb_y),

  1059.         s->new_picture.data[1] + 8*(mb_x + uvstride*mb_y),

  1060.         s->new_picture.data[2] + 8*(mb_x + uvstride*mb_y)

  1061.     };

  1062.     uint8_t *ref_data[3]= {

  1063.         s->last_picture.data[0] + 16*(mb_x + stride*mb_y),

  1064.         s->last_picture.data[1] + 8*(mb_x + uvstride*mb_y),

  1065.         s->last_picture.data[2] + 8*(mb_x + uvstride*mb_y)

  1066.     };

  1067. 

  1068.     assert(s->quarter_sample==0 || s->quarter_sample==1);

  1069. 

  1070.     s->me.penalty_factor    = get_penalty_factor(s, s->avctx->me_cmp);

  1071.     s->me.sub_penalty_factor= get_penalty_factor(s, s->avctx->me_sub_cmp);

  1072.     s->me.mb_penalty_factor = get_penalty_factor(s, s->avctx->mb_cmp);

  1073. 

  1074.     get_limits(s, 16*mb_x, 16*mb_y);

  1075.     s->me.skip=0;

  1076. 

  1077.     switch(s->me_method) {

  1078.     case ME_ZERO:

  1079.     default:

  1080. 	no_motion_search(s, &mx, &my);

  1081.         mx-= mb_x*16;

  1082.         my-= mb_y*16;

  1083.         dmin = 0;

  1084.         break;

  1085. #if 0

  1086.     case ME_FULL:

  1087. 	dmin = full_motion_search(s, &mx, &my, range, ref_picture);

  1088.         mx-= mb_x*16;

  1089.         my-= mb_y*16;

  1090.         break;

  1091.     case ME_LOG:

  1092. 	dmin = log_motion_search(s, &mx, &my, range / 2, ref_picture);

  1093.         mx-= mb_x*16;

  1094.         my-= mb_y*16;

  1095.         break;

  1096.     case ME_PHODS:

  1097. 	dmin = phods_motion_search(s, &mx, &my, range / 2, ref_picture);

  1098.         mx-= mb_x*16;

  1099.         my-= mb_y*16;

  1100.         break;

  1101. #endif

  1102.     case ME_X1:

  1103.     case ME_EPZS:

  1104.        {

  1105.             const int mot_stride = s->block_wrap[0];

  1106.             const int mot_xy = s->block_index[0];

  1107. 

  1108.             P_LEFT[0]       = s->current_picture.motion_val[0][mot_xy - 1][0];

  1109.             P_LEFT[1]       = s->current_picture.motion_val[0][mot_xy - 1][1];

  1110. 

  1111.             if(P_LEFT[0]       > (s->me.xmax<<shift)) P_LEFT[0]       = (s->me.xmax<<shift);

  1112. 

  1113.             if(!s->first_slice_line) {

  1114.                 P_TOP[0]      = s->current_picture.motion_val[0][mot_xy - mot_stride    ][0];

  1115.                 P_TOP[1]      = s->current_picture.motion_val[0][mot_xy - mot_stride    ][1];

  1116.                 P_TOPRIGHT[0] = s->current_picture.motion_val[0][mot_xy - mot_stride + 2][0];

  1117.                 P_TOPRIGHT[1] = s->current_picture.motion_val[0][mot_xy - mot_stride + 2][1];

  1118.                 if(P_TOP[1]      > (s->me.ymax<<shift)) P_TOP[1]     = (s->me.ymax<<shift);

  1119.                 if(P_TOPRIGHT[0] < (s->me.xmin<<shift)) P_TOPRIGHT[0]= (s->me.xmin<<shift);

  1120.                 if(P_TOPRIGHT[1] > (s->me.ymax<<shift)) P_TOPRIGHT[1]= (s->me.ymax<<shift);

  1121.         

  1122.                 P_MEDIAN[0]= mid_pred(P_LEFT[0], P_TOP[0], P_TOPRIGHT[0]);

  1123.                 P_MEDIAN[1]= mid_pred(P_LEFT[1], P_TOP[1], P_TOPRIGHT[1]);

  1124. 

  1125.                 if(s->out_format == FMT_H263){

  1126.                     pred_x = P_MEDIAN[0];

  1127.                     pred_y = P_MEDIAN[1];

  1128.                 }else { /* mpeg1 at least */

  1129.                     pred_x= P_LEFT[0];

  1130.                     pred_y= P_LEFT[1];

  1131.                 }

  1132.             }else{

  1133.                 pred_x= P_LEFT[0];

  1134.                 pred_y= P_LEFT[1];

  1135.             }

  1136. 

  1137.         }

  1138.         dmin = s->me.motion_search[0](s, &mx, &my, P, pred_x, pred_y, 

  1139.                                       src_data, ref_data, stride, uvstride, s->p_mv_table, (1<<16)>>shift, mv_penalty);

  1140.  

  1141.         break;

  1142.     }

  1143. 

  1144.     /* intra / predictive decision */

  1145.     xx = mb_x * 16;

  1146.     yy = mb_y * 16;

  1147. 

  1148.     pix = src_data[0];

  1149.     /* At this point (mx,my) are full-pell and the relative displacement */

  1150.     ppix = ref_data[0] + (my * s->linesize) + mx;

  1151.     

  1152.     sum = s->dsp.pix_sum(pix, s->linesize);

  1153.     

  1154.     varc = (s->dsp.pix_norm1(pix, s->linesize) - (((unsigned)(sum*sum))>>8) + 500 + 128)>>8;

  1155.     vard = (s->dsp.sse[0](NULL, pix, ppix, s->linesize, 16)+128)>>8;

  1156. 

  1157. //printf("%d %d %d %X %X %X\n", s->mb_width, mb_x, mb_y,(int)s, (int)s->mb_var, (int)s->mc_mb_var); fflush(stdout);

  1158.     pic->mb_var   [s->mb_stride * mb_y + mb_x] = varc;

  1159.     pic->mc_mb_var[s->mb_stride * mb_y + mb_x] = vard;

  1160.     pic->mb_mean  [s->mb_stride * mb_y + mb_x] = (sum+128)>>8;

  1161. //    pic->mb_cmp_score[s->mb_stride * mb_y + mb_x] = dmin; 

  1162.     s->mb_var_sum_temp    += varc;

  1163.     s->mc_mb_var_sum_temp += vard;

  1164. //printf("E%d %d %d %X %X %X\n", s->mb_width, mb_x, mb_y,(int)s, (int)s->mb_var, (int)s->mc_mb_var); fflush(stdout);

  1165.     

  1166. #if 0

  1167.     printf("varc=%4d avg_var=%4d (sum=%4d) vard=%4d mx=%2d my=%2d\n",

  1168. 	   varc, s->avg_mb_var, sum, vard, mx - xx, my - yy);

  1169. #endif

  1170.     if(s->avctx->mb_decision > FF_MB_DECISION_SIMPLE){

  1171.         if (vard <= 64 || vard < varc)

  1172.             s->scene_change_score+= ff_sqrt(vard) - ff_sqrt(varc);

  1173.         else

  1174.             s->scene_change_score+= s->qscale;

  1175. 

  1176.         if (vard*2 + 200 > varc)

  1177.             mb_type|= CANDIDATE_MB_TYPE_INTRA;

  1178.         if (varc*2 + 200 > vard){

  1179.             mb_type|= CANDIDATE_MB_TYPE_INTER;

  1180.             s->me.sub_motion_search(s, &mx, &my, dmin,

  1181. 				   pred_x, pred_y, src_data, ref_data, stride, uvstride, 0, 16, mv_penalty);

  1182.             if(s->flags&CODEC_FLAG_MV0)

  1183.                 if(mx || my)

  1184.                     mb_type |= CANDIDATE_MB_TYPE_SKIPED; //FIXME check difference

  1185.         }else{

  1186.             mx <<=shift;

  1187.             my <<=shift;

  1188.         }

  1189.         if((s->flags&CODEC_FLAG_4MV)

  1190.            && !s->me.skip && varc>50 && vard>10){

  1191.             if(h263_mv4_search(s, mx, my, shift) < INT_MAX)

  1192.                 mb_type|=CANDIDATE_MB_TYPE_INTER4V;

  1193. 

  1194.             set_p_mv_tables(s, mx, my, 0);

  1195.         }else

  1196.             set_p_mv_tables(s, mx, my, 1);

  1197.         if((s->flags&CODEC_FLAG_INTERLACED_ME)

  1198.            && !s->me.skip){ //FIXME varc/d checks

  1199.             if(interlaced_search(s, src_data, ref_data, s->p_field_mv_table, s->p_field_select_table, s->f_code, mx, my) < INT_MAX)

  1200.                 mb_type |= CANDIDATE_MB_TYPE_INTER_I;

  1201.         }

  1202.     }else{

  1203.         int intra_score, i;

  1204.         mb_type= CANDIDATE_MB_TYPE_INTER;

  1205. 

  1206.         dmin= s->me.sub_motion_search(s, &mx, &my, dmin,

  1207.                                     pred_x, pred_y, src_data, ref_data, stride, uvstride, 0, 16, mv_penalty);

  1208.         if(s->avctx->me_sub_cmp != s->avctx->mb_cmp && !s->me.skip)

  1209.             dmin= s->me.get_mb_score(s, mx, my, pred_x, pred_y, src_data, ref_data, stride, uvstride, mv_penalty);

  1210. 

  1211.         if((s->flags&CODEC_FLAG_4MV)

  1212.            && !s->me.skip && varc>50 && vard>10){

  1213.             int dmin4= h263_mv4_search(s, mx, my, shift);

  1214.             if(dmin4 < dmin){

  1215.                 mb_type= CANDIDATE_MB_TYPE_INTER4V;

  1216.                 dmin=dmin4;

  1217.             }

  1218.         }

  1219.         if((s->flags&CODEC_FLAG_INTERLACED_ME)

  1220.            && !s->me.skip){ //FIXME varc/d checks

  1221.             int dmin_i= interlaced_search(s, src_data, ref_data, s->p_field_mv_table, s->p_field_select_table, s->f_code, mx, my);

  1222.             if(dmin_i < dmin){

  1223.                 mb_type = CANDIDATE_MB_TYPE_INTER_I;

  1224.                 dmin= dmin_i;

  1225.             }

  1226.         }

  1227.                 

  1228. //        pic->mb_cmp_score[s->mb_stride * mb_y + mb_x] = dmin; 

  1229.         set_p_mv_tables(s, mx, my, mb_type!=CANDIDATE_MB_TYPE_INTER4V);

  1230. 

  1231.         /* get intra luma score */

  1232.         if((s->avctx->mb_cmp&0xFF)==FF_CMP_SSE){

  1233.             intra_score= (varc<<8) - 500; //FIXME dont scale it down so we dont have to fix it

  1234.         }else{

  1235.             int mean= (sum+128)>>8;

  1236.             mean*= 0x01010101;

  1237.             

  1238.             for(i=0; i<16; i++){

  1239.                 *(uint32_t*)(&s->me.scratchpad[i*s->linesize+ 0]) = mean;

  1240.                 *(uint32_t*)(&s->me.scratchpad[i*s->linesize+ 4]) = mean;

  1241.                 *(uint32_t*)(&s->me.scratchpad[i*s->linesize+ 8]) = mean;

  1242.                 *(uint32_t*)(&s->me.scratchpad[i*s->linesize+12]) = mean;

  1243.             }

  1244. 

  1245.             intra_score= s->dsp.mb_cmp[0](s, s->me.scratchpad, pix, s->linesize, 16);

  1246.         }

  1247. #if 0 //FIXME

  1248.         /* get chroma score */

  1249.         if(s->avctx->mb_cmp&FF_CMP_CHROMA){

  1250.             for(i=1; i<3; i++){

  1251.                 uint8_t *dest_c;

  1252.                 int mean;

  1253.                 

  1254.                 if(s->out_format == FMT_H263){

  1255.                     mean= (s->dc_val[i][mb_x + (mb_y+1)*(s->mb_width+2)] + 4)>>3; //FIXME not exact but simple ;)

  1256.                 }else{

  1257.                     mean= (s->last_dc[i] + 4)>>3;

  1258.                 }

  1259.                 dest_c = s->new_picture.data[i] + (mb_y * 8  * (s->uvlinesize)) + mb_x * 8;

  1260.                 

  1261.                 mean*= 0x01010101;

  1262.                 for(i=0; i<8; i++){

  1263.                     *(uint32_t*)(&s->me.scratchpad[i*s->uvlinesize+ 0]) = mean;

  1264.                     *(uint32_t*)(&s->me.scratchpad[i*s->uvlinesize+ 4]) = mean;

  1265.                 }

  1266.                 

  1267.                 intra_score+= s->dsp.mb_cmp[1](s, s->me.scratchpad, dest_c, s->uvlinesize);

  1268.             }                

  1269.         }

  1270. #endif

  1271.         intra_score += s->me.mb_penalty_factor*16;

  1272.         

  1273.         if(intra_score < dmin){

  1274.             mb_type= CANDIDATE_MB_TYPE_INTRA;

  1275.             s->current_picture.mb_type[mb_y*s->mb_stride + mb_x]= CANDIDATE_MB_TYPE_INTRA; //FIXME cleanup

  1276.         }else

  1277.             s->current_picture.mb_type[mb_y*s->mb_stride + mb_x]= 0;

  1278.         

  1279.         if (vard <= 64 || vard < varc) { //FIXME

  1280.             s->scene_change_score+= ff_sqrt(vard) - ff_sqrt(varc);

  1281.         }else{

  1282.             s->scene_change_score+= s->qscale;

  1283.         }

  1284.     }

  1285. 

  1286.     s->mb_type[mb_y*s->mb_stride + mb_x]= mb_type;

  1287. }

  1288. 

  1289. int ff_pre_estimate_p_frame_motion(MpegEncContext * s,

  1290.                                     int mb_x, int mb_y)

  1291. {

  1292.     int mx, my, dmin;

  1293.     int pred_x=0, pred_y=0;

  1294.     int P[10][2];

  1295.     const int shift= 1+s->quarter_sample;

  1296.     uint8_t * const mv_penalty= s->me.mv_penalty[s->f_code] + MAX_MV;

  1297.     const int xy= mb_x + mb_y*s->mb_stride;

  1298.     const int stride= s->linesize;

  1299.     const int uvstride= s->uvlinesize;

  1300.     uint8_t *src_data[3]= {

  1301.         s->new_picture.data[0] + 16*(mb_x + stride*mb_y),

  1302.         s->new_picture.data[1] + 8*(mb_x + uvstride*mb_y),

  1303.         s->new_picture.data[2] + 8*(mb_x + uvstride*mb_y)

  1304.     };

  1305.     uint8_t *ref_data[3]= {

  1306.         s->last_picture.data[0] + 16*(mb_x + stride*mb_y),

  1307.         s->last_picture.data[1] + 8*(mb_x + uvstride*mb_y),

  1308.         s->last_picture.data[2] + 8*(mb_x + uvstride*mb_y)

  1309.     };

  1310.     

  1311.     assert(s->quarter_sample==0 || s->quarter_sample==1);

  1312. 

  1313.     s->me.pre_penalty_factor    = get_penalty_factor(s, s->avctx->me_pre_cmp);

  1314. 

  1315.     get_limits(s, 16*mb_x, 16*mb_y);

  1316.     s->me.skip=0;

  1317. 

  1318.     P_LEFT[0]       = s->p_mv_table[xy + 1][0];

  1319.     P_LEFT[1]       = s->p_mv_table[xy + 1][1];

  1320. 

  1321.     if(P_LEFT[0]       < (s->me.xmin<<shift)) P_LEFT[0]       = (s->me.xmin<<shift);

  1322. 

  1323.     /* special case for first line */

  1324.     if (s->first_slice_line) {

  1325.         pred_x= P_LEFT[0];

  1326.         pred_y= P_LEFT[1];

  1327.         P_TOP[0]= P_TOPRIGHT[0]= P_MEDIAN[0]=

  1328.         P_TOP[1]= P_TOPRIGHT[1]= P_MEDIAN[1]= 0; //FIXME 

  1329.     } else {

  1330.         P_TOP[0]      = s->p_mv_table[xy + s->mb_stride    ][0];

  1331.         P_TOP[1]      = s->p_mv_table[xy + s->mb_stride    ][1];

  1332.         P_TOPRIGHT[0] = s->p_mv_table[xy + s->mb_stride - 1][0];

  1333.         P_TOPRIGHT[1] = s->p_mv_table[xy + s->mb_stride - 1][1];

  1334.         if(P_TOP[1]      < (s->me.ymin<<shift)) P_TOP[1]     = (s->me.ymin<<shift);

  1335.         if(P_TOPRIGHT[0] > (s->me.xmax<<shift)) P_TOPRIGHT[0]= (s->me.xmax<<shift);

  1336.         if(P_TOPRIGHT[1] < (s->me.ymin<<shift)) P_TOPRIGHT[1]= (s->me.ymin<<shift);

  1337.     

  1338.         P_MEDIAN[0]= mid_pred(P_LEFT[0], P_TOP[0], P_TOPRIGHT[0]);

  1339.         P_MEDIAN[1]= mid_pred(P_LEFT[1], P_TOP[1], P_TOPRIGHT[1]);

  1340. 

  1341.         pred_x = P_MEDIAN[0];

  1342.         pred_y = P_MEDIAN[1];

  1343.     }

  1344.     dmin = s->me.pre_motion_search(s, &mx, &my, P, pred_x, pred_y, 

  1345.                                    src_data, ref_data, stride, uvstride, s->p_mv_table, (1<<16)>>shift, mv_penalty);

  1346. 

  1347.     s->p_mv_table[xy][0] = mx<<shift;

  1348.     s->p_mv_table[xy][1] = my<<shift;

  1349.     

  1350.     return dmin;

  1351. }

  1352. 

  1353. static int ff_estimate_motion_b(MpegEncContext * s,

  1354.                        int mb_x, int mb_y, int16_t (*mv_table)[2], uint8_t *src_data[3],

  1355.                        uint8_t *ref_data[3], int stride, int uvstride, int f_code)

  1356. {

  1357.     int mx, my, dmin;

  1358.     int pred_x=0, pred_y=0;

  1359.     int P[10][2];

  1360.     const int shift= 1+s->quarter_sample;

  1361.     const int mot_stride = s->mb_stride;

  1362.     const int mot_xy = mb_y*mot_stride + mb_x;

  1363.     uint8_t * const ref_picture= ref_data[0] - 16*s->mb_x - 16*s->mb_y*s->linesize; //FIXME ugly

  1364.     uint8_t * const mv_penalty= s->me.mv_penalty[f_code] + MAX_MV;

  1365.     int mv_scale;

  1366.         

  1367.     s->me.penalty_factor    = get_penalty_factor(s, s->avctx->me_cmp);

  1368.     s->me.sub_penalty_factor= get_penalty_factor(s, s->avctx->me_sub_cmp);

  1369.     s->me.mb_penalty_factor = get_penalty_factor(s, s->avctx->mb_cmp);

  1370. 

  1371.     get_limits(s, 16*mb_x, 16*mb_y);

  1372. 

  1373.     switch(s->me_method) {

  1374.     case ME_ZERO:

  1375.     default:

  1376. 	no_motion_search(s, &mx, &my);

  1377.         dmin = 0;

  1378.         mx-= mb_x*16;

  1379.         my-= mb_y*16;

  1380.         break;

  1381. #if 0

  1382.     case ME_FULL:

  1383. 	dmin = full_motion_search(s, &mx, &my, range, ref_picture);

  1384.         mx-= mb_x*16;

  1385.         my-= mb_y*16;

  1386.         break;

  1387.     case ME_LOG:

  1388. 	dmin = log_motion_search(s, &mx, &my, range / 2, ref_picture);

  1389.         mx-= mb_x*16;

  1390.         my-= mb_y*16;

  1391.         break;

  1392.     case ME_PHODS:

  1393. 	dmin = phods_motion_search(s, &mx, &my, range / 2, ref_picture);

  1394.         mx-= mb_x*16;

  1395.         my-= mb_y*16;

  1396.         break;

  1397. #endif

  1398.     case ME_X1:

  1399.     case ME_EPZS:

  1400.        {

  1401.             P_LEFT[0]        = mv_table[mot_xy - 1][0];

  1402.             P_LEFT[1]        = mv_table[mot_xy - 1][1];

  1403. 

  1404.             if(P_LEFT[0]       > (s->me.xmax<<shift)) P_LEFT[0]       = (s->me.xmax<<shift);

  1405. 

  1406.             /* special case for first line */

  1407.             if (!s->first_slice_line) {

  1408.                 P_TOP[0] = mv_table[mot_xy - mot_stride             ][0];

  1409.                 P_TOP[1] = mv_table[mot_xy - mot_stride             ][1];

  1410.                 P_TOPRIGHT[0] = mv_table[mot_xy - mot_stride + 1         ][0];

  1411.                 P_TOPRIGHT[1] = mv_table[mot_xy - mot_stride + 1         ][1];

  1412.                 if(P_TOP[1] > (s->me.ymax<<shift)) P_TOP[1]= (s->me.ymax<<shift);

  1413.                 if(P_TOPRIGHT[0] < (s->me.xmin<<shift)) P_TOPRIGHT[0]= (s->me.xmin<<shift);

  1414.                 if(P_TOPRIGHT[1] > (s->me.ymax<<shift)) P_TOPRIGHT[1]= (s->me.ymax<<shift);

  1415.         

  1416.                 P_MEDIAN[0]= mid_pred(P_LEFT[0], P_TOP[0], P_TOPRIGHT[0]);

  1417.                 P_MEDIAN[1]= mid_pred(P_LEFT[1], P_TOP[1], P_TOPRIGHT[1]);

  1418.             }

  1419.             pred_x= P_LEFT[0];

  1420.             pred_y= P_LEFT[1];

  1421.         }

  1422.         

  1423.         if(mv_table == s->b_forw_mv_table){

  1424.             mv_scale= (s->pb_time<<16) / (s->pp_time<<shift);

  1425.         }else{

  1426.             mv_scale= ((s->pb_time - s->pp_time)<<16) / (s->pp_time<<shift);

  1427.         }

  1428.         

  1429.         dmin = s->me.motion_search[0](s, &mx, &my, P, pred_x, pred_y, 

  1430.                                       src_data, ref_data, stride, uvstride, s->p_mv_table, mv_scale, mv_penalty);

  1431.  

  1432.         break;

  1433.     }

  1434.     

  1435.     dmin= s->me.sub_motion_search(s, &mx, &my, dmin,

  1436. 				   pred_x, pred_y, src_data, ref_data, stride, uvstride, 0, 16, mv_penalty);

  1437.                                    

  1438.     if(s->avctx->me_sub_cmp != s->avctx->mb_cmp && !s->me.skip)

  1439.         dmin= s->me.get_mb_score(s, mx, my, pred_x, pred_y, src_data, ref_data, stride, uvstride, mv_penalty);

  1440. 

  1441. //printf("%d %d %d %d//", s->mb_x, s->mb_y, mx, my);

  1442. //    s->mb_type[mb_y*s->mb_width + mb_x]= mb_type;

  1443.     mv_table[mot_xy][0]= mx;

  1444.     mv_table[mot_xy][1]= my;

  1445. 

  1446.     return dmin;

  1447. }

  1448. 

  1449. static inline int check_bidir_mv(MpegEncContext * s, uint8_t *src_data[3], uint8_t *ref_data[6],

  1450.                    int stride, int uvstride,

  1451.                    int motion_fx, int motion_fy,

  1452.                    int motion_bx, int motion_by,

  1453.                    int pred_fx, int pred_fy,

  1454.                    int pred_bx, int pred_by,

  1455.                    int size, int h)

  1456. {

  1457.     //FIXME optimize?

  1458.     //FIXME move into template?

  1459.     //FIXME better f_code prediction (max mv & distance)

  1460.     //FIXME pointers

  1461.     uint8_t * const mv_penalty= s->me.mv_penalty[s->f_code] + MAX_MV; // f_code of the prev frame

  1462.     uint8_t *dest_y = s->me.scratchpad;

  1463.     uint8_t *ptr;

  1464.     int dxy;

  1465.     int src_x, src_y;

  1466.     int fbmin;

  1467. 

  1468.     if(s->quarter_sample){

  1469.         dxy = ((motion_fy & 3) << 2) | (motion_fx & 3);

  1470.         src_x = motion_fx >> 2;

  1471.         src_y = motion_fy >> 2;

  1472. 

  1473.         ptr = ref_data[0] + (src_y * stride) + src_x;

  1474.         s->dsp.put_qpel_pixels_tab[0][dxy](dest_y    , ptr    , stride);

  1475. 

  1476.         dxy = ((motion_by & 3) << 2) | (motion_bx & 3);

  1477.         src_x = motion_bx >> 2;

  1478.         src_y = motion_by >> 2;

  1479.     

  1480.         ptr = ref_data[3] + (src_y * stride) + src_x;

  1481.         s->dsp.avg_qpel_pixels_tab[size][dxy](dest_y    , ptr    , stride);

  1482.     }else{

  1483.         dxy = ((motion_fy & 1) << 1) | (motion_fx & 1);

  1484.         src_x = motion_fx >> 1;

  1485.         src_y = motion_fy >> 1;

  1486. 

  1487.         ptr = ref_data[0] + (src_y * stride) + src_x;

  1488.         s->dsp.put_pixels_tab[size][dxy](dest_y    , ptr    , stride, h);

  1489. 

  1490.         dxy = ((motion_by & 1) << 1) | (motion_bx & 1);

  1491.         src_x = motion_bx >> 1;

  1492.         src_y = motion_by >> 1;

  1493.     

  1494.         ptr = ref_data[3] + (src_y * stride) + src_x;

  1495.         s->dsp.avg_pixels_tab[size][dxy](dest_y    , ptr    , stride, h);

  1496.     }

  1497. 

  1498.     fbmin = (mv_penalty[motion_fx-pred_fx] + mv_penalty[motion_fy-pred_fy])*s->me.mb_penalty_factor

  1499.            +(mv_penalty[motion_bx-pred_bx] + mv_penalty[motion_by-pred_by])*s->me.mb_penalty_factor

  1500.            + s->dsp.mb_cmp[size](s, src_data[0], dest_y, stride, h); //FIXME new_pic

  1501.            

  1502.     if(s->avctx->mb_cmp&FF_CMP_CHROMA){

  1503.     }

  1504.     //FIXME CHROMA !!!

  1505.            

  1506.     return fbmin;

  1507. }

  1508. 

  1509. /* refine the bidir vectors in hq mode and return the score in both lq & hq mode*/

  1510. static inline int bidir_refine(MpegEncContext * s, uint8_t *src_data[3], uint8_t *ref_data[6],

  1511.                                   int stride, int uvstride,

  1512.                                   int mb_x, int mb_y)

  1513. {

  1514.     const int mot_stride = s->mb_stride;

  1515.     const int xy = mb_y *mot_stride + mb_x;

  1516.     int fbmin;

  1517.     int pred_fx= s->b_bidir_forw_mv_table[xy-1][0];

  1518.     int pred_fy= s->b_bidir_forw_mv_table[xy-1][1];

  1519.     int pred_bx= s->b_bidir_back_mv_table[xy-1][0];

  1520.     int pred_by= s->b_bidir_back_mv_table[xy-1][1];

  1521.     int motion_fx= s->b_bidir_forw_mv_table[xy][0]= s->b_forw_mv_table[xy][0];

  1522.     int motion_fy= s->b_bidir_forw_mv_table[xy][1]= s->b_forw_mv_table[xy][1];

  1523.     int motion_bx= s->b_bidir_back_mv_table[xy][0]= s->b_back_mv_table[xy][0];

  1524.     int motion_by= s->b_bidir_back_mv_table[xy][1]= s->b_back_mv_table[xy][1];

  1525. 

  1526.     //FIXME do refinement and add flag

  1527.     

  1528.     fbmin= check_bidir_mv(s, src_data, ref_data, stride, uvstride,

  1529.                           motion_fx, motion_fy,

  1530.                           motion_bx, motion_by,

  1531.                           pred_fx, pred_fy,

  1532.                           pred_bx, pred_by,

  1533.                           0, 16);

  1534. 

  1535.    return fbmin;

  1536. }

  1537. 

  1538. static inline int direct_search(MpegEncContext * s, uint8_t *src_data[3], uint8_t *ref_data[6],

  1539.                                 int stride, int uvstride,

  1540.                                 int mb_x, int mb_y)

  1541. {

  1542.     int P[10][2];

  1543.     const int mot_stride = s->mb_stride;

  1544.     const int mot_xy = mb_y*mot_stride + mb_x;

  1545.     const int shift= 1+s->quarter_sample;

  1546.     int dmin, i;

  1547.     const int time_pp= s->pp_time;

  1548.     const int time_pb= s->pb_time;

  1549.     int mx, my, xmin, xmax, ymin, ymax;

  1550.     int16_t (*mv_table)[2]= s->b_direct_mv_table;

  1551.     uint8_t * const mv_penalty= s->me.mv_penalty[1] + MAX_MV;

  1552.     

  1553.     ymin= xmin=(-32)>>shift;

  1554.     ymax= xmax=   31>>shift;

  1555. 

  1556.     if(IS_8X8(s->next_picture.mb_type[mot_xy])){

  1557.         s->mv_type= MV_TYPE_8X8;

  1558.     }else{

  1559.         s->mv_type= MV_TYPE_16X16;

  1560.     }

  1561. 

  1562.     for(i=0; i<4; i++){

  1563.         int index= s->block_index[i];

  1564.         int min, max;

  1565.     

  1566.         s->me.co_located_mv[i][0]= s->next_picture.motion_val[0][index][0];

  1567.         s->me.co_located_mv[i][1]= s->next_picture.motion_val[0][index][1];

  1568.         s->me.direct_basis_mv[i][0]= s->me.co_located_mv[i][0]*time_pb/time_pp + ((i& 1)<<(shift+3));

  1569.         s->me.direct_basis_mv[i][1]= s->me.co_located_mv[i][1]*time_pb/time_pp + ((i>>1)<<(shift+3));

  1570. //        s->me.direct_basis_mv[1][i][0]= s->me.co_located_mv[i][0]*(time_pb - time_pp)/time_pp + ((i &1)<<(shift+3);

  1571. //        s->me.direct_basis_mv[1][i][1]= s->me.co_located_mv[i][1]*(time_pb - time_pp)/time_pp + ((i>>1)<<(shift+3);

  1572. 

  1573.         max= FFMAX(s->me.direct_basis_mv[i][0], s->me.direct_basis_mv[i][0] - s->me.co_located_mv[i][0])>>shift;

  1574.         min= FFMIN(s->me.direct_basis_mv[i][0], s->me.direct_basis_mv[i][0] - s->me.co_located_mv[i][0])>>shift;

  1575.         max+= 16*mb_x + 1; // +-1 is for the simpler rounding

  1576.         min+= 16*mb_x - 1;

  1577.         xmax= FFMIN(xmax, s->width - max);

  1578.         xmin= FFMAX(xmin, - 16     - min);

  1579. 

  1580.         max= FFMAX(s->me.direct_basis_mv[i][1], s->me.direct_basis_mv[i][1] - s->me.co_located_mv[i][1])>>shift;

  1581.         min= FFMIN(s->me.direct_basis_mv[i][1], s->me.direct_basis_mv[i][1] - s->me.co_located_mv[i][1])>>shift;

  1582.         max+= 16*mb_y + 1; // +-1 is for the simpler rounding

  1583.         min+= 16*mb_y - 1;

  1584.         ymax= FFMIN(ymax, s->height - max);

  1585.         ymin= FFMAX(ymin, - 16      - min);

  1586.         

  1587.         if(s->mv_type == MV_TYPE_16X16) break;

  1588.     }

  1589.     

  1590.     assert(xmax <= 15 && ymax <= 15 && xmin >= -16 && ymin >= -16);

  1591.     

  1592.     if(xmax < 0 || xmin >0 || ymax < 0 || ymin > 0){

  1593.         s->b_direct_mv_table[mot_xy][0]= 0;

  1594.         s->b_direct_mv_table[mot_xy][1]= 0;

  1595. 

  1596.         return 256*256*256*64;

  1597.     }

  1598.     

  1599.     s->me.xmin= xmin;

  1600.     s->me.ymin= ymin;

  1601.     s->me.xmax= xmax;

  1602.     s->me.ymax= ymax;

  1603. 

  1604.     P_LEFT[0]        = clip(mv_table[mot_xy - 1][0], xmin<<shift, xmax<<shift);

  1605.     P_LEFT[1]        = clip(mv_table[mot_xy - 1][1], ymin<<shift, ymax<<shift);

  1606. 

  1607.     /* special case for first line */

  1608.     if (!s->first_slice_line) { //FIXME maybe allow this over thread boundary as its cliped

  1609.         P_TOP[0]      = clip(mv_table[mot_xy - mot_stride             ][0], xmin<<shift, xmax<<shift);

  1610.         P_TOP[1]      = clip(mv_table[mot_xy - mot_stride             ][1], ymin<<shift, ymax<<shift);

  1611.         P_TOPRIGHT[0] = clip(mv_table[mot_xy - mot_stride + 1         ][0], xmin<<shift, xmax<<shift);

  1612.         P_TOPRIGHT[1] = clip(mv_table[mot_xy - mot_stride + 1         ][1], ymin<<shift, ymax<<shift);

  1613.     

  1614.         P_MEDIAN[0]= mid_pred(P_LEFT[0], P_TOP[0], P_TOPRIGHT[0]);

  1615.         P_MEDIAN[1]= mid_pred(P_LEFT[1], P_TOP[1], P_TOPRIGHT[1]);

  1616.     }

  1617.  

  1618.     //FIXME direct_search  ptr in context!!! (needed for chroma anyway or this will get messy)   

  1619.     if(s->flags&CODEC_FLAG_QPEL){

  1620.         dmin = simple_direct_qpel_epzs_motion_search(s, &mx, &my, P, 0, 0, 

  1621.                                                      src_data, ref_data, stride, uvstride, mv_table, 1<<14, mv_penalty);

  1622.         dmin = simple_direct_qpel_qpel_motion_search(s, &mx, &my, dmin,

  1623.                                                 0, 0, src_data, ref_data, stride, uvstride, 0, 16, mv_penalty);

  1624.         

  1625.         if(s->avctx->me_sub_cmp != s->avctx->mb_cmp && !s->me.skip)

  1626.             dmin= simple_direct_qpel_qpel_get_mb_score(s, mx, my, 0, 0, src_data, ref_data, stride, uvstride, mv_penalty);

  1627.     }else{

  1628.         dmin = simple_direct_hpel_epzs_motion_search(s, &mx, &my, P, 0, 0, 

  1629.                                                      src_data, ref_data, stride, uvstride, mv_table, 1<<15, mv_penalty);

  1630.         dmin = simple_direct_hpel_hpel_motion_search(s, &mx, &my, dmin,

  1631.                                                 0, 0, src_data, ref_data, stride, uvstride, 0, 16, mv_penalty);

  1632. 

  1633.         if(s->avctx->me_sub_cmp != s->avctx->mb_cmp && !s->me.skip)

  1634.             dmin= simple_direct_hpel_hpel_get_mb_score(s, mx, my, 0, 0, src_data, ref_data, stride, uvstride, mv_penalty);

  1635.     }

  1636.     

  1637.     get_limits(s, 16*mb_x, 16*mb_y); //restore s->me.?min/max, maybe not needed

  1638. 

  1639.     s->b_direct_mv_table[mot_xy][0]= mx;

  1640.     s->b_direct_mv_table[mot_xy][1]= my;

  1641.     return dmin;

  1642. }

  1643. 

  1644. void ff_estimate_b_frame_motion(MpegEncContext * s,

  1645.                              int mb_x, int mb_y)

  1646. {

  1647.     const int penalty_factor= s->me.mb_penalty_factor;

  1648.     int fmin, bmin, dmin, fbmin, bimin, fimin;

  1649.     int type=0;

  1650.     const int stride= s->linesize;

  1651.     const int uvstride= s->uvlinesize;

  1652.     uint8_t *src_data[3]= {

  1653.         s->new_picture.data[0] + 16*(s->mb_x + stride*s->mb_y),

  1654.         s->new_picture.data[1] + 8*(s->mb_x + uvstride*s->mb_y),

  1655.         s->new_picture.data[2] + 8*(s->mb_x + uvstride*s->mb_y)

  1656.     };

  1657.     uint8_t *ref_data[6]= {

  1658.         s->last_picture.data[0] + 16*(s->mb_x + stride*s->mb_y),

  1659.         s->last_picture.data[1] + 8*(s->mb_x + uvstride*s->mb_y),

  1660.         s->last_picture.data[2] + 8*(s->mb_x + uvstride*s->mb_y),

  1661.         s->next_picture.data[0] + 16*(s->mb_x + stride*s->mb_y),

  1662.         s->next_picture.data[1] + 8*(s->mb_x + uvstride*s->mb_y),

  1663.         s->next_picture.data[2] + 8*(s->mb_x + uvstride*s->mb_y)

  1664.     };

  1665.     

  1666.     s->me.skip=0;

  1667.     if (s->codec_id == CODEC_ID_MPEG4)

  1668.         dmin= direct_search(s, src_data, ref_data, stride, uvstride, mb_x, mb_y);

  1669.     else

  1670.         dmin= INT_MAX;

  1671. //FIXME penalty stuff for non mpeg4

  1672.     s->me.skip=0;

  1673.     fmin= ff_estimate_motion_b(s, mb_x, mb_y, s->b_forw_mv_table, src_data, 

  1674.                                ref_data, stride, uvstride, s->f_code) + 3*penalty_factor;

  1675.     

  1676.     s->me.skip=0;

  1677.     bmin= ff_estimate_motion_b(s, mb_x, mb_y, s->b_back_mv_table, src_data, 

  1678.                                ref_data+3, stride, uvstride, s->b_code) + 2*penalty_factor;

  1679. //printf(" %d %d ", s->b_forw_mv_table[xy][0], s->b_forw_mv_table[xy][1]);

  1680. 

  1681.     s->me.skip=0;

  1682.     fbmin= bidir_refine(s, src_data, ref_data, stride, uvstride, mb_x, mb_y) + penalty_factor;

  1683. //printf("%d %d %d %d\n", dmin, fmin, bmin, fbmin);

  1684.     

  1685.     if(s->flags & CODEC_FLAG_INTERLACED_ME){

  1686.         const int xy = mb_y*s->mb_stride + mb_x;

  1687. 

  1688. //FIXME mb type penalty

  1689.         s->me.skip=0;

  1690.         fimin= interlaced_search(s, src_data, ref_data  , 

  1691.                                  s->b_field_mv_table[0], s->b_field_select_table[0], s->f_code,

  1692.                                  s->b_forw_mv_table[xy][0], s->b_forw_mv_table[xy][1]);

  1693.         bimin= interlaced_search(s, src_data, ref_data+3, 

  1694.                                  s->b_field_mv_table[1], s->b_field_select_table[1], s->b_code,

  1695.                                  s->b_back_mv_table[xy][0], s->b_back_mv_table[xy][1]);

  1696.     }else

  1697.         fimin= bimin= INT_MAX;

  1698. 

  1699.     {

  1700.         int score= fmin;

  1701.         type = CANDIDATE_MB_TYPE_FORWARD;

  1702.         

  1703.         if (dmin <= score){

  1704.             score = dmin;

  1705.             type = CANDIDATE_MB_TYPE_DIRECT;

  1706.         }

  1707.         if(bmin<score){

  1708.             score=bmin;

  1709.             type= CANDIDATE_MB_TYPE_BACKWARD; 

  1710.         }

  1711.         if(fbmin<score){

  1712.             score=fbmin;

  1713.             type= CANDIDATE_MB_TYPE_BIDIR;

  1714.         }

  1715.         if(fimin<score){

  1716.             score=fimin;

  1717.             type= CANDIDATE_MB_TYPE_FORWARD_I;

  1718.         }

  1719.         if(bimin<score){

  1720.             score=bimin;

  1721.             type= CANDIDATE_MB_TYPE_BACKWARD_I;

  1722.         }

  1723.         

  1724.         score= ((unsigned)(score*score + 128*256))>>16;

  1725.         s->mc_mb_var_sum_temp += score;

  1726.         s->current_picture.mc_mb_var[mb_y*s->mb_stride + mb_x] = score; //FIXME use SSE

  1727.     }

  1728. 

  1729.     if(s->avctx->mb_decision > FF_MB_DECISION_SIMPLE){

  1730.         type= CANDIDATE_MB_TYPE_FORWARD | CANDIDATE_MB_TYPE_BACKWARD | CANDIDATE_MB_TYPE_BIDIR | CANDIDATE_MB_TYPE_DIRECT;

  1731.         if(fimin < INT_MAX)

  1732.             type |= CANDIDATE_MB_TYPE_FORWARD_I;

  1733.         if(bimin < INT_MAX)

  1734.             type |= CANDIDATE_MB_TYPE_BACKWARD_I;

  1735.         if(fimin < INT_MAX && bimin < INT_MAX){

  1736.             type |= CANDIDATE_MB_TYPE_BIDIR_I;

  1737.         }

  1738.          //FIXME something smarter

  1739.         if(dmin>256*256*16) type&= ~CANDIDATE_MB_TYPE_DIRECT; //dont try direct mode if its invalid for this MB

  1740. #if 0        

  1741.         if(s->out_format == FMT_MPEG1)

  1742.             type |= CANDIDATE_MB_TYPE_INTRA;

  1743. #endif

  1744.     }

  1745. 

  1746.     s->mb_type[mb_y*s->mb_stride + mb_x]= type;

  1747. }

  1748. 

  1749. /* find best f_code for ME which do unlimited searches */

  1750. int ff_get_best_fcode(MpegEncContext * s, int16_t (*mv_table)[2], int type)

  1751. {

  1752.     if(s->me_method>=ME_EPZS){

  1753.         int score[8];

  1754.         int i, y;

  1755.         uint8_t * fcode_tab= s->fcode_tab;

  1756.         int best_fcode=-1;

  1757.         int best_score=-10000000;

  1758. 

  1759.         for(i=0; i<8; i++) score[i]= s->mb_num*(8-i);

  1760. 

  1761.         for(y=0; y<s->mb_height; y++){

  1762.             int x;

  1763.             int xy= y*s->mb_stride;

  1764.             for(x=0; x<s->mb_width; x++){

  1765.                 if(s->mb_type[xy] & type){

  1766.                     int fcode= FFMAX(fcode_tab[mv_table[xy][0] + MAX_MV],

  1767.                                      fcode_tab[mv_table[xy][1] + MAX_MV]);

  1768.                     int j;

  1769.                     

  1770.                     for(j=0; j<fcode && j<8; j++){

  1771.                         if(s->pict_type==B_TYPE || s->current_picture.mc_mb_var[xy] < s->current_picture.mb_var[xy])

  1772.                             score[j]-= 170;

  1773.                     }

  1774.                 }

  1775.                 xy++;

  1776.             }

  1777.         }

  1778.         

  1779.         for(i=1; i<8; i++){

  1780.             if(score[i] > best_score){

  1781.                 best_score= score[i];

  1782.                 best_fcode= i;

  1783.             }

  1784. //            printf("%d %d\n", i, score[i]);

  1785.         }

  1786. 

  1787. //    printf("fcode: %d type: %d\n", i, s->pict_type);

  1788.         return best_fcode;

  1789. /*        for(i=0; i<=MAX_FCODE; i++){

  1790.             printf("%d ", mv_num[i]);

  1791.         }

  1792.         printf("\n");*/

  1793.     }else{

  1794.         return 1;

  1795.     }

  1796. }

  1797. 

  1798. void ff_fix_long_p_mvs(MpegEncContext * s)

  1799. {

  1800.     const int f_code= s->f_code;

  1801.     int y, range;

  1802.     assert(s->pict_type==P_TYPE);

  1803. 

  1804.     range = (((s->out_format == FMT_MPEG1) ? 8 : 16) << f_code);

  1805.     

  1806.     if(s->msmpeg4_version) range= 16;

  1807.     

  1808.     if(s->avctx->me_range && range > s->avctx->me_range) range= s->avctx->me_range;

  1809.     

  1810. //printf("%d no:%d %d//\n", clip, noclip, f_code);

  1811.     if(s->flags&CODEC_FLAG_4MV){

  1812.         const int wrap= 2+ s->mb_width*2;

  1813. 

  1814.         /* clip / convert to intra 8x8 type MVs */

  1815.         for(y=0; y<s->mb_height; y++){

  1816.             int xy= (y*2 + 1)*wrap + 1;

  1817.             int i= y*s->mb_stride;

  1818.             int x;

  1819. 

  1820.             for(x=0; x<s->mb_width; x++){

  1821.                 if(s->mb_type[i]&CANDIDATE_MB_TYPE_INTER4V){

  1822.                     int block;

  1823.                     for(block=0; block<4; block++){

  1824.                         int off= (block& 1) + (block>>1)*wrap;

  1825.                         int mx= s->current_picture.motion_val[0][ xy + off ][0];

  1826.                         int my= s->current_picture.motion_val[0][ xy + off ][1];

  1827. 

  1828.                         if(   mx >=range || mx <-range

  1829.                            || my >=range || my <-range){

  1830.                             s->mb_type[i] &= ~CANDIDATE_MB_TYPE_INTER4V;

  1831.                             s->mb_type[i] |= CANDIDATE_MB_TYPE_INTRA;

  1832.                             s->current_picture.mb_type[i]= CANDIDATE_MB_TYPE_INTRA;

  1833.                         }

  1834.                     }

  1835.                 }

  1836.                 xy+=2;

  1837.                 i++;

  1838.             }

  1839.         }

  1840.     }

  1841. }

  1842. 

  1843. /**

  1844.  *

  1845.  * @param truncate 1 for truncation, 0 for using intra

  1846.  */

  1847. void ff_fix_long_mvs(MpegEncContext * s, uint8_t *field_select_table, int field_select, 

  1848.                      int16_t (*mv_table)[2], int f_code, int type, int truncate)

  1849. {

  1850.     int y, h_range, v_range;

  1851. 

  1852.     // RAL: 8 in MPEG-1, 16 in MPEG-4

  1853.     int range = (((s->out_format == FMT_MPEG1) ? 8 : 16) << f_code);

  1854. 

  1855.     if(s->msmpeg4_version) range= 16;

  1856.     if(s->avctx->me_range && range > s->avctx->me_range) range= s->avctx->me_range;

  1857. 

  1858.     h_range= range;

  1859.     v_range= field_select_table ? range>>1 : range;

  1860. 

  1861.     /* clip / convert to intra 16x16 type MVs */

  1862.     for(y=0; y<s->mb_height; y++){

  1863.         int x;

  1864.         int xy= y*s->mb_stride;

  1865.         for(x=0; x<s->mb_width; x++){

  1866.             if (s->mb_type[xy] & type){    // RAL: "type" test added...

  1867.                 if(field_select_table==NULL || field_select_table[xy] == field_select){

  1868.                     if(   mv_table[xy][0] >=h_range || mv_table[xy][0] <-h_range

  1869.                        || mv_table[xy][1] >=v_range || mv_table[xy][1] <-v_range){

  1870. 

  1871.                         if(truncate){

  1872.                             if     (mv_table[xy][0] > h_range-1) mv_table[xy][0]=  h_range-1;

  1873.                             else if(mv_table[xy][0] < -h_range ) mv_table[xy][0]= -h_range;

  1874.                             if     (mv_table[xy][1] > v_range-1) mv_table[xy][1]=  v_range-1;

  1875.                             else if(mv_table[xy][1] < -v_range ) mv_table[xy][1]= -v_range;

  1876.                         }else{

  1877.                             s->mb_type[xy] &= ~type;

  1878.                             s->mb_type[xy] |= CANDIDATE_MB_TYPE_INTRA;

  1879.                             mv_table[xy][0]=

  1880.                             mv_table[xy][1]= 0;

  1881.                         }

  1882.                     }

  1883.                 }

  1884.             }

  1885.             xy++;

  1886.         }

  1887.     }

  1888. }