• Main Page
  • Related Pages
  • Modules
  • Data Structures
  • Files
  • Examples
  • File List
  • Globals

libavcodec/cavsdsp.c

Go to the documentation of this file.
00001 /*
00002  * Chinese AVS video (AVS1-P2, JiZhun profile) decoder.
00003  *
00004  * DSP functions
00005  *
00006  * Copyright (c) 2006  Stefan Gehrer <stefan.gehrer@gmx.de>
00007  *
00008  * This file is part of Libav.
00009  *
00010  * Libav is free software; you can redistribute it and/or
00011  * modify it under the terms of the GNU Lesser General Public
00012  * License as published by the Free Software Foundation; either
00013  * version 2.1 of the License, or (at your option) any later version.
00014  *
00015  * Libav is distributed in the hope that it will be useful,
00016  * but WITHOUT ANY WARRANTY; without even the implied warranty of
00017  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
00018  * Lesser General Public License for more details.
00019  *
00020  * You should have received a copy of the GNU Lesser General Public
00021  * License along with Libav; if not, write to the Free Software
00022  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
00023  */
00024 
00025 #include <stdio.h>
00026 #include "dsputil.h"
00027 #include "cavsdsp.h"
00028 
00029 /*****************************************************************************
00030  *
00031  * in-loop deblocking filter
00032  *
00033  ****************************************************************************/
00034 
00035 #define P2 p0_p[-3*stride]
00036 #define P1 p0_p[-2*stride]
00037 #define P0 p0_p[-1*stride]
00038 #define Q0 p0_p[ 0*stride]
00039 #define Q1 p0_p[ 1*stride]
00040 #define Q2 p0_p[ 2*stride]
00041 
00042 static inline void loop_filter_l2(uint8_t *p0_p,int stride,int alpha, int beta) {
00043     int p0 = P0;
00044     int q0 = Q0;
00045 
00046     if(abs(p0-q0)<alpha && abs(P1-p0)<beta && abs(Q1-q0)<beta) {
00047         int s = p0 + q0 + 2;
00048         alpha = (alpha>>2) + 2;
00049         if(abs(P2-p0) < beta && abs(p0-q0) < alpha) {
00050             P0 = (P1 + p0 + s) >> 2;
00051             P1 = (2*P1 + s) >> 2;
00052         } else
00053             P0 = (2*P1 + s) >> 2;
00054         if(abs(Q2-q0) < beta && abs(q0-p0) < alpha) {
00055             Q0 = (Q1 + q0 + s) >> 2;
00056             Q1 = (2*Q1 + s) >> 2;
00057         } else
00058             Q0 = (2*Q1 + s) >> 2;
00059     }
00060 }
00061 
00062 static inline void loop_filter_l1(uint8_t *p0_p, int stride, int alpha, int beta, int tc) {
00063     int p0 = P0;
00064     int q0 = Q0;
00065 
00066     if(abs(p0-q0)<alpha && abs(P1-p0)<beta && abs(Q1-q0)<beta) {
00067         int delta = av_clip(((q0-p0)*3+P1-Q1+4)>>3,-tc, tc);
00068         P0 = av_clip_uint8(p0+delta);
00069         Q0 = av_clip_uint8(q0-delta);
00070         if(abs(P2-p0)<beta) {
00071             delta = av_clip(((P0-P1)*3+P2-Q0+4)>>3, -tc, tc);
00072             P1 = av_clip_uint8(P1+delta);
00073         }
00074         if(abs(Q2-q0)<beta) {
00075             delta = av_clip(((Q1-Q0)*3+P0-Q2+4)>>3, -tc, tc);
00076             Q1 = av_clip_uint8(Q1-delta);
00077         }
00078     }
00079 }
00080 
00081 static inline void loop_filter_c2(uint8_t *p0_p,int stride,int alpha, int beta) {
00082     int p0 = P0;
00083     int q0 = Q0;
00084 
00085     if(abs(p0-q0)<alpha && abs(P1-p0)<beta && abs(Q1-q0)<beta) {
00086         int s = p0 + q0 + 2;
00087         alpha = (alpha>>2) + 2;
00088         if(abs(P2-p0) < beta && abs(p0-q0) < alpha) {
00089             P0 = (P1 + p0 + s) >> 2;
00090         } else
00091             P0 = (2*P1 + s) >> 2;
00092         if(abs(Q2-q0) < beta && abs(q0-p0) < alpha) {
00093             Q0 = (Q1 + q0 + s) >> 2;
00094         } else
00095             Q0 = (2*Q1 + s) >> 2;
00096     }
00097 }
00098 
00099 static inline void loop_filter_c1(uint8_t *p0_p,int stride,int alpha, int beta,
00100                                   int tc) {
00101     if(abs(P0-Q0)<alpha && abs(P1-P0)<beta && abs(Q1-Q0)<beta) {
00102         int delta = av_clip(((Q0-P0)*3+P1-Q1+4)>>3, -tc, tc);
00103         P0 = av_clip_uint8(P0+delta);
00104         Q0 = av_clip_uint8(Q0-delta);
00105     }
00106 }
00107 
00108 #undef P0
00109 #undef P1
00110 #undef P2
00111 #undef Q0
00112 #undef Q1
00113 #undef Q2
00114 
00115 static void cavs_filter_lv_c(uint8_t *d, int stride, int alpha, int beta, int tc,
00116                            int bs1, int bs2) {
00117     int i;
00118     if(bs1==2)
00119         for(i=0;i<16;i++)
00120             loop_filter_l2(d + i*stride,1,alpha,beta);
00121     else {
00122         if(bs1)
00123             for(i=0;i<8;i++)
00124                 loop_filter_l1(d + i*stride,1,alpha,beta,tc);
00125         if (bs2)
00126             for(i=8;i<16;i++)
00127                 loop_filter_l1(d + i*stride,1,alpha,beta,tc);
00128     }
00129 }
00130 
00131 static void cavs_filter_lh_c(uint8_t *d, int stride, int alpha, int beta, int tc,
00132                            int bs1, int bs2) {
00133     int i;
00134     if(bs1==2)
00135         for(i=0;i<16;i++)
00136             loop_filter_l2(d + i,stride,alpha,beta);
00137     else {
00138         if(bs1)
00139             for(i=0;i<8;i++)
00140                 loop_filter_l1(d + i,stride,alpha,beta,tc);
00141         if (bs2)
00142             for(i=8;i<16;i++)
00143                 loop_filter_l1(d + i,stride,alpha,beta,tc);
00144     }
00145 }
00146 
00147 static void cavs_filter_cv_c(uint8_t *d, int stride, int alpha, int beta, int tc,
00148                            int bs1, int bs2) {
00149     int i;
00150     if(bs1==2)
00151         for(i=0;i<8;i++)
00152             loop_filter_c2(d + i*stride,1,alpha,beta);
00153     else {
00154         if(bs1)
00155             for(i=0;i<4;i++)
00156                 loop_filter_c1(d + i*stride,1,alpha,beta,tc);
00157         if (bs2)
00158             for(i=4;i<8;i++)
00159                 loop_filter_c1(d + i*stride,1,alpha,beta,tc);
00160     }
00161 }
00162 
00163 static void cavs_filter_ch_c(uint8_t *d, int stride, int alpha, int beta, int tc,
00164                            int bs1, int bs2) {
00165     int i;
00166     if(bs1==2)
00167         for(i=0;i<8;i++)
00168             loop_filter_c2(d + i,stride,alpha,beta);
00169     else {
00170         if(bs1)
00171             for(i=0;i<4;i++)
00172                 loop_filter_c1(d + i,stride,alpha,beta,tc);
00173         if (bs2)
00174             for(i=4;i<8;i++)
00175                 loop_filter_c1(d + i,stride,alpha,beta,tc);
00176     }
00177 }
00178 
00179 /*****************************************************************************
00180  *
00181  * inverse transform
00182  *
00183  ****************************************************************************/
00184 
00185 static void cavs_idct8_add_c(uint8_t *dst, DCTELEM *block, int stride) {
00186     int i;
00187     DCTELEM (*src)[8] = (DCTELEM(*)[8])block;
00188     uint8_t *cm = ff_cropTbl + MAX_NEG_CROP;
00189 
00190     src[0][0] += 8;
00191 
00192     for( i = 0; i < 8; i++ ) {
00193         const int a0 =  3*src[i][1] - (src[i][7]<<1);
00194         const int a1 =  3*src[i][3] + (src[i][5]<<1);
00195         const int a2 =  (src[i][3]<<1) - 3*src[i][5];
00196         const int a3 =  (src[i][1]<<1) + 3*src[i][7];
00197 
00198         const int b4 = ((a0 + a1 + a3)<<1) + a1;
00199         const int b5 = ((a0 - a1 + a2)<<1) + a0;
00200         const int b6 = ((a3 - a2 - a1)<<1) + a3;
00201         const int b7 = ((a0 - a2 - a3)<<1) - a2;
00202 
00203         const int a7 = (src[i][2]<<2) - 10*src[i][6];
00204         const int a6 = (src[i][6]<<2) + 10*src[i][2];
00205         const int a5 = ((src[i][0] - src[i][4]) << 3) + 4;
00206         const int a4 = ((src[i][0] + src[i][4]) << 3) + 4;
00207 
00208         const int b0 = a4 + a6;
00209         const int b1 = a5 + a7;
00210         const int b2 = a5 - a7;
00211         const int b3 = a4 - a6;
00212 
00213         src[i][0] = (b0 + b4) >> 3;
00214         src[i][1] = (b1 + b5) >> 3;
00215         src[i][2] = (b2 + b6) >> 3;
00216         src[i][3] = (b3 + b7) >> 3;
00217         src[i][4] = (b3 - b7) >> 3;
00218         src[i][5] = (b2 - b6) >> 3;
00219         src[i][6] = (b1 - b5) >> 3;
00220         src[i][7] = (b0 - b4) >> 3;
00221     }
00222     for( i = 0; i < 8; i++ ) {
00223         const int a0 =  3*src[1][i] - (src[7][i]<<1);
00224         const int a1 =  3*src[3][i] + (src[5][i]<<1);
00225         const int a2 =  (src[3][i]<<1) - 3*src[5][i];
00226         const int a3 =  (src[1][i]<<1) + 3*src[7][i];
00227 
00228         const int b4 = ((a0 + a1 + a3)<<1) + a1;
00229         const int b5 = ((a0 - a1 + a2)<<1) + a0;
00230         const int b6 = ((a3 - a2 - a1)<<1) + a3;
00231         const int b7 = ((a0 - a2 - a3)<<1) - a2;
00232 
00233         const int a7 = (src[2][i]<<2) - 10*src[6][i];
00234         const int a6 = (src[6][i]<<2) + 10*src[2][i];
00235         const int a5 = (src[0][i] - src[4][i]) << 3;
00236         const int a4 = (src[0][i] + src[4][i]) << 3;
00237 
00238         const int b0 = a4 + a6;
00239         const int b1 = a5 + a7;
00240         const int b2 = a5 - a7;
00241         const int b3 = a4 - a6;
00242 
00243         dst[i + 0*stride] = cm[ dst[i + 0*stride] + ((b0 + b4) >> 7)];
00244         dst[i + 1*stride] = cm[ dst[i + 1*stride] + ((b1 + b5) >> 7)];
00245         dst[i + 2*stride] = cm[ dst[i + 2*stride] + ((b2 + b6) >> 7)];
00246         dst[i + 3*stride] = cm[ dst[i + 3*stride] + ((b3 + b7) >> 7)];
00247         dst[i + 4*stride] = cm[ dst[i + 4*stride] + ((b3 - b7) >> 7)];
00248         dst[i + 5*stride] = cm[ dst[i + 5*stride] + ((b2 - b6) >> 7)];
00249         dst[i + 6*stride] = cm[ dst[i + 6*stride] + ((b1 - b5) >> 7)];
00250         dst[i + 7*stride] = cm[ dst[i + 7*stride] + ((b0 - b4) >> 7)];
00251     }
00252 }
00253 
00254 /*****************************************************************************
00255  *
00256  * motion compensation
00257  *
00258  ****************************************************************************/
00259 
00260 #define CAVS_SUBPIX(OPNAME, OP, NAME, A, B, C, D, E, F) \
00261 static void OPNAME ## cavs_filt8_h_ ## NAME(uint8_t *dst, uint8_t *src, int dstStride, int srcStride){\
00262     const int h=8;\
00263     uint8_t *cm = ff_cropTbl + MAX_NEG_CROP;\
00264     int i;\
00265     for(i=0; i<h; i++)\
00266     {\
00267         OP(dst[0], A*src[-2] + B*src[-1] + C*src[0] + D*src[1] + E*src[2] + F*src[3]);\
00268         OP(dst[1], A*src[-1] + B*src[ 0] + C*src[1] + D*src[2] + E*src[3] + F*src[4]);\
00269         OP(dst[2], A*src[ 0] + B*src[ 1] + C*src[2] + D*src[3] + E*src[4] + F*src[5]);\
00270         OP(dst[3], A*src[ 1] + B*src[ 2] + C*src[3] + D*src[4] + E*src[5] + F*src[6]);\
00271         OP(dst[4], A*src[ 2] + B*src[ 3] + C*src[4] + D*src[5] + E*src[6] + F*src[7]);\
00272         OP(dst[5], A*src[ 3] + B*src[ 4] + C*src[5] + D*src[6] + E*src[7] + F*src[8]);\
00273         OP(dst[6], A*src[ 4] + B*src[ 5] + C*src[6] + D*src[7] + E*src[8] + F*src[9]);\
00274         OP(dst[7], A*src[ 5] + B*src[ 6] + C*src[7] + D*src[8] + E*src[9] + F*src[10]);\
00275         dst+=dstStride;\
00276         src+=srcStride;\
00277     }\
00278 }\
00279 \
00280 static void OPNAME ## cavs_filt8_v_  ## NAME(uint8_t *dst, uint8_t *src, int dstStride, int srcStride){\
00281     const int w=8;\
00282     uint8_t *cm = ff_cropTbl + MAX_NEG_CROP;\
00283     int i;\
00284     for(i=0; i<w; i++)\
00285     {\
00286         const int srcB= src[-2*srcStride];\
00287         const int srcA= src[-1*srcStride];\
00288         const int src0= src[0 *srcStride];\
00289         const int src1= src[1 *srcStride];\
00290         const int src2= src[2 *srcStride];\
00291         const int src3= src[3 *srcStride];\
00292         const int src4= src[4 *srcStride];\
00293         const int src5= src[5 *srcStride];\
00294         const int src6= src[6 *srcStride];\
00295         const int src7= src[7 *srcStride];\
00296         const int src8= src[8 *srcStride];\
00297         const int src9= src[9 *srcStride];\
00298         const int src10= src[10 *srcStride];\
00299         OP(dst[0*dstStride], A*srcB + B*srcA + C*src0 + D*src1 + E*src2 + F*src3);\
00300         OP(dst[1*dstStride], A*srcA + B*src0 + C*src1 + D*src2 + E*src3 + F*src4);\
00301         OP(dst[2*dstStride], A*src0 + B*src1 + C*src2 + D*src3 + E*src4 + F*src5);\
00302         OP(dst[3*dstStride], A*src1 + B*src2 + C*src3 + D*src4 + E*src5 + F*src6);\
00303         OP(dst[4*dstStride], A*src2 + B*src3 + C*src4 + D*src5 + E*src6 + F*src7);\
00304         OP(dst[5*dstStride], A*src3 + B*src4 + C*src5 + D*src6 + E*src7 + F*src8);\
00305         OP(dst[6*dstStride], A*src4 + B*src5 + C*src6 + D*src7 + E*src8 + F*src9);\
00306         OP(dst[7*dstStride], A*src5 + B*src6 + C*src7 + D*src8 + E*src9 + F*src10);\
00307         dst++;\
00308         src++;\
00309     }\
00310 }\
00311 \
00312 static void OPNAME ## cavs_filt16_v_ ## NAME(uint8_t *dst, uint8_t *src, int dstStride, int srcStride){\
00313     OPNAME ## cavs_filt8_v_ ## NAME(dst  , src  , dstStride, srcStride);\
00314     OPNAME ## cavs_filt8_v_ ## NAME(dst+8, src+8, dstStride, srcStride);\
00315     src += 8*srcStride;\
00316     dst += 8*dstStride;\
00317     OPNAME ## cavs_filt8_v_ ## NAME(dst  , src  , dstStride, srcStride);\
00318     OPNAME ## cavs_filt8_v_ ## NAME(dst+8, src+8, dstStride, srcStride);\
00319 }\
00320 \
00321 static void OPNAME ## cavs_filt16_h_ ## NAME(uint8_t *dst, uint8_t *src, int dstStride, int srcStride){\
00322     OPNAME ## cavs_filt8_h_ ## NAME(dst  , src  , dstStride, srcStride);\
00323     OPNAME ## cavs_filt8_h_ ## NAME(dst+8, src+8, dstStride, srcStride);\
00324     src += 8*srcStride;\
00325     dst += 8*dstStride;\
00326     OPNAME ## cavs_filt8_h_ ## NAME(dst  , src  , dstStride, srcStride);\
00327     OPNAME ## cavs_filt8_h_ ## NAME(dst+8, src+8, dstStride, srcStride);\
00328 }\
00329 
00330 #define CAVS_SUBPIX_HV(OPNAME, OP, NAME, AH, BH, CH, DH, EH, FH, AV, BV, CV, DV, EV, FV, FULL) \
00331 static void OPNAME ## cavs_filt8_hv_ ## NAME(uint8_t *dst, uint8_t *src1, uint8_t *src2, int dstStride, int srcStride){\
00332     int16_t temp[8*(8+5)];\
00333     int16_t *tmp = temp;\
00334     const int h=8;\
00335     const int w=8;\
00336     uint8_t *cm = ff_cropTbl + MAX_NEG_CROP;\
00337     int i;\
00338     src1 -= 2*srcStride;\
00339     for(i=0; i<h+5; i++)\
00340     {\
00341         tmp[0]= AH*src1[-2] + BH*src1[-1] + CH*src1[0] + DH*src1[1] + EH*src1[2] + FH*src1[3];\
00342         tmp[1]= AH*src1[-1] + BH*src1[ 0] + CH*src1[1] + DH*src1[2] + EH*src1[3] + FH*src1[4];\
00343         tmp[2]= AH*src1[ 0] + BH*src1[ 1] + CH*src1[2] + DH*src1[3] + EH*src1[4] + FH*src1[5];\
00344         tmp[3]= AH*src1[ 1] + BH*src1[ 2] + CH*src1[3] + DH*src1[4] + EH*src1[5] + FH*src1[6];\
00345         tmp[4]= AH*src1[ 2] + BH*src1[ 3] + CH*src1[4] + DH*src1[5] + EH*src1[6] + FH*src1[7];\
00346         tmp[5]= AH*src1[ 3] + BH*src1[ 4] + CH*src1[5] + DH*src1[6] + EH*src1[7] + FH*src1[8];\
00347         tmp[6]= AH*src1[ 4] + BH*src1[ 5] + CH*src1[6] + DH*src1[7] + EH*src1[8] + FH*src1[9];\
00348         tmp[7]= AH*src1[ 5] + BH*src1[ 6] + CH*src1[7] + DH*src1[8] + EH*src1[9] + FH*src1[10];\
00349         tmp+=8;\
00350         src1+=srcStride;\
00351     }\
00352     if(FULL) {\
00353       tmp = temp+8*2;                           \
00354       for(i=0; i<w; i++)                        \
00355         {                                       \
00356           const int tmpB= tmp[-2*8];    \
00357           const int tmpA= tmp[-1*8];    \
00358           const int tmp0= tmp[0 *8];    \
00359           const int tmp1= tmp[1 *8];    \
00360           const int tmp2= tmp[2 *8];    \
00361           const int tmp3= tmp[3 *8];    \
00362           const int tmp4= tmp[4 *8];    \
00363           const int tmp5= tmp[5 *8];    \
00364           const int tmp6= tmp[6 *8];    \
00365           const int tmp7= tmp[7 *8];    \
00366           const int tmp8= tmp[8 *8];    \
00367           const int tmp9= tmp[9 *8];    \
00368           const int tmp10=tmp[10*8];                            \
00369           OP(dst[0*dstStride], AV*tmpB + BV*tmpA + CV*tmp0 + DV*tmp1 + EV*tmp2 + FV*tmp3 + 64*src2[0*srcStride]); \
00370           OP(dst[1*dstStride], AV*tmpA + BV*tmp0 + CV*tmp1 + DV*tmp2 + EV*tmp3 + FV*tmp4 + 64*src2[1*srcStride]); \
00371           OP(dst[2*dstStride], AV*tmp0 + BV*tmp1 + CV*tmp2 + DV*tmp3 + EV*tmp4 + FV*tmp5 + 64*src2[2*srcStride]); \
00372           OP(dst[3*dstStride], AV*tmp1 + BV*tmp2 + CV*tmp3 + DV*tmp4 + EV*tmp5 + FV*tmp6 + 64*src2[3*srcStride]); \
00373           OP(dst[4*dstStride], AV*tmp2 + BV*tmp3 + CV*tmp4 + DV*tmp5 + EV*tmp6 + FV*tmp7 + 64*src2[4*srcStride]); \
00374           OP(dst[5*dstStride], AV*tmp3 + BV*tmp4 + CV*tmp5 + DV*tmp6 + EV*tmp7 + FV*tmp8 + 64*src2[5*srcStride]); \
00375           OP(dst[6*dstStride], AV*tmp4 + BV*tmp5 + CV*tmp6 + DV*tmp7 + EV*tmp8 + FV*tmp9 + 64*src2[6*srcStride]); \
00376           OP(dst[7*dstStride], AV*tmp5 + BV*tmp6 + CV*tmp7 + DV*tmp8 + EV*tmp9 + FV*tmp10 + 64*src2[7*srcStride]); \
00377           dst++;                                                        \
00378           tmp++;                                                        \
00379           src2++;                                                       \
00380         }                                                               \
00381     } else {\
00382       tmp = temp+8*2;                           \
00383       for(i=0; i<w; i++)                        \
00384         {                                       \
00385           const int tmpB= tmp[-2*8];    \
00386           const int tmpA= tmp[-1*8];    \
00387           const int tmp0= tmp[0 *8];    \
00388           const int tmp1= tmp[1 *8];    \
00389           const int tmp2= tmp[2 *8];    \
00390           const int tmp3= tmp[3 *8];    \
00391           const int tmp4= tmp[4 *8];    \
00392           const int tmp5= tmp[5 *8];    \
00393           const int tmp6= tmp[6 *8];    \
00394           const int tmp7= tmp[7 *8];    \
00395           const int tmp8= tmp[8 *8];    \
00396           const int tmp9= tmp[9 *8];    \
00397           const int tmp10=tmp[10*8];                            \
00398           OP(dst[0*dstStride], AV*tmpB + BV*tmpA + CV*tmp0 + DV*tmp1 + EV*tmp2 + FV*tmp3); \
00399           OP(dst[1*dstStride], AV*tmpA + BV*tmp0 + CV*tmp1 + DV*tmp2 + EV*tmp3 + FV*tmp4); \
00400           OP(dst[2*dstStride], AV*tmp0 + BV*tmp1 + CV*tmp2 + DV*tmp3 + EV*tmp4 + FV*tmp5); \
00401           OP(dst[3*dstStride], AV*tmp1 + BV*tmp2 + CV*tmp3 + DV*tmp4 + EV*tmp5 + FV*tmp6); \
00402           OP(dst[4*dstStride], AV*tmp2 + BV*tmp3 + CV*tmp4 + DV*tmp5 + EV*tmp6 + FV*tmp7); \
00403           OP(dst[5*dstStride], AV*tmp3 + BV*tmp4 + CV*tmp5 + DV*tmp6 + EV*tmp7 + FV*tmp8); \
00404           OP(dst[6*dstStride], AV*tmp4 + BV*tmp5 + CV*tmp6 + DV*tmp7 + EV*tmp8 + FV*tmp9); \
00405           OP(dst[7*dstStride], AV*tmp5 + BV*tmp6 + CV*tmp7 + DV*tmp8 + EV*tmp9 + FV*tmp10); \
00406           dst++;                                                        \
00407           tmp++;                                                        \
00408         }                                                               \
00409     }\
00410 }\
00411 \
00412 static void OPNAME ## cavs_filt16_hv_ ## NAME(uint8_t *dst, uint8_t *src1, uint8_t *src2, int dstStride, int srcStride){ \
00413     OPNAME ## cavs_filt8_hv_ ## NAME(dst  , src1,   src2  , dstStride, srcStride); \
00414     OPNAME ## cavs_filt8_hv_ ## NAME(dst+8, src1+8, src2+8, dstStride, srcStride); \
00415     src1 += 8*srcStride;\
00416     src2 += 8*srcStride;\
00417     dst += 8*dstStride;\
00418     OPNAME ## cavs_filt8_hv_ ## NAME(dst  , src1,   src2  , dstStride, srcStride); \
00419     OPNAME ## cavs_filt8_hv_ ## NAME(dst+8, src1+8, src2+8, dstStride, srcStride); \
00420 }\
00421 
00422 #define CAVS_MC(OPNAME, SIZE) \
00423 static void ff_ ## OPNAME ## cavs_qpel ## SIZE ## _mc10_c(uint8_t *dst, uint8_t *src, int stride){\
00424     OPNAME ## cavs_filt ## SIZE ## _h_qpel_l(dst, src, stride, stride);\
00425 }\
00426 \
00427 static void ff_ ## OPNAME ## cavs_qpel ## SIZE ## _mc20_c(uint8_t *dst, uint8_t *src, int stride){\
00428     OPNAME ## cavs_filt ## SIZE ## _h_hpel(dst, src, stride, stride);\
00429 }\
00430 \
00431 static void ff_ ## OPNAME ## cavs_qpel ## SIZE ## _mc30_c(uint8_t *dst, uint8_t *src, int stride){\
00432     OPNAME ## cavs_filt ## SIZE ## _h_qpel_r(dst, src, stride, stride);\
00433 }\
00434 \
00435 static void ff_ ## OPNAME ## cavs_qpel ## SIZE ## _mc01_c(uint8_t *dst, uint8_t *src, int stride){\
00436     OPNAME ## cavs_filt ## SIZE ## _v_qpel_l(dst, src, stride, stride);\
00437 }\
00438 \
00439 static void ff_ ## OPNAME ## cavs_qpel ## SIZE ## _mc02_c(uint8_t *dst, uint8_t *src, int stride){\
00440     OPNAME ## cavs_filt ## SIZE ## _v_hpel(dst, src, stride, stride);\
00441 }\
00442 \
00443 static void ff_ ## OPNAME ## cavs_qpel ## SIZE ## _mc03_c(uint8_t *dst, uint8_t *src, int stride){\
00444     OPNAME ## cavs_filt ## SIZE ## _v_qpel_r(dst, src, stride, stride);\
00445 }\
00446 \
00447 static void ff_ ## OPNAME ## cavs_qpel ## SIZE ## _mc22_c(uint8_t *dst, uint8_t *src, int stride){\
00448   OPNAME ## cavs_filt ## SIZE ## _hv_jj(dst, src, NULL, stride, stride); \
00449 }\
00450 \
00451 static void ff_ ## OPNAME ## cavs_qpel ## SIZE ## _mc11_c(uint8_t *dst, uint8_t *src, int stride){\
00452   OPNAME ## cavs_filt ## SIZE ## _hv_egpr(dst, src, src, stride, stride); \
00453 }\
00454 \
00455 static void ff_ ## OPNAME ## cavs_qpel ## SIZE ## _mc13_c(uint8_t *dst, uint8_t *src, int stride){\
00456   OPNAME ## cavs_filt ## SIZE ## _hv_egpr(dst, src, src+stride, stride, stride); \
00457 }\
00458 \
00459 static void ff_ ## OPNAME ## cavs_qpel ## SIZE ## _mc31_c(uint8_t *dst, uint8_t *src, int stride){\
00460   OPNAME ## cavs_filt ## SIZE ## _hv_egpr(dst, src, src+1, stride, stride); \
00461 }\
00462 \
00463 static void ff_ ## OPNAME ## cavs_qpel ## SIZE ## _mc33_c(uint8_t *dst, uint8_t *src, int stride){\
00464   OPNAME ## cavs_filt ## SIZE ## _hv_egpr(dst, src, src+stride+1,stride, stride); \
00465 }\
00466 \
00467 static void ff_ ## OPNAME ## cavs_qpel ## SIZE ## _mc21_c(uint8_t *dst, uint8_t *src, int stride){\
00468   OPNAME ## cavs_filt ## SIZE ## _hv_ff(dst, src, src+stride+1,stride, stride); \
00469 }\
00470 \
00471 static void ff_ ## OPNAME ## cavs_qpel ## SIZE ## _mc12_c(uint8_t *dst, uint8_t *src, int stride){\
00472   OPNAME ## cavs_filt ## SIZE ## _hv_ii(dst, src, src+stride+1,stride, stride); \
00473 }\
00474 \
00475 static void ff_ ## OPNAME ## cavs_qpel ## SIZE ## _mc32_c(uint8_t *dst, uint8_t *src, int stride){\
00476   OPNAME ## cavs_filt ## SIZE ## _hv_kk(dst, src, src+stride+1,stride, stride); \
00477 }\
00478 \
00479 static void ff_ ## OPNAME ## cavs_qpel ## SIZE ## _mc23_c(uint8_t *dst, uint8_t *src, int stride){\
00480   OPNAME ## cavs_filt ## SIZE ## _hv_qq(dst, src, src+stride+1,stride, stride); \
00481 }\
00482 
00483 #define op_put1(a, b)  a = cm[((b)+4)>>3]
00484 #define op_put2(a, b)  a = cm[((b)+64)>>7]
00485 #define op_put3(a, b)  a = cm[((b)+32)>>6]
00486 #define op_put4(a, b)  a = cm[((b)+512)>>10]
00487 #define op_avg1(a, b)  a = ((a)+cm[((b)+4)>>3]   +1)>>1
00488 #define op_avg2(a, b)  a = ((a)+cm[((b)+64)>>7]  +1)>>1
00489 #define op_avg3(a, b)  a = ((a)+cm[((b)+32)>>6]  +1)>>1
00490 #define op_avg4(a, b)  a = ((a)+cm[((b)+512)>>10]+1)>>1
00491 CAVS_SUBPIX(put_   , op_put1, hpel,    0, -1,  5,  5, -1,  0)
00492 CAVS_SUBPIX(put_   , op_put2, qpel_l, -1, -2, 96, 42, -7,  0)
00493 CAVS_SUBPIX(put_   , op_put2, qpel_r,  0, -7, 42, 96, -2, -1)
00494 CAVS_SUBPIX_HV(put_, op_put3, jj,      0, -1,  5,  5, -1,  0,  0, -1,  5,  5, -1, 0, 0)
00495 CAVS_SUBPIX_HV(put_, op_put4, ff,      0, -1,  5,  5, -1,  0, -1, -2, 96, 42, -7, 0, 0)
00496 CAVS_SUBPIX_HV(put_, op_put4, ii,     -1, -2, 96, 42, -7,  0,  0, -1,  5,  5, -1, 0, 0)
00497 CAVS_SUBPIX_HV(put_, op_put4, kk,      0, -7, 42, 96, -2, -1,  0, -1,  5,  5, -1, 0, 0)
00498 CAVS_SUBPIX_HV(put_, op_put4, qq,      0, -1,  5,  5, -1,  0,  0, -7, 42, 96, -2,-1, 0)
00499 CAVS_SUBPIX_HV(put_, op_put2, egpr,    0, -1,  5,  5, -1,  0,  0, -1,  5,  5, -1, 0, 1)
00500 CAVS_SUBPIX(avg_   , op_avg1, hpel,    0, -1,  5,  5, -1,  0)
00501 CAVS_SUBPIX(avg_   , op_avg2, qpel_l, -1, -2, 96, 42, -7,  0)
00502 CAVS_SUBPIX(avg_   , op_avg2, qpel_r,  0, -7, 42, 96, -2, -1)
00503 CAVS_SUBPIX_HV(avg_, op_avg3, jj,      0, -1,  5,  5, -1,  0,  0, -1,  5,  5, -1, 0, 0)
00504 CAVS_SUBPIX_HV(avg_, op_avg4, ff,      0, -1,  5,  5, -1,  0, -1, -2, 96, 42, -7, 0, 0)
00505 CAVS_SUBPIX_HV(avg_, op_avg4, ii,     -1, -2, 96, 42, -7,  0,  0, -1,  5,  5, -1, 0, 0)
00506 CAVS_SUBPIX_HV(avg_, op_avg4, kk,      0, -7, 42, 96, -2, -1,  0, -1,  5,  5, -1, 0, 0)
00507 CAVS_SUBPIX_HV(avg_, op_avg4, qq,      0, -1,  5,  5, -1,  0,  0, -7, 42, 96, -2,-1, 0)
00508 CAVS_SUBPIX_HV(avg_, op_avg2, egpr,    0, -1,  5,  5, -1,  0,  0, -1,  5,  5, -1, 0, 1)
00509 CAVS_MC(put_, 8)
00510 CAVS_MC(put_, 16)
00511 CAVS_MC(avg_, 8)
00512 CAVS_MC(avg_, 16)
00513 
00514 #define ff_put_cavs_qpel8_mc00_c  ff_put_pixels8x8_c
00515 #define ff_avg_cavs_qpel8_mc00_c  ff_avg_pixels8x8_c
00516 #define ff_put_cavs_qpel16_mc00_c ff_put_pixels16x16_c
00517 #define ff_avg_cavs_qpel16_mc00_c ff_avg_pixels16x16_c
00518 
00519 av_cold void ff_cavsdsp_init(CAVSDSPContext* c, AVCodecContext *avctx) {
00520 #define dspfunc(PFX, IDX, NUM) \
00521     c->PFX ## _pixels_tab[IDX][ 0] = ff_ ## PFX ## NUM ## _mc00_c; \
00522     c->PFX ## _pixels_tab[IDX][ 1] = ff_ ## PFX ## NUM ## _mc10_c; \
00523     c->PFX ## _pixels_tab[IDX][ 2] = ff_ ## PFX ## NUM ## _mc20_c; \
00524     c->PFX ## _pixels_tab[IDX][ 3] = ff_ ## PFX ## NUM ## _mc30_c; \
00525     c->PFX ## _pixels_tab[IDX][ 4] = ff_ ## PFX ## NUM ## _mc01_c; \
00526     c->PFX ## _pixels_tab[IDX][ 5] = ff_ ## PFX ## NUM ## _mc11_c; \
00527     c->PFX ## _pixels_tab[IDX][ 6] = ff_ ## PFX ## NUM ## _mc21_c; \
00528     c->PFX ## _pixels_tab[IDX][ 7] = ff_ ## PFX ## NUM ## _mc31_c; \
00529     c->PFX ## _pixels_tab[IDX][ 8] = ff_ ## PFX ## NUM ## _mc02_c; \
00530     c->PFX ## _pixels_tab[IDX][ 9] = ff_ ## PFX ## NUM ## _mc12_c; \
00531     c->PFX ## _pixels_tab[IDX][10] = ff_ ## PFX ## NUM ## _mc22_c; \
00532     c->PFX ## _pixels_tab[IDX][11] = ff_ ## PFX ## NUM ## _mc32_c; \
00533     c->PFX ## _pixels_tab[IDX][12] = ff_ ## PFX ## NUM ## _mc03_c; \
00534     c->PFX ## _pixels_tab[IDX][13] = ff_ ## PFX ## NUM ## _mc13_c; \
00535     c->PFX ## _pixels_tab[IDX][14] = ff_ ## PFX ## NUM ## _mc23_c; \
00536     c->PFX ## _pixels_tab[IDX][15] = ff_ ## PFX ## NUM ## _mc33_c
00537     dspfunc(put_cavs_qpel, 0, 16);
00538     dspfunc(put_cavs_qpel, 1, 8);
00539     dspfunc(avg_cavs_qpel, 0, 16);
00540     dspfunc(avg_cavs_qpel, 1, 8);
00541     c->cavs_filter_lv = cavs_filter_lv_c;
00542     c->cavs_filter_lh = cavs_filter_lh_c;
00543     c->cavs_filter_cv = cavs_filter_cv_c;
00544     c->cavs_filter_ch = cavs_filter_ch_c;
00545     c->cavs_idct8_add = cavs_idct8_add_c;
00546 
00547     if (HAVE_MMX) ff_cavsdsp_init_mmx(c, avctx);
00548 }
Generated on Thu Jul 11 2013 15:38:18 for Libav by doxygen 1.7.1