]> git.sesse.net Git - ffmpeg/blobdiff - postproc/swscale_template.c
1000l
[ffmpeg] / postproc / swscale_template.c
index 4e106d6e98d50964c684c75b746e83b17217cce5..ebdb6c33d11d2dba274fb39bbc988d753b038b0d 100644 (file)
@@ -1,9 +1,20 @@
+/*
+    Copyright (C) 2001-2003 Michael Niedermayer <michaelni@gmx.at>
+
+    This program is free software; you can redistribute it and/or modify
+    it under the terms of the GNU General Public License as published by
+    the Free Software Foundation; either version 2 of the License, or
+    (at your option) any later version.
 
-// Software scaling and colorspace conversion routines for MPlayer
+    This program is distributed in the hope that it will be useful,
+    but WITHOUT ANY WARRANTY; without even the implied warranty of
+    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+    GNU General Public License for more details.
 
-// Orginal C implementation by A'rpi/ESP-team <arpi@thot.banki.hu>
-// current version mostly by Michael Niedermayer (michaelni@gmx.at)
-// the parts written by michael are under GNU GPL
+    You should have received a copy of the GNU General Public License
+    along with this program; if not, write to the Free Software
+    Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
+*/
 
 #undef MOVNTQ
 #undef PAVGB
 #define MOVNTQ(a,b) "movq " #a ", " #b " \n\t"
 #endif
 
-
-#define YSCALEYUV2YV12X(x) \
+#define YSCALEYUV2YV12X(x, offset) \
                        "xorl %%eax, %%eax              \n\t"\
-                       "pxor %%mm3, %%mm3              \n\t"\
-                       "pxor %%mm4, %%mm4              \n\t"\
-                       "movl %0, %%edx                 \n\t"\
+                       "movq "VROUNDER_OFFSET"(%0), %%mm3\n\t"\
+                       "movq %%mm3, %%mm4              \n\t"\
+                       "leal " offset "(%0), %%edx     \n\t"\
+                       "movl (%%edx), %%esi            \n\t"\
                        ".balign 16                     \n\t" /* FIXME Unroll? */\
                        "1:                             \n\t"\
-                       "movl (%1, %%edx, 4), %%esi     \n\t"\
-                       "movq (%2, %%edx, 8), %%mm0     \n\t" /* filterCoeff */\
+                       "movq 8(%%edx), %%mm0           \n\t" /* filterCoeff */\
                        "movq " #x "(%%esi, %%eax, 2), %%mm2    \n\t" /* srcData */\
                        "movq 8+" #x "(%%esi, %%eax, 2), %%mm5  \n\t" /* srcData */\
+                       "addl $16, %%edx                \n\t"\
+                       "movl (%%edx), %%esi            \n\t"\
+                       "testl %%esi, %%esi             \n\t"\
                        "pmulhw %%mm0, %%mm2            \n\t"\
                        "pmulhw %%mm0, %%mm5            \n\t"\
                        "paddw %%mm2, %%mm3             \n\t"\
                        "paddw %%mm5, %%mm4             \n\t"\
-                       "addl $1, %%edx                 \n\t"\
                        " jnz 1b                        \n\t"\
                        "psraw $3, %%mm3                \n\t"\
                        "psraw $3, %%mm4                \n\t"\
                        "packuswb %%mm4, %%mm3          \n\t"\
-                       MOVNTQ(%%mm3, (%3, %%eax))\
+                       MOVNTQ(%%mm3, (%1, %%eax))\
                        "addl $8, %%eax                 \n\t"\
-                       "cmpl %4, %%eax                 \n\t"\
-                       "pxor %%mm3, %%mm3              \n\t"\
-                       "pxor %%mm4, %%mm4              \n\t"\
-                       "movl %0, %%edx                 \n\t"\
+                       "cmpl %2, %%eax                 \n\t"\
+                       "movq "VROUNDER_OFFSET"(%0), %%mm3\n\t"\
+                       "movq %%mm3, %%mm4              \n\t"\
+                       "leal " offset "(%0), %%edx     \n\t"\
+                       "movl (%%edx), %%esi            \n\t"\
                        "jb 1b                          \n\t"
 
 #define YSCALEYUV2YV121 \
                           "m" (lumSrc+lumFilterSize), "m" (chrSrc+chrFilterSize)
                        : "%eax", "%ebx", "%ecx", "%edx", "%esi"
 */
-#define YSCALEYUV2RGBX \
+#define YSCALEYUV2PACKEDX \
                "xorl %%eax, %%eax              \n\t"\
                ".balign 16                     \n\t"\
+               "nop                            \n\t"\
                "1:                             \n\t"\
-               "movl %1, %%edx                 \n\t" /* -chrFilterSize */\
-               "movl %3, %%ebx                 \n\t" /* chrMmxFilter+lumFilterSize */\
-               "movl %7, %%ecx                 \n\t" /* chrSrc+lumFilterSize */\
-               "pxor %%mm3, %%mm3              \n\t"\
-               "pxor %%mm4, %%mm4              \n\t"\
+               "leal "CHR_MMX_FILTER_OFFSET"(%0), %%edx        \n\t"\
+               "movl (%%edx), %%esi            \n\t"\
+               "movq "VROUNDER_OFFSET"(%0), %%mm3\n\t"\
+               "movq %%mm3, %%mm4              \n\t"\
+               ".balign 16                     \n\t"\
                "2:                             \n\t"\
-               "movl (%%ecx, %%edx, 4), %%esi  \n\t"\
-               "movq (%%ebx, %%edx, 8), %%mm0  \n\t" /* filterCoeff */\
+               "movq 8(%%edx), %%mm0           \n\t" /* filterCoeff */\
                "movq (%%esi, %%eax), %%mm2     \n\t" /* UsrcData */\
                "movq 4096(%%esi, %%eax), %%mm5 \n\t" /* VsrcData */\
+               "addl $16, %%edx                \n\t"\
+               "movl (%%edx), %%esi            \n\t"\
                "pmulhw %%mm0, %%mm2            \n\t"\
                "pmulhw %%mm0, %%mm5            \n\t"\
                "paddw %%mm2, %%mm3             \n\t"\
                "paddw %%mm5, %%mm4             \n\t"\
-               "addl $1, %%edx                 \n\t"\
+               "testl %%esi, %%esi             \n\t"\
                " jnz 2b                        \n\t"\
 \
-               "movl %0, %%edx                 \n\t" /* -lumFilterSize */\
-               "movl %2, %%ebx                 \n\t" /* lumMmxFilter+lumFilterSize */\
-               "movl %6, %%ecx                 \n\t" /* lumSrc+lumFilterSize */\
-               "pxor %%mm1, %%mm1              \n\t"\
-               "pxor %%mm7, %%mm7              \n\t"\
+               "leal "LUM_MMX_FILTER_OFFSET"(%0), %%edx        \n\t"\
+               "movl (%%edx), %%esi            \n\t"\
+               "movq "VROUNDER_OFFSET"(%0), %%mm1\n\t"\
+               "movq %%mm1, %%mm7              \n\t"\
+               ".balign 16                     \n\t"\
                "2:                             \n\t"\
-               "movl (%%ecx, %%edx, 4), %%esi  \n\t"\
-               "movq (%%ebx, %%edx, 8), %%mm0  \n\t" /* filterCoeff */\
+               "movq 8(%%edx), %%mm0           \n\t" /* filterCoeff */\
                "movq (%%esi, %%eax, 2), %%mm2  \n\t" /* Y1srcData */\
                "movq 8(%%esi, %%eax, 2), %%mm5 \n\t" /* Y2srcData */\
+               "addl $16, %%edx                \n\t"\
+               "movl (%%edx), %%esi            \n\t"\
                "pmulhw %%mm0, %%mm2            \n\t"\
                "pmulhw %%mm0, %%mm5            \n\t"\
                "paddw %%mm2, %%mm1             \n\t"\
                "paddw %%mm5, %%mm7             \n\t"\
-               "addl $1, %%edx                 \n\t"\
+               "testl %%esi, %%esi             \n\t"\
                " jnz 2b                        \n\t"\
-\
-               "psubw w400, %%mm3              \n\t" /* (U-128)8*/\
-               "psubw w400, %%mm4              \n\t" /* (V-128)8*/\
+
+
+#define YSCALEYUV2RGBX \
+               YSCALEYUV2PACKEDX\
+               "psubw "U_OFFSET"(%0), %%mm3    \n\t" /* (U-128)8*/\
+               "psubw "V_OFFSET"(%0), %%mm4    \n\t" /* (V-128)8*/\
                "movq %%mm3, %%mm2              \n\t" /* (U-128)8*/\
                "movq %%mm4, %%mm5              \n\t" /* (V-128)8*/\
-               "pmulhw ugCoeff, %%mm3          \n\t"\
-               "pmulhw vgCoeff, %%mm4          \n\t"\
+               "pmulhw "UG_COEFF"(%0), %%mm3   \n\t"\
+               "pmulhw "VG_COEFF"(%0), %%mm4   \n\t"\
        /* mm2=(U-128)8, mm3=ug, mm4=vg mm5=(V-128)8 */\
-               "pmulhw ubCoeff, %%mm2          \n\t"\
-               "pmulhw vrCoeff, %%mm5          \n\t"\
-               "psubw w80, %%mm1               \n\t" /* 8(Y-16)*/\
-               "psubw w80, %%mm7               \n\t" /* 8(Y-16)*/\
-               "pmulhw yCoeff, %%mm1           \n\t"\
-               "pmulhw yCoeff, %%mm7           \n\t"\
+               "pmulhw "UB_COEFF"(%0), %%mm2   \n\t"\
+               "pmulhw "VR_COEFF"(%0), %%mm5   \n\t"\
+               "psubw "Y_OFFSET"(%0), %%mm1    \n\t" /* 8(Y-16)*/\
+               "psubw "Y_OFFSET"(%0), %%mm7    \n\t" /* 8(Y-16)*/\
+               "pmulhw "Y_COEFF"(%0), %%mm1    \n\t"\
+               "pmulhw "Y_COEFF"(%0), %%mm7    \n\t"\
        /* mm1= Y1, mm2=ub, mm3=ug, mm4=vg mm5=vr, mm7=Y2 */\
                "paddw %%mm3, %%mm4             \n\t"\
                "movq %%mm2, %%mm0              \n\t"\
                "packuswb %%mm6, %%mm5          \n\t"\
                "packuswb %%mm3, %%mm4          \n\t"\
                "pxor %%mm7, %%mm7              \n\t"
-
+#if 0
 #define FULL_YSCALEYUV2RGB \
                "pxor %%mm7, %%mm7              \n\t"\
                "movd %6, %%mm6                 \n\t" /*yalpha1*/\
                "movq 4096(%3, %%eax,2), %%mm0  \n\t" /* uvbuf1[eax+2048]*/\
                "paddw %%mm2, %%mm3             \n\t" /* uvbuf0[eax]uvalpha1 - uvbuf1[eax](1-uvalpha1)*/\
                "psubw %%mm0, %%mm4             \n\t" /* uvbuf0[eax+2048] - uvbuf1[eax+2048]*/\
-               "psubw w80, %%mm1               \n\t" /* 8(Y-16)*/\
-               "psubw w400, %%mm3              \n\t" /* 8(U-128)*/\
-               "pmulhw yCoeff, %%mm1           \n\t"\
+               "psubw "MANGLE(w80)", %%mm1     \n\t" /* 8(Y-16)*/\
+               "psubw "MANGLE(w400)", %%mm3    \n\t" /* 8(U-128)*/\
+               "pmulhw "MANGLE(yCoeff)", %%mm1 \n\t"\
 \
 \
                "pmulhw %%mm5, %%mm4            \n\t" /* (uvbuf0[eax+2048] - uvbuf1[eax+2048])uvalpha1>>16*/\
                "movq %%mm3, %%mm2              \n\t" /* (U-128)8*/\
-               "pmulhw ubCoeff, %%mm3          \n\t"\
+               "pmulhw "MANGLE(ubCoeff)", %%mm3\n\t"\
                "psraw $4, %%mm0                \n\t" /* uvbuf0[eax+2048] - uvbuf1[eax+2048] >>4*/\
-               "pmulhw ugCoeff, %%mm2          \n\t"\
+               "pmulhw "MANGLE(ugCoeff)", %%mm2\n\t"\
                "paddw %%mm4, %%mm0             \n\t" /* uvbuf0[eax+2048]uvalpha1 - uvbuf1[eax+2048](1-uvalpha1)*/\
-               "psubw w400, %%mm0              \n\t" /* (V-128)8*/\
+               "psubw "MANGLE(w400)", %%mm0    \n\t" /* (V-128)8*/\
 \
 \
                "movq %%mm0, %%mm4              \n\t" /* (V-128)8*/\
-               "pmulhw vrCoeff, %%mm0          \n\t"\
-               "pmulhw vgCoeff, %%mm4          \n\t"\
+               "pmulhw "MANGLE(vrCoeff)", %%mm0\n\t"\
+               "pmulhw "MANGLE(vgCoeff)", %%mm4\n\t"\
                "paddw %%mm1, %%mm3             \n\t" /* B*/\
                "paddw %%mm1, %%mm0             \n\t" /* R*/\
                "packuswb %%mm3, %%mm3          \n\t"\
                "paddw %%mm2, %%mm1             \n\t" /* G*/\
 \
                "packuswb %%mm1, %%mm1          \n\t"
+#endif
 
-#define YSCALEYUV2RGB \
-               "movd %6, %%mm6                 \n\t" /*yalpha1*/\
-               "punpcklwd %%mm6, %%mm6         \n\t"\
-               "punpcklwd %%mm6, %%mm6         \n\t"\
-               "movq %%mm6, asm_yalpha1        \n\t"\
-               "movd %7, %%mm5                 \n\t" /*uvalpha1*/\
-               "punpcklwd %%mm5, %%mm5         \n\t"\
-               "punpcklwd %%mm5, %%mm5         \n\t"\
-               "movq %%mm5, asm_uvalpha1       \n\t"\
-               "xorl %%eax, %%eax              \n\t"\
+#define YSCALEYUV2PACKED(index, c) \
+               "movq "CHR_MMX_FILTER_OFFSET"+8("#c"), %%mm0\n\t"\
+               "movq "LUM_MMX_FILTER_OFFSET"+8("#c"), %%mm1\n\t"\
+               "psraw $3, %%mm0                \n\t"\
+               "psraw $3, %%mm1                \n\t"\
+               "movq %%mm0, "CHR_MMX_FILTER_OFFSET"+8("#c")\n\t"\
+               "movq %%mm1, "LUM_MMX_FILTER_OFFSET"+8("#c")\n\t"\
+               "xorl "#index", "#index"                \n\t"\
+               ".balign 16                     \n\t"\
+               "1:                             \n\t"\
+               "movq (%2, "#index"), %%mm2     \n\t" /* uvbuf0[eax]*/\
+               "movq (%3, "#index"), %%mm3     \n\t" /* uvbuf1[eax]*/\
+               "movq 4096(%2, "#index"), %%mm5 \n\t" /* uvbuf0[eax+2048]*/\
+               "movq 4096(%3, "#index"), %%mm4 \n\t" /* uvbuf1[eax+2048]*/\
+               "psubw %%mm3, %%mm2             \n\t" /* uvbuf0[eax] - uvbuf1[eax]*/\
+               "psubw %%mm4, %%mm5             \n\t" /* uvbuf0[eax+2048] - uvbuf1[eax+2048]*/\
+               "movq "CHR_MMX_FILTER_OFFSET"+8("#c"), %%mm0\n\t"\
+               "pmulhw %%mm0, %%mm2            \n\t" /* (uvbuf0[eax] - uvbuf1[eax])uvalpha1>>16*/\
+               "pmulhw %%mm0, %%mm5            \n\t" /* (uvbuf0[eax+2048] - uvbuf1[eax+2048])uvalpha1>>16*/\
+               "psraw $7, %%mm3                \n\t" /* uvbuf0[eax] - uvbuf1[eax] >>4*/\
+               "psraw $7, %%mm4                \n\t" /* uvbuf0[eax+2048] - uvbuf1[eax+2048] >>4*/\
+               "paddw %%mm2, %%mm3             \n\t" /* uvbuf0[eax]uvalpha1 - uvbuf1[eax](1-uvalpha1)*/\
+               "paddw %%mm5, %%mm4             \n\t" /* uvbuf0[eax+2048]uvalpha1 - uvbuf1[eax+2048](1-uvalpha1)*/\
+               "movq (%0, "#index", 2), %%mm0  \n\t" /*buf0[eax]*/\
+               "movq (%1, "#index", 2), %%mm1  \n\t" /*buf1[eax]*/\
+               "movq 8(%0, "#index", 2), %%mm6 \n\t" /*buf0[eax]*/\
+               "movq 8(%1, "#index", 2), %%mm7 \n\t" /*buf1[eax]*/\
+               "psubw %%mm1, %%mm0             \n\t" /* buf0[eax] - buf1[eax]*/\
+               "psubw %%mm7, %%mm6             \n\t" /* buf0[eax] - buf1[eax]*/\
+               "pmulhw "LUM_MMX_FILTER_OFFSET"+8("#c"), %%mm0\n\t" /* (buf0[eax] - buf1[eax])yalpha1>>16*/\
+               "pmulhw "LUM_MMX_FILTER_OFFSET"+8("#c"), %%mm6\n\t" /* (buf0[eax] - buf1[eax])yalpha1>>16*/\
+               "psraw $7, %%mm1                \n\t" /* buf0[eax] - buf1[eax] >>4*/\
+               "psraw $7, %%mm7                \n\t" /* buf0[eax] - buf1[eax] >>4*/\
+               "paddw %%mm0, %%mm1             \n\t" /* buf0[eax]yalpha1 + buf1[eax](1-yalpha1) >>16*/\
+               "paddw %%mm6, %%mm7             \n\t" /* buf0[eax]yalpha1 + buf1[eax](1-yalpha1) >>16*/\
+                
+#define YSCALEYUV2RGB(index, c) \
+               "xorl "#index", "#index"        \n\t"\
                ".balign 16                     \n\t"\
                "1:                             \n\t"\
-               "movq (%2, %%eax), %%mm2        \n\t" /* uvbuf0[eax]*/\
-               "movq (%3, %%eax), %%mm3        \n\t" /* uvbuf1[eax]*/\
-               "movq 4096(%2, %%eax), %%mm5    \n\t" /* uvbuf0[eax+2048]*/\
-               "movq 4096(%3, %%eax), %%mm4    \n\t" /* uvbuf1[eax+2048]*/\
+               "movq (%2, "#index"), %%mm2     \n\t" /* uvbuf0[eax]*/\
+               "movq (%3, "#index"), %%mm3     \n\t" /* uvbuf1[eax]*/\
+               "movq 4096(%2, "#index"), %%mm5\n\t" /* uvbuf0[eax+2048]*/\
+               "movq 4096(%3, "#index"), %%mm4\n\t" /* uvbuf1[eax+2048]*/\
                "psubw %%mm3, %%mm2             \n\t" /* uvbuf0[eax] - uvbuf1[eax]*/\
                "psubw %%mm4, %%mm5             \n\t" /* uvbuf0[eax+2048] - uvbuf1[eax+2048]*/\
-               "movq asm_uvalpha1, %%mm0       \n\t"\
+               "movq "CHR_MMX_FILTER_OFFSET"+8("#c"), %%mm0\n\t"\
                "pmulhw %%mm0, %%mm2            \n\t" /* (uvbuf0[eax] - uvbuf1[eax])uvalpha1>>16*/\
                "pmulhw %%mm0, %%mm5            \n\t" /* (uvbuf0[eax+2048] - uvbuf1[eax+2048])uvalpha1>>16*/\
                "psraw $4, %%mm3                \n\t" /* uvbuf0[eax] - uvbuf1[eax] >>4*/\
                "psraw $4, %%mm4                \n\t" /* uvbuf0[eax+2048] - uvbuf1[eax+2048] >>4*/\
                "paddw %%mm2, %%mm3             \n\t" /* uvbuf0[eax]uvalpha1 - uvbuf1[eax](1-uvalpha1)*/\
                "paddw %%mm5, %%mm4             \n\t" /* uvbuf0[eax+2048]uvalpha1 - uvbuf1[eax+2048](1-uvalpha1)*/\
-               "psubw w400, %%mm3              \n\t" /* (U-128)8*/\
-               "psubw w400, %%mm4              \n\t" /* (V-128)8*/\
+               "psubw "U_OFFSET"("#c"), %%mm3  \n\t" /* (U-128)8*/\
+               "psubw "V_OFFSET"("#c"), %%mm4  \n\t" /* (V-128)8*/\
                "movq %%mm3, %%mm2              \n\t" /* (U-128)8*/\
                "movq %%mm4, %%mm5              \n\t" /* (V-128)8*/\
-               "pmulhw ugCoeff, %%mm3          \n\t"\
-               "pmulhw vgCoeff, %%mm4          \n\t"\
+               "pmulhw "UG_COEFF"("#c"), %%mm3\n\t"\
+               "pmulhw "VG_COEFF"("#c"), %%mm4\n\t"\
        /* mm2=(U-128)8, mm3=ug, mm4=vg mm5=(V-128)8 */\
-               "movq (%0, %%eax, 2), %%mm0     \n\t" /*buf0[eax]*/\
-               "movq (%1, %%eax, 2), %%mm1     \n\t" /*buf1[eax]*/\
-               "movq 8(%0, %%eax, 2), %%mm6    \n\t" /*buf0[eax]*/\
-               "movq 8(%1, %%eax, 2), %%mm7    \n\t" /*buf1[eax]*/\
+               "movq (%0, "#index", 2), %%mm0  \n\t" /*buf0[eax]*/\
+               "movq (%1, "#index", 2), %%mm1  \n\t" /*buf1[eax]*/\
+               "movq 8(%0, "#index", 2), %%mm6\n\t" /*buf0[eax]*/\
+               "movq 8(%1, "#index", 2), %%mm7\n\t" /*buf1[eax]*/\
                "psubw %%mm1, %%mm0             \n\t" /* buf0[eax] - buf1[eax]*/\
                "psubw %%mm7, %%mm6             \n\t" /* buf0[eax] - buf1[eax]*/\
-               "pmulhw asm_yalpha1, %%mm0      \n\t" /* (buf0[eax] - buf1[eax])yalpha1>>16*/\
-               "pmulhw asm_yalpha1, %%mm6      \n\t" /* (buf0[eax] - buf1[eax])yalpha1>>16*/\
+               "pmulhw "LUM_MMX_FILTER_OFFSET"+8("#c"), %%mm0\n\t" /* (buf0[eax] - buf1[eax])yalpha1>>16*/\
+               "pmulhw "LUM_MMX_FILTER_OFFSET"+8("#c"), %%mm6\n\t" /* (buf0[eax] - buf1[eax])yalpha1>>16*/\
                "psraw $4, %%mm1                \n\t" /* buf0[eax] - buf1[eax] >>4*/\
                "psraw $4, %%mm7                \n\t" /* buf0[eax] - buf1[eax] >>4*/\
                "paddw %%mm0, %%mm1             \n\t" /* buf0[eax]yalpha1 + buf1[eax](1-yalpha1) >>16*/\
                "paddw %%mm6, %%mm7             \n\t" /* buf0[eax]yalpha1 + buf1[eax](1-yalpha1) >>16*/\
-               "pmulhw ubCoeff, %%mm2          \n\t"\
-               "pmulhw vrCoeff, %%mm5          \n\t"\
-               "psubw w80, %%mm1               \n\t" /* 8(Y-16)*/\
-               "psubw w80, %%mm7               \n\t" /* 8(Y-16)*/\
-               "pmulhw yCoeff, %%mm1           \n\t"\
-               "pmulhw yCoeff, %%mm7           \n\t"\
+               "pmulhw "UB_COEFF"("#c"), %%mm2\n\t"\
+               "pmulhw "VR_COEFF"("#c"), %%mm5\n\t"\
+               "psubw "Y_OFFSET"("#c"), %%mm1  \n\t" /* 8(Y-16)*/\
+               "psubw "Y_OFFSET"("#c"), %%mm7  \n\t" /* 8(Y-16)*/\
+               "pmulhw "Y_COEFF"("#c"), %%mm1  \n\t"\
+               "pmulhw "Y_COEFF"("#c"), %%mm7  \n\t"\
        /* mm1= Y1, mm2=ub, mm3=ug, mm4=vg mm5=vr, mm7=Y2 */\
                "paddw %%mm3, %%mm4             \n\t"\
                "movq %%mm2, %%mm0              \n\t"\
                "packuswb %%mm6, %%mm5          \n\t"\
                "packuswb %%mm3, %%mm4          \n\t"\
                "pxor %%mm7, %%mm7              \n\t"
-
-#define YSCALEYUV2RGB1 \
-               "xorl %%eax, %%eax              \n\t"\
+                
+#define YSCALEYUV2PACKED1(index, c) \
+               "xorl "#index", "#index"                \n\t"\
+               ".balign 16                     \n\t"\
+               "1:                             \n\t"\
+               "movq (%2, "#index"), %%mm3     \n\t" /* uvbuf0[eax]*/\
+               "movq 4096(%2, "#index"), %%mm4 \n\t" /* uvbuf0[eax+2048]*/\
+               "psraw $7, %%mm3                \n\t" \
+               "psraw $7, %%mm4                \n\t" \
+               "movq (%0, "#index", 2), %%mm1  \n\t" /*buf0[eax]*/\
+               "movq 8(%0, "#index", 2), %%mm7 \n\t" /*buf0[eax]*/\
+               "psraw $7, %%mm1                \n\t" \
+               "psraw $7, %%mm7                \n\t" \
+                
+#define YSCALEYUV2RGB1(index, c) \
+               "xorl "#index", "#index"        \n\t"\
                ".balign 16                     \n\t"\
                "1:                             \n\t"\
-               "movq (%2, %%eax), %%mm3        \n\t" /* uvbuf0[eax]*/\
-               "movq 4096(%2, %%eax), %%mm4    \n\t" /* uvbuf0[eax+2048]*/\
+               "movq (%2, "#index"), %%mm3     \n\t" /* uvbuf0[eax]*/\
+               "movq 4096(%2, "#index"), %%mm4 \n\t" /* uvbuf0[eax+2048]*/\
                "psraw $4, %%mm3                \n\t" /* uvbuf0[eax] - uvbuf1[eax] >>4*/\
                "psraw $4, %%mm4                \n\t" /* uvbuf0[eax+2048] - uvbuf1[eax+2048] >>4*/\
-               "psubw w400, %%mm3              \n\t" /* (U-128)8*/\
-               "psubw w400, %%mm4              \n\t" /* (V-128)8*/\
+               "psubw "U_OFFSET"("#c"), %%mm3  \n\t" /* (U-128)8*/\
+               "psubw "V_OFFSET"("#c"), %%mm4  \n\t" /* (V-128)8*/\
                "movq %%mm3, %%mm2              \n\t" /* (U-128)8*/\
                "movq %%mm4, %%mm5              \n\t" /* (V-128)8*/\
-               "pmulhw ugCoeff, %%mm3          \n\t"\
-               "pmulhw vgCoeff, %%mm4          \n\t"\
+               "pmulhw "UG_COEFF"("#c"), %%mm3\n\t"\
+               "pmulhw "VG_COEFF"("#c"), %%mm4\n\t"\
        /* mm2=(U-128)8, mm3=ug, mm4=vg mm5=(V-128)8 */\
-               "movq (%0, %%eax, 2), %%mm1     \n\t" /*buf0[eax]*/\
-               "movq 8(%0, %%eax, 2), %%mm7    \n\t" /*buf0[eax]*/\
+               "movq (%0, "#index", 2), %%mm1  \n\t" /*buf0[eax]*/\
+               "movq 8(%0, "#index", 2), %%mm7 \n\t" /*buf0[eax]*/\
                "psraw $4, %%mm1                \n\t" /* buf0[eax] - buf1[eax] >>4*/\
                "psraw $4, %%mm7                \n\t" /* buf0[eax] - buf1[eax] >>4*/\
-               "pmulhw ubCoeff, %%mm2          \n\t"\
-               "pmulhw vrCoeff, %%mm5          \n\t"\
-               "psubw w80, %%mm1               \n\t" /* 8(Y-16)*/\
-               "psubw w80, %%mm7               \n\t" /* 8(Y-16)*/\
-               "pmulhw yCoeff, %%mm1           \n\t"\
-               "pmulhw yCoeff, %%mm7           \n\t"\
+               "pmulhw "UB_COEFF"("#c"), %%mm2\n\t"\
+               "pmulhw "VR_COEFF"("#c"), %%mm5\n\t"\
+               "psubw "Y_OFFSET"("#c"), %%mm1  \n\t" /* 8(Y-16)*/\
+               "psubw "Y_OFFSET"("#c"), %%mm7  \n\t" /* 8(Y-16)*/\
+               "pmulhw "Y_COEFF"("#c"), %%mm1  \n\t"\
+               "pmulhw "Y_COEFF"("#c"), %%mm7  \n\t"\
        /* mm1= Y1, mm2=ub, mm3=ug, mm4=vg mm5=vr, mm7=Y2 */\
                "paddw %%mm3, %%mm4             \n\t"\
                "movq %%mm2, %%mm0              \n\t"\
                "packuswb %%mm3, %%mm4          \n\t"\
                "pxor %%mm7, %%mm7              \n\t"
 
+#define YSCALEYUV2PACKED1b(index, c) \
+               "xorl "#index", "#index"                \n\t"\
+               ".balign 16                     \n\t"\
+               "1:                             \n\t"\
+               "movq (%2, "#index"), %%mm2     \n\t" /* uvbuf0[eax]*/\
+               "movq (%3, "#index"), %%mm3     \n\t" /* uvbuf1[eax]*/\
+               "movq 4096(%2, "#index"), %%mm5 \n\t" /* uvbuf0[eax+2048]*/\
+               "movq 4096(%3, "#index"), %%mm4 \n\t" /* uvbuf1[eax+2048]*/\
+               "paddw %%mm2, %%mm3             \n\t" /* uvbuf0[eax] + uvbuf1[eax]*/\
+               "paddw %%mm5, %%mm4             \n\t" /* uvbuf0[eax+2048] + uvbuf1[eax+2048]*/\
+               "psrlw $8, %%mm3                \n\t" \
+               "psrlw $8, %%mm4                \n\t" \
+               "movq (%0, "#index", 2), %%mm1  \n\t" /*buf0[eax]*/\
+               "movq 8(%0, "#index", 2), %%mm7 \n\t" /*buf0[eax]*/\
+               "psraw $7, %%mm1                \n\t" \
+               "psraw $7, %%mm7                \n\t" 
+                
 // do vertical chrominance interpolation
-#define YSCALEYUV2RGB1b \
-               "xorl %%eax, %%eax              \n\t"\
+#define YSCALEYUV2RGB1b(index, c) \
+               "xorl "#index", "#index"                \n\t"\
                ".balign 16                     \n\t"\
                "1:                             \n\t"\
-               "movq (%2, %%eax), %%mm2        \n\t" /* uvbuf0[eax]*/\
-               "movq (%3, %%eax), %%mm3        \n\t" /* uvbuf1[eax]*/\
-               "movq 4096(%2, %%eax), %%mm5    \n\t" /* uvbuf0[eax+2048]*/\
-               "movq 4096(%3, %%eax), %%mm4    \n\t" /* uvbuf1[eax+2048]*/\
+               "movq (%2, "#index"), %%mm2     \n\t" /* uvbuf0[eax]*/\
+               "movq (%3, "#index"), %%mm3     \n\t" /* uvbuf1[eax]*/\
+               "movq 4096(%2, "#index"), %%mm5 \n\t" /* uvbuf0[eax+2048]*/\
+               "movq 4096(%3, "#index"), %%mm4 \n\t" /* uvbuf1[eax+2048]*/\
                "paddw %%mm2, %%mm3             \n\t" /* uvbuf0[eax] + uvbuf1[eax]*/\
                "paddw %%mm5, %%mm4             \n\t" /* uvbuf0[eax+2048] + uvbuf1[eax+2048]*/\
                "psrlw $5, %%mm3                \n\t" /*FIXME might overflow*/\
                "psrlw $5, %%mm4                \n\t" /*FIXME might overflow*/\
-               "psubw w400, %%mm3              \n\t" /* (U-128)8*/\
-               "psubw w400, %%mm4              \n\t" /* (V-128)8*/\
+               "psubw "U_OFFSET"("#c"), %%mm3  \n\t" /* (U-128)8*/\
+               "psubw "V_OFFSET"("#c"), %%mm4  \n\t" /* (V-128)8*/\
                "movq %%mm3, %%mm2              \n\t" /* (U-128)8*/\
                "movq %%mm4, %%mm5              \n\t" /* (V-128)8*/\
-               "pmulhw ugCoeff, %%mm3          \n\t"\
-               "pmulhw vgCoeff, %%mm4          \n\t"\
+               "pmulhw "UG_COEFF"("#c"), %%mm3\n\t"\
+               "pmulhw "VG_COEFF"("#c"), %%mm4\n\t"\
        /* mm2=(U-128)8, mm3=ug, mm4=vg mm5=(V-128)8 */\
-               "movq (%0, %%eax, 2), %%mm1     \n\t" /*buf0[eax]*/\
-               "movq 8(%0, %%eax, 2), %%mm7    \n\t" /*buf0[eax]*/\
+               "movq (%0, "#index", 2), %%mm1  \n\t" /*buf0[eax]*/\
+               "movq 8(%0, "#index", 2), %%mm7 \n\t" /*buf0[eax]*/\
                "psraw $4, %%mm1                \n\t" /* buf0[eax] - buf1[eax] >>4*/\
                "psraw $4, %%mm7                \n\t" /* buf0[eax] - buf1[eax] >>4*/\
-               "pmulhw ubCoeff, %%mm2          \n\t"\
-               "pmulhw vrCoeff, %%mm5          \n\t"\
-               "psubw w80, %%mm1               \n\t" /* 8(Y-16)*/\
-               "psubw w80, %%mm7               \n\t" /* 8(Y-16)*/\
-               "pmulhw yCoeff, %%mm1           \n\t"\
-               "pmulhw yCoeff, %%mm7           \n\t"\
+               "pmulhw "UB_COEFF"("#c"), %%mm2\n\t"\
+               "pmulhw "VR_COEFF"("#c"), %%mm5\n\t"\
+               "psubw "Y_OFFSET"("#c"), %%mm1  \n\t" /* 8(Y-16)*/\
+               "psubw "Y_OFFSET"("#c"), %%mm7  \n\t" /* 8(Y-16)*/\
+               "pmulhw "Y_COEFF"("#c"), %%mm1  \n\t"\
+               "pmulhw "Y_COEFF"("#c"), %%mm7  \n\t"\
        /* mm1= Y1, mm2=ub, mm3=ug, mm4=vg mm5=vr, mm7=Y2 */\
                "paddw %%mm3, %%mm4             \n\t"\
                "movq %%mm2, %%mm0              \n\t"\
                "packuswb %%mm3, %%mm4          \n\t"\
                "pxor %%mm7, %%mm7              \n\t"
 
-#define WRITEBGR32 \
+#define WRITEBGR32(dst, dstw, index) \
                /* mm2=B, %%mm4=G, %%mm5=R, %%mm7=0 */\
                        "movq %%mm2, %%mm1              \n\t" /* B */\
                        "movq %%mm5, %%mm6              \n\t" /* R */\
                        "punpcklwd %%mm6, %%mm1         \n\t" /* 0RGB0RGB 2 */\
                        "punpckhwd %%mm6, %%mm3         \n\t" /* 0RGB0RGB 3 */\
 \
-                       MOVNTQ(%%mm0, (%4, %%eax, 4))\
-                       MOVNTQ(%%mm2, 8(%4, %%eax, 4))\
-                       MOVNTQ(%%mm1, 16(%4, %%eax, 4))\
-                       MOVNTQ(%%mm3, 24(%4, %%eax, 4))\
+                       MOVNTQ(%%mm0, (dst, index, 4))\
+                       MOVNTQ(%%mm2, 8(dst, index, 4))\
+                       MOVNTQ(%%mm1, 16(dst, index, 4))\
+                       MOVNTQ(%%mm3, 24(dst, index, 4))\
 \
-                       "addl $8, %%eax                 \n\t"\
-                       "cmpl %5, %%eax                 \n\t"\
+                       "addl $8, "#index"              \n\t"\
+                       "cmpl "#dstw", "#index"         \n\t"\
                        " jb 1b                         \n\t"
 
-#define WRITEBGR16 \
-                       "pand bF8, %%mm2                \n\t" /* B */\
-                       "pand bFC, %%mm4                \n\t" /* G */\
-                       "pand bF8, %%mm5                \n\t" /* R */\
+#define WRITEBGR16(dst, dstw, index) \
+                       "pand "MANGLE(bF8)", %%mm2      \n\t" /* B */\
+                       "pand "MANGLE(bFC)", %%mm4      \n\t" /* G */\
+                       "pand "MANGLE(bF8)", %%mm5      \n\t" /* R */\
                        "psrlq $3, %%mm2                \n\t"\
 \
                        "movq %%mm2, %%mm1              \n\t"\
                        "por %%mm3, %%mm2               \n\t"\
                        "por %%mm4, %%mm1               \n\t"\
 \
-                       MOVNTQ(%%mm2, (%4, %%eax, 2))\
-                       MOVNTQ(%%mm1, 8(%4, %%eax, 2))\
+                       MOVNTQ(%%mm2, (dst, index, 2))\
+                       MOVNTQ(%%mm1, 8(dst, index, 2))\
 \
-                       "addl $8, %%eax                 \n\t"\
-                       "cmpl %5, %%eax                 \n\t"\
+                       "addl $8, "#index"              \n\t"\
+                       "cmpl "#dstw", "#index"         \n\t"\
                        " jb 1b                         \n\t"
 
-#define WRITEBGR15 \
-                       "pand bF8, %%mm2                \n\t" /* B */\
-                       "pand bF8, %%mm4                \n\t" /* G */\
-                       "pand bF8, %%mm5                \n\t" /* R */\
+#define WRITEBGR15(dst, dstw, index) \
+                       "pand "MANGLE(bF8)", %%mm2      \n\t" /* B */\
+                       "pand "MANGLE(bF8)", %%mm4      \n\t" /* G */\
+                       "pand "MANGLE(bF8)", %%mm5      \n\t" /* R */\
                        "psrlq $3, %%mm2                \n\t"\
                        "psrlq $1, %%mm5                \n\t"\
 \
                        "por %%mm3, %%mm2               \n\t"\
                        "por %%mm4, %%mm1               \n\t"\
 \
-                       MOVNTQ(%%mm2, (%4, %%eax, 2))\
-                       MOVNTQ(%%mm1, 8(%4, %%eax, 2))\
+                       MOVNTQ(%%mm2, (dst, index, 2))\
+                       MOVNTQ(%%mm1, 8(dst, index, 2))\
 \
-                       "addl $8, %%eax                 \n\t"\
-                       "cmpl %5, %%eax                 \n\t"\
+                       "addl $8, "#index"              \n\t"\
+                       "cmpl "#dstw", "#index"         \n\t"\
                        " jb 1b                         \n\t"
 
-#define WRITEBGR24OLD \
+#define WRITEBGR24OLD(dst, dstw, index) \
                /* mm2=B, %%mm4=G, %%mm5=R, %%mm7=0 */\
                        "movq %%mm2, %%mm1              \n\t" /* B */\
                        "movq %%mm5, %%mm6              \n\t" /* R */\
 \
                        "movq %%mm0, %%mm4              \n\t" /* 0RGB0RGB 0 */\
                        "psrlq $8, %%mm0                \n\t" /* 00RGB0RG 0 */\
-                       "pand bm00000111, %%mm4         \n\t" /* 00000RGB 0 */\
-                       "pand bm11111000, %%mm0         \n\t" /* 00RGB000 0.5 */\
+                       "pand "MANGLE(bm00000111)", %%mm4\n\t" /* 00000RGB 0 */\
+                       "pand "MANGLE(bm11111000)", %%mm0\n\t" /* 00RGB000 0.5 */\
                        "por %%mm4, %%mm0               \n\t" /* 00RGBRGB 0 */\
                        "movq %%mm2, %%mm4              \n\t" /* 0RGB0RGB 1 */\
                        "psllq $48, %%mm2               \n\t" /* GB000000 1 */\
                        "psrld $16, %%mm4               \n\t" /* 000R000R 1 */\
                        "psrlq $24, %%mm2               \n\t" /* 0000RGB0 1.5 */\
                        "por %%mm4, %%mm2               \n\t" /* 000RRGBR 1 */\
-                       "pand bm00001111, %%mm2         \n\t" /* 0000RGBR 1 */\
+                       "pand "MANGLE(bm00001111)", %%mm2\n\t" /* 0000RGBR 1 */\
                        "movq %%mm1, %%mm4              \n\t" /* 0RGB0RGB 2 */\
                        "psrlq $8, %%mm1                \n\t" /* 00RGB0RG 2 */\
-                       "pand bm00000111, %%mm4         \n\t" /* 00000RGB 2 */\
-                       "pand bm11111000, %%mm1         \n\t" /* 00RGB000 2.5 */\
+                       "pand "MANGLE(bm00000111)", %%mm4\n\t" /* 00000RGB 2 */\
+                       "pand "MANGLE(bm11111000)", %%mm1\n\t" /* 00RGB000 2.5 */\
                        "por %%mm4, %%mm1               \n\t" /* 00RGBRGB 2 */\
                        "movq %%mm1, %%mm4              \n\t" /* 00RGBRGB 2 */\
                        "psllq $32, %%mm1               \n\t" /* BRGB0000 2 */\
                        "psrlq $32, %%mm4               \n\t" /* 000000RG 2.5 */\
                        "movq %%mm3, %%mm5              \n\t" /* 0RGB0RGB 3 */\
                        "psrlq $8, %%mm3                \n\t" /* 00RGB0RG 3 */\
-                       "pand bm00000111, %%mm5         \n\t" /* 00000RGB 3 */\
-                       "pand bm11111000, %%mm3         \n\t" /* 00RGB000 3.5 */\
+                       "pand "MANGLE(bm00000111)", %%mm5\n\t" /* 00000RGB 3 */\
+                       "pand "MANGLE(bm11111000)", %%mm3\n\t" /* 00RGB000 3.5 */\
                        "por %%mm5, %%mm3               \n\t" /* 00RGBRGB 3 */\
                        "psllq $16, %%mm3               \n\t" /* RGBRGB00 3 */\
                        "por %%mm4, %%mm3               \n\t" /* RGBRGBRG 2.5 */\
 \
-                       MOVNTQ(%%mm0, (%%ebx))\
-                       MOVNTQ(%%mm2, 8(%%ebx))\
-                       MOVNTQ(%%mm3, 16(%%ebx))\
-                       "addl $24, %%ebx                \n\t"\
+                       MOVNTQ(%%mm0, (dst))\
+                       MOVNTQ(%%mm2, 8(dst))\
+                       MOVNTQ(%%mm3, 16(dst))\
+                       "addl $24, "#dst"               \n\t"\
 \
-                       "addl $8, %%eax                 \n\t"\
-                       "cmpl %5, %%eax                 \n\t"\
+                       "addl $8, "#index"              \n\t"\
+                       "cmpl "#dstw", "#index"         \n\t"\
                        " jb 1b                         \n\t"
 
-#define WRITEBGR24MMX \
+#define WRITEBGR24MMX(dst, dstw, index) \
                /* mm2=B, %%mm4=G, %%mm5=R, %%mm7=0 */\
                        "movq %%mm2, %%mm1              \n\t" /* B */\
                        "movq %%mm5, %%mm6              \n\t" /* R */\
                        "movq %%mm2, %%mm6              \n\t" /* 0RGBRGB0 1 */\
                        "psllq $40, %%mm2               \n\t" /* GB000000 1 */\
                        "por %%mm2, %%mm0               \n\t" /* GBRGBRGB 0 */\
-                       MOVNTQ(%%mm0, (%%ebx))\
+                       MOVNTQ(%%mm0, (dst))\
 \
                        "psrlq $24, %%mm6               \n\t" /* 0000RGBR 1 */\
                        "movq %%mm1, %%mm5              \n\t" /* 0RGBRGB0 2 */\
                        "psllq $24, %%mm1               \n\t" /* BRGB0000 2 */\
                        "por %%mm1, %%mm6               \n\t" /* BRGBRGBR 1 */\
-                       MOVNTQ(%%mm6, 8(%%ebx))\
+                       MOVNTQ(%%mm6, 8(dst))\
 \
                        "psrlq $40, %%mm5               \n\t" /* 000000RG 2 */\
                        "psllq $8, %%mm3                \n\t" /* RGBRGB00 3 */\
                        "por %%mm3, %%mm5               \n\t" /* RGBRGBRG 2 */\
-                       MOVNTQ(%%mm5, 16(%%ebx))\
+                       MOVNTQ(%%mm5, 16(dst))\
 \
-                       "addl $24, %%ebx                \n\t"\
+                       "addl $24, "#dst"               \n\t"\
 \
-                       "addl $8, %%eax                 \n\t"\
-                       "cmpl %5, %%eax                 \n\t"\
+                       "addl $8, "#index"                      \n\t"\
+                       "cmpl "#dstw", "#index"                 \n\t"\
                        " jb 1b                         \n\t"
 
-#define WRITEBGR24MMX2 \
+#define WRITEBGR24MMX2(dst, dstw, index) \
                /* mm2=B, %%mm4=G, %%mm5=R, %%mm7=0 */\
-                       "movq M24A, %%mm0               \n\t"\
-                       "movq M24C, %%mm7               \n\t"\
+                       "movq "MANGLE(M24A)", %%mm0     \n\t"\
+                       "movq "MANGLE(M24C)", %%mm7     \n\t"\
                        "pshufw $0x50, %%mm2, %%mm1     \n\t" /* B3 B2 B3 B2  B1 B0 B1 B0 */\
                        "pshufw $0x50, %%mm4, %%mm3     \n\t" /* G3 G2 G3 G2  G1 G0 G1 G0 */\
                        "pshufw $0x00, %%mm5, %%mm6     \n\t" /* R1 R0 R1 R0  R1 R0 R1 R0 */\
                        "psllq $8, %%mm3                \n\t" /* G2        G1       G0    */\
                        "por %%mm1, %%mm6               \n\t"\
                        "por %%mm3, %%mm6               \n\t"\
-                       MOVNTQ(%%mm6, (%%ebx))\
+                       MOVNTQ(%%mm6, (dst))\
 \
                        "psrlq $8, %%mm4                \n\t" /* 00 G7 G6 G5  G4 G3 G2 G1 */\
                        "pshufw $0xA5, %%mm2, %%mm1     \n\t" /* B5 B4 B5 B4  B3 B2 B3 B2 */\
                        "pshufw $0x55, %%mm4, %%mm3     \n\t" /* G4 G3 G4 G3  G4 G3 G4 G3 */\
                        "pshufw $0xA5, %%mm5, %%mm6     \n\t" /* R5 R4 R5 R4  R3 R2 R3 R2 */\
 \
-                       "pand M24B, %%mm1               \n\t" /* B5       B4        B3    */\
+                       "pand "MANGLE(M24B)", %%mm1     \n\t" /* B5       B4        B3    */\
                        "pand %%mm7, %%mm3              \n\t" /*       G4        G3       */\
                        "pand %%mm0, %%mm6              \n\t" /*    R4        R3       R2 */\
 \
                        "por %%mm1, %%mm3               \n\t" /* B5    G4 B4     G3 B3    */\
                        "por %%mm3, %%mm6               \n\t"\
-                       MOVNTQ(%%mm6, 8(%%ebx))\
+                       MOVNTQ(%%mm6, 8(dst))\
 \
                        "pshufw $0xFF, %%mm2, %%mm1     \n\t" /* B7 B6 B7 B6  B7 B6 B6 B7 */\
                        "pshufw $0xFA, %%mm4, %%mm3     \n\t" /* 00 G7 00 G7  G6 G5 G6 G5 */\
 \
                        "pand %%mm7, %%mm1              \n\t" /*       B7        B6       */\
                        "pand %%mm0, %%mm3              \n\t" /*    G7        G6       G5 */\
-                       "pand M24B, %%mm6               \n\t" /* R7       R6        R5    */\
+                       "pand "MANGLE(M24B)", %%mm6     \n\t" /* R7       R6        R5    */\
 \
                        "por %%mm1, %%mm3               \n\t"\
                        "por %%mm3, %%mm6               \n\t"\
-                       MOVNTQ(%%mm6, 16(%%ebx))\
+                       MOVNTQ(%%mm6, 16(dst))\
 \
-                       "addl $24, %%ebx                \n\t"\
+                       "addl $24, "#dst"               \n\t"\
 \
-                       "addl $8, %%eax                 \n\t"\
-                       "cmpl %5, %%eax                 \n\t"\
+                       "addl $8, "#index"              \n\t"\
+                       "cmpl "#dstw", "#index"         \n\t"\
                        " jb 1b                         \n\t"
 
 #ifdef HAVE_MMX2
 #define WRITEBGR24 WRITEBGR24MMX
 #endif
 
-static inline void RENAME(yuv2yuvX)(int16_t *lumFilter, int16_t **lumSrc, int lumFilterSize,
+#define WRITEYUY2(dst, dstw, index) \
+                       "packuswb %%mm3, %%mm3          \n\t"\
+                       "packuswb %%mm4, %%mm4          \n\t"\
+                       "packuswb %%mm7, %%mm1          \n\t"\
+                       "punpcklbw %%mm4, %%mm3         \n\t"\
+                       "movq %%mm1, %%mm7              \n\t"\
+                       "punpcklbw %%mm3, %%mm1         \n\t"\
+                       "punpckhbw %%mm3, %%mm7         \n\t"\
+\
+                       MOVNTQ(%%mm1, (dst, index, 2))\
+                       MOVNTQ(%%mm7, 8(dst, index, 2))\
+\
+                       "addl $8, "#index"              \n\t"\
+                       "cmpl "#dstw", "#index"         \n\t"\
+                       " jb 1b                         \n\t"
+
+
+static inline void RENAME(yuv2yuvX)(SwsContext *c, 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,
-                                   int16_t * lumMmxFilter, int16_t * chrMmxFilter)
+                                   uint8_t *dest, uint8_t *uDest, uint8_t *vDest, int dstW, int chrDstW)
 {
 #ifdef HAVE_MMX
        if(uDest != NULL)
        {
                asm volatile(
-                               YSCALEYUV2YV12X(0)
-                               :: "m" (-chrFilterSize), "r" (chrSrc+chrFilterSize),
-                               "r" (chrMmxFilter+chrFilterSize*4), "r" (uDest), "m" (dstW>>1)
+                               YSCALEYUV2YV12X(0, CHR_MMX_FILTER_OFFSET)
+                               :: "r" (&c->redDither),
+                               "r" (uDest), "m" (chrDstW)
                                : "%eax", "%edx", "%esi"
                        );
 
                asm volatile(
-                               YSCALEYUV2YV12X(4096)
-                               :: "m" (-chrFilterSize), "r" (chrSrc+chrFilterSize),
-                               "r" (chrMmxFilter+chrFilterSize*4), "r" (vDest), "m" (dstW>>1)
+                               YSCALEYUV2YV12X(4096, CHR_MMX_FILTER_OFFSET)
+                               :: "r" (&c->redDither),
+                               "r" (vDest), "m" (chrDstW)
                                : "%eax", "%edx", "%esi"
                        );
        }
 
        asm volatile(
-                       YSCALEYUV2YV12X(0)
-                       :: "m" (-lumFilterSize), "r" (lumSrc+lumFilterSize),
-                          "r" (lumMmxFilter+lumFilterSize*4), "r" (dest), "m" (dstW)
+                       YSCALEYUV2YV12X(0, LUM_MMX_FILTER_OFFSET)
+                       :: "r" (&c->redDither),
+                          "r" (dest), "m" (dstW)
                        : "%eax", "%edx", "%esi"
                );
 #else
 yuv2yuvXinC(lumFilter, lumSrc, lumFilterSize,
            chrFilter, chrSrc, chrFilterSize,
-           dest, uDest, vDest, dstW);
+           dest, uDest, vDest, dstW, chrDstW);
 #endif
 }
 
 static inline void RENAME(yuv2yuv1)(int16_t *lumSrc, int16_t *chrSrc,
-                                   uint8_t *dest, uint8_t *uDest, uint8_t *vDest, int dstW)
+                                   uint8_t *dest, uint8_t *uDest, uint8_t *vDest, int dstW, int chrDstW)
 {
 #ifdef HAVE_MMX
        if(uDest != NULL)
        {
                asm volatile(
                                YSCALEYUV2YV121
-                               :: "r" (chrSrc + (dstW>>1)), "r" (uDest + (dstW>>1)),
-                               "g" (-(dstW>>1))
+                               :: "r" (chrSrc + chrDstW), "r" (uDest + chrDstW),
+                               "g" (-chrDstW)
                                : "%eax"
                        );
 
                asm volatile(
                                YSCALEYUV2YV121
-                               :: "r" (chrSrc + 2048 + (dstW>>1)), "r" (vDest + (dstW>>1)),
-                               "g" (-(dstW>>1))
+                               :: "r" (chrSrc + 2048 + chrDstW), "r" (vDest + chrDstW),
+                               "g" (-chrDstW)
                                : "%eax"
                        );
        }
@@ -706,24 +800,34 @@ static inline void RENAME(yuv2yuv1)(int16_t *lumSrc, int16_t *chrSrc,
                : "%eax"
        );
 #else
-       //FIXME Optimize (just quickly writen not opti..)
-       //FIXME replace MINMAX with LUTs
        int i;
        for(i=0; i<dstW; i++)
        {
                int val= lumSrc[i]>>7;
+               
+               if(val&256){
+                       if(val<0) val=0;
+                       else      val=255;
+               }
 
-               dest[i]= MIN(MAX(val>>19, 0), 255);
+               dest[i]= val;
        }
 
        if(uDest != NULL)
-               for(i=0; i<(dstW>>1); i++)
+               for(i=0; i<chrDstW; i++)
                {
                        int u=chrSrc[i]>>7;
                        int v=chrSrc[i + 2048]>>7;
 
-                       uDest[i]= MIN(MAX(u>>19, 0), 255);
-                       vDest[i]= MIN(MAX(v>>19, 0), 255);
+                       if((u|v)&256){
+                               if(u<0)         u=0;
+                               else if (u>255) u=255;
+                               if(v<0)         v=0;
+                               else if (v>255) v=255;
+                       }
+
+                       uDest[i]= u;
+                       vDest[i]= v;
                }
 #endif
 }
@@ -732,110 +836,127 @@ static inline void RENAME(yuv2yuv1)(int16_t *lumSrc, int16_t *chrSrc,
 /**
  * vertical scale YV12 to RGB
  */
-static inline void RENAME(yuv2rgbX)(int16_t *lumFilter, int16_t **lumSrc, int lumFilterSize,
+static inline void RENAME(yuv2packedX)(SwsContext *c, int16_t *lumFilter, int16_t **lumSrc, int lumFilterSize,
                                    int16_t *chrFilter, int16_t **chrSrc, int chrFilterSize,
-                           uint8_t *dest, int dstW, int dstbpp, int16_t * lumMmxFilter, int16_t * chrMmxFilter)
+                           uint8_t *dest, int dstW, int dstY)
 {
-       if(fullUVIpol)
-       {
-//FIXME
-       }//FULL_UV_IPOL
-       else
+       int dummy=0;
+       switch(c->dstFormat)
        {
 #ifdef HAVE_MMX
-               if(dstbpp == 32) //FIXME untested
+       case IMGFMT_BGR32:
                {
                        asm volatile(
                                YSCALEYUV2RGBX
-                               WRITEBGR32
+                               WRITEBGR32(%4, %5, %%eax)
 
-                       :: "m" (-lumFilterSize), "m" (-chrFilterSize),
-                          "m" (lumMmxFilter+lumFilterSize*4), "m" (chrMmxFilter+chrFilterSize*4),
-                          "r" (dest), "m" (dstW),
-                          "m" (lumSrc+lumFilterSize), "m" (chrSrc+chrFilterSize)
-                       : "%eax", "%ebx", "%ecx", "%edx", "%esi"
+                       :: "r" (&c->redDither), 
+                          "m" (dummy), "m" (dummy), "m" (dummy),
+                          "r" (dest), "m" (dstW)
+                       : "%eax", "%edx", "%esi"
                        );
                }
-               else if(dstbpp==24) //FIXME untested
+               break;
+       case IMGFMT_BGR24:
                {
                        asm volatile(
                                YSCALEYUV2RGBX
                                "leal (%%eax, %%eax, 2), %%ebx  \n\t" //FIXME optimize
                                "addl %4, %%ebx                 \n\t"
-                               WRITEBGR24
+                               WRITEBGR24(%%ebx, %5, %%eax)
 
-                       :: "m" (-lumFilterSize), "m" (-chrFilterSize),
-                          "m" (lumMmxFilter+lumFilterSize*4), "m" (chrMmxFilter+chrFilterSize*4),
-                          "r" (dest), "m" (dstW),
-                          "m" (lumSrc+lumFilterSize), "m" (chrSrc+chrFilterSize)
-                       : "%eax", "%ebx", "%ecx", "%edx", "%esi"
+                       :: "r" (&c->redDither), 
+                          "m" (dummy), "m" (dummy), "m" (dummy),
+                          "r" (dest), "m" (dstW)
+                       : "%eax", "%ebx", "%edx", "%esi" //FIXME ebx
                        );
                }
-               else if(dstbpp==15)
+               break;
+       case IMGFMT_BGR15:
                {
                        asm volatile(
                                YSCALEYUV2RGBX
                /* mm2=B, %%mm4=G, %%mm5=R, %%mm7=0 */
 #ifdef DITHER1XBPP
-                               "paddusb b5Dither, %%mm2        \n\t"
-                               "paddusb g5Dither, %%mm4        \n\t"
-                               "paddusb r5Dither, %%mm5        \n\t"
+                               "paddusb "MANGLE(b5Dither)", %%mm2\n\t"
+                               "paddusb "MANGLE(g5Dither)", %%mm4\n\t"
+                               "paddusb "MANGLE(r5Dither)", %%mm5\n\t"
 #endif
 
-                               WRITEBGR15
+                               WRITEBGR15(%4, %5, %%eax)
 
-                       :: "m" (-lumFilterSize), "m" (-chrFilterSize),
-                          "m" (lumMmxFilter+lumFilterSize*4), "m" (chrMmxFilter+chrFilterSize*4),
-                          "r" (dest), "m" (dstW),
-                          "m" (lumSrc+lumFilterSize), "m" (chrSrc+chrFilterSize)
-                       : "%eax", "%ebx", "%ecx", "%edx", "%esi"
+                       :: "r" (&c->redDither), 
+                          "m" (dummy), "m" (dummy), "m" (dummy),
+                          "r" (dest), "m" (dstW)
+                       : "%eax", "%edx", "%esi"
                        );
                }
-               else if(dstbpp==16)
+               break;
+       case IMGFMT_BGR16:
                {
                        asm volatile(
                                YSCALEYUV2RGBX
                /* mm2=B, %%mm4=G, %%mm5=R, %%mm7=0 */
 #ifdef DITHER1XBPP
-                               "paddusb b5Dither, %%mm2        \n\t"
-                               "paddusb g6Dither, %%mm4        \n\t"
-                               "paddusb r5Dither, %%mm5        \n\t"
+                               "paddusb "MANGLE(b5Dither)", %%mm2\n\t"
+                               "paddusb "MANGLE(g6Dither)", %%mm4\n\t"
+                               "paddusb "MANGLE(r5Dither)", %%mm5\n\t"
 #endif
 
-                               WRITEBGR16
+                               WRITEBGR16(%4, %5, %%eax)
 
-                       :: "m" (-lumFilterSize), "m" (-chrFilterSize),
-                          "m" (lumMmxFilter+lumFilterSize*4), "m" (chrMmxFilter+chrFilterSize*4),
-                          "r" (dest), "m" (dstW),
-                          "m" (lumSrc+lumFilterSize), "m" (chrSrc+chrFilterSize)
-                       : "%eax", "%ebx", "%ecx", "%edx", "%esi"
+                       :: "r" (&c->redDither), 
+                          "m" (dummy), "m" (dummy), "m" (dummy),
+                          "r" (dest), "m" (dstW)
+                       : "%eax", "%edx", "%esi"
                        );
                }
-#else
-yuv2rgbXinC(lumFilter, lumSrc, lumFilterSize,
-           chrFilter, chrSrc, chrFilterSize,
-           dest, dstW, dstbpp);
+               break;
+       case IMGFMT_YUY2:
+               {
+                       asm volatile(
+                               YSCALEYUV2PACKEDX
+               /* mm2=B, %%mm4=G, %%mm5=R, %%mm7=0 */
+
+                               "psraw $3, %%mm3                \n\t"
+                               "psraw $3, %%mm4                \n\t"
+                               "psraw $3, %%mm1                \n\t"
+                               "psraw $3, %%mm7                \n\t"
+                               WRITEYUY2(%4, %5, %%eax)
 
+                       :: "r" (&c->redDither), 
+                          "m" (dummy), "m" (dummy), "m" (dummy),
+                          "r" (dest), "m" (dstW)
+                       : "%eax", "%edx", "%esi"
+                       );
+               }
+               break;
 #endif
-       } //!FULL_UV_IPOL
+       default:
+               yuv2packedXinC(c, lumFilter, lumSrc, lumFilterSize,
+                           chrFilter, chrSrc, chrFilterSize,
+                           dest, dstW, dstY);
+               break;
+       }
 }
 
-
 /**
  * vertical bilinear scale YV12 to RGB
  */
-static inline void RENAME(yuv2rgb2)(uint16_t *buf0, uint16_t *buf1, uint16_t *uvbuf0, uint16_t *uvbuf1,
-                           uint8_t *dest, int dstW, int yalpha, int uvalpha, int dstbpp)
+static inline void RENAME(yuv2packed2)(SwsContext *c, uint16_t *buf0, uint16_t *buf1, uint16_t *uvbuf0, uint16_t *uvbuf1,
+                           uint8_t *dest, int dstW, int yalpha, int uvalpha, int y)
 {
        int yalpha1=yalpha^4095;
        int uvalpha1=uvalpha^4095;
+       int i;
 
-       if(fullUVIpol)
+#if 0 //isn't used
+       if(flags&SWS_FULL_CHR_H_INT)
        {
-
-#ifdef HAVE_MMX
-               if(dstbpp == 32)
+               switch(dstFormat)
                {
+#ifdef HAVE_MMX
+               case IMGFMT_BGR32:
                        asm volatile(
 
 
@@ -859,9 +980,8 @@ FULL_YSCALEYUV2RGB
                        "m" (yalpha1), "m" (uvalpha1)
                        : "%eax"
                        );
-               }
-               else if(dstbpp==24)
-               {
+                       break;
+               case IMGFMT_BGR24:
                        asm volatile(
 
 FULL_YSCALEYUV2RGB
@@ -876,8 +996,8 @@ FULL_YSCALEYUV2RGB
 
                        "movq %%mm3, %%mm2              \n\t" // BGR0BGR0
                        "psrlq $8, %%mm3                \n\t" // GR0BGR00
-                       "pand bm00000111, %%mm2         \n\t" // BGR00000
-                       "pand bm11111000, %%mm3         \n\t" // 000BGR00
+                       "pand "MANGLE(bm00000111)", %%mm2\n\t" // BGR00000
+                       "pand "MANGLE(bm11111000)", %%mm3\n\t" // 000BGR00
                        "por %%mm2, %%mm3               \n\t" // BGRBGR00
                        "movq %%mm1, %%mm2              \n\t"
                        "psllq $48, %%mm1               \n\t" // 000000BG
@@ -909,16 +1029,15 @@ FULL_YSCALEYUV2RGB
                        "m" (yalpha1), "m" (uvalpha1)
                        : "%eax", "%ebx"
                        );
-               }
-               else if(dstbpp==15)
-               {
+                       break;
+               case IMGFMT_BGR15:
                        asm volatile(
 
 FULL_YSCALEYUV2RGB
 #ifdef DITHER1XBPP
-                       "paddusb g5Dither, %%mm1        \n\t"
-                       "paddusb r5Dither, %%mm0        \n\t"
-                       "paddusb b5Dither, %%mm3        \n\t"
+                       "paddusb "MANGLE(g5Dither)", %%mm1\n\t"
+                       "paddusb "MANGLE(r5Dither)", %%mm0\n\t"
+                       "paddusb "MANGLE(b5Dither)", %%mm3\n\t"
 #endif
                        "punpcklbw %%mm7, %%mm1         \n\t" // 0G0G0G0G
                        "punpcklbw %%mm7, %%mm3         \n\t" // 0B0B0B0B
@@ -927,8 +1046,8 @@ FULL_YSCALEYUV2RGB
                        "psrlw $3, %%mm3                \n\t"
                        "psllw $2, %%mm1                \n\t"
                        "psllw $7, %%mm0                \n\t"
-                       "pand g15Mask, %%mm1            \n\t"
-                       "pand r15Mask, %%mm0            \n\t"
+                       "pand "MANGLE(g15Mask)", %%mm1  \n\t"
+                       "pand "MANGLE(r15Mask)", %%mm0  \n\t"
 
                        "por %%mm3, %%mm1               \n\t"
                        "por %%mm1, %%mm0               \n\t"
@@ -943,16 +1062,15 @@ FULL_YSCALEYUV2RGB
                        "m" (yalpha1), "m" (uvalpha1)
                        : "%eax"
                        );
-               }
-               else if(dstbpp==16)
-               {
+                       break;
+               case IMGFMT_BGR16:
                        asm volatile(
 
 FULL_YSCALEYUV2RGB
 #ifdef DITHER1XBPP
-                       "paddusb g6Dither, %%mm1        \n\t"
-                       "paddusb r5Dither, %%mm0        \n\t"
-                       "paddusb b5Dither, %%mm3        \n\t"
+                       "paddusb "MANGLE(g6Dither)", %%mm1\n\t"
+                       "paddusb "MANGLE(r5Dither)", %%mm0\n\t"
+                       "paddusb "MANGLE(b5Dither)", %%mm3\n\t"
 #endif
                        "punpcklbw %%mm7, %%mm1         \n\t" // 0G0G0G0G
                        "punpcklbw %%mm7, %%mm3         \n\t" // 0B0B0B0B
@@ -961,8 +1079,8 @@ FULL_YSCALEYUV2RGB
                        "psrlw $3, %%mm3                \n\t"
                        "psllw $3, %%mm1                \n\t"
                        "psllw $8, %%mm0                \n\t"
-                       "pand g16Mask, %%mm1            \n\t"
-                       "pand r16Mask, %%mm0            \n\t"
+                       "pand "MANGLE(g16Mask)", %%mm1  \n\t"
+                       "pand "MANGLE(r16Mask)", %%mm0  \n\t"
 
                        "por %%mm3, %%mm1               \n\t"
                        "por %%mm1, %%mm0               \n\t"
@@ -977,9 +1095,30 @@ FULL_YSCALEYUV2RGB
                        "m" (yalpha1), "m" (uvalpha1)
                        : "%eax"
                        );
+               break;
+#endif
+               case IMGFMT_RGB32:
+#ifndef HAVE_MMX
+               case IMGFMT_BGR32:
+#endif
+               if(dstFormat==IMGFMT_BGR32)
+               {
+                       int i;
+#ifdef WORDS_BIGENDIAN
+                       dest++;
+#endif
+                       for(i=0;i<dstW;i++){
+                               // vertical linear interpolation && yuv2rgb in a single step:
+                               int Y=yuvtab_2568[((buf0[i]*yalpha1+buf1[i]*yalpha)>>19)];
+                               int U=((uvbuf0[i]*uvalpha1+uvbuf1[i]*uvalpha)>>19);
+                               int V=((uvbuf0[i+2048]*uvalpha1+uvbuf1[i+2048]*uvalpha)>>19);
+                               dest[0]=clip_table[((Y + yuvtab_40cf[U]) >>13)];
+                               dest[1]=clip_table[((Y + yuvtab_1a1e[V] + yuvtab_0c92[U]) >>13)];
+                               dest[2]=clip_table[((Y + yuvtab_3343[V]) >>13)];
+                               dest+= 4;
+                       }
                }
-#else
-               if(dstbpp==32 || dstbpp==24)
+               else if(dstFormat==IMGFMT_BGR24)
                {
                        int i;
                        for(i=0;i<dstW;i++){
@@ -990,10 +1129,10 @@ FULL_YSCALEYUV2RGB
                                dest[0]=clip_table[((Y + yuvtab_40cf[U]) >>13)];
                                dest[1]=clip_table[((Y + yuvtab_1a1e[V] + yuvtab_0c92[U]) >>13)];
                                dest[2]=clip_table[((Y + yuvtab_3343[V]) >>13)];
-                               dest+=dstbpp>>3;
+                               dest+= 3;
                        }
                }
-               else if(dstbpp==16)
+               else if(dstFormat==IMGFMT_BGR16)
                {
                        int i;
                        for(i=0;i<dstW;i++){
@@ -1008,7 +1147,7 @@ FULL_YSCALEYUV2RGB
                                        clip_table16r[(Y + yuvtab_3343[V]) >>13];
                        }
                }
-               else if(dstbpp==15)
+               else if(dstFormat==IMGFMT_BGR15)
                {
                        int i;
                        for(i=0;i<dstW;i++){
@@ -1023,419 +1162,859 @@ FULL_YSCALEYUV2RGB
                                        clip_table15r[(Y + yuvtab_3343[V]) >>13];
                        }
                }
-#endif
        }//FULL_UV_IPOL
        else
        {
+#endif // if 0
 #ifdef HAVE_MMX
-               if(dstbpp == 32)
-               {
+       switch(c->dstFormat)
+       {
+//Note 8280 == DSTW_OFFSET but the preprocessor can't handle that there :(
+       case IMGFMT_BGR32:
                        asm volatile(
-                               YSCALEYUV2RGB
-                               WRITEBGR32
-
-                       :: "r" (buf0), "r" (buf1), "r" (uvbuf0), "r" (uvbuf1), "r" (dest), "m" (dstW),
-                       "m" (yalpha1), "m" (uvalpha1)
+                               "movl %%esp, "ESP_OFFSET"(%5)           \n\t"
+                               "movl %4, %%esp                         \n\t"
+                               YSCALEYUV2RGB(%%eax, %5)
+                               WRITEBGR32(%%esp, 8280(%5), %%eax)
+                               "movl "ESP_OFFSET"(%5), %%esp           \n\t"
+
+                       :: "r" (buf0), "r" (buf1), "r" (uvbuf0), "r" (uvbuf1), "m" (dest),
+                       "r" (&c->redDither)
                        : "%eax"
                        );
-               }
-               else if(dstbpp==24)
-               {
+                       return;
+       case IMGFMT_BGR24:
                        asm volatile(
-                               "movl %4, %%ebx                 \n\t"
-                               YSCALEYUV2RGB
-                               WRITEBGR24
-
-                       :: "r" (buf0), "r" (buf1), "r" (uvbuf0), "r" (uvbuf1), "m" (dest), "m" (dstW),
-                       "m" (yalpha1), "m" (uvalpha1)
-                       : "%eax", "%ebx"
+                               "movl %%esp, "ESP_OFFSET"(%5)           \n\t"
+                               "movl %4, %%esp                 \n\t"
+                               YSCALEYUV2RGB(%%eax, %5)
+                               WRITEBGR24(%%esp, 8280(%5), %%eax)
+                               "movl "ESP_OFFSET"(%5), %%esp           \n\t"
+                       :: "r" (buf0), "r" (buf1), "r" (uvbuf0), "r" (uvbuf1), "m" (dest),
+                       "r" (&c->redDither)
+                       : "%eax"
                        );
-               }
-               else if(dstbpp==15)
-               {
+                       return;
+       case IMGFMT_BGR15:
                        asm volatile(
-                               YSCALEYUV2RGB
+                               "movl %%esp, "ESP_OFFSET"(%5)           \n\t"
+                               "movl %4, %%esp                         \n\t"
+                               YSCALEYUV2RGB(%%eax, %5)
                /* mm2=B, %%mm4=G, %%mm5=R, %%mm7=0 */
 #ifdef DITHER1XBPP
-                               "paddusb b5Dither, %%mm2        \n\t"
-                               "paddusb g5Dither, %%mm4        \n\t"
-                               "paddusb r5Dither, %%mm5        \n\t"
+                               "paddusb "MANGLE(b5Dither)", %%mm2\n\t"
+                               "paddusb "MANGLE(g5Dither)", %%mm4\n\t"
+                               "paddusb "MANGLE(r5Dither)", %%mm5\n\t"
 #endif
 
-                               WRITEBGR15
+                               WRITEBGR15(%%esp, 8280(%5), %%eax)
+                               "movl "ESP_OFFSET"(%5), %%esp           \n\t"
 
-                       :: "r" (buf0), "r" (buf1), "r" (uvbuf0), "r" (uvbuf1), "r" (dest), "m" (dstW),
-                       "m" (yalpha1), "m" (uvalpha1)
+                       :: "r" (buf0), "r" (buf1), "r" (uvbuf0), "r" (uvbuf1), "m" (dest),
+                       "r" (&c->redDither)
                        : "%eax"
                        );
-               }
-               else if(dstbpp==16)
-               {
+                       return;
+       case IMGFMT_BGR16:
                        asm volatile(
-                               YSCALEYUV2RGB
+                               "movl %%esp, "ESP_OFFSET"(%5)           \n\t"
+                               "movl %4, %%esp                         \n\t"
+                               YSCALEYUV2RGB(%%eax, %5)
                /* mm2=B, %%mm4=G, %%mm5=R, %%mm7=0 */
 #ifdef DITHER1XBPP
-                               "paddusb b5Dither, %%mm2        \n\t"
-                               "paddusb g6Dither, %%mm4        \n\t"
-                               "paddusb r5Dither, %%mm5        \n\t"
+                               "paddusb "MANGLE(b5Dither)", %%mm2\n\t"
+                               "paddusb "MANGLE(g6Dither)", %%mm4\n\t"
+                               "paddusb "MANGLE(r5Dither)", %%mm5\n\t"
 #endif
 
-                               WRITEBGR16
-
-                       :: "r" (buf0), "r" (buf1), "r" (uvbuf0), "r" (uvbuf1), "r" (dest), "m" (dstW),
-                       "m" (yalpha1), "m" (uvalpha1)
+                               WRITEBGR16(%%esp, 8280(%5), %%eax)
+                               "movl "ESP_OFFSET"(%5), %%esp           \n\t"
+                       :: "r" (buf0), "r" (buf1), "r" (uvbuf0), "r" (uvbuf1), "m" (dest),
+                       "r" (&c->redDither)
                        : "%eax"
                        );
-               }
-#else
-               if(dstbpp==32)
-               {
-                       int i;
-                       for(i=0; i<dstW-1; i+=2){
-                               // vertical linear interpolation && yuv2rgb in a single step:
-                               int Y1=yuvtab_2568[((buf0[i]*yalpha1+buf1[i]*yalpha)>>19)];
-                               int Y2=yuvtab_2568[((buf0[i+1]*yalpha1+buf1[i+1]*yalpha)>>19)];
-                               int U=((uvbuf0[i>>1]*uvalpha1+uvbuf1[i>>1]*uvalpha)>>19);
-                               int V=((uvbuf0[(i>>1)+2048]*uvalpha1+uvbuf1[(i>>1)+2048]*uvalpha)>>19);
-
-                               int Cb= yuvtab_40cf[U];
-                               int Cg= yuvtab_1a1e[V] + yuvtab_0c92[U];
-                               int Cr= yuvtab_3343[V];
-
-                               dest[4*i+0]=clip_table[((Y1 + Cb) >>13)];
-                               dest[4*i+1]=clip_table[((Y1 + Cg) >>13)];
-                               dest[4*i+2]=clip_table[((Y1 + Cr) >>13)];
-
-                               dest[4*i+4]=clip_table[((Y2 + Cb) >>13)];
-                               dest[4*i+5]=clip_table[((Y2 + Cg) >>13)];
-                               dest[4*i+6]=clip_table[((Y2 + Cr) >>13)];
-                       }
-               }
-               else if(dstbpp==24)
-               {
-                       int i;
-                       for(i=0; i<dstW-1; i+=2){
-                               // vertical linear interpolation && yuv2rgb in a single step:
-                               int Y1=yuvtab_2568[((buf0[i]*yalpha1+buf1[i]*yalpha)>>19)];
-                               int Y2=yuvtab_2568[((buf0[i+1]*yalpha1+buf1[i+1]*yalpha)>>19)];
-                               int U=((uvbuf0[i>>1]*uvalpha1+uvbuf1[i>>1]*uvalpha)>>19);
-                               int V=((uvbuf0[(i>>1)+2048]*uvalpha1+uvbuf1[(i>>1)+2048]*uvalpha)>>19);
-
-                               int Cb= yuvtab_40cf[U];
-                               int Cg= yuvtab_1a1e[V] + yuvtab_0c92[U];
-                               int Cr= yuvtab_3343[V];
-
-                               dest[0]=clip_table[((Y1 + Cb) >>13)];
-                               dest[1]=clip_table[((Y1 + Cg) >>13)];
-                               dest[2]=clip_table[((Y1 + Cr) >>13)];
-
-                               dest[3]=clip_table[((Y2 + Cb) >>13)];
-                               dest[4]=clip_table[((Y2 + Cg) >>13)];
-                               dest[5]=clip_table[((Y2 + Cr) >>13)];
-                               dest+=6;
-                       }
-               }
-               else if(dstbpp==16)
-               {
-                       int i;
-                       for(i=0; i<dstW-1; i+=2){
-                               // vertical linear interpolation && yuv2rgb in a single step:
-                               int Y1=yuvtab_2568[((buf0[i]*yalpha1+buf1[i]*yalpha)>>19)];
-                               int Y2=yuvtab_2568[((buf0[i+1]*yalpha1+buf1[i+1]*yalpha)>>19)];
-                               int U=((uvbuf0[i>>1]*uvalpha1+uvbuf1[i>>1]*uvalpha)>>19);
-                               int V=((uvbuf0[(i>>1)+2048]*uvalpha1+uvbuf1[(i>>1)+2048]*uvalpha)>>19);
-
-                               int Cb= yuvtab_40cf[U];
-                               int Cg= yuvtab_1a1e[V] + yuvtab_0c92[U];
-                               int Cr= yuvtab_3343[V];
-
-                               ((uint16_t*)dest)[i] =
-                                       clip_table16b[(Y1 + Cb) >>13] |
-                                       clip_table16g[(Y1 + Cg) >>13] |
-                                       clip_table16r[(Y1 + Cr) >>13];
-
-                               ((uint16_t*)dest)[i+1] =
-                                       clip_table16b[(Y2 + Cb) >>13] |
-                                       clip_table16g[(Y2 + Cg) >>13] |
-                                       clip_table16r[(Y2 + Cr) >>13];
-                       }
-               }
-               else if(dstbpp==15)
-               {
-                       int i;
-                       for(i=0; i<dstW-1; i+=2){
-                               // vertical linear interpolation && yuv2rgb in a single step:
-                               int Y1=yuvtab_2568[((buf0[i]*yalpha1+buf1[i]*yalpha)>>19)];
-                               int Y2=yuvtab_2568[((buf0[i+1]*yalpha1+buf1[i+1]*yalpha)>>19)];
-                               int U=((uvbuf0[i>>1]*uvalpha1+uvbuf1[i>>1]*uvalpha)>>19);
-                               int V=((uvbuf0[(i>>1)+2048]*uvalpha1+uvbuf1[(i>>1)+2048]*uvalpha)>>19);
-
-                               int Cb= yuvtab_40cf[U];
-                               int Cg= yuvtab_1a1e[V] + yuvtab_0c92[U];
-                               int Cr= yuvtab_3343[V];
-
-                               ((uint16_t*)dest)[i] =
-                                       clip_table15b[(Y1 + Cb) >>13] |
-                                       clip_table15g[(Y1 + Cg) >>13] |
-                                       clip_table15r[(Y1 + Cr) >>13];
-
-                               ((uint16_t*)dest)[i+1] =
-                                       clip_table15b[(Y2 + Cb) >>13] |
-                                       clip_table15g[(Y2 + Cg) >>13] |
-                                       clip_table15r[(Y2 + Cr) >>13];
-                       }
-               }
-#endif
-       } //!FULL_UV_IPOL
+                       return;
+       case IMGFMT_YUY2:
+                       asm volatile(
+                               "movl %%esp, "ESP_OFFSET"(%5)           \n\t"
+                               "movl %4, %%esp                         \n\t"
+                               YSCALEYUV2PACKED(%%eax, %5)
+                               WRITEYUY2(%%esp, 8280(%5), %%eax)
+                               "movl "ESP_OFFSET"(%5), %%esp           \n\t"
+                       :: "r" (buf0), "r" (buf1), "r" (uvbuf0), "r" (uvbuf1), "m" (dest),
+                       "r" (&c->redDither)
+                       : "%eax"
+                       );
+                       return;
+       default: break;
+       }
+#endif //HAVE_MMX
+YSCALE_YUV_2_ANYRGB_C(YSCALE_YUV_2_RGB2_C, YSCALE_YUV_2_PACKED2_C)
 }
 
 /**
  * YV12 to RGB without scaling or interpolating
  */
-static inline void RENAME(yuv2rgb1)(uint16_t *buf0, uint16_t *uvbuf0, uint16_t *uvbuf1,
-                           uint8_t *dest, int dstW, int uvalpha, int dstbpp)
+static inline void RENAME(yuv2packed1)(SwsContext *c, uint16_t *buf0, uint16_t *uvbuf0, uint16_t *uvbuf1,
+                           uint8_t *dest, int dstW, int uvalpha, int dstFormat, int flags, int y)
 {
-       int uvalpha1=uvalpha^4095;
        const int yalpha1=0;
+       int i;
+       
+       uint16_t *buf1= buf0; //FIXME needed for the rgb1/bgr1
+       const int yalpha= 4096; //FIXME ...
 
-       if(fullUVIpol || allwaysIpol)
+       if(flags&SWS_FULL_CHR_H_INT)
        {
-               RENAME(yuv2rgb2)(buf0, buf0, uvbuf0, uvbuf1, dest, dstW, 0, uvalpha, dstbpp);
+               RENAME(yuv2packed2)(c, buf0, buf0, uvbuf0, uvbuf1, dest, dstW, 0, uvalpha, y);
                return;
        }
 
 #ifdef HAVE_MMX
        if( uvalpha < 2048 ) // note this is not correct (shifts chrominance by 0.5 pixels) but its a bit faster
        {
-               if(dstbpp == 32)
+               switch(dstFormat)
                {
+               case IMGFMT_BGR32:
                        asm volatile(
-                               YSCALEYUV2RGB1
-                               WRITEBGR32
-                       :: "r" (buf0), "r" (buf0), "r" (uvbuf0), "r" (uvbuf1), "r" (dest), "m" (dstW),
-                       "m" (yalpha1), "m" (uvalpha1)
+                               "movl %%esp, "ESP_OFFSET"(%5)           \n\t"
+                               "movl %4, %%esp                         \n\t"
+                               YSCALEYUV2RGB1(%%eax, %5)
+                               WRITEBGR32(%%esp, 8280(%5), %%eax)
+                               "movl "ESP_OFFSET"(%5), %%esp           \n\t"
+
+                       :: "r" (buf0), "r" (buf1), "r" (uvbuf0), "r" (uvbuf1), "m" (dest),
+                       "r" (&c->redDither)
                        : "%eax"
                        );
-               }
-               else if(dstbpp==24)
-               {
+                       return;
+               case IMGFMT_BGR24:
                        asm volatile(
-                               "movl %4, %%ebx                 \n\t"
-                               YSCALEYUV2RGB1
-                               WRITEBGR24
-                       :: "r" (buf0), "r" (buf0), "r" (uvbuf0), "r" (uvbuf1), "m" (dest), "m" (dstW),
-                       "m" (yalpha1), "m" (uvalpha1)
-                       : "%eax", "%ebx"
+                               "movl %%esp, "ESP_OFFSET"(%5)           \n\t"
+                               "movl %4, %%esp                         \n\t"
+                               YSCALEYUV2RGB1(%%eax, %5)
+                               WRITEBGR24(%%esp, 8280(%5), %%eax)
+                               "movl "ESP_OFFSET"(%5), %%esp           \n\t"
+
+                       :: "r" (buf0), "r" (buf1), "r" (uvbuf0), "r" (uvbuf1), "m" (dest),
+                       "r" (&c->redDither)
+                       : "%eax"
                        );
-               }
-               else if(dstbpp==15)
-               {
+                       return;
+               case IMGFMT_BGR15:
                        asm volatile(
-                               YSCALEYUV2RGB1
+                               "movl %%esp, "ESP_OFFSET"(%5)           \n\t"
+                               "movl %4, %%esp                         \n\t"
+                               YSCALEYUV2RGB1(%%eax, %5)
                /* mm2=B, %%mm4=G, %%mm5=R, %%mm7=0 */
 #ifdef DITHER1XBPP
-                               "paddusb b5Dither, %%mm2        \n\t"
-                               "paddusb g5Dither, %%mm4        \n\t"
-                               "paddusb r5Dither, %%mm5        \n\t"
+                               "paddusb "MANGLE(b5Dither)", %%mm2\n\t"
+                               "paddusb "MANGLE(g5Dither)", %%mm4\n\t"
+                               "paddusb "MANGLE(r5Dither)", %%mm5\n\t"
 #endif
-                               WRITEBGR15
-                       :: "r" (buf0), "r" (buf0), "r" (uvbuf0), "r" (uvbuf1), "r" (dest), "m" (dstW),
-                       "m" (yalpha1), "m" (uvalpha1)
+                               WRITEBGR15(%%esp, 8280(%5), %%eax)
+                               "movl "ESP_OFFSET"(%5), %%esp           \n\t"
+
+                       :: "r" (buf0), "r" (buf1), "r" (uvbuf0), "r" (uvbuf1), "m" (dest),
+                       "r" (&c->redDither)
                        : "%eax"
                        );
-               }
-               else if(dstbpp==16)
-               {
+                       return;
+               case IMGFMT_BGR16:
                        asm volatile(
-                               YSCALEYUV2RGB1
+                               "movl %%esp, "ESP_OFFSET"(%5)           \n\t"
+                               "movl %4, %%esp                         \n\t"
+                               YSCALEYUV2RGB1(%%eax, %5)
                /* mm2=B, %%mm4=G, %%mm5=R, %%mm7=0 */
 #ifdef DITHER1XBPP
-                               "paddusb b5Dither, %%mm2        \n\t"
-                               "paddusb g6Dither, %%mm4        \n\t"
-                               "paddusb r5Dither, %%mm5        \n\t"
+                               "paddusb "MANGLE(b5Dither)", %%mm2\n\t"
+                               "paddusb "MANGLE(g6Dither)", %%mm4\n\t"
+                               "paddusb "MANGLE(r5Dither)", %%mm5\n\t"
 #endif
 
-                               WRITEBGR16
-                       :: "r" (buf0), "r" (buf0), "r" (uvbuf0), "r" (uvbuf1), "r" (dest), "m" (dstW),
-                       "m" (yalpha1), "m" (uvalpha1)
+                               WRITEBGR16(%%esp, 8280(%5), %%eax)
+                               "movl "ESP_OFFSET"(%5), %%esp           \n\t"
+
+                       :: "r" (buf0), "r" (buf1), "r" (uvbuf0), "r" (uvbuf1), "m" (dest),
+                       "r" (&c->redDither)
+                       : "%eax"
+                       );
+                       return;
+               case IMGFMT_YUY2:
+                       asm volatile(
+                               "movl %%esp, "ESP_OFFSET"(%5)           \n\t"
+                               "movl %4, %%esp                         \n\t"
+                               YSCALEYUV2PACKED1(%%eax, %5)
+                               WRITEYUY2(%%esp, 8280(%5), %%eax)
+                               "movl "ESP_OFFSET"(%5), %%esp           \n\t"
+
+                       :: "r" (buf0), "r" (buf1), "r" (uvbuf0), "r" (uvbuf1), "m" (dest),
+                       "r" (&c->redDither)
                        : "%eax"
                        );
+                       return;
                }
        }
        else
        {
-               if(dstbpp == 32)
+               switch(dstFormat)
                {
+               case IMGFMT_BGR32:
                        asm volatile(
-                               YSCALEYUV2RGB1b
-                               WRITEBGR32
-                       :: "r" (buf0), "r" (buf0), "r" (uvbuf0), "r" (uvbuf1), "r" (dest), "m" (dstW),
-                       "m" (yalpha1), "m" (uvalpha1)
+                               "movl %%esp, "ESP_OFFSET"(%5)           \n\t"
+                               "movl %4, %%esp                         \n\t"
+                               YSCALEYUV2RGB1b(%%eax, %5)
+                               WRITEBGR32(%%esp, 8280(%5), %%eax)
+                               "movl "ESP_OFFSET"(%5), %%esp           \n\t"
+
+                       :: "r" (buf0), "r" (buf1), "r" (uvbuf0), "r" (uvbuf1), "m" (dest),
+                       "r" (&c->redDither)
                        : "%eax"
                        );
-               }
-               else if(dstbpp==24)
-               {
+                       return;
+               case IMGFMT_BGR24:
                        asm volatile(
-                               "movl %4, %%ebx                 \n\t"
-                               YSCALEYUV2RGB1b
-                               WRITEBGR24
-                       :: "r" (buf0), "r" (buf0), "r" (uvbuf0), "r" (uvbuf1), "m" (dest), "m" (dstW),
-                       "m" (yalpha1), "m" (uvalpha1)
-                       : "%eax", "%ebx"
+                               "movl %%esp, "ESP_OFFSET"(%5)           \n\t"
+                               "movl %4, %%esp                         \n\t"
+                               YSCALEYUV2RGB1b(%%eax, %5)
+                               WRITEBGR24(%%esp, 8280(%5), %%eax)
+                               "movl "ESP_OFFSET"(%5), %%esp           \n\t"
+
+                       :: "r" (buf0), "r" (buf1), "r" (uvbuf0), "r" (uvbuf1), "m" (dest),
+                       "r" (&c->redDither)
+                       : "%eax"
                        );
-               }
-               else if(dstbpp==15)
-               {
+                       return;
+               case IMGFMT_BGR15:
                        asm volatile(
-                               YSCALEYUV2RGB1b
+                               "movl %%esp, "ESP_OFFSET"(%5)           \n\t"
+                               "movl %4, %%esp                         \n\t"
+                               YSCALEYUV2RGB1b(%%eax, %5)
                /* mm2=B, %%mm4=G, %%mm5=R, %%mm7=0 */
 #ifdef DITHER1XBPP
-                               "paddusb b5Dither, %%mm2        \n\t"
-                               "paddusb g5Dither, %%mm4        \n\t"
-                               "paddusb r5Dither, %%mm5        \n\t"
+                               "paddusb "MANGLE(b5Dither)", %%mm2\n\t"
+                               "paddusb "MANGLE(g5Dither)", %%mm4\n\t"
+                               "paddusb "MANGLE(r5Dither)", %%mm5\n\t"
 #endif
-                               WRITEBGR15
-                       :: "r" (buf0), "r" (buf0), "r" (uvbuf0), "r" (uvbuf1), "r" (dest), "m" (dstW),
-                       "m" (yalpha1), "m" (uvalpha1)
+                               WRITEBGR15(%%esp, 8280(%5), %%eax)
+                               "movl "ESP_OFFSET"(%5), %%esp           \n\t"
+
+                       :: "r" (buf0), "r" (buf1), "r" (uvbuf0), "r" (uvbuf1), "m" (dest),
+                       "r" (&c->redDither)
                        : "%eax"
                        );
-               }
-               else if(dstbpp==16)
-               {
+                       return;
+               case IMGFMT_BGR16:
                        asm volatile(
-                               YSCALEYUV2RGB1b
+                               "movl %%esp, "ESP_OFFSET"(%5)           \n\t"
+                               "movl %4, %%esp                         \n\t"
+                               YSCALEYUV2RGB1b(%%eax, %5)
                /* mm2=B, %%mm4=G, %%mm5=R, %%mm7=0 */
 #ifdef DITHER1XBPP
-                               "paddusb b5Dither, %%mm2        \n\t"
-                               "paddusb g6Dither, %%mm4        \n\t"
-                               "paddusb r5Dither, %%mm5        \n\t"
+                               "paddusb "MANGLE(b5Dither)", %%mm2\n\t"
+                               "paddusb "MANGLE(g6Dither)", %%mm4\n\t"
+                               "paddusb "MANGLE(r5Dither)", %%mm5\n\t"
 #endif
 
-                               WRITEBGR16
-                       :: "r" (buf0), "r" (buf0), "r" (uvbuf0), "r" (uvbuf1), "r" (dest), "m" (dstW),
-                       "m" (yalpha1), "m" (uvalpha1)
+                               WRITEBGR16(%%esp, 8280(%5), %%eax)
+                               "movl "ESP_OFFSET"(%5), %%esp           \n\t"
+
+                       :: "r" (buf0), "r" (buf1), "r" (uvbuf0), "r" (uvbuf1), "m" (dest),
+                       "r" (&c->redDither)
                        : "%eax"
                        );
+                       return;
+               case IMGFMT_YUY2:
+                       asm volatile(
+                               "movl %%esp, "ESP_OFFSET"(%5)           \n\t"
+                               "movl %4, %%esp                         \n\t"
+                               YSCALEYUV2PACKED1b(%%eax, %5)
+                               WRITEYUY2(%%esp, 8280(%5), %%eax)
+                               "movl "ESP_OFFSET"(%5), %%esp           \n\t"
+
+                       :: "r" (buf0), "r" (buf1), "r" (uvbuf0), "r" (uvbuf1), "m" (dest),
+                       "r" (&c->redDither)
+                       : "%eax"
+                       );
+                       return;
                }
        }
-#else
-//FIXME write 2 versions (for even & odd lines)
-
-       if(dstbpp==32)
-       {
-               int i;
-               for(i=0; i<dstW-1; i+=2){
-                       // vertical linear interpolation && yuv2rgb in a single step:
-                       int Y1=yuvtab_2568[buf0[i]>>7];
-                       int Y2=yuvtab_2568[buf0[i+1]>>7];
-                       int U=((uvbuf0[i>>1]*uvalpha1+uvbuf1[i>>1]*uvalpha)>>19);
-                       int V=((uvbuf0[(i>>1)+2048]*uvalpha1+uvbuf1[(i>>1)+2048]*uvalpha)>>19);
-
-                       int Cb= yuvtab_40cf[U];
-                       int Cg= yuvtab_1a1e[V] + yuvtab_0c92[U];
-                       int Cr= yuvtab_3343[V];
-
-                       dest[4*i+0]=clip_table[((Y1 + Cb) >>13)];
-                       dest[4*i+1]=clip_table[((Y1 + Cg) >>13)];
-                       dest[4*i+2]=clip_table[((Y1 + Cr) >>13)];
-
-                       dest[4*i+4]=clip_table[((Y2 + Cb) >>13)];
-                       dest[4*i+5]=clip_table[((Y2 + Cg) >>13)];
-                       dest[4*i+6]=clip_table[((Y2 + Cr) >>13)];
-               }
-       }
-       else if(dstbpp==24)
-       {
-               int i;
-               for(i=0; i<dstW-1; i+=2){
-                       // vertical linear interpolation && yuv2rgb in a single step:
-                       int Y1=yuvtab_2568[buf0[i]>>7];
-                       int Y2=yuvtab_2568[buf0[i+1]>>7];
-                       int U=((uvbuf0[i>>1]*uvalpha1+uvbuf1[i>>1]*uvalpha)>>19);
-                       int V=((uvbuf0[(i>>1)+2048]*uvalpha1+uvbuf1[(i>>1)+2048]*uvalpha)>>19);
-
-                       int Cb= yuvtab_40cf[U];
-                       int Cg= yuvtab_1a1e[V] + yuvtab_0c92[U];
-                       int Cr= yuvtab_3343[V];
-
-                       dest[0]=clip_table[((Y1 + Cb) >>13)];
-                       dest[1]=clip_table[((Y1 + Cg) >>13)];
-                       dest[2]=clip_table[((Y1 + Cr) >>13)];
-
-                       dest[3]=clip_table[((Y2 + Cb) >>13)];
-                       dest[4]=clip_table[((Y2 + Cg) >>13)];
-                       dest[5]=clip_table[((Y2 + Cr) >>13)];
-                       dest+=6;
-               }
-       }
-       else if(dstbpp==16)
+#endif
+       if( uvalpha < 2048 )
        {
-               int i;
-               for(i=0; i<dstW-1; i+=2){
-                       // vertical linear interpolation && yuv2rgb in a single step:
-                       int Y1=yuvtab_2568[buf0[i]>>7];
-                       int Y2=yuvtab_2568[buf0[i+1]>>7];
-                       int U=((uvbuf0[i>>1]*uvalpha1+uvbuf1[i>>1]*uvalpha)>>19);
-                       int V=((uvbuf0[(i>>1)+2048]*uvalpha1+uvbuf1[(i>>1)+2048]*uvalpha)>>19);
-
-                       int Cb= yuvtab_40cf[U];
-                       int Cg= yuvtab_1a1e[V] + yuvtab_0c92[U];
-                       int Cr= yuvtab_3343[V];
-
-                       ((uint16_t*)dest)[i] =
-                               clip_table16b[(Y1 + Cb) >>13] |
-                               clip_table16g[(Y1 + Cg) >>13] |
-                               clip_table16r[(Y1 + Cr) >>13];
-
-                       ((uint16_t*)dest)[i+1] =
-                               clip_table16b[(Y2 + Cb) >>13] |
-                               clip_table16g[(Y2 + Cg) >>13] |
-                               clip_table16r[(Y2 + Cr) >>13];
-               }
+               YSCALE_YUV_2_ANYRGB_C(YSCALE_YUV_2_RGB1_C, YSCALE_YUV_2_PACKED1_C)
+       }else{
+               YSCALE_YUV_2_ANYRGB_C(YSCALE_YUV_2_RGB1B_C, YSCALE_YUV_2_PACKED1B_C)
        }
-       else if(dstbpp==15)
+}
+
+//FIXME yuy2* can read upto 7 samples to much
+
+static inline void RENAME(yuy2ToY)(uint8_t *dst, uint8_t *src, int width)
+{
+#ifdef HAVE_MMX
+       asm volatile(
+               "movq "MANGLE(bm01010101)", %%mm2\n\t"
+               "movl %0, %%eax                 \n\t"
+               "1:                             \n\t"
+               "movq (%1, %%eax,2), %%mm0      \n\t"
+               "movq 8(%1, %%eax,2), %%mm1     \n\t"
+               "pand %%mm2, %%mm0              \n\t"
+               "pand %%mm2, %%mm1              \n\t"
+               "packuswb %%mm1, %%mm0          \n\t"
+               "movq %%mm0, (%2, %%eax)        \n\t"
+               "addl $8, %%eax                 \n\t"
+               " js 1b                         \n\t"
+               : : "g" (-width), "r" (src+width*2), "r" (dst+width)
+               : "%eax"
+       );
+#else
+       int i;
+       for(i=0; i<width; i++)
+               dst[i]= src[2*i];
+#endif
+}
+
+static inline void RENAME(yuy2ToUV)(uint8_t *dstU, uint8_t *dstV, uint8_t *src1, uint8_t *src2, int width)
+{
+#if defined (HAVE_MMX2) || defined (HAVE_3DNOW)
+       asm volatile(
+               "movq "MANGLE(bm01010101)", %%mm4\n\t"
+               "movl %0, %%eax                 \n\t"
+               "1:                             \n\t"
+               "movq (%1, %%eax,4), %%mm0      \n\t"
+               "movq 8(%1, %%eax,4), %%mm1     \n\t"
+               "movq (%2, %%eax,4), %%mm2      \n\t"
+               "movq 8(%2, %%eax,4), %%mm3     \n\t"
+               PAVGB(%%mm2, %%mm0)
+               PAVGB(%%mm3, %%mm1)
+               "psrlw $8, %%mm0                \n\t"
+               "psrlw $8, %%mm1                \n\t"
+               "packuswb %%mm1, %%mm0          \n\t"
+               "movq %%mm0, %%mm1              \n\t"
+               "psrlw $8, %%mm0                \n\t"
+               "pand %%mm4, %%mm1              \n\t"
+               "packuswb %%mm0, %%mm0          \n\t"
+               "packuswb %%mm1, %%mm1          \n\t"
+               "movd %%mm0, (%4, %%eax)        \n\t"
+               "movd %%mm1, (%3, %%eax)        \n\t"
+               "addl $4, %%eax                 \n\t"
+               " js 1b                         \n\t"
+               : : "g" (-width), "r" (src1+width*4), "r" (src2+width*4), "r" (dstU+width), "r" (dstV+width)
+               : "%eax"
+       );
+#else
+       int i;
+       for(i=0; i<width; i++)
        {
-               int i;
-               for(i=0; i<dstW-1; i+=2){
-                       // vertical linear interpolation && yuv2rgb in a single step:
-                       int Y1=yuvtab_2568[buf0[i]>>7];
-                       int Y2=yuvtab_2568[buf0[i+1]>>7];
-                       int U=((uvbuf0[i>>1]*uvalpha1+uvbuf1[i>>1]*uvalpha)>>19);
-                       int V=((uvbuf0[(i>>1)+2048]*uvalpha1+uvbuf1[(i>>1)+2048]*uvalpha)>>19);
-
-                       int Cb= yuvtab_40cf[U];
-                       int Cg= yuvtab_1a1e[V] + yuvtab_0c92[U];
-                       int Cr= yuvtab_3343[V];
-
-                       ((uint16_t*)dest)[i] =
-                               clip_table15b[(Y1 + Cb) >>13] |
-                               clip_table15g[(Y1 + Cg) >>13] |
-                               clip_table15r[(Y1 + Cr) >>13];
-
-                       ((uint16_t*)dest)[i+1] =
-                               clip_table15b[(Y2 + Cb) >>13] |
-                               clip_table15g[(Y2 + Cg) >>13] |
-                               clip_table15r[(Y2 + Cr) >>13];
-               }
+               dstU[i]= (src1[4*i + 1] + src2[4*i + 1])>>1;
+               dstV[i]= (src1[4*i + 3] + src2[4*i + 3])>>1;
        }
 #endif
 }
 
-// Bilinear / Bicubic scaling
-static inline void RENAME(hScale)(int16_t *dst, int dstW, uint8_t *src, int srcW, int xInc,
-                                 int16_t *filter, int16_t *filterPos, int filterSize)
+//this is allmost identical to the previous, end exists only cuz yuy2ToY/UV)(dst, src+1, ...) would have 100% unaligned accesses
+static inline void RENAME(uyvyToY)(uint8_t *dst, uint8_t *src, int width)
 {
 #ifdef HAVE_MMX
-       if(filterSize==4) // allways true for upscaling, sometimes for down too
-       {
-               int counter= -2*dstW;
+       asm volatile(
+               "movl %0, %%eax                 \n\t"
+               "1:                             \n\t"
+               "movq (%1, %%eax,2), %%mm0      \n\t"
+               "movq 8(%1, %%eax,2), %%mm1     \n\t"
+               "psrlw $8, %%mm0                \n\t"
+               "psrlw $8, %%mm1                \n\t"
+               "packuswb %%mm1, %%mm0          \n\t"
+               "movq %%mm0, (%2, %%eax)        \n\t"
+               "addl $8, %%eax                 \n\t"
+               " js 1b                         \n\t"
+               : : "g" (-width), "r" (src+width*2), "r" (dst+width)
+               : "%eax"
+       );
+#else
+       int i;
+       for(i=0; i<width; i++)
+               dst[i]= src[2*i+1];
+#endif
+}
+
+static inline void RENAME(uyvyToUV)(uint8_t *dstU, uint8_t *dstV, uint8_t *src1, uint8_t *src2, int width)
+{
+#if defined (HAVE_MMX2) || defined (HAVE_3DNOW)
+       asm volatile(
+               "movq "MANGLE(bm01010101)", %%mm4\n\t"
+               "movl %0, %%eax                 \n\t"
+               "1:                             \n\t"
+               "movq (%1, %%eax,4), %%mm0      \n\t"
+               "movq 8(%1, %%eax,4), %%mm1     \n\t"
+               "movq (%2, %%eax,4), %%mm2      \n\t"
+               "movq 8(%2, %%eax,4), %%mm3     \n\t"
+               PAVGB(%%mm2, %%mm0)
+               PAVGB(%%mm3, %%mm1)
+               "pand %%mm4, %%mm0              \n\t"
+               "pand %%mm4, %%mm1              \n\t"
+               "packuswb %%mm1, %%mm0          \n\t"
+               "movq %%mm0, %%mm1              \n\t"
+               "psrlw $8, %%mm0                \n\t"
+               "pand %%mm4, %%mm1              \n\t"
+               "packuswb %%mm0, %%mm0          \n\t"
+               "packuswb %%mm1, %%mm1          \n\t"
+               "movd %%mm0, (%4, %%eax)        \n\t"
+               "movd %%mm1, (%3, %%eax)        \n\t"
+               "addl $4, %%eax                 \n\t"
+               " js 1b                         \n\t"
+               : : "g" (-width), "r" (src1+width*4), "r" (src2+width*4), "r" (dstU+width), "r" (dstV+width)
+               : "%eax"
+       );
+#else
+       int i;
+       for(i=0; i<width; i++)
+       {
+               dstU[i]= (src1[4*i + 0] + src2[4*i + 0])>>1;
+               dstV[i]= (src1[4*i + 2] + src2[4*i + 2])>>1;
+       }
+#endif
+}
+
+static inline void RENAME(bgr32ToY)(uint8_t *dst, uint8_t *src, int width)
+{
+#ifdef HAVE_MMXFIXME
+#else
+       int i;
+       for(i=0; i<width; i++)
+       {
+               int b=  ((uint32_t*)src)[i]&0xFF;
+               int g= (((uint32_t*)src)[i]>>8)&0xFF;
+               int r= (((uint32_t*)src)[i]>>16)&0xFF;
+
+               dst[i]= ((RY*r + GY*g + BY*b + (33<<(RGB2YUV_SHIFT-1)) )>>RGB2YUV_SHIFT);
+       }
+#endif
+}
+
+static inline void RENAME(bgr32ToUV)(uint8_t *dstU, uint8_t *dstV, uint8_t *src1, uint8_t *src2, int width)
+{
+#ifdef HAVE_MMXFIXME
+#else
+       int i;
+       for(i=0; i<width; i++)
+       {
+               const int a= ((uint32_t*)src1)[2*i+0];
+               const int e= ((uint32_t*)src1)[2*i+1];
+               const int c= ((uint32_t*)src2)[2*i+0];
+               const int d= ((uint32_t*)src2)[2*i+1];
+               const int l= (a&0xFF00FF) + (e&0xFF00FF) + (c&0xFF00FF) + (d&0xFF00FF);
+               const int h= (a&0x00FF00) + (e&0x00FF00) + (c&0x00FF00) + (d&0x00FF00);
+               const int b=  l&0x3FF;
+               const int g=  h>>8;
+               const int r=  l>>16;
+
+               dstU[i]= ((RU*r + GU*g + BU*b)>>(RGB2YUV_SHIFT+2)) + 128;
+               dstV[i]= ((RV*r + GV*g + BV*b)>>(RGB2YUV_SHIFT+2)) + 128;
+       }
+#endif
+}
+
+static inline void RENAME(bgr24ToY)(uint8_t *dst, uint8_t *src, int width)
+{
+#ifdef HAVE_MMX
+       asm volatile(
+               "movl %2, %%eax                 \n\t"
+               "movq "MANGLE(bgr2YCoeff)", %%mm6               \n\t"
+               "movq "MANGLE(w1111)", %%mm5            \n\t"
+               "pxor %%mm7, %%mm7              \n\t"
+               "leal (%%eax, %%eax, 2), %%ebx  \n\t"
+               ".balign 16                     \n\t"
+               "1:                             \n\t"
+               PREFETCH" 64(%0, %%ebx)         \n\t"
+               "movd (%0, %%ebx), %%mm0        \n\t"
+               "movd 3(%0, %%ebx), %%mm1       \n\t"
+               "punpcklbw %%mm7, %%mm0         \n\t"
+               "punpcklbw %%mm7, %%mm1         \n\t"
+               "movd 6(%0, %%ebx), %%mm2       \n\t"
+               "movd 9(%0, %%ebx), %%mm3       \n\t"
+               "punpcklbw %%mm7, %%mm2         \n\t"
+               "punpcklbw %%mm7, %%mm3         \n\t"
+               "pmaddwd %%mm6, %%mm0           \n\t"
+               "pmaddwd %%mm6, %%mm1           \n\t"
+               "pmaddwd %%mm6, %%mm2           \n\t"
+               "pmaddwd %%mm6, %%mm3           \n\t"
+#ifndef FAST_BGR2YV12
+               "psrad $8, %%mm0                \n\t"
+               "psrad $8, %%mm1                \n\t"
+               "psrad $8, %%mm2                \n\t"
+               "psrad $8, %%mm3                \n\t"
+#endif
+               "packssdw %%mm1, %%mm0          \n\t"
+               "packssdw %%mm3, %%mm2          \n\t"
+               "pmaddwd %%mm5, %%mm0           \n\t"
+               "pmaddwd %%mm5, %%mm2           \n\t"
+               "packssdw %%mm2, %%mm0          \n\t"
+               "psraw $7, %%mm0                \n\t"
+
+               "movd 12(%0, %%ebx), %%mm4      \n\t"
+               "movd 15(%0, %%ebx), %%mm1      \n\t"
+               "punpcklbw %%mm7, %%mm4         \n\t"
+               "punpcklbw %%mm7, %%mm1         \n\t"
+               "movd 18(%0, %%ebx), %%mm2      \n\t"
+               "movd 21(%0, %%ebx), %%mm3      \n\t"
+               "punpcklbw %%mm7, %%mm2         \n\t"
+               "punpcklbw %%mm7, %%mm3         \n\t"
+               "pmaddwd %%mm6, %%mm4           \n\t"
+               "pmaddwd %%mm6, %%mm1           \n\t"
+               "pmaddwd %%mm6, %%mm2           \n\t"
+               "pmaddwd %%mm6, %%mm3           \n\t"
+#ifndef FAST_BGR2YV12
+               "psrad $8, %%mm4                \n\t"
+               "psrad $8, %%mm1                \n\t"
+               "psrad $8, %%mm2                \n\t"
+               "psrad $8, %%mm3                \n\t"
+#endif
+               "packssdw %%mm1, %%mm4          \n\t"
+               "packssdw %%mm3, %%mm2          \n\t"
+               "pmaddwd %%mm5, %%mm4           \n\t"
+               "pmaddwd %%mm5, %%mm2           \n\t"
+               "addl $24, %%ebx                \n\t"
+               "packssdw %%mm2, %%mm4          \n\t"
+               "psraw $7, %%mm4                \n\t"
+
+               "packuswb %%mm4, %%mm0          \n\t"
+               "paddusb "MANGLE(bgr2YOffset)", %%mm0   \n\t"
+
+               "movq %%mm0, (%1, %%eax)        \n\t"
+               "addl $8, %%eax                 \n\t"
+               " js 1b                         \n\t"
+               : : "r" (src+width*3), "r" (dst+width), "g" (-width)
+               : "%eax", "%ebx"
+       );
+#else
+       int i;
+       for(i=0; i<width; i++)
+       {
+               int b= src[i*3+0];
+               int g= src[i*3+1];
+               int r= src[i*3+2];
+
+               dst[i]= ((RY*r + GY*g + BY*b + (33<<(RGB2YUV_SHIFT-1)) )>>RGB2YUV_SHIFT);
+       }
+#endif
+}
+
+static inline void RENAME(bgr24ToUV)(uint8_t *dstU, uint8_t *dstV, uint8_t *src1, uint8_t *src2, int width)
+{
+#ifdef HAVE_MMX
+       asm volatile(
+               "movl %4, %%eax                 \n\t"
+               "movq "MANGLE(w1111)", %%mm5            \n\t"
+               "movq "MANGLE(bgr2UCoeff)", %%mm6               \n\t"
+               "pxor %%mm7, %%mm7              \n\t"
+               "leal (%%eax, %%eax, 2), %%ebx  \n\t"
+               "addl %%ebx, %%ebx              \n\t"
+               ".balign 16                     \n\t"
+               "1:                             \n\t"
+               PREFETCH" 64(%0, %%ebx)         \n\t"
+               PREFETCH" 64(%1, %%ebx)         \n\t"
+#if defined (HAVE_MMX2) || defined (HAVE_3DNOW)
+               "movq (%0, %%ebx), %%mm0        \n\t"
+               "movq (%1, %%ebx), %%mm1        \n\t"
+               "movq 6(%0, %%ebx), %%mm2       \n\t"
+               "movq 6(%1, %%ebx), %%mm3       \n\t"
+               PAVGB(%%mm1, %%mm0)
+               PAVGB(%%mm3, %%mm2)
+               "movq %%mm0, %%mm1              \n\t"
+               "movq %%mm2, %%mm3              \n\t"
+               "psrlq $24, %%mm0               \n\t"
+               "psrlq $24, %%mm2               \n\t"
+               PAVGB(%%mm1, %%mm0)
+               PAVGB(%%mm3, %%mm2)
+               "punpcklbw %%mm7, %%mm0         \n\t"
+               "punpcklbw %%mm7, %%mm2         \n\t"
+#else
+               "movd (%0, %%ebx), %%mm0        \n\t"
+               "movd (%1, %%ebx), %%mm1        \n\t"
+               "movd 3(%0, %%ebx), %%mm2       \n\t"
+               "movd 3(%1, %%ebx), %%mm3       \n\t"
+               "punpcklbw %%mm7, %%mm0         \n\t"
+               "punpcklbw %%mm7, %%mm1         \n\t"
+               "punpcklbw %%mm7, %%mm2         \n\t"
+               "punpcklbw %%mm7, %%mm3         \n\t"
+               "paddw %%mm1, %%mm0             \n\t"
+               "paddw %%mm3, %%mm2             \n\t"
+               "paddw %%mm2, %%mm0             \n\t"
+               "movd 6(%0, %%ebx), %%mm4       \n\t"
+               "movd 6(%1, %%ebx), %%mm1       \n\t"
+               "movd 9(%0, %%ebx), %%mm2       \n\t"
+               "movd 9(%1, %%ebx), %%mm3       \n\t"
+               "punpcklbw %%mm7, %%mm4         \n\t"
+               "punpcklbw %%mm7, %%mm1         \n\t"
+               "punpcklbw %%mm7, %%mm2         \n\t"
+               "punpcklbw %%mm7, %%mm3         \n\t"
+               "paddw %%mm1, %%mm4             \n\t"
+               "paddw %%mm3, %%mm2             \n\t"
+               "paddw %%mm4, %%mm2             \n\t"
+               "psrlw $2, %%mm0                \n\t"
+               "psrlw $2, %%mm2                \n\t"
+#endif
+               "movq "MANGLE(bgr2VCoeff)", %%mm1               \n\t"
+               "movq "MANGLE(bgr2VCoeff)", %%mm3               \n\t"
+               
+               "pmaddwd %%mm0, %%mm1           \n\t"
+               "pmaddwd %%mm2, %%mm3           \n\t"
+               "pmaddwd %%mm6, %%mm0           \n\t"
+               "pmaddwd %%mm6, %%mm2           \n\t"
+#ifndef FAST_BGR2YV12
+               "psrad $8, %%mm0                \n\t"
+               "psrad $8, %%mm1                \n\t"
+               "psrad $8, %%mm2                \n\t"
+               "psrad $8, %%mm3                \n\t"
+#endif
+               "packssdw %%mm2, %%mm0          \n\t"
+               "packssdw %%mm3, %%mm1          \n\t"
+               "pmaddwd %%mm5, %%mm0           \n\t"
+               "pmaddwd %%mm5, %%mm1           \n\t"
+               "packssdw %%mm1, %%mm0          \n\t" // V1 V0 U1 U0
+               "psraw $7, %%mm0                \n\t"
+
+#if defined (HAVE_MMX2) || defined (HAVE_3DNOW)
+               "movq 12(%0, %%ebx), %%mm4      \n\t"
+               "movq 12(%1, %%ebx), %%mm1      \n\t"
+               "movq 18(%0, %%ebx), %%mm2      \n\t"
+               "movq 18(%1, %%ebx), %%mm3      \n\t"
+               PAVGB(%%mm1, %%mm4)
+               PAVGB(%%mm3, %%mm2)
+               "movq %%mm4, %%mm1              \n\t"
+               "movq %%mm2, %%mm3              \n\t"
+               "psrlq $24, %%mm4               \n\t"
+               "psrlq $24, %%mm2               \n\t"
+               PAVGB(%%mm1, %%mm4)
+               PAVGB(%%mm3, %%mm2)
+               "punpcklbw %%mm7, %%mm4         \n\t"
+               "punpcklbw %%mm7, %%mm2         \n\t"
+#else
+               "movd 12(%0, %%ebx), %%mm4      \n\t"
+               "movd 12(%1, %%ebx), %%mm1      \n\t"
+               "movd 15(%0, %%ebx), %%mm2      \n\t"
+               "movd 15(%1, %%ebx), %%mm3      \n\t"
+               "punpcklbw %%mm7, %%mm4         \n\t"
+               "punpcklbw %%mm7, %%mm1         \n\t"
+               "punpcklbw %%mm7, %%mm2         \n\t"
+               "punpcklbw %%mm7, %%mm3         \n\t"
+               "paddw %%mm1, %%mm4             \n\t"
+               "paddw %%mm3, %%mm2             \n\t"
+               "paddw %%mm2, %%mm4             \n\t"
+               "movd 18(%0, %%ebx), %%mm5      \n\t"
+               "movd 18(%1, %%ebx), %%mm1      \n\t"
+               "movd 21(%0, %%ebx), %%mm2      \n\t"
+               "movd 21(%1, %%ebx), %%mm3      \n\t"
+               "punpcklbw %%mm7, %%mm5         \n\t"
+               "punpcklbw %%mm7, %%mm1         \n\t"
+               "punpcklbw %%mm7, %%mm2         \n\t"
+               "punpcklbw %%mm7, %%mm3         \n\t"
+               "paddw %%mm1, %%mm5             \n\t"
+               "paddw %%mm3, %%mm2             \n\t"
+               "paddw %%mm5, %%mm2             \n\t"
+               "movq "MANGLE(w1111)", %%mm5            \n\t"
+               "psrlw $2, %%mm4                \n\t"
+               "psrlw $2, %%mm2                \n\t"
+#endif
+               "movq "MANGLE(bgr2VCoeff)", %%mm1               \n\t"
+               "movq "MANGLE(bgr2VCoeff)", %%mm3               \n\t"
+               
+               "pmaddwd %%mm4, %%mm1           \n\t"
+               "pmaddwd %%mm2, %%mm3           \n\t"
+               "pmaddwd %%mm6, %%mm4           \n\t"
+               "pmaddwd %%mm6, %%mm2           \n\t"
+#ifndef FAST_BGR2YV12
+               "psrad $8, %%mm4                \n\t"
+               "psrad $8, %%mm1                \n\t"
+               "psrad $8, %%mm2                \n\t"
+               "psrad $8, %%mm3                \n\t"
+#endif
+               "packssdw %%mm2, %%mm4          \n\t"
+               "packssdw %%mm3, %%mm1          \n\t"
+               "pmaddwd %%mm5, %%mm4           \n\t"
+               "pmaddwd %%mm5, %%mm1           \n\t"
+               "addl $24, %%ebx                \n\t"
+               "packssdw %%mm1, %%mm4          \n\t" // V3 V2 U3 U2
+               "psraw $7, %%mm4                \n\t"
+               
+               "movq %%mm0, %%mm1              \n\t"
+               "punpckldq %%mm4, %%mm0         \n\t"
+               "punpckhdq %%mm4, %%mm1         \n\t"
+               "packsswb %%mm1, %%mm0          \n\t"
+               "paddb "MANGLE(bgr2UVOffset)", %%mm0    \n\t"
+
+               "movd %%mm0, (%2, %%eax)        \n\t"
+               "punpckhdq %%mm0, %%mm0         \n\t"
+               "movd %%mm0, (%3, %%eax)        \n\t"
+               "addl $4, %%eax                 \n\t"
+               " js 1b                         \n\t"
+               : : "r" (src1+width*6), "r" (src2+width*6), "r" (dstU+width), "r" (dstV+width), "g" (-width)
+               : "%eax", "%ebx"
+       );
+#else
+       int i;
+       for(i=0; i<width; i++)
+       {
+               int b= src1[6*i + 0] + src1[6*i + 3] + src2[6*i + 0] + src2[6*i + 3];
+               int g= src1[6*i + 1] + src1[6*i + 4] + src2[6*i + 1] + src2[6*i + 4];
+               int r= src1[6*i + 2] + src1[6*i + 5] + src2[6*i + 2] + src2[6*i + 5];
+
+               dstU[i]= ((RU*r + GU*g + BU*b)>>(RGB2YUV_SHIFT+2)) + 128;
+               dstV[i]= ((RV*r + GV*g + BV*b)>>(RGB2YUV_SHIFT+2)) + 128;
+       }
+#endif
+}
+
+static inline void RENAME(bgr16ToY)(uint8_t *dst, uint8_t *src, int width)
+{
+       int i;
+       for(i=0; i<width; i++)
+       {
+               int d= ((uint16_t*)src)[i];
+               int b= d&0x1F;
+               int g= (d>>5)&0x3F;
+               int r= (d>>11)&0x1F;
+
+               dst[i]= ((2*RY*r + GY*g + 2*BY*b)>>(RGB2YUV_SHIFT-2)) + 16;
+       }
+}
+
+static inline void RENAME(bgr16ToUV)(uint8_t *dstU, uint8_t *dstV, uint8_t *src1, uint8_t *src2, int width)
+{
+       int i;
+       for(i=0; i<width; i++)
+       {
+               int d0= ((uint32_t*)src1)[i];
+               int d1= ((uint32_t*)src2)[i];
+               
+               int dl= (d0&0x07E0F81F) + (d1&0x07E0F81F);
+               int dh= ((d0>>5)&0x07C0F83F) + ((d1>>5)&0x07C0F83F);
+
+               int dh2= (dh>>11) + (dh<<21);
+               int d= dh2 + dl;
+
+               int b= d&0x7F;
+               int r= (d>>11)&0x7F;
+               int g= d>>21;
+               dstU[i]= ((2*RU*r + GU*g + 2*BU*b)>>(RGB2YUV_SHIFT+2-2)) + 128;
+               dstV[i]= ((2*RV*r + GV*g + 2*BV*b)>>(RGB2YUV_SHIFT+2-2)) + 128;
+       }
+}
+
+static inline void RENAME(bgr15ToY)(uint8_t *dst, uint8_t *src, int width)
+{
+       int i;
+       for(i=0; i<width; i++)
+       {
+               int d= ((uint16_t*)src)[i];
+               int b= d&0x1F;
+               int g= (d>>5)&0x1F;
+               int r= (d>>10)&0x1F;
+
+               dst[i]= ((RY*r + GY*g + BY*b)>>(RGB2YUV_SHIFT-3)) + 16;
+       }
+}
+
+static inline void RENAME(bgr15ToUV)(uint8_t *dstU, uint8_t *dstV, uint8_t *src1, uint8_t *src2, int width)
+{
+       int i;
+       for(i=0; i<width; i++)
+       {
+               int d0= ((uint32_t*)src1)[i];
+               int d1= ((uint32_t*)src2)[i];
+               
+               int dl= (d0&0x03E07C1F) + (d1&0x03E07C1F);
+               int dh= ((d0>>5)&0x03E0F81F) + ((d1>>5)&0x03E0F81F);
+
+               int dh2= (dh>>11) + (dh<<21);
+               int d= dh2 + dl;
+
+               int b= d&0x7F;
+               int r= (d>>10)&0x7F;
+               int g= d>>21;
+               dstU[i]= ((RU*r + GU*g + BU*b)>>(RGB2YUV_SHIFT+2-3)) + 128;
+               dstV[i]= ((RV*r + GV*g + BV*b)>>(RGB2YUV_SHIFT+2-3)) + 128;
+       }
+}
+
+
+static inline void RENAME(rgb32ToY)(uint8_t *dst, uint8_t *src, int width)
+{
+       int i;
+       for(i=0; i<width; i++)
+       {
+               int r=  ((uint32_t*)src)[i]&0xFF;
+               int g= (((uint32_t*)src)[i]>>8)&0xFF;
+               int b= (((uint32_t*)src)[i]>>16)&0xFF;
+
+               dst[i]= ((RY*r + GY*g + BY*b + (33<<(RGB2YUV_SHIFT-1)) )>>RGB2YUV_SHIFT);
+       }
+}
+
+static inline void RENAME(rgb32ToUV)(uint8_t *dstU, uint8_t *dstV, uint8_t *src1, uint8_t *src2, int width)
+{
+       int i;
+       for(i=0; i<width; i++)
+       {
+               const int a= ((uint32_t*)src1)[2*i+0];
+               const int e= ((uint32_t*)src1)[2*i+1];
+               const int c= ((uint32_t*)src2)[2*i+0];
+               const int d= ((uint32_t*)src2)[2*i+1];
+               const int l= (a&0xFF00FF) + (e&0xFF00FF) + (c&0xFF00FF) + (d&0xFF00FF);
+               const int h= (a&0x00FF00) + (e&0x00FF00) + (c&0x00FF00) + (d&0x00FF00);
+               const int r=  l&0x3FF;
+               const int g=  h>>8;
+               const int b=  l>>16;
+
+               dstU[i]= ((RU*r + GU*g + BU*b)>>(RGB2YUV_SHIFT+2)) + 128;
+               dstV[i]= ((RV*r + GV*g + BV*b)>>(RGB2YUV_SHIFT+2)) + 128;
+       }
+}
+
+static inline void RENAME(rgb24ToY)(uint8_t *dst, uint8_t *src, int width)
+{
+       int i;
+       for(i=0; i<width; i++)
+       {
+               int r= src[i*3+0];
+               int g= src[i*3+1];
+               int b= src[i*3+2];
+
+               dst[i]= ((RY*r + GY*g + BY*b + (33<<(RGB2YUV_SHIFT-1)) )>>RGB2YUV_SHIFT);
+       }
+}
+
+static inline void RENAME(rgb24ToUV)(uint8_t *dstU, uint8_t *dstV, uint8_t *src1, uint8_t *src2, int width)
+{
+       int i;
+       for(i=0; i<width; i++)
+       {
+               int r= src1[6*i + 0] + src1[6*i + 3] + src2[6*i + 0] + src2[6*i + 3];
+               int g= src1[6*i + 1] + src1[6*i + 4] + src2[6*i + 1] + src2[6*i + 4];
+               int b= src1[6*i + 2] + src1[6*i + 5] + src2[6*i + 2] + src2[6*i + 5];
+
+               dstU[i]= ((RU*r + GU*g + BU*b)>>(RGB2YUV_SHIFT+2)) + 128;
+               dstV[i]= ((RV*r + GV*g + BV*b)>>(RGB2YUV_SHIFT+2)) + 128;
+       }
+}
+
+
+// Bilinear / Bicubic scaling
+static inline void RENAME(hScale)(int16_t *dst, int dstW, uint8_t *src, int srcW, int xInc,
+                                 int16_t *filter, int16_t *filterPos, int filterSize)
+{
+#ifdef HAVE_MMX
+       assert(filterSize % 4 == 0 && filterSize>0);
+       if(filterSize==4) // allways true for upscaling, sometimes for down too
+       {
+               int counter= -2*dstW;
                filter-= counter*2;
                filterPos-= counter/2;
                dst-= counter/2;
                asm volatile(
                        "pxor %%mm7, %%mm7              \n\t"
-                       "movq w02, %%mm6                \n\t"
+                       "movq "MANGLE(w02)", %%mm6      \n\t"
                        "pushl %%ebp                    \n\t" // we use 7 regs here ...
                        "movl %%eax, %%ebp              \n\t"
                        ".balign 16                     \n\t"
@@ -1473,7 +2052,7 @@ static inline void RENAME(hScale)(int16_t *dst, int dstW, uint8_t *src, int srcW
                dst-= counter/2;
                asm volatile(
                        "pxor %%mm7, %%mm7              \n\t"
-                       "movq w02, %%mm6                \n\t"
+                       "movq "MANGLE(w02)", %%mm6      \n\t"
                        "pushl %%ebp                    \n\t" // we use 7 regs here ...
                        "movl %%eax, %%ebp              \n\t"
                        ".balign 16                     \n\t"
@@ -1523,7 +2102,7 @@ static inline void RENAME(hScale)(int16_t *dst, int dstW, uint8_t *src, int srcW
                dst-= counter/2;
                asm volatile(
                        "pxor %%mm7, %%mm7              \n\t"
-                       "movq w02, %%mm6                \n\t"
+                       "movq "MANGLE(w02)", %%mm6      \n\t"
                        ".balign 16                     \n\t"
                        "1:                             \n\t"
                        "movl %2, %%ecx                 \n\t"
@@ -1584,13 +2163,58 @@ static inline void RENAME(hScale)(int16_t *dst, int dstW, uint8_t *src, int srcW
 #endif
 }
       // *** horizontal scale Y line to temp buffer
-static inline void RENAME(hyscale)(uint16_t *dst, int dstWidth, uint8_t *src, int srcW, int xInc)
+static inline void RENAME(hyscale)(uint16_t *dst, int dstWidth, uint8_t *src, int srcW, int xInc,
+                                  int flags, int canMMX2BeUsed, int16_t *hLumFilter,
+                                  int16_t *hLumFilterPos, int hLumFilterSize, void *funnyYCode, 
+                                  int srcFormat, uint8_t *formatConvBuffer, int16_t *mmx2Filter,
+                                  int32_t *mmx2FilterPos)
 {
+    if(srcFormat==IMGFMT_YUY2)
+    {
+       RENAME(yuy2ToY)(formatConvBuffer, src, srcW);
+       src= formatConvBuffer;
+    }
+    else if(srcFormat==IMGFMT_UYVY)
+    {
+       RENAME(uyvyToY)(formatConvBuffer, src, srcW);
+       src= formatConvBuffer;
+    }
+    else if(srcFormat==IMGFMT_BGR32)
+    {
+       RENAME(bgr32ToY)(formatConvBuffer, src, srcW);
+       src= formatConvBuffer;
+    }
+    else if(srcFormat==IMGFMT_BGR24)
+    {
+       RENAME(bgr24ToY)(formatConvBuffer, src, srcW);
+       src= formatConvBuffer;
+    }
+    else if(srcFormat==IMGFMT_BGR16)
+    {
+       RENAME(bgr16ToY)(formatConvBuffer, src, srcW);
+       src= formatConvBuffer;
+    }
+    else if(srcFormat==IMGFMT_BGR15)
+    {
+       RENAME(bgr15ToY)(formatConvBuffer, src, srcW);
+       src= formatConvBuffer;
+    }
+    else if(srcFormat==IMGFMT_RGB32)
+    {
+       RENAME(rgb32ToY)(formatConvBuffer, src, srcW);
+       src= formatConvBuffer;
+    }
+    else if(srcFormat==IMGFMT_RGB24)
+    {
+       RENAME(rgb24ToY)(formatConvBuffer, src, srcW);
+       src= formatConvBuffer;
+    }
+
 #ifdef HAVE_MMX
-       // use the new MMX scaler if th mmx2 cant be used (its faster than the x86asm one)
-    if(sws_flags != SWS_FAST_BILINEAR || (!canMMX2BeUsed))
+       // use the new MMX scaler if the mmx2 can't be used (its faster than the x86asm one)
+    if(!(flags&SWS_FAST_BILINEAR) || (!canMMX2BeUsed))
 #else
-    if(sws_flags != SWS_FAST_BILINEAR)
+    if(!(flags&SWS_FAST_BILINEAR))
 #endif
     {
        RENAME(hScale)(dst, dstWidth, src, srcW, xInc, hLumFilter, hLumFilterPos, hLumFilterSize);
@@ -1604,35 +2228,21 @@ static inline void RENAME(hyscale)(uint16_t *dst, int dstWidth, uint8_t *src, in
        {
                asm volatile(
                        "pxor %%mm7, %%mm7              \n\t"
-                       "pxor %%mm2, %%mm2              \n\t" // 2*xalpha
-                       "movd %5, %%mm6                 \n\t" // xInc&0xFFFF
-                       "punpcklwd %%mm6, %%mm6         \n\t"
-                       "punpcklwd %%mm6, %%mm6         \n\t"
-                       "movq %%mm6, %%mm2              \n\t"
-                       "psllq $16, %%mm2               \n\t"
-                       "paddw %%mm6, %%mm2             \n\t"
-                       "psllq $16, %%mm2               \n\t"
-                       "paddw %%mm6, %%mm2             \n\t"
-                       "psllq $16, %%mm2               \n\t" //0,t,2t,3t               t=xInc&0xFF
-                       "movq %%mm2, temp0              \n\t"
-                       "movd %4, %%mm6                 \n\t" //(xInc*4)&0xFFFF
-                       "punpcklwd %%mm6, %%mm6         \n\t"
-                       "punpcklwd %%mm6, %%mm6         \n\t"
+                       "movl %0, %%ecx                 \n\t"
+                       "movl %1, %%edi                 \n\t"
+                       "movl %2, %%edx                 \n\t"
+                       "movl %3, %%ebx                 \n\t"
                        "xorl %%eax, %%eax              \n\t" // i
-                       "movl %0, %%esi                 \n\t" // src
-                       "movl %1, %%edi                 \n\t" // buf1
-                       "movl %3, %%edx                 \n\t" // (xInc*4)>>16
-                       "xorl %%ecx, %%ecx              \n\t"
-                       "xorl %%ebx, %%ebx              \n\t"
-                       "movw %4, %%bx                  \n\t" // (xInc*4)&0xFFFF
+                       PREFETCH" (%%ecx)               \n\t"
+                       PREFETCH" 32(%%ecx)             \n\t"
+                       PREFETCH" 64(%%ecx)             \n\t"
 
 #define FUNNY_Y_CODE \
-                       PREFETCH" 1024(%%esi)           \n\t"\
-                       PREFETCH" 1056(%%esi)           \n\t"\
-                       PREFETCH" 1088(%%esi)           \n\t"\
-                       "call funnyYCode                \n\t"\
-                       "movq temp0, %%mm2              \n\t"\
-                       "xorl %%ecx, %%ecx              \n\t"
+                       "movl (%%ebx), %%esi            \n\t"\
+                       "call *%4                       \n\t"\
+                       "addl (%%ebx, %%eax), %%ecx     \n\t"\
+                       "addl %%eax, %%edi              \n\t"\
+                       "xorl %%eax, %%eax              \n\t"\
 
 FUNNY_Y_CODE
 FUNNY_Y_CODE
@@ -1643,8 +2253,8 @@ FUNNY_Y_CODE
 FUNNY_Y_CODE
 FUNNY_Y_CODE
 
-                       :: "m" (src), "m" (dst), "m" (dstWidth), "m" ((xInc*4)>>16),
-                       "m" ((xInc*4)&0xFFFF), "m" (xInc&0xFFFF)
+                       :: "m" (src), "m" (dst), "m" (mmx2Filter), "m" (mmx2FilterPos),
+                       "m" (funnyYCode)
                        : "%eax", "%ebx", "%ecx", "%edx", "%esi", "%edi"
                );
                for(i=dstWidth-1; (i*xInc)>>16 >=srcW-1; i--) dst[i] = src[srcW-1]*128;
@@ -1693,7 +2303,7 @@ FUNNY_Y_CODE
                : "%eax", "%ebx", "%ecx", "%edi", "%esi"
                );
 #ifdef HAVE_MMX2
-       } //if MMX2 cant be used
+       } //if MMX2 can't be used
 #endif
 #else
        int i;
@@ -1709,14 +2319,70 @@ FUNNY_Y_CODE
     }
 }
 
-inline static void RENAME(hcscale)(uint16_t *dst, int dstWidth,
-                               uint8_t *src1, uint8_t *src2, int srcW, int xInc)
+inline static void RENAME(hcscale)(uint16_t *dst, int dstWidth, uint8_t *src1, uint8_t *src2,
+                                  int srcW, int xInc, int flags, int canMMX2BeUsed, int16_t *hChrFilter,
+                                  int16_t *hChrFilterPos, int hChrFilterSize, void *funnyUVCode,
+                                  int srcFormat, uint8_t *formatConvBuffer, int16_t *mmx2Filter,
+                                  int32_t *mmx2FilterPos)
 {
+    if(srcFormat==IMGFMT_YUY2)
+    {
+       RENAME(yuy2ToUV)(formatConvBuffer, formatConvBuffer+2048, src1, src2, srcW);
+       src1= formatConvBuffer;
+       src2= formatConvBuffer+2048;
+    }
+    else if(srcFormat==IMGFMT_UYVY)
+    {
+       RENAME(uyvyToUV)(formatConvBuffer, formatConvBuffer+2048, src1, src2, srcW);
+       src1= formatConvBuffer;
+       src2= formatConvBuffer+2048;
+    }
+    else if(srcFormat==IMGFMT_BGR32)
+    {
+       RENAME(bgr32ToUV)(formatConvBuffer, formatConvBuffer+2048, src1, src2, srcW);
+       src1= formatConvBuffer;
+       src2= formatConvBuffer+2048;
+    }
+    else if(srcFormat==IMGFMT_BGR24)
+    {
+       RENAME(bgr24ToUV)(formatConvBuffer, formatConvBuffer+2048, src1, src2, srcW);
+       src1= formatConvBuffer;
+       src2= formatConvBuffer+2048;
+    }
+    else if(srcFormat==IMGFMT_BGR16)
+    {
+       RENAME(bgr16ToUV)(formatConvBuffer, formatConvBuffer+2048, src1, src2, srcW);
+       src1= formatConvBuffer;
+       src2= formatConvBuffer+2048;
+    }
+    else if(srcFormat==IMGFMT_BGR15)
+    {
+       RENAME(bgr15ToUV)(formatConvBuffer, formatConvBuffer+2048, src1, src2, srcW);
+       src1= formatConvBuffer;
+       src2= formatConvBuffer+2048;
+    }
+    else if(srcFormat==IMGFMT_RGB32)
+    {
+       RENAME(rgb32ToUV)(formatConvBuffer, formatConvBuffer+2048, src1, src2, srcW);
+       src1= formatConvBuffer;
+       src2= formatConvBuffer+2048;
+    }
+    else if(srcFormat==IMGFMT_RGB24)
+    {
+       RENAME(rgb24ToUV)(formatConvBuffer, formatConvBuffer+2048, src1, src2, srcW);
+       src1= formatConvBuffer;
+       src2= formatConvBuffer+2048;
+    }
+    else if(isGray(srcFormat))
+    {
+       return;
+    }
+
 #ifdef HAVE_MMX
-       // use the new MMX scaler if th mmx2 cant be used (its faster than the x86asm one)
-    if(sws_flags != SWS_FAST_BILINEAR || (!canMMX2BeUsed))
+       // use the new MMX scaler if the mmx2 can't be used (its faster than the x86asm one)
+    if(!(flags&SWS_FAST_BILINEAR) || (!canMMX2BeUsed))
 #else
-    if(sws_flags != SWS_FAST_BILINEAR)
+    if(!(flags&SWS_FAST_BILINEAR))
 #endif
     {
        RENAME(hScale)(dst     , dstWidth, src1, srcW, xInc, hChrFilter, hChrFilterPos, hChrFilterSize);
@@ -1730,65 +2396,44 @@ inline static void RENAME(hcscale)(uint16_t *dst, int dstWidth,
        if(canMMX2BeUsed)
        {
                asm volatile(
-               "pxor %%mm7, %%mm7              \n\t"
-               "pxor %%mm2, %%mm2              \n\t" // 2*xalpha
-               "movd %5, %%mm6                 \n\t" // xInc&0xFFFF
-               "punpcklwd %%mm6, %%mm6         \n\t"
-               "punpcklwd %%mm6, %%mm6         \n\t"
-               "movq %%mm6, %%mm2              \n\t"
-               "psllq $16, %%mm2               \n\t"
-               "paddw %%mm6, %%mm2             \n\t"
-               "psllq $16, %%mm2               \n\t"
-               "paddw %%mm6, %%mm2             \n\t"
-               "psllq $16, %%mm2               \n\t" //0,t,2t,3t               t=xInc&0xFFFF
-               "movq %%mm2, temp0              \n\t"
-               "movd %4, %%mm6                 \n\t" //(xInc*4)&0xFFFF
-               "punpcklwd %%mm6, %%mm6         \n\t"
-               "punpcklwd %%mm6, %%mm6         \n\t"
-               "xorl %%eax, %%eax              \n\t" // i
-               "movl %0, %%esi                 \n\t" // src
-               "movl %1, %%edi                 \n\t" // buf1
-               "movl %3, %%edx                 \n\t" // (xInc*4)>>16
-               "xorl %%ecx, %%ecx              \n\t"
-               "xorl %%ebx, %%ebx              \n\t"
-               "movw %4, %%bx                  \n\t" // (xInc*4)&0xFFFF
-
-#define FUNNYUVCODE \
-                       PREFETCH" 1024(%%esi)           \n\t"\
-                       PREFETCH" 1056(%%esi)           \n\t"\
-                       PREFETCH" 1088(%%esi)           \n\t"\
-                       "call funnyUVCode               \n\t"\
-                       "movq temp0, %%mm2              \n\t"\
-                       "xorl %%ecx, %%ecx              \n\t"
-
-FUNNYUVCODE
-FUNNYUVCODE
-FUNNYUVCODE
-FUNNYUVCODE
-
-FUNNYUVCODE
-FUNNYUVCODE
-FUNNYUVCODE
-FUNNYUVCODE
-               "xorl %%eax, %%eax              \n\t" // i
-               "movl %6, %%esi                 \n\t" // src
-               "movl %1, %%edi                 \n\t" // buf1
-               "addl $4096, %%edi              \n\t"
-
-FUNNYUVCODE
-FUNNYUVCODE
-FUNNYUVCODE
-FUNNYUVCODE
-
-FUNNYUVCODE
-FUNNYUVCODE
-FUNNYUVCODE
-FUNNYUVCODE
-
-               :: "m" (src1), "m" (dst), "m" (dstWidth), "m" ((xInc*4)>>16),
-                 "m" ((xInc*4)&0xFFFF), "m" (xInc&0xFFFF), "m" (src2)
-               : "%eax", "%ebx", "%ecx", "%edx", "%esi", "%edi"
-       );
+                       "pxor %%mm7, %%mm7              \n\t"
+                       "movl %0, %%ecx                 \n\t"
+                       "movl %1, %%edi                 \n\t"
+                       "movl %2, %%edx                 \n\t"
+                       "movl %3, %%ebx                 \n\t"
+                       "xorl %%eax, %%eax              \n\t" // i
+                       PREFETCH" (%%ecx)               \n\t"
+                       PREFETCH" 32(%%ecx)             \n\t"
+                       PREFETCH" 64(%%ecx)             \n\t"
+
+#define FUNNY_UV_CODE \
+                       "movl (%%ebx), %%esi            \n\t"\
+                       "call *%4                       \n\t"\
+                       "addl (%%ebx, %%eax), %%ecx     \n\t"\
+                       "addl %%eax, %%edi              \n\t"\
+                       "xorl %%eax, %%eax              \n\t"\
+
+FUNNY_UV_CODE
+FUNNY_UV_CODE
+FUNNY_UV_CODE
+FUNNY_UV_CODE
+                       "xorl %%eax, %%eax              \n\t" // i
+                       "movl %5, %%ecx                 \n\t" // src
+                       "movl %1, %%edi                 \n\t" // buf1
+                       "addl $4096, %%edi              \n\t"
+                       PREFETCH" (%%ecx)               \n\t"
+                       PREFETCH" 32(%%ecx)             \n\t"
+                       PREFETCH" 64(%%ecx)             \n\t"
+
+FUNNY_UV_CODE
+FUNNY_UV_CODE
+FUNNY_UV_CODE
+FUNNY_UV_CODE
+
+                       :: "m" (src1), "m" (dst), "m" (mmx2Filter), "m" (mmx2FilterPos),
+                       "m" (funnyUVCode), "m" (src2)
+                       : "%eax", "%ebx", "%ecx", "%edx", "%esi", "%edi"
+               );
                for(i=dstWidth-1; (i*xInc)>>16 >=srcW-1; i--)
                {
 //                     printf("%d %d %d\n", dstWidth, i, srcW);
@@ -1837,7 +2482,7 @@ FUNNYUVCODE
                : "%eax", "%ebx", "%ecx", "%edi", "%esi"
                );
 #ifdef HAVE_MMX2
-       } //if MMX2 cant be used
+       } //if MMX2 can't be used
 #endif
 #else
        int i;
@@ -1858,618 +2503,157 @@ FUNNYUVCODE
    }
 }
 
-static inline void RENAME(initFilter)(int16_t *dstFilter, int16_t *filterPos, int *filterSize, int xInc,
-                                     int srcW, int dstW, int filterAlign, int one)
-{
-       int i;
-       double filter[8000];
-#ifdef HAVE_MMX
-       asm volatile("emms\n\t"::: "memory"); //FIXME this shouldnt be required but it IS (even for non mmx versions)
-#endif
-
-       if(ABS(xInc - 0x10000) <10) // unscaled
-       {
-               int i;
-               *filterSize= (1 +(filterAlign-1)) & (~(filterAlign-1)); // 1 or 4 normaly
-               for(i=0; i<dstW*(*filterSize); i++) filter[i]=0;
-
-               for(i=0; i<dstW; i++)
-               {
-                       filter[i*(*filterSize)]=1;
-                       filterPos[i]=i;
-               }
-
+static int RENAME(swScale)(SwsContext *c, uint8_t* src[], int srcStride[], int srcSliceY,
+             int srcSliceH, uint8_t* dst[], int dstStride[]){
+
+       /* load a few things into local vars to make the code more readable? and faster */
+       const int srcW= c->srcW;
+       const int dstW= c->dstW;
+       const int dstH= c->dstH;
+       const int chrDstW= c->chrDstW;
+       const int chrSrcW= c->chrSrcW;
+       const int lumXInc= c->lumXInc;
+       const int chrXInc= c->chrXInc;
+       const int dstFormat= c->dstFormat;
+       const int srcFormat= c->srcFormat;
+       const int flags= c->flags;
+       const int canMMX2BeUsed= c->canMMX2BeUsed;
+       int16_t *vLumFilterPos= c->vLumFilterPos;
+       int16_t *vChrFilterPos= c->vChrFilterPos;
+       int16_t *hLumFilterPos= c->hLumFilterPos;
+       int16_t *hChrFilterPos= c->hChrFilterPos;
+       int16_t *vLumFilter= c->vLumFilter;
+       int16_t *vChrFilter= c->vChrFilter;
+       int16_t *hLumFilter= c->hLumFilter;
+       int16_t *hChrFilter= c->hChrFilter;
+       int32_t *lumMmxFilter= c->lumMmxFilter;
+       int32_t *chrMmxFilter= c->chrMmxFilter;
+       const int vLumFilterSize= c->vLumFilterSize;
+       const int vChrFilterSize= c->vChrFilterSize;
+       const int hLumFilterSize= c->hLumFilterSize;
+       const int hChrFilterSize= c->hChrFilterSize;
+       int16_t **lumPixBuf= c->lumPixBuf;
+       int16_t **chrPixBuf= c->chrPixBuf;
+       const int vLumBufSize= c->vLumBufSize;
+       const int vChrBufSize= c->vChrBufSize;
+       uint8_t *funnyYCode= c->funnyYCode;
+       uint8_t *funnyUVCode= c->funnyUVCode;
+       uint8_t *formatConvBuffer= c->formatConvBuffer;
+       const int chrSrcSliceY= srcSliceY >> c->chrSrcVSubSample;
+       const int chrSrcSliceH= -((-srcSliceH) >> c->chrSrcVSubSample);
+       int lastDstY;
+
+       /* vars whch will change and which we need to storw back in the context */
+       int dstY= c->dstY;
+       int lumBufIndex= c->lumBufIndex;
+       int chrBufIndex= c->chrBufIndex;
+       int lastInLumBuf= c->lastInLumBuf;
+       int lastInChrBuf= c->lastInChrBuf;
+       
+       if(isPacked(c->srcFormat)){
+               src[0]=
+               src[1]=
+               src[2]= src[0];
+               srcStride[0]=
+               srcStride[1]=
+               srcStride[2]= srcStride[0];
        }
-       else if(xInc <= (1<<16) || sws_flags==SWS_FAST_BILINEAR) // upscale
-       {
-               int i;
-               int xDstInSrc;
-               if(sws_flags==SWS_BICUBIC) *filterSize= 4;
-               else                       *filterSize= 2;
-//             printf("%d %d %d\n", filterSize, srcW, dstW);
-               *filterSize= (*filterSize +(filterAlign-1)) & (~(filterAlign-1));
-
-               xDstInSrc= xInc - 0x8000;
-               for(i=0; i<dstW; i++)
-               {
-                       int xx= (xDstInSrc>>16) - (*filterSize>>1) + 1;
-                       int j;
-
-                       filterPos[i]= xx;
-                       if(sws_flags == SWS_BICUBIC)
-                       {
-                               double d= ABS(((xx+1)<<16) - xDstInSrc)/(double)(1<<16);
-                               double y1,y2,y3,y4;
-                               double A= -0.75;
-                                       // Equation is from VirtualDub
-                               y1 = (        +     A*d -       2.0*A*d*d +       A*d*d*d);
-                               y2 = (+ 1.0             -     (A+3.0)*d*d + (A+2.0)*d*d*d);
-                               y3 = (        -     A*d + (2.0*A+3.0)*d*d - (A+2.0)*d*d*d);
-                               y4 = (                  +           A*d*d -       A*d*d*d);
-
-//                             printf("%d %d %d \n", coeff, (int)d, xDstInSrc);
-                               filter[i*(*filterSize) + 0]= y1;
-                               filter[i*(*filterSize) + 1]= y2;
-                               filter[i*(*filterSize) + 2]= y3;
-                               filter[i*(*filterSize) + 3]= y4;
-//                             printf("%1.3f %d, %d, %d, %d\n",d , y1, y2, y3, y4);
-                       }
-                       else
-                       {
-                               for(j=0; j<*filterSize; j++)
-                               {
-                                       double d= ABS((xx<<16) - xDstInSrc)/(double)(1<<16);
-                                       double coeff= 1.0 - d;
-                                       if(coeff<0) coeff=0;
-       //                              printf("%d %d %d \n", coeff, (int)d, xDstInSrc);
-                                       filter[i*(*filterSize) + j]= coeff;
-                                       xx++;
-                               }
-                       }
-                       xDstInSrc+= xInc;
-               }
-       }
-       else // downscale
-       {
-               int xDstInSrc;
-               if(sws_flags==SWS_BICUBIC) *filterSize= (int)ceil(1 + 4.0*srcW / (double)dstW);
-               else                       *filterSize= (int)ceil(1 + 2.0*srcW / (double)dstW);
-//             printf("%d %d %d\n", *filterSize, srcW, dstW);
-               *filterSize= (*filterSize +(filterAlign-1)) & (~(filterAlign-1));
-
-               xDstInSrc= xInc - 0x8000;
-               for(i=0; i<dstW; i++)
-               {
-                       int xx= (int)((double)xDstInSrc/(double)(1<<16) - *filterSize*0.5 + 0.5);
-                       int j;
+       srcStride[1]<<= c->vChrDrop;
+       srcStride[2]<<= c->vChrDrop;
 
-                       filterPos[i]= xx;
-                       for(j=0; j<*filterSize; j++)
-                       {
-                               double d= ABS((xx<<16) - xDstInSrc)/(double)xInc;
-                               double coeff;
-                               if(sws_flags == SWS_BICUBIC)
-                               {
-                                       double A= -0.75;
-//                                     d*=2;
-                                       // Equation is from VirtualDub
-                                       if(d<1.0)
-                                               coeff = (1.0 - (A+3.0)*d*d + (A+2.0)*d*d*d);
-                                       else if(d<2.0)
-                                               coeff = (-4.0*A + 8.0*A*d - 5.0*A*d*d + A*d*d*d);
-                                       else
-                                               coeff=0.0;
-                               }
-                               else
-                               {
-                                       coeff= 1.0 - d;
-                                       if(coeff<0) coeff=0;
-                               }
-//                             if(filterAlign==1) printf("%d %d %d \n", coeff, (int)d, xDstInSrc);
-                               filter[i*(*filterSize) + j]= coeff;
-                               xx++;
-                       }
-                       xDstInSrc+= xInc;
-               }
-       }
-
-       //fix borders
-       for(i=0; i<dstW; i++)
-       {
-               int j;
-               if(filterPos[i] < 0)
-               {
-                       // Move filter coeffs left to compensate for filterPos
-                       for(j=1; j<*filterSize; j++)
-                       {
-                               int left= MAX(j + filterPos[i], 0);
-                               filter[i*(*filterSize) + left] += filter[i*(*filterSize) + j];
-                               filter[i*(*filterSize) + j]=0;
-                       }
-                       filterPos[i]= 0;
-               }
-
-               if(filterPos[i] + (*filterSize) > srcW)
-               {
-                       int shift= filterPos[i] + (*filterSize) - srcW;
-                       // Move filter coeffs right to compensate for filterPos
-                       for(j=(*filterSize)-2; j>=0; j--)
-                       {
-                               int right= MIN(j + shift, (*filterSize)-1);
-                               filter[i*(*filterSize) +right] += filter[i*(*filterSize) +j];
-                               filter[i*(*filterSize) +j]=0;
-                       }
-                       filterPos[i]= srcW - (*filterSize);
-               }
-       }
-
-       //FIXME try to align filterpos if possible / try to shift filterpos to put zeros at the end
-       // and skip these than later
-
-       //Normalize
-       for(i=0; i<dstW; i++)
-       {
-               int j;
-               double sum=0;
-               double scale= one;
-               for(j=0; j<*filterSize; j++)
-               {
-                       sum+= filter[i*(*filterSize) + j];
-               }
-               scale/= sum;
-               for(j=0; j<*filterSize; j++)
-               {
-                       dstFilter[i*(*filterSize) + j]= (int)(filter[i*(*filterSize) + j]*scale);
-               }
-       }
-}
+//     printf("swscale %X %X %X -> %X %X %X\n", (int)src[0], (int)src[1], (int)src[2],
+//             (int)dst[0], (int)dst[1], (int)dst[2]);
 
-#ifdef HAVE_MMX2
-static void initMMX2HScaler(int dstW, int xInc, uint8_t *funnyCode)
+#if 0 //self test FIXME move to a vfilter or something
 {
-       uint8_t *fragment;
-       int imm8OfPShufW1;
-       int imm8OfPShufW2;
-       int fragmentLength;
-
-       int xpos, i;
-
-       // create an optimized horizontal scaling routine
-
-       //code fragment
-
-       asm volatile(
-               "jmp 9f                         \n\t"
-       // Begin
-               "0:                             \n\t"
-               "movq (%%esi), %%mm0            \n\t" //FIXME Alignment
-               "movq %%mm0, %%mm1              \n\t"
-               "psrlq $8, %%mm0                \n\t"
-               "punpcklbw %%mm7, %%mm1 \n\t"
-               "movq %%mm2, %%mm3              \n\t"
-               "punpcklbw %%mm7, %%mm0 \n\t"
-               "addw %%bx, %%cx                \n\t" //2*xalpha += (4*lumXInc)&0xFFFF
-               "pshufw $0xFF, %%mm1, %%mm1     \n\t"
-               "1:                             \n\t"
-               "adcl %%edx, %%esi              \n\t" //xx+= (4*lumXInc)>>16 + carry
-               "pshufw $0xFF, %%mm0, %%mm0     \n\t"
-               "2:                             \n\t"
-               "psrlw $9, %%mm3                \n\t"
-               "psubw %%mm1, %%mm0             \n\t"
-               "pmullw %%mm3, %%mm0            \n\t"
-               "paddw %%mm6, %%mm2             \n\t" // 2*alpha += xpos&0xFFFF
-               "psllw $7, %%mm1                \n\t"
-               "paddw %%mm1, %%mm0             \n\t"
-
-               "movq %%mm0, (%%edi, %%eax)     \n\t"
-
-               "addl $8, %%eax                 \n\t"
-       // End
-               "9:                             \n\t"
-//             "int $3\n\t"
-               "leal 0b, %0                    \n\t"
-               "leal 1b, %1                    \n\t"
-               "leal 2b, %2                    \n\t"
-               "decl %1                        \n\t"
-               "decl %2                        \n\t"
-               "subl %0, %1                    \n\t"
-               "subl %0, %2                    \n\t"
-               "leal 9b, %3                    \n\t"
-               "subl %0, %3                    \n\t"
-               :"=r" (fragment), "=r" (imm8OfPShufW1), "=r" (imm8OfPShufW2),
-               "=r" (fragmentLength)
-       );
-
-       xpos= 0; //lumXInc/2 - 0x8000; // difference between pixel centers
-
-       for(i=0; i<dstW/8; i++)
-       {
-               int xx=xpos>>16;
-
-               if((i&3) == 0)
-               {
-                       int a=0;
-                       int b=((xpos+xInc)>>16) - xx;
-                       int c=((xpos+xInc*2)>>16) - xx;
-                       int d=((xpos+xInc*3)>>16) - xx;
-
-                       memcpy(funnyCode + fragmentLength*i/4, fragment, fragmentLength);
-
-                       funnyCode[fragmentLength*i/4 + imm8OfPShufW1]=
-                       funnyCode[fragmentLength*i/4 + imm8OfPShufW2]=
-                               a | (b<<2) | (c<<4) | (d<<6);
-
-                       // if we dont need to read 8 bytes than dont :), reduces the chance of
-                       // crossing a cache line
-                       if(d<3) funnyCode[fragmentLength*i/4 + 1]= 0x6E;
-
-                       funnyCode[fragmentLength*(i+4)/4]= RET;
-               }
-               xpos+=xInc;
-       }
-/*
-       xpos= 0; //chrXInc/2 - 0x10000; // difference between centers of chrom samples
-       for(i=0; i<dstUVw/8; i++)
-       {
-               int xx=xpos>>16;
-
-               if((i&3) == 0)
-               {
-                       int a=0;
-                       int b=((xpos+chrXInc)>>16) - xx;
-                       int c=((xpos+chrXInc*2)>>16) - xx;
-                       int d=((xpos+chrXInc*3)>>16) - xx;
-
-                       memcpy(funnyUVCode + fragmentLength*i/4, fragment, fragmentLength);
-
-                       funnyUVCode[fragmentLength*i/4 + imm8OfPShufW1]=
-                       funnyUVCode[fragmentLength*i/4 + imm8OfPShufW2]=
-                               a | (b<<2) | (c<<4) | (d<<6);
-
-                       // if we dont need to read 8 bytes than dont :), reduces the chance of
-                       // crossing a cache line
-                       if(d<3) funnyUVCode[fragmentLength*i/4 + 1]= 0x6E;
-
-                       funnyUVCode[fragmentLength*(i+4)/4]= RET;
-               }
-               xpos+=chrXInc;
-       }
-*/
-//     funnyCode[0]= RET;
+static volatile int i=0;
+i++;
+if(srcFormat==IMGFMT_YV12 && i==1 && srcSliceH>= c->srcH)
+       selfTest(src, srcStride, c->srcW, c->srcH);
+i--;
 }
-#endif // HAVE_MMX2
-
-static void RENAME(SwScale_YV12slice)(unsigned char* srcptr[],int stride[], int srcSliceY ,
-                            int srcSliceH, uint8_t* dstptr[], int dststride, int dstbpp,
-                            int srcW, int srcH, int dstW, int dstH){
-
-
-unsigned int lumXInc= (srcW << 16) / dstW;
-unsigned int lumYInc= (srcH << 16) / dstH;
-unsigned int chrXInc;
-unsigned int chrYInc;
-
-static int dstY;
-
-// used to detect a size change
-static int oldDstW= -1;
-static int oldSrcW= -1;
-static int oldDstH= -1;
-static int oldSrcH= -1;
-static int oldFlags=-1;
-
-static int lastInLumBuf;
-static int lastInChrBuf;
-
-int chrDstW, chrDstH;
-
-static int lumBufIndex=0;
-static int chrBufIndex=0;
-
-static int firstTime=1;
-
-const int widthAlign= dstbpp==12 ? 16 : 8;
-const int bytespp= (dstbpp+1)/8; //(12->1, 15&16->2, 24->3, 32->4)
-const int over= dstbpp==12 ?     (((dstW+15)&(~15))) - dststride
-                               : (((dstW+7)&(~7)))*bytespp - dststride;
-if(dststride%widthAlign !=0 )
-{
-       if(firstTime)
-               fprintf(stderr, "SwScaler: Warning: dstStride is not a multiple of %d!\n"
-                               "SwScaler:          ->cannot do aligned memory acesses anymore\n",
-                               widthAlign);
-}
-
-if(over>0)
-{
-       if(firstTime)
-               fprintf(stderr, "SwScaler: Warning: output width is not a multiple of 8 (16 for YV12)\n"
-                               "SwScaler:          and dststride is not large enough to handle %d extra bytes\n"
-                               "SwScaler:          ->using unoptimized C version for last line(s)\n",
-                               over);
-}
-
-
-
-//printf("%d %d %d %d\n", srcW, srcH, dstW, dstH);
-//printf("%d %d %d %d\n", lumXInc, lumYInc, srcSliceY, srcSliceH);
-
-#ifdef HAVE_MMX2
-canMMX2BeUsed= (lumXInc <= 0x10000 && (dstW&31)==0 && (srcW&15)==0) ? 1 : 0;
-if(!canMMX2BeUsed && lumXInc <= 0x10000 && (srcW&15)==0 && sws_flags==SWS_FAST_BILINEAR)
-{
-       if(firstTime)
-               fprintf(stderr, "SwScaler: output Width is not a multiple of 32 -> no MMX2 scaler\n");
-}
-#else
-canMMX2BeUsed=0; // should be 0 anyway but ...
 #endif
 
-if(firstTime)
-{
-#if defined (DITHER1XBPP) && defined (HAVE_MMX)
-       char *dither= " dithered";
-#else
-       char *dither= "";
-#endif
-       if(sws_flags==SWS_FAST_BILINEAR)
-               fprintf(stderr, "SwScaler: FAST_BILINEAR scaler ");
-       else if(sws_flags==SWS_BILINEAR)
-               fprintf(stderr, "SwScaler: BILINEAR scaler ");
-       else if(sws_flags==SWS_BICUBIC)
-               fprintf(stderr, "SwScaler: BICUBIC scaler ");
-       else
-               fprintf(stderr, "SwScaler: ehh flags invalid?! ");
-
-       if(dstbpp==15)
-               fprintf(stderr, "with%s BGR15 output ", dither);
-       else if(dstbpp==16)
-               fprintf(stderr, "with%s BGR16 output ", dither);
-       else if(dstbpp==24)
-               fprintf(stderr, "with BGR24 output ");
-       else if(dstbpp==32)
-               fprintf(stderr, "with BGR32 output ");
-       else if(dstbpp==12)
-               fprintf(stderr, "with YV12 output ");
-       else
-               fprintf(stderr, "without output ");
+//printf("sws Strides:%d %d %d -> %d %d %d\n", srcStride[0],srcStride[1],srcStride[2],
+//dstStride[0],dstStride[1],dstStride[2]);
 
-#ifdef HAVE_MMX2
-               fprintf(stderr, "using MMX2\n");
-#elif defined (HAVE_3DNOW)
-               fprintf(stderr, "using 3DNOW\n");
-#elif defined (HAVE_MMX)
-               fprintf(stderr, "using MMX\n");
-#elif defined (ARCH_X86)
-               fprintf(stderr, "using X86 ASM\n");
-#else
-               fprintf(stderr, "using C\n");
-#endif
-}
-
-
-// match pixel 0 of the src to pixel 0 of dst and match pixel n-2 of src to pixel n-2 of dst
-// n-2 is the last chrominance sample available
-// this is not perfect, but noone shuld notice the difference, the more correct variant
-// would be like the vertical one, but that would require some special code for the
-// first and last pixel
-if(sws_flags==SWS_FAST_BILINEAR)
-{
-       if(canMMX2BeUsed)       lumXInc+= 20;
-#ifndef HAVE_MMX //we dont use the x86asm scaler if mmx is available
-       else                    lumXInc = ((srcW-2)<<16)/(dstW-2) - 20;
-#endif
-}
-
-if(fullUVIpol && !(dstbpp==12))        chrXInc= lumXInc>>1, chrDstW= dstW;
-else                                   chrXInc= lumXInc,    chrDstW= (dstW+1)>>1;
-
-if(dstbpp==12) chrYInc= lumYInc,    chrDstH= (dstH+1)>>1;
-else           chrYInc= lumYInc>>1, chrDstH= dstH;
-
-  // force calculation of the horizontal interpolation of the first line
-
-  if(srcSliceY ==0){
-//     printf("dstW %d, srcw %d, mmx2 %d\n", dstW, srcW, canMMX2BeUsed);
-       lumBufIndex=0;
-       chrBufIndex=0;
-       dstY=0;
-
-       //precalculate horizontal scaler filter coefficients
-       if(oldDstW!=dstW || oldSrcW!=srcW || oldFlags!=sws_flags)
+       if(dstStride[0]%8 !=0 || dstStride[1]%8 !=0 || dstStride[2]%8 !=0)
        {
-#ifdef HAVE_MMX
-               const int filterAlign=4;
-#else
-               const int filterAlign=1;
-#endif
-               oldDstW= dstW; oldSrcW= srcW; oldFlags= sws_flags;
-
-               RENAME(initFilter)(hLumFilter, hLumFilterPos, &hLumFilterSize, lumXInc,
-                               srcW   , dstW   , filterAlign, 1<<14);
-               RENAME(initFilter)(hChrFilter, hChrFilterPos, &hChrFilterSize, chrXInc,
-                               (srcW+1)>>1, chrDstW, filterAlign, 1<<14);
-
-#ifdef HAVE_MMX2
-// cant downscale !!!
-               if(canMMX2BeUsed && sws_flags == SWS_FAST_BILINEAR)
+               static int firstTime=1; //FIXME move this into the context perhaps
+               if(flags & SWS_PRINT_INFO && firstTime)
                {
-                       initMMX2HScaler(dstW   , lumXInc, funnyYCode);
-                       initMMX2HScaler(chrDstW, chrXInc, funnyUVCode);
+                       MSG_WARN("SwScaler: Warning: dstStride is not aligned!\n"
+                                       "SwScaler:          ->cannot do aligned memory acesses anymore\n");
+                       firstTime=0;
                }
-#endif
-       } // Init Horizontal stuff
-
-       if(oldDstH!=dstH || oldSrcH!=srcH || oldFlags!=sws_flags)
-       {
-               int i;
-               oldDstH= dstH; oldSrcH= srcH; oldFlags= sws_flags; //FIXME swsflags conflict with x check
-
-               // deallocate pixbufs
-               for(i=0; i<vLumBufSize; i++) free(lumPixBuf[i]);
-               for(i=0; i<vChrBufSize; i++) free(chrPixBuf[i]);
-
-               RENAME(initFilter)(vLumFilter, vLumFilterPos, &vLumFilterSize, lumYInc,
-                               srcH   , dstH,    1, (1<<12)-4);
-               RENAME(initFilter)(vChrFilter, vChrFilterPos, &vChrFilterSize, chrYInc,
-                               (srcH+1)>>1, chrDstH, 1, (1<<12)-4);
-
-               // Calculate Buffer Sizes so that they wont run out while handling these damn slices
-               vLumBufSize= vLumFilterSize; vChrBufSize= vChrFilterSize;
-               for(i=0; i<dstH; i++)
-               {
-                       int chrI= i*chrDstH / dstH;
-                       int nextSlice= MAX(vLumFilterPos[i   ] + vLumFilterSize - 1,
-                                        ((vChrFilterPos[chrI] + vChrFilterSize - 1)<<1));
-                       nextSlice&= ~1; // Slices start at even boundaries
-                       if(vLumFilterPos[i   ] + vLumBufSize < nextSlice)
-                               vLumBufSize= nextSlice - vLumFilterPos[i   ];
-                       if(vChrFilterPos[chrI] + vChrBufSize < (nextSlice>>1))
-                               vChrBufSize= (nextSlice>>1) - vChrFilterPos[chrI];
-               }
-
-               // allocate pixbufs (we use dynamic allocation because otherwise we would need to
-               // allocate several megabytes to handle all possible cases)
-               for(i=0; i<vLumBufSize; i++)
-                       lumPixBuf[i]= lumPixBuf[i+vLumBufSize]= (uint16_t*)memalign(8, 4000);
-               for(i=0; i<vChrBufSize; i++)
-                       chrPixBuf[i]= chrPixBuf[i+vChrBufSize]= (uint16_t*)memalign(8, 8000);
-
-               //try to avoid drawing green stuff between the right end and the stride end
-               for(i=0; i<vLumBufSize; i++) memset(lumPixBuf[i], 0, 4000);
-               for(i=0; i<vChrBufSize; i++) memset(chrPixBuf[i], 64, 8000);
-
-               ASSERT(chrDstH<=dstH)
-               ASSERT(vLumFilterSize*dstH*4<16000)
-               ASSERT(vChrFilterSize*chrDstH*4<16000)
-#ifdef HAVE_MMX
-               // pack filter data for mmx code
-               for(i=0; i<vLumFilterSize*dstH; i++)
-                       lumMmxFilter[4*i]=lumMmxFilter[4*i+1]=lumMmxFilter[4*i+2]=lumMmxFilter[4*i+3]=
-                               vLumFilter[i];
-               for(i=0; i<vChrFilterSize*chrDstH; i++)
-                       chrMmxFilter[4*i]=chrMmxFilter[4*i+1]=chrMmxFilter[4*i+2]=chrMmxFilter[4*i+3]=
-                               vChrFilter[i];
-#endif
        }
 
-       if(firstTime && verbose)
-       {
-#ifdef HAVE_MMX2
-               int mmx2=1;
-#else
-               int mmx2=0;
-#endif
-#ifdef HAVE_MMX
-               int mmx=1;
-#else
-               int mmx=0;
-#endif
-
-#ifdef HAVE_MMX
-               if(canMMX2BeUsed && sws_flags==SWS_FAST_BILINEAR)
-                       printf("SwScaler: using FAST_BILINEAR MMX2 scaler for horizontal scaling\n");
-               else
-               {
-                       if(hLumFilterSize==4)
-                               printf("SwScaler: using 4-tap MMX scaler for horizontal luminance scaling\n");
-                       else if(hLumFilterSize==8)
-                               printf("SwScaler: using 8-tap MMX scaler for horizontal luminance scaling\n");
-                       else
-                               printf("SwScaler: using n-tap MMX scaler for horizontal luminance scaling\n");
-
-                       if(hChrFilterSize==4)
-                               printf("SwScaler: using 4-tap MMX scaler for horizontal chrominance scaling\n");
-                       else if(hChrFilterSize==8)
-                               printf("SwScaler: using 8-tap MMX scaler for horizontal chrominance scaling\n");
-                       else
-                               printf("SwScaler: using n-tap MMX scaler for horizontal chrominance scaling\n");
-               }
-#elif defined (ARCH_X86)
-               printf("SwScaler: using X86-Asm scaler for horizontal scaling\n");
-#else
-               if(sws_flags==SWS_FAST_BILINEAR)
-                       printf("SwScaler: using FAST_BILINEAR C scaler for horizontal scaling\n");
-               else
-                       printf("SwScaler: using C scaler for horizontal scaling\n");
-#endif
-
-               if(dstbpp==12)
-               {
-                       if(vLumFilterSize==1)
-                               printf("SwScaler: using 1-tap %s \"scaler\" for vertical scaling (YV12)\n", mmx ? "MMX" : "C");
-                       else
-                               printf("SwScaler: using n-tap %s scaler for vertical scaling (YV12)\n", mmx ? "MMX" : "C");
-               }
-               else
-               {
-                       if(vLumFilterSize==1 && vChrFilterSize==2)
-                               printf("SwScaler: using 1-tap %s \"scaler\" for vertical luminance scaling (BGR)\n"
-                                      "SwScaler:       2-tap scaler for vertical chrominance scaling (BGR)\n", mmx ? "MMX" : "C");
-                       else if(vLumFilterSize==2 && vChrFilterSize==2)
-                               printf("SwScaler: using 2-tap linear %s scaler for vertical scaling (BGR)\n", mmx ? "MMX" : "C");
-                       else
-                               printf("SwScaler: using n-tap %s scaler for vertical scaling (BGR)\n", mmx ? "MMX" : "C");
-               }
-
-               if(dstbpp==24)
-                       printf("SwScaler: using %s YV12->BGR24 Converter\n",
-                               mmx2 ? "MMX2" : (mmx ? "MMX" : "C"));
-               else
-                       printf("SwScaler: using %s YV12->BGR%d Converter\n", mmx ? "MMX" : "C", dstbpp);
-
-               printf("SwScaler: %dx%d -> %dx%d\n", srcW, srcH, dstW, dstH);
+       /* Note the user might start scaling the picture in the middle so this will not get executed
+          this is not really intended but works currently, so ppl might do it */
+       if(srcSliceY ==0){
+               lumBufIndex=0;
+               chrBufIndex=0;
+               dstY=0; 
+               lastInLumBuf= -1;
+               lastInChrBuf= -1;
        }
 
-       lastInLumBuf= -1;
-       lastInChrBuf= -1;
-  } // if(firstLine)
+       lastDstY= dstY;
 
        for(;dstY < dstH; dstY++){
-               unsigned char *dest =dstptr[0]+dststride*dstY;
-               unsigned char *uDest=dstptr[1]+(dststride>>1)*(dstY>>1);
-               unsigned char *vDest=dstptr[2]+(dststride>>1)*(dstY>>1);
-               const int chrDstY= dstbpp==12 ? (dstY>>1) : dstY;
+               unsigned char *dest =dst[0]+dstStride[0]*dstY;
+               const int chrDstY= dstY>>c->chrDstVSubSample;
+               unsigned char *uDest=dst[1]+dstStride[1]*chrDstY;
+               unsigned char *vDest=dst[2]+dstStride[2]*chrDstY;
 
                const int firstLumSrcY= vLumFilterPos[dstY]; //First line needed as input
                const int firstChrSrcY= vChrFilterPos[chrDstY]; //First line needed as input
                const int lastLumSrcY= firstLumSrcY + vLumFilterSize -1; // Last line needed as input
                const int lastChrSrcY= firstChrSrcY + vChrFilterSize -1; // Last line needed as input
 
-               if(sws_flags == SWS_FAST_BILINEAR)
-               {
-                       //handle holes
-                       if(firstLumSrcY > lastInLumBuf) lastInLumBuf= firstLumSrcY-1;
-                       if(firstChrSrcY > lastInChrBuf) lastInChrBuf= firstChrSrcY-1;
-               }
-
+//printf("dstY:%d dstH:%d firstLumSrcY:%d lastInLumBuf:%d vLumBufSize: %d vChrBufSize: %d slice: %d %d vLumFilterSize: %d firstChrSrcY: %d vChrFilterSize: %d c->chrSrcVSubSample: %d\n",
+// dstY, dstH, firstLumSrcY, lastInLumBuf, vLumBufSize, vChrBufSize, srcSliceY, srcSliceH, vLumFilterSize, firstChrSrcY, vChrFilterSize,  c->chrSrcVSubSample);
+               //handle holes (FAST_BILINEAR & weird filters)
+               if(firstLumSrcY > lastInLumBuf) lastInLumBuf= firstLumSrcY-1;
+               if(firstChrSrcY > lastInChrBuf) lastInChrBuf= firstChrSrcY-1;
+//printf("%d %d %d\n", firstChrSrcY, lastInChrBuf, vChrBufSize);
                ASSERT(firstLumSrcY >= lastInLumBuf - vLumBufSize + 1)
                ASSERT(firstChrSrcY >= lastInChrBuf - vChrBufSize + 1)
 
                // Do we have enough lines in this slice to output the dstY line
-               if(lastLumSrcY < srcSliceY + srcSliceH && lastChrSrcY < ((srcSliceY + srcSliceH)>>1))
+               if(lastLumSrcY < srcSliceY + srcSliceH && lastChrSrcY < -((-srcSliceY - srcSliceH)>>c->chrSrcVSubSample))
                {
                        //Do horizontal scaling
                        while(lastInLumBuf < lastLumSrcY)
                        {
-                               uint8_t *src= srcptr[0]+(lastInLumBuf + 1 - srcSliceY)*stride[0];
+                               uint8_t *s= src[0]+(lastInLumBuf + 1 - srcSliceY)*srcStride[0];
                                lumBufIndex++;
+//                             printf("%d %d %d %d\n", lumBufIndex, vLumBufSize, lastInLumBuf,  lastLumSrcY);
                                ASSERT(lumBufIndex < 2*vLumBufSize)
                                ASSERT(lastInLumBuf + 1 - srcSliceY < srcSliceH)
                                ASSERT(lastInLumBuf + 1 - srcSliceY >= 0)
 //                             printf("%d %d\n", lumBufIndex, vLumBufSize);
-                               RENAME(hyscale)(lumPixBuf[ lumBufIndex ], dstW, src, srcW, lumXInc);
+                               RENAME(hyscale)(lumPixBuf[ lumBufIndex ], dstW, s, srcW, lumXInc,
+                                               flags, canMMX2BeUsed, hLumFilter, hLumFilterPos, hLumFilterSize,
+                                               funnyYCode, c->srcFormat, formatConvBuffer, 
+                                               c->lumMmx2Filter, c->lumMmx2FilterPos);
                                lastInLumBuf++;
                        }
                        while(lastInChrBuf < lastChrSrcY)
                        {
-                               uint8_t *src1= srcptr[1]+(lastInChrBuf + 1 - (srcSliceY>>1))*stride[1];
-                               uint8_t *src2= srcptr[2]+(lastInChrBuf + 1 - (srcSliceY>>1))*stride[2];
+                               uint8_t *src1= src[1]+(lastInChrBuf + 1 - chrSrcSliceY)*srcStride[1];
+                               uint8_t *src2= src[2]+(lastInChrBuf + 1 - chrSrcSliceY)*srcStride[2];
                                chrBufIndex++;
                                ASSERT(chrBufIndex < 2*vChrBufSize)
-                               ASSERT(lastInChrBuf + 1 - (srcSliceY>>1) < (srcSliceH>>1))
-                               ASSERT(lastInChrBuf + 1 - (srcSliceY>>1) >= 0)
-                               RENAME(hcscale)(chrPixBuf[ chrBufIndex ], chrDstW, src1, src2, (srcW+1)>>1, chrXInc);
+                               ASSERT(lastInChrBuf + 1 - chrSrcSliceY < (chrSrcSliceH))
+                               ASSERT(lastInChrBuf + 1 - chrSrcSliceY >= 0)
+                               //FIXME replace parameters through context struct (some at least)
+
+                               if(!(isGray(srcFormat) || isGray(dstFormat)))
+                                       RENAME(hcscale)(chrPixBuf[ chrBufIndex ], chrDstW, src1, src2, chrSrcW, chrXInc,
+                                               flags, canMMX2BeUsed, hChrFilter, hChrFilterPos, hChrFilterSize,
+                                               funnyUVCode, c->srcFormat, formatConvBuffer, 
+                                               c->chrMmx2Filter, c->chrMmx2FilterPos);
                                lastInChrBuf++;
                        }
                        //wrap buf index around to stay inside the ring buffer
@@ -2481,34 +2665,42 @@ else            chrYInc= lumYInc>>1, chrDstH= dstH;
 /*             printf("%d %d Last:%d %d LastInBuf:%d %d Index:%d %d Y:%d FSize: %d %d BSize: %d %d\n",
                        firstChrSrcY,firstLumSrcY,lastChrSrcY,lastLumSrcY,
                        lastInChrBuf,lastInLumBuf,chrBufIndex,lumBufIndex,dstY,vChrFilterSize,vLumFilterSize,
-                       vChrBufSize, vLumBufSize);
-*/
+                       vChrBufSize, vLumBufSize);*/
+
                        //Do horizontal scaling
                        while(lastInLumBuf+1 < srcSliceY + srcSliceH)
                        {
-                               uint8_t *src= srcptr[0]+(lastInLumBuf + 1 - srcSliceY)*stride[0];
+                               uint8_t *s= src[0]+(lastInLumBuf + 1 - srcSliceY)*srcStride[0];
                                lumBufIndex++;
                                ASSERT(lumBufIndex < 2*vLumBufSize)
                                ASSERT(lastInLumBuf + 1 - srcSliceY < srcSliceH)
                                ASSERT(lastInLumBuf + 1 - srcSliceY >= 0)
-                               RENAME(hyscale)(lumPixBuf[ lumBufIndex ], dstW, src, srcW, lumXInc);
+                               RENAME(hyscale)(lumPixBuf[ lumBufIndex ], dstW, s, srcW, lumXInc,
+                                               flags, canMMX2BeUsed, hLumFilter, hLumFilterPos, hLumFilterSize,
+                                               funnyYCode, c->srcFormat, formatConvBuffer, 
+                                               c->lumMmx2Filter, c->lumMmx2FilterPos);
                                lastInLumBuf++;
                        }
-                       while(lastInChrBuf+1 < ((srcSliceY + srcSliceH)>>1))
+                       while(lastInChrBuf+1 < (chrSrcSliceY + chrSrcSliceH))
                        {
-                               uint8_t *src1= srcptr[1]+(lastInChrBuf + 1 - (srcSliceY>>1))*stride[1];
-                               uint8_t *src2= srcptr[2]+(lastInChrBuf + 1 - (srcSliceY>>1))*stride[2];
+                               uint8_t *src1= src[1]+(lastInChrBuf + 1 - chrSrcSliceY)*srcStride[1];
+                               uint8_t *src2= src[2]+(lastInChrBuf + 1 - chrSrcSliceY)*srcStride[2];
                                chrBufIndex++;
                                ASSERT(chrBufIndex < 2*vChrBufSize)
-                               ASSERT(lastInChrBuf + 1 - (srcSliceY>>1) < (srcSliceH>>1))
-                               ASSERT(lastInChrBuf + 1 - (srcSliceY>>1) >= 0)
-                               RENAME(hcscale)(chrPixBuf[ chrBufIndex ], chrDstW, src1, src2, (srcW+1)>>1, chrXInc);
+                               ASSERT(lastInChrBuf + 1 - chrSrcSliceY < chrSrcSliceH)
+                               ASSERT(lastInChrBuf + 1 - chrSrcSliceY >= 0)
+
+                               if(!(isGray(srcFormat) || isGray(dstFormat)))
+                                       RENAME(hcscale)(chrPixBuf[ chrBufIndex ], chrDstW, src1, src2, chrSrcW, chrXInc,
+                                               flags, canMMX2BeUsed, hChrFilter, hChrFilterPos, hChrFilterSize,
+                                               funnyUVCode, c->srcFormat, formatConvBuffer, 
+                                               c->chrMmx2Filter, c->chrMmx2FilterPos);
                                lastInChrBuf++;
                        }
                        //wrap buf index around to stay inside the ring buffer
                        if(lumBufIndex >= vLumBufSize ) lumBufIndex-= vLumBufSize;
                        if(chrBufIndex >= vChrBufSize ) chrBufIndex-= vChrBufSize;
-                       break; //we cant output a dstY line so lets try with the next slice
+                       break; //we can't output a dstY line so let's try with the next slice
                }
 
 #ifdef HAVE_MMX
@@ -2517,80 +2709,92 @@ else            chrYInc= lumYInc>>1, chrDstH= dstH;
                g5Dither= dither8[dstY&1];
                r5Dither= dither8[(dstY+1)&1];
 #endif
-           if(dstY < dstH-2 || over<=0)
+           if(dstY < dstH-2)
            {
-               if(dstbpp==12) //YV12
+               int16_t **lumSrcPtr= lumPixBuf + lumBufIndex + firstLumSrcY - lastInLumBuf + vLumBufSize;
+               int16_t **chrSrcPtr= chrPixBuf + chrBufIndex + firstChrSrcY - lastInChrBuf + vChrBufSize;
+#ifdef HAVE_MMX
+               int i;
+               for(i=0; i<vLumFilterSize; i++)
+               {
+                       lumMmxFilter[4*i+0]= (int32_t)lumSrcPtr[i];
+                       lumMmxFilter[4*i+2]= 
+                       lumMmxFilter[4*i+3]= 
+                               ((uint16_t)vLumFilter[dstY*vLumFilterSize + i])*0x10001;
+               }
+               for(i=0; i<vChrFilterSize; i++)
                {
-                       if(dstY&1) uDest=vDest= NULL; //FIXME split functions in lumi / chromi
+                       chrMmxFilter[4*i+0]= (int32_t)chrSrcPtr[i];
+                       chrMmxFilter[4*i+2]= 
+                       chrMmxFilter[4*i+3]= 
+                               ((uint16_t)vChrFilter[chrDstY*vChrFilterSize + i])*0x10001;
+               }
+#endif
+               if(isPlanarYUV(dstFormat) || isGray(dstFormat)) //YV12 like
+               {
+                       const int chrSkipMask= (1<<c->chrDstVSubSample)-1;
+                       if((dstY&chrSkipMask) || isGray(dstFormat)) uDest=vDest= NULL; //FIXME split functions in lumi / chromi
                        if(vLumFilterSize == 1 && vChrFilterSize == 1) // Unscaled YV12
                        {
                                int16_t *lumBuf = lumPixBuf[0];
                                int16_t *chrBuf= chrPixBuf[0];
-                               RENAME(yuv2yuv1)(lumBuf, chrBuf, dest, uDest, vDest, dstW);
+                               RENAME(yuv2yuv1)(lumBuf, chrBuf, dest, uDest, vDest, dstW, chrDstW);
                        }
                        else //General YV12
                        {
-                               int16_t **lumSrcPtr= lumPixBuf + lumBufIndex + firstLumSrcY - lastInLumBuf + vLumBufSize;
-                               int16_t **chrSrcPtr= chrPixBuf + chrBufIndex + firstChrSrcY - lastInChrBuf + vChrBufSize;
-                               RENAME(yuv2yuvX)(
-                                       vLumFilter+dstY*vLumFilterSize     , lumSrcPtr, vLumFilterSize,
-                                       vChrFilter+(dstY>>1)*vChrFilterSize, chrSrcPtr, vChrFilterSize,
-                                       dest, uDest, vDest, dstW,
-                                       lumMmxFilter+dstY*vLumFilterSize*4, chrMmxFilter+(dstY>>1)*vChrFilterSize*4);
+                               RENAME(yuv2yuvX)(c,
+                                       vLumFilter+dstY*vLumFilterSize   , lumSrcPtr, vLumFilterSize,
+                                       vChrFilter+chrDstY*vChrFilterSize, chrSrcPtr, vChrFilterSize,
+                                       dest, uDest, vDest, dstW, chrDstW);
                        }
                }
                else
                {
-                       int16_t **lumSrcPtr= lumPixBuf + lumBufIndex + firstLumSrcY - lastInLumBuf + vLumBufSize;
-                       int16_t **chrSrcPtr= chrPixBuf + chrBufIndex + firstChrSrcY - lastInChrBuf + vChrBufSize;
-
                        ASSERT(lumSrcPtr + vLumFilterSize - 1 < lumPixBuf + vLumBufSize*2);
                        ASSERT(chrSrcPtr + vChrFilterSize - 1 < chrPixBuf + vChrBufSize*2);
                        if(vLumFilterSize == 1 && vChrFilterSize == 2) //Unscaled RGB
                        {
                                int chrAlpha= vChrFilter[2*dstY+1];
-
-                               RENAME(yuv2rgb1)(*lumSrcPtr, *chrSrcPtr, *(chrSrcPtr+1),
-                                                dest, dstW, chrAlpha, dstbpp);
+                               RENAME(yuv2packed1)(c, *lumSrcPtr, *chrSrcPtr, *(chrSrcPtr+1),
+                                                dest, dstW, chrAlpha, dstFormat, flags, dstY);
                        }
                        else if(vLumFilterSize == 2 && vChrFilterSize == 2) //BiLinear Upscale RGB
                        {
                                int lumAlpha= vLumFilter[2*dstY+1];
                                int chrAlpha= vChrFilter[2*dstY+1];
-
-                               RENAME(yuv2rgb2)(*lumSrcPtr, *(lumSrcPtr+1), *chrSrcPtr, *(chrSrcPtr+1),
-                                                dest, dstW, lumAlpha, chrAlpha, dstbpp);
+                               RENAME(yuv2packed2)(c, *lumSrcPtr, *(lumSrcPtr+1), *chrSrcPtr, *(chrSrcPtr+1),
+                                                dest, dstW, lumAlpha, chrAlpha, dstY);
                        }
                        else //General RGB
                        {
-                               RENAME(yuv2rgbX)(
+                               RENAME(yuv2packedX)(c,
                                        vLumFilter+dstY*vLumFilterSize, lumSrcPtr, vLumFilterSize,
                                        vChrFilter+dstY*vChrFilterSize, chrSrcPtr, vChrFilterSize,
-                                       dest, dstW, dstbpp,
-                                       lumMmxFilter+dstY*vLumFilterSize*4, chrMmxFilter+dstY*vChrFilterSize*4);
+                                       dest, dstW, dstY);
                        }
                }
             }
-           else // hmm looks like we cant use MMX here without overwriting this arrays tail
+           else // hmm looks like we can't use MMX here without overwriting this array's tail
            {
                int16_t **lumSrcPtr= lumPixBuf + lumBufIndex + firstLumSrcY - lastInLumBuf + vLumBufSize;
                int16_t **chrSrcPtr= chrPixBuf + chrBufIndex + firstChrSrcY - lastInChrBuf + vChrBufSize;
-               if(dstbpp==12) //YV12
+               if(isPlanarYUV(dstFormat) || isGray(dstFormat)) //YV12
                {
-                       if(dstY&1) uDest=vDest= NULL; //FIXME split functions in lumi / chromi
+                       const int chrSkipMask= (1<<c->chrDstVSubSample)-1;
+                       if((dstY&chrSkipMask) || isGray(dstFormat)) uDest=vDest= NULL; //FIXME split functions in lumi / chromi
                        yuv2yuvXinC(
-                               vLumFilter+dstY*vLumFilterSize     , lumSrcPtr, vLumFilterSize,
-                               vChrFilter+(dstY>>1)*vChrFilterSize, chrSrcPtr, vChrFilterSize,
-                               dest, uDest, vDest, dstW);
+                               vLumFilter+dstY*vLumFilterSize   , lumSrcPtr, vLumFilterSize,
+                               vChrFilter+chrDstY*vChrFilterSize, chrSrcPtr, vChrFilterSize,
+                               dest, uDest, vDest, dstW, chrDstW);
                }
                else
                {
                        ASSERT(lumSrcPtr + vLumFilterSize - 1 < lumPixBuf + vLumBufSize*2);
                        ASSERT(chrSrcPtr + vChrFilterSize - 1 < chrPixBuf + vChrBufSize*2);
-                       yuv2rgbXinC(
+                       yuv2packedXinC(c, 
                                vLumFilter+dstY*vLumFilterSize, lumSrcPtr, vLumFilterSize,
                                vChrFilter+dstY*vChrFilterSize, chrSrcPtr, vChrFilterSize,
-                               dest, dstW, dstbpp);
+                               dest, dstW, dstY);
                }
            }
        }
@@ -2599,5 +2803,12 @@ else             chrYInc= lumYInc>>1, chrDstH= dstH;
        __asm __volatile(SFENCE:::"memory");
        __asm __volatile(EMMS:::"memory");
 #endif
-       firstTime=0;
+       /* store changed local vars back in the context */
+       c->dstY= dstY;
+       c->lumBufIndex= lumBufIndex;
+       c->chrBufIndex= chrBufIndex;
+       c->lastInLumBuf= lastInLumBuf;
+       c->lastInChrBuf= lastInChrBuf;
+
+       return dstY - lastDstY;
 }