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

源码开发语言/平台
当前位置：首页 > 源码/资料 > 多媒体 > 流媒体/Mpeg4/MP4 > 查看源码
dct.c
资源名称：X264CODEC.rar　[点击查看]
上传用户：lctgjx
上传日期：2022-06-04
资源大小：8887k
文件大小：21k
源码类别：

流媒体/Mpeg4/MP4

开发平台：

Visual C++

dct.c：源码内容
  1. /*****************************************************************************
  2.  * dct.c: h264 encoder library
  3.  *****************************************************************************
  4.  * Copyright (C) 2003-2008 x264 project
  5.  *
  6.  * Authors: Guillaume Poirier <gpoirier@mplayerhq.hu>
  7.  *          Eric Petit <eric.petit@lapsus.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., 51 Franklin Street, Fifth Floor, Boston, MA  02111, USA.
  22.  *****************************************************************************/
  24. #include "common/common.h"
  25. #include "ppccommon.h"
  27. #define VEC_DCT(a0,a1,a2,a3,b0,b1,b2,b3) 
  28.     b1 = vec_add( a0, a3 );              
  29.     b3 = vec_add( a1, a2 );              
  30.     b0 = vec_add( b1, b3 );              
  31.     b2 = vec_sub( b1, b3 );              
  32.     a0 = vec_sub( a0, a3 );              
  33.     a1 = vec_sub( a1, a2 );              
  34.     b1 = vec_add( a0, a0 );              
  35.     b1 = vec_add( b1, a1 );              
  36.     b3 = vec_sub( a0, a1 );              
  37.     b3 = vec_sub( b3, a1 )
  39. void x264_sub4x4_dct_altivec( int16_t dct[4][4],
  40.         uint8_t *pix1, uint8_t *pix2 )
  41. {
  42.     PREP_DIFF_8BYTEALIGNED;
  43.     vec_s16_t dct0v, dct1v, dct2v, dct3v;
  44.     vec_s16_t tmp0v, tmp1v, tmp2v, tmp3v;
  46.     vec_u8_t permHighv;
  48.     VEC_DIFF_H_8BYTE_ALIGNED( pix1, FENC_STRIDE, pix2, FDEC_STRIDE, 4, dct0v );
  49.     VEC_DIFF_H_8BYTE_ALIGNED( pix1, FENC_STRIDE, pix2, FDEC_STRIDE, 4, dct1v );
  50.     VEC_DIFF_H_8BYTE_ALIGNED( pix1, FENC_STRIDE, pix2, FDEC_STRIDE, 4, dct2v );
  51.     VEC_DIFF_H_8BYTE_ALIGNED( pix1, FENC_STRIDE, pix2, FDEC_STRIDE, 4, dct3v );
  52.     VEC_DCT( dct0v, dct1v, dct2v, dct3v, tmp0v, tmp1v, tmp2v, tmp3v );
  53.     VEC_TRANSPOSE_4( tmp0v, tmp1v, tmp2v, tmp3v,
  54.                      dct0v, dct1v, dct2v, dct3v );
  55.     permHighv = (vec_u8_t) CV(0x00, 0x01, 0x02, 0x03, 0x04, 0x05, 0x06, 0x07, 0x10, 0x11, 0x12, 0x13, 0x14, 0x15, 0x16, 0x17);
  56.     VEC_DCT( dct0v, dct1v, dct2v, dct3v, tmp0v, tmp1v, tmp2v, tmp3v );
  58.     vec_st(vec_perm(tmp0v, tmp1v, permHighv), 0,  (int16_t*)dct);
  59.     vec_st(vec_perm(tmp2v, tmp3v, permHighv), 16, (int16_t*)dct);
  60. }
  62. void x264_sub8x8_dct_altivec( int16_t dct[4][4][4],
  63.         uint8_t *pix1, uint8_t *pix2 )
  64. {
  65.     PREP_DIFF_8BYTEALIGNED;
  66.     vec_s16_t dct0v, dct1v, dct2v, dct3v, dct4v, dct5v, dct6v, dct7v;
  67.     vec_s16_t tmp0v, tmp1v, tmp2v, tmp3v, tmp4v, tmp5v, tmp6v, tmp7v;
  69.     vec_u8_t permHighv, permLowv;
  71.     VEC_DIFF_H_8BYTE_ALIGNED( pix1, FENC_STRIDE, pix2, FDEC_STRIDE, 8, dct0v );
  72.     VEC_DIFF_H_8BYTE_ALIGNED( pix1, FENC_STRIDE, pix2, FDEC_STRIDE, 8, dct1v );
  73.     VEC_DIFF_H_8BYTE_ALIGNED( pix1, FENC_STRIDE, pix2, FDEC_STRIDE, 8, dct2v );
  74.     VEC_DIFF_H_8BYTE_ALIGNED( pix1, FENC_STRIDE, pix2, FDEC_STRIDE, 8, dct3v );
  75.     VEC_DIFF_H_8BYTE_ALIGNED( pix1, FENC_STRIDE, pix2, FDEC_STRIDE, 8, dct4v );
  76.     VEC_DIFF_H_8BYTE_ALIGNED( pix1, FENC_STRIDE, pix2, FDEC_STRIDE, 8, dct5v );
  77.     VEC_DIFF_H_8BYTE_ALIGNED( pix1, FENC_STRIDE, pix2, FDEC_STRIDE, 8, dct6v );
  78.     VEC_DIFF_H_8BYTE_ALIGNED( pix1, FENC_STRIDE, pix2, FDEC_STRIDE, 8, dct7v );
  79.     VEC_DCT( dct0v, dct1v, dct2v, dct3v, tmp0v, tmp1v, tmp2v, tmp3v );
  80.     VEC_DCT( dct4v, dct5v, dct6v, dct7v, tmp4v, tmp5v, tmp6v, tmp7v );
  81.     VEC_TRANSPOSE_8( tmp0v, tmp1v, tmp2v, tmp3v,
  82.                      tmp4v, tmp5v, tmp6v, tmp7v,
  83.                      dct0v, dct1v, dct2v, dct3v,
  84.                      dct4v, dct5v, dct6v, dct7v );
  86.     permHighv = (vec_u8_t) CV(0x00, 0x01, 0x02, 0x03, 0x04, 0x05, 0x06, 0x07, 0x10, 0x11, 0x12, 0x13, 0x14, 0x15, 0x16, 0x17);
  87.     permLowv  = (vec_u8_t) CV(0x08, 0x09, 0x0A, 0x0B, 0x0C, 0x0D, 0x0E, 0x0F, 0x18, 0x19, 0x1A, 0x1B, 0x1C, 0x1D, 0x1E, 0x1F);
  89.     VEC_DCT( dct0v, dct1v, dct2v, dct3v, tmp0v, tmp1v, tmp2v, tmp3v );
  90.     VEC_DCT( dct4v, dct5v, dct6v, dct7v, tmp4v, tmp5v, tmp6v, tmp7v );
  92.     vec_st(vec_perm(tmp0v, tmp1v, permHighv), 0,   (int16_t*)dct);
  93.     vec_st(vec_perm(tmp2v, tmp3v, permHighv), 16,  (int16_t*)dct);
  94.     vec_st(vec_perm(tmp4v, tmp5v, permHighv), 32,  (int16_t*)dct);
  95.     vec_st(vec_perm(tmp6v, tmp7v, permHighv), 48,  (int16_t*)dct);
  96.     vec_st(vec_perm(tmp0v, tmp1v, permLowv),  64,  (int16_t*)dct);
  97.     vec_st(vec_perm(tmp2v, tmp3v, permLowv),  80,  (int16_t*)dct);
  98.     vec_st(vec_perm(tmp4v, tmp5v, permLowv),  96,  (int16_t*)dct);
  99.     vec_st(vec_perm(tmp6v, tmp7v, permLowv),  112, (int16_t*)dct);
  100. }
  102. void x264_sub16x16_dct_altivec( int16_t dct[16][4][4],
  103.         uint8_t *pix1, uint8_t *pix2 ) 
  104. {
  105.     x264_sub8x8_dct_altivec( &dct[ 0], &pix1[0], &pix2[0] );
  106.     x264_sub8x8_dct_altivec( &dct[ 4], &pix1[8], &pix2[8] );
  107.     x264_sub8x8_dct_altivec( &dct[ 8], &pix1[8*FENC_STRIDE+0], &pix2[8*FDEC_STRIDE+0] );
  108.     x264_sub8x8_dct_altivec( &dct[12], &pix1[8*FENC_STRIDE+8], &pix2[8*FDEC_STRIDE+8] );
  109. }
  111. /***************************************************************************
  112.  * 8x8 transform:
  113.  ***************************************************************************/
  115. /* DCT8_1D unrolled by 8 in Altivec */
  116. #define DCT8_1D_ALTIVEC( dct0v, dct1v, dct2v, dct3v, dct4v, dct5v, dct6v, dct7v ) 
  117. { 
  118.     /* int s07 = SRC(0) + SRC(7);         */ 
  119.     vec_s16_t s07v = vec_add( dct0v, dct7v); 
  120.     /* int s16 = SRC(1) + SRC(6);         */ 
  121.     vec_s16_t s16v = vec_add( dct1v, dct6v); 
  122.     /* int s25 = SRC(2) + SRC(5);         */ 
  123.     vec_s16_t s25v = vec_add( dct2v, dct5v); 
  124.     /* int s34 = SRC(3) + SRC(4);         */ 
  125.     vec_s16_t s34v = vec_add( dct3v, dct4v); 
  127.     /* int a0 = s07 + s34;                */ 
  128.     vec_s16_t a0v = vec_add(s07v, s34v);     
  129.     /* int a1 = s16 + s25;                */ 
  130.     vec_s16_t a1v = vec_add(s16v, s25v);     
  131.     /* int a2 = s07 - s34;                */ 
  132.     vec_s16_t a2v = vec_sub(s07v, s34v);     
  133.     /* int a3 = s16 - s25;                */ 
  134.     vec_s16_t a3v = vec_sub(s16v, s25v);     
  136.     /* int d07 = SRC(0) - SRC(7);         */ 
  137.     vec_s16_t d07v = vec_sub( dct0v, dct7v); 
  138.     /* int d16 = SRC(1) - SRC(6);         */ 
  139.     vec_s16_t d16v = vec_sub( dct1v, dct6v); 
  140.     /* int d25 = SRC(2) - SRC(5);         */ 
  141.     vec_s16_t d25v = vec_sub( dct2v, dct5v); 
  142.     /* int d34 = SRC(3) - SRC(4);         */ 
  143.     vec_s16_t d34v = vec_sub( dct3v, dct4v); 
  145.     /* int a4 = d16 + d25 + (d07 + (d07>>1)); */ 
  146.     vec_s16_t a4v = vec_add( vec_add(d16v, d25v), vec_add(d07v, vec_sra(d07v, onev)) );
  147.     /* int a5 = d07 - d34 - (d25 + (d25>>1)); */ 
  148.     vec_s16_t a5v = vec_sub( vec_sub(d07v, d34v), vec_add(d25v, vec_sra(d25v, onev)) );
  149.     /* int a6 = d07 + d34 - (d16 + (d16>>1)); */ 
  150.     vec_s16_t a6v = vec_sub( vec_add(d07v, d34v), vec_add(d16v, vec_sra(d16v, onev)) );
  151.     /* int a7 = d16 - d25 + (d34 + (d34>>1)); */ 
  152.     vec_s16_t a7v = vec_add( vec_sub(d16v, d25v), vec_add(d34v, vec_sra(d34v, onev)) );
  154.     /* DST(0) =  a0 + a1;                    */ 
  155.     dct0v = vec_add( a0v, a1v );                
  156.     /* DST(1) =  a4 + (a7>>2);               */ 
  157.     dct1v = vec_add( a4v, vec_sra(a7v, twov) ); 
  158.     /* DST(2) =  a2 + (a3>>1);               */ 
  159.     dct2v = vec_add( a2v, vec_sra(a3v, onev) ); 
  160.     /* DST(3) =  a5 + (a6>>2);               */ 
  161.     dct3v = vec_add( a5v, vec_sra(a6v, twov) ); 
  162.     /* DST(4) =  a0 - a1;                    */ 
  163.     dct4v = vec_sub( a0v, a1v );                
  164.     /* DST(5) =  a6 - (a5>>2);               */ 
  165.     dct5v = vec_sub( a6v, vec_sra(a5v, twov) ); 
  166.     /* DST(6) = (a2>>1) - a3 ;               */ 
  167.     dct6v = vec_sub( vec_sra(a2v, onev), a3v ); 
  168.     /* DST(7) = (a4>>2) - a7 ;               */ 
  169.     dct7v = vec_sub( vec_sra(a4v, twov), a7v ); 
  170. }
  173. void x264_sub8x8_dct8_altivec( int16_t dct[8][8], uint8_t *pix1, uint8_t *pix2 )
  174. {
  175.     vec_u16_t onev = vec_splat_u16(1);
  176.     vec_u16_t twov = vec_add( onev, onev );
  178.     PREP_DIFF_8BYTEALIGNED;
  180.     vec_s16_t dct0v, dct1v, dct2v, dct3v,
  181.               dct4v, dct5v, dct6v, dct7v;
  183.     VEC_DIFF_H_8BYTE_ALIGNED( pix1, FENC_STRIDE, pix2, FDEC_STRIDE, 8, dct0v );
  184.     VEC_DIFF_H_8BYTE_ALIGNED( pix1, FENC_STRIDE, pix2, FDEC_STRIDE, 8, dct1v );
  185.     VEC_DIFF_H_8BYTE_ALIGNED( pix1, FENC_STRIDE, pix2, FDEC_STRIDE, 8, dct2v );
  186.     VEC_DIFF_H_8BYTE_ALIGNED( pix1, FENC_STRIDE, pix2, FDEC_STRIDE, 8, dct3v );
  188.     VEC_DIFF_H_8BYTE_ALIGNED( pix1, FENC_STRIDE, pix2, FDEC_STRIDE, 8, dct4v );
  189.     VEC_DIFF_H_8BYTE_ALIGNED( pix1, FENC_STRIDE, pix2, FDEC_STRIDE, 8, dct5v );
  190.     VEC_DIFF_H_8BYTE_ALIGNED( pix1, FENC_STRIDE, pix2, FDEC_STRIDE, 8, dct6v );
  191.     VEC_DIFF_H_8BYTE_ALIGNED( pix1, FENC_STRIDE, pix2, FDEC_STRIDE, 8, dct7v );
  193.     DCT8_1D_ALTIVEC( dct0v, dct1v, dct2v, dct3v,
  194.                      dct4v, dct5v, dct6v, dct7v );
  196.     vec_s16_t dct_tr0v, dct_tr1v, dct_tr2v, dct_tr3v,
  197.         dct_tr4v, dct_tr5v, dct_tr6v, dct_tr7v;
  199.     VEC_TRANSPOSE_8(dct0v, dct1v, dct2v, dct3v,
  200.                     dct4v, dct5v, dct6v, dct7v,
  201.                     dct_tr0v, dct_tr1v, dct_tr2v, dct_tr3v,
  202.                     dct_tr4v, dct_tr5v, dct_tr6v, dct_tr7v );
  204.     DCT8_1D_ALTIVEC( dct_tr0v, dct_tr1v, dct_tr2v, dct_tr3v,
  205.                      dct_tr4v, dct_tr5v, dct_tr6v, dct_tr7v );
  207.     vec_st( dct_tr0v,  0,  (signed short *)dct );
  208.     vec_st( dct_tr1v, 16,  (signed short *)dct );
  209.     vec_st( dct_tr2v, 32,  (signed short *)dct );
  210.     vec_st( dct_tr3v, 48,  (signed short *)dct );
  211.     
  212.     vec_st( dct_tr4v, 64,  (signed short *)dct );
  213.     vec_st( dct_tr5v, 80,  (signed short *)dct );
  214.     vec_st( dct_tr6v, 96,  (signed short *)dct );
  215.     vec_st( dct_tr7v, 112, (signed short *)dct );
  216. }
  218. void x264_sub16x16_dct8_altivec( int16_t dct[4][8][8], uint8_t *pix1, uint8_t *pix2 )
  219. {
  220.     x264_sub8x8_dct8_altivec( dct[0], &pix1[0],               &pix2[0] );
  221.     x264_sub8x8_dct8_altivec( dct[1], &pix1[8],               &pix2[8] );
  222.     x264_sub8x8_dct8_altivec( dct[2], &pix1[8*FENC_STRIDE+0], &pix2[8*FDEC_STRIDE+0] );
  223.     x264_sub8x8_dct8_altivec( dct[3], &pix1[8*FENC_STRIDE+8], &pix2[8*FDEC_STRIDE+8] );
  224. }
  227. /****************************************************************************
  228.  * IDCT transform:
  229.  ****************************************************************************/
  231. #define IDCT_1D_ALTIVEC(s0, s1, s2, s3,  d0, d1, d2, d3) 
  232. {                                                        
  233.     /*        a0  = SRC(0) + SRC(2); */                  
  234.     vec_s16_t a0v = vec_add(s0, s2);                     
  235.     /*        a1  = SRC(0) - SRC(2); */                  
  236.     vec_s16_t a1v = vec_sub(s0, s2);                     
  237.     /*        a2  =           (SRC(1)>>1) - SRC(3); */   
  238.     vec_s16_t a2v = vec_sub(vec_sra(s1, onev), s3);      
  239.     /*        a3  =           (SRC(3)>>1) + SRC(1); */   
  240.     vec_s16_t a3v = vec_add(vec_sra(s3, onev), s1);      
  241.     /* DST(0,    a0 + a3); */                            
  242.     d0 = vec_add(a0v, a3v);                              
  243.     /* DST(1,    a1 + a2); */                            
  244.     d1 = vec_add(a1v, a2v);                              
  245.     /* DST(2,    a1 - a2); */                            
  246.     d2 = vec_sub(a1v, a2v);                              
  247.     /* DST(3,    a0 - a3); */                            
  248.     d3 = vec_sub(a0v, a3v);                              
  249. }
  251. #define VEC_LOAD_U8_ADD_S16_STORE_U8(va)             
  252.     vdst_orig = vec_ld(0, dst);                      
  253.     vdst = vec_perm(vdst_orig, zero_u8v, vdst_mask); 
  254.     vdst_ss = (vec_s16_t)vec_mergeh(zero_u8v, vdst); 
  255.     va = vec_add(va, vdst_ss);                       
  256.     va_u8 = vec_s16_to_u8(va);                       
  257.     va_u32 = vec_splat((vec_u32_t)va_u8, 0);         
  258.     vec_ste(va_u32, element, (uint32_t*)dst);
  260. #define ALTIVEC_STORE4_SUM_CLIP(dest, idctv, perm_ldv)          
  261. {                                                               
  262.     /* unaligned load */                                        
  263.     vec_u8_t lv = vec_ld(0, dest);                              
  264.     vec_u8_t dstv = vec_perm(lv, zero_u8v, (vec_u8_t)perm_ldv); 
  265.     vec_s16_t idct_sh6 = vec_sra(idctv, sixv);                  
  266.     vec_u16_t dst16 = (vec_u16_t)vec_mergeh(zero_u8v, dstv);    
  267.     vec_s16_t idstsum = vec_adds(idct_sh6, (vec_s16_t)dst16);   
  268.     vec_u8_t idstsum8 = vec_s16_to_u8(idstsum);                 
  269.     /* unaligned store */                                       
  270.     vec_u32_t bodyv = vec_splat((vec_u32_t)idstsum8, 0);        
  271.     int element = ((unsigned long)dest & 0xf) >> 2;             
  272.     vec_ste(bodyv, element, (uint32_t *)dest);                  
  273. }
  275. void x264_add4x4_idct_altivec( uint8_t *dst, int16_t dct[4][4] )
  276. {
  277.     vec_u16_t onev = vec_splat_u16(1);
  279.     dct[0][0] += 32; // rounding for the >>6 at the end
  281.     vec_s16_t s0, s1, s2, s3;
  283.     s0 = vec_ld( 0x00, (int16_t*)dct );
  284.     s1 = vec_sld( s0, s0, 8 );
  285.     s2 = vec_ld( 0x10, (int16_t*)dct );
  286.     s3 = vec_sld( s2, s2, 8 );
  288.     vec_s16_t d0, d1, d2, d3;
  289.     IDCT_1D_ALTIVEC( s0, s1, s2, s3, d0, d1, d2, d3 );
  291.     vec_s16_t tr0, tr1, tr2, tr3;
  293.     VEC_TRANSPOSE_4( d0, d1, d2, d3, tr0, tr1, tr2, tr3 );
  295.     vec_s16_t idct0, idct1, idct2, idct3;
  296.     IDCT_1D_ALTIVEC( tr0, tr1, tr2, tr3, idct0, idct1, idct2, idct3 );
  298.     vec_u8_t perm_ldv = vec_lvsl( 0, dst );
  299.     vec_u16_t sixv = vec_splat_u16(6);
  300.     LOAD_ZERO;
  302.     ALTIVEC_STORE4_SUM_CLIP( &dst[0*FDEC_STRIDE], idct0, perm_ldv );
  303.     ALTIVEC_STORE4_SUM_CLIP( &dst[1*FDEC_STRIDE], idct1, perm_ldv );
  304.     ALTIVEC_STORE4_SUM_CLIP( &dst[2*FDEC_STRIDE], idct2, perm_ldv );
  305.     ALTIVEC_STORE4_SUM_CLIP( &dst[3*FDEC_STRIDE], idct3, perm_ldv );
  306. }
  308. void x264_add8x8_idct_altivec( uint8_t *p_dst, int16_t dct[4][4][4] )
  309. {
  310.     x264_add4x4_idct_altivec( &p_dst[0],               dct[0] );
  311.     x264_add4x4_idct_altivec( &p_dst[4],               dct[1] );
  312.     x264_add4x4_idct_altivec( &p_dst[4*FDEC_STRIDE+0], dct[2] );
  313.     x264_add4x4_idct_altivec( &p_dst[4*FDEC_STRIDE+4], dct[3] );
  314. }
  316. void x264_add16x16_idct_altivec( uint8_t *p_dst, int16_t dct[16][4][4] )
  317. {
  318.     x264_add8x8_idct_altivec( &p_dst[0],               &dct[0] );
  319.     x264_add8x8_idct_altivec( &p_dst[8],               &dct[4] );
  320.     x264_add8x8_idct_altivec( &p_dst[8*FDEC_STRIDE+0], &dct[8] );
  321.     x264_add8x8_idct_altivec( &p_dst[8*FDEC_STRIDE+8], &dct[12] );
  322. }
  324. #define IDCT8_1D_ALTIVEC(s0, s1, s2, s3, s4, s5, s6, s7,  d0, d1, d2, d3, d4, d5, d6, d7)
  325. {
  326.     /*        a0  = SRC(0) + SRC(4); */ 
  327.     vec_s16_t a0v = vec_add(s0, s4);    
  328.     /*        a2  = SRC(0) - SRC(4); */ 
  329.     vec_s16_t a2v = vec_sub(s0, s4);    
  330.     /*        a4  =           (SRC(2)>>1) - SRC(6); */ 
  331.     vec_s16_t a4v = vec_sub(vec_sra(s2, onev), s6);    
  332.     /*        a6  =           (SRC(6)>>1) + SRC(2); */ 
  333.     vec_s16_t a6v = vec_add(vec_sra(s6, onev), s2);    
  334.     /*        b0  =         a0 + a6; */ 
  335.     vec_s16_t b0v = vec_add(a0v, a6v);  
  336.     /*        b2  =         a2 + a4; */ 
  337.     vec_s16_t b2v = vec_add(a2v, a4v);  
  338.     /*        b4  =         a2 - a4; */ 
  339.     vec_s16_t b4v = vec_sub(a2v, a4v);  
  340.     /*        b6  =         a0 - a6; */ 
  341.     vec_s16_t b6v = vec_sub(a0v, a6v);  
  342.     /* a1 =  SRC(5) - SRC(3) - SRC(7) - (SRC(7)>>1); */ 
  343.     /*        a1 =             (SRC(5)-SRC(3)) -  (SRC(7)  +  (SRC(7)>>1)); */ 
  344.     vec_s16_t a1v = vec_sub( vec_sub(s5, s3), vec_add(s7, vec_sra(s7, onev)) );
  345.     /* a3 =  SRC(7) + SRC(1) - SRC(3) - (SRC(3)>>1); */ 
  346.     /*        a3 =             (SRC(7)+SRC(1)) -  (SRC(3)  +  (SRC(3)>>1)); */ 
  347.     vec_s16_t a3v = vec_sub( vec_add(s7, s1), vec_add(s3, vec_sra(s3, onev)) );
  348.     /* a5 =  SRC(7) - SRC(1) + SRC(5) + (SRC(5)>>1); */ 
  349.     /*        a5 =             (SRC(7)-SRC(1)) +   SRC(5) +   (SRC(5)>>1); */  
  350.     vec_s16_t a5v = vec_add( vec_sub(s7, s1), vec_add(s5, vec_sra(s5, onev)) );
  351.     /*        a7 =                SRC(5)+SRC(3) +  SRC(1) +   (SRC(1)>>1); */  
  352.     vec_s16_t a7v = vec_add( vec_add(s5, s3), vec_add(s1, vec_sra(s1, onev)) );
  353.     /*        b1 =                  (a7>>2)  +  a1; */  
  354.     vec_s16_t b1v = vec_add( vec_sra(a7v, twov), a1v);  
  355.     /*        b3 =          a3 +        (a5>>2); */     
  356.     vec_s16_t b3v = vec_add(a3v, vec_sra(a5v, twov));   
  357.     /*        b5 =                  (a3>>2)  -   a5; */ 
  358.     vec_s16_t b5v = vec_sub( vec_sra(a3v, twov), a5v);  
  359.     /*        b7 =           a7 -        (a1>>2); */    
  360.     vec_s16_t b7v = vec_sub( a7v, vec_sra(a1v, twov));  
  361.     /* DST(0,    b0 + b7); */ 
  362.     d0 = vec_add(b0v, b7v); 
  363.     /* DST(1,    b2 + b5); */ 
  364.     d1 = vec_add(b2v, b5v); 
  365.     /* DST(2,    b4 + b3); */ 
  366.     d2 = vec_add(b4v, b3v); 
  367.     /* DST(3,    b6 + b1); */ 
  368.     d3 = vec_add(b6v, b1v); 
  369.     /* DST(4,    b6 - b1); */ 
  370.     d4 = vec_sub(b6v, b1v); 
  371.     /* DST(5,    b4 - b3); */ 
  372.     d5 = vec_sub(b4v, b3v); 
  373.     /* DST(6,    b2 - b5); */ 
  374.     d6 = vec_sub(b2v, b5v); 
  375.     /* DST(7,    b0 - b7); */ 
  376.     d7 = vec_sub(b0v, b7v); 
  377. }
  379. #define ALTIVEC_STORE_SUM_CLIP(dest, idctv, perm_ldv, perm_stv, sel)
  380. {
  381.     /* unaligned load */                                       
  382.     vec_u8_t hv = vec_ld( 0, dest );                           
  383.     vec_u8_t lv = vec_ld( 7, dest );                           
  384.     vec_u8_t dstv   = vec_perm( hv, lv, (vec_u8_t)perm_ldv );  
  385.     vec_s16_t idct_sh6 = vec_sra(idctv, sixv);                 
  386.     vec_u16_t dst16 = (vec_u16_t)vec_mergeh(zero_u8v, dstv);   
  387.     vec_s16_t idstsum = vec_adds(idct_sh6, (vec_s16_t)dst16);  
  388.     vec_u8_t idstsum8 = vec_packsu(zero_s16v, idstsum);        
  389.     /* unaligned store */                                      
  390.     vec_u8_t bodyv  = vec_perm( idstsum8, idstsum8, perm_stv );
  391.     vec_u8_t edgelv = vec_perm( sel, zero_u8v, perm_stv );     
  392.     lv    = vec_sel( lv, bodyv, edgelv );                      
  393.     vec_st( lv, 7, dest );                                     
  394.     hv    = vec_ld( 0, dest );                                 
  395.     vec_u8_t edgehv = vec_perm( zero_u8v, sel, perm_stv );     
  396.     hv    = vec_sel( hv, bodyv, edgehv );                      
  397.     vec_st( hv, 0, dest );                                     
  398. }
  400. void x264_add8x8_idct8_altivec( uint8_t *dst, int16_t dct[8][8] )
  401. {
  402.     vec_u16_t onev = vec_splat_u16(1);
  403.     vec_u16_t twov = vec_splat_u16(2);
  405.     dct[0][0] += 32; // rounding for the >>6 at the end
  407.     vec_s16_t s0, s1, s2, s3, s4, s5, s6, s7;
  409.     s0 = vec_ld(0x00, (int16_t*)dct);
  410.     s1 = vec_ld(0x10, (int16_t*)dct);
  411.     s2 = vec_ld(0x20, (int16_t*)dct);
  412.     s3 = vec_ld(0x30, (int16_t*)dct);
  413.     s4 = vec_ld(0x40, (int16_t*)dct);
  414.     s5 = vec_ld(0x50, (int16_t*)dct);
  415.     s6 = vec_ld(0x60, (int16_t*)dct);
  416.     s7 = vec_ld(0x70, (int16_t*)dct);
  418.     vec_s16_t d0, d1, d2, d3, d4, d5, d6, d7;
  419.     IDCT8_1D_ALTIVEC(s0, s1, s2, s3, s4, s5, s6, s7,  d0, d1, d2, d3, d4, d5, d6, d7);
  421.     vec_s16_t tr0, tr1, tr2, tr3, tr4, tr5, tr6, tr7;
  423.     VEC_TRANSPOSE_8( d0,  d1,  d2,  d3,  d4,  d5,  d6, d7,
  424.                     tr0, tr1, tr2, tr3, tr4, tr5, tr6, tr7);
  426.     vec_s16_t idct0, idct1, idct2, idct3, idct4, idct5, idct6, idct7;
  427.     IDCT8_1D_ALTIVEC(tr0,     tr1,   tr2,   tr3,   tr4,   tr5,   tr6,   tr7,
  428.                      idct0, idct1, idct2, idct3, idct4, idct5, idct6, idct7);
  430.     vec_u8_t perm_ldv = vec_lvsl(0, dst);
  431.     vec_u8_t perm_stv = vec_lvsr(8, dst);
  432.     vec_u16_t sixv = vec_splat_u16(6);
  433.     const vec_u8_t sel = (vec_u8_t) CV(0,0,0,0,0,0,0,0,-1,-1,-1,-1,-1,-1,-1,-1);
  434.     LOAD_ZERO;
  436.     ALTIVEC_STORE_SUM_CLIP(&dst[0*FDEC_STRIDE], idct0, perm_ldv, perm_stv, sel);
  437.     ALTIVEC_STORE_SUM_CLIP(&dst[1*FDEC_STRIDE], idct1, perm_ldv, perm_stv, sel);
  438.     ALTIVEC_STORE_SUM_CLIP(&dst[2*FDEC_STRIDE], idct2, perm_ldv, perm_stv, sel);
  439.     ALTIVEC_STORE_SUM_CLIP(&dst[3*FDEC_STRIDE], idct3, perm_ldv, perm_stv, sel);
  440.     ALTIVEC_STORE_SUM_CLIP(&dst[4*FDEC_STRIDE], idct4, perm_ldv, perm_stv, sel);
  441.     ALTIVEC_STORE_SUM_CLIP(&dst[5*FDEC_STRIDE], idct5, perm_ldv, perm_stv, sel);
  442.     ALTIVEC_STORE_SUM_CLIP(&dst[6*FDEC_STRIDE], idct6, perm_ldv, perm_stv, sel);
  443.     ALTIVEC_STORE_SUM_CLIP(&dst[7*FDEC_STRIDE], idct7, perm_ldv, perm_stv, sel);
  444. }
  446. void x264_add16x16_idct8_altivec( uint8_t *dst, int16_t dct[4][8][8] )
  447. {
  448.     x264_add8x8_idct8_altivec( &dst[0],               dct[0] );
  449.     x264_add8x8_idct8_altivec( &dst[8],               dct[1] );
  450.     x264_add8x8_idct8_altivec( &dst[8*FDEC_STRIDE+0], dct[2] );
  451.     x264_add8x8_idct8_altivec( &dst[8*FDEC_STRIDE+8], dct[3] );
  452. }
  454. void x264_zigzag_scan_4x4_frame_altivec( int16_t level[16], int16_t dct[4][4] )
  455. {
  456.     vec_s16_t dct0v, dct1v;
  457.     vec_s16_t tmp0v, tmp1v;
  459.     dct0v = vec_ld(0x00, (int16_t*)dct);
  460.     dct1v = vec_ld(0x10, (int16_t*)dct);
  462.     const vec_u8_t sel0 = (vec_u8_t) CV(0,1,8,9,2,3,4,5,10,11,16,17,24,25,18,19);
  463.     const vec_u8_t sel1 = (vec_u8_t) CV(12,13,6,7,14,15,20,21,26,27,28,29,22,23,30,31);
  465.     tmp0v = vec_perm( dct0v, dct1v, sel0 );
  466.     tmp1v = vec_perm( dct0v, dct1v, sel1 );
  468.     vec_st( tmp0v, 0x00, level );
  469.     vec_st( tmp1v, 0x10, level );
  470. }
  472. void x264_zigzag_scan_4x4_field_altivec( int16_t level[16], int16_t dct[4][4] )
  473. {
  474.     vec_s16_t dct0v, dct1v;
  475.     vec_s16_t tmp0v, tmp1v;
  477.     dct0v = vec_ld(0x00, (int16_t*)dct);
  478.     dct1v = vec_ld(0x10, (int16_t*)dct);
  480.     const vec_u8_t sel0 = (vec_u8_t) CV(0,1,2,3,8,9,4,5,6,7,10,11,12,13,14,15);
  482.     tmp0v = vec_perm( dct0v, dct1v, sel0 );
  483.     tmp1v = dct1v;
  485.     vec_st( tmp0v, 0x00, level );
  486.     vec_st( tmp1v, 0x10, level );
  487. }