开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 多媒体 > Audio > 查看源码
pixel.c
资源名称：chapter15.rar　[点击查看]
上传用户：hjq518
上传日期：2021-12-09
资源大小：5084k
文件大小：25k
源码类别：

Audio

开发平台：

Visual C++

pixel.c：源码内容
  1. /*****************************************************************************
  2.  * pixel.c: h264 encoder
  3.  *****************************************************************************
  4.  * Copyright (C) 2003 Laurent Aimar
  5.  * $Id: pixel.c,v 1.1 2004/06/03 19:27:07 fenrir Exp $
  6.  *
  7.  * Authors: Eric Petit <titer@m0k.org>
  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.  *****************************************************************************/
  24. #ifdef SYS_LINUX
  25. #include <altivec.h>
  26. #endif
  28. #include "common/common.h"
  29. #include "ppccommon.h"
  31. /***********************************************************************
  32.  * SAD routines
  33.  **********************************************************************/
  35. #define PIXEL_SAD_ALTIVEC( name, lx, ly, a, b )        
  36. static int name( uint8_t *pix1, int i_pix1,            
  37.                  uint8_t *pix2, int i_pix2 )           
  38. {                                                      
  39.     int y;                                             
  40.     DECLARE_ALIGNED( int, sum, 16 );                   
  41.                                                        
  42.     LOAD_ZERO;                                         
  43.     PREP_LOAD;                                         
  44.     vec_u8_t  pix1v, pix2v;                            
  45.     vec_s32_t sumv = zero_s32v;                        
  46.     for( y = 0; y < ly; y++ )                          
  47.     {                                                  
  48.         VEC_LOAD( pix1, pix1v, lx, vec_u8_t );         
  49.         VEC_LOAD( pix2, pix2v, lx, vec_u8_t );         
  50.         sumv = (vec_s32_t) vec_sum4s(                  
  51.                    vec_sub( vec_max( pix1v, pix2v ),   
  52.                             vec_min( pix1v, pix2v ) ), 
  53.                    (vec_u32_t) sumv );                 
  54.         pix1 += i_pix1;                                
  55.         pix2 += i_pix2;                                
  56.     }                                                  
  57.     sumv = vec_sum##a( sumv, zero_s32v );              
  58.     sumv = vec_splat( sumv, b );                       
  59.     vec_ste( sumv, 0, &sum );                          
  60.     return sum;                                        
  61. }
  63. PIXEL_SAD_ALTIVEC( pixel_sad_16x16_altivec, 16, 16, s,  3 )
  64. PIXEL_SAD_ALTIVEC( pixel_sad_8x16_altivec,  8,  16, 2s, 1 )
  65. PIXEL_SAD_ALTIVEC( pixel_sad_16x8_altivec,  16, 8,  s,  3 )
  66. PIXEL_SAD_ALTIVEC( pixel_sad_8x8_altivec,   8,  8,  2s, 1 )
  68. /***********************************************************************
  69.  * SATD routines
  70.  **********************************************************************/
  72. /***********************************************************************
  73.  * VEC_HADAMAR
  74.  ***********************************************************************
  75.  * b[0] = a[0] + a[1] + a[2] + a[3]
  76.  * b[1] = a[0] + a[1] - a[2] - a[3]
  77.  * b[2] = a[0] - a[1] - a[2] + a[3]
  78.  * b[3] = a[0] - a[1] + a[2] - a[3]
  79.  **********************************************************************/
  80. #define VEC_HADAMAR(a0,a1,a2,a3,b0,b1,b2,b3) 
  81.     b2 = vec_add( a0, a1 ); 
  82.     b3 = vec_add( a2, a3 ); 
  83.     a0 = vec_sub( a0, a1 ); 
  84.     a2 = vec_sub( a2, a3 ); 
  85.     b0 = vec_add( b2, b3 ); 
  86.     b1 = vec_sub( b2, b3 ); 
  87.     b2 = vec_sub( a0, a2 ); 
  88.     b3 = vec_add( a0, a2 )
  90. /***********************************************************************
  91.  * VEC_ABS
  92.  ***********************************************************************
  93.  * a: s16v
  94.  *
  95.  * a = abs(a)
  96.  * 
  97.  * Call vec_sub()/vec_max() instead of vec_abs() because vec_abs()
  98.  * actually also calls vec_splat(0), but we already have a null vector.
  99.  **********************************************************************/
  100. #define VEC_ABS(a) 
  101.     pix1v = vec_sub( zero_s16v, a ); 
  102.     a     = vec_max( a, pix1v ); 
  104. /***********************************************************************
  105.  * VEC_ADD_ABS
  106.  ***********************************************************************
  107.  * a:    s16v
  108.  * b, c: s32v
  109.  *
  110.  * c[i] = abs(a[2*i]) + abs(a[2*i+1]) + [bi]
  111.  **********************************************************************/
  112. #define VEC_ADD_ABS(a,b,c) 
  113.     VEC_ABS( a ); 
  114.     c = vec_sum4s( a, b )
  116. /***********************************************************************
  117.  * SATD 4x4
  118.  **********************************************************************/
  119. static int pixel_satd_4x4_altivec( uint8_t *pix1, int i_pix1,
  120.                                    uint8_t *pix2, int i_pix2 )
  121. {
  122.     DECLARE_ALIGNED( int, i_satd, 16 );
  124.     PREP_DIFF;
  125.     vec_s16_t diff0v, diff1v, diff2v, diff3v;
  126.     vec_s16_t temp0v, temp1v, temp2v, temp3v;
  127.     vec_s32_t satdv;
  129.     VEC_DIFF_H( pix1, i_pix1, pix2, i_pix2, 4, diff0v );
  130.     VEC_DIFF_H( pix1, i_pix1, pix2, i_pix2, 4, diff1v );
  131.     VEC_DIFF_H( pix1, i_pix1, pix2, i_pix2, 4, diff2v );
  132.     VEC_DIFF_H( pix1, i_pix1, pix2, i_pix2, 4, diff3v );
  134.     /* Hadamar H */
  135.     VEC_HADAMAR( diff0v, diff1v, diff2v, diff3v,
  136.                  temp0v, temp1v, temp2v, temp3v );
  137.     
  138.     VEC_TRANSPOSE_4( temp0v, temp1v, temp2v, temp3v,
  139.                      diff0v, diff1v, diff2v, diff3v );
  140.     /* Hadamar V */
  141.     VEC_HADAMAR( diff0v, diff1v, diff2v, diff3v,
  142.                  temp0v, temp1v, temp2v, temp3v );
  144.     VEC_ADD_ABS( temp0v, zero_s32v, satdv );
  145.     VEC_ADD_ABS( temp1v, satdv,     satdv );
  146.     VEC_ADD_ABS( temp2v, satdv,     satdv );
  147.     VEC_ADD_ABS( temp3v, satdv,     satdv );
  149.     satdv = vec_sum2s( satdv, zero_s32v );
  150.     satdv = vec_splat( satdv, 1 );
  151.     vec_ste( satdv, 0, &i_satd );
  153.     return i_satd / 2;
  154. }
  156. /***********************************************************************
  157.  * SATD 4x8
  158.  **********************************************************************/
  159. static int pixel_satd_4x8_altivec( uint8_t *pix1, int i_pix1,
  160.                                    uint8_t *pix2, int i_pix2 )
  161. {
  162.     DECLARE_ALIGNED( int, i_satd, 16 );
  164.     PREP_DIFF;
  165.     vec_s16_t diff0v, diff1v, diff2v, diff3v;
  166.     vec_s16_t temp0v, temp1v, temp2v, temp3v;
  167.     vec_s32_t satdv;
  169.     VEC_DIFF_H( pix1, i_pix1, pix2, i_pix2, 4, diff0v );
  170.     VEC_DIFF_H( pix1, i_pix1, pix2, i_pix2, 4, diff1v );
  171.     VEC_DIFF_H( pix1, i_pix1, pix2, i_pix2, 4, diff2v );
  172.     VEC_DIFF_H( pix1, i_pix1, pix2, i_pix2, 4, diff3v );
  173.     VEC_HADAMAR( diff0v, diff1v, diff2v, diff3v,
  174.                  temp0v, temp1v, temp2v, temp3v );
  175.     VEC_TRANSPOSE_4( temp0v, temp1v, temp2v, temp3v,
  176.                      diff0v, diff1v, diff2v, diff3v );
  177.     VEC_HADAMAR( diff0v, diff1v, diff2v, diff3v,
  178.                  temp0v, temp1v, temp2v, temp3v );
  179.     VEC_ADD_ABS( temp0v, zero_s32v, satdv );
  180.     VEC_ADD_ABS( temp1v, satdv,     satdv );
  181.     VEC_ADD_ABS( temp2v, satdv,     satdv );
  182.     VEC_ADD_ABS( temp3v, satdv,     satdv );
  184.     VEC_DIFF_H( pix1, i_pix1, pix2, i_pix2, 4, diff0v );
  185.     VEC_DIFF_H( pix1, i_pix1, pix2, i_pix2, 4, diff1v );
  186.     VEC_DIFF_H( pix1, i_pix1, pix2, i_pix2, 4, diff2v );
  187.     VEC_DIFF_H( pix1, i_pix1, pix2, i_pix2, 4, diff3v );
  188.     VEC_HADAMAR( diff0v, diff1v, diff2v, diff3v,
  189.                  temp0v, temp1v, temp2v, temp3v );
  190.     VEC_TRANSPOSE_4( temp0v, temp1v, temp2v, temp3v,
  191.                      diff0v, diff1v, diff2v, diff3v );
  192.     VEC_HADAMAR( diff0v, diff1v, diff2v, diff3v,
  193.                  temp0v, temp1v, temp2v, temp3v );
  194.     VEC_ADD_ABS( temp0v, satdv,     satdv );
  195.     VEC_ADD_ABS( temp1v, satdv,     satdv );
  196.     VEC_ADD_ABS( temp2v, satdv,     satdv );
  197.     VEC_ADD_ABS( temp3v, satdv,     satdv );
  199.     satdv = vec_sum2s( satdv, zero_s32v );
  200.     satdv = vec_splat( satdv, 1 );
  201.     vec_ste( satdv, 0, &i_satd );
  203.     return i_satd / 2;
  204. }
  206. /***********************************************************************
  207.  * SATD 8x4
  208.  **********************************************************************/
  209. static int pixel_satd_8x4_altivec( uint8_t *pix1, int i_pix1,
  210.                                    uint8_t *pix2, int i_pix2 )
  211. {
  212.     DECLARE_ALIGNED( int, i_satd, 16 );
  214.     PREP_DIFF;
  215.     vec_s16_t diff0v, diff1v, diff2v, diff3v,
  216.               diff4v, diff5v, diff6v, diff7v;
  217.     vec_s16_t temp0v, temp1v, temp2v, temp3v,
  218.               temp4v, temp5v, temp6v, temp7v;
  219.     vec_s32_t satdv;
  221.     VEC_DIFF_H( pix1, i_pix1, pix2, i_pix2, 8, diff0v );
  222.     VEC_DIFF_H( pix1, i_pix1, pix2, i_pix2, 8, diff1v );
  223.     VEC_DIFF_H( pix1, i_pix1, pix2, i_pix2, 8, diff2v );
  224.     VEC_DIFF_H( pix1, i_pix1, pix2, i_pix2, 8, diff3v );
  226.     VEC_HADAMAR( diff0v, diff1v, diff2v, diff3v,
  227.                  temp0v, temp1v, temp2v, temp3v );
  228.     /* This causes warnings because temp4v...temp7v haven't be set,
  229.        but we don't care */
  230.     VEC_TRANSPOSE_8( temp0v, temp1v, temp2v, temp3v,
  231.                      temp4v, temp5v, temp6v, temp7v,
  232.                      diff0v, diff1v, diff2v, diff3v,
  233.                      diff4v, diff5v, diff6v, diff7v );
  234.     VEC_HADAMAR( diff0v, diff1v, diff2v, diff3v,
  235.                  temp0v, temp1v, temp2v, temp3v );
  236.     VEC_HADAMAR( diff4v, diff5v, diff6v, diff7v,
  237.                  temp4v, temp5v, temp6v, temp7v );
  239.     VEC_ADD_ABS( temp0v, zero_s32v, satdv );
  240.     VEC_ADD_ABS( temp1v, satdv,     satdv );
  241.     VEC_ADD_ABS( temp2v, satdv,     satdv );
  242.     VEC_ADD_ABS( temp3v, satdv,     satdv );
  243.     VEC_ADD_ABS( temp4v, satdv,     satdv );
  244.     VEC_ADD_ABS( temp5v, satdv,     satdv );
  245.     VEC_ADD_ABS( temp6v, satdv,     satdv );
  246.     VEC_ADD_ABS( temp7v, satdv,     satdv );
  248.     satdv = vec_sum2s( satdv, zero_s32v );
  249.     satdv = vec_splat( satdv, 1 );
  250.     vec_ste( satdv, 0, &i_satd );
  252.     return i_satd / 2;
  253. }
  255. /***********************************************************************
  256.  * SATD 8x8
  257.  **********************************************************************/
  258. static int pixel_satd_8x8_altivec( uint8_t *pix1, int i_pix1,
  259.                                    uint8_t *pix2, int i_pix2 )
  260. {
  261.     DECLARE_ALIGNED( int, i_satd, 16 );
  263.     PREP_DIFF;
  264.     vec_s16_t diff0v, diff1v, diff2v, diff3v,
  265.               diff4v, diff5v, diff6v, diff7v;
  266.     vec_s16_t temp0v, temp1v, temp2v, temp3v,
  267.               temp4v, temp5v, temp6v, temp7v;
  268.     vec_s32_t satdv;
  270.     VEC_DIFF_H( pix1, i_pix1, pix2, i_pix2, 8, diff0v );
  271.     VEC_DIFF_H( pix1, i_pix1, pix2, i_pix2, 8, diff1v );
  272.     VEC_DIFF_H( pix1, i_pix1, pix2, i_pix2, 8, diff2v );
  273.     VEC_DIFF_H( pix1, i_pix1, pix2, i_pix2, 8, diff3v );
  274.     VEC_DIFF_H( pix1, i_pix1, pix2, i_pix2, 8, diff4v );
  275.     VEC_DIFF_H( pix1, i_pix1, pix2, i_pix2, 8, diff5v );
  276.     VEC_DIFF_H( pix1, i_pix1, pix2, i_pix2, 8, diff6v );
  277.     VEC_DIFF_H( pix1, i_pix1, pix2, i_pix2, 8, diff7v );
  279.     VEC_HADAMAR( diff0v, diff1v, diff2v, diff3v,
  280.                  temp0v, temp1v, temp2v, temp3v );
  281.     VEC_HADAMAR( diff4v, diff5v, diff6v, diff7v,
  282.                  temp4v, temp5v, temp6v, temp7v );
  284.     VEC_TRANSPOSE_8( temp0v, temp1v, temp2v, temp3v,
  285.                      temp4v, temp5v, temp6v, temp7v,
  286.                      diff0v, diff1v, diff2v, diff3v,
  287.                      diff4v, diff5v, diff6v, diff7v );
  289.     VEC_HADAMAR( diff0v, diff1v, diff2v, diff3v,
  290.                  temp0v, temp1v, temp2v, temp3v );
  291.     VEC_HADAMAR( diff4v, diff5v, diff6v, diff7v,
  292.                  temp4v, temp5v, temp6v, temp7v );
  294.     VEC_ADD_ABS( temp0v, zero_s32v, satdv );
  295.     VEC_ADD_ABS( temp1v, satdv,     satdv );
  296.     VEC_ADD_ABS( temp2v, satdv,     satdv );
  297.     VEC_ADD_ABS( temp3v, satdv,     satdv );
  298.     VEC_ADD_ABS( temp4v, satdv,     satdv );
  299.     VEC_ADD_ABS( temp5v, satdv,     satdv );
  300.     VEC_ADD_ABS( temp6v, satdv,     satdv );
  301.     VEC_ADD_ABS( temp7v, satdv,     satdv );
  303.     satdv = vec_sums( satdv, zero_s32v );
  304.     satdv = vec_splat( satdv, 3 );
  305.     vec_ste( satdv, 0, &i_satd );
  307.     return i_satd / 2;
  308. }
  310. /***********************************************************************
  311.  * SATD 8x16
  312.  **********************************************************************/
  313. static int pixel_satd_8x16_altivec( uint8_t *pix1, int i_pix1,
  314.                                     uint8_t *pix2, int i_pix2 )
  315. {
  316.     DECLARE_ALIGNED( int, i_satd, 16 );
  318.     PREP_DIFF;
  319.     vec_s16_t diff0v, diff1v, diff2v, diff3v,
  320.               diff4v, diff5v, diff6v, diff7v;
  321.     vec_s16_t temp0v, temp1v, temp2v, temp3v,
  322.               temp4v, temp5v, temp6v, temp7v;
  323.     vec_s32_t satdv;
  325.     VEC_DIFF_H( pix1, i_pix1, pix2, i_pix2, 8, diff0v );
  326.     VEC_DIFF_H( pix1, i_pix1, pix2, i_pix2, 8, diff1v );
  327.     VEC_DIFF_H( pix1, i_pix1, pix2, i_pix2, 8, diff2v );
  328.     VEC_DIFF_H( pix1, i_pix1, pix2, i_pix2, 8, diff3v );
  329.     VEC_DIFF_H( pix1, i_pix1, pix2, i_pix2, 8, diff4v );
  330.     VEC_DIFF_H( pix1, i_pix1, pix2, i_pix2, 8, diff5v );
  331.     VEC_DIFF_H( pix1, i_pix1, pix2, i_pix2, 8, diff6v );
  332.     VEC_DIFF_H( pix1, i_pix1, pix2, i_pix2, 8, diff7v );
  333.     VEC_HADAMAR( diff0v, diff1v, diff2v, diff3v,
  334.                  temp0v, temp1v, temp2v, temp3v );
  335.     VEC_HADAMAR( diff4v, diff5v, diff6v, diff7v,
  336.                  temp4v, temp5v, temp6v, temp7v );
  337.     VEC_TRANSPOSE_8( temp0v, temp1v, temp2v, temp3v,
  338.                      temp4v, temp5v, temp6v, temp7v,
  339.                      diff0v, diff1v, diff2v, diff3v,
  340.                      diff4v, diff5v, diff6v, diff7v );
  341.     VEC_HADAMAR( diff0v, diff1v, diff2v, diff3v,
  342.                  temp0v, temp1v, temp2v, temp3v );
  343.     VEC_HADAMAR( diff4v, diff5v, diff6v, diff7v,
  344.                  temp4v, temp5v, temp6v, temp7v );
  345.     VEC_ADD_ABS( temp0v, zero_s32v, satdv );
  346.     VEC_ADD_ABS( temp1v, satdv,     satdv );
  347.     VEC_ADD_ABS( temp2v, satdv,     satdv );
  348.     VEC_ADD_ABS( temp3v, satdv,     satdv );
  349.     VEC_ADD_ABS( temp4v, satdv,     satdv );
  350.     VEC_ADD_ABS( temp5v, satdv,     satdv );
  351.     VEC_ADD_ABS( temp6v, satdv,     satdv );
  352.     VEC_ADD_ABS( temp7v, satdv,     satdv );
  354.     VEC_DIFF_H( pix1, i_pix1, pix2, i_pix2, 8, diff0v );
  355.     VEC_DIFF_H( pix1, i_pix1, pix2, i_pix2, 8, diff1v );
  356.     VEC_DIFF_H( pix1, i_pix1, pix2, i_pix2, 8, diff2v );
  357.     VEC_DIFF_H( pix1, i_pix1, pix2, i_pix2, 8, diff3v );
  358.     VEC_DIFF_H( pix1, i_pix1, pix2, i_pix2, 8, diff4v );
  359.     VEC_DIFF_H( pix1, i_pix1, pix2, i_pix2, 8, diff5v );
  360.     VEC_DIFF_H( pix1, i_pix1, pix2, i_pix2, 8, diff6v );
  361.     VEC_DIFF_H( pix1, i_pix1, pix2, i_pix2, 8, diff7v );
  362.     VEC_HADAMAR( diff0v, diff1v, diff2v, diff3v,
  363.                  temp0v, temp1v, temp2v, temp3v );
  364.     VEC_HADAMAR( diff4v, diff5v, diff6v, diff7v,
  365.                  temp4v, temp5v, temp6v, temp7v );
  366.     VEC_TRANSPOSE_8( temp0v, temp1v, temp2v, temp3v,
  367.                      temp4v, temp5v, temp6v, temp7v,
  368.                      diff0v, diff1v, diff2v, diff3v,
  369.                      diff4v, diff5v, diff6v, diff7v );
  370.     VEC_HADAMAR( diff0v, diff1v, diff2v, diff3v,
  371.                  temp0v, temp1v, temp2v, temp3v );
  372.     VEC_HADAMAR( diff4v, diff5v, diff6v, diff7v,
  373.                  temp4v, temp5v, temp6v, temp7v );
  374.     VEC_ADD_ABS( temp0v, satdv,     satdv );
  375.     VEC_ADD_ABS( temp1v, satdv,     satdv );
  376.     VEC_ADD_ABS( temp2v, satdv,     satdv );
  377.     VEC_ADD_ABS( temp3v, satdv,     satdv );
  378.     VEC_ADD_ABS( temp4v, satdv,     satdv );
  379.     VEC_ADD_ABS( temp5v, satdv,     satdv );
  380.     VEC_ADD_ABS( temp6v, satdv,     satdv );
  381.     VEC_ADD_ABS( temp7v, satdv,     satdv );
  383.     satdv = vec_sums( satdv, zero_s32v );
  384.     satdv = vec_splat( satdv, 3 );
  385.     vec_ste( satdv, 0, &i_satd );
  387.     return i_satd / 2;
  388. }
  390. /***********************************************************************
  391.  * SATD 16x8
  392.  **********************************************************************/
  393. static int pixel_satd_16x8_altivec( uint8_t *pix1, int i_pix1,
  394.                                     uint8_t *pix2, int i_pix2 )
  395. {
  396.     DECLARE_ALIGNED( int, i_satd, 16 );
  398.     LOAD_ZERO;
  399.     PREP_LOAD;
  400.     vec_s32_t satdv;
  401.     vec_s16_t pix1v, pix2v;
  402.     vec_s16_t diffh0v, diffh1v, diffh2v, diffh3v,
  403.               diffh4v, diffh5v, diffh6v, diffh7v;
  404.     vec_s16_t diffl0v, diffl1v, diffl2v, diffl3v,
  405.               diffl4v, diffl5v, diffl6v, diffl7v;
  406.     vec_s16_t temp0v, temp1v, temp2v, temp3v,
  407.               temp4v, temp5v, temp6v, temp7v;
  409.     VEC_DIFF_HL( pix1, i_pix1, pix2, i_pix2, diffh0v, diffl0v );
  410.     VEC_DIFF_HL( pix1, i_pix1, pix2, i_pix2, diffh1v, diffl1v );
  411.     VEC_DIFF_HL( pix1, i_pix1, pix2, i_pix2, diffh2v, diffl2v );
  412.     VEC_DIFF_HL( pix1, i_pix1, pix2, i_pix2, diffh3v, diffl3v );
  413.     VEC_DIFF_HL( pix1, i_pix1, pix2, i_pix2, diffh4v, diffl4v );
  414.     VEC_DIFF_HL( pix1, i_pix1, pix2, i_pix2, diffh5v, diffl5v );
  415.     VEC_DIFF_HL( pix1, i_pix1, pix2, i_pix2, diffh6v, diffl6v );
  416.     VEC_DIFF_HL( pix1, i_pix1, pix2, i_pix2, diffh7v, diffl7v );
  418.     VEC_HADAMAR( diffh0v, diffh1v, diffh2v, diffh3v,
  419.                  temp0v, temp1v, temp2v, temp3v );
  420.     VEC_HADAMAR( diffh4v, diffh5v, diffh6v, diffh7v,
  421.                  temp4v, temp5v, temp6v, temp7v );
  423.     VEC_TRANSPOSE_8( temp0v, temp1v, temp2v, temp3v,
  424.                      temp4v, temp5v, temp6v, temp7v,
  425.                      diffh0v, diffh1v, diffh2v, diffh3v,
  426.                      diffh4v, diffh5v, diffh6v, diffh7v );
  428.     VEC_HADAMAR( diffh0v, diffh1v, diffh2v, diffh3v,
  429.                  temp0v, temp1v, temp2v, temp3v );
  430.     VEC_HADAMAR( diffh4v, diffh5v, diffh6v, diffh7v,
  431.                  temp4v, temp5v, temp6v, temp7v );
  433.     VEC_ADD_ABS( temp0v, zero_s32v, satdv );
  434.     VEC_ADD_ABS( temp1v, satdv,     satdv );
  435.     VEC_ADD_ABS( temp2v, satdv,     satdv );
  436.     VEC_ADD_ABS( temp3v, satdv,     satdv );
  437.     VEC_ADD_ABS( temp4v, satdv,     satdv );
  438.     VEC_ADD_ABS( temp5v, satdv,     satdv );
  439.     VEC_ADD_ABS( temp6v, satdv,     satdv );
  440.     VEC_ADD_ABS( temp7v, satdv,     satdv );
  442.     VEC_HADAMAR( diffl0v, diffl1v, diffl2v, diffl3v,
  443.                  temp0v, temp1v, temp2v, temp3v );
  444.     VEC_HADAMAR( diffl4v, diffl5v, diffl6v, diffl7v,
  445.                  temp4v, temp5v, temp6v, temp7v );
  447.     VEC_TRANSPOSE_8( temp0v, temp1v, temp2v, temp3v,
  448.                      temp4v, temp5v, temp6v, temp7v,
  449.                      diffl0v, diffl1v, diffl2v, diffl3v,
  450.                      diffl4v, diffl5v, diffl6v, diffl7v );
  452.     VEC_HADAMAR( diffl0v, diffl1v, diffl2v, diffl3v,
  453.                  temp0v, temp1v, temp2v, temp3v );
  454.     VEC_HADAMAR( diffl4v, diffl5v, diffl6v, diffl7v,
  455.                  temp4v, temp5v, temp6v, temp7v );
  457.     VEC_ADD_ABS( temp0v, satdv,     satdv );
  458.     VEC_ADD_ABS( temp1v, satdv,     satdv );
  459.     VEC_ADD_ABS( temp2v, satdv,     satdv );
  460.     VEC_ADD_ABS( temp3v, satdv,     satdv );
  461.     VEC_ADD_ABS( temp4v, satdv,     satdv );
  462.     VEC_ADD_ABS( temp5v, satdv,     satdv );
  463.     VEC_ADD_ABS( temp6v, satdv,     satdv );
  464.     VEC_ADD_ABS( temp7v, satdv,     satdv );
  466.     satdv = vec_sums( satdv, zero_s32v );
  467.     satdv = vec_splat( satdv, 3 );
  468.     vec_ste( satdv, 0, &i_satd );
  470.     return i_satd / 2;
  471. }
  473. /***********************************************************************
  474.  * SATD 16x16
  475.  **********************************************************************/
  476. static int pixel_satd_16x16_altivec( uint8_t *pix1, int i_pix1,
  477.                                      uint8_t *pix2, int i_pix2 )
  478. {
  479.     DECLARE_ALIGNED( int, i_satd, 16 );
  481.     LOAD_ZERO;
  482.     PREP_LOAD;
  483.     vec_s32_t satdv;
  484.     vec_s16_t pix1v, pix2v;
  485.     vec_s16_t diffh0v, diffh1v, diffh2v, diffh3v,
  486.               diffh4v, diffh5v, diffh6v, diffh7v;
  487.     vec_s16_t diffl0v, diffl1v, diffl2v, diffl3v,
  488.               diffl4v, diffl5v, diffl6v, diffl7v;
  489.     vec_s16_t temp0v, temp1v, temp2v, temp3v,
  490.               temp4v, temp5v, temp6v, temp7v;
  492.     VEC_DIFF_HL( pix1, i_pix1, pix2, i_pix2, diffh0v, diffl0v );
  493.     VEC_DIFF_HL( pix1, i_pix1, pix2, i_pix2, diffh1v, diffl1v );
  494.     VEC_DIFF_HL( pix1, i_pix1, pix2, i_pix2, diffh2v, diffl2v );
  495.     VEC_DIFF_HL( pix1, i_pix1, pix2, i_pix2, diffh3v, diffl3v );
  496.     VEC_DIFF_HL( pix1, i_pix1, pix2, i_pix2, diffh4v, diffl4v );
  497.     VEC_DIFF_HL( pix1, i_pix1, pix2, i_pix2, diffh5v, diffl5v );
  498.     VEC_DIFF_HL( pix1, i_pix1, pix2, i_pix2, diffh6v, diffl6v );
  499.     VEC_DIFF_HL( pix1, i_pix1, pix2, i_pix2, diffh7v, diffl7v );
  500.     VEC_HADAMAR( diffh0v, diffh1v, diffh2v, diffh3v,
  501.                  temp0v, temp1v, temp2v, temp3v );
  502.     VEC_HADAMAR( diffh4v, diffh5v, diffh6v, diffh7v,
  503.                  temp4v, temp5v, temp6v, temp7v );
  504.     VEC_TRANSPOSE_8( temp0v, temp1v, temp2v, temp3v,
  505.                      temp4v, temp5v, temp6v, temp7v,
  506.                      diffh0v, diffh1v, diffh2v, diffh3v,
  507.                      diffh4v, diffh5v, diffh6v, diffh7v );
  508.     VEC_HADAMAR( diffh0v, diffh1v, diffh2v, diffh3v,
  509.                  temp0v, temp1v, temp2v, temp3v );
  510.     VEC_HADAMAR( diffh4v, diffh5v, diffh6v, diffh7v,
  511.                  temp4v, temp5v, temp6v, temp7v );
  512.     VEC_ADD_ABS( temp0v, zero_s32v, satdv );
  513.     VEC_ADD_ABS( temp1v, satdv,     satdv );
  514.     VEC_ADD_ABS( temp2v, satdv,     satdv );
  515.     VEC_ADD_ABS( temp3v, satdv,     satdv );
  516.     VEC_ADD_ABS( temp4v, satdv,     satdv );
  517.     VEC_ADD_ABS( temp5v, satdv,     satdv );
  518.     VEC_ADD_ABS( temp6v, satdv,     satdv );
  519.     VEC_ADD_ABS( temp7v, satdv,     satdv );
  520.     VEC_HADAMAR( diffl0v, diffl1v, diffl2v, diffl3v,
  521.                  temp0v, temp1v, temp2v, temp3v );
  522.     VEC_HADAMAR( diffl4v, diffl5v, diffl6v, diffl7v,
  523.                  temp4v, temp5v, temp6v, temp7v );
  524.     VEC_TRANSPOSE_8( temp0v, temp1v, temp2v, temp3v,
  525.                      temp4v, temp5v, temp6v, temp7v,
  526.                      diffl0v, diffl1v, diffl2v, diffl3v,
  527.                      diffl4v, diffl5v, diffl6v, diffl7v );
  528.     VEC_HADAMAR( diffl0v, diffl1v, diffl2v, diffl3v,
  529.                  temp0v, temp1v, temp2v, temp3v );
  530.     VEC_HADAMAR( diffl4v, diffl5v, diffl6v, diffl7v,
  531.                  temp4v, temp5v, temp6v, temp7v );
  532.     VEC_ADD_ABS( temp0v, satdv,     satdv );
  533.     VEC_ADD_ABS( temp1v, satdv,     satdv );
  534.     VEC_ADD_ABS( temp2v, satdv,     satdv );
  535.     VEC_ADD_ABS( temp3v, satdv,     satdv );
  536.     VEC_ADD_ABS( temp4v, satdv,     satdv );
  537.     VEC_ADD_ABS( temp5v, satdv,     satdv );
  538.     VEC_ADD_ABS( temp6v, satdv,     satdv );
  539.     VEC_ADD_ABS( temp7v, satdv,     satdv );
  541.     VEC_DIFF_HL( pix1, i_pix1, pix2, i_pix2, diffh0v, diffl0v );
  542.     VEC_DIFF_HL( pix1, i_pix1, pix2, i_pix2, diffh1v, diffl1v );
  543.     VEC_DIFF_HL( pix1, i_pix1, pix2, i_pix2, diffh2v, diffl2v );
  544.     VEC_DIFF_HL( pix1, i_pix1, pix2, i_pix2, diffh3v, diffl3v );
  545.     VEC_DIFF_HL( pix1, i_pix1, pix2, i_pix2, diffh4v, diffl4v );
  546.     VEC_DIFF_HL( pix1, i_pix1, pix2, i_pix2, diffh5v, diffl5v );
  547.     VEC_DIFF_HL( pix1, i_pix1, pix2, i_pix2, diffh6v, diffl6v );
  548.     VEC_DIFF_HL( pix1, i_pix1, pix2, i_pix2, diffh7v, diffl7v );
  549.     VEC_HADAMAR( diffh0v, diffh1v, diffh2v, diffh3v,
  550.                  temp0v, temp1v, temp2v, temp3v );
  551.     VEC_HADAMAR( diffh4v, diffh5v, diffh6v, diffh7v,
  552.                  temp4v, temp5v, temp6v, temp7v );
  553.     VEC_TRANSPOSE_8( temp0v, temp1v, temp2v, temp3v,
  554.                      temp4v, temp5v, temp6v, temp7v,
  555.                      diffh0v, diffh1v, diffh2v, diffh3v,
  556.                      diffh4v, diffh5v, diffh6v, diffh7v );
  557.     VEC_HADAMAR( diffh0v, diffh1v, diffh2v, diffh3v,
  558.                  temp0v, temp1v, temp2v, temp3v );
  559.     VEC_HADAMAR( diffh4v, diffh5v, diffh6v, diffh7v,
  560.                  temp4v, temp5v, temp6v, temp7v );
  561.     VEC_ADD_ABS( temp0v, satdv,     satdv );
  562.     VEC_ADD_ABS( temp1v, satdv,     satdv );
  563.     VEC_ADD_ABS( temp2v, satdv,     satdv );
  564.     VEC_ADD_ABS( temp3v, satdv,     satdv );
  565.     VEC_ADD_ABS( temp4v, satdv,     satdv );
  566.     VEC_ADD_ABS( temp5v, satdv,     satdv );
  567.     VEC_ADD_ABS( temp6v, satdv,     satdv );
  568.     VEC_ADD_ABS( temp7v, satdv,     satdv );
  569.     VEC_HADAMAR( diffl0v, diffl1v, diffl2v, diffl3v,
  570.                  temp0v, temp1v, temp2v, temp3v );
  571.     VEC_HADAMAR( diffl4v, diffl5v, diffl6v, diffl7v,
  572.                  temp4v, temp5v, temp6v, temp7v );
  573.     VEC_TRANSPOSE_8( temp0v, temp1v, temp2v, temp3v,
  574.                      temp4v, temp5v, temp6v, temp7v,
  575.                      diffl0v, diffl1v, diffl2v, diffl3v,
  576.                      diffl4v, diffl5v, diffl6v, diffl7v );
  577.     VEC_HADAMAR( diffl0v, diffl1v, diffl2v, diffl3v,
  578.                  temp0v, temp1v, temp2v, temp3v );
  579.     VEC_HADAMAR( diffl4v, diffl5v, diffl6v, diffl7v,
  580.                  temp4v, temp5v, temp6v, temp7v );
  581.     VEC_ADD_ABS( temp0v, satdv,     satdv );
  582.     VEC_ADD_ABS( temp1v, satdv,     satdv );
  583.     VEC_ADD_ABS( temp2v, satdv,     satdv );
  584.     VEC_ADD_ABS( temp3v, satdv,     satdv );
  585.     VEC_ADD_ABS( temp4v, satdv,     satdv );
  586.     VEC_ADD_ABS( temp5v, satdv,     satdv );
  587.     VEC_ADD_ABS( temp6v, satdv,     satdv );
  588.     VEC_ADD_ABS( temp7v, satdv,     satdv );
  590.     satdv = vec_sums( satdv, zero_s32v );
  591.     satdv = vec_splat( satdv, 3 );
  592.     vec_ste( satdv, 0, &i_satd );
  594.     return i_satd / 2;
  595. }
  597. /****************************************************************************
  598.  * x264_pixel_init:
  599.  ****************************************************************************/
  600. void x264_pixel_altivec_init( x264_pixel_function_t *pixf )
  601. {
  602.     pixf->sad[PIXEL_16x16]  = pixel_sad_16x16_altivec;
  603.     pixf->sad[PIXEL_8x16]   = pixel_sad_8x16_altivec;
  604.     pixf->sad[PIXEL_16x8]   = pixel_sad_16x8_altivec;
  605.     pixf->sad[PIXEL_8x8]    = pixel_sad_8x8_altivec;
  607.     pixf->satd[PIXEL_16x16] = pixel_satd_16x16_altivec;
  608.     pixf->satd[PIXEL_8x16]  = pixel_satd_8x16_altivec;
  609.     pixf->satd[PIXEL_16x8]  = pixel_satd_16x8_altivec;
  610.     pixf->satd[PIXEL_8x8]   = pixel_satd_8x8_altivec;
  611.     pixf->satd[PIXEL_8x4]   = pixel_satd_8x4_altivec;
  612.     pixf->satd[PIXEL_4x8]   = pixel_satd_4x8_altivec;
  613.     pixf->satd[PIXEL_4x4]   = pixel_satd_4x4_altivec;
  614. }