frame.c
上传用户:hjq518
上传日期:2021-12-09
资源大小:5084k
文件大小:26k
源码类别:

Audio

开发平台:

Visual C++

  1. /*****************************************************************************
  2.  * frame.c: h264 encoder library
  3.  *****************************************************************************
  4.  * Copyright (C) 2003 Laurent Aimar
  5.  * $Id: frame.c,v 1.1 2004/06/03 19:27:06 fenrir Exp $
  6.  *
  7.  * Authors: Laurent Aimar <fenrir@via.ecp.fr>
  8.  *
  9.  * This program is free software; you can redistribute it and/or modify
  10.  * it under the terms of the GNU General Public License as published by
  11.  * the Free Software Foundation; either version 2 of the License, or
  12.  * (at your option) any later version.
  13.  *
  14.  * This program is distributed in the hope that it will be useful,
  15.  * but WITHOUT ANY WARRANTY; without even the implied warranty of
  16.  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  17.  * GNU General Public License for more details.
  18.  *
  19.  * You should have received a copy of the GNU General Public License
  20.  * along with this program; if not, write to the Free Software
  21.  * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA  02111, USA.
  22.  *****************************************************************************/
  23. #include <stdio.h>
  24. #include <string.h>
  25. #include "common.h"
  26. x264_frame_t *x264_frame_new( x264_t *h )
  27. {
  28.     x264_frame_t   *frame = x264_malloc( sizeof( x264_frame_t ) );
  29.     int i;
  30.     int i_mb_count = h->mb.i_mb_count;
  31.     int i_stride;
  32.     int i_lines;
  33.     memset( frame, 0, sizeof(x264_frame_t) );
  34.     /* allocate frame data (+64 for extra data for me) */
  35.     i_stride = ( ( h->param.i_width  + 15 )&0xfffff0 )+ 64;
  36.     i_lines  = ( ( h->param.i_height + 15 )&0xfffff0 );
  37.     frame->i_plane = 3;
  38.     for( i = 0; i < 3; i++ )
  39.     {
  40.         int i_divh = 1;
  41.         int i_divw = 1;
  42.         if( i > 0 )
  43.         {
  44.             if( h->param.i_csp == X264_CSP_I420 )
  45.                 i_divh = i_divw = 2;
  46.             else if( h->param.i_csp == X264_CSP_I422 )
  47.                 i_divw = 2;
  48.         }
  49.         frame->i_stride[i] = i_stride / i_divw;
  50.         frame->i_lines[i] = i_lines / i_divh;
  51.         frame->buffer[i] = x264_malloc( frame->i_stride[i] *
  52.                                         ( frame->i_lines[i] + 64 / i_divh ) );
  53.         frame->plane[i] = ((uint8_t*)frame->buffer[i]) +
  54.                           frame->i_stride[i] * 32 / i_divh + 32 / i_divw;
  55.     }
  56.     frame->i_stride[3] = 0;
  57.     frame->i_lines[3] = 0;
  58.     frame->buffer[3] = NULL;
  59.     frame->plane[3] = NULL;
  60.     frame->filtered[0] = frame->plane[0];
  61.     for( i = 0; i < 3; i++ )
  62.     {
  63.         frame->buffer[4+i] = x264_malloc( frame->i_stride[0] *
  64.                                         ( frame->i_lines[0] + 64 ) );
  65.         frame->filtered[i+1] = ((uint8_t*)frame->buffer[4+i]) +
  66.                                 frame->i_stride[0] * 32 + 32;
  67.     }
  68.     if( h->frames.b_have_lowres )
  69.     {
  70.         frame->i_stride_lowres = frame->i_stride[0]/2 + 32;
  71.         frame->i_lines_lowres = frame->i_lines[0]/2;
  72.         for( i = 0; i < 4; i++ )
  73.         {
  74.             frame->buffer[7+i] = x264_malloc( frame->i_stride_lowres *
  75.                                             ( frame->i_lines[0]/2 + 64 ) );
  76.             frame->lowres[i] = ((uint8_t*)frame->buffer[7+i]) +
  77.                                 frame->i_stride_lowres * 32 + 32;
  78.         }
  79.     }
  80.     if( h->param.analyse.i_me_method == X264_ME_ESA )
  81.     {
  82.         frame->buffer[11] = x264_malloc( frame->i_stride[0] * (frame->i_lines[0] + 64) * sizeof(uint16_t) );
  83.         frame->integral = (uint16_t*)frame->buffer[11] + frame->i_stride[0] * 32 + 32;
  84.     }
  85.     frame->i_poc = -1;
  86.     frame->i_type = X264_TYPE_AUTO;
  87.     frame->i_qpplus1 = 0;
  88.     frame->i_pts = -1;
  89.     frame->i_frame = -1;
  90.     frame->i_frame_num = -1;
  91.     frame->mb_type= x264_malloc( i_mb_count * sizeof( int8_t) );
  92.     frame->mv[0]  = x264_malloc( 2*16 * i_mb_count * sizeof( int16_t ) );
  93.     frame->ref[0] = x264_malloc( 4 * i_mb_count * sizeof( int8_t ) );
  94.     if( h->param.i_bframe )
  95.     {
  96.         frame->mv[1]  = x264_malloc( 2*16 * i_mb_count * sizeof( int16_t ) );
  97.         frame->ref[1] = x264_malloc( 4 * i_mb_count * sizeof( int8_t ) );
  98.     }
  99.     else
  100.     {
  101.         frame->mv[1]  = NULL;
  102.         frame->ref[1] = NULL;
  103.     }
  104.     return frame;
  105. }
  106. void x264_frame_delete( x264_frame_t *frame )
  107. {
  108.     int i;
  109.     for( i = 0; i < frame->i_plane; i++ )
  110.     {
  111.         x264_free( frame->buffer[i] );
  112.     }
  113.     for( i = 4; i < 12; i++ ) /* filtered planes */
  114.     {
  115.         x264_free( frame->buffer[i] );
  116.     }
  117.     x264_free( frame->mb_type );
  118.     x264_free( frame->mv[0] );
  119.     x264_free( frame->mv[1] );
  120.     x264_free( frame->ref[0] );
  121.     x264_free( frame->ref[1] );
  122.     x264_free( frame );
  123. }
  124. void x264_frame_copy_picture( x264_t *h, x264_frame_t *dst, x264_picture_t *src )
  125. {
  126.     dst->i_type     = src->i_type;
  127.     dst->i_qpplus1  = src->i_qpplus1;
  128.     dst->i_pts      = src->i_pts;
  129.     switch( src->img.i_csp & X264_CSP_MASK )
  130.     {
  131.         case X264_CSP_I420:
  132.             h->csp.i420( dst, &src->img, h->param.i_width, h->param.i_height );
  133.             break;
  134.         case X264_CSP_YV12:
  135.             h->csp.yv12( dst, &src->img, h->param.i_width, h->param.i_height );
  136.             break;
  137.         case X264_CSP_I422:
  138.             h->csp.i422( dst, &src->img, h->param.i_width, h->param.i_height );
  139.             break;
  140.         case X264_CSP_I444:
  141.             h->csp.i444( dst, &src->img, h->param.i_width, h->param.i_height );
  142.             break;
  143.         case X264_CSP_YUYV:
  144.             h->csp.yuyv( dst, &src->img, h->param.i_width, h->param.i_height );
  145.             break;
  146.         case X264_CSP_RGB:
  147.             h->csp.rgb( dst, &src->img, h->param.i_width, h->param.i_height );
  148.             break;
  149.         case X264_CSP_BGR:
  150.             h->csp.bgr( dst, &src->img, h->param.i_width, h->param.i_height );
  151.             break;
  152.         case X264_CSP_BGRA:
  153.             h->csp.bgra( dst, &src->img, h->param.i_width, h->param.i_height );
  154.             break;
  155.         default:
  156.             x264_log( h, X264_LOG_ERROR, "Arg invalid CSPn" );
  157.             break;
  158.     }
  159. }
  160. static void plane_expand_border( uint8_t *pix, int i_stride, int i_height, int i_pad )
  161. {
  162. #define PPIXEL(x, y) ( pix + (x) + (y)*i_stride )
  163.     const int i_width = i_stride - 2*i_pad;
  164.     int y;
  165.     for( y = 0; y < i_height; y++ )
  166.     {
  167.         /* left band */
  168.         memset( PPIXEL(-i_pad, y), PPIXEL(0, y)[0], i_pad );
  169.         /* right band */
  170.         memset( PPIXEL(i_width, y), PPIXEL(i_width-1, y)[0], i_pad );
  171.     }
  172.     /* upper band */
  173.     for( y = 0; y < i_pad; y++ )
  174.         memcpy( PPIXEL(-i_pad, -y-1), PPIXEL(-i_pad, 0), i_stride );
  175.     /* lower band */
  176.     for( y = 0; y < i_pad; y++ )
  177.         memcpy( PPIXEL(-i_pad, i_height+y), PPIXEL(-i_pad, i_height-1), i_stride );
  178. #undef PPIXEL
  179. }
  180. void x264_frame_expand_border( x264_frame_t *frame )
  181. {
  182.     int i;
  183.     for( i = 0; i < frame->i_plane; i++ )
  184.     {
  185.         int i_pad = i ? 16 : 32;
  186.         plane_expand_border( frame->plane[i], frame->i_stride[i], frame->i_lines[i], i_pad );
  187.     }
  188. }
  189. void x264_frame_expand_border_filtered( x264_frame_t *frame )
  190. {
  191.     /* during filtering, 8 extra pixels were filtered on each edge. 
  192.        we want to expand border from the last filtered pixel */
  193.     int i;
  194.     for( i = 1; i < 4; i++ )
  195.         plane_expand_border( frame->filtered[i] - 8*frame->i_stride[0] - 8, frame->i_stride[0], frame->i_lines[0]+2*8, 24 );
  196. }
  197. void x264_frame_expand_border_lowres( x264_frame_t *frame )
  198. {
  199.     int i;
  200.     for( i = 0; i < 4; i++ )
  201.         plane_expand_border( frame->lowres[i], frame->i_stride_lowres, frame->i_lines_lowres, 32 );
  202. }
  203. void x264_frame_expand_border_mod16( x264_t *h, x264_frame_t *frame )
  204. {
  205.     int i, y;
  206.     for( i = 0; i < frame->i_plane; i++ )
  207.     {
  208.         int i_subsample = i ? 1 : 0;
  209.         int i_width = h->param.i_width >> i_subsample;
  210.         int i_height = h->param.i_height >> i_subsample;
  211.         int i_padx = ( h->sps->i_mb_width * 16 - h->param.i_width ) >> i_subsample;
  212.         int i_pady = ( h->sps->i_mb_height * 16 - h->param.i_height ) >> i_subsample;
  213.         if( i_padx )
  214.         {
  215.             for( y = 0; y < i_height; y++ )
  216.                 memset( &frame->plane[i][y*frame->i_stride[i] + i_width],
  217.                          frame->plane[i][y*frame->i_stride[i] + i_width - 1],
  218.                          i_padx );
  219.         }
  220.         if( i_pady )
  221.         {
  222.             for( y = i_height; y < i_height + i_pady; y++ );
  223.                 memcpy( &frame->plane[i][y*frame->i_stride[i]],
  224.                         &frame->plane[i][(i_height-1)*frame->i_stride[i]],
  225.                         i_width + i_padx );
  226.         }
  227.     }
  228. }
  229. /* Deblocking filter */
  230. static const int i_alpha_table[52] =
  231. {
  232.      0,  0,  0,  0,  0,  0,  0,  0,  0,  0,
  233.      0,  0,  0,  0,  0,  0,  4,  4,  5,  6,
  234.      7,  8,  9, 10, 12, 13, 15, 17, 20, 22,
  235.     25, 28, 32, 36, 40, 45, 50, 56, 63, 71,
  236.     80, 90,101,113,127,144,162,182,203,226,
  237.     255, 255
  238. };
  239. static const int i_beta_table[52] =
  240. {
  241.      0,  0,  0,  0,  0,  0,  0,  0,  0,  0,
  242.      0,  0,  0,  0,  0,  0,  2,  2,  2,  3,
  243.      3,  3,  3,  4,  4,  4,  6,  6,  7,  7,
  244.      8,  8,  9,  9, 10, 10, 11, 11, 12, 12,
  245.     13, 13, 14, 14, 15, 15, 16, 16, 17, 17,
  246.     18, 18
  247. };
  248. static const int i_tc0_table[52][3] =
  249. {
  250.     { 0, 0, 0 }, { 0, 0, 0 }, { 0, 0, 0 }, { 0, 0, 0 }, { 0, 0, 0 }, { 0, 0, 0 },
  251.     { 0, 0, 0 }, { 0, 0, 0 }, { 0, 0, 0 }, { 0, 0, 0 }, { 0, 0, 0 }, { 0, 0, 0 },
  252.     { 0, 0, 0 }, { 0, 0, 0 }, { 0, 0, 0 }, { 0, 0, 0 }, { 0, 0, 0 }, { 0, 0, 1 },
  253.     { 0, 0, 1 }, { 0, 0, 1 }, { 0, 0, 1 }, { 0, 1, 1 }, { 0, 1, 1 }, { 1, 1, 1 },
  254.     { 1, 1, 1 }, { 1, 1, 1 }, { 1, 1, 1 }, { 1, 1, 2 }, { 1, 1, 2 }, { 1, 1, 2 },
  255.     { 1, 1, 2 }, { 1, 2, 3 }, { 1, 2, 3 }, { 2, 2, 3 }, { 2, 2, 4 }, { 2, 3, 4 },
  256.     { 2, 3, 4 }, { 3, 3, 5 }, { 3, 4, 6 }, { 3, 4, 6 }, { 4, 5, 7 }, { 4, 5, 8 },
  257.     { 4, 6, 9 }, { 5, 7,10 }, { 6, 8,11 }, { 6, 8,13 }, { 7,10,14 }, { 8,11,16 },
  258.     { 9,12,18 }, {10,13,20 }, {11,15,23 }, {13,17,25 }
  259. };
  260. /* From ffmpeg */
  261. static inline int clip_uint8( int a )
  262. {
  263.     if (a&(~255))
  264.         return (-a)>>31;
  265.     else
  266.         return a;
  267. }
  268. static inline void deblock_luma_c( uint8_t *pix, int xstride, int ystride, int alpha, int beta, int8_t *tc0 )
  269. {
  270.     int i, d;
  271.     for( i = 0; i < 4; i++ ) {
  272.         if( tc0[i] < 0 ) {
  273.             pix += 4*ystride;
  274.             continue;
  275.         }
  276.         for( d = 0; d < 4; d++ ) {
  277.             const int p2 = pix[-3*xstride];
  278.             const int p1 = pix[-2*xstride];
  279.             const int p0 = pix[-1*xstride];
  280.             const int q0 = pix[ 0*xstride];
  281.             const int q1 = pix[ 1*xstride];
  282.             const int q2 = pix[ 2*xstride];
  283.    
  284.             if( abs( p0 - q0 ) < alpha &&
  285.                 abs( p1 - p0 ) < beta &&
  286.                 abs( q1 - q0 ) < beta ) {
  287.    
  288.                 int tc = tc0[i];
  289.                 int delta;
  290.    
  291.                 if( abs( p2 - p0 ) < beta ) {
  292.                     pix[-2*xstride] = p1 + x264_clip3( (( p2 + ((p0 + q0 + 1) >> 1)) >> 1) - p1, -tc0[i], tc0[i] );
  293.                     tc++; 
  294.                 }
  295.                 if( abs( q2 - q0 ) < beta ) {
  296.                     pix[ 1*xstride] = q1 + x264_clip3( (( q2 + ((p0 + q0 + 1) >> 1)) >> 1) - q1, -tc0[i], tc0[i] );
  297.                     tc++;
  298.                 }
  299.     
  300.                 delta = x264_clip3( (((q0 - p0 ) << 2) + (p1 - q1) + 4) >> 3, -tc, tc );
  301.                 pix[-1*xstride] = clip_uint8( p0 + delta );    /* p0' */
  302.                 pix[ 0*xstride] = clip_uint8( q0 - delta );    /* q0' */
  303.             }
  304.             pix += ystride;
  305.         }
  306.     }
  307. }
  308. static void deblock_v_luma_c( uint8_t *pix, int stride, int alpha, int beta, int8_t *tc0 )
  309. {
  310.     deblock_luma_c( pix, stride, 1, alpha, beta, tc0 ); 
  311. }
  312. static void deblock_h_luma_c( uint8_t *pix, int stride, int alpha, int beta, int8_t *tc0 )
  313. {
  314.     deblock_luma_c( pix, 1, stride, alpha, beta, tc0 );
  315. }
  316. static inline void deblock_chroma_c( uint8_t *pix, int xstride, int ystride, int alpha, int beta, int8_t *tc0 )
  317. {
  318.     int i, d;
  319.     for( i = 0; i < 4; i++ ) {
  320.         const int tc = tc0[i];
  321.         if( tc <= 0 ) {
  322.             pix += 2*ystride;
  323.             continue;
  324.         }
  325.         for( d = 0; d < 2; d++ ) {
  326.             const int p1 = pix[-2*xstride];
  327.             const int p0 = pix[-1*xstride];
  328.             const int q0 = pix[ 0*xstride];
  329.             const int q1 = pix[ 1*xstride];
  330.             if( abs( p0 - q0 ) < alpha &&
  331.                 abs( p1 - p0 ) < beta &&
  332.                 abs( q1 - q0 ) < beta ) {
  333.                 int delta = x264_clip3( (((q0 - p0 ) << 2) + (p1 - q1) + 4) >> 3, -tc, tc );
  334.                 pix[-1*xstride] = clip_uint8( p0 + delta );    /* p0' */
  335.                 pix[ 0*xstride] = clip_uint8( q0 - delta );    /* q0' */
  336.             }
  337.             pix += ystride;
  338.         }
  339.     }
  340. }
  341. static void deblock_v_chroma_c( uint8_t *pix, int stride, int alpha, int beta, int8_t *tc0 )
  342. {   
  343.     deblock_chroma_c( pix, stride, 1, alpha, beta, tc0 );
  344. }
  345. static void deblock_h_chroma_c( uint8_t *pix, int stride, int alpha, int beta, int8_t *tc0 )
  346. {   
  347.     deblock_chroma_c( pix, 1, stride, alpha, beta, tc0 );
  348. }
  349. static inline void deblock_luma_intra_c( uint8_t *pix, int xstride, int ystride, int alpha, int beta )
  350. {
  351.     int d;
  352.     for( d = 0; d < 16; d++ ) {
  353.         const int p2 = pix[-3*xstride];
  354.         const int p1 = pix[-2*xstride];
  355.         const int p0 = pix[-1*xstride];
  356.         const int q0 = pix[ 0*xstride];
  357.         const int q1 = pix[ 1*xstride];
  358.         const int q2 = pix[ 2*xstride];
  359.         if( abs( p0 - q0 ) < alpha &&
  360.             abs( p1 - p0 ) < beta &&
  361.             abs( q1 - q0 ) < beta ) {
  362.             if(abs( p0 - q0 ) < ((alpha >> 2) + 2) ){
  363.                 if( abs( p2 - p0 ) < beta)
  364.                 {
  365.                     const int p3 = pix[-4*xstride];
  366.                     /* p0', p1', p2' */
  367.                     pix[-1*xstride] = ( p2 + 2*p1 + 2*p0 + 2*q0 + q1 + 4 ) >> 3;
  368.                     pix[-2*xstride] = ( p2 + p1 + p0 + q0 + 2 ) >> 2;
  369.                     pix[-3*xstride] = ( 2*p3 + 3*p2 + p1 + p0 + q0 + 4 ) >> 3;
  370.                 } else {
  371.                     /* p0' */
  372.                     pix[-1*xstride] = ( 2*p1 + p0 + q1 + 2 ) >> 2;
  373.                 }
  374.                 if( abs( q2 - q0 ) < beta)
  375.                 {
  376.                     const int q3 = pix[3*xstride];
  377.                     /* q0', q1', q2' */
  378.                     pix[0*xstride] = ( p1 + 2*p0 + 2*q0 + 2*q1 + q2 + 4 ) >> 3;
  379.                     pix[1*xstride] = ( p0 + q0 + q1 + q2 + 2 ) >> 2;
  380.                     pix[2*xstride] = ( 2*q3 + 3*q2 + q1 + q0 + p0 + 4 ) >> 3;
  381.                 } else {
  382.                     /* q0' */
  383.                     pix[0*xstride] = ( 2*q1 + q0 + p1 + 2 ) >> 2;
  384.                 }
  385.             }else{
  386.                 /* p0', q0' */
  387.                 pix[-1*xstride] = ( 2*p1 + p0 + q1 + 2 ) >> 2;
  388.                 pix[ 0*xstride] = ( 2*q1 + q0 + p1 + 2 ) >> 2;
  389.             }
  390.         }
  391.         pix += ystride;
  392.     }
  393. }
  394. static void deblock_v_luma_intra_c( uint8_t *pix, int stride, int alpha, int beta )
  395. {   
  396.     deblock_luma_intra_c( pix, stride, 1, alpha, beta );
  397. }
  398. static void deblock_h_luma_intra_c( uint8_t *pix, int stride, int alpha, int beta )
  399. {   
  400.     deblock_luma_intra_c( pix, 1, stride, alpha, beta );
  401. }
  402. static inline void deblock_chroma_intra_c( uint8_t *pix, int xstride, int ystride, int alpha, int beta )
  403. {   
  404.     int d; 
  405.     for( d = 0; d < 8; d++ ) {
  406.         const int p1 = pix[-2*xstride];
  407.         const int p0 = pix[-1*xstride];
  408.         const int q0 = pix[ 0*xstride];
  409.         const int q1 = pix[ 1*xstride];
  410.         if( abs( p0 - q0 ) < alpha &&
  411.             abs( p1 - p0 ) < beta &&
  412.             abs( q1 - q0 ) < beta ) {
  413.             pix[-1*xstride] = (2*p1 + p0 + q1 + 2) >> 2;   /* p0' */
  414.             pix[ 0*xstride] = (2*q1 + q0 + p1 + 2) >> 2;   /* q0' */
  415.         }
  416.         pix += ystride;
  417.     }
  418. }
  419. static void deblock_v_chroma_intra_c( uint8_t *pix, int stride, int alpha, int beta )
  420. {   
  421.     deblock_chroma_intra_c( pix, stride, 1, alpha, beta );
  422. }
  423. static void deblock_h_chroma_intra_c( uint8_t *pix, int stride, int alpha, int beta )
  424. {   
  425.     deblock_chroma_intra_c( pix, 1, stride, alpha, beta );
  426. }
  427. static inline void deblock_edge( x264_t *h, uint8_t *pix, int i_stride, int bS[4], int i_qp, int b_chroma,
  428.                                  x264_deblock_inter_t pf_inter, x264_deblock_intra_t pf_intra )
  429. {
  430.     int i;
  431.     const int index_a = x264_clip3( i_qp + h->sh.i_alpha_c0_offset, 0, 51 );
  432.     const int alpha = i_alpha_table[index_a];
  433.     const int beta  = i_beta_table[x264_clip3( i_qp + h->sh.i_beta_offset, 0, 51 )];
  434.     if( bS[0] < 4 ) {
  435.         int8_t tc[4]; 
  436.         for(i=0; i<4; i++)
  437.             tc[i] = (bS[i] ? i_tc0_table[index_a][bS[i] - 1] : -1) + b_chroma;
  438.         pf_inter( pix, i_stride, alpha, beta, tc );
  439.     } else {
  440.         pf_intra( pix, i_stride, alpha, beta );
  441.     }
  442. }
  443. void x264_frame_deblocking_filter( x264_t *h, int i_slice_type )
  444. {
  445.     const int s8x8 = 2 * h->mb.i_mb_stride;
  446.     const int s4x4 = 4 * h->mb.i_mb_stride;
  447.     int mb_y, mb_x;
  448.     for( mb_y = 0, mb_x = 0; mb_y < h->sps->i_mb_height; )
  449.     {
  450.         const int mb_xy  = mb_y * h->mb.i_mb_stride + mb_x;
  451.         const int mb_8x8 = 2 * s8x8 * mb_y + 2 * mb_x;
  452.         const int mb_4x4 = 4 * s4x4 * mb_y + 4 * mb_x;
  453.         int i_edge;
  454.         int i_dir;
  455.         const int b_8x8_transform = h->mb.mb_transform_size[mb_xy];
  456.         /* cavlc + 8x8 transform stores nnz per 16 coeffs for the purpose of
  457.          * entropy coding, but per 64 coeffs for the purpose of deblocking */
  458.         if( !h->param.b_cabac && b_8x8_transform )
  459.         {
  460.             uint32_t *nnz = (uint32_t*)h->mb.non_zero_count[mb_xy];
  461.             if( nnz[0] ) nnz[0] = 0x01010101;
  462.             if( nnz[1] ) nnz[1] = 0x01010101;
  463.             if( nnz[2] ) nnz[2] = 0x01010101;
  464.             if( nnz[3] ) nnz[3] = 0x01010101;
  465.         }
  466.         /* i_dir == 0 -> vertical edge
  467.          * i_dir == 1 -> horizontal edge */
  468.         for( i_dir = 0; i_dir < 2; i_dir++ )
  469.         {
  470.             int i_start;
  471.             int i_qp, i_qpn;
  472.             i_start = (( i_dir == 0 && mb_x != 0 ) || ( i_dir == 1 && mb_y != 0 ) ) ? 0 : 1;
  473.             for( i_edge = i_start; i_edge < 4; i_edge++ )
  474.             {
  475.                 int mbn_xy  = i_edge > 0 ? mb_xy  : ( i_dir == 0 ? mb_xy  - 1 : mb_xy - h->mb.i_mb_stride );
  476.                 int mbn_8x8 = i_edge > 0 ? mb_8x8 : ( i_dir == 0 ? mb_8x8 - 2 : mb_8x8 - 2 * s8x8 );
  477.                 int mbn_4x4 = i_edge > 0 ? mb_4x4 : ( i_dir == 0 ? mb_4x4 - 4 : mb_4x4 - 4 * s4x4 );
  478.                 int bS[4];  /* filtering strength */
  479.                 /* *** Get bS for each 4px for the current edge *** */
  480.                 if( IS_INTRA( h->mb.type[mb_xy] ) || IS_INTRA( h->mb.type[mbn_xy] ) )
  481.                 {
  482.                     bS[0] = bS[1] = bS[2] = bS[3] = ( i_edge == 0 ? 4 : 3 );
  483.                 }
  484.                 else
  485.                 {
  486.                     int i;
  487.                     for( i = 0; i < 4; i++ )
  488.                     {
  489.                         int x  = i_dir == 0 ? i_edge : i;
  490.                         int y  = i_dir == 0 ? i      : i_edge;
  491.                         int xn = (x - (i_dir == 0 ? 1 : 0 ))&0x03;
  492.                         int yn = (y - (i_dir == 0 ? 0 : 1 ))&0x03;
  493.                         if( h->mb.non_zero_count[mb_xy][block_idx_xy[x][y]] != 0 ||
  494.                             h->mb.non_zero_count[mbn_xy][block_idx_xy[xn][yn]] != 0 )
  495.                         {
  496.                             bS[i] = 2;
  497.                         }
  498.                         else
  499.                         {
  500.                             /* FIXME: A given frame may occupy more than one position in
  501.                              * the reference list. So we should compare the frame numbers,
  502.                              * not the indices in the ref list.
  503.                              * No harm yet, as we don't generate that case.*/
  504.                             int i8p= mb_8x8+(x/2)+(y/2)*s8x8;
  505.                             int i8q= mbn_8x8+(xn/2)+(yn/2)*s8x8;
  506.                             int i4p= mb_4x4+x+y*s4x4;
  507.                             int i4q= mbn_4x4+xn+yn*s4x4;
  508.                             int l;
  509.                             bS[i] = 0;
  510.                             for( l = 0; l < 1 + (i_slice_type == SLICE_TYPE_B); l++ )
  511.                             {
  512.                                 if( h->mb.ref[l][i8p] != h->mb.ref[l][i8q] ||
  513.                                     abs( h->mb.mv[l][i4p][0] - h->mb.mv[l][i4q][0] ) >= 4 ||
  514.                                     abs( h->mb.mv[l][i4p][1] - h->mb.mv[l][i4q][1] ) >= 4 )
  515.                                 {
  516.                                     bS[i] = 1;
  517.                                     break;
  518.                                 }
  519.                             }
  520.                         }
  521.                     }
  522.                 }
  523.                 /* *** filter *** */
  524.                 /* Y plane */
  525.                 i_qp = h->mb.qp[mb_xy];
  526.                 i_qpn= h->mb.qp[mbn_xy];
  527.                 if( i_dir == 0 )
  528.                 {
  529.                     /* vertical edge */
  530.                     if( !b_8x8_transform || !(i_edge & 1) )
  531.                     {
  532.                         deblock_edge( h, &h->fdec->plane[0][16*mb_y * h->fdec->i_stride[0] + 16*mb_x + 4*i_edge],
  533.                                       h->fdec->i_stride[0], bS, (i_qp+i_qpn+1) >> 1, 0,
  534.                                       h->loopf.deblock_h_luma, h->loopf.deblock_h_luma_intra );
  535.                     }
  536.                     if( !(i_edge & 1) )
  537.                     {
  538.                         /* U/V planes */
  539.                         int i_qpc = ( i_chroma_qp_table[x264_clip3( i_qp + h->pps->i_chroma_qp_index_offset, 0, 51 )] +
  540.                                       i_chroma_qp_table[x264_clip3( i_qpn + h->pps->i_chroma_qp_index_offset, 0, 51 )] + 1 ) >> 1;
  541.                         deblock_edge( h, &h->fdec->plane[1][8*(mb_y*h->fdec->i_stride[1]+mb_x)+2*i_edge],
  542.                                       h->fdec->i_stride[1], bS, i_qpc, 1,
  543.                                       h->loopf.deblock_h_chroma, h->loopf.deblock_h_chroma_intra );
  544.                         deblock_edge( h, &h->fdec->plane[2][8*(mb_y*h->fdec->i_stride[2]+mb_x)+2*i_edge],
  545.                                       h->fdec->i_stride[2], bS, i_qpc, 1,
  546.                                       h->loopf.deblock_h_chroma, h->loopf.deblock_h_chroma_intra );
  547.                     }
  548.                 }
  549.                 else
  550.                 {
  551.                     /* horizontal edge */
  552.                     if( !b_8x8_transform || !(i_edge & 1) )
  553.                     {
  554.                         deblock_edge( h, &h->fdec->plane[0][(16*mb_y + 4*i_edge) * h->fdec->i_stride[0] + 16*mb_x],
  555.                                       h->fdec->i_stride[0], bS, (i_qp+i_qpn+1) >> 1, 0,
  556.                                       h->loopf.deblock_v_luma, h->loopf.deblock_v_luma_intra );
  557.                     }
  558.                     /* U/V planes */
  559.                     if( !(i_edge & 1) )
  560.                     {
  561.                         int i_qpc = ( i_chroma_qp_table[x264_clip3( i_qp + h->pps->i_chroma_qp_index_offset, 0, 51 )] +
  562.                                       i_chroma_qp_table[x264_clip3( i_qpn + h->pps->i_chroma_qp_index_offset, 0, 51 )] + 1 ) >> 1;
  563.                         deblock_edge( h, &h->fdec->plane[1][8*(mb_y*h->fdec->i_stride[1]+mb_x)+2*i_edge*h->fdec->i_stride[1]],
  564.                                       h->fdec->i_stride[1], bS, i_qpc, 1,
  565.                                       h->loopf.deblock_v_chroma, h->loopf.deblock_v_chroma_intra );
  566.                         deblock_edge( h, &h->fdec->plane[2][8*(mb_y*h->fdec->i_stride[2]+mb_x)+2*i_edge*h->fdec->i_stride[2]],
  567.                                       h->fdec->i_stride[2], bS, i_qpc, 1,
  568.                                       h->loopf.deblock_v_chroma, h->loopf.deblock_v_chroma_intra );
  569.                     }
  570.                 }
  571.             }
  572.         }
  573.         /* newt mb */
  574.         mb_x++;
  575.         if( mb_x >= h->sps->i_mb_width )
  576.         {
  577.             mb_x = 0;
  578.             mb_y++;
  579.         }
  580.     }
  581. }
  582. #ifdef HAVE_MMXEXT
  583. void x264_deblock_v_chroma_mmxext( uint8_t *pix, int stride, int alpha, int beta, int8_t *tc0 );
  584. void x264_deblock_h_chroma_mmxext( uint8_t *pix, int stride, int alpha, int beta, int8_t *tc0 );
  585. void x264_deblock_v_chroma_intra_mmxext( uint8_t *pix, int stride, int alpha, int beta );
  586. void x264_deblock_h_chroma_intra_mmxext( uint8_t *pix, int stride, int alpha, int beta );
  587. #endif
  588. #ifdef ARCH_X86_64
  589. void x264_deblock_v_luma_sse2( uint8_t *pix, int stride, int alpha, int beta, int8_t *tc0 );
  590. void x264_deblock_h_luma_sse2( uint8_t *pix, int stride, int alpha, int beta, int8_t *tc0 );
  591. #elif defined( HAVE_MMXEXT )
  592. void x264_deblock_h_luma_mmxext( uint8_t *pix, int stride, int alpha, int beta, int8_t *tc0 );
  593. void x264_deblock_v8_luma_mmxext( uint8_t *pix, int stride, int alpha, int beta, int8_t *tc0 );
  594. void x264_deblock_v_luma_mmxext( uint8_t *pix, int stride, int alpha, int beta, int8_t *tc0 )
  595. {
  596.     x264_deblock_v8_luma_mmxext( pix,   stride, alpha, beta, tc0   );
  597.     x264_deblock_v8_luma_mmxext( pix+8, stride, alpha, beta, tc0+2 );
  598. }
  599. #endif
  600. void x264_deblock_init( int cpu, x264_deblock_function_t *pf )
  601. {
  602.     pf->deblock_v_luma = deblock_v_luma_c;
  603.     pf->deblock_h_luma = deblock_h_luma_c;
  604.     pf->deblock_v_chroma = deblock_v_chroma_c;
  605.     pf->deblock_h_chroma = deblock_h_chroma_c;
  606.     pf->deblock_v_luma_intra = deblock_v_luma_intra_c;
  607.     pf->deblock_h_luma_intra = deblock_h_luma_intra_c;
  608.     pf->deblock_v_chroma_intra = deblock_v_chroma_intra_c;
  609.     pf->deblock_h_chroma_intra = deblock_h_chroma_intra_c;
  610. #ifdef HAVE_MMXEXT
  611.     if( cpu&X264_CPU_MMXEXT )
  612.     {
  613.         pf->deblock_v_chroma = x264_deblock_v_chroma_mmxext;
  614.         pf->deblock_h_chroma = x264_deblock_h_chroma_mmxext;
  615.         pf->deblock_v_chroma_intra = x264_deblock_v_chroma_intra_mmxext;
  616.         pf->deblock_h_chroma_intra = x264_deblock_h_chroma_intra_mmxext;
  617. #ifdef ARCH_X86_64
  618.         if( cpu&X264_CPU_SSE2 )
  619.         {
  620.             pf->deblock_v_luma = x264_deblock_v_luma_sse2;
  621.             pf->deblock_h_luma = x264_deblock_h_luma_sse2;
  622.         }
  623. #else
  624.         pf->deblock_v_luma = x264_deblock_v_luma_mmxext;
  625.         pf->deblock_h_luma = x264_deblock_h_luma_mmxext;
  626. #endif
  627.     }
  628. #endif
  629. }