dsputil.c
上传用户:chinavct
上传日期:2022-06-20
资源大小:330k
文件大小:174k
源码类别:

流媒体/Mpeg4/MP4

开发平台:

Visual C++

  1.     }
  2. }
  3. static void OPNAME ## h264_qpel8_v_lowpass(uint8_t *dst, uint8_t *src, int dstStride, int srcStride){
  4.     const int w=8;
  5.     uint8_t *cm = ff_cropTbl + MAX_NEG_CROP;
  6.     int i;
  7.     for(i=0; i<w; i++)
  8.     {
  9.         const int srcB= src[-2*srcStride];
  10.         const int srcA= src[-1*srcStride];
  11.         const int src0= src[0 *srcStride];
  12.         const int src1= src[1 *srcStride];
  13.         const int src2= src[2 *srcStride];
  14.         const int src3= src[3 *srcStride];
  15.         const int src4= src[4 *srcStride];
  16.         const int src5= src[5 *srcStride];
  17.         const int src6= src[6 *srcStride];
  18.         const int src7= src[7 *srcStride];
  19.         const int src8= src[8 *srcStride];
  20.         const int src9= src[9 *srcStride];
  21.         const int src10=src[10*srcStride];
  22.         OP(dst[0*dstStride], (src0+src1)*20 - (srcA+src2)*5 + (srcB+src3));
  23.         OP(dst[1*dstStride], (src1+src2)*20 - (src0+src3)*5 + (srcA+src4));
  24.         OP(dst[2*dstStride], (src2+src3)*20 - (src1+src4)*5 + (src0+src5));
  25.         OP(dst[3*dstStride], (src3+src4)*20 - (src2+src5)*5 + (src1+src6));
  26.         OP(dst[4*dstStride], (src4+src5)*20 - (src3+src6)*5 + (src2+src7));
  27.         OP(dst[5*dstStride], (src5+src6)*20 - (src4+src7)*5 + (src3+src8));
  28.         OP(dst[6*dstStride], (src6+src7)*20 - (src5+src8)*5 + (src4+src9));
  29.         OP(dst[7*dstStride], (src7+src8)*20 - (src6+src9)*5 + (src5+src10));
  30.         dst++;
  31.         src++;
  32.     }
  33. }
  34. static void OPNAME ## h264_qpel8_hv_lowpass(uint8_t *dst, int16_t *tmp, uint8_t *src, int dstStride, int tmpStride, int srcStride){
  35.     const int h=8;
  36.     const int w=8;
  37.     uint8_t *cm = ff_cropTbl + MAX_NEG_CROP;
  38.     int i;
  39.     src -= 2*srcStride;
  40.     for(i=0; i<h+5; i++)
  41.     {
  42.         tmp[0]= (src[0]+src[1])*20 - (src[-1]+src[2])*5 + (src[-2]+src[3 ]);
  43.         tmp[1]= (src[1]+src[2])*20 - (src[0 ]+src[3])*5 + (src[-1]+src[4 ]);
  44.         tmp[2]= (src[2]+src[3])*20 - (src[1 ]+src[4])*5 + (src[0 ]+src[5 ]);
  45.         tmp[3]= (src[3]+src[4])*20 - (src[2 ]+src[5])*5 + (src[1 ]+src[6 ]);
  46.         tmp[4]= (src[4]+src[5])*20 - (src[3 ]+src[6])*5 + (src[2 ]+src[7 ]);
  47.         tmp[5]= (src[5]+src[6])*20 - (src[4 ]+src[7])*5 + (src[3 ]+src[8 ]);
  48.         tmp[6]= (src[6]+src[7])*20 - (src[5 ]+src[8])*5 + (src[4 ]+src[9 ]);
  49.         tmp[7]= (src[7]+src[8])*20 - (src[6 ]+src[9])*5 + (src[5 ]+src[10]);
  50.         tmp+=tmpStride;
  51.         src+=srcStride;
  52.     }
  53.     tmp -= tmpStride*(h+5-2);
  54.     for(i=0; i<w; i++)
  55.     {
  56.         const int tmpB= tmp[-2*tmpStride];
  57.         const int tmpA= tmp[-1*tmpStride];
  58.         const int tmp0= tmp[0 *tmpStride];
  59.         const int tmp1= tmp[1 *tmpStride];
  60.         const int tmp2= tmp[2 *tmpStride];
  61.         const int tmp3= tmp[3 *tmpStride];
  62.         const int tmp4= tmp[4 *tmpStride];
  63.         const int tmp5= tmp[5 *tmpStride];
  64.         const int tmp6= tmp[6 *tmpStride];
  65.         const int tmp7= tmp[7 *tmpStride];
  66.         const int tmp8= tmp[8 *tmpStride];
  67.         const int tmp9= tmp[9 *tmpStride];
  68.         const int tmp10=tmp[10*tmpStride];
  69.         OP2(dst[0*dstStride], (tmp0+tmp1)*20 - (tmpA+tmp2)*5 + (tmpB+tmp3));
  70.         OP2(dst[1*dstStride], (tmp1+tmp2)*20 - (tmp0+tmp3)*5 + (tmpA+tmp4));
  71.         OP2(dst[2*dstStride], (tmp2+tmp3)*20 - (tmp1+tmp4)*5 + (tmp0+tmp5));
  72.         OP2(dst[3*dstStride], (tmp3+tmp4)*20 - (tmp2+tmp5)*5 + (tmp1+tmp6));
  73.         OP2(dst[4*dstStride], (tmp4+tmp5)*20 - (tmp3+tmp6)*5 + (tmp2+tmp7));
  74.         OP2(dst[5*dstStride], (tmp5+tmp6)*20 - (tmp4+tmp7)*5 + (tmp3+tmp8));
  75.         OP2(dst[6*dstStride], (tmp6+tmp7)*20 - (tmp5+tmp8)*5 + (tmp4+tmp9));
  76.         OP2(dst[7*dstStride], (tmp7+tmp8)*20 - (tmp6+tmp9)*5 + (tmp5+tmp10));
  77.         dst++;
  78.         tmp++;
  79.     }
  80. }
  81. static void OPNAME ## h264_qpel16_v_lowpass(uint8_t *dst, uint8_t *src, int dstStride, int srcStride){
  82.     OPNAME ## h264_qpel8_v_lowpass(dst  , src  , dstStride, srcStride);
  83.     OPNAME ## h264_qpel8_v_lowpass(dst+8, src+8, dstStride, srcStride);
  84.     src += 8*srcStride;
  85.     dst += 8*dstStride;
  86.     OPNAME ## h264_qpel8_v_lowpass(dst  , src  , dstStride, srcStride);
  87.     OPNAME ## h264_qpel8_v_lowpass(dst+8, src+8, dstStride, srcStride);
  88. }
  89. static void OPNAME ## h264_qpel16_h_lowpass(uint8_t *dst, uint8_t *src, int dstStride, int srcStride){
  90.     OPNAME ## h264_qpel8_h_lowpass(dst  , src  , dstStride, srcStride);
  91.     OPNAME ## h264_qpel8_h_lowpass(dst+8, src+8, dstStride, srcStride);
  92.     src += 8*srcStride;
  93.     dst += 8*dstStride;
  94.     OPNAME ## h264_qpel8_h_lowpass(dst  , src  , dstStride, srcStride);
  95.     OPNAME ## h264_qpel8_h_lowpass(dst+8, src+8, dstStride, srcStride);
  96. }
  97. static void OPNAME ## h264_qpel16_hv_lowpass(uint8_t *dst, int16_t *tmp, uint8_t *src, int dstStride, int tmpStride, int srcStride){
  98.     OPNAME ## h264_qpel8_hv_lowpass(dst  , tmp  , src  , dstStride, tmpStride, srcStride);
  99.     OPNAME ## h264_qpel8_hv_lowpass(dst+8, tmp+8, src+8, dstStride, tmpStride, srcStride);
  100.     src += 8*srcStride;
  101.     dst += 8*dstStride;
  102.     OPNAME ## h264_qpel8_hv_lowpass(dst  , tmp  , src  , dstStride, tmpStride, srcStride);
  103.     OPNAME ## h264_qpel8_hv_lowpass(dst+8, tmp+8, src+8, dstStride, tmpStride, srcStride);
  104. }
  105. #define H264_MC(OPNAME, SIZE) 
  106. static void OPNAME ## h264_qpel ## SIZE ## _mc00_c (uint8_t *dst, uint8_t *src, int stride){
  107.     OPNAME ## pixels ## SIZE ## _c(dst, src, stride, SIZE);
  108. }
  109. static void OPNAME ## h264_qpel ## SIZE ## _mc10_c(uint8_t *dst, uint8_t *src, int stride){
  110.     uint8_t half[SIZE*SIZE];
  111.     put_h264_qpel ## SIZE ## _h_lowpass(half, src, SIZE, stride);
  112.     OPNAME ## pixels ## SIZE ## _l2(dst, src, half, stride, stride, SIZE, SIZE);
  113. }
  114. static void OPNAME ## h264_qpel ## SIZE ## _mc20_c(uint8_t *dst, uint8_t *src, int stride){
  115.     OPNAME ## h264_qpel ## SIZE ## _h_lowpass(dst, src, stride, stride);
  116. }
  117. static void OPNAME ## h264_qpel ## SIZE ## _mc30_c(uint8_t *dst, uint8_t *src, int stride){
  118.     uint8_t half[SIZE*SIZE];
  119.     put_h264_qpel ## SIZE ## _h_lowpass(half, src, SIZE, stride);
  120.     OPNAME ## pixels ## SIZE ## _l2(dst, src+1, half, stride, stride, SIZE, SIZE);
  121. }
  122. static void OPNAME ## h264_qpel ## SIZE ## _mc01_c(uint8_t *dst, uint8_t *src, int stride){
  123.     uint8_t full[SIZE*(SIZE+5)];
  124.     uint8_t * const full_mid= full + SIZE*2;
  125.     uint8_t half[SIZE*SIZE];
  126.     copy_block ## SIZE (full, src - stride*2, SIZE,  stride, SIZE + 5);
  127.     put_h264_qpel ## SIZE ## _v_lowpass(half, full_mid, SIZE, SIZE);
  128.     OPNAME ## pixels ## SIZE ## _l2(dst, full_mid, half, stride, SIZE, SIZE, SIZE);
  129. }
  130. static void OPNAME ## h264_qpel ## SIZE ## _mc02_c(uint8_t *dst, uint8_t *src, int stride){
  131.     uint8_t full[SIZE*(SIZE+5)];
  132.     uint8_t * const full_mid= full + SIZE*2;
  133.     copy_block ## SIZE (full, src - stride*2, SIZE,  stride, SIZE + 5);
  134.     OPNAME ## h264_qpel ## SIZE ## _v_lowpass(dst, full_mid, stride, SIZE);
  135. }
  136. static void OPNAME ## h264_qpel ## SIZE ## _mc03_c(uint8_t *dst, uint8_t *src, int stride){
  137.     uint8_t full[SIZE*(SIZE+5)];
  138.     uint8_t * const full_mid= full + SIZE*2;
  139.     uint8_t half[SIZE*SIZE];
  140.     copy_block ## SIZE (full, src - stride*2, SIZE,  stride, SIZE + 5);
  141.     put_h264_qpel ## SIZE ## _v_lowpass(half, full_mid, SIZE, SIZE);
  142.     OPNAME ## pixels ## SIZE ## _l2(dst, full_mid+SIZE, half, stride, SIZE, SIZE, SIZE);
  143. }
  144. static void OPNAME ## h264_qpel ## SIZE ## _mc11_c(uint8_t *dst, uint8_t *src, int stride){
  145.     uint8_t full[SIZE*(SIZE+5)];
  146.     uint8_t * const full_mid= full + SIZE*2;
  147.     uint8_t halfH[SIZE*SIZE];
  148.     uint8_t halfV[SIZE*SIZE];
  149.     put_h264_qpel ## SIZE ## _h_lowpass(halfH, src, SIZE, stride);
  150.     copy_block ## SIZE (full, src - stride*2, SIZE,  stride, SIZE + 5);
  151.     put_h264_qpel ## SIZE ## _v_lowpass(halfV, full_mid, SIZE, SIZE);
  152.     OPNAME ## pixels ## SIZE ## _l2(dst, halfH, halfV, stride, SIZE, SIZE, SIZE);
  153. }
  154. static void OPNAME ## h264_qpel ## SIZE ## _mc31_c(uint8_t *dst, uint8_t *src, int stride){
  155.     uint8_t full[SIZE*(SIZE+5)];
  156.     uint8_t * const full_mid= full + SIZE*2;
  157.     uint8_t halfH[SIZE*SIZE];
  158.     uint8_t halfV[SIZE*SIZE];
  159.     put_h264_qpel ## SIZE ## _h_lowpass(halfH, src, SIZE, stride);
  160.     copy_block ## SIZE (full, src - stride*2 + 1, SIZE,  stride, SIZE + 5);
  161.     put_h264_qpel ## SIZE ## _v_lowpass(halfV, full_mid, SIZE, SIZE);
  162.     OPNAME ## pixels ## SIZE ## _l2(dst, halfH, halfV, stride, SIZE, SIZE, SIZE);
  163. }
  164. static void OPNAME ## h264_qpel ## SIZE ## _mc13_c(uint8_t *dst, uint8_t *src, int stride){
  165.     uint8_t full[SIZE*(SIZE+5)];
  166.     uint8_t * const full_mid= full + SIZE*2;
  167.     uint8_t halfH[SIZE*SIZE];
  168.     uint8_t halfV[SIZE*SIZE];
  169.     put_h264_qpel ## SIZE ## _h_lowpass(halfH, src + stride, SIZE, stride);
  170.     copy_block ## SIZE (full, src - stride*2, SIZE,  stride, SIZE + 5);
  171.     put_h264_qpel ## SIZE ## _v_lowpass(halfV, full_mid, SIZE, SIZE);
  172.     OPNAME ## pixels ## SIZE ## _l2(dst, halfH, halfV, stride, SIZE, SIZE, SIZE);
  173. }
  174. static void OPNAME ## h264_qpel ## SIZE ## _mc33_c(uint8_t *dst, uint8_t *src, int stride){
  175.     uint8_t full[SIZE*(SIZE+5)];
  176.     uint8_t * const full_mid= full + SIZE*2;
  177.     uint8_t halfH[SIZE*SIZE];
  178.     uint8_t halfV[SIZE*SIZE];
  179.     put_h264_qpel ## SIZE ## _h_lowpass(halfH, src + stride, SIZE, stride);
  180.     copy_block ## SIZE (full, src - stride*2 + 1, SIZE,  stride, SIZE + 5);
  181.     put_h264_qpel ## SIZE ## _v_lowpass(halfV, full_mid, SIZE, SIZE);
  182.     OPNAME ## pixels ## SIZE ## _l2(dst, halfH, halfV, stride, SIZE, SIZE, SIZE);
  183. }
  184. static void OPNAME ## h264_qpel ## SIZE ## _mc22_c(uint8_t *dst, uint8_t *src, int stride){
  185.     int16_t tmp[SIZE*(SIZE+5)];
  186.     OPNAME ## h264_qpel ## SIZE ## _hv_lowpass(dst, tmp, src, stride, SIZE, stride);
  187. }
  188. static void OPNAME ## h264_qpel ## SIZE ## _mc21_c(uint8_t *dst, uint8_t *src, int stride){
  189.     int16_t tmp[SIZE*(SIZE+5)];
  190.     uint8_t halfH[SIZE*SIZE];
  191.     uint8_t halfHV[SIZE*SIZE];
  192.     put_h264_qpel ## SIZE ## _h_lowpass(halfH, src, SIZE, stride);
  193.     put_h264_qpel ## SIZE ## _hv_lowpass(halfHV, tmp, src, SIZE, SIZE, stride);
  194.     OPNAME ## pixels ## SIZE ## _l2(dst, halfH, halfHV, stride, SIZE, SIZE, SIZE);
  195. }
  196. static void OPNAME ## h264_qpel ## SIZE ## _mc23_c(uint8_t *dst, uint8_t *src, int stride){
  197.     int16_t tmp[SIZE*(SIZE+5)];
  198.     uint8_t halfH[SIZE*SIZE];
  199.     uint8_t halfHV[SIZE*SIZE];
  200.     put_h264_qpel ## SIZE ## _h_lowpass(halfH, src + stride, SIZE, stride);
  201.     put_h264_qpel ## SIZE ## _hv_lowpass(halfHV, tmp, src, SIZE, SIZE, stride);
  202.     OPNAME ## pixels ## SIZE ## _l2(dst, halfH, halfHV, stride, SIZE, SIZE, SIZE);
  203. }
  204. static void OPNAME ## h264_qpel ## SIZE ## _mc12_c(uint8_t *dst, uint8_t *src, int stride){
  205.     uint8_t full[SIZE*(SIZE+5)];
  206.     uint8_t * const full_mid= full + SIZE*2;
  207.     int16_t tmp[SIZE*(SIZE+5)];
  208.     uint8_t halfV[SIZE*SIZE];
  209.     uint8_t halfHV[SIZE*SIZE];
  210.     copy_block ## SIZE (full, src - stride*2, SIZE,  stride, SIZE + 5);
  211.     put_h264_qpel ## SIZE ## _v_lowpass(halfV, full_mid, SIZE, SIZE);
  212.     put_h264_qpel ## SIZE ## _hv_lowpass(halfHV, tmp, src, SIZE, SIZE, stride);
  213.     OPNAME ## pixels ## SIZE ## _l2(dst, halfV, halfHV, stride, SIZE, SIZE, SIZE);
  214. }
  215. static void OPNAME ## h264_qpel ## SIZE ## _mc32_c(uint8_t *dst, uint8_t *src, int stride){
  216.     uint8_t full[SIZE*(SIZE+5)];
  217.     uint8_t * const full_mid= full + SIZE*2;
  218.     int16_t tmp[SIZE*(SIZE+5)];
  219.     uint8_t halfV[SIZE*SIZE];
  220.     uint8_t halfHV[SIZE*SIZE];
  221.     copy_block ## SIZE (full, src - stride*2 + 1, SIZE,  stride, SIZE + 5);
  222.     put_h264_qpel ## SIZE ## _v_lowpass(halfV, full_mid, SIZE, SIZE);
  223.     put_h264_qpel ## SIZE ## _hv_lowpass(halfHV, tmp, src, SIZE, SIZE, stride);
  224.     OPNAME ## pixels ## SIZE ## _l2(dst, halfV, halfHV, stride, SIZE, SIZE, SIZE);
  225. }
  226. #define op_avg(a, b)   a = (((a)+cm[((b) + 16)>>5]+1)>>1)
  227. #define op2_avg(a, b)  a = (((a)+cm[((b) + 512)>>10]+1)>>1)
  228. #define op_put(a, b)   a = cm[((b) + 16)>>5]
  229. #define op2_put(a, b)  a = cm[((b) + 512)>>10]
  230. H264_LOWPASS(avg_       , op_avg, op2_avg)
  231. H264_LOWPASS(put_       , op_put, op2_put)
  232. static av_unused void put_h264_qpel2_h_lowpass(uint8_t *dst, uint8_t *src, int dstStride, int srcStride)
  233. {
  234.     const int h=2;
  235.     uint8_t *cm = ff_cropTbl + MAX_NEG_CROP;
  236.     int i;
  237.     for(i=0; i<h; i++)
  238.     {
  239.         /*OP(dst[0], (src[0]+src[1])*20 - (src[-1]+src[2])*5 + (src[-2]+src[3]));*/
  240. dst[0] = (((dst[0])+cm[(((src[0]+src[1])*20 - (src[-1]+src[2])*5 + (src[-2]+src[3])) + 16)>>5]+1)>>1);
  241.         /*OP(dst[1], (src[1]+src[2])*20 - (src[0 ]+src[3])*5 + (src[-1]+src[4]));*/
  242. dst[1] = (((dst[1])+cm[(((src[1]+src[2])*20 - (src[0 ]+src[3])*5 + (src[-1]+src[4])) + 16)>>5]+1)>>1);
  243.         dst+=dstStride;
  244.         src+=srcStride;
  245.     }
  246. }
  247. static av_unused void put_h264_qpel2_v_lowpass(uint8_t *dst, uint8_t *src, int dstStride, int srcStride)
  248. {
  249.     const int w=2;
  250.     uint8_t *cm = ff_cropTbl + MAX_NEG_CROP;
  251.     int i;
  252.     for(i=0; i<w; i++)
  253.     {
  254.         const int srcB= src[-2*srcStride];
  255.         const int srcA= src[-1*srcStride];
  256.         const int src0= src[0 *srcStride];
  257.         const int src1= src[1 *srcStride];
  258.         const int src2= src[2 *srcStride];
  259.         const int src3= src[3 *srcStride];
  260.         const int src4= src[4 *srcStride];
  261.         /*OP(dst[0*dstStride], (src0+src1)*20 - (srcA+src2)*5 + (srcB+src3));*/
  262. dst[0*dstStride] = (((dst[0*dstStride])+cm[(((src0+src1)*20 - (srcA+src2)*5 + (srcB+src3)) + 16)>>5]+1)>>1);
  263.         /*OP(dst[1*dstStride], (src1+src2)*20 - (src0+src3)*5 + (srcA+src4));*/
  264. dst[1*dstStride] = (((dst[1*dstStride])+cm[(((src1+src2)*20 - (src0+src3)*5 + (srcA+src4)) + 16)>>5]+1)>>1);
  265.         dst++;
  266.         src++;
  267.     }
  268. }
  269. static av_unused void put_h264_qpel2_hv_lowpass(uint8_t *dst, int16_t *tmp, uint8_t *src, int dstStride, int tmpStride, int srcStride)
  270. {
  271.     const int h=2;
  272.     const int w=2;
  273.     uint8_t *cm = ff_cropTbl + MAX_NEG_CROP;
  274.     int i;
  275.     src -= 2*srcStride;
  276.     for(i=0; i<h+5; i++)
  277.     {
  278.         tmp[0]= (src[0]+src[1])*20 - (src[-1]+src[2])*5 + (src[-2]+src[3]);
  279.         tmp[1]= (src[1]+src[2])*20 - (src[0 ]+src[3])*5 + (src[-1]+src[4]);
  280.         tmp+=tmpStride;
  281.         src+=srcStride;
  282.     }
  283.     tmp -= tmpStride*(h+5-2);
  284.     for(i=0; i<w; i++)
  285.     {
  286.         const int tmpB= tmp[-2*tmpStride];
  287.         const int tmpA= tmp[-1*tmpStride];
  288.         const int tmp0= tmp[0 *tmpStride];
  289.         const int tmp1= tmp[1 *tmpStride];
  290.         const int tmp2= tmp[2 *tmpStride];
  291.         const int tmp3= tmp[3 *tmpStride];
  292.         const int tmp4= tmp[4 *tmpStride];
  293.         /*OP2(dst[0*dstStride], (tmp0+tmp1)*20 - (tmpA+tmp2)*5 + (tmpB+tmp3));*/
  294. dst[0*dstStride] = (((dst[0*dstStride])+cm[(((tmp0+tmp1)*20 - (tmpA+tmp2)*5 + (tmpB+tmp3)) + 512)>>10]+1)>>1);
  295.         /*OP2(dst[1*dstStride], (tmp1+tmp2)*20 - (tmp0+tmp3)*5 + (tmpA+tmp4));*/
  296. dst[1*dstStride] = (((dst[1*dstStride])+cm[(((tmp1+tmp2)*20 - (tmp0+tmp3)*5 + (tmpA+tmp4)) + 512)>>10]+1)>>1);
  297.         dst++;
  298.         tmp++;
  299.     }
  300. }
  301. H264_MC(put_, 2)
  302. H264_MC(put_, 4)
  303. H264_MC(put_, 8)
  304. H264_MC(put_, 16)
  305. H264_MC(avg_, 4)
  306. H264_MC(avg_, 8)
  307. H264_MC(avg_, 16)
  308. #undef op_avg
  309. #undef op_put
  310. #undef op2_avg
  311. #undef op2_put
  312. #endif
  313. #define op_scale1(x)  block[x] = av_clip_uint8( (block[x]*weight + offset) >> log2_denom )
  314. #define op_scale2(x)  dst[x] = av_clip_uint8( (src[x]*weights + dst[x]*weightd + offset) >> (log2_denom+1))
  315. #define H264_WEIGHT(W,H) 
  316. static void weight_h264_pixels ## W ## x ## H ## _c(uint8_t *block, int stride, int log2_denom, int weight, int offset){ 
  317.     int y; 
  318.     offset <<= log2_denom; 
  319.     if(log2_denom) offset += 1<<(log2_denom-1); 
  320.     for(y=0; y<H; y++, block += stride){ 
  321.         op_scale1(0); 
  322.         op_scale1(1); 
  323.         if(W==2) continue; 
  324.         op_scale1(2); 
  325.         op_scale1(3); 
  326.         if(W==4) continue; 
  327.         op_scale1(4); 
  328.         op_scale1(5); 
  329.         op_scale1(6); 
  330.         op_scale1(7); 
  331.         if(W==8) continue; 
  332.         op_scale1(8); 
  333.         op_scale1(9); 
  334.         op_scale1(10); 
  335.         op_scale1(11); 
  336.         op_scale1(12); 
  337.         op_scale1(13); 
  338.         op_scale1(14); 
  339.         op_scale1(15); 
  340.     } 
  341. static void biweight_h264_pixels ## W ## x ## H ## _c(uint8_t *dst, uint8_t *src, int stride, int log2_denom, int weightd, int weights, int offset){ 
  342.     int y; 
  343.     offset = ((offset + 1) | 1) << log2_denom; 
  344.     for(y=0; y<H; y++, dst += stride, src += stride){ 
  345.         op_scale2(0); 
  346.         op_scale2(1); 
  347.         if(W==2) continue; 
  348.         op_scale2(2); 
  349.         op_scale2(3); 
  350.         if(W==4) continue; 
  351.         op_scale2(4); 
  352.         op_scale2(5); 
  353.         op_scale2(6); 
  354.         op_scale2(7); 
  355.         if(W==8) continue; 
  356.         op_scale2(8); 
  357.         op_scale2(9); 
  358.         op_scale2(10); 
  359.         op_scale2(11); 
  360.         op_scale2(12); 
  361.         op_scale2(13); 
  362.         op_scale2(14); 
  363.         op_scale2(15); 
  364.     } 
  365. }
  366. H264_WEIGHT(16,16)
  367. H264_WEIGHT(16,8)
  368. H264_WEIGHT(8,16)
  369. H264_WEIGHT(8,8)
  370. H264_WEIGHT(8,4)
  371. H264_WEIGHT(4,8)
  372. H264_WEIGHT(4,4)
  373. H264_WEIGHT(4,2)
  374. H264_WEIGHT(2,4)
  375. H264_WEIGHT(2,2)
  376. #undef op_scale1
  377. #undef op_scale2
  378. #undef H264_WEIGHT
  379. static void wmv2_mspel8_h_lowpass(uint8_t *dst, uint8_t *src, int dstStride, int srcStride, int h){
  380.     uint8_t *cm = ff_cropTbl + MAX_NEG_CROP;
  381.     int i;
  382.     for(i=0; i<h; i++){
  383.         dst[0]= cm[(9*(src[0] + src[1]) - (src[-1] + src[2]) + 8)>>4];
  384.         dst[1]= cm[(9*(src[1] + src[2]) - (src[ 0] + src[3]) + 8)>>4];
  385.         dst[2]= cm[(9*(src[2] + src[3]) - (src[ 1] + src[4]) + 8)>>4];
  386.         dst[3]= cm[(9*(src[3] + src[4]) - (src[ 2] + src[5]) + 8)>>4];
  387.         dst[4]= cm[(9*(src[4] + src[5]) - (src[ 3] + src[6]) + 8)>>4];
  388.         dst[5]= cm[(9*(src[5] + src[6]) - (src[ 4] + src[7]) + 8)>>4];
  389.         dst[6]= cm[(9*(src[6] + src[7]) - (src[ 5] + src[8]) + 8)>>4];
  390.         dst[7]= cm[(9*(src[7] + src[8]) - (src[ 6] + src[9]) + 8)>>4];
  391.         dst+=dstStride;
  392.         src+=srcStride;
  393.     }
  394. }
  395. #ifdef CONFIG_CAVS_DECODER
  396. /* AVS specific */
  397. void ff_cavsdsp_init(DSPContext* c, AVCodecContext *avctx);
  398. void ff_put_cavs_qpel8_mc00_c(uint8_t *dst, uint8_t *src, int stride) {
  399.     put_pixels8_c(dst, src, stride, 8);
  400. }
  401. void ff_avg_cavs_qpel8_mc00_c(uint8_t *dst, uint8_t *src, int stride) {
  402.     avg_pixels8_c(dst, src, stride, 8);
  403. }
  404. void ff_put_cavs_qpel16_mc00_c(uint8_t *dst, uint8_t *src, int stride) {
  405.     put_pixels16_c(dst, src, stride, 16);
  406. }
  407. void ff_avg_cavs_qpel16_mc00_c(uint8_t *dst, uint8_t *src, int stride) {
  408.     avg_pixels16_c(dst, src, stride, 16);
  409. }
  410. #endif /* CONFIG_CAVS_DECODER */
  411. #if defined(CONFIG_VC1_DECODER) || defined(CONFIG_WMV3_DECODER)
  412. /* VC-1 specific */
  413. void ff_vc1dsp_init(DSPContext* c, AVCodecContext *avctx);
  414. void ff_put_vc1_mspel_mc00_c(uint8_t *dst, uint8_t *src, int stride, int rnd) {
  415.     put_pixels8_c(dst, src, stride, 8);
  416. }
  417. #endif /* CONFIG_VC1_DECODER||CONFIG_WMV3_DECODER */
  418. void ff_intrax8dsp_init(DSPContext* c, AVCodecContext *avctx);
  419. /* H264 specific */
  420. void ff_h264dspenc_init(DSPContext* c, AVCodecContext *avctx);
  421. #if defined(CONFIG_RV30_DECODER)
  422. void ff_rv30dsp_init(DSPContext* c, AVCodecContext *avctx);
  423. #endif /* CONFIG_RV30_DECODER */
  424. #if defined(CONFIG_RV40_DECODER)
  425. static void put_rv40_qpel16_mc33_c(uint8_t *dst, uint8_t *src, int stride){
  426.     put_pixels16_xy2_c(dst, src, stride, 16);
  427. }
  428. static void avg_rv40_qpel16_mc33_c(uint8_t *dst, uint8_t *src, int stride){
  429.     avg_pixels16_xy2_c(dst, src, stride, 16);
  430. }
  431. static void put_rv40_qpel8_mc33_c(uint8_t *dst, uint8_t *src, int stride){
  432.     put_pixels8_xy2_c(dst, src, stride, 8);
  433. }
  434. static void avg_rv40_qpel8_mc33_c(uint8_t *dst, uint8_t *src, int stride){
  435.     avg_pixels8_xy2_c(dst, src, stride, 8);
  436. }
  437. void ff_rv40dsp_init(DSPContext* c, AVCodecContext *avctx);
  438. #endif /* CONFIG_RV40_DECODER */
  439. static void wmv2_mspel8_v_lowpass(uint8_t *dst, uint8_t *src, int dstStride, int srcStride, int w){
  440.     uint8_t *cm = ff_cropTbl + MAX_NEG_CROP;
  441.     int i;
  442.     for(i=0; i<w; i++){
  443.         const int src_1= src[ -srcStride];
  444.         const int src0 = src[0          ];
  445.         const int src1 = src[  srcStride];
  446.         const int src2 = src[2*srcStride];
  447.         const int src3 = src[3*srcStride];
  448.         const int src4 = src[4*srcStride];
  449.         const int src5 = src[5*srcStride];
  450.         const int src6 = src[6*srcStride];
  451.         const int src7 = src[7*srcStride];
  452.         const int src8 = src[8*srcStride];
  453.         const int src9 = src[9*srcStride];
  454.         dst[0*dstStride]= cm[(9*(src0 + src1) - (src_1 + src2) + 8)>>4];
  455.         dst[1*dstStride]= cm[(9*(src1 + src2) - (src0  + src3) + 8)>>4];
  456.         dst[2*dstStride]= cm[(9*(src2 + src3) - (src1  + src4) + 8)>>4];
  457.         dst[3*dstStride]= cm[(9*(src3 + src4) - (src2  + src5) + 8)>>4];
  458.         dst[4*dstStride]= cm[(9*(src4 + src5) - (src3  + src6) + 8)>>4];
  459.         dst[5*dstStride]= cm[(9*(src5 + src6) - (src4  + src7) + 8)>>4];
  460.         dst[6*dstStride]= cm[(9*(src6 + src7) - (src5  + src8) + 8)>>4];
  461.         dst[7*dstStride]= cm[(9*(src7 + src8) - (src6  + src9) + 8)>>4];
  462.         src++;
  463.         dst++;
  464.     }
  465. }
  466. static void put_mspel8_mc00_c (uint8_t *dst, uint8_t *src, int stride){
  467.     put_pixels8_c(dst, src, stride, 8);
  468. }
  469. static void put_mspel8_mc10_c(uint8_t *dst, uint8_t *src, int stride){
  470.     uint8_t half[64];
  471.     wmv2_mspel8_h_lowpass(half, src, 8, stride, 8);
  472.     put_pixels8_l2(dst, src, half, stride, stride, 8, 8);
  473. }
  474. static void put_mspel8_mc20_c(uint8_t *dst, uint8_t *src, int stride){
  475.     wmv2_mspel8_h_lowpass(dst, src, stride, stride, 8);
  476. }
  477. static void put_mspel8_mc30_c(uint8_t *dst, uint8_t *src, int stride){
  478.     uint8_t half[64];
  479.     wmv2_mspel8_h_lowpass(half, src, 8, stride, 8);
  480.     put_pixels8_l2(dst, src+1, half, stride, stride, 8, 8);
  481. }
  482. static void put_mspel8_mc02_c(uint8_t *dst, uint8_t *src, int stride){
  483.     wmv2_mspel8_v_lowpass(dst, src, stride, stride, 8);
  484. }
  485. static void put_mspel8_mc12_c(uint8_t *dst, uint8_t *src, int stride){
  486.     uint8_t halfH[88];
  487.     uint8_t halfV[64];
  488.     uint8_t halfHV[64];
  489.     wmv2_mspel8_h_lowpass(halfH, src-stride, 8, stride, 11);
  490.     wmv2_mspel8_v_lowpass(halfV, src, 8, stride, 8);
  491.     wmv2_mspel8_v_lowpass(halfHV, halfH+8, 8, 8, 8);
  492.     put_pixels8_l2(dst, halfV, halfHV, stride, 8, 8, 8);
  493. }
  494. static void put_mspel8_mc32_c(uint8_t *dst, uint8_t *src, int stride){
  495.     uint8_t halfH[88];
  496.     uint8_t halfV[64];
  497.     uint8_t halfHV[64];
  498.     wmv2_mspel8_h_lowpass(halfH, src-stride, 8, stride, 11);
  499.     wmv2_mspel8_v_lowpass(halfV, src+1, 8, stride, 8);
  500.     wmv2_mspel8_v_lowpass(halfHV, halfH+8, 8, 8, 8);
  501.     put_pixels8_l2(dst, halfV, halfHV, stride, 8, 8, 8);
  502. }
  503. static void put_mspel8_mc22_c(uint8_t *dst, uint8_t *src, int stride){
  504.     uint8_t halfH[88];
  505.     wmv2_mspel8_h_lowpass(halfH, src-stride, 8, stride, 11);
  506.     wmv2_mspel8_v_lowpass(dst, halfH+8, stride, 8, 8);
  507. }
  508. /*
  509. static void h263_v_loop_filter_c(uint8_t *src, int stride, int qscale){
  510.     if(ENABLE_ANY_H263) {
  511.     int x;
  512.     const int strength= ff_h263_loop_filter_strength[qscale];
  513.     for(x=0; x<8; x++){
  514.         int d1, d2, ad1;
  515.         int p0= src[x-2*stride];
  516.         int p1= src[x-1*stride];
  517.         int p2= src[x+0*stride];
  518.         int p3= src[x+1*stride];
  519.         int d = (p0 - p3 + 4*(p2 - p1)) / 8;
  520.         if     (d<-2*strength) d1= 0;
  521.         else if(d<-  strength) d1=-2*strength - d;
  522.         else if(d<   strength) d1= d;
  523.         else if(d< 2*strength) d1= 2*strength - d;
  524.         else                   d1= 0;
  525.         p1 += d1;
  526.         p2 -= d1;
  527.         if(p1&256) p1= ~(p1>>31);
  528.         if(p2&256) p2= ~(p2>>31);
  529.         src[x-1*stride] = p1;
  530.         src[x+0*stride] = p2;
  531.         ad1= FFABS(d1)>>1;
  532.         d2= av_clip((p0-p3)/4, -ad1, ad1);
  533.         src[x-2*stride] = p0 - d2;
  534.         src[x+  stride] = p3 + d2;
  535.     }
  536.     }
  537. }
  538. static void h263_h_loop_filter_c(uint8_t *src, int stride, int qscale){
  539.     if(ENABLE_ANY_H263) {
  540.     int y;
  541.     const int strength= ff_h263_loop_filter_strength[qscale];
  542.     for(y=0; y<8; y++){
  543.         int d1, d2, ad1;
  544.         int p0= src[y*stride-2];
  545.         int p1= src[y*stride-1];
  546.         int p2= src[y*stride+0];
  547.         int p3= src[y*stride+1];
  548.         int d = (p0 - p3 + 4*(p2 - p1)) / 8;
  549.         if     (d<-2*strength) d1= 0;
  550.         else if(d<-  strength) d1=-2*strength - d;
  551.         else if(d<   strength) d1= d;
  552.         else if(d< 2*strength) d1= 2*strength - d;
  553.         else                   d1= 0;
  554.         p1 += d1;
  555.         p2 -= d1;
  556.         if(p1&256) p1= ~(p1>>31);
  557.         if(p2&256) p2= ~(p2>>31);
  558.         src[y*stride-1] = p1;
  559.         src[y*stride+0] = p2;
  560.         ad1= FFABS(d1)>>1;
  561.         d2= av_clip((p0-p3)/4, -ad1, ad1);
  562.         src[y*stride-2] = p0 - d2;
  563.         src[y*stride+1] = p3 + d2;
  564.     }
  565.     }
  566. }
  567. */
  568. static void h261_loop_filter_c(uint8_t *src, int stride){
  569.     int x,y,xy,yz;
  570.     int temp[64];
  571.     for(x=0; x<8; x++){
  572.         temp[x      ] = 4*src[x           ];
  573.         temp[x + 7*8] = 4*src[x + 7*stride];
  574.     }
  575.     for(y=1; y<7; y++){
  576.         for(x=0; x<8; x++){
  577.             xy = y * stride + x;
  578.             yz = y * 8 + x;
  579.             temp[yz] = src[xy - stride] + 2*src[xy] + src[xy + stride];
  580.         }
  581.     }
  582.     for(y=0; y<8; y++){
  583.         src[  y*stride] = (temp[  y*8] + 2)>>2;
  584.         src[7+y*stride] = (temp[7+y*8] + 2)>>2;
  585.         for(x=1; x<7; x++){
  586.             xy = y * stride + x;
  587.             yz = y * 8 + x;
  588.             src[xy] = (temp[yz-1] + 2*temp[yz] + temp[yz+1] + 8)>>4;
  589.         }
  590.     }
  591. }
  592. static inline void h264_loop_filter_luma_c(uint8_t *pix, int xstride, int ystride, int alpha, int beta, int8_t *tc0)
  593. {
  594.     int i, d;
  595.     for( i = 0; i < 4; i++ ) {
  596.         if( tc0[i] < 0 ) {
  597.             pix += 4*ystride;
  598.             continue;
  599.         }
  600.         for( d = 0; d < 4; d++ ) {
  601.             const int p0 = pix[-1*xstride];
  602.             const int p1 = pix[-2*xstride];
  603.             const int p2 = pix[-3*xstride];
  604.             const int q0 = pix[0];
  605.             const int q1 = pix[1*xstride];
  606.             const int q2 = pix[2*xstride];
  607.             if( FFABS( p0 - q0 ) < alpha &&
  608.                 FFABS( p1 - p0 ) < beta &&
  609.                 FFABS( q1 - q0 ) < beta ) {
  610.                 int tc = tc0[i];
  611.                 int i_delta;
  612.                 if( FFABS( p2 - p0 ) < beta ) {
  613.                     pix[-2*xstride] = p1 + av_clip( (( p2 + ( ( p0 + q0 + 1 ) >> 1 ) ) >> 1) - p1, -tc0[i], tc0[i] );
  614.                     tc++;
  615.                 }
  616.                 if( FFABS( q2 - q0 ) < beta ) {
  617.                     pix[   xstride] = q1 + av_clip( (( q2 + ( ( p0 + q0 + 1 ) >> 1 ) ) >> 1) - q1, -tc0[i], tc0[i] );
  618.                     tc++;
  619.                 }
  620.                 i_delta = av_clip( (((q0 - p0 ) << 2) + (p1 - q1) + 4) >> 3, -tc, tc );
  621.                 pix[-xstride] = av_clip_uint8( p0 + i_delta );    /* p0' */
  622.                 pix[0]        = av_clip_uint8( q0 - i_delta );    /* q0' */
  623.             }
  624.             pix += ystride;
  625.         }
  626.     }
  627. }
  628. static void h264_v_loop_filter_luma_c(uint8_t *pix, int stride, int alpha, int beta, int8_t *tc0)
  629. {
  630.     h264_loop_filter_luma_c(pix, stride, 1, alpha, beta, tc0);
  631. }
  632. static void h264_h_loop_filter_luma_c(uint8_t *pix, int stride, int alpha, int beta, int8_t *tc0)
  633. {
  634.     h264_loop_filter_luma_c(pix, 1, stride, alpha, beta, tc0);
  635. }
  636. static inline void h264_loop_filter_luma_intra_c(uint8_t *pix, int xstride, int ystride, int alpha, int beta)
  637. {
  638.     int d;
  639.     for( d = 0; d < 16; d++ ) {
  640.         const int p2 = pix[-3*xstride];
  641.         const int p1 = pix[-2*xstride];
  642.         const int p0 = pix[-1*xstride];
  643.         const int q0 = pix[ 0*xstride];
  644.         const int q1 = pix[ 1*xstride];
  645.         const int q2 = pix[ 2*xstride];
  646.         if( FFABS( p0 - q0 ) < alpha &&
  647.             FFABS( p1 - p0 ) < beta &&
  648.             FFABS( q1 - q0 ) < beta ) {
  649.             if(FFABS( p0 - q0 ) < (( alpha >> 2 ) + 2 )){
  650.                 if( FFABS( p2 - p0 ) < beta)
  651.                 {
  652.                     const int p3 = pix[-4*xstride];
  653.                     /* p0', p1', p2' */
  654.                     pix[-1*xstride] = ( p2 + 2*p1 + 2*p0 + 2*q0 + q1 + 4 ) >> 3;
  655.                     pix[-2*xstride] = ( p2 + p1 + p0 + q0 + 2 ) >> 2;
  656.                     pix[-3*xstride] = ( 2*p3 + 3*p2 + p1 + p0 + q0 + 4 ) >> 3;
  657.                 } else {
  658.                     /* p0' */
  659.                     pix[-1*xstride] = ( 2*p1 + p0 + q1 + 2 ) >> 2;
  660.                 }
  661.                 if( FFABS( q2 - q0 ) < beta)
  662.                 {
  663.                     const int q3 = pix[3*xstride];
  664.                     /* q0', q1', q2' */
  665.                     pix[0*xstride] = ( p1 + 2*p0 + 2*q0 + 2*q1 + q2 + 4 ) >> 3;
  666.                     pix[1*xstride] = ( p0 + q0 + q1 + q2 + 2 ) >> 2;
  667.                     pix[2*xstride] = ( 2*q3 + 3*q2 + q1 + q0 + p0 + 4 ) >> 3;
  668.                 } else {
  669.                     /* q0' */
  670.                     pix[0*xstride] = ( 2*q1 + q0 + p1 + 2 ) >> 2;
  671.                 }
  672.             }else{
  673.                 /* p0', q0' */
  674.                 pix[-1*xstride] = ( 2*p1 + p0 + q1 + 2 ) >> 2;
  675.                 pix[ 0*xstride] = ( 2*q1 + q0 + p1 + 2 ) >> 2;
  676.             }
  677.         }
  678.         pix += ystride;
  679.     }
  680. }
  681. static void h264_v_loop_filter_luma_intra_c(uint8_t *pix, int stride, int alpha, int beta)
  682. {
  683.     h264_loop_filter_luma_intra_c(pix, stride, 1, alpha, beta);
  684. }
  685. static void h264_h_loop_filter_luma_intra_c(uint8_t *pix, int stride, int alpha, int beta)
  686. {
  687.     h264_loop_filter_luma_intra_c(pix, 1, stride, alpha, beta);
  688. }
  689. static inline void h264_loop_filter_chroma_c(uint8_t *pix, int xstride, int ystride, int alpha, int beta, int8_t *tc0)
  690. {
  691.     int i, d;
  692.     for( i = 0; i < 4; i++ ) {
  693.         const int tc = tc0[i];
  694.         if( tc <= 0 ) {
  695.             pix += 2*ystride;
  696.             continue;
  697.         }
  698.         for( d = 0; d < 2; d++ ) {
  699.             const int p0 = pix[-1*xstride];
  700.             const int p1 = pix[-2*xstride];
  701.             const int q0 = pix[0];
  702.             const int q1 = pix[1*xstride];
  703.             if( FFABS( p0 - q0 ) < alpha &&
  704.                 FFABS( p1 - p0 ) < beta &&
  705.                 FFABS( q1 - q0 ) < beta ) {
  706.                 int delta = av_clip( (((q0 - p0 ) << 2) + (p1 - q1) + 4) >> 3, -tc, tc );
  707.                 pix[-xstride] = av_clip_uint8( p0 + delta );    /* p0' */
  708.                 pix[0]        = av_clip_uint8( q0 - delta );    /* q0' */
  709.             }
  710.             pix += ystride;
  711.         }
  712.     }
  713. }
  714. static void h264_v_loop_filter_chroma_c(uint8_t *pix, int stride, int alpha, int beta, int8_t *tc0)
  715. {
  716.     h264_loop_filter_chroma_c(pix, stride, 1, alpha, beta, tc0);
  717. }
  718. static void h264_h_loop_filter_chroma_c(uint8_t *pix, int stride, int alpha, int beta, int8_t *tc0)
  719. {
  720.     h264_loop_filter_chroma_c(pix, 1, stride, alpha, beta, tc0);
  721. }
  722. static inline void h264_loop_filter_chroma_intra_c(uint8_t *pix, int xstride, int ystride, int alpha, int beta)
  723. {
  724.     int d;
  725.     for( d = 0; d < 8; d++ ) {
  726.         const int p0 = pix[-1*xstride];
  727.         const int p1 = pix[-2*xstride];
  728.         const int q0 = pix[0];
  729.         const int q1 = pix[1*xstride];
  730.         if( FFABS( p0 - q0 ) < alpha &&
  731.             FFABS( p1 - p0 ) < beta &&
  732.             FFABS( q1 - q0 ) < beta ) {
  733.             pix[-xstride] = ( 2*p1 + p0 + q1 + 2 ) >> 2;   /* p0' */
  734.             pix[0]        = ( 2*q1 + q0 + p1 + 2 ) >> 2;   /* q0' */
  735.         }
  736.         pix += ystride;
  737.     }
  738. }
  739. static void h264_v_loop_filter_chroma_intra_c(uint8_t *pix, int stride, int alpha, int beta)
  740. {
  741.     h264_loop_filter_chroma_intra_c(pix, stride, 1, alpha, beta);
  742. }
  743. static void h264_h_loop_filter_chroma_intra_c(uint8_t *pix, int stride, int alpha, int beta)
  744. {
  745.     h264_loop_filter_chroma_intra_c(pix, 1, stride, alpha, beta);
  746. }
  747. static inline int pix_abs16_c(void *v, uint8_t *pix1, uint8_t *pix2, int line_size, int h)
  748. {
  749.     int s, i;
  750.     s = 0;
  751.     for(i=0;i<h;i++) {
  752.         s += abs(pix1[0] - pix2[0]);
  753.         s += abs(pix1[1] - pix2[1]);
  754.         s += abs(pix1[2] - pix2[2]);
  755.         s += abs(pix1[3] - pix2[3]);
  756.         s += abs(pix1[4] - pix2[4]);
  757.         s += abs(pix1[5] - pix2[5]);
  758.         s += abs(pix1[6] - pix2[6]);
  759.         s += abs(pix1[7] - pix2[7]);
  760.         s += abs(pix1[8] - pix2[8]);
  761.         s += abs(pix1[9] - pix2[9]);
  762.         s += abs(pix1[10] - pix2[10]);
  763.         s += abs(pix1[11] - pix2[11]);
  764.         s += abs(pix1[12] - pix2[12]);
  765.         s += abs(pix1[13] - pix2[13]);
  766.         s += abs(pix1[14] - pix2[14]);
  767.         s += abs(pix1[15] - pix2[15]);
  768.         pix1 += line_size;
  769.         pix2 += line_size;
  770.     }
  771.     return s;
  772. }
  773. static int pix_abs16_x2_c(void *v, uint8_t *pix1, uint8_t *pix2, int line_size, int h)
  774. {
  775.     int s, i;
  776.     s = 0;
  777.     for(i=0;i<h;i++) {
  778.         s += abs(pix1[0] - avg2(pix2[0], pix2[1]));
  779.         s += abs(pix1[1] - avg2(pix2[1], pix2[2]));
  780.         s += abs(pix1[2] - avg2(pix2[2], pix2[3]));
  781.         s += abs(pix1[3] - avg2(pix2[3], pix2[4]));
  782.         s += abs(pix1[4] - avg2(pix2[4], pix2[5]));
  783.         s += abs(pix1[5] - avg2(pix2[5], pix2[6]));
  784.         s += abs(pix1[6] - avg2(pix2[6], pix2[7]));
  785.         s += abs(pix1[7] - avg2(pix2[7], pix2[8]));
  786.         s += abs(pix1[8] - avg2(pix2[8], pix2[9]));
  787.         s += abs(pix1[9] - avg2(pix2[9], pix2[10]));
  788.         s += abs(pix1[10] - avg2(pix2[10], pix2[11]));
  789.         s += abs(pix1[11] - avg2(pix2[11], pix2[12]));
  790.         s += abs(pix1[12] - avg2(pix2[12], pix2[13]));
  791.         s += abs(pix1[13] - avg2(pix2[13], pix2[14]));
  792.         s += abs(pix1[14] - avg2(pix2[14], pix2[15]));
  793.         s += abs(pix1[15] - avg2(pix2[15], pix2[16]));
  794.         pix1 += line_size;
  795.         pix2 += line_size;
  796.     }
  797.     return s;
  798. }
  799. static int pix_abs16_y2_c(void *v, uint8_t *pix1, uint8_t *pix2, int line_size, int h)
  800. {
  801.     int s, i;
  802.     uint8_t *pix3 = pix2 + line_size;
  803.     s = 0;
  804.     for(i=0;i<h;i++) {
  805.         s += abs(pix1[0] - avg2(pix2[0], pix3[0]));
  806.         s += abs(pix1[1] - avg2(pix2[1], pix3[1]));
  807.         s += abs(pix1[2] - avg2(pix2[2], pix3[2]));
  808.         s += abs(pix1[3] - avg2(pix2[3], pix3[3]));
  809.         s += abs(pix1[4] - avg2(pix2[4], pix3[4]));
  810.         s += abs(pix1[5] - avg2(pix2[5], pix3[5]));
  811.         s += abs(pix1[6] - avg2(pix2[6], pix3[6]));
  812.         s += abs(pix1[7] - avg2(pix2[7], pix3[7]));
  813.         s += abs(pix1[8] - avg2(pix2[8], pix3[8]));
  814.         s += abs(pix1[9] - avg2(pix2[9], pix3[9]));
  815.         s += abs(pix1[10] - avg2(pix2[10], pix3[10]));
  816.         s += abs(pix1[11] - avg2(pix2[11], pix3[11]));
  817.         s += abs(pix1[12] - avg2(pix2[12], pix3[12]));
  818.         s += abs(pix1[13] - avg2(pix2[13], pix3[13]));
  819.         s += abs(pix1[14] - avg2(pix2[14], pix3[14]));
  820.         s += abs(pix1[15] - avg2(pix2[15], pix3[15]));
  821.         pix1 += line_size;
  822.         pix2 += line_size;
  823.         pix3 += line_size;
  824.     }
  825.     return s;
  826. }
  827. static int pix_abs16_xy2_c(void *v, uint8_t *pix1, uint8_t *pix2, int line_size, int h)
  828. {
  829.     int s, i;
  830.     uint8_t *pix3 = pix2 + line_size;
  831.     s = 0;
  832.     for(i=0;i<h;i++) {
  833.         s += abs(pix1[0] - avg4(pix2[0], pix2[1], pix3[0], pix3[1]));
  834.         s += abs(pix1[1] - avg4(pix2[1], pix2[2], pix3[1], pix3[2]));
  835.         s += abs(pix1[2] - avg4(pix2[2], pix2[3], pix3[2], pix3[3]));
  836.         s += abs(pix1[3] - avg4(pix2[3], pix2[4], pix3[3], pix3[4]));
  837.         s += abs(pix1[4] - avg4(pix2[4], pix2[5], pix3[4], pix3[5]));
  838.         s += abs(pix1[5] - avg4(pix2[5], pix2[6], pix3[5], pix3[6]));
  839.         s += abs(pix1[6] - avg4(pix2[6], pix2[7], pix3[6], pix3[7]));
  840.         s += abs(pix1[7] - avg4(pix2[7], pix2[8], pix3[7], pix3[8]));
  841.         s += abs(pix1[8] - avg4(pix2[8], pix2[9], pix3[8], pix3[9]));
  842.         s += abs(pix1[9] - avg4(pix2[9], pix2[10], pix3[9], pix3[10]));
  843.         s += abs(pix1[10] - avg4(pix2[10], pix2[11], pix3[10], pix3[11]));
  844.         s += abs(pix1[11] - avg4(pix2[11], pix2[12], pix3[11], pix3[12]));
  845.         s += abs(pix1[12] - avg4(pix2[12], pix2[13], pix3[12], pix3[13]));
  846.         s += abs(pix1[13] - avg4(pix2[13], pix2[14], pix3[13], pix3[14]));
  847.         s += abs(pix1[14] - avg4(pix2[14], pix2[15], pix3[14], pix3[15]));
  848.         s += abs(pix1[15] - avg4(pix2[15], pix2[16], pix3[15], pix3[16]));
  849.         pix1 += line_size;
  850.         pix2 += line_size;
  851.         pix3 += line_size;
  852.     }
  853.     return s;
  854. }
  855. static inline int pix_abs8_c(void *v, uint8_t *pix1, uint8_t *pix2, int line_size, int h)
  856. {
  857.     int s, i;
  858.     s = 0;
  859.     for(i=0;i<h;i++) {
  860.         s += abs(pix1[0] - pix2[0]);
  861.         s += abs(pix1[1] - pix2[1]);
  862.         s += abs(pix1[2] - pix2[2]);
  863.         s += abs(pix1[3] - pix2[3]);
  864.         s += abs(pix1[4] - pix2[4]);
  865.         s += abs(pix1[5] - pix2[5]);
  866.         s += abs(pix1[6] - pix2[6]);
  867.         s += abs(pix1[7] - pix2[7]);
  868.         pix1 += line_size;
  869.         pix2 += line_size;
  870.     }
  871.     return s;
  872. }
  873. static int pix_abs8_x2_c(void *v, uint8_t *pix1, uint8_t *pix2, int line_size, int h)
  874. {
  875.     int s, i;
  876.     s = 0;
  877.     for(i=0;i<h;i++) {
  878.         s += abs(pix1[0] - avg2(pix2[0], pix2[1]));
  879.         s += abs(pix1[1] - avg2(pix2[1], pix2[2]));
  880.         s += abs(pix1[2] - avg2(pix2[2], pix2[3]));
  881.         s += abs(pix1[3] - avg2(pix2[3], pix2[4]));
  882.         s += abs(pix1[4] - avg2(pix2[4], pix2[5]));
  883.         s += abs(pix1[5] - avg2(pix2[5], pix2[6]));
  884.         s += abs(pix1[6] - avg2(pix2[6], pix2[7]));
  885.         s += abs(pix1[7] - avg2(pix2[7], pix2[8]));
  886.         pix1 += line_size;
  887.         pix2 += line_size;
  888.     }
  889.     return s;
  890. }
  891. static int pix_abs8_y2_c(void *v, uint8_t *pix1, uint8_t *pix2, int line_size, int h)
  892. {
  893.     int s, i;
  894.     uint8_t *pix3 = pix2 + line_size;
  895.     s = 0;
  896.     for(i=0;i<h;i++) {
  897.         s += abs(pix1[0] - avg2(pix2[0], pix3[0]));
  898.         s += abs(pix1[1] - avg2(pix2[1], pix3[1]));
  899.         s += abs(pix1[2] - avg2(pix2[2], pix3[2]));
  900.         s += abs(pix1[3] - avg2(pix2[3], pix3[3]));
  901.         s += abs(pix1[4] - avg2(pix2[4], pix3[4]));
  902.         s += abs(pix1[5] - avg2(pix2[5], pix3[5]));
  903.         s += abs(pix1[6] - avg2(pix2[6], pix3[6]));
  904.         s += abs(pix1[7] - avg2(pix2[7], pix3[7]));
  905.         pix1 += line_size;
  906.         pix2 += line_size;
  907.         pix3 += line_size;
  908.     }
  909.     return s;
  910. }
  911. static int pix_abs8_xy2_c(void *v, uint8_t *pix1, uint8_t *pix2, int line_size, int h)
  912. {
  913.     int s, i;
  914.     uint8_t *pix3 = pix2 + line_size;
  915.     s = 0;
  916.     for(i=0;i<h;i++) {
  917.         s += abs(pix1[0] - avg4(pix2[0], pix2[1], pix3[0], pix3[1]));
  918.         s += abs(pix1[1] - avg4(pix2[1], pix2[2], pix3[1], pix3[2]));
  919.         s += abs(pix1[2] - avg4(pix2[2], pix2[3], pix3[2], pix3[3]));
  920.         s += abs(pix1[3] - avg4(pix2[3], pix2[4], pix3[3], pix3[4]));
  921.         s += abs(pix1[4] - avg4(pix2[4], pix2[5], pix3[4], pix3[5]));
  922.         s += abs(pix1[5] - avg4(pix2[5], pix2[6], pix3[5], pix3[6]));
  923.         s += abs(pix1[6] - avg4(pix2[6], pix2[7], pix3[6], pix3[7]));
  924.         s += abs(pix1[7] - avg4(pix2[7], pix2[8], pix3[7], pix3[8]));
  925.         pix1 += line_size;
  926.         pix2 += line_size;
  927.         pix3 += line_size;
  928.     }
  929.     return s;
  930. }
  931. static int nsse16_c(void *v, uint8_t *s1, uint8_t *s2, int stride, int h){
  932.     MpegEncContext *c = v;
  933.     int score1=0;
  934.     int score2=0;
  935.     int x,y;
  936.     for(y=0; y<h; y++){
  937.         for(x=0; x<16; x++){
  938.             score1+= (s1[x  ] - s2[x ])*(s1[x  ] - s2[x ]);
  939.         }
  940.         if(y+1<h){
  941.             for(x=0; x<15; x++){
  942.                 score2+= FFABS(  s1[x  ] - s1[x  +stride]
  943.                              - s1[x+1] + s1[x+1+stride])
  944.                         -FFABS(  s2[x  ] - s2[x  +stride]
  945.                              - s2[x+1] + s2[x+1+stride]);
  946.             }
  947.         }
  948.         s1+= stride;
  949.         s2+= stride;
  950.     }
  951.     if(c) return score1 + FFABS(score2)*c->avctx->nsse_weight;
  952.     else  return score1 + FFABS(score2)*8;
  953. }
  954. static int nsse8_c(void *v, uint8_t *s1, uint8_t *s2, int stride, int h){
  955.     MpegEncContext *c = v;
  956.     int score1=0;
  957.     int score2=0;
  958.     int x,y;
  959.     for(y=0; y<h; y++){
  960.         for(x=0; x<8; x++){
  961.             score1+= (s1[x  ] - s2[x ])*(s1[x  ] - s2[x ]);
  962.         }
  963.         if(y+1<h){
  964.             for(x=0; x<7; x++){
  965.                 score2+= FFABS(  s1[x  ] - s1[x  +stride]
  966.                              - s1[x+1] + s1[x+1+stride])
  967.                         -FFABS(  s2[x  ] - s2[x  +stride]
  968.                              - s2[x+1] + s2[x+1+stride]);
  969.             }
  970.         }
  971.         s1+= stride;
  972.         s2+= stride;
  973.     }
  974.     if(c) return score1 + FFABS(score2)*c->avctx->nsse_weight;
  975.     else  return score1 + FFABS(score2)*8;
  976. }
  977. static int try_8x8basis_c(int16_t rem[64], int16_t weight[64], int16_t basis[64], int scale){
  978.     int i;
  979.     unsigned int sum=0;
  980.     for(i=0; i<8*8; i++){
  981.         int b= rem[i] + ((basis[i]*scale + (1<<(BASIS_SHIFT - RECON_SHIFT-1)))>>(BASIS_SHIFT - RECON_SHIFT));
  982.         int w= weight[i];
  983.         b>>= RECON_SHIFT;
  984.         assert(-512<b && b<512);
  985.         sum += (w*b)*(w*b)>>4;
  986.     }
  987.     return sum>>2;
  988. }
  989. static void add_8x8basis_c(int16_t rem[64], int16_t basis[64], int scale){
  990.     int i;
  991.     for(i=0; i<8*8; i++){
  992.         rem[i] += (basis[i]*scale + (1<<(BASIS_SHIFT - RECON_SHIFT-1)))>>(BASIS_SHIFT - RECON_SHIFT);
  993.     }
  994. }
  995. /**
  996.  * permutes an 8x8 block.
  997.  * @param block the block which will be permuted according to the given permutation vector
  998.  * @param permutation the permutation vector
  999.  * @param last the last non zero coefficient in scantable order, used to speed the permutation up
  1000.  * @param scantable the used scantable, this is only used to speed the permutation up, the block is not
  1001.  *                  (inverse) permutated to scantable order!
  1002.  */
  1003. void ff_block_permute(DCTELEM *block, uint8_t *permutation, const uint8_t *scantable, int last)
  1004. {
  1005.     int i;
  1006.     DCTELEM temp[64];
  1007.     if(last<=0) return;
  1008.     //if(permutation[1]==1) return; //FIXME it is ok but not clean and might fail for some permutations
  1009.     for(i=0; i<=last; i++){
  1010.         const int j= scantable[i];
  1011.         temp[j]= block[j];
  1012.         block[j]=0;
  1013.     }
  1014.     for(i=0; i<=last; i++){
  1015.         const int j= scantable[i];
  1016.         const int perm_j= permutation[j];
  1017.         block[perm_j]= temp[j];
  1018.     }
  1019. }
  1020. static int zero_cmp(void *s, uint8_t *a, uint8_t *b, int stride, int h){
  1021.     return 0;
  1022. }
  1023. void ff_set_cmp(DSPContext* c, me_cmp_func *cmp, int type){
  1024.     int i;
  1025.     memset(cmp, 0, sizeof(void*)*5);
  1026.     for(i=0; i<5; i++){
  1027.         switch(type&0xFF){
  1028.         case FF_CMP_SAD:
  1029.             cmp[i]= c->sad[i];
  1030.             break;
  1031.         case FF_CMP_SATD:
  1032.             cmp[i]= c->hadamard8_diff[i];
  1033.             break;
  1034.         case FF_CMP_SSE:
  1035.             cmp[i]= c->sse[i];
  1036.             break;
  1037.         case FF_CMP_DCT:
  1038.             cmp[i]= c->dct_sad[i];
  1039.             break;
  1040.         case FF_CMP_DCT264:
  1041.             cmp[i]= c->dct264_sad[i];
  1042.             break;
  1043.         case FF_CMP_DCTMAX:
  1044.             cmp[i]= c->dct_max[i];
  1045.             break;
  1046.         case FF_CMP_PSNR:
  1047.             cmp[i]= c->quant_psnr[i];
  1048.             break;
  1049.         case FF_CMP_BIT:
  1050.             cmp[i]= c->bit[i];
  1051.             break;
  1052.         case FF_CMP_RD:
  1053.             cmp[i]= c->rd[i];
  1054.             break;
  1055.         case FF_CMP_VSAD:
  1056.             cmp[i]= c->vsad[i];
  1057.             break;
  1058.         case FF_CMP_VSSE:
  1059.             cmp[i]= c->vsse[i];
  1060.             break;
  1061.         case FF_CMP_ZERO:
  1062.             cmp[i]= zero_cmp;
  1063.             break;
  1064.         case FF_CMP_NSSE:
  1065.             cmp[i]= c->nsse[i];
  1066.             break;
  1067. #ifdef CONFIG_SNOW_ENCODER
  1068.         case FF_CMP_W53:
  1069.             cmp[i]= c->w53[i];
  1070.             break;
  1071.         case FF_CMP_W97:
  1072.             cmp[i]= c->w97[i];
  1073.             break;
  1074. #endif
  1075.         default:
  1076.             av_log(NULL, AV_LOG_ERROR,"internal error in cmp function selectionn");
  1077.         }
  1078.     }
  1079. }
  1080. static void clear_block_c(DCTELEM *block)
  1081. {
  1082.     memset(block, 0, sizeof(DCTELEM)*64);
  1083. }
  1084. /**
  1085.  * memset(blocks, 0, sizeof(DCTELEM)*6*64)
  1086.  */
  1087. static void clear_blocks_c(DCTELEM *blocks)
  1088. {
  1089.     memset(blocks, 0, sizeof(DCTELEM)*6*64);
  1090. }
  1091. static void add_bytes_c(uint8_t *dst, uint8_t *src, int w){
  1092.     long i;
  1093.     for(i=0; i<=w-(long)sizeof(long); i+=(long)sizeof(long)){
  1094.         long a = *(long*)(src+i);
  1095.         long b = *(long*)(dst+i);
  1096.         *(long*)(dst+i) = ((a&pb_7f) + (b&pb_7f)) ^ ((a^b)&pb_80);
  1097.     }
  1098.     for(; i<w; i++)
  1099.         dst[i+0] += src[i+0];
  1100. }
  1101. static void add_bytes_l2_c(uint8_t *dst, uint8_t *src1, uint8_t *src2, int w){
  1102.     long i;
  1103.     for(i=0; i<=w-(long)sizeof(long); i+=(long)sizeof(long)){
  1104.         long a = *(long*)(src1+i);
  1105.         long b = *(long*)(src2+i);
  1106.         *(long*)(dst+i) = ((a&pb_7f) + (b&pb_7f)) ^ ((a^b)&pb_80);
  1107.     }
  1108.     for(; i<w; i++)
  1109.         dst[i] = src1[i]+src2[i];
  1110. }
  1111. static void diff_bytes_c(uint8_t *dst, uint8_t *src1, uint8_t *src2, int w){
  1112.     long i;
  1113. #ifndef HAVE_FAST_UNALIGNED
  1114.     if((long)src2 & (sizeof(long)-1)){
  1115.         for(i=0; i+7<w; i+=8){
  1116.             dst[i+0] = src1[i+0]-src2[i+0];
  1117.             dst[i+1] = src1[i+1]-src2[i+1];
  1118.             dst[i+2] = src1[i+2]-src2[i+2];
  1119.             dst[i+3] = src1[i+3]-src2[i+3];
  1120.             dst[i+4] = src1[i+4]-src2[i+4];
  1121.             dst[i+5] = src1[i+5]-src2[i+5];
  1122.             dst[i+6] = src1[i+6]-src2[i+6];
  1123.             dst[i+7] = src1[i+7]-src2[i+7];
  1124.         }
  1125.     }else
  1126. #endif
  1127.     for(i=0; i<=w-(long)sizeof(long); i+=(long)sizeof(long)){
  1128.         long a = *(long*)(src1+i);
  1129.         long b = *(long*)(src2+i);
  1130.         *(long*)(dst+i) = ((a|pb_80) - (b&pb_7f)) ^ ((a^b^pb_80)&pb_80);
  1131.     }
  1132.     for(; i<w; i++)
  1133.         dst[i+0] = src1[i+0]-src2[i+0];
  1134. }
  1135. static void sub_hfyu_median_prediction_c(uint8_t *dst, uint8_t *src1, uint8_t *src2, int w, int *left, int *left_top){
  1136.     int i;
  1137.     uint8_t l, lt;
  1138.     l= *left;
  1139.     lt= *left_top;
  1140.     for(i=0; i<w; i++){
  1141.         const int pred= mid_pred(l, src1[i], (l + src1[i] - lt)&0xFF);
  1142.         lt= src1[i];
  1143.         l= src2[i];
  1144.         dst[i]= l - pred;
  1145.     }
  1146.     *left= l;
  1147.     *left_top= lt;
  1148. }
  1149. #define BUTTERFLY2(o1,o2,i1,i2) 
  1150. o1= (i1)+(i2);
  1151. o2= (i1)-(i2);
  1152. #define BUTTERFLY1(x,y) 
  1153. {
  1154.     int a,b;
  1155.     a= x;
  1156.     b= y;
  1157.     x= a+b;
  1158.     y= a-b;
  1159. }
  1160. #define BUTTERFLYA(x,y) (FFABS((x)+(y)) + FFABS((x)-(y)))
  1161. static int hadamard8_diff8x8_c(/*MpegEncContext*/ void *s, uint8_t *dst, uint8_t *src, int stride, int h){
  1162.     int i;
  1163.     int temp[64];
  1164.     int sum=0;
  1165.     assert(h==8);
  1166.     for(i=0; i<8; i++){
  1167.         //FIXME try pointer walks
  1168.         BUTTERFLY2(temp[8*i+0], temp[8*i+1], src[stride*i+0]-dst[stride*i+0],src[stride*i+1]-dst[stride*i+1]);
  1169.         BUTTERFLY2(temp[8*i+2], temp[8*i+3], src[stride*i+2]-dst[stride*i+2],src[stride*i+3]-dst[stride*i+3]);
  1170.         BUTTERFLY2(temp[8*i+4], temp[8*i+5], src[stride*i+4]-dst[stride*i+4],src[stride*i+5]-dst[stride*i+5]);
  1171.         BUTTERFLY2(temp[8*i+6], temp[8*i+7], src[stride*i+6]-dst[stride*i+6],src[stride*i+7]-dst[stride*i+7]);
  1172.         BUTTERFLY1(temp[8*i+0], temp[8*i+2]);
  1173.         BUTTERFLY1(temp[8*i+1], temp[8*i+3]);
  1174.         BUTTERFLY1(temp[8*i+4], temp[8*i+6]);
  1175.         BUTTERFLY1(temp[8*i+5], temp[8*i+7]);
  1176.         BUTTERFLY1(temp[8*i+0], temp[8*i+4]);
  1177.         BUTTERFLY1(temp[8*i+1], temp[8*i+5]);
  1178.         BUTTERFLY1(temp[8*i+2], temp[8*i+6]);
  1179.         BUTTERFLY1(temp[8*i+3], temp[8*i+7]);
  1180.     }
  1181.     for(i=0; i<8; i++){
  1182.         BUTTERFLY1(temp[8*0+i], temp[8*1+i]);
  1183.         BUTTERFLY1(temp[8*2+i], temp[8*3+i]);
  1184.         BUTTERFLY1(temp[8*4+i], temp[8*5+i]);
  1185.         BUTTERFLY1(temp[8*6+i], temp[8*7+i]);
  1186.         BUTTERFLY1(temp[8*0+i], temp[8*2+i]);
  1187.         BUTTERFLY1(temp[8*1+i], temp[8*3+i]);
  1188.         BUTTERFLY1(temp[8*4+i], temp[8*6+i]);
  1189.         BUTTERFLY1(temp[8*5+i], temp[8*7+i]);
  1190.         sum +=
  1191.              BUTTERFLYA(temp[8*0+i], temp[8*4+i])
  1192.             +BUTTERFLYA(temp[8*1+i], temp[8*5+i])
  1193.             +BUTTERFLYA(temp[8*2+i], temp[8*6+i])
  1194.             +BUTTERFLYA(temp[8*3+i], temp[8*7+i]);
  1195.     }
  1196. #if 0
  1197. static int maxi=0;
  1198. if(sum>maxi){
  1199.     maxi=sum;
  1200.     printf("MAX:%dn", maxi);
  1201. }
  1202. #endif
  1203.     return sum;
  1204. }
  1205. static int hadamard8_intra8x8_c(/*MpegEncContext*/ void *s, uint8_t *src, uint8_t *dummy, int stride, int h){
  1206.     int i;
  1207.     int temp[64];
  1208.     int sum=0;
  1209.     assert(h==8);
  1210.     for(i=0; i<8; i++){
  1211.         //FIXME try pointer walks
  1212.         BUTTERFLY2(temp[8*i+0], temp[8*i+1], src[stride*i+0],src[stride*i+1]);
  1213.         BUTTERFLY2(temp[8*i+2], temp[8*i+3], src[stride*i+2],src[stride*i+3]);
  1214.         BUTTERFLY2(temp[8*i+4], temp[8*i+5], src[stride*i+4],src[stride*i+5]);
  1215.         BUTTERFLY2(temp[8*i+6], temp[8*i+7], src[stride*i+6],src[stride*i+7]);
  1216.         BUTTERFLY1(temp[8*i+0], temp[8*i+2]);
  1217.         BUTTERFLY1(temp[8*i+1], temp[8*i+3]);
  1218.         BUTTERFLY1(temp[8*i+4], temp[8*i+6]);
  1219.         BUTTERFLY1(temp[8*i+5], temp[8*i+7]);
  1220.         BUTTERFLY1(temp[8*i+0], temp[8*i+4]);
  1221.         BUTTERFLY1(temp[8*i+1], temp[8*i+5]);
  1222.         BUTTERFLY1(temp[8*i+2], temp[8*i+6]);
  1223.         BUTTERFLY1(temp[8*i+3], temp[8*i+7]);
  1224.     }
  1225.     for(i=0; i<8; i++){
  1226.         BUTTERFLY1(temp[8*0+i], temp[8*1+i]);
  1227.         BUTTERFLY1(temp[8*2+i], temp[8*3+i]);
  1228.         BUTTERFLY1(temp[8*4+i], temp[8*5+i]);
  1229.         BUTTERFLY1(temp[8*6+i], temp[8*7+i]);
  1230.         BUTTERFLY1(temp[8*0+i], temp[8*2+i]);
  1231.         BUTTERFLY1(temp[8*1+i], temp[8*3+i]);
  1232.         BUTTERFLY1(temp[8*4+i], temp[8*6+i]);
  1233.         BUTTERFLY1(temp[8*5+i], temp[8*7+i]);
  1234.         sum +=
  1235.              BUTTERFLYA(temp[8*0+i], temp[8*4+i])
  1236.             +BUTTERFLYA(temp[8*1+i], temp[8*5+i])
  1237.             +BUTTERFLYA(temp[8*2+i], temp[8*6+i])
  1238.             +BUTTERFLYA(temp[8*3+i], temp[8*7+i]);
  1239.     }
  1240.     sum -= FFABS(temp[8*0] + temp[8*4]); // -mean
  1241.     return sum;
  1242. }
  1243. static int dct_sad8x8_c(/*MpegEncContext*/ void *c, uint8_t *src1, uint8_t *src2, int stride, int h){
  1244.     MpegEncContext * const s= (MpegEncContext *)c;
  1245.     DECLARE_ALIGNED_16(uint64_t, aligned_temp[sizeof(DCTELEM)*64/8]);
  1246.     DCTELEM * const temp= (DCTELEM*)aligned_temp;
  1247.     assert(h==8);
  1248.     s->dsp.diff_pixels(temp, src1, src2, stride);
  1249.     s->dsp.fdct(temp);
  1250.     return s->dsp.sum_abs_dctelem(temp);
  1251. }
  1252. #ifdef CONFIG_GPL
  1253. #define DCT8_1D {
  1254.     const int s07 = SRC(0) + SRC(7);
  1255.     const int s16 = SRC(1) + SRC(6);
  1256.     const int s25 = SRC(2) + SRC(5);
  1257.     const int s34 = SRC(3) + SRC(4);
  1258.     const int a0 = s07 + s34;
  1259.     const int a1 = s16 + s25;
  1260.     const int a2 = s07 - s34;
  1261.     const int a3 = s16 - s25;
  1262.     const int d07 = SRC(0) - SRC(7);
  1263.     const int d16 = SRC(1) - SRC(6);
  1264.     const int d25 = SRC(2) - SRC(5);
  1265.     const int d34 = SRC(3) - SRC(4);
  1266.     const int a4 = d16 + d25 + (d07 + (d07>>1));
  1267.     const int a5 = d07 - d34 - (d25 + (d25>>1));
  1268.     const int a6 = d07 + d34 - (d16 + (d16>>1));
  1269.     const int a7 = d16 - d25 + (d34 + (d34>>1));
  1270.     DST(0,  a0 + a1     ) ;
  1271.     DST(1,  a4 + (a7>>2)) ;
  1272.     DST(2,  a2 + (a3>>1)) ;
  1273.     DST(3,  a5 + (a6>>2)) ;
  1274.     DST(4,  a0 - a1     ) ;
  1275.     DST(5,  a6 - (a5>>2)) ;
  1276.     DST(6, (a2>>1) - a3 ) ;
  1277.     DST(7, (a4>>2) - a7 ) ;
  1278. }
  1279. static int dct264_sad8x8_c(/*MpegEncContext*/ void *c, uint8_t *src1, uint8_t *src2, int stride, int h){
  1280.     MpegEncContext * const s= (MpegEncContext *)c;
  1281.     DCTELEM dct[8][8];
  1282.     int i;
  1283.     int sum=0;
  1284.     s->dsp.diff_pixels(dct[0], src1, src2, stride);
  1285. #define SRC(x) dct[i][x]
  1286. #define DST(x,v) dct[i][x]= v
  1287.     for( i = 0; i < 8; i++ )
  1288.         DCT8_1D
  1289. #undef SRC
  1290. #undef DST
  1291. #define SRC(x) dct[x][i]
  1292. #define DST(x,v) sum += FFABS(v)
  1293.     for( i = 0; i < 8; i++ )
  1294.         DCT8_1D
  1295. #undef SRC
  1296. #undef DST
  1297.     return sum;
  1298. }
  1299. #endif
  1300. static int dct_max8x8_c(/*MpegEncContext*/ void *c, uint8_t *src1, uint8_t *src2, int stride, int h){
  1301.     MpegEncContext * const s= (MpegEncContext *)c;
  1302.     DECLARE_ALIGNED_8(uint64_t, aligned_temp[sizeof(DCTELEM)*64/8]);
  1303.     DCTELEM * const temp= (DCTELEM*)aligned_temp;
  1304.     int sum=0, i;
  1305.     assert(h==8);
  1306.     s->dsp.diff_pixels(temp, src1, src2, stride);
  1307.     s->dsp.fdct(temp);
  1308.     for(i=0; i<64; i++)
  1309.         sum= FFMAX(sum, FFABS(temp[i]));
  1310.     return sum;
  1311. }
  1312. static int quant_psnr8x8_c(/*MpegEncContext*/ void *c, uint8_t *src1, uint8_t *src2, int stride, int h){
  1313.     MpegEncContext * const s= (MpegEncContext *)c;
  1314.     DECLARE_ALIGNED_8 (uint64_t, aligned_temp[sizeof(DCTELEM)*64*2/8]);
  1315.     DCTELEM * const temp= (DCTELEM*)aligned_temp;
  1316.     DCTELEM * const bak = ((DCTELEM*)aligned_temp)+64;
  1317.     int sum=0, i;
  1318.     assert(h==8);
  1319.     s->mb_intra=0;
  1320.     s->dsp.diff_pixels(temp, src1, src2, stride);
  1321.     memcpy(bak, temp, 64*sizeof(DCTELEM));
  1322.     s->block_last_index[0/*FIXME*/]= s->fast_dct_quantize(s, temp, 0/*FIXME*/, s->qscale, &i);
  1323.     s->dct_unquantize_inter(s, temp, 0, s->qscale);
  1324.     ff_simple_idct(temp); //FIXME
  1325.     for(i=0; i<64; i++)
  1326.         sum+= (temp[i]-bak[i])*(temp[i]-bak[i]);
  1327.     return sum;
  1328. }
  1329. static int rd8x8_c(/*MpegEncContext*/ void *c, uint8_t *src1, uint8_t *src2, int stride, int h){
  1330.     MpegEncContext * const s= (MpegEncContext *)c;
  1331.     const uint8_t *scantable= s->intra_scantable.permutated;
  1332.     DECLARE_ALIGNED_8 (uint64_t, aligned_temp[sizeof(DCTELEM)*64/8]);
  1333.     //DECLARE_ALIGNED_8 (uint64_t, aligned_bak[stride]);
  1334. uint64_t* aligned_bak = (uint64_t*)av_malloc(stride*sizeof(uint64_t));
  1335.     DCTELEM * const temp= (DCTELEM*)aligned_temp;
  1336.     uint8_t * const bak= (uint8_t*)aligned_bak;
  1337.     int i, last, run, bits, level, distortion, start_i;
  1338.     const int esc_length= s->ac_esc_length;
  1339.     uint8_t * length;
  1340.     uint8_t * last_length;
  1341.     assert(h==8);
  1342.     for(i=0; i<8; i++){
  1343.         ((uint32_t*)(bak + i*stride))[0]= ((uint32_t*)(src2 + i*stride))[0];
  1344.         ((uint32_t*)(bak + i*stride))[1]= ((uint32_t*)(src2 + i*stride))[1];
  1345.     }
  1346.     s->dsp.diff_pixels(temp, src1, src2, stride);
  1347.     s->block_last_index[0/*FIXME*/]= last= s->fast_dct_quantize(s, temp, 0/*FIXME*/, s->qscale, &i);
  1348.     bits=0;
  1349.     if (s->mb_intra) {
  1350.         start_i = 1;
  1351.         length     = s->intra_ac_vlc_length;
  1352.         last_length= s->intra_ac_vlc_last_length;
  1353.         bits+= s->luma_dc_vlc_length[temp[0] + 256]; //FIXME chroma
  1354.     } else {
  1355.         start_i = 0;
  1356.         length     = s->inter_ac_vlc_length;
  1357.         last_length= s->inter_ac_vlc_last_length;
  1358.     }
  1359.     if(last>=start_i){
  1360.         run=0;
  1361.         for(i=start_i; i<last; i++){
  1362.             int j= scantable[i];
  1363.             level= temp[j];
  1364.             if(level){
  1365.                 level+=64;
  1366.                 if((level&(~127)) == 0){
  1367.                     bits+= length[UNI_AC_ENC_INDEX(run, level)];
  1368.                 }else
  1369.                     bits+= esc_length;
  1370.                 run=0;
  1371.             }else
  1372.                 run++;
  1373.         }
  1374.         i= scantable[last];
  1375.         level= temp[i] + 64;
  1376.         assert(level - 64);
  1377.         if((level&(~127)) == 0){
  1378.             bits+= last_length[UNI_AC_ENC_INDEX(run, level)];
  1379.         }else
  1380.             bits+= esc_length;
  1381.     }
  1382.     if(last>=0){
  1383.         if(s->mb_intra)
  1384.             s->dct_unquantize_intra(s, temp, 0, s->qscale);
  1385.         else
  1386.             s->dct_unquantize_inter(s, temp, 0, s->qscale);
  1387.     }
  1388.     s->dsp.idct_add(bak, stride, temp);
  1389.     distortion= s->dsp.sse[1](NULL, bak, src1, stride, 8);
  1390. av_free(aligned_bak);
  1391.     return distortion + ((bits*s->qscale*s->qscale*109 + 64)>>7);
  1392. }
  1393. static int bit8x8_c(/*MpegEncContext*/ void *c, uint8_t *src1, uint8_t *src2, int stride, int h){
  1394.     MpegEncContext * const s= (MpegEncContext *)c;
  1395.     const uint8_t *scantable= s->intra_scantable.permutated;
  1396.     DECLARE_ALIGNED_8 (uint64_t, aligned_temp[sizeof(DCTELEM)*64/8]);
  1397.     DCTELEM * const temp= (DCTELEM*)aligned_temp;
  1398.     int i, last, run, bits, level, start_i;
  1399.     const int esc_length= s->ac_esc_length;
  1400.     uint8_t * length;
  1401.     uint8_t * last_length;
  1402.     assert(h==8);
  1403.     s->dsp.diff_pixels(temp, src1, src2, stride);
  1404.     s->block_last_index[0/*FIXME*/]= last= s->fast_dct_quantize(s, temp, 0/*FIXME*/, s->qscale, &i);
  1405.     bits=0;
  1406.     if (s->mb_intra) {
  1407.         start_i = 1;
  1408.         length     = s->intra_ac_vlc_length;
  1409.         last_length= s->intra_ac_vlc_last_length;
  1410.         bits+= s->luma_dc_vlc_length[temp[0] + 256]; //FIXME chroma
  1411.     } else {
  1412.         start_i = 0;
  1413.         length     = s->inter_ac_vlc_length;
  1414.         last_length= s->inter_ac_vlc_last_length;
  1415.     }
  1416.     if(last>=start_i){
  1417.         run=0;
  1418.         for(i=start_i; i<last; i++){
  1419.             int j= scantable[i];
  1420.             level= temp[j];
  1421.             if(level){
  1422.                 level+=64;
  1423.                 if((level&(~127)) == 0){
  1424.                     bits+= length[UNI_AC_ENC_INDEX(run, level)];
  1425.                 }else
  1426.                     bits+= esc_length;
  1427.                 run=0;
  1428.             }else
  1429.                 run++;
  1430.         }
  1431.         i= scantable[last];
  1432.         level= temp[i] + 64;
  1433.         assert(level - 64);
  1434.         if((level&(~127)) == 0){
  1435.             bits+= last_length[UNI_AC_ENC_INDEX(run, level)];
  1436.         }else
  1437.             bits+= esc_length;
  1438.     }
  1439.     return bits;
  1440. }
  1441. static int vsad_intra16_c(/*MpegEncContext*/ void *c, uint8_t *s, uint8_t *dummy, int stride, int h){
  1442.     int score=0;
  1443.     int x,y;
  1444.     for(y=1; y<h; y++){
  1445.         for(x=0; x<16; x+=4){
  1446.             score+= FFABS(s[x  ] - s[x  +stride]) + FFABS(s[x+1] - s[x+1+stride])
  1447.                    +FFABS(s[x+2] - s[x+2+stride]) + FFABS(s[x+3] - s[x+3+stride]);
  1448.         }
  1449.         s+= stride;
  1450.     }
  1451.     return score;
  1452. }
  1453. static int vsad16_c(/*MpegEncContext*/ void *c, uint8_t *s1, uint8_t *s2, int stride, int h){
  1454.     int score=0;
  1455.     int x,y;
  1456.     for(y=1; y<h; y++){
  1457.         for(x=0; x<16; x++){
  1458.             score+= FFABS(s1[x  ] - s2[x ] - s1[x  +stride] + s2[x +stride]);
  1459.         }
  1460.         s1+= stride;
  1461.         s2+= stride;
  1462.     }
  1463.     return score;
  1464. }
  1465. #define SQ(a) ((a)*(a))
  1466. static int vsse_intra16_c(/*MpegEncContext*/ void *c, uint8_t *s, uint8_t *dummy, int stride, int h){
  1467.     int score=0;
  1468.     int x,y;
  1469.     for(y=1; y<h; y++){
  1470.         for(x=0; x<16; x+=4){
  1471.             score+= SQ(s[x  ] - s[x  +stride]) + SQ(s[x+1] - s[x+1+stride])
  1472.                    +SQ(s[x+2] - s[x+2+stride]) + SQ(s[x+3] - s[x+3+stride]);
  1473.         }
  1474.         s+= stride;
  1475.     }
  1476.     return score;
  1477. }
  1478. static int vsse16_c(/*MpegEncContext*/ void *c, uint8_t *s1, uint8_t *s2, int stride, int h){
  1479.     int score=0;
  1480.     int x,y;
  1481.     for(y=1; y<h; y++){
  1482.         for(x=0; x<16; x++){
  1483.             score+= SQ(s1[x  ] - s2[x ] - s1[x  +stride] + s2[x +stride]);
  1484.         }
  1485.         s1+= stride;
  1486.         s2+= stride;
  1487.     }
  1488.     return score;
  1489. }
  1490. static int ssd_int8_vs_int16_c(const int8_t *pix1, const int16_t *pix2,
  1491.                                int size){
  1492.     int score=0;
  1493.     int i;
  1494.     for(i=0; i<size; i++)
  1495.         score += (pix1[i]-pix2[i])*(pix1[i]-pix2[i]);
  1496.     return score;
  1497. }
  1498. WRAPPER8_16_SQ(hadamard8_diff8x8_c, hadamard8_diff16_c)
  1499. WRAPPER8_16_SQ(hadamard8_intra8x8_c, hadamard8_intra16_c)
  1500. WRAPPER8_16_SQ(dct_sad8x8_c, dct_sad16_c)
  1501. #ifdef CONFIG_GPL
  1502. WRAPPER8_16_SQ(dct264_sad8x8_c, dct264_sad16_c)
  1503. #endif
  1504. WRAPPER8_16_SQ(dct_max8x8_c, dct_max16_c)
  1505. WRAPPER8_16_SQ(quant_psnr8x8_c, quant_psnr16_c)
  1506. WRAPPER8_16_SQ(rd8x8_c, rd16_c)
  1507. WRAPPER8_16_SQ(bit8x8_c, bit16_c)
  1508. static void vector_fmul_c(float *dst, const float *src, int len){
  1509.     int i;
  1510.     for(i=0; i<len; i++)
  1511.         dst[i] *= src[i];
  1512. }
  1513. static void vector_fmul_reverse_c(float *dst, const float *src0, const float *src1, int len){
  1514.     int i;
  1515.     src1 += len-1;
  1516.     for(i=0; i<len; i++)
  1517.         dst[i] = src0[i] * src1[-i];
  1518. }
  1519. void ff_vector_fmul_add_add_c(float *dst, const float *src0, const float *src1, const float *src2, int src3, int len, int step){
  1520.     int i;
  1521.     for(i=0; i<len; i++)
  1522.         dst[i*step] = src0[i] * src1[i] + src2[i] + src3;
  1523. }
  1524. void ff_vector_fmul_window_c(float *dst, const float *src0, const float *src1, const float *win, float add_bias, int len){
  1525.     int i,j;
  1526.     dst += len;
  1527.     win += len;
  1528.     src0+= len;
  1529.     for(i=-len, j=len-1; i<0; i++, j--) {
  1530.         float s0 = src0[i];
  1531.         float s1 = src1[j];
  1532.         float wi = win[i];
  1533.         float wj = win[j];
  1534.         dst[i] = s0*wj - s1*wi + add_bias;
  1535.         dst[j] = s0*wi + s1*wj + add_bias;
  1536.     }
  1537. }
  1538. static void int32_to_float_fmul_scalar_c(float *dst, const int *src, float mul, int len){
  1539.     int i;
  1540.     for(i=0; i<len; i++)
  1541.         dst[i] = src[i] * mul;
  1542. }
  1543. static av_always_inline int float_to_int16_one(const float *src){
  1544.     int_fast32_t tmp = *(const int32_t*)src;
  1545.     if(tmp & 0xf0000){
  1546.         tmp = (0x43c0ffff - tmp)>>31;
  1547.         // is this faster on some gcc/cpu combinations?
  1548. //      if(tmp > 0x43c0ffff) tmp = 0xFFFF;
  1549. //      else                 tmp = 0;
  1550.     }
  1551.     return tmp - 0x8000;
  1552. }
  1553. void ff_float_to_int16_c(int16_t *dst, const float *src, long len){
  1554.     int i;
  1555.     for(i=0; i<len; i++)
  1556.         dst[i] = float_to_int16_one(src+i);
  1557. }
  1558. void ff_float_to_int16_interleave_c(int16_t *dst, const float **src, long len, int channels){
  1559.     int i,j,c;
  1560.     if(channels==2){
  1561.         for(i=0; i<len; i++){
  1562.             dst[2*i]   = float_to_int16_one(src[0]+i);
  1563.             dst[2*i+1] = float_to_int16_one(src[1]+i);
  1564.         }
  1565.     }else{
  1566.         for(c=0; c<channels; c++)
  1567.             for(i=0, j=c; i<len; i++, j+=channels)
  1568.                 dst[j] = float_to_int16_one(src[c]+i);
  1569.     }
  1570. }
  1571. static void add_int16_c(int16_t * v1, int16_t * v2, int order)
  1572. {
  1573.     while (order--)
  1574.        *v1++ += *v2++;
  1575. }
  1576. static void sub_int16_c(int16_t * v1, int16_t * v2, int order)
  1577. {
  1578.     while (order--)
  1579.         *v1++ -= *v2++;
  1580. }
  1581. static int32_t scalarproduct_int16_c(int16_t * v1, int16_t * v2, int order, int shift)
  1582. {
  1583.     int res = 0;
  1584.     while (order--)
  1585.         res += (*v1++ * *v2++) >> shift;
  1586.     return res;
  1587. }
  1588. #define W0 2048
  1589. #define W1 2841 /* 2048*sqrt (2)*cos (1*pi/16) */
  1590. #define W2 2676 /* 2048*sqrt (2)*cos (2*pi/16) */
  1591. #define W3 2408 /* 2048*sqrt (2)*cos (3*pi/16) */
  1592. #define W4 2048 /* 2048*sqrt (2)*cos (4*pi/16) */
  1593. #define W5 1609 /* 2048*sqrt (2)*cos (5*pi/16) */
  1594. #define W6 1108 /* 2048*sqrt (2)*cos (6*pi/16) */
  1595. #define W7 565  /* 2048*sqrt (2)*cos (7*pi/16) */
  1596. static void wmv2_idct_row(short * b)
  1597. {
  1598.     int s1,s2;
  1599.     int a0,a1,a2,a3,a4,a5,a6,a7;
  1600.     /*step 1*/
  1601.     a1 = W1*b[1]+W7*b[7];
  1602.     a7 = W7*b[1]-W1*b[7];
  1603.     a5 = W5*b[5]+W3*b[3];
  1604.     a3 = W3*b[5]-W5*b[3];
  1605.     a2 = W2*b[2]+W6*b[6];
  1606.     a6 = W6*b[2]-W2*b[6];
  1607.     a0 = W0*b[0]+W0*b[4];
  1608.     a4 = W0*b[0]-W0*b[4];
  1609.     /*step 2*/
  1610.     s1 = (181*(a1-a5+a7-a3)+128)>>8;//1,3,5,7,
  1611.     s2 = (181*(a1-a5-a7+a3)+128)>>8;
  1612.     /*step 3*/
  1613.     b[0] = (a0+a2+a1+a5 + (1<<7))>>8;
  1614.     b[1] = (a4+a6 +s1   + (1<<7))>>8;
  1615.     b[2] = (a4-a6 +s2   + (1<<7))>>8;
  1616.     b[3] = (a0-a2+a7+a3 + (1<<7))>>8;
  1617.     b[4] = (a0-a2-a7-a3 + (1<<7))>>8;
  1618.     b[5] = (a4-a6 -s2   + (1<<7))>>8;
  1619.     b[6] = (a4+a6 -s1   + (1<<7))>>8;
  1620.     b[7] = (a0+a2-a1-a5 + (1<<7))>>8;
  1621. }
  1622. static void wmv2_idct_col(short * b)
  1623. {
  1624.     int s1,s2;
  1625.     int a0,a1,a2,a3,a4,a5,a6,a7;
  1626.     /*step 1, with extended precision*/
  1627.     a1 = (W1*b[8*1]+W7*b[8*7] + 4)>>3;
  1628.     a7 = (W7*b[8*1]-W1*b[8*7] + 4)>>3;
  1629.     a5 = (W5*b[8*5]+W3*b[8*3] + 4)>>3;
  1630.     a3 = (W3*b[8*5]-W5*b[8*3] + 4)>>3;
  1631.     a2 = (W2*b[8*2]+W6*b[8*6] + 4)>>3;
  1632.     a6 = (W6*b[8*2]-W2*b[8*6] + 4)>>3;
  1633.     a0 = (W0*b[8*0]+W0*b[8*4]    )>>3;
  1634.     a4 = (W0*b[8*0]-W0*b[8*4]    )>>3;
  1635.     /*step 2*/
  1636.     s1 = (181*(a1-a5+a7-a3)+128)>>8;
  1637.     s2 = (181*(a1-a5-a7+a3)+128)>>8;
  1638.     /*step 3*/
  1639.     b[8*0] = (a0+a2+a1+a5 + (1<<13))>>14;
  1640.     b[8*1] = (a4+a6 +s1   + (1<<13))>>14;
  1641.     b[8*2] = (a4-a6 +s2   + (1<<13))>>14;
  1642.     b[8*3] = (a0-a2+a7+a3 + (1<<13))>>14;
  1643.     b[8*4] = (a0-a2-a7-a3 + (1<<13))>>14;
  1644.     b[8*5] = (a4-a6 -s2   + (1<<13))>>14;
  1645.     b[8*6] = (a4+a6 -s1   + (1<<13))>>14;
  1646.     b[8*7] = (a0+a2-a1-a5 + (1<<13))>>14;
  1647. }
  1648. void ff_wmv2_idct_c(short * block){
  1649.     int i;
  1650.     for(i=0;i<64;i+=8){
  1651.         wmv2_idct_row(block+i);
  1652.     }
  1653.     for(i=0;i<8;i++){
  1654.         wmv2_idct_col(block+i);
  1655.     }
  1656. }
  1657. /* XXX: those functions should be suppressed ASAP when all IDCTs are
  1658.  converted */
  1659. /*
  1660. static void ff_wmv2_idct_put_c(uint8_t *dest, int line_size, DCTELEM *block)
  1661. {
  1662.     ff_wmv2_idct_c(block);
  1663.     put_pixels_clamped_c(block, dest, line_size);
  1664. }
  1665. static void ff_wmv2_idct_add_c(uint8_t *dest, int line_size, DCTELEM *block)
  1666. {
  1667.     ff_wmv2_idct_c(block);
  1668.     add_pixels_clamped_c(block, dest, line_size);
  1669. }
  1670. */
  1671. static void ff_jref_idct_put(uint8_t *dest, int line_size, DCTELEM *block)
  1672. {
  1673.     j_rev_dct (block);
  1674.     put_pixels_clamped_c(block, dest, line_size);
  1675. }
  1676. static void ff_jref_idct_add(uint8_t *dest, int line_size, DCTELEM *block)
  1677. {
  1678.     j_rev_dct (block);
  1679.     add_pixels_clamped_c(block, dest, line_size);
  1680. }
  1681. static void ff_jref_idct4_put(uint8_t *dest, int line_size, DCTELEM *block)
  1682. {
  1683.     j_rev_dct4 (block);
  1684.     put_pixels_clamped4_c(block, dest, line_size);
  1685. }
  1686. static void ff_jref_idct4_add(uint8_t *dest, int line_size, DCTELEM *block)
  1687. {
  1688.     j_rev_dct4 (block);
  1689.     add_pixels_clamped4_c(block, dest, line_size);
  1690. }
  1691. static void ff_jref_idct2_put(uint8_t *dest, int line_size, DCTELEM *block)
  1692. {
  1693.     j_rev_dct2 (block);
  1694.     put_pixels_clamped2_c(block, dest, line_size);
  1695. }
  1696. static void ff_jref_idct2_add(uint8_t *dest, int line_size, DCTELEM *block)
  1697. {
  1698.     j_rev_dct2 (block);
  1699.     add_pixels_clamped2_c(block, dest, line_size);
  1700. }
  1701. static void ff_jref_idct1_put(uint8_t *dest, int line_size, DCTELEM *block)
  1702. {
  1703.     uint8_t *cm = ff_cropTbl + MAX_NEG_CROP;
  1704.     dest[0] = cm[(block[0] + 4)>>3];
  1705. }
  1706. static void ff_jref_idct1_add(uint8_t *dest, int line_size, DCTELEM *block)
  1707. {
  1708.     uint8_t *cm = ff_cropTbl + MAX_NEG_CROP;
  1709.     dest[0] = cm[dest[0] + ((block[0] + 4)>>3)];
  1710. }
  1711. static void just_return(void *mem av_unused, int stride av_unused, int h av_unused) { return; }
  1712. /* init static data */
  1713. void dsputil_static_init(void)
  1714. {
  1715.     int i;
  1716.     for(i=0;i<256;i++) ff_cropTbl[i + MAX_NEG_CROP] = i;
  1717.     for(i=0;i<MAX_NEG_CROP;i++) {
  1718.         ff_cropTbl[i] = 0;
  1719.         ff_cropTbl[i + MAX_NEG_CROP + 256] = 255;
  1720.     }
  1721.     for(i=0;i<512;i++) {
  1722.         ff_squareTbl[i] = (i - 256) * (i - 256);
  1723.     }
  1724.     for(i=0; i<64; i++) inv_zigzag_direct16[ff_zigzag_direct[i]]= i+1;
  1725. }
  1726. int ff_check_alignment(void){
  1727.     static int did_fail=0;
  1728.     DECLARE_ALIGNED_16(int, aligned);
  1729.     if((long)&aligned & 15){
  1730.         if(!did_fail){
  1731. #if defined(HAVE_MMX) || defined(HAVE_ALTIVEC)
  1732.             av_log(NULL, AV_LOG_ERROR,
  1733.                 "Compiler did not align stack variables. Libavcodec has been miscompiledn"
  1734.                 "and may be very slow or crash. This is not a bug in libavcodec,n"
  1735.                 "but in the compiler. You may try recompiling using gcc >= 4.2.n"
  1736.                 "Do not report crashes to FFmpeg developers.n");
  1737. #endif
  1738.             did_fail=1;
  1739.         }
  1740.         return -1;
  1741.     }
  1742.     return 0;
  1743. }
  1744. void dsputil_init(DSPContext* c, AVCodecContext *avctx)
  1745. {
  1746.     int i;
  1747.     ff_check_alignment();
  1748. #ifdef CONFIG_ENCODERS
  1749.     if(avctx->dct_algo==FF_DCT_FASTINT) {
  1750.         c->fdct = fdct_ifast;
  1751.         c->fdct248 = fdct_ifast248;
  1752.     }
  1753.     else if(avctx->dct_algo==FF_DCT_FAAN) {
  1754.         c->fdct = ff_faandct;
  1755.         c->fdct248 = ff_faandct248;
  1756.     }
  1757.     else {
  1758.         c->fdct = ff_jpeg_fdct_islow; //slow/accurate/default
  1759.         c->fdct248 = ff_fdct248_islow;
  1760.     }
  1761. #endif //CONFIG_ENCODERS
  1762.     if(avctx->lowres==1){
  1763.         if(avctx->idct_algo==FF_IDCT_INT || avctx->idct_algo==FF_IDCT_AUTO || !ENABLE_H264_DECODER){
  1764.             c->idct_put= ff_jref_idct4_put;
  1765.             c->idct_add= ff_jref_idct4_add;
  1766.         }else{
  1767.             c->idct_put= ff_h264_lowres_idct_put_c;
  1768.             c->idct_add= ff_h264_lowres_idct_add_c;
  1769.         }
  1770.         c->idct    = j_rev_dct4;
  1771.         c->idct_permutation_type= FF_NO_IDCT_PERM;
  1772.     }else if(avctx->lowres==2){
  1773.         c->idct_put= ff_jref_idct2_put;
  1774.         c->idct_add= ff_jref_idct2_add;
  1775.         c->idct    = j_rev_dct2;
  1776.         c->idct_permutation_type= FF_NO_IDCT_PERM;
  1777.     }else if(avctx->lowres==3){
  1778.         c->idct_put= ff_jref_idct1_put;
  1779.         c->idct_add= ff_jref_idct1_add;
  1780.         c->idct    = j_rev_dct1;
  1781.         c->idct_permutation_type= FF_NO_IDCT_PERM;
  1782.     }
  1783. else
  1784. {
  1785.         if(avctx->idct_algo==FF_IDCT_INT)
  1786. {
  1787.             c->idct_put= ff_jref_idct_put;
  1788.             c->idct_add= ff_jref_idct_add;
  1789.             c->idct    = j_rev_dct;
  1790.             c->idct_permutation_type= FF_LIBMPEG2_IDCT_PERM;
  1791.         }
  1792. /*else if((ENABLE_VP3_DECODER || ENABLE_VP5_DECODER || ENABLE_VP6_DECODER || ENABLE_THEORA_DECODER ) && avctx->idct_algo==FF_IDCT_VP3)
  1793. {
  1794.             c->idct_put= ff_vp3_idct_put_c;
  1795.             c->idct_add= ff_vp3_idct_add_c;
  1796.             c->idct    = ff_vp3_idct_c;
  1797.             c->idct_permutation_type= FF_NO_IDCT_PERM;
  1798.         }
  1799. else if(avctx->idct_algo==FF_IDCT_WMV2)
  1800. {
  1801.             c->idct_put= ff_wmv2_idct_put_c;
  1802.             c->idct_add= ff_wmv2_idct_add_c;
  1803.             c->idct    = ff_wmv2_idct_c;
  1804.             c->idct_permutation_type= FF_NO_IDCT_PERM;
  1805.         }
  1806. else if(avctx->idct_algo==FF_IDCT_FAAN)
  1807. {
  1808.             c->idct_put= ff_faanidct_put;
  1809.             c->idct_add= ff_faanidct_add;
  1810.             c->idct    = ff_faanidct;
  1811.             c->idct_permutation_type= FF_NO_IDCT_PERM;
  1812.         }
  1813. else if(ENABLE_EATGQ_DECODER && avctx->idct_algo==FF_IDCT_EA) 
  1814. {
  1815.             c->idct_put= ff_ea_idct_put_c;
  1816.             c->idct_permutation_type= FF_NO_IDCT_PERM;
  1817.         }*/
  1818. else
  1819. { //accurate/default
  1820.             c->idct_put= ff_simple_idct_put;
  1821.             c->idct_add= ff_simple_idct_add;
  1822.             c->idct    = ff_simple_idct;
  1823.             c->idct_permutation_type= FF_NO_IDCT_PERM;
  1824.         }
  1825.     }
  1826.     if (ENABLE_H264_DECODER) {
  1827.         c->h264_idct_add= ff_h264_idct_add_c;
  1828.         c->h264_idct8_add= ff_h264_idct8_add_c;
  1829.         c->h264_idct_dc_add= ff_h264_idct_dc_add_c;
  1830.         c->h264_idct8_dc_add= ff_h264_idct8_dc_add_c;
  1831.         c->h264_idct_add16     = ff_h264_idct_add16_c;
  1832.         c->h264_idct8_add4     = ff_h264_idct8_add4_c;
  1833.         c->h264_idct_add8      = ff_h264_idct_add8_c;
  1834.         c->h264_idct_add16intra= ff_h264_idct_add16intra_c;
  1835.     }
  1836.     c->get_pixels = get_pixels_c;
  1837.     c->diff_pixels = diff_pixels_c;
  1838.     c->put_pixels_clamped = put_pixels_clamped_c;
  1839.     c->put_signed_pixels_clamped = put_signed_pixels_clamped_c;
  1840.     c->add_pixels_clamped = add_pixels_clamped_c;
  1841.     c->add_pixels8 = add_pixels8_c;
  1842.     c->add_pixels4 = add_pixels4_c;
  1843.     c->sum_abs_dctelem = sum_abs_dctelem_c;
  1844.     c->gmc1 = gmc1_c;
  1845.     c->gmc = ff_gmc_c;
  1846.     c->clear_block = clear_block_c;
  1847.     c->clear_blocks = clear_blocks_c;
  1848.     c->pix_sum = pix_sum_c;
  1849.     c->pix_norm1 = pix_norm1_c;
  1850.     /* TODO [0] 16  [1] 8 */
  1851.     c->pix_abs[0][0] = pix_abs16_c;
  1852.     c->pix_abs[0][1] = pix_abs16_x2_c;
  1853.     c->pix_abs[0][2] = pix_abs16_y2_c;
  1854.     c->pix_abs[0][3] = pix_abs16_xy2_c;
  1855.     c->pix_abs[1][0] = pix_abs8_c;
  1856.     c->pix_abs[1][1] = pix_abs8_x2_c;
  1857.     c->pix_abs[1][2] = pix_abs8_y2_c;
  1858.     c->pix_abs[1][3] = pix_abs8_xy2_c;
  1859. #define dspfunc(PFX, IDX, NUM) 
  1860.     c->PFX ## _pixels_tab[IDX][0] = PFX ## _pixels ## NUM ## _c;     
  1861.     c->PFX ## _pixels_tab[IDX][1] = PFX ## _pixels ## NUM ## _x2_c;  
  1862.     c->PFX ## _pixels_tab[IDX][2] = PFX ## _pixels ## NUM ## _y2_c;  
  1863.     c->PFX ## _pixels_tab[IDX][3] = PFX ## _pixels ## NUM ## _xy2_c
  1864.     dspfunc(put, 0, 16);
  1865.     dspfunc(put_no_rnd, 0, 16);
  1866.     dspfunc(put, 1, 8);
  1867.     dspfunc(put_no_rnd, 1, 8);
  1868.     dspfunc(put, 2, 4);
  1869.     dspfunc(put, 3, 2);
  1870.     dspfunc(avg, 0, 16);
  1871.     dspfunc(avg_no_rnd, 0, 16);
  1872.     dspfunc(avg, 1, 8);
  1873.     dspfunc(avg_no_rnd, 1, 8);
  1874.     dspfunc(avg, 2, 4);
  1875.     dspfunc(avg, 3, 2);
  1876. #undef dspfunc
  1877.     c->put_no_rnd_pixels_l2[0]= put_no_rnd_pixels16_l2_c;
  1878.     c->put_no_rnd_pixels_l2[1]= put_no_rnd_pixels8_l2_c;
  1879.     c->put_tpel_pixels_tab[ 0] = put_tpel_pixels_mc00_c;
  1880.     c->put_tpel_pixels_tab[ 1] = put_tpel_pixels_mc10_c;
  1881.     c->put_tpel_pixels_tab[ 2] = put_tpel_pixels_mc20_c;
  1882.     c->put_tpel_pixels_tab[ 4] = put_tpel_pixels_mc01_c;
  1883.     c->put_tpel_pixels_tab[ 5] = put_tpel_pixels_mc11_c;
  1884.     c->put_tpel_pixels_tab[ 6] = put_tpel_pixels_mc21_c;
  1885.     c->put_tpel_pixels_tab[ 8] = put_tpel_pixels_mc02_c;
  1886.     c->put_tpel_pixels_tab[ 9] = put_tpel_pixels_mc12_c;
  1887.     c->put_tpel_pixels_tab[10] = put_tpel_pixels_mc22_c;
  1888.     c->avg_tpel_pixels_tab[ 0] = avg_tpel_pixels_mc00_c;
  1889.     c->avg_tpel_pixels_tab[ 1] = avg_tpel_pixels_mc10_c;
  1890.     c->avg_tpel_pixels_tab[ 2] = avg_tpel_pixels_mc20_c;
  1891.     c->avg_tpel_pixels_tab[ 4] = avg_tpel_pixels_mc01_c;
  1892.     c->avg_tpel_pixels_tab[ 5] = avg_tpel_pixels_mc11_c;
  1893.     c->avg_tpel_pixels_tab[ 6] = avg_tpel_pixels_mc21_c;
  1894.     c->avg_tpel_pixels_tab[ 8] = avg_tpel_pixels_mc02_c;
  1895.     c->avg_tpel_pixels_tab[ 9] = avg_tpel_pixels_mc12_c;
  1896.     c->avg_tpel_pixels_tab[10] = avg_tpel_pixels_mc22_c;
  1897. #define dspfunc(PFX, IDX, NUM) 
  1898.     c->PFX ## _pixels_tab[IDX][ 0] = PFX ## NUM ## _mc00_c; 
  1899.     c->PFX ## _pixels_tab[IDX][ 1] = PFX ## NUM ## _mc10_c; 
  1900.     c->PFX ## _pixels_tab[IDX][ 2] = PFX ## NUM ## _mc20_c; 
  1901.     c->PFX ## _pixels_tab[IDX][ 3] = PFX ## NUM ## _mc30_c; 
  1902.     c->PFX ## _pixels_tab[IDX][ 4] = PFX ## NUM ## _mc01_c; 
  1903.     c->PFX ## _pixels_tab[IDX][ 5] = PFX ## NUM ## _mc11_c; 
  1904.     c->PFX ## _pixels_tab[IDX][ 6] = PFX ## NUM ## _mc21_c; 
  1905.     c->PFX ## _pixels_tab[IDX][ 7] = PFX ## NUM ## _mc31_c; 
  1906.     c->PFX ## _pixels_tab[IDX][ 8] = PFX ## NUM ## _mc02_c; 
  1907.     c->PFX ## _pixels_tab[IDX][ 9] = PFX ## NUM ## _mc12_c; 
  1908.     c->PFX ## _pixels_tab[IDX][10] = PFX ## NUM ## _mc22_c; 
  1909.     c->PFX ## _pixels_tab[IDX][11] = PFX ## NUM ## _mc32_c; 
  1910.     c->PFX ## _pixels_tab[IDX][12] = PFX ## NUM ## _mc03_c; 
  1911.     c->PFX ## _pixels_tab[IDX][13] = PFX ## NUM ## _mc13_c; 
  1912.     c->PFX ## _pixels_tab[IDX][14] = PFX ## NUM ## _mc23_c; 
  1913.     c->PFX ## _pixels_tab[IDX][15] = PFX ## NUM ## _mc33_c
  1914.     dspfunc(put_qpel, 0, 16);
  1915.     dspfunc(put_no_rnd_qpel, 0, 16);
  1916.     dspfunc(avg_qpel, 0, 16);
  1917.     /* dspfunc(avg_no_rnd_qpel, 0, 16); */
  1918.     dspfunc(put_qpel, 1, 8);
  1919.     dspfunc(put_no_rnd_qpel, 1, 8);
  1920.     dspfunc(avg_qpel, 1, 8);
  1921.     /* dspfunc(avg_no_rnd_qpel, 1, 8); */
  1922.     dspfunc(put_h264_qpel, 0, 16);
  1923.     dspfunc(put_h264_qpel, 1, 8);
  1924.     dspfunc(put_h264_qpel, 2, 4);
  1925.     dspfunc(put_h264_qpel, 3, 2);
  1926.     dspfunc(avg_h264_qpel, 0, 16);
  1927.     dspfunc(avg_h264_qpel, 1, 8);
  1928.     dspfunc(avg_h264_qpel, 2, 4);
  1929. #undef dspfunc
  1930.     c->put_h264_chroma_pixels_tab[0]= put_h264_chroma_mc8_c;
  1931.     c->put_h264_chroma_pixels_tab[1]= put_h264_chroma_mc4_c;
  1932.     c->put_h264_chroma_pixels_tab[2]= put_h264_chroma_mc2_c;
  1933.     c->avg_h264_chroma_pixels_tab[0]= avg_h264_chroma_mc8_c;
  1934.     c->avg_h264_chroma_pixels_tab[1]= avg_h264_chroma_mc4_c;
  1935.     c->avg_h264_chroma_pixels_tab[2]= avg_h264_chroma_mc2_c;
  1936.     c->put_no_rnd_h264_chroma_pixels_tab[0]= put_no_rnd_h264_chroma_mc8_c;
  1937.     c->weight_h264_pixels_tab[0]= weight_h264_pixels16x16_c;
  1938.     c->weight_h264_pixels_tab[1]= weight_h264_pixels16x8_c;
  1939.     c->weight_h264_pixels_tab[2]= weight_h264_pixels8x16_c;
  1940.     c->weight_h264_pixels_tab[3]= weight_h264_pixels8x8_c;
  1941.     c->weight_h264_pixels_tab[4]= weight_h264_pixels8x4_c;
  1942.     c->weight_h264_pixels_tab[5]= weight_h264_pixels4x8_c;
  1943.     c->weight_h264_pixels_tab[6]= weight_h264_pixels4x4_c;
  1944.     c->weight_h264_pixels_tab[7]= weight_h264_pixels4x2_c;
  1945.     c->weight_h264_pixels_tab[8]= weight_h264_pixels2x4_c;
  1946.     c->weight_h264_pixels_tab[9]= weight_h264_pixels2x2_c;
  1947.     c->biweight_h264_pixels_tab[0]= biweight_h264_pixels16x16_c;
  1948.     c->biweight_h264_pixels_tab[1]= biweight_h264_pixels16x8_c;
  1949.     c->biweight_h264_pixels_tab[2]= biweight_h264_pixels8x16_c;
  1950.     c->biweight_h264_pixels_tab[3]= biweight_h264_pixels8x8_c;
  1951.     c->biweight_h264_pixels_tab[4]= biweight_h264_pixels8x4_c;
  1952.     c->biweight_h264_pixels_tab[5]= biweight_h264_pixels4x8_c;
  1953.     c->biweight_h264_pixels_tab[6]= biweight_h264_pixels4x4_c;
  1954.     c->biweight_h264_pixels_tab[7]= biweight_h264_pixels4x2_c;
  1955.     c->biweight_h264_pixels_tab[8]= biweight_h264_pixels2x4_c;
  1956.     c->biweight_h264_pixels_tab[9]= biweight_h264_pixels2x2_c;
  1957.     c->draw_edges = draw_edges_c;
  1958. #ifdef CONFIG_CAVS_DECODER
  1959.     ff_cavsdsp_init(c,avctx);
  1960. #endif
  1961. #if defined(CONFIG_VC1_DECODER) || defined(CONFIG_WMV3_DECODER)
  1962.     ff_vc1dsp_init(c,avctx);
  1963. #endif
  1964. #if defined(CONFIG_WMV2_DECODER) || defined(CONFIG_VC1_DECODER) || defined(CONFIG_WMV3_DECODER)
  1965.     ff_intrax8dsp_init(c,avctx);
  1966. #endif
  1967. #if defined(CONFIG_H264_ENCODER)
  1968.     ff_h264dspenc_init(c,avctx);
  1969. #endif
  1970. #if defined(CONFIG_RV30_DECODER)
  1971.     ff_rv30dsp_init(c,avctx);
  1972. #endif
  1973. #if defined(CONFIG_RV40_DECODER)
  1974.     ff_rv40dsp_init(c,avctx);
  1975.     c->put_rv40_qpel_pixels_tab[0][15] = put_rv40_qpel16_mc33_c;
  1976.     c->avg_rv40_qpel_pixels_tab[0][15] = avg_rv40_qpel16_mc33_c;
  1977.     c->put_rv40_qpel_pixels_tab[1][15] = put_rv40_qpel8_mc33_c;
  1978.     c->avg_rv40_qpel_pixels_tab[1][15] = avg_rv40_qpel8_mc33_c;
  1979. #endif
  1980.     c->put_mspel_pixels_tab[0]= put_mspel8_mc00_c;
  1981.     c->put_mspel_pixels_tab[1]= put_mspel8_mc10_c;
  1982.     c->put_mspel_pixels_tab[2]= put_mspel8_mc20_c;
  1983.     c->put_mspel_pixels_tab[3]= put_mspel8_mc30_c;
  1984.     c->put_mspel_pixels_tab[4]= put_mspel8_mc02_c;
  1985.     c->put_mspel_pixels_tab[5]= put_mspel8_mc12_c;
  1986.     c->put_mspel_pixels_tab[6]= put_mspel8_mc22_c;
  1987.     c->put_mspel_pixels_tab[7]= put_mspel8_mc32_c;
  1988. #define SET_CMP_FUNC(name) 
  1989.     c->name[0]= name ## 16_c;
  1990.     c->name[1]= name ## 8x8_c;
  1991.     SET_CMP_FUNC(hadamard8_diff)
  1992.     c->hadamard8_diff[4]= hadamard8_intra16_c;
  1993.     SET_CMP_FUNC(dct_sad)
  1994.     SET_CMP_FUNC(dct_max)
  1995. #ifdef CONFIG_GPL
  1996.     SET_CMP_FUNC(dct264_sad)
  1997. #endif
  1998.     c->sad[0]= pix_abs16_c;
  1999.     c->sad[1]= pix_abs8_c;
  2000.     c->sse[0]= sse16_c;
  2001.     c->sse[1]= sse8_c;
  2002.     c->sse[2]= sse4_c;
  2003.     SET_CMP_FUNC(quant_psnr)
  2004.     SET_CMP_FUNC(rd)
  2005.     SET_CMP_FUNC(bit)
  2006.     c->vsad[0]= vsad16_c;
  2007.     c->vsad[4]= vsad_intra16_c;
  2008.     c->vsse[0]= vsse16_c;
  2009.     c->vsse[4]= vsse_intra16_c;
  2010.     c->nsse[0]= nsse16_c;
  2011.     c->nsse[1]= nsse8_c;
  2012. #ifdef CONFIG_SNOW_ENCODER
  2013.     c->w53[0]= w53_16_c;
  2014.     c->w53[1]= w53_8_c;
  2015.     c->w97[0]= w97_16_c;
  2016.     c->w97[1]= w97_8_c;
  2017. #endif
  2018.     c->ssd_int8_vs_int16 = ssd_int8_vs_int16_c;
  2019.     c->add_bytes= add_bytes_c;
  2020.     c->add_bytes_l2= add_bytes_l2_c;
  2021.     c->diff_bytes= diff_bytes_c;
  2022.     c->sub_hfyu_median_prediction= sub_hfyu_median_prediction_c;
  2023.     c->bswap_buf= bswap_buf;
  2024. #ifdef CONFIG_PNG_DECODER
  2025.     c->add_png_paeth_prediction= ff_add_png_paeth_prediction;
  2026. #endif
  2027.     c->h264_v_loop_filter_luma= h264_v_loop_filter_luma_c;
  2028.     c->h264_h_loop_filter_luma= h264_h_loop_filter_luma_c;
  2029.     c->h264_v_loop_filter_luma_intra= h264_v_loop_filter_luma_intra_c;
  2030.     c->h264_h_loop_filter_luma_intra= h264_h_loop_filter_luma_intra_c;
  2031.     c->h264_v_loop_filter_chroma= h264_v_loop_filter_chroma_c;
  2032.     c->h264_h_loop_filter_chroma= h264_h_loop_filter_chroma_c;
  2033.     c->h264_v_loop_filter_chroma_intra= h264_v_loop_filter_chroma_intra_c;
  2034.     c->h264_h_loop_filter_chroma_intra= h264_h_loop_filter_chroma_intra_c;
  2035.     c->h264_loop_filter_strength= NULL;
  2036.     /*if (ENABLE_ANY_H263) {
  2037.         c->h263_h_loop_filter= h263_h_loop_filter_c;
  2038.         c->h263_v_loop_filter= h263_v_loop_filter_c;
  2039.     }
  2040.     if (ENABLE_VP3_DECODER || ENABLE_THEORA_DECODER) {
  2041.         c->vp3_h_loop_filter= ff_vp3_h_loop_filter_c;
  2042.         c->vp3_v_loop_filter= ff_vp3_v_loop_filter_c;
  2043.     }*/
  2044.     c->h261_loop_filter= h261_loop_filter_c;
  2045.     c->try_8x8basis= try_8x8basis_c;
  2046.     c->add_8x8basis= add_8x8basis_c;
  2047. #ifdef CONFIG_SNOW_DECODER
  2048.     c->vertical_compose97i = ff_snow_vertical_compose97i;
  2049.     c->horizontal_compose97i = ff_snow_horizontal_compose97i;
  2050.     c->inner_add_yblock = ff_snow_inner_add_yblock;
  2051. #endif
  2052. #ifdef CONFIG_VORBIS_DECODER
  2053.     c->vorbis_inverse_coupling = vorbis_inverse_coupling;
  2054. #endif
  2055. #ifdef CONFIG_AC3_DECODER
  2056.     c->ac3_downmix = ff_ac3_downmix_c;
  2057. #endif
  2058. #ifdef CONFIG_FLAC_ENCODER
  2059.     c->flac_compute_autocorr = ff_flac_compute_autocorr;
  2060. #endif
  2061.     c->vector_fmul = vector_fmul_c;
  2062.     c->vector_fmul_reverse = vector_fmul_reverse_c;
  2063.     c->vector_fmul_add_add = ff_vector_fmul_add_add_c;
  2064.     c->vector_fmul_window = ff_vector_fmul_window_c;
  2065.     c->int32_to_float_fmul_scalar = int32_to_float_fmul_scalar_c;
  2066.     c->float_to_int16 = ff_float_to_int16_c;
  2067.     c->float_to_int16_interleave = ff_float_to_int16_interleave_c;
  2068.     c->add_int16 = add_int16_c;
  2069.     c->sub_int16 = sub_int16_c;
  2070.     c->scalarproduct_int16 = scalarproduct_int16_c;
  2071.     c->shrink[0]= ff_img_copy_plane;
  2072.     c->shrink[1]= ff_shrink22;
  2073.     c->shrink[2]= ff_shrink44;
  2074.     c->shrink[3]= ff_shrink88;
  2075.     c->prefetch= just_return;
  2076.     memset(c->put_2tap_qpel_pixels_tab, 0, sizeof(c->put_2tap_qpel_pixels_tab));
  2077.     memset(c->avg_2tap_qpel_pixels_tab, 0, sizeof(c->avg_2tap_qpel_pixels_tab));
  2078.     /*if (ENABLE_MMX)      dsputil_init_mmx   (c, avctx);
  2079.     if (ENABLE_ARM)      dsputil_init_arm   (c, avctx);
  2080.     if (ENABLE_MLIB)     dsputil_init_mlib  (c, avctx);
  2081.     if (ENABLE_VIS)      dsputil_init_vis   (c, avctx);
  2082.     if (ENABLE_ALPHA)    dsputil_init_alpha (c, avctx);
  2083.     if (ENABLE_PPC)      dsputil_init_ppc   (c, avctx);
  2084.     if (ENABLE_MMI)      dsputil_init_mmi   (c, avctx);
  2085.     if (ENABLE_SH4)      dsputil_init_sh4   (c, avctx);
  2086.     if (ENABLE_BFIN)     dsputil_init_bfin  (c, avctx);*/
  2087.     for(i=0; i<64; i++){
  2088.         if(!c->put_2tap_qpel_pixels_tab[0][i])
  2089.             c->put_2tap_qpel_pixels_tab[0][i]= c->put_h264_qpel_pixels_tab[0][i];
  2090.         if(!c->avg_2tap_qpel_pixels_tab[0][i])
  2091.             c->avg_2tap_qpel_pixels_tab[0][i]= c->avg_h264_qpel_pixels_tab[0][i];
  2092.     }
  2093.     switch(c->idct_permutation_type){
  2094.     case FF_NO_IDCT_PERM:
  2095.         for(i=0; i<64; i++)
  2096.             c->idct_permutation[i]= i;
  2097.         break;
  2098.     case FF_LIBMPEG2_IDCT_PERM:
  2099.         for(i=0; i<64; i++)
  2100.             c->idct_permutation[i]= (i & 0x38) | ((i & 6) >> 1) | ((i & 1) << 2);
  2101.         break;
  2102.     case FF_SIMPLE_IDCT_PERM:
  2103.         for(i=0; i<64; i++)
  2104.             c->idct_permutation[i]= simple_mmx_permutation[i];
  2105.         break;
  2106.     case FF_TRANSPOSE_IDCT_PERM:
  2107.         for(i=0; i<64; i++)
  2108.             c->idct_permutation[i]= ((i&7)<<3) | (i>>3);
  2109.         break;
  2110.     case FF_PARTTRANS_IDCT_PERM:
  2111.         for(i=0; i<64; i++)
  2112.             c->idct_permutation[i]= (i&0x24) | ((i&3)<<3) | ((i>>3)&3);
  2113.         break;
  2114.     case FF_SSE2_IDCT_PERM:
  2115.         for(i=0; i<64; i++)
  2116.             c->idct_permutation[i]= (i&0x38) | idct_sse2_row_perm[i&7];
  2117.         break;
  2118.     default:
  2119.         av_log(avctx, AV_LOG_ERROR, "Internal error, IDCT permutation not setn");
  2120.     }
  2121. }