]> git.sesse.net Git - ffmpeg/blobdiff - postproc/swscale.c
AltiVec hScale, all size patch by (Romain Dolbeau <dolbeaur at club-internet dot...
[ffmpeg] / postproc / swscale.c
index ee230b7a26d4ec56d502f2a824338172c94a3afd..dd10521582ceaeee8204cc500b519b645a0cd474 100644 (file)
@@ -18,7 +18,7 @@
 
 /*
   supported Input formats: YV12, I420/IYUV, YUY2, UYVY, BGR32, BGR24, BGR16, BGR15, RGB32, RGB24, Y8/Y800, YVU9/IF09
-  supported output formats: YV12, I420/IYUV, YUY2, {BGR,RGB}{1,4,8,15,16,24,32}, Y8/Y800, YVU9/IF09
+  supported output formats: YV12, I420/IYUV, YUY2, UYVY, {BGR,RGB}{1,4,8,15,16,24,32}, Y8/Y800, YVU9/IF09
   {BGR,RGB}{1,4,8,15,16} support dithering
   
   unscaled special converters (YV12=I420=IYUV, Y800=Y8)
@@ -61,6 +61,9 @@ untested special converters
 #else
 #include <stdlib.h>
 #endif
+#ifdef HAVE_ALTIVEC_H
+#include <altivec.h>
+#endif
 #include "swscale.h"
 #include "swscale_internal.h"
 #include "../cpudetect.h"
@@ -68,14 +71,6 @@ untested special converters
 #include "../libvo/img_format.h"
 #include "rgb2rgb.h"
 #include "../libvo/fastmemcpy.h"
-#include "../mp_msg.h"
-
-#define MSG_WARN(args...) mp_msg(MSGT_SWS,MSGL_WARN, ##args )
-#define MSG_FATAL(args...) mp_msg(MSGT_SWS,MSGL_FATAL, ##args )
-#define MSG_ERR(args...) mp_msg(MSGT_SWS,MSGL_ERR, ##args )
-#define MSG_V(args...) mp_msg(MSGT_SWS,MSGL_V, ##args )
-#define MSG_DBG2(args...) mp_msg(MSGT_SWS,MSGL_DBG2, ##args )
-#define MSG_INFO(args...) mp_msg(MSGT_SWS,MSGL_INFO, ##args )
 
 #undef MOVNTQ
 #undef PAVGB
@@ -115,7 +110,7 @@ untested special converters
                        || (x)==IMGFMT_RGB32|| (x)==IMGFMT_RGB24\
                        || (x)==IMGFMT_Y800 || (x)==IMGFMT_YVU9\
                        || (x)==IMGFMT_444P || (x)==IMGFMT_422P || (x)==IMGFMT_411P)
-#define isSupportedOut(x) ((x)==IMGFMT_YV12 || (x)==IMGFMT_YUY2\
+#define isSupportedOut(x) ((x)==IMGFMT_YV12 || (x)==IMGFMT_YUY2 || (x)==IMGFMT_UYVY\
                        || (x)==IMGFMT_444P || (x)==IMGFMT_422P || (x)==IMGFMT_411P\
                        || isRGB(x) || isBGR(x)\
                        || (x)==IMGFMT_Y800 || (x)==IMGFMT_YVU9)
@@ -132,7 +127,6 @@ untested special converters
 #define RV ((int)( 0.439*(1<<RGB2YUV_SHIFT)+0.5))
 #define RU ((int)(-0.148*(1<<RGB2YUV_SHIFT)+0.5))
 
-extern int verbose; // defined in mplayer.c
 extern const int32_t Inverse_Table_6_9[8][4];
 
 /*
@@ -155,10 +149,6 @@ write special BGR->BGR scaler
 #define MAX(a,b) ((a) < (b) ? (b) : (a))
 
 #ifdef ARCH_X86
-#define CAN_COMPILE_X86_ASM
-#endif
-
-#ifdef CAN_COMPILE_X86_ASM
 static uint64_t __attribute__((aligned(8))) bF8=       0xF8F8F8F8F8F8F8F8LL;
 static uint64_t __attribute__((aligned(8))) bFC=       0xFCFCFCFCFCFCFCFCLL;
 static uint64_t __attribute__((aligned(8))) w10=       0x0010001000100010LL;
@@ -209,26 +199,6 @@ static const uint64_t w1111       __attribute__((aligned(8))) = 0x00010001000100
 // clipping helper table for C implementations:
 static unsigned char clip_table[768];
 
-//global sws_flags from the command line
-int sws_flags=2;
-
-//global srcFilter
-SwsFilter src_filter= {NULL, NULL, NULL, NULL};
-
-float sws_lum_gblur= 0.0;
-float sws_chr_gblur= 0.0;
-int sws_chr_vshift= 0;
-int sws_chr_hshift= 0;
-float sws_chr_sharpen= 0.0;
-float sws_lum_sharpen= 0.0;
-
-/* cpuCaps combined from cpudetect and whats actually compiled in
-   (if there is no support for something compiled in it wont appear here) */
-static CpuCaps cpuCaps;
-
-int (*swScale)(SwsContext *context, uint8_t* src[], int srcStride[], int srcSliceY,
-             int srcSliceH, uint8_t* dst[], int dstStride[])=NULL;
-
 static SwsVector *sws_getConvVec(SwsVector *a, SwsVector *b);
                  
 extern const uint8_t dither_2x2_4[2][8];
@@ -237,7 +207,7 @@ extern const uint8_t dither_8x8_32[8][8];
 extern const uint8_t dither_8x8_73[8][8];
 extern const uint8_t dither_8x8_220[8][8];
 
-#ifdef CAN_COMPILE_X86_ASM
+#ifdef ARCH_X86
 void in_asm_used_var_warning_killer()
 {
  volatile int i= bF8+bFC+w10+
@@ -247,139 +217,6 @@ void in_asm_used_var_warning_killer()
 }
 #endif
 
-static int testFormat[]={
-IMGFMT_YVU9,
-IMGFMT_YV12,
-//IMGFMT_IYUV,
-IMGFMT_I420,
-IMGFMT_BGR15,
-IMGFMT_BGR16,
-IMGFMT_BGR24,
-IMGFMT_BGR32,
-IMGFMT_RGB24,
-IMGFMT_RGB32,
-//IMGFMT_Y8,
-IMGFMT_Y800,
-//IMGFMT_YUY2,
-0
-};
-
-static uint64_t getSSD(uint8_t *src1, uint8_t *src2, int stride1, int stride2, int w, int h){
-       int x,y;
-       uint64_t ssd=0;
-
-       for(y=0; y<h; y++){
-               for(x=0; x<w; x++){
-                       int d= src1[x + y*stride1] - src2[x + y*stride2];
-                       ssd+= d*d;
-               }
-       }
-       return ssd;
-}
-
-// test by ref -> src -> dst -> out & compare out against ref
-// ref & out are YV12
-static void doTest(uint8_t *ref[3], int refStride[3], int w, int h, int srcFormat, int dstFormat, 
-                   int srcW, int srcH, int dstW, int dstH, int flags){
-       uint8_t *src[3];
-       uint8_t *dst[3];
-       uint8_t *out[3];
-       int srcStride[3], dstStride[3];
-       int i;
-       uint64_t ssdY, ssdU, ssdV;
-       SwsContext *srcContext, *dstContext, *outContext;
-       
-       for(i=0; i<3; i++){
-               // avoid stride % bpp != 0
-               if(srcFormat==IMGFMT_RGB24 || srcFormat==IMGFMT_BGR24)
-                       srcStride[i]= srcW*3;
-               else
-                       srcStride[i]= srcW*4;
-               
-               if(dstFormat==IMGFMT_RGB24 || dstFormat==IMGFMT_BGR24)
-                       dstStride[i]= dstW*3;
-               else
-                       dstStride[i]= dstW*4;
-       
-               src[i]= malloc(srcStride[i]*srcH);
-               dst[i]= malloc(dstStride[i]*dstH);
-               out[i]= malloc(refStride[i]*h);
-       }
-
-       srcContext= sws_getContext(w, h, IMGFMT_YV12, srcW, srcH, srcFormat, flags, NULL, NULL);
-       dstContext= sws_getContext(srcW, srcH, srcFormat, dstW, dstH, dstFormat, flags, NULL, NULL);
-       outContext= sws_getContext(dstW, dstH, dstFormat, w, h, IMGFMT_YV12, flags, NULL, NULL);
-       if(srcContext==NULL ||dstContext==NULL ||outContext==NULL){
-               printf("Failed allocating swsContext\n");
-               goto end;
-       }
-//     printf("test %X %X %X -> %X %X %X\n", (int)ref[0], (int)ref[1], (int)ref[2],
-//             (int)src[0], (int)src[1], (int)src[2]);
-
-       srcContext->swScale(srcContext, ref, refStride, 0, h   , src, srcStride);
-       dstContext->swScale(dstContext, src, srcStride, 0, srcH, dst, dstStride);
-       outContext->swScale(outContext, dst, dstStride, 0, dstH, out, refStride);
-            
-       ssdY= getSSD(ref[0], out[0], refStride[0], refStride[0], w, h);
-       ssdU= getSSD(ref[1], out[1], refStride[1], refStride[1], (w+1)>>1, (h+1)>>1);
-       ssdV= getSSD(ref[2], out[2], refStride[2], refStride[2], (w+1)>>1, (h+1)>>1);
-       
-       if(isGray(srcFormat) || isGray(dstFormat)) ssdU=ssdV=0; //FIXME check that output is really gray
-       
-       ssdY/= w*h;
-       ssdU/= w*h/4;
-       ssdV/= w*h/4;
-       
-       if(ssdY>100 || ssdU>50 || ssdV>50){
-               printf(" %s %dx%d -> %s %4dx%4d flags=%2d SSD=%5lld,%5lld,%5lld\n", 
-                       vo_format_name(srcFormat), srcW, srcH, 
-                       vo_format_name(dstFormat), dstW, dstH,
-                       flags,
-                       ssdY, ssdU, ssdV);
-       }
-
-       end:
-       
-       sws_freeContext(srcContext);
-       sws_freeContext(dstContext);
-       sws_freeContext(outContext);
-
-       for(i=0; i<3; i++){
-               free(src[i]);
-               free(dst[i]);
-               free(out[i]);
-       }
-}
-
-static void selfTest(uint8_t *src[3], int stride[3], int w, int h){
-       int srcFormat, dstFormat, srcFormatIndex, dstFormatIndex;
-       int srcW, srcH, dstW, dstH;
-       int flags;
-
-       for(srcFormatIndex=0; ;srcFormatIndex++){
-               srcFormat= testFormat[srcFormatIndex];
-               if(!srcFormat) break;
-               for(dstFormatIndex=0; ;dstFormatIndex++){
-                       dstFormat= testFormat[dstFormatIndex];
-                       if(!dstFormat) break;
-                       if(!isSupportedOut(dstFormat)) continue;
-printf("%s -> %s\n", 
-       vo_format_name(srcFormat),
-       vo_format_name(dstFormat));
-
-                       srcW= w+w/3;
-                       srcH= h+h/3;
-                       for(dstW=w; dstW<w*2; dstW+= dstW/3){
-                               for(dstH=h; dstH<h*2; dstH+= dstH/3){
-                                       for(flags=1; flags<33; flags*=2)
-                                               doTest(src, stride, w, h, srcFormat, dstFormat,
-                                                       srcW, srcH, dstW, dstH, flags);
-                               }
-                       }
-               }
-       }
-}
-
 static inline void yuv2yuvXinC(int16_t *lumFilter, int16_t **lumSrc, int lumFilterSize,
                                    int16_t *chrFilter, int16_t **chrSrc, int chrFilterSize,
                                    uint8_t *dest, uint8_t *uDest, uint8_t *vDest, int dstW, int chrDstW)
@@ -388,7 +225,7 @@ static inline void yuv2yuvXinC(int16_t *lumFilter, int16_t **lumSrc, int lumFilt
        int i;
        for(i=0; i<dstW; i++)
        {
-               int val=0;
+               int val=1<<18;
                int j;
                for(j=0; j<lumFilterSize; j++)
                        val += lumSrc[j][i] * lumFilter[j];
@@ -399,8 +236,8 @@ static inline void yuv2yuvXinC(int16_t *lumFilter, int16_t **lumSrc, int lumFilt
        if(uDest != NULL)
                for(i=0; i<chrDstW; i++)
                {
-                       int u=0;
-                       int v=0;
+                       int u=1<<18;
+                       int v=1<<18;
                        int j;
                        for(j=0; j<chrFilterSize; j++)
                        {
@@ -417,10 +254,10 @@ static inline void yuv2yuvXinC(int16_t *lumFilter, int16_t **lumSrc, int lumFilt
 #define YSCALE_YUV_2_PACKEDX_C(type) \
                for(i=0; i<(dstW>>1); i++){\
                        int j;\
-                       int Y1=0;\
-                       int Y2=0;\
-                       int U=0;\
-                       int V=0;\
+                       int Y1=1<<18;\
+                       int Y2=1<<18;\
+                       int U=1<<18;\
+                       int V=1<<18;\
                        type *r, *b, *g;\
                        const int i2= 2*i;\
                        \
@@ -669,6 +506,14 @@ static inline void yuv2yuvXinC(int16_t *lumFilter, int16_t **lumSrc, int lumFilt
                        ((uint8_t*)dest)[2*i2+3]= V;\
                }               \
                break;\
+       case IMGFMT_UYVY:\
+               func2\
+                       ((uint8_t*)dest)[2*i2+0]= U;\
+                       ((uint8_t*)dest)[2*i2+1]= Y1;\
+                       ((uint8_t*)dest)[2*i2+2]= V;\
+                       ((uint8_t*)dest)[2*i2+3]= Y2;\
+               }               \
+               break;\
        }\
 
 
@@ -779,8 +624,8 @@ static inline void yuv2packedXinC(SwsContext *c, int16_t *lumFilter, int16_t **l
                        int acc=0;
                        for(i=0; i<dstW-1; i+=2){
                                int j;
-                               int Y1=0;
-                               int Y2=0;
+                               int Y1=1<<18;
+                               int Y2=1<<18;
 
                                for(j=0; j<lumFilterSize; j++)
                                {
@@ -813,6 +658,14 @@ static inline void yuv2packedXinC(SwsContext *c, int16_t *lumFilter, int16_t **l
                        ((uint8_t*)dest)[2*i2+3]= V;
                }
                 break;
+       case IMGFMT_UYVY:
+               YSCALE_YUV_2_PACKEDX_C(void)
+                       ((uint8_t*)dest)[2*i2+0]= U;
+                       ((uint8_t*)dest)[2*i2+1]= Y1;
+                       ((uint8_t*)dest)[2*i2+2]= V;
+                       ((uint8_t*)dest)[2*i2+3]= Y2;
+               }
+                break;
        }
 }
 
@@ -823,7 +676,13 @@ static inline void yuv2packedXinC(SwsContext *c, int16_t *lumFilter, int16_t **l
 #define COMPILE_C
 #endif
 
-#ifdef CAN_COMPILE_X86_ASM
+#ifdef ARCH_POWERPC
+#ifdef HAVE_ALTIVEC
+#define COMPILE_ALTIVEC
+#endif //HAVE_ALTIVEC
+#endif //ARCH_POWERPC
+
+#ifdef ARCH_X86
 
 #if (defined (HAVE_MMX) && !defined (HAVE_3DNOW) && !defined (HAVE_MMX2)) || defined (RUNTIME_CPUDETECT)
 #define COMPILE_MMX
@@ -836,7 +695,7 @@ static inline void yuv2packedXinC(SwsContext *c, int16_t *lumFilter, int16_t **l
 #if (defined (HAVE_3DNOW) && !defined (HAVE_MMX2)) || defined (RUNTIME_CPUDETECT)
 #define COMPILE_3DNOW
 #endif
-#endif //CAN_COMPILE_X86_ASM
+#endif //ARCH_X86
 
 #undef HAVE_MMX
 #undef HAVE_MMX2
@@ -846,11 +705,21 @@ static inline void yuv2packedXinC(SwsContext *c, int16_t *lumFilter, int16_t **l
 #undef HAVE_MMX
 #undef HAVE_MMX2
 #undef HAVE_3DNOW
+#undef HAVE_ALTIVEC
 #define RENAME(a) a ## _C
 #include "swscale_template.c"
 #endif
 
-#ifdef CAN_COMPILE_X86_ASM
+#ifdef ARCH_POWERPC
+#ifdef COMPILE_ALTIVEC
+#undef RENAME
+#define HAVE_ALTIVEC
+#define RENAME(a) a ## _altivec
+#include "swscale_template.c"
+#endif
+#endif //ARCH_POWERPC
+
+#ifdef ARCH_X86
 
 //X86 versions
 /*
@@ -892,139 +761,9 @@ static inline void yuv2packedXinC(SwsContext *c, int16_t *lumFilter, int16_t **l
 #include "swscale_template.c"
 #endif
 
-#endif //CAN_COMPILE_X86_ASM
-
-// minor note: the HAVE_xyz is messed up after that line so dont use it
-
-
-// old global scaler, dont use for new code
-// will use sws_flags from the command line
-void SwScale_YV12slice(unsigned char* src[], int srcStride[], int srcSliceY ,
-                            int srcSliceH, uint8_t* dst[], int dstStride, int dstbpp,
-                            int srcW, int srcH, int dstW, int dstH){
-
-       static SwsContext *context=NULL;
-       int dstFormat;
-       int dstStride3[3]= {dstStride, dstStride>>1, dstStride>>1};
-
-       switch(dstbpp)
-       {
-               case 8 : dstFormat= IMGFMT_Y8;          break;
-               case 12: dstFormat= IMGFMT_YV12;        break;
-               case 15: dstFormat= IMGFMT_BGR15;       break;
-               case 16: dstFormat= IMGFMT_BGR16;       break;
-               case 24: dstFormat= IMGFMT_BGR24;       break;
-               case 32: dstFormat= IMGFMT_BGR32;       break;
-               default: return;
-       }
-
-       if(!context) context=sws_getContextFromCmdLine(srcW, srcH, IMGFMT_YV12, dstW, dstH, dstFormat);
-
-       context->swScale(context, src, srcStride, srcSliceY, srcSliceH, dst, dstStride3);
-}
-
-void sws_getFlagsAndFilterFromCmdLine(int *flags, SwsFilter **srcFilterParam, SwsFilter **dstFilterParam)
-{
-       static int firstTime=1;
-       *flags=0;
-
-#ifdef ARCH_X86
-       if(gCpuCaps.hasMMX)
-               asm volatile("emms\n\t"::: "memory"); //FIXME this shouldnt be required but it IS (even for non mmx versions)
-#endif
-       if(firstTime)
-       {
-               firstTime=0;
-               *flags= SWS_PRINT_INFO;
-       }
-       else if(verbose>1) *flags= SWS_PRINT_INFO;
-
-       if(src_filter.lumH) sws_freeVec(src_filter.lumH);
-       if(src_filter.lumV) sws_freeVec(src_filter.lumV);
-       if(src_filter.chrH) sws_freeVec(src_filter.chrH);
-       if(src_filter.chrV) sws_freeVec(src_filter.chrV);
+#endif //ARCH_X86
 
-       if(sws_lum_gblur!=0.0){
-               src_filter.lumH= sws_getGaussianVec(sws_lum_gblur, 3.0);
-               src_filter.lumV= sws_getGaussianVec(sws_lum_gblur, 3.0);
-       }else{
-               src_filter.lumH= sws_getIdentityVec();
-               src_filter.lumV= sws_getIdentityVec();
-       }
-
-       if(sws_chr_gblur!=0.0){
-               src_filter.chrH= sws_getGaussianVec(sws_chr_gblur, 3.0);
-               src_filter.chrV= sws_getGaussianVec(sws_chr_gblur, 3.0);
-       }else{
-               src_filter.chrH= sws_getIdentityVec();
-               src_filter.chrV= sws_getIdentityVec();
-       }
-
-       if(sws_chr_sharpen!=0.0){
-               SwsVector *g= sws_getConstVec(-1.0, 3);
-               SwsVector *id= sws_getConstVec(10.0/sws_chr_sharpen, 1);
-               g->coeff[1]=2.0;
-               sws_addVec(id, g);
-               sws_convVec(src_filter.chrH, id);
-               sws_convVec(src_filter.chrV, id);
-               sws_freeVec(g);
-               sws_freeVec(id);
-       }
-
-       if(sws_lum_sharpen!=0.0){
-               SwsVector *g= sws_getConstVec(-1.0, 3);
-               SwsVector *id= sws_getConstVec(10.0/sws_lum_sharpen, 1);
-               g->coeff[1]=2.0;
-               sws_addVec(id, g);
-               sws_convVec(src_filter.lumH, id);
-               sws_convVec(src_filter.lumV, id);
-               sws_freeVec(g);
-               sws_freeVec(id);
-       }
-
-       if(sws_chr_hshift)
-               sws_shiftVec(src_filter.chrH, sws_chr_hshift);
-
-       if(sws_chr_vshift)
-               sws_shiftVec(src_filter.chrV, sws_chr_vshift);
-
-       sws_normalizeVec(src_filter.chrH, 1.0);
-       sws_normalizeVec(src_filter.chrV, 1.0);
-       sws_normalizeVec(src_filter.lumH, 1.0);
-       sws_normalizeVec(src_filter.lumV, 1.0);
-
-       if(verbose > 1) sws_printVec(src_filter.chrH);
-       if(verbose > 1) sws_printVec(src_filter.lumH);
-
-       switch(sws_flags)
-       {
-               case 0: *flags|= SWS_FAST_BILINEAR; break;
-               case 1: *flags|= SWS_BILINEAR; break;
-               case 2: *flags|= SWS_BICUBIC; break;
-               case 3: *flags|= SWS_X; break;
-               case 4: *flags|= SWS_POINT; break;
-               case 5: *flags|= SWS_AREA; break;
-               case 6: *flags|= SWS_BICUBLIN; break;
-               case 7: *flags|= SWS_GAUSS; break;
-               case 8: *flags|= SWS_SINC; break;
-               case 9: *flags|= SWS_LANCZOS; break;
-               case 10:*flags|= SWS_SPLINE; break;
-               default:*flags|= SWS_BILINEAR; break;
-       }
-       
-       *srcFilterParam= &src_filter;
-       *dstFilterParam= NULL;
-}
-
-// will use sws_flags & src_filter (from cmd line)
-SwsContext *sws_getContextFromCmdLine(int srcW, int srcH, int srcFormat, int dstW, int dstH, int dstFormat)
-{
-       int flags;
-       SwsFilter *dstFilterParam, *srcFilterParam;
-       sws_getFlagsAndFilterFromCmdLine(&flags, &srcFilterParam, &dstFilterParam);
-
-       return sws_getContext(srcW, srcH, srcFormat, dstW, dstH, dstFormat, flags, srcFilterParam, dstFilterParam);
-}
+// minor note: the HAVE_xyz is messed up after that line so don't use it
 
 static double getSplineCoeff(double a, double b, double c, double d, double dist)
 {
@@ -1048,7 +787,7 @@ static inline void initFilter(int16_t **outFilter, int16_t **filterPos, int *out
        double *filter=NULL;
        double *filter2=NULL;
 #ifdef ARCH_X86
-       if(gCpuCaps.hasMMX)
+       if(flags & SWS_CPU_CAPS_MMX)
                asm volatile("emms\n\t"::: "memory"); //FIXME this shouldnt be required but it IS (even for non mmx versions)
 #endif
 
@@ -1279,7 +1018,7 @@ static inline void initFilter(int16_t **outFilter, int16_t **filterPos, int *out
 
                        if(cutOff > SWS_MAX_REDUCE_CUTOFF) break;
 
-                       /* preserve Monotonicity because the core cant handle the filter otherwise */
+                       /* preserve Monotonicity because the core can't handle the filter otherwise */
                        if(i<dstW-1 && (*filterPos)[i] >= (*filterPos)[i+1]) break;
 
                        // Move filter coeffs left
@@ -1302,6 +1041,21 @@ static inline void initFilter(int16_t **outFilter, int16_t **filterPos, int *out
                if(min>minFilterSize) minFilterSize= min;
        }
 
+        if (flags & SWS_CPU_CAPS_ALTIVEC) {
+          // we can handle the special case 4,
+          // so we don't want to go to the full 8
+          if (minFilterSize < 5)
+            filterAlign = 4;
+
+          // we really don't want to waste our time
+          // doing useless computation, so fall-back on
+          // the scalar C code for very small filter.
+          // vectorizing is worth it only if you have
+          // decent-sized vector.
+          if (minFilterSize < 3)
+            filterAlign = 1;
+        }
+
        ASSERT(minFilterSize > 0)
        filterSize= (minFilterSize +(filterAlign-1)) & (~(filterAlign-1));
        ASSERT(filterSize > 0)
@@ -1364,8 +1118,10 @@ static inline void initFilter(int16_t **outFilter, int16_t **filterPos, int *out
        for(i=0; i<dstW; i++)
        {
                int j;
+               double error=0;
                double sum=0;
                double scale= one;
+
                for(j=0; j<filterSize; j++)
                {
                        sum+= filter[i*filterSize + j];
@@ -1373,7 +1129,10 @@ static inline void initFilter(int16_t **outFilter, int16_t **filterPos, int *out
                scale/= sum;
                for(j=0; j<*outFilterSize; j++)
                {
-                       (*outFilter)[i*(*outFilterSize) + j]= (int)(filter[i*filterSize + j]*scale);
+                       double v= filter[i*filterSize + j]*scale + error;
+                       int intV= floor(v + 0.5);
+                       (*outFilter)[i*(*outFilterSize) + j]= intV;
+                       error = v - intV;
                }
        }
        
@@ -1560,10 +1319,6 @@ static void initMMX2HScaler(int dstW, int xInc, uint8_t *funnyCode, int16_t *fil
 }
 #endif // ARCH_X86
 
-//FIXME remove
-void SwScale_Init(){
-}
-
 static void globalInit(){
     // generating tables:
     int i;
@@ -1571,38 +1326,42 @@ static void globalInit(){
        int c= MIN(MAX(i-256, 0), 255);
        clip_table[i]=c;
     }
+}
 
-cpuCaps= gCpuCaps;
-
+static SwsFunc getSwsFunc(int flags){
+    
 #ifdef RUNTIME_CPUDETECT
-#ifdef CAN_COMPILE_X86_ASM
+#ifdef ARCH_X86
        // ordered per speed fasterst first
-       if(gCpuCaps.hasMMX2)
-               swScale= swScale_MMX2;
-       else if(gCpuCaps.has3DNow)
-               swScale= swScale_3DNow;
-       else if(gCpuCaps.hasMMX)
-               swScale= swScale_MMX;
+       if(flags & SWS_CPU_CAPS_MMX2)
+               return swScale_MMX2;
+       else if(flags & SWS_CPU_CAPS_3DNOW)
+               return swScale_3DNow;
+       else if(flags & SWS_CPU_CAPS_MMX)
+               return swScale_MMX;
        else
-               swScale= swScale_C;
+               return swScale_C;
 
 #else
-       swScale= swScale_C;
-       cpuCaps.hasMMX2 = cpuCaps.hasMMX = cpuCaps.has3DNow = 0;
+#ifdef ARCH_POWERPC
+       if(flags & SWS_CPU_CAPS_ALTIVEC)
+         return swScale_altivec;
+       else
+         return swScale_C;
+#endif
+       return swScale_C;
 #endif
 #else //RUNTIME_CPUDETECT
 #ifdef HAVE_MMX2
-       swScale= swScale_MMX2;
-       cpuCaps.has3DNow = 0;
+       return swScale_MMX2;
 #elif defined (HAVE_3DNOW)
-       swScale= swScale_3DNow;
-       cpuCaps.hasMMX2 = 0;
+       return swScale_3DNow;
 #elif defined (HAVE_MMX)
-       swScale= swScale_MMX;
-       cpuCaps.hasMMX2 = cpuCaps.has3DNow = 0;
+       return swScale_MMX;
+#elif defined (HAVE_ALTIVEC)
+       return swScale_altivec;
 #else
-       swScale= swScale_C;
-       cpuCaps.hasMMX2 = cpuCaps.hasMMX = cpuCaps.has3DNow = 0;
+       return swScale_C;
 #endif
 #endif //!RUNTIME_CPUDETECT
 }
@@ -1640,6 +1399,15 @@ static int PlanarToYuy2Wrapper(SwsContext *c, uint8_t* src[], int srcStride[], i
        return srcSliceH;
 }
 
+static int PlanarToUyvyWrapper(SwsContext *c, uint8_t* src[], int srcStride[], int srcSliceY,
+             int srcSliceH, uint8_t* dstParam[], int dstStride[]){
+       uint8_t *dst=dstParam[0] + dstStride[0]*srcSliceY;
+
+       yv12touyvy( src[0],src[1],src[2],dst,c->srcW,srcSliceH,srcStride[0],srcStride[1],dstStride[0] );
+
+       return srcSliceH;
+}
+
 /* {RGB,BGR}{15,16,24,32} -> {RGB,BGR}{15,16,24,32} */
 static int rgb2rgbWrapper(SwsContext *c, uint8_t* src[], int srcStride[], int srcSliceY,
                           int srcSliceH, uint8_t* dst[], int dstStride[]){
@@ -1770,7 +1538,7 @@ static inline void sws_orderYUV(int format, uint8_t * sortedP[], int sortedStrid
                sortedStride[1]= stride[2];
                sortedStride[2]= stride[1];
        }
-       else if(isPacked(format) || isGray(format))
+       else if(isPacked(format) || isGray(format) || format == IMGFMT_Y8)
        {
                sortedP[0]= p[0];
                sortedP[1]= 
@@ -1984,17 +1752,31 @@ SwsContext *sws_getContext(int srcW, int srcH, int origSrcFormat, int dstW, int
 
        SwsContext *c;
        int i;
-       int usesFilter;
+       int usesVFilter, usesHFilter;
        int unscaled, needsDither;
        int srcFormat, dstFormat;
        SwsFilter dummyFilter= {NULL, NULL, NULL, NULL};
 #ifdef ARCH_X86
-       if(gCpuCaps.hasMMX)
+       if(flags & SWS_CPU_CAPS_MMX)
                asm volatile("emms\n\t"::: "memory");
 #endif
-       if(swScale==NULL) globalInit();
 
-       /* avoid dupplicate Formats, so we dont need to check to much */
+#ifndef RUNTIME_CPUDETECT //ensure that the flags match the compiled variant if cpudetect is off
+       flags &= ~(SWS_CPU_CAPS_MMX|SWS_CPU_CAPS_MMX2|SWS_CPU_CAPS_3DNOW|SWS_CPU_CAPS_ALTIVEC);
+#ifdef HAVE_MMX2
+       flags |= SWS_CPU_CAPS_MMX|SWS_CPU_CAPS_MMX2;
+#elif defined (HAVE_3DNOW)
+       flags |= SWS_CPU_CAPS_MMX|SWS_CPU_CAPS_3DNOW;
+#elif defined (HAVE_MMX)
+       flags |= SWS_CPU_CAPS_MMX;
+#elif defined (HAVE_ALTIVEC)
+       flags |= SWS_CPU_CAPS_ALTIVEC;
+#endif
+#endif
+       if(clip_table[512] != 255) globalInit();
+       if(rgb15to16 == NULL) sws_rgb2rgb_init(flags);
+
+       /* avoid duplicate Formats, so we don't need to check to much */
        srcFormat = remove_dup_fourcc(origSrcFormat);
        dstFormat = remove_dup_fourcc(origDstFormat);
 
@@ -2039,16 +1821,17 @@ SwsContext *sws_getContext(int srcW, int srcH, int origSrcFormat, int dstW, int
        c->srcFormat= srcFormat;
        c->origDstFormat= origDstFormat;
        c->origSrcFormat= origSrcFormat;
-
-       usesFilter=0;
-       if(dstFilter->lumV!=NULL && dstFilter->lumV->length>1) usesFilter=1;
-       if(dstFilter->lumH!=NULL && dstFilter->lumH->length>1) usesFilter=1;
-       if(dstFilter->chrV!=NULL && dstFilter->chrV->length>1) usesFilter=1;
-       if(dstFilter->chrH!=NULL && dstFilter->chrH->length>1) usesFilter=1;
-       if(srcFilter->lumV!=NULL && srcFilter->lumV->length>1) usesFilter=1;
-       if(srcFilter->lumH!=NULL && srcFilter->lumH->length>1) usesFilter=1;
-       if(srcFilter->chrV!=NULL && srcFilter->chrV->length>1) usesFilter=1;
-       if(srcFilter->chrH!=NULL && srcFilter->chrH->length>1) usesFilter=1;
+        c->vRounder= 4* 0x0001000100010001ULL;
+
+       usesHFilter= usesVFilter= 0;
+       if(dstFilter->lumV!=NULL && dstFilter->lumV->length>1) usesVFilter=1;
+       if(dstFilter->lumH!=NULL && dstFilter->lumH->length>1) usesHFilter=1;
+       if(dstFilter->chrV!=NULL && dstFilter->chrV->length>1) usesVFilter=1;
+       if(dstFilter->chrH!=NULL && dstFilter->chrH->length>1) usesHFilter=1;
+       if(srcFilter->lumV!=NULL && srcFilter->lumV->length>1) usesVFilter=1;
+       if(srcFilter->lumH!=NULL && srcFilter->lumH->length>1) usesHFilter=1;
+       if(srcFilter->chrV!=NULL && srcFilter->chrV->length>1) usesVFilter=1;
+       if(srcFilter->chrH!=NULL && srcFilter->chrH->length>1) usesHFilter=1;
 
        getSubSampleFactors(&c->chrSrcHSubSample, &c->chrSrcVSubSample, srcFormat);
        getSubSampleFactors(&c->chrDstHSubSample, &c->chrDstVSubSample, dstFormat);
@@ -2076,7 +1859,7 @@ SwsContext *sws_getContext(int srcW, int srcH, int origSrcFormat, int dstW, int
        sws_setColorspaceDetails(c, Inverse_Table_6_9[SWS_CS_DEFAULT], 0, Inverse_Table_6_9[SWS_CS_DEFAULT] /* FIXME*/, 0, 0, 1<<16, 1<<16); 
 
        /* unscaled special Cases */
-       if(unscaled && !usesFilter)
+       if(unscaled && !usesHFilter && !usesVFilter)
        {
                /* yv12_to_nv12 */
                if(srcFormat == IMGFMT_YV12 && dstFormat == IMGFMT_NV12)
@@ -2113,9 +1896,13 @@ SwsContext *sws_getContext(int srcW, int srcH, int origSrcFormat, int dstW, int
                                c->swScale= rgb2rgbWrapper;
 
                        /* yv12_to_yuy2 */
-                       if(srcFormat == IMGFMT_YV12 && dstFormat == IMGFMT_YUY2)
+                       if(srcFormat == IMGFMT_YV12 && 
+                           (dstFormat == IMGFMT_YUY2 || dstFormat == IMGFMT_UYVY))
                        {
-                               c->swScale= PlanarToYuy2Wrapper;
+                               if (dstFormat == IMGFMT_YUY2)
+                                   c->swScale= PlanarToYuy2Wrapper;
+                               else
+                                   c->swScale= PlanarToUyvyWrapper;
                        }
                }
 
@@ -2136,7 +1923,7 @@ SwsContext *sws_getContext(int srcW, int srcH, int origSrcFormat, int dstW, int
                }
        }
 
-       if(cpuCaps.hasMMX2)
+       if(flags & SWS_CPU_CAPS_MMX2)
        {
                c->canMMX2BeUsed= (dstW >=srcW && (dstW&31)==0 && (srcW&15)==0) ? 1 : 0;
                if(!c->canMMX2BeUsed && dstW >=srcW && (srcW&15)==0 && (flags&SWS_FAST_BILINEAR))
@@ -2144,6 +1931,7 @@ SwsContext *sws_getContext(int srcW, int srcH, int origSrcFormat, int dstW, int
                        if(flags&SWS_PRINT_INFO)
                                MSG_INFO("SwScaler: output Width is not a multiple of 32 -> no MMX2 scaler\n");
                }
+               if(usesHFilter) c->canMMX2BeUsed=0;
        }
        else
                c->canMMX2BeUsed=0;
@@ -2164,8 +1952,8 @@ SwsContext *sws_getContext(int srcW, int srcH, int origSrcFormat, int dstW, int
                        c->lumXInc+= 20;
                        c->chrXInc+= 20;
                }
-               //we dont use the x86asm scaler if mmx is available
-               else if(cpuCaps.hasMMX)
+               //we don't use the x86asm scaler if mmx is available
+               else if(flags & SWS_CPU_CAPS_MMX)
                {
                        c->lumXInc = ((srcW-2)<<16)/(dstW-2) - 20;
                        c->chrXInc = ((c->chrSrcW-2)<<16)/(c->chrDstW-2) - 20;
@@ -2174,7 +1962,10 @@ SwsContext *sws_getContext(int srcW, int srcH, int origSrcFormat, int dstW, int
 
        /* precalculate horizontal scaler filter coefficients */
        {
-               const int filterAlign= cpuCaps.hasMMX ? 4 : 1;
+               const int filterAlign=
+                 (flags & SWS_CPU_CAPS_MMX) ? 4 :
+                 (flags & SWS_CPU_CAPS_ALTIVEC) ? 8 :
+                 1;
 
                initFilter(&c->hLumFilter, &c->hLumFilterPos, &c->hLumFilterSize, c->lumXInc,
                                 srcW      ,       dstW, filterAlign, 1<<14,
@@ -2186,7 +1977,7 @@ SwsContext *sws_getContext(int srcW, int srcH, int origSrcFormat, int dstW, int
                                 srcFilter->chrH, dstFilter->chrH);
 
 #ifdef ARCH_X86
-// cant downscale !!!
+// can't downscale !!!
                if(c->canMMX2BeUsed && (flags & SWS_FAST_BILINEAR))
                {
                        c->lumMmx2Filter   = (int16_t*)memalign(8, (dstW        /8+8)*sizeof(int16_t));
@@ -2203,16 +1994,22 @@ SwsContext *sws_getContext(int srcW, int srcH, int origSrcFormat, int dstW, int
 
 
        /* precalculate vertical scaler filter coefficients */
-       initFilter(&c->vLumFilter, &c->vLumFilterPos, &c->vLumFilterSize, c->lumYInc,
-                       srcH      ,        dstH, 1, (1<<12)-4,
-                       (flags&SWS_BICUBLIN) ? (flags|SWS_BICUBIC)  : flags,
-                       srcFilter->lumV, dstFilter->lumV);
-       initFilter(&c->vChrFilter, &c->vChrFilterPos, &c->vChrFilterSize, c->chrYInc,
-                       c->chrSrcH, c->chrDstH, 1, (1<<12)-4,
-                       (flags&SWS_BICUBLIN) ? (flags|SWS_BILINEAR) : flags,
-                       srcFilter->chrV, dstFilter->chrV);
-
-       // Calculate Buffer Sizes so that they wont run out while handling these damn slices
+       {
+               const int filterAlign=
+                 (flags & SWS_CPU_CAPS_ALTIVEC) ? 8 :
+                 1;
+
+               initFilter(&c->vLumFilter, &c->vLumFilterPos, &c->vLumFilterSize, c->lumYInc,
+                               srcH      ,        dstH, filterAlign, (1<<12)-4,
+                               (flags&SWS_BICUBLIN) ? (flags|SWS_BICUBIC)  : flags,
+                               srcFilter->lumV, dstFilter->lumV);
+               initFilter(&c->vChrFilter, &c->vChrFilterPos, &c->vChrFilterSize, c->chrYInc,
+                               c->chrSrcH, c->chrDstH, filterAlign, (1<<12)-4,
+                               (flags&SWS_BICUBLIN) ? (flags|SWS_BILINEAR) : flags,
+                               srcFilter->chrV, dstFilter->chrV);
+       }
+
+       // Calculate Buffer Sizes so that they won't run out while handling these damn slices
        c->vLumBufSize= c->vLumFilterSize;
        c->vChrBufSize= c->vChrFilterSize;
        for(i=0; i<dstH; i++)
@@ -2220,7 +2017,9 @@ SwsContext *sws_getContext(int srcW, int srcH, int origSrcFormat, int dstW, int
                int chrI= i*c->chrDstH / dstH;
                int nextSlice= MAX(c->vLumFilterPos[i   ] + c->vLumFilterSize - 1,
                                 ((c->vChrFilterPos[chrI] + c->vChrFilterSize - 1)<<c->chrSrcVSubSample));
-               nextSlice&= ~3; // Slices start at boundaries which are divisable through 4
+
+               nextSlice>>= c->chrSrcVSubSample;
+               nextSlice<<= c->chrSrcVSubSample;
                if(c->vLumFilterPos[i   ] + c->vLumBufSize < nextSlice)
                        c->vLumBufSize= nextSlice - c->vLumFilterPos[i   ];
                if(c->vChrFilterPos[chrI] + c->vChrBufSize < (nextSlice>>c->chrSrcVSubSample))
@@ -2281,19 +2080,21 @@ SwsContext *sws_getContext(int srcW, int srcH, int origSrcFormat, int dstW, int
                        MSG_INFO("from %s to %s ", 
                                vo_format_name(srcFormat), vo_format_name(dstFormat));
 
-               if(cpuCaps.hasMMX2)
+               if(flags & SWS_CPU_CAPS_MMX2)
                        MSG_INFO("using MMX2\n");
-               else if(cpuCaps.has3DNow)
+               else if(flags & SWS_CPU_CAPS_3DNOW)
                        MSG_INFO("using 3DNOW\n");
-               else if(cpuCaps.hasMMX)
+               else if(flags & SWS_CPU_CAPS_MMX)
                        MSG_INFO("using MMX\n");
-               else
+               else if(flags & SWS_CPU_CAPS_ALTIVEC)
+                       MSG_INFO("using AltiVec\n");
+               else 
                        MSG_INFO("using C\n");
        }
 
-       if((flags & SWS_PRINT_INFO) && verbose>0)
+       if(flags & SWS_PRINT_INFO)
        {
-               if(cpuCaps.hasMMX)
+               if(flags & SWS_CPU_CAPS_MMX)
                {
                        if(c->canMMX2BeUsed && (flags&SWS_FAST_BILINEAR))
                                MSG_V("SwScaler: using FAST_BILINEAR MMX2 scaler for horizontal scaling\n");
@@ -2328,34 +2129,34 @@ SwsContext *sws_getContext(int srcW, int srcH, int origSrcFormat, int dstW, int
                if(isPlanarYUV(dstFormat))
                {
                        if(c->vLumFilterSize==1)
-                               MSG_V("SwScaler: using 1-tap %s \"scaler\" for vertical scaling (YV12 like)\n", cpuCaps.hasMMX ? "MMX" : "C");
+                               MSG_V("SwScaler: using 1-tap %s \"scaler\" for vertical scaling (YV12 like)\n", (flags & SWS_CPU_CAPS_MMX) ? "MMX" : "C");
                        else
-                               MSG_V("SwScaler: using n-tap %s scaler for vertical scaling (YV12 like)\n", cpuCaps.hasMMX ? "MMX" : "C");
+                               MSG_V("SwScaler: using n-tap %s scaler for vertical scaling (YV12 like)\n", (flags & SWS_CPU_CAPS_MMX) ? "MMX" : "C");
                }
                else
                {
                        if(c->vLumFilterSize==1 && c->vChrFilterSize==2)
                                MSG_V("SwScaler: using 1-tap %s \"scaler\" for vertical luminance scaling (BGR)\n"
-                                      "SwScaler:       2-tap scaler for vertical chrominance scaling (BGR)\n",cpuCaps.hasMMX ? "MMX" : "C");
+                                      "SwScaler:       2-tap scaler for vertical chrominance scaling (BGR)\n",(flags & SWS_CPU_CAPS_MMX) ? "MMX" : "C");
                        else if(c->vLumFilterSize==2 && c->vChrFilterSize==2)
-                               MSG_V("SwScaler: using 2-tap linear %s scaler for vertical scaling (BGR)\n", cpuCaps.hasMMX ? "MMX" : "C");
+                               MSG_V("SwScaler: using 2-tap linear %s scaler for vertical scaling (BGR)\n", (flags & SWS_CPU_CAPS_MMX) ? "MMX" : "C");
                        else
-                               MSG_V("SwScaler: using n-tap %s scaler for vertical scaling (BGR)\n", cpuCaps.hasMMX ? "MMX" : "C");
+                               MSG_V("SwScaler: using n-tap %s scaler for vertical scaling (BGR)\n", (flags & SWS_CPU_CAPS_MMX) ? "MMX" : "C");
                }
 
                if(dstFormat==IMGFMT_BGR24)
                        MSG_V("SwScaler: using %s YV12->BGR24 Converter\n",
-                               cpuCaps.hasMMX2 ? "MMX2" : (cpuCaps.hasMMX ? "MMX" : "C"));
+                               (flags & SWS_CPU_CAPS_MMX2) ? "MMX2" : ((flags & SWS_CPU_CAPS_MMX) ? "MMX" : "C"));
                else if(dstFormat==IMGFMT_BGR32)
-                       MSG_V("SwScaler: using %s YV12->BGR32 Converter\n", cpuCaps.hasMMX ? "MMX" : "C");
+                       MSG_V("SwScaler: using %s YV12->BGR32 Converter\n", (flags & SWS_CPU_CAPS_MMX) ? "MMX" : "C");
                else if(dstFormat==IMGFMT_BGR16)
-                       MSG_V("SwScaler: using %s YV12->BGR16 Converter\n", cpuCaps.hasMMX ? "MMX" : "C");
+                       MSG_V("SwScaler: using %s YV12->BGR16 Converter\n", (flags & SWS_CPU_CAPS_MMX) ? "MMX" : "C");
                else if(dstFormat==IMGFMT_BGR15)
-                       MSG_V("SwScaler: using %s YV12->BGR15 Converter\n", cpuCaps.hasMMX ? "MMX" : "C");
+                       MSG_V("SwScaler: using %s YV12->BGR15 Converter\n", (flags & SWS_CPU_CAPS_MMX) ? "MMX" : "C");
 
                MSG_V("SwScaler: %dx%d -> %dx%d\n", srcW, srcH, dstW, dstH);
        }
-       if((flags & SWS_PRINT_INFO) && verbose>1)
+       if(flags & SWS_PRINT_INFO)
        {
                MSG_DBG2("SwScaler:Lum srcW=%d srcH=%d dstW=%d dstH=%d xInc=%d yInc=%d\n",
                        c->srcW, c->srcH, c->dstW, c->dstH, c->lumXInc, c->lumYInc);
@@ -2363,21 +2164,24 @@ SwsContext *sws_getContext(int srcW, int srcH, int origSrcFormat, int dstW, int
                        c->chrSrcW, c->chrSrcH, c->chrDstW, c->chrDstH, c->chrXInc, c->chrYInc);
        }
 
-       c->swScale= swScale;
+       c->swScale= getSwsFunc(flags);
        return c;
 }
 
 /**
- * swscale warper, so we dont need to export the SwsContext.
+ * swscale warper, so we don't need to export the SwsContext.
  * assumes planar YUV to be in YUV order instead of YVU
  */
 int sws_scale_ordered(SwsContext *c, uint8_t* src[], int srcStride[], int srcSliceY,
                            int srcSliceH, uint8_t* dst[], int dstStride[]){
-       c->swScale(c, src, srcStride, srcSliceY, srcSliceH, dst, dstStride);
+       //copy strides, so they can safely be modified
+       int srcStride2[3]= {srcStride[0], srcStride[1], srcStride[2]};
+       int dstStride2[3]= {dstStride[0], dstStride[1], dstStride[2]};
+       return c->swScale(c, src, srcStride2, srcSliceY, srcSliceH, dst, dstStride2);
 }
 
 /**
- * swscale warper, so we dont need to export the SwsContext
+ * swscale warper, so we don't need to export the SwsContext
  */
 int sws_scale(SwsContext *c, uint8_t* srcParam[], int srcStrideParam[], int srcSliceY,
                            int srcSliceH, uint8_t* dstParam[], int dstStrideParam[]){
@@ -2388,14 +2192,76 @@ int sws_scale(SwsContext *c, uint8_t* srcParam[], int srcStrideParam[], int srcS
        sws_orderYUV(c->origSrcFormat, src, srcStride, srcParam, srcStrideParam);
        sws_orderYUV(c->origDstFormat, dst, dstStride, dstParam, dstStrideParam);
 //printf("sws: slice %d %d\n", srcSliceY, srcSliceH);
-       c->swScale(c, src, srcStride, srcSliceY, srcSliceH, dst, dstStride);
+
+       return c->swScale(c, src, srcStride, srcSliceY, srcSliceH, dst, dstStride);
+}
+
+SwsFilter *sws_getDefaultFilter(float lumaGBlur, float chromaGBlur, 
+                               float lumaSharpen, float chromaSharpen,
+                               float chromaHShift, float chromaVShift,
+                               int verbose)
+{
+       SwsFilter *filter= malloc(sizeof(SwsFilter));
+
+       if(lumaGBlur!=0.0){
+               filter->lumH= sws_getGaussianVec(lumaGBlur, 3.0);
+               filter->lumV= sws_getGaussianVec(lumaGBlur, 3.0);
+       }else{
+               filter->lumH= sws_getIdentityVec();
+               filter->lumV= sws_getIdentityVec();
+       }
+
+       if(chromaGBlur!=0.0){
+               filter->chrH= sws_getGaussianVec(chromaGBlur, 3.0);
+               filter->chrV= sws_getGaussianVec(chromaGBlur, 3.0);
+       }else{
+               filter->chrH= sws_getIdentityVec();
+               filter->chrV= sws_getIdentityVec();
+       }
+
+       if(chromaSharpen!=0.0){
+               SwsVector *g= sws_getConstVec(-1.0, 3);
+               SwsVector *id= sws_getConstVec(10.0/chromaSharpen, 1);
+               g->coeff[1]=2.0;
+               sws_addVec(id, g);
+               sws_convVec(filter->chrH, id);
+               sws_convVec(filter->chrV, id);
+               sws_freeVec(g);
+               sws_freeVec(id);
+       }
+
+       if(lumaSharpen!=0.0){
+               SwsVector *g= sws_getConstVec(-1.0, 3);
+               SwsVector *id= sws_getConstVec(10.0/lumaSharpen, 1);
+               g->coeff[1]=2.0;
+               sws_addVec(id, g);
+               sws_convVec(filter->lumH, id);
+               sws_convVec(filter->lumV, id);
+               sws_freeVec(g);
+               sws_freeVec(id);
+       }
+
+       if(chromaHShift != 0.0)
+               sws_shiftVec(filter->chrH, (int)(chromaHShift+0.5));
+
+       if(chromaVShift != 0.0)
+               sws_shiftVec(filter->chrV, (int)(chromaVShift+0.5));
+
+       sws_normalizeVec(filter->chrH, 1.0);
+       sws_normalizeVec(filter->chrV, 1.0);
+       sws_normalizeVec(filter->lumH, 1.0);
+       sws_normalizeVec(filter->lumV, 1.0);
+
+       if(verbose) sws_printVec(filter->chrH);
+       if(verbose) sws_printVec(filter->lumH);
+
+        return filter;
 }
 
 /**
  * returns a normalized gaussian curve used to filter stuff
  * quality=3 is high quality, lowwer is lowwer quality
  */
-
 SwsVector *sws_getGaussianVec(double variance, double quality){
        const int length= (int)(variance*quality + 0.5) | 1;
        int i;
@@ -2616,6 +2482,17 @@ void sws_freeVec(SwsVector *a){
        free(a);
 }
 
+void sws_freeFilter(SwsFilter *filter){
+       if(!filter) return;
+
+       if(filter->lumH) sws_freeVec(filter->lumH);
+       if(filter->lumV) sws_freeVec(filter->lumV);
+       if(filter->chrH) sws_freeVec(filter->chrH);
+       if(filter->chrV) sws_freeVec(filter->chrV);
+       free(filter);
+}
+
+
 void sws_freeContext(SwsContext *c){
        int i;
        if(!c) return;