]> git.sesse.net Git - ffmpeg/blobdiff - postproc/rgb2rgb_template.c
fixing "internal error in orderYUV" bug
[ffmpeg] / postproc / rgb2rgb_template.c
index b873462e18147f04e0934225284071a306bd59d0..85c8f440e3968067cf10d4e122013e0289bd0cae 100644 (file)
 /*
  *
  *  rgb2rgb.c, Software RGB to RGB convertor
+ *  pluralize by Software PAL8 to RGB convertor
+ *               Software YUV to YUV convertor
+ *               Software YUV to RGB convertor
  *  Written by Nick Kurshev.
- *  palette stuff & yuv stuff by Michael
+ *  palette & yuv & runtime cpu stuff by Michael (michaelni@gmx.at) (under GPL)
  */
-#include <inttypes.h>
-#include "../config.h"
-#include "rgb2rgb.h"
-#include "../mmx_defs.h"
 
-#ifdef HAVE_MMX
-static const uint64_t mask32   __attribute__((aligned(8))) = 0x00FFFFFF00FFFFFFULL;
-static const uint64_t mask24l  __attribute__((aligned(8))) = 0x0000000000FFFFFFULL;
-static const uint64_t mask24h  __attribute__((aligned(8))) = 0x0000FFFFFF000000ULL;
-static const uint64_t mask15b  __attribute__((aligned(8))) = 0x001F001F001F001FULL; /* 00000000 00011111  xxB */
-static const uint64_t mask15rg __attribute__((aligned(8))) = 0x7FE07FE07FE07FE0ULL; /* 01111111 11100000  RGx */
-static const uint64_t mask15s  __attribute__((aligned(8))) = 0xFFE0FFE0FFE0FFE0ULL;
+#include <stddef.h>
+#include <inttypes.h> /* for __WORDSIZE */
+
+#ifndef __WORDSIZE
+// #warning You have misconfigured system and probably will lose performance!
+#define __WORDSIZE MP_WORDSIZE
+#endif
+
+#undef PREFETCH
+#undef MOVNTQ
+#undef EMMS
+#undef SFENCE
+#undef MMREG_SIZE
+#undef PREFETCHW
+#undef PAVGB
+
+#ifdef HAVE_SSE2
+#define MMREG_SIZE 16
+#else
+#define MMREG_SIZE 8
+#endif
+
+#ifdef HAVE_3DNOW
+#define PREFETCH  "prefetch"
+#define PREFETCHW "prefetchw"
+#define PAVGB    "pavgusb"
+#elif defined ( HAVE_MMX2 )
+#define PREFETCH "prefetchnta"
+#define PREFETCHW "prefetcht0"
+#define PAVGB    "pavgb"
+#else
+#define PREFETCH "/nop"
+#define PREFETCHW "/nop"
+#endif
+
+#ifdef HAVE_3DNOW
+/* On K6 femms is faster of emms. On K7 femms is directly mapped on emms. */
+#define EMMS     "femms"
+#else
+#define EMMS     "emms"
+#endif
+
+#ifdef HAVE_MMX2
+#define MOVNTQ "movntq"
+#define SFENCE "sfence"
+#else
+#define MOVNTQ "movq"
+#define SFENCE "/nop"
 #endif
 
-void rgb24to32(const uint8_t *src,uint8_t *dst,uint32_t src_size)
+static inline void RENAME(rgb24to32)(const uint8_t *src,uint8_t *dst,unsigned src_size)
 {
   uint8_t *dest = dst;
   const uint8_t *s = src;
   const uint8_t *end;
 #ifdef HAVE_MMX
-  uint8_t *mm_end;
+  const uint8_t *mm_end;
 #endif
   end = s + src_size;
 #ifdef HAVE_MMX
   __asm __volatile(PREFETCH"   %0"::"m"(*s):"memory");
-  mm_end = (uint8_t*)((((unsigned long)end)/(MMREG_SIZE*2))*(MMREG_SIZE*2));
+  mm_end = end - 23;
   __asm __volatile("movq       %0, %%mm7"::"m"(mask32):"memory");
-  if(mm_end == end) mm_end -= MMREG_SIZE*2;
   while(s < mm_end)
   {
     __asm __volatile(
        PREFETCH"       32%1\n\t"
        "movd   %1, %%mm0\n\t"
-       "movd   3%1, %%mm1\n\t"
-       "movd   6%1, %%mm2\n\t"
-       "movd   9%1, %%mm3\n\t"
-       "punpckldq %%mm1, %%mm0\n\t"
-       "punpckldq %%mm3, %%mm2\n\t"
+       "punpckldq 3%1, %%mm0\n\t"
+       "movd   6%1, %%mm1\n\t"
+       "punpckldq 9%1, %%mm1\n\t"
+       "movd   12%1, %%mm2\n\t"
+       "punpckldq 15%1, %%mm2\n\t"
+       "movd   18%1, %%mm3\n\t"
+       "punpckldq 21%1, %%mm3\n\t"
        "pand   %%mm7, %%mm0\n\t"
+       "pand   %%mm7, %%mm1\n\t"
        "pand   %%mm7, %%mm2\n\t"
+       "pand   %%mm7, %%mm3\n\t"
        MOVNTQ" %%mm0, %0\n\t"
-       MOVNTQ" %%mm2, 8%0"
+       MOVNTQ" %%mm1, 8%0\n\t"
+       MOVNTQ" %%mm2, 16%0\n\t"
+       MOVNTQ" %%mm3, 24%0"
        :"=m"(*dest)
        :"m"(*s)
        :"memory");
-    dest += 16;
-    s += 12;
+    dest += 32;
+    s += 24;
   }
   __asm __volatile(SFENCE:::"memory");
   __asm __volatile(EMMS:::"memory");
@@ -64,46 +109,70 @@ void rgb24to32(const uint8_t *src,uint8_t *dst,uint32_t src_size)
   }
 }
 
-void rgb32to24(const uint8_t *src,uint8_t *dst,uint32_t src_size)
+static inline void RENAME(rgb32to24)(const uint8_t *src,uint8_t *dst,unsigned src_size)
 {
   uint8_t *dest = dst;
   const uint8_t *s = src;
   const uint8_t *end;
 #ifdef HAVE_MMX
-  uint8_t *mm_end;
+  const uint8_t *mm_end;
 #endif
   end = s + src_size;
 #ifdef HAVE_MMX
   __asm __volatile(PREFETCH"   %0"::"m"(*s):"memory");
-  mm_end = (uint8_t*)((((unsigned long)end)/(MMREG_SIZE*2))*(MMREG_SIZE*2));
-  __asm __volatile(
-       "movq   %0, %%mm7\n\t"
-       "movq   %1, %%mm6"
-       ::"m"(mask24l),"m"(mask24h):"memory");
-  if(mm_end == end) mm_end -= MMREG_SIZE*2;
+  mm_end = end - 31;
   while(s < mm_end)
   {
     __asm __volatile(
        PREFETCH"       32%1\n\t"
        "movq   %1, %%mm0\n\t"
        "movq   8%1, %%mm1\n\t"
+       "movq   16%1, %%mm4\n\t"
+       "movq   24%1, %%mm5\n\t"
        "movq   %%mm0, %%mm2\n\t"
        "movq   %%mm1, %%mm3\n\t"
+       "movq   %%mm4, %%mm6\n\t"
+       "movq   %%mm5, %%mm7\n\t"
        "psrlq  $8, %%mm2\n\t"
        "psrlq  $8, %%mm3\n\t"
-       "pand   %%mm7, %%mm0\n\t"
-       "pand   %%mm7, %%mm1\n\t"
-       "pand   %%mm6, %%mm2\n\t"
-       "pand   %%mm6, %%mm3\n\t"
+       "psrlq  $8, %%mm6\n\t"
+       "psrlq  $8, %%mm7\n\t"
+       "pand   %2, %%mm0\n\t"
+       "pand   %2, %%mm1\n\t"
+       "pand   %2, %%mm4\n\t"
+       "pand   %2, %%mm5\n\t"
+       "pand   %3, %%mm2\n\t"
+       "pand   %3, %%mm3\n\t"
+       "pand   %3, %%mm6\n\t"
+       "pand   %3, %%mm7\n\t"
+       "por    %%mm2, %%mm0\n\t"
+       "por    %%mm3, %%mm1\n\t"
+       "por    %%mm6, %%mm4\n\t"
+       "por    %%mm7, %%mm5\n\t"
+
+       "movq   %%mm1, %%mm2\n\t"
+       "movq   %%mm4, %%mm3\n\t"
+       "psllq  $48, %%mm2\n\t"
+       "psllq  $32, %%mm3\n\t"
+       "pand   %4, %%mm2\n\t"
+       "pand   %5, %%mm3\n\t"
        "por    %%mm2, %%mm0\n\t"
+       "psrlq  $16, %%mm1\n\t"
+       "psrlq  $32, %%mm4\n\t"
+       "psllq  $16, %%mm5\n\t"
        "por    %%mm3, %%mm1\n\t"
+       "pand   %6, %%mm5\n\t"
+       "por    %%mm5, %%mm4\n\t"
+
        MOVNTQ" %%mm0, %0\n\t"
-       MOVNTQ" %%mm1, 6%0"
+       MOVNTQ" %%mm1, 8%0\n\t"
+       MOVNTQ" %%mm4, 16%0"
        :"=m"(*dest)
-       :"m"(*s)
+       :"m"(*s),"m"(mask24l),
+        "m"(mask24h),"m"(mask24hh),"m"(mask24hhh),"m"(mask24hhhh)
        :"memory");
-    dest += 12;
-    s += 16;
+    dest += 24;
+    s += 32;
   }
   __asm __volatile(SFENCE:::"memory");
   __asm __volatile(EMMS:::"memory");
@@ -123,17 +192,18 @@ void rgb32to24(const uint8_t *src,uint8_t *dst,uint32_t src_size)
  MMX2, 3DNOW optimization by Nick Kurshev
  32bit c version, and and&add trick by Michael Niedermayer
 */
-void rgb15to16(const uint8_t *src,uint8_t *dst,uint32_t src_size)
+static inline void RENAME(rgb15to16)(const uint8_t *src,uint8_t *dst,unsigned src_size)
 {
+  register const uint8_t* s=src;
+  register uint8_t* d=dst;
+  register const uint8_t *end;
+  const uint8_t *mm_end;
+  end = s + src_size;
 #ifdef HAVE_MMX
-  register const char* s=src+src_size;
-  register char* d=dst+src_size;
-  register int offs=-src_size;
-  __asm __volatile(PREFETCH"   %0"::"m"(*(s+offs)));
-  __asm __volatile(
-       "movq   %0, %%mm4\n\t"
-       ::"m"(mask15s));
-  while(offs<0)
+  __asm __volatile(PREFETCH"   %0"::"m"(*s));
+  __asm __volatile("movq       %0, %%mm4"::"m"(mask15s));
+  mm_end = end - 15;
+  while(s<mm_end)
   {
        __asm __volatile(
                PREFETCH"       32%1\n\t"
@@ -147,238 +217,2306 @@ void rgb15to16(const uint8_t *src,uint8_t *dst,uint32_t src_size)
                "paddw  %%mm3, %%mm2\n\t"
                MOVNTQ" %%mm0, %0\n\t"
                MOVNTQ" %%mm2, 8%0"
-               :"=m"(*(d+offs))
-               :"m"(*(s+offs))
+               :"=m"(*d)
+               :"m"(*s)
                );
-       offs+=16;
+       d+=16;
+       s+=16;
   }
   __asm __volatile(SFENCE:::"memory");
   __asm __volatile(EMMS:::"memory");
+#endif
+    mm_end = end - 3;
+    while(s < mm_end)
+    {
+       register unsigned x= *((uint32_t *)s);
+       *((uint32_t *)d) = (x&0x7FFF7FFF) + (x&0x7FE07FE0);
+       d+=4;
+       s+=4;
+    }
+    if(s < end)
+    {
+       register unsigned short x= *((uint16_t *)s);
+       *((uint16_t *)d) = (x&0x7FFF) + (x&0x7FE0);
+    }
+}
+
+static inline void RENAME(rgb16to15)(const uint8_t *src,uint8_t *dst,unsigned src_size)
+{
+  register const uint8_t* s=src;
+  register uint8_t* d=dst;
+  register const uint8_t *end;
+  const uint8_t *mm_end;
+  end = s + src_size;
+#ifdef HAVE_MMX
+  __asm __volatile(PREFETCH"   %0"::"m"(*s));
+  __asm __volatile("movq       %0, %%mm7"::"m"(mask15rg));
+  __asm __volatile("movq       %0, %%mm6"::"m"(mask15b));
+  mm_end = end - 15;
+  while(s<mm_end)
+  {
+       __asm __volatile(
+               PREFETCH"       32%1\n\t"
+               "movq   %1, %%mm0\n\t"
+               "movq   8%1, %%mm2\n\t"
+               "movq   %%mm0, %%mm1\n\t"
+               "movq   %%mm2, %%mm3\n\t"
+               "psrlq  $1, %%mm0\n\t"
+               "psrlq  $1, %%mm2\n\t"
+               "pand   %%mm7, %%mm0\n\t"
+               "pand   %%mm7, %%mm2\n\t"
+               "pand   %%mm6, %%mm1\n\t"
+               "pand   %%mm6, %%mm3\n\t"
+               "por    %%mm1, %%mm0\n\t"
+               "por    %%mm3, %%mm2\n\t"
+               MOVNTQ" %%mm0, %0\n\t"
+               MOVNTQ" %%mm2, 8%0"
+               :"=m"(*d)
+               :"m"(*s)
+               );
+       d+=16;
+       s+=16;
+  }
+  __asm __volatile(SFENCE:::"memory");
+  __asm __volatile(EMMS:::"memory");
+#endif
+    mm_end = end - 3;
+    while(s < mm_end)
+    {
+       register uint32_t x= *((uint32_t *)s);
+       *((uint32_t *)d) = ((x>>1)&0x7FE07FE0) | (x&0x001F001F);
+       s+=4;
+       d+=4;
+    }
+    if(s < end)
+    {
+       register uint16_t x= *((uint16_t *)s);
+       *((uint16_t *)d) = ((x>>1)&0x7FE0) | (x&0x001F);
+       s+=2;
+       d+=2;
+    }
+}
+
+static inline void RENAME(rgb32to16)(const uint8_t *src, uint8_t *dst, unsigned src_size)
+{
+       const uint8_t *s = src;
+       const uint8_t *end;
+#ifdef HAVE_MMX
+       const uint8_t *mm_end;
+#endif
+       uint16_t *d = (uint16_t *)dst;
+       end = s + src_size;
+#ifdef HAVE_MMX
+       mm_end = end - 15;
+#if 1 //is faster only if multiplies are reasonable fast (FIXME figure out on which cpus this is faster, on Athlon its slightly faster)
+       asm volatile(
+               "movq %3, %%mm5                 \n\t"
+               "movq %4, %%mm6                 \n\t"
+               "movq %5, %%mm7                 \n\t"
+               ".balign 16                     \n\t"
+               "1:                             \n\t"
+               PREFETCH" 32(%1)                \n\t"
+               "movd   (%1), %%mm0             \n\t"
+               "movd   4(%1), %%mm3            \n\t"
+               "punpckldq 8(%1), %%mm0         \n\t"
+               "punpckldq 12(%1), %%mm3        \n\t"
+               "movq %%mm0, %%mm1              \n\t"
+               "movq %%mm3, %%mm4              \n\t"
+               "pand %%mm6, %%mm0              \n\t"
+               "pand %%mm6, %%mm3              \n\t"
+               "pmaddwd %%mm7, %%mm0           \n\t"
+               "pmaddwd %%mm7, %%mm3           \n\t"
+               "pand %%mm5, %%mm1              \n\t"
+               "pand %%mm5, %%mm4              \n\t"
+               "por %%mm1, %%mm0               \n\t"   
+               "por %%mm4, %%mm3               \n\t"
+               "psrld $5, %%mm0                \n\t"
+               "pslld $11, %%mm3               \n\t"
+               "por %%mm3, %%mm0               \n\t"
+               MOVNTQ" %%mm0, (%0)             \n\t"
+               "addl $16, %1                   \n\t"
+               "addl $8, %0                    \n\t"
+               "cmpl %2, %1                    \n\t"
+               " jb 1b                         \n\t"
+               : "+r" (d), "+r"(s)
+               : "r" (mm_end), "m" (mask3216g), "m" (mask3216br), "m" (mul3216)
+       );
 #else
-#if 0
-   const uint16_t *s1=( uint16_t * )src;
-   uint16_t *d1=( uint16_t * )dst;
-   uint16_t *e=((uint8_t *)s1)+src_size;
-   while( s1<e ){
-     register int x=*( s1++ );
-     /* rrrrrggggggbbbbb
-        0rrrrrgggggbbbbb
-        0111 1111 1110 0000=0x7FE0
-        00000000000001 1111=0x001F */
-     *( d1++ )=( x&0x001F )|( ( x&0x7FE0 )<<1 );
-   }
+       __asm __volatile(PREFETCH"      %0"::"m"(*src):"memory");
+       __asm __volatile(
+           "movq       %0, %%mm7\n\t"
+           "movq       %1, %%mm6\n\t"
+           ::"m"(red_16mask),"m"(green_16mask));
+       while(s < mm_end)
+       {
+           __asm __volatile(
+               PREFETCH" 32%1\n\t"
+               "movd   %1, %%mm0\n\t"
+               "movd   4%1, %%mm3\n\t"
+               "punpckldq 8%1, %%mm0\n\t"
+               "punpckldq 12%1, %%mm3\n\t"
+               "movq   %%mm0, %%mm1\n\t"
+               "movq   %%mm0, %%mm2\n\t"
+               "movq   %%mm3, %%mm4\n\t"
+               "movq   %%mm3, %%mm5\n\t"
+               "psrlq  $3, %%mm0\n\t"
+               "psrlq  $3, %%mm3\n\t"
+               "pand   %2, %%mm0\n\t"
+               "pand   %2, %%mm3\n\t"
+               "psrlq  $5, %%mm1\n\t"
+               "psrlq  $5, %%mm4\n\t"
+               "pand   %%mm6, %%mm1\n\t"
+               "pand   %%mm6, %%mm4\n\t"
+               "psrlq  $8, %%mm2\n\t"
+               "psrlq  $8, %%mm5\n\t"
+               "pand   %%mm7, %%mm2\n\t"
+               "pand   %%mm7, %%mm5\n\t"
+               "por    %%mm1, %%mm0\n\t"
+               "por    %%mm4, %%mm3\n\t"
+               "por    %%mm2, %%mm0\n\t"
+               "por    %%mm5, %%mm3\n\t"
+               "psllq  $16, %%mm3\n\t"
+               "por    %%mm3, %%mm0\n\t"
+               MOVNTQ" %%mm0, %0\n\t"
+               :"=m"(*d):"m"(*s),"m"(blue_16mask):"memory");
+               d += 4;
+               s += 16;
+       }
+#endif
+       __asm __volatile(SFENCE:::"memory");
+       __asm __volatile(EMMS:::"memory");
+#endif
+       while(s < end)
+       {
+               const int src= *((uint32_t*)s)++;
+               *d++ = ((src&0xFF)>>3) + ((src&0xFC00)>>5) + ((src&0xF80000)>>8);
+//             *d++ = ((src>>3)&0x1F) + ((src>>5)&0x7E0) + ((src>>8)&0xF800);
+       }
+}
+
+static inline void RENAME(rgb32tobgr16)(const uint8_t *src, uint8_t *dst, unsigned int src_size)
+{
+       const uint8_t *s = src;
+       const uint8_t *end;
+#ifdef HAVE_MMX
+       const uint8_t *mm_end;
+#endif
+       uint16_t *d = (uint16_t *)dst;
+       end = s + src_size;
+#ifdef HAVE_MMX
+       __asm __volatile(PREFETCH"      %0"::"m"(*src):"memory");
+       __asm __volatile(
+           "movq       %0, %%mm7\n\t"
+           "movq       %1, %%mm6\n\t"
+           ::"m"(red_16mask),"m"(green_16mask));
+       mm_end = end - 15;
+       while(s < mm_end)
+       {
+           __asm __volatile(
+               PREFETCH" 32%1\n\t"
+               "movd   %1, %%mm0\n\t"
+               "movd   4%1, %%mm3\n\t"
+               "punpckldq 8%1, %%mm0\n\t"
+               "punpckldq 12%1, %%mm3\n\t"
+               "movq   %%mm0, %%mm1\n\t"
+               "movq   %%mm0, %%mm2\n\t"
+               "movq   %%mm3, %%mm4\n\t"
+               "movq   %%mm3, %%mm5\n\t"
+               "psllq  $8, %%mm0\n\t"
+               "psllq  $8, %%mm3\n\t"
+               "pand   %%mm7, %%mm0\n\t"
+               "pand   %%mm7, %%mm3\n\t"
+               "psrlq  $5, %%mm1\n\t"
+               "psrlq  $5, %%mm4\n\t"
+               "pand   %%mm6, %%mm1\n\t"
+               "pand   %%mm6, %%mm4\n\t"
+               "psrlq  $19, %%mm2\n\t"
+               "psrlq  $19, %%mm5\n\t"
+               "pand   %2, %%mm2\n\t"
+               "pand   %2, %%mm5\n\t"
+               "por    %%mm1, %%mm0\n\t"
+               "por    %%mm4, %%mm3\n\t"
+               "por    %%mm2, %%mm0\n\t"
+               "por    %%mm5, %%mm3\n\t"
+               "psllq  $16, %%mm3\n\t"
+               "por    %%mm3, %%mm0\n\t"
+               MOVNTQ" %%mm0, %0\n\t"
+               :"=m"(*d):"m"(*s),"m"(blue_16mask):"memory");
+               d += 4;
+               s += 16;
+       }
+       __asm __volatile(SFENCE:::"memory");
+       __asm __volatile(EMMS:::"memory");
+#endif
+       while(s < end)
+       {
+               const int src= *((uint32_t*)s)++;
+               *d++ = ((src&0xF8)<<8) + ((src&0xFC00)>>5) + ((src&0xF80000)>>19);
+       }
+}
+
+static inline void RENAME(rgb32to15)(const uint8_t *src, uint8_t *dst, unsigned src_size)
+{
+       const uint8_t *s = src;
+       const uint8_t *end;
+#ifdef HAVE_MMX
+       const uint8_t *mm_end;
+#endif
+       uint16_t *d = (uint16_t *)dst;
+       end = s + src_size;
+#ifdef HAVE_MMX
+       mm_end = end - 15;
+#if 1 //is faster only if multiplies are reasonable fast (FIXME figure out on which cpus this is faster, on Athlon its slightly faster)
+       asm volatile(
+               "movq %3, %%mm5                 \n\t"
+               "movq %4, %%mm6                 \n\t"
+               "movq %5, %%mm7                 \n\t"
+               ".balign 16                     \n\t"
+               "1:                             \n\t"
+               PREFETCH" 32(%1)                \n\t"
+               "movd   (%1), %%mm0             \n\t"
+               "movd   4(%1), %%mm3            \n\t"
+               "punpckldq 8(%1), %%mm0         \n\t"
+               "punpckldq 12(%1), %%mm3        \n\t"
+               "movq %%mm0, %%mm1              \n\t"
+               "movq %%mm3, %%mm4              \n\t"
+               "pand %%mm6, %%mm0              \n\t"
+               "pand %%mm6, %%mm3              \n\t"
+               "pmaddwd %%mm7, %%mm0           \n\t"
+               "pmaddwd %%mm7, %%mm3           \n\t"
+               "pand %%mm5, %%mm1              \n\t"
+               "pand %%mm5, %%mm4              \n\t"
+               "por %%mm1, %%mm0               \n\t"   
+               "por %%mm4, %%mm3               \n\t"
+               "psrld $6, %%mm0                \n\t"
+               "pslld $10, %%mm3               \n\t"
+               "por %%mm3, %%mm0               \n\t"
+               MOVNTQ" %%mm0, (%0)             \n\t"
+               "addl $16, %1                   \n\t"
+               "addl $8, %0                    \n\t"
+               "cmpl %2, %1                    \n\t"
+               " jb 1b                         \n\t"
+               : "+r" (d), "+r"(s)
+               : "r" (mm_end), "m" (mask3215g), "m" (mask3216br), "m" (mul3215)
+       );
 #else
-       const uint32_t *s1=( uint32_t * )src;
-       uint32_t *d1=( uint32_t * )dst;
-       int i;
-       int size= src_size>>2;
-       for(i=0; i<size; i++)
+       __asm __volatile(PREFETCH"      %0"::"m"(*src):"memory");
+       __asm __volatile(
+           "movq       %0, %%mm7\n\t"
+           "movq       %1, %%mm6\n\t"
+           ::"m"(red_15mask),"m"(green_15mask));
+       while(s < mm_end)
+       {
+           __asm __volatile(
+               PREFETCH" 32%1\n\t"
+               "movd   %1, %%mm0\n\t"
+               "movd   4%1, %%mm3\n\t"
+               "punpckldq 8%1, %%mm0\n\t"
+               "punpckldq 12%1, %%mm3\n\t"
+               "movq   %%mm0, %%mm1\n\t"
+               "movq   %%mm0, %%mm2\n\t"
+               "movq   %%mm3, %%mm4\n\t"
+               "movq   %%mm3, %%mm5\n\t"
+               "psrlq  $3, %%mm0\n\t"
+               "psrlq  $3, %%mm3\n\t"
+               "pand   %2, %%mm0\n\t"
+               "pand   %2, %%mm3\n\t"
+               "psrlq  $6, %%mm1\n\t"
+               "psrlq  $6, %%mm4\n\t"
+               "pand   %%mm6, %%mm1\n\t"
+               "pand   %%mm6, %%mm4\n\t"
+               "psrlq  $9, %%mm2\n\t"
+               "psrlq  $9, %%mm5\n\t"
+               "pand   %%mm7, %%mm2\n\t"
+               "pand   %%mm7, %%mm5\n\t"
+               "por    %%mm1, %%mm0\n\t"
+               "por    %%mm4, %%mm3\n\t"
+               "por    %%mm2, %%mm0\n\t"
+               "por    %%mm5, %%mm3\n\t"
+               "psllq  $16, %%mm3\n\t"
+               "por    %%mm3, %%mm0\n\t"
+               MOVNTQ" %%mm0, %0\n\t"
+               :"=m"(*d):"m"(*s),"m"(blue_15mask):"memory");
+               d += 4;
+               s += 16;
+       }
+#endif
+       __asm __volatile(SFENCE:::"memory");
+       __asm __volatile(EMMS:::"memory");
+#endif
+       while(s < end)
        {
-               register int x= s1[i];
-//             d1[i] = x + (x&0x7FE07FE0); //faster but need msbit =0 which might not allways be true
-               d1[i] = (x&0x7FFF7FFF) + (x&0x7FE07FE0);
+               const int src= *((uint32_t*)s)++;
+               *d++ = ((src&0xFF)>>3) + ((src&0xF800)>>6) + ((src&0xF80000)>>9);
+       }
+}
+
+static inline void RENAME(rgb32tobgr15)(const uint8_t *src, uint8_t *dst, unsigned src_size)
+{
+       const uint8_t *s = src;
+       const uint8_t *end;
+#ifdef HAVE_MMX
+       const uint8_t *mm_end;
+#endif
+       uint16_t *d = (uint16_t *)dst;
+       end = s + src_size;
+#ifdef HAVE_MMX
+       __asm __volatile(PREFETCH"      %0"::"m"(*src):"memory");
+       __asm __volatile(
+           "movq       %0, %%mm7\n\t"
+           "movq       %1, %%mm6\n\t"
+           ::"m"(red_15mask),"m"(green_15mask));
+       mm_end = end - 15;
+       while(s < mm_end)
+       {
+           __asm __volatile(
+               PREFETCH" 32%1\n\t"
+               "movd   %1, %%mm0\n\t"
+               "movd   4%1, %%mm3\n\t"
+               "punpckldq 8%1, %%mm0\n\t"
+               "punpckldq 12%1, %%mm3\n\t"
+               "movq   %%mm0, %%mm1\n\t"
+               "movq   %%mm0, %%mm2\n\t"
+               "movq   %%mm3, %%mm4\n\t"
+               "movq   %%mm3, %%mm5\n\t"
+               "psllq  $7, %%mm0\n\t"
+               "psllq  $7, %%mm3\n\t"
+               "pand   %%mm7, %%mm0\n\t"
+               "pand   %%mm7, %%mm3\n\t"
+               "psrlq  $6, %%mm1\n\t"
+               "psrlq  $6, %%mm4\n\t"
+               "pand   %%mm6, %%mm1\n\t"
+               "pand   %%mm6, %%mm4\n\t"
+               "psrlq  $19, %%mm2\n\t"
+               "psrlq  $19, %%mm5\n\t"
+               "pand   %2, %%mm2\n\t"
+               "pand   %2, %%mm5\n\t"
+               "por    %%mm1, %%mm0\n\t"
+               "por    %%mm4, %%mm3\n\t"
+               "por    %%mm2, %%mm0\n\t"
+               "por    %%mm5, %%mm3\n\t"
+               "psllq  $16, %%mm3\n\t"
+               "por    %%mm3, %%mm0\n\t"
+               MOVNTQ" %%mm0, %0\n\t"
+               :"=m"(*d):"m"(*s),"m"(blue_15mask):"memory");
+               d += 4;
+               s += 16;
+       }
+       __asm __volatile(SFENCE:::"memory");
+       __asm __volatile(EMMS:::"memory");
+#endif
+       while(s < end)
+       {
+               const int src= *((uint32_t*)s)++;
+               *d++ = ((src&0xF8)<<7) + ((src&0xF800)>>6) + ((src&0xF80000)>>19);
+       }
+}
 
+static inline void RENAME(rgb24to16)(const uint8_t *src, uint8_t *dst, unsigned src_size)
+{
+       const uint8_t *s = src;
+       const uint8_t *end;
+#ifdef HAVE_MMX
+       const uint8_t *mm_end;
+#endif
+       uint16_t *d = (uint16_t *)dst;
+       end = s + src_size;
+#ifdef HAVE_MMX
+       __asm __volatile(PREFETCH"      %0"::"m"(*src):"memory");
+       __asm __volatile(
+           "movq       %0, %%mm7\n\t"
+           "movq       %1, %%mm6\n\t"
+           ::"m"(red_16mask),"m"(green_16mask));
+       mm_end = end - 11;
+       while(s < mm_end)
+       {
+           __asm __volatile(
+               PREFETCH" 32%1\n\t"
+               "movd   %1, %%mm0\n\t"
+               "movd   3%1, %%mm3\n\t"
+               "punpckldq 6%1, %%mm0\n\t"
+               "punpckldq 9%1, %%mm3\n\t"
+               "movq   %%mm0, %%mm1\n\t"
+               "movq   %%mm0, %%mm2\n\t"
+               "movq   %%mm3, %%mm4\n\t"
+               "movq   %%mm3, %%mm5\n\t"
+               "psrlq  $3, %%mm0\n\t"
+               "psrlq  $3, %%mm3\n\t"
+               "pand   %2, %%mm0\n\t"
+               "pand   %2, %%mm3\n\t"
+               "psrlq  $5, %%mm1\n\t"
+               "psrlq  $5, %%mm4\n\t"
+               "pand   %%mm6, %%mm1\n\t"
+               "pand   %%mm6, %%mm4\n\t"
+               "psrlq  $8, %%mm2\n\t"
+               "psrlq  $8, %%mm5\n\t"
+               "pand   %%mm7, %%mm2\n\t"
+               "pand   %%mm7, %%mm5\n\t"
+               "por    %%mm1, %%mm0\n\t"
+               "por    %%mm4, %%mm3\n\t"
+               "por    %%mm2, %%mm0\n\t"
+               "por    %%mm5, %%mm3\n\t"
+               "psllq  $16, %%mm3\n\t"
+               "por    %%mm3, %%mm0\n\t"
+               MOVNTQ" %%mm0, %0\n\t"
+               :"=m"(*d):"m"(*s),"m"(blue_16mask):"memory");
+               d += 4;
+               s += 12;
+       }
+       __asm __volatile(SFENCE:::"memory");
+       __asm __volatile(EMMS:::"memory");
+#endif
+       while(s < end)
+       {
+               const int b= *s++;
+               const int g= *s++;
+               const int r= *s++;
+               *d++ = (b>>3) | ((g&0xFC)<<3) | ((r&0xF8)<<8);
        }
+}
+
+static inline void RENAME(rgb24tobgr16)(const uint8_t *src, uint8_t *dst, unsigned int src_size)
+{
+       const uint8_t *s = src;
+       const uint8_t *end;
+#ifdef HAVE_MMX
+       const uint8_t *mm_end;
 #endif
+       uint16_t *d = (uint16_t *)dst;
+       end = s + src_size;
+#ifdef HAVE_MMX
+       __asm __volatile(PREFETCH"      %0"::"m"(*src):"memory");
+       __asm __volatile(
+           "movq       %0, %%mm7\n\t"
+           "movq       %1, %%mm6\n\t"
+           ::"m"(red_16mask),"m"(green_16mask));
+       mm_end = end - 15;
+       while(s < mm_end)
+       {
+           __asm __volatile(
+               PREFETCH" 32%1\n\t"
+               "movd   %1, %%mm0\n\t"
+               "movd   3%1, %%mm3\n\t"
+               "punpckldq 6%1, %%mm0\n\t"
+               "punpckldq 9%1, %%mm3\n\t"
+               "movq   %%mm0, %%mm1\n\t"
+               "movq   %%mm0, %%mm2\n\t"
+               "movq   %%mm3, %%mm4\n\t"
+               "movq   %%mm3, %%mm5\n\t"
+               "psllq  $8, %%mm0\n\t"
+               "psllq  $8, %%mm3\n\t"
+               "pand   %%mm7, %%mm0\n\t"
+               "pand   %%mm7, %%mm3\n\t"
+               "psrlq  $5, %%mm1\n\t"
+               "psrlq  $5, %%mm4\n\t"
+               "pand   %%mm6, %%mm1\n\t"
+               "pand   %%mm6, %%mm4\n\t"
+               "psrlq  $19, %%mm2\n\t"
+               "psrlq  $19, %%mm5\n\t"
+               "pand   %2, %%mm2\n\t"
+               "pand   %2, %%mm5\n\t"
+               "por    %%mm1, %%mm0\n\t"
+               "por    %%mm4, %%mm3\n\t"
+               "por    %%mm2, %%mm0\n\t"
+               "por    %%mm5, %%mm3\n\t"
+               "psllq  $16, %%mm3\n\t"
+               "por    %%mm3, %%mm0\n\t"
+               MOVNTQ" %%mm0, %0\n\t"
+               :"=m"(*d):"m"(*s),"m"(blue_16mask):"memory");
+               d += 4;
+               s += 12;
+       }
+       __asm __volatile(SFENCE:::"memory");
+       __asm __volatile(EMMS:::"memory");
 #endif
+       while(s < end)
+       {
+               const int r= *s++;
+               const int g= *s++;
+               const int b= *s++;
+               *d++ = (b>>3) | ((g&0xFC)<<3) | ((r&0xF8)<<8);
+       }
 }
 
-/**
- * Pallete is assumed to contain bgr32
- */
-void palette8torgb32(const uint8_t *src, uint8_t *dst, uint32_t num_pixels, const uint8_t *palette)
+static inline void RENAME(rgb24to15)(const uint8_t *src, uint8_t *dst, unsigned src_size)
 {
-       uint32_t i;
-       for(i=0; i<num_pixels; i++)
-               ((uint32_t *)dst)[i] = ((uint32_t *)palette)[ src[i] ];
+       const uint8_t *s = src;
+       const uint8_t *end;
+#ifdef HAVE_MMX
+       const uint8_t *mm_end;
+#endif
+       uint16_t *d = (uint16_t *)dst;
+       end = s + src_size;
+#ifdef HAVE_MMX
+       __asm __volatile(PREFETCH"      %0"::"m"(*src):"memory");
+       __asm __volatile(
+           "movq       %0, %%mm7\n\t"
+           "movq       %1, %%mm6\n\t"
+           ::"m"(red_15mask),"m"(green_15mask));
+       mm_end = end - 11;
+       while(s < mm_end)
+       {
+           __asm __volatile(
+               PREFETCH" 32%1\n\t"
+               "movd   %1, %%mm0\n\t"
+               "movd   3%1, %%mm3\n\t"
+               "punpckldq 6%1, %%mm0\n\t"
+               "punpckldq 9%1, %%mm3\n\t"
+               "movq   %%mm0, %%mm1\n\t"
+               "movq   %%mm0, %%mm2\n\t"
+               "movq   %%mm3, %%mm4\n\t"
+               "movq   %%mm3, %%mm5\n\t"
+               "psrlq  $3, %%mm0\n\t"
+               "psrlq  $3, %%mm3\n\t"
+               "pand   %2, %%mm0\n\t"
+               "pand   %2, %%mm3\n\t"
+               "psrlq  $6, %%mm1\n\t"
+               "psrlq  $6, %%mm4\n\t"
+               "pand   %%mm6, %%mm1\n\t"
+               "pand   %%mm6, %%mm4\n\t"
+               "psrlq  $9, %%mm2\n\t"
+               "psrlq  $9, %%mm5\n\t"
+               "pand   %%mm7, %%mm2\n\t"
+               "pand   %%mm7, %%mm5\n\t"
+               "por    %%mm1, %%mm0\n\t"
+               "por    %%mm4, %%mm3\n\t"
+               "por    %%mm2, %%mm0\n\t"
+               "por    %%mm5, %%mm3\n\t"
+               "psllq  $16, %%mm3\n\t"
+               "por    %%mm3, %%mm0\n\t"
+               MOVNTQ" %%mm0, %0\n\t"
+               :"=m"(*d):"m"(*s),"m"(blue_15mask):"memory");
+               d += 4;
+               s += 12;
+       }
+       __asm __volatile(SFENCE:::"memory");
+       __asm __volatile(EMMS:::"memory");
+#endif
+       while(s < end)
+       {
+               const int b= *s++;
+               const int g= *s++;
+               const int r= *s++;
+               *d++ = (b>>3) | ((g&0xF8)<<2) | ((r&0xF8)<<7);
+       }
 }
 
-/**
- * Pallete is assumed to contain bgr32
- */
-void palette8torgb24(const uint8_t *src, uint8_t *dst, uint32_t num_pixels, const uint8_t *palette)
+static inline void RENAME(rgb24tobgr15)(const uint8_t *src, uint8_t *dst, unsigned src_size)
 {
-       uint32_t i;
+       const uint8_t *s = src;
+       const uint8_t *end;
+#ifdef HAVE_MMX
+       const uint8_t *mm_end;
+#endif
+       uint16_t *d = (uint16_t *)dst;
+       end = s + src_size;
+#ifdef HAVE_MMX
+       __asm __volatile(PREFETCH"      %0"::"m"(*src):"memory");
+       __asm __volatile(
+           "movq       %0, %%mm7\n\t"
+           "movq       %1, %%mm6\n\t"
+           ::"m"(red_15mask),"m"(green_15mask));
+       mm_end = end - 15;
+       while(s < mm_end)
+       {
+           __asm __volatile(
+               PREFETCH" 32%1\n\t"
+               "movd   %1, %%mm0\n\t"
+               "movd   3%1, %%mm3\n\t"
+               "punpckldq 6%1, %%mm0\n\t"
+               "punpckldq 9%1, %%mm3\n\t"
+               "movq   %%mm0, %%mm1\n\t"
+               "movq   %%mm0, %%mm2\n\t"
+               "movq   %%mm3, %%mm4\n\t"
+               "movq   %%mm3, %%mm5\n\t"
+               "psllq  $7, %%mm0\n\t"
+               "psllq  $7, %%mm3\n\t"
+               "pand   %%mm7, %%mm0\n\t"
+               "pand   %%mm7, %%mm3\n\t"
+               "psrlq  $6, %%mm1\n\t"
+               "psrlq  $6, %%mm4\n\t"
+               "pand   %%mm6, %%mm1\n\t"
+               "pand   %%mm6, %%mm4\n\t"
+               "psrlq  $19, %%mm2\n\t"
+               "psrlq  $19, %%mm5\n\t"
+               "pand   %2, %%mm2\n\t"
+               "pand   %2, %%mm5\n\t"
+               "por    %%mm1, %%mm0\n\t"
+               "por    %%mm4, %%mm3\n\t"
+               "por    %%mm2, %%mm0\n\t"
+               "por    %%mm5, %%mm3\n\t"
+               "psllq  $16, %%mm3\n\t"
+               "por    %%mm3, %%mm0\n\t"
+               MOVNTQ" %%mm0, %0\n\t"
+               :"=m"(*d):"m"(*s),"m"(blue_15mask):"memory");
+               d += 4;
+               s += 12;
+       }
+       __asm __volatile(SFENCE:::"memory");
+       __asm __volatile(EMMS:::"memory");
+#endif
+       while(s < end)
+       {
+               const int r= *s++;
+               const int g= *s++;
+               const int b= *s++;
+               *d++ = (b>>3) | ((g&0xF8)<<2) | ((r&0xF8)<<7);
+       }
+}
+
 /*
-       writes 1 byte o much and might cause alignment issues on some architectures?
-       for(i=0; i<num_pixels; i++)
-               ((uint32_t *)(&dst[i*3])) = ((uint32_t *)palette)[ src[i] ];
+  I use here less accurate approximation by simply
+ left-shifting the input
+  value and filling the low order bits with
+ zeroes. This method improves png's
+  compression but this scheme cannot reproduce white exactly, since it does not
+  generate an all-ones maximum value; the net effect is to darken the
+  image slightly.
+
+  The better method should be "left bit replication":
+
+   4 3 2 1 0
+   ---------
+   1 1 0 1 1
+
+   7 6 5 4 3  2 1 0
+   ----------------
+   1 1 0 1 1  1 1 0
+   |=======|  |===|
+       |      Leftmost Bits Repeated to Fill Open Bits
+       |
+   Original Bits
 */
-       for(i=0; i<num_pixels; i++)
+static inline void RENAME(rgb15to24)(const uint8_t *src, uint8_t *dst, unsigned src_size)
+{
+       const uint16_t *end;
+#ifdef HAVE_MMX
+       const uint16_t *mm_end;
+#endif
+       uint8_t *d = (uint8_t *)dst;
+       const uint16_t *s = (uint16_t *)src;
+       end = s + src_size/2;
+#ifdef HAVE_MMX
+       __asm __volatile(PREFETCH"      %0"::"m"(*s):"memory");
+       mm_end = end - 7;
+       while(s < mm_end)
+       {
+           __asm __volatile(
+               PREFETCH" 32%1\n\t"
+               "movq   %1, %%mm0\n\t"
+               "movq   %1, %%mm1\n\t"
+               "movq   %1, %%mm2\n\t"
+               "pand   %2, %%mm0\n\t"
+               "pand   %3, %%mm1\n\t"
+               "pand   %4, %%mm2\n\t"
+               "psllq  $3, %%mm0\n\t"
+               "psrlq  $2, %%mm1\n\t"
+               "psrlq  $7, %%mm2\n\t"
+               "movq   %%mm0, %%mm3\n\t"
+               "movq   %%mm1, %%mm4\n\t"
+               "movq   %%mm2, %%mm5\n\t"
+               "punpcklwd %5, %%mm0\n\t"
+               "punpcklwd %5, %%mm1\n\t"
+               "punpcklwd %5, %%mm2\n\t"
+               "punpckhwd %5, %%mm3\n\t"
+               "punpckhwd %5, %%mm4\n\t"
+               "punpckhwd %5, %%mm5\n\t"
+               "psllq  $8, %%mm1\n\t"
+               "psllq  $16, %%mm2\n\t"
+               "por    %%mm1, %%mm0\n\t"
+               "por    %%mm2, %%mm0\n\t"
+               "psllq  $8, %%mm4\n\t"
+               "psllq  $16, %%mm5\n\t"
+               "por    %%mm4, %%mm3\n\t"
+               "por    %%mm5, %%mm3\n\t"
+
+               "movq   %%mm0, %%mm6\n\t"
+               "movq   %%mm3, %%mm7\n\t"
+               
+               "movq   8%1, %%mm0\n\t"
+               "movq   8%1, %%mm1\n\t"
+               "movq   8%1, %%mm2\n\t"
+               "pand   %2, %%mm0\n\t"
+               "pand   %3, %%mm1\n\t"
+               "pand   %4, %%mm2\n\t"
+               "psllq  $3, %%mm0\n\t"
+               "psrlq  $2, %%mm1\n\t"
+               "psrlq  $7, %%mm2\n\t"
+               "movq   %%mm0, %%mm3\n\t"
+               "movq   %%mm1, %%mm4\n\t"
+               "movq   %%mm2, %%mm5\n\t"
+               "punpcklwd %5, %%mm0\n\t"
+               "punpcklwd %5, %%mm1\n\t"
+               "punpcklwd %5, %%mm2\n\t"
+               "punpckhwd %5, %%mm3\n\t"
+               "punpckhwd %5, %%mm4\n\t"
+               "punpckhwd %5, %%mm5\n\t"
+               "psllq  $8, %%mm1\n\t"
+               "psllq  $16, %%mm2\n\t"
+               "por    %%mm1, %%mm0\n\t"
+               "por    %%mm2, %%mm0\n\t"
+               "psllq  $8, %%mm4\n\t"
+               "psllq  $16, %%mm5\n\t"
+               "por    %%mm4, %%mm3\n\t"
+               "por    %%mm5, %%mm3\n\t"
+
+               :"=m"(*d)
+               :"m"(*s),"m"(mask15b),"m"(mask15g),"m"(mask15r), "m"(mmx_null)
+               :"memory");
+           /* Borrowed 32 to 24 */
+           __asm __volatile(
+               "movq   %%mm0, %%mm4\n\t"
+               "movq   %%mm3, %%mm5\n\t"
+               "movq   %%mm6, %%mm0\n\t"
+               "movq   %%mm7, %%mm1\n\t"
+               
+               "movq   %%mm4, %%mm6\n\t"
+               "movq   %%mm5, %%mm7\n\t"
+               "movq   %%mm0, %%mm2\n\t"
+               "movq   %%mm1, %%mm3\n\t"
+
+               "psrlq  $8, %%mm2\n\t"
+               "psrlq  $8, %%mm3\n\t"
+               "psrlq  $8, %%mm6\n\t"
+               "psrlq  $8, %%mm7\n\t"
+               "pand   %2, %%mm0\n\t"
+               "pand   %2, %%mm1\n\t"
+               "pand   %2, %%mm4\n\t"
+               "pand   %2, %%mm5\n\t"
+               "pand   %3, %%mm2\n\t"
+               "pand   %3, %%mm3\n\t"
+               "pand   %3, %%mm6\n\t"
+               "pand   %3, %%mm7\n\t"
+               "por    %%mm2, %%mm0\n\t"
+               "por    %%mm3, %%mm1\n\t"
+               "por    %%mm6, %%mm4\n\t"
+               "por    %%mm7, %%mm5\n\t"
+
+               "movq   %%mm1, %%mm2\n\t"
+               "movq   %%mm4, %%mm3\n\t"
+               "psllq  $48, %%mm2\n\t"
+               "psllq  $32, %%mm3\n\t"
+               "pand   %4, %%mm2\n\t"
+               "pand   %5, %%mm3\n\t"
+               "por    %%mm2, %%mm0\n\t"
+               "psrlq  $16, %%mm1\n\t"
+               "psrlq  $32, %%mm4\n\t"
+               "psllq  $16, %%mm5\n\t"
+               "por    %%mm3, %%mm1\n\t"
+               "pand   %6, %%mm5\n\t"
+               "por    %%mm5, %%mm4\n\t"
+
+               MOVNTQ" %%mm0, %0\n\t"
+               MOVNTQ" %%mm1, 8%0\n\t"
+               MOVNTQ" %%mm4, 16%0"
+
+               :"=m"(*d)
+               :"m"(*s),"m"(mask24l),"m"(mask24h),"m"(mask24hh),"m"(mask24hhh),"m"(mask24hhhh)
+               :"memory");
+               d += 24;
+               s += 8;
+       }
+       __asm __volatile(SFENCE:::"memory");
+       __asm __volatile(EMMS:::"memory");
+#endif
+       while(s < end)
        {
-               //FIXME slow?
-               dst[0]= palette[ src[i]*4+0 ];
-               dst[1]= palette[ src[i]*4+1 ];
-               dst[2]= palette[ src[i]*4+2 ];
-               dst+= 3;
+               register uint16_t bgr;
+               bgr = *s++;
+               *d++ = (bgr&0x1F)<<3;
+               *d++ = (bgr&0x3E0)>>2;
+               *d++ = (bgr&0x7C00)>>7;
        }
 }
 
-void rgb32to16(const uint8_t *src, uint8_t *dst, uint32_t num_pixels)
+static inline void RENAME(rgb16to24)(const uint8_t *src, uint8_t *dst, unsigned src_size)
 {
-       uint32_t i;
-       for(i=0; i<num_pixels; i+=4)
+       const uint16_t *end;
+#ifdef HAVE_MMX
+       const uint16_t *mm_end;
+#endif
+       uint8_t *d = (uint8_t *)dst;
+       const uint16_t *s = (const uint16_t *)src;
+       end = s + src_size/2;
+#ifdef HAVE_MMX
+       __asm __volatile(PREFETCH"      %0"::"m"(*s):"memory");
+       mm_end = end - 7;
+       while(s < mm_end)
        {
-               const int b= src[i+0];
-               const int g= src[i+1];
-               const int r= src[i+2];
+           __asm __volatile(
+               PREFETCH" 32%1\n\t"
+               "movq   %1, %%mm0\n\t"
+               "movq   %1, %%mm1\n\t"
+               "movq   %1, %%mm2\n\t"
+               "pand   %2, %%mm0\n\t"
+               "pand   %3, %%mm1\n\t"
+               "pand   %4, %%mm2\n\t"
+               "psllq  $3, %%mm0\n\t"
+               "psrlq  $3, %%mm1\n\t"
+               "psrlq  $8, %%mm2\n\t"
+               "movq   %%mm0, %%mm3\n\t"
+               "movq   %%mm1, %%mm4\n\t"
+               "movq   %%mm2, %%mm5\n\t"
+               "punpcklwd %5, %%mm0\n\t"
+               "punpcklwd %5, %%mm1\n\t"
+               "punpcklwd %5, %%mm2\n\t"
+               "punpckhwd %5, %%mm3\n\t"
+               "punpckhwd %5, %%mm4\n\t"
+               "punpckhwd %5, %%mm5\n\t"
+               "psllq  $8, %%mm1\n\t"
+               "psllq  $16, %%mm2\n\t"
+               "por    %%mm1, %%mm0\n\t"
+               "por    %%mm2, %%mm0\n\t"
+               "psllq  $8, %%mm4\n\t"
+               "psllq  $16, %%mm5\n\t"
+               "por    %%mm4, %%mm3\n\t"
+               "por    %%mm5, %%mm3\n\t"
+               
+               "movq   %%mm0, %%mm6\n\t"
+               "movq   %%mm3, %%mm7\n\t"
+
+               "movq   8%1, %%mm0\n\t"
+               "movq   8%1, %%mm1\n\t"
+               "movq   8%1, %%mm2\n\t"
+               "pand   %2, %%mm0\n\t"
+               "pand   %3, %%mm1\n\t"
+               "pand   %4, %%mm2\n\t"
+               "psllq  $3, %%mm0\n\t"
+               "psrlq  $3, %%mm1\n\t"
+               "psrlq  $8, %%mm2\n\t"
+               "movq   %%mm0, %%mm3\n\t"
+               "movq   %%mm1, %%mm4\n\t"
+               "movq   %%mm2, %%mm5\n\t"
+               "punpcklwd %5, %%mm0\n\t"
+               "punpcklwd %5, %%mm1\n\t"
+               "punpcklwd %5, %%mm2\n\t"
+               "punpckhwd %5, %%mm3\n\t"
+               "punpckhwd %5, %%mm4\n\t"
+               "punpckhwd %5, %%mm5\n\t"
+               "psllq  $8, %%mm1\n\t"
+               "psllq  $16, %%mm2\n\t"
+               "por    %%mm1, %%mm0\n\t"
+               "por    %%mm2, %%mm0\n\t"
+               "psllq  $8, %%mm4\n\t"
+               "psllq  $16, %%mm5\n\t"
+               "por    %%mm4, %%mm3\n\t"
+               "por    %%mm5, %%mm3\n\t"
+               :"=m"(*d)
+               :"m"(*s),"m"(mask16b),"m"(mask16g),"m"(mask16r),"m"(mmx_null)           
+               :"memory");
+           /* Borrowed 32 to 24 */
+           __asm __volatile(
+               "movq   %%mm0, %%mm4\n\t"
+               "movq   %%mm3, %%mm5\n\t"
+               "movq   %%mm6, %%mm0\n\t"
+               "movq   %%mm7, %%mm1\n\t"
+               
+               "movq   %%mm4, %%mm6\n\t"
+               "movq   %%mm5, %%mm7\n\t"
+               "movq   %%mm0, %%mm2\n\t"
+               "movq   %%mm1, %%mm3\n\t"
+
+               "psrlq  $8, %%mm2\n\t"
+               "psrlq  $8, %%mm3\n\t"
+               "psrlq  $8, %%mm6\n\t"
+               "psrlq  $8, %%mm7\n\t"
+               "pand   %2, %%mm0\n\t"
+               "pand   %2, %%mm1\n\t"
+               "pand   %2, %%mm4\n\t"
+               "pand   %2, %%mm5\n\t"
+               "pand   %3, %%mm2\n\t"
+               "pand   %3, %%mm3\n\t"
+               "pand   %3, %%mm6\n\t"
+               "pand   %3, %%mm7\n\t"
+               "por    %%mm2, %%mm0\n\t"
+               "por    %%mm3, %%mm1\n\t"
+               "por    %%mm6, %%mm4\n\t"
+               "por    %%mm7, %%mm5\n\t"
+
+               "movq   %%mm1, %%mm2\n\t"
+               "movq   %%mm4, %%mm3\n\t"
+               "psllq  $48, %%mm2\n\t"
+               "psllq  $32, %%mm3\n\t"
+               "pand   %4, %%mm2\n\t"
+               "pand   %5, %%mm3\n\t"
+               "por    %%mm2, %%mm0\n\t"
+               "psrlq  $16, %%mm1\n\t"
+               "psrlq  $32, %%mm4\n\t"
+               "psllq  $16, %%mm5\n\t"
+               "por    %%mm3, %%mm1\n\t"
+               "pand   %6, %%mm5\n\t"
+               "por    %%mm5, %%mm4\n\t"
+
+               MOVNTQ" %%mm0, %0\n\t"
+               MOVNTQ" %%mm1, 8%0\n\t"
+               MOVNTQ" %%mm4, 16%0"
 
-               ((uint16_t *)dst)[i]= (b>>3) | ((g&0xFC)<<3) | ((r&0xF8)<<8);
+               :"=m"(*d)
+               :"m"(*s),"m"(mask24l),"m"(mask24h),"m"(mask24hh),"m"(mask24hhh),"m"(mask24hhhh)
+               :"memory");
+               d += 24;
+               s += 8;
+       }
+       __asm __volatile(SFENCE:::"memory");
+       __asm __volatile(EMMS:::"memory");
+#endif
+       while(s < end)
+       {
+               register uint16_t bgr;
+               bgr = *s++;
+               *d++ = (bgr&0x1F)<<3;
+               *d++ = (bgr&0x7E0)>>3;
+               *d++ = (bgr&0xF800)>>8;
+       }
+}
+
+static inline void RENAME(rgb15to32)(const uint8_t *src, uint8_t *dst, unsigned src_size)
+{
+       const uint16_t *end;
+#ifdef HAVE_MMX
+       const uint16_t *mm_end;
+#endif
+       uint8_t *d = (uint8_t *)dst;
+       const uint16_t *s = (const uint16_t *)src;
+       end = s + src_size/2;
+#ifdef HAVE_MMX
+       __asm __volatile(PREFETCH"      %0"::"m"(*s):"memory");
+       __asm __volatile("pxor  %%mm7,%%mm7\n\t":::"memory");
+       mm_end = end - 3;
+       while(s < mm_end)
+       {
+           __asm __volatile(
+               PREFETCH" 32%1\n\t"
+               "movq   %1, %%mm0\n\t"
+               "movq   %1, %%mm1\n\t"
+               "movq   %1, %%mm2\n\t"
+               "pand   %2, %%mm0\n\t"
+               "pand   %3, %%mm1\n\t"
+               "pand   %4, %%mm2\n\t"
+               "psllq  $3, %%mm0\n\t"
+               "psrlq  $2, %%mm1\n\t"
+               "psrlq  $7, %%mm2\n\t"
+               "movq   %%mm0, %%mm3\n\t"
+               "movq   %%mm1, %%mm4\n\t"
+               "movq   %%mm2, %%mm5\n\t"
+               "punpcklwd %%mm7, %%mm0\n\t"
+               "punpcklwd %%mm7, %%mm1\n\t"
+               "punpcklwd %%mm7, %%mm2\n\t"
+               "punpckhwd %%mm7, %%mm3\n\t"
+               "punpckhwd %%mm7, %%mm4\n\t"
+               "punpckhwd %%mm7, %%mm5\n\t"
+               "psllq  $8, %%mm1\n\t"
+               "psllq  $16, %%mm2\n\t"
+               "por    %%mm1, %%mm0\n\t"
+               "por    %%mm2, %%mm0\n\t"
+               "psllq  $8, %%mm4\n\t"
+               "psllq  $16, %%mm5\n\t"
+               "por    %%mm4, %%mm3\n\t"
+               "por    %%mm5, %%mm3\n\t"
+               MOVNTQ" %%mm0, %0\n\t"
+               MOVNTQ" %%mm3, 8%0\n\t"
+               :"=m"(*d)
+               :"m"(*s),"m"(mask15b),"m"(mask15g),"m"(mask15r)
+               :"memory");
+               d += 16;
+               s += 4;
+       }
+       __asm __volatile(SFENCE:::"memory");
+       __asm __volatile(EMMS:::"memory");
+#endif
+       while(s < end)
+       {
+#if 0 //slightly slower on athlon
+               int bgr= *s++;
+               *((uint32_t*)d)++ = ((bgr&0x1F)<<3) + ((bgr&0x3E0)<<6) + ((bgr&0x7C00)<<9);
+#else
+//FIXME this is very likely wrong for bigendian (and the following converters too)
+               register uint16_t bgr;
+               bgr = *s++;
+               *d++ = (bgr&0x1F)<<3;
+               *d++ = (bgr&0x3E0)>>2;
+               *d++ = (bgr&0x7C00)>>7;
+               *d++ = 0;
+#endif
+       }
+}
+
+static inline void RENAME(rgb16to32)(const uint8_t *src, uint8_t *dst, unsigned src_size)
+{
+       const uint16_t *end;
+#ifdef HAVE_MMX
+       const uint16_t *mm_end;
+#endif
+       uint8_t *d = (uint8_t *)dst;
+       const uint16_t *s = (uint16_t *)src;
+       end = s + src_size/2;
+#ifdef HAVE_MMX
+       __asm __volatile(PREFETCH"      %0"::"m"(*s):"memory");
+       __asm __volatile("pxor  %%mm7,%%mm7\n\t":::"memory");
+       mm_end = end - 3;
+       while(s < mm_end)
+       {
+           __asm __volatile(
+               PREFETCH" 32%1\n\t"
+               "movq   %1, %%mm0\n\t"
+               "movq   %1, %%mm1\n\t"
+               "movq   %1, %%mm2\n\t"
+               "pand   %2, %%mm0\n\t"
+               "pand   %3, %%mm1\n\t"
+               "pand   %4, %%mm2\n\t"
+               "psllq  $3, %%mm0\n\t"
+               "psrlq  $3, %%mm1\n\t"
+               "psrlq  $8, %%mm2\n\t"
+               "movq   %%mm0, %%mm3\n\t"
+               "movq   %%mm1, %%mm4\n\t"
+               "movq   %%mm2, %%mm5\n\t"
+               "punpcklwd %%mm7, %%mm0\n\t"
+               "punpcklwd %%mm7, %%mm1\n\t"
+               "punpcklwd %%mm7, %%mm2\n\t"
+               "punpckhwd %%mm7, %%mm3\n\t"
+               "punpckhwd %%mm7, %%mm4\n\t"
+               "punpckhwd %%mm7, %%mm5\n\t"
+               "psllq  $8, %%mm1\n\t"
+               "psllq  $16, %%mm2\n\t"
+               "por    %%mm1, %%mm0\n\t"
+               "por    %%mm2, %%mm0\n\t"
+               "psllq  $8, %%mm4\n\t"
+               "psllq  $16, %%mm5\n\t"
+               "por    %%mm4, %%mm3\n\t"
+               "por    %%mm5, %%mm3\n\t"
+               MOVNTQ" %%mm0, %0\n\t"
+               MOVNTQ" %%mm3, 8%0\n\t"
+               :"=m"(*d)
+               :"m"(*s),"m"(mask16b),"m"(mask16g),"m"(mask16r)
+               :"memory");
+               d += 16;
+               s += 4;
+       }
+       __asm __volatile(SFENCE:::"memory");
+       __asm __volatile(EMMS:::"memory");
+#endif
+       while(s < end)
+       {
+               register uint16_t bgr;
+               bgr = *s++;
+               *d++ = (bgr&0x1F)<<3;
+               *d++ = (bgr&0x7E0)>>3;
+               *d++ = (bgr&0xF800)>>8;
+               *d++ = 0;
        }
 }
 
-void rgb32to15(const uint8_t *src, uint8_t *dst, uint32_t num_pixels)
+static inline void RENAME(rgb32tobgr32)(const uint8_t *src, uint8_t *dst, unsigned int src_size)
 {
-       uint32_t i;
-       for(i=0; i<num_pixels; i+=4)
+#ifdef HAVE_MMX
+/* TODO: unroll this loop */
+       asm volatile (
+               "xorl %%eax, %%eax              \n\t"
+               ".balign 16                     \n\t"
+               "1:                             \n\t"
+               PREFETCH" 32(%0, %%eax)         \n\t"
+               "movq (%0, %%eax), %%mm0        \n\t"
+               "movq %%mm0, %%mm1              \n\t"
+               "movq %%mm0, %%mm2              \n\t"
+               "pslld $16, %%mm0               \n\t"
+               "psrld $16, %%mm1               \n\t"
+               "pand "MANGLE(mask32r)", %%mm0  \n\t"
+               "pand "MANGLE(mask32g)", %%mm2  \n\t"
+               "pand "MANGLE(mask32b)", %%mm1  \n\t"
+               "por %%mm0, %%mm2               \n\t"
+               "por %%mm1, %%mm2               \n\t"
+               MOVNTQ" %%mm2, (%1, %%eax)      \n\t"
+               "addl $8, %%eax                 \n\t"
+               "cmpl %2, %%eax                 \n\t"
+               " jb 1b                         \n\t"
+               :: "r" (src), "r"(dst), "r" (src_size-7)
+               : "%eax"
+       );
+
+       __asm __volatile(SFENCE:::"memory");
+       __asm __volatile(EMMS:::"memory");
+#else
+       unsigned i;
+       unsigned num_pixels = src_size >> 2;
+       for(i=0; i<num_pixels; i++)
        {
-               const int b= src[i+0];
-               const int g= src[i+1];
-               const int r= src[i+2];
+#ifdef WORDS_BIGENDIAN  
+         dst[4*i + 1] = src[4*i + 3];
+         dst[4*i + 2] = src[4*i + 2];
+         dst[4*i + 3] = src[4*i + 1];
+#else
+         dst[4*i + 0] = src[4*i + 2];
+         dst[4*i + 1] = src[4*i + 1];
+         dst[4*i + 2] = src[4*i + 0];
+#endif
+       }
+#endif
+}
+
+static inline void RENAME(rgb24tobgr24)(const uint8_t *src, uint8_t *dst, unsigned int src_size)
+{
+       unsigned i;
+#ifdef HAVE_MMX
+       int mmx_size= 23 - src_size;
+       asm volatile (
+               "movq "MANGLE(mask24r)", %%mm5  \n\t"
+               "movq "MANGLE(mask24g)", %%mm6  \n\t"
+               "movq "MANGLE(mask24b)", %%mm7  \n\t"
+               ".balign 16                     \n\t"
+               "1:                             \n\t"
+               PREFETCH" 32(%1, %%eax)         \n\t"
+               "movq   (%1, %%eax), %%mm0      \n\t" // BGR BGR BG
+               "movq   (%1, %%eax), %%mm1      \n\t" // BGR BGR BG
+               "movq  2(%1, %%eax), %%mm2      \n\t" // R BGR BGR B
+               "psllq $16, %%mm0               \n\t" // 00 BGR BGR
+               "pand %%mm5, %%mm0              \n\t"
+               "pand %%mm6, %%mm1              \n\t"
+               "pand %%mm7, %%mm2              \n\t"
+               "por %%mm0, %%mm1               \n\t"
+               "por %%mm2, %%mm1               \n\t"                
+               "movq  6(%1, %%eax), %%mm0      \n\t" // BGR BGR BG
+               MOVNTQ" %%mm1,   (%2, %%eax)    \n\t" // RGB RGB RG
+               "movq  8(%1, %%eax), %%mm1      \n\t" // R BGR BGR B
+               "movq 10(%1, %%eax), %%mm2      \n\t" // GR BGR BGR
+               "pand %%mm7, %%mm0              \n\t"
+               "pand %%mm5, %%mm1              \n\t"
+               "pand %%mm6, %%mm2              \n\t"
+               "por %%mm0, %%mm1               \n\t"
+               "por %%mm2, %%mm1               \n\t"                
+               "movq 14(%1, %%eax), %%mm0      \n\t" // R BGR BGR B
+               MOVNTQ" %%mm1,  8(%2, %%eax)    \n\t" // B RGB RGB R
+               "movq 16(%1, %%eax), %%mm1      \n\t" // GR BGR BGR
+               "movq 18(%1, %%eax), %%mm2      \n\t" // BGR BGR BG
+               "pand %%mm6, %%mm0              \n\t"
+               "pand %%mm7, %%mm1              \n\t"
+               "pand %%mm5, %%mm2              \n\t"
+               "por %%mm0, %%mm1               \n\t"
+               "por %%mm2, %%mm1               \n\t"                
+               MOVNTQ" %%mm1, 16(%2, %%eax)    \n\t"
+               "addl $24, %%eax                \n\t"
+               " js 1b                         \n\t"
+               : "+a" (mmx_size)
+               : "r" (src-mmx_size), "r"(dst-mmx_size)
+       );
+
+       __asm __volatile(SFENCE:::"memory");
+       __asm __volatile(EMMS:::"memory");
+
+       if(mmx_size==23) return; //finihsed, was multiple of 8
 
-               ((uint16_t *)dst)[i]= (b>>3) | ((g&0xF8)<<3) | ((r&0xF8)<<7);
+       src+= src_size;
+       dst+= src_size;
+       src_size= 23-mmx_size;
+       src-= src_size;
+       dst-= src_size;
+#endif
+       for(i=0; i<src_size; i+=3)
+       {
+               register uint8_t x;
+               x          = src[i + 2];
+               dst[i + 1] = src[i + 1];
+               dst[i + 2] = src[i + 0];
+               dst[i + 0] = x;
        }
 }
 
+static inline void RENAME(yuvPlanartoyuy2)(const uint8_t *ysrc, const uint8_t *usrc, const uint8_t *vsrc, uint8_t *dst,
+       unsigned int width, unsigned int height,
+       int lumStride, int chromStride, int dstStride, int vertLumPerChroma)
+{
+       unsigned y;
+       const unsigned chromWidth= width>>1;
+       for(y=0; y<height; y++)
+       {
+#ifdef HAVE_MMX
+//FIXME handle 2 lines a once (fewer prefetch, reuse some chrom, but very likely limited by mem anyway)
+               asm volatile(
+                       "xorl %%eax, %%eax              \n\t"
+                       ".balign 16                     \n\t"
+                       "1:                             \n\t"
+                       PREFETCH" 32(%1, %%eax, 2)      \n\t"
+                       PREFETCH" 32(%2, %%eax)         \n\t"
+                       PREFETCH" 32(%3, %%eax)         \n\t"
+                       "movq (%2, %%eax), %%mm0        \n\t" // U(0)
+                       "movq %%mm0, %%mm2              \n\t" // U(0)
+                       "movq (%3, %%eax), %%mm1        \n\t" // V(0)
+                       "punpcklbw %%mm1, %%mm0         \n\t" // UVUV UVUV(0)
+                       "punpckhbw %%mm1, %%mm2         \n\t" // UVUV UVUV(8)
+
+                       "movq (%1, %%eax,2), %%mm3      \n\t" // Y(0)
+                       "movq 8(%1, %%eax,2), %%mm5     \n\t" // Y(8)
+                       "movq %%mm3, %%mm4              \n\t" // Y(0)
+                       "movq %%mm5, %%mm6              \n\t" // Y(8)
+                       "punpcklbw %%mm0, %%mm3         \n\t" // YUYV YUYV(0)
+                       "punpckhbw %%mm0, %%mm4         \n\t" // YUYV YUYV(4)
+                       "punpcklbw %%mm2, %%mm5         \n\t" // YUYV YUYV(8)
+                       "punpckhbw %%mm2, %%mm6         \n\t" // YUYV YUYV(12)
+
+                       MOVNTQ" %%mm3, (%0, %%eax, 4)   \n\t"
+                       MOVNTQ" %%mm4, 8(%0, %%eax, 4)  \n\t"
+                       MOVNTQ" %%mm5, 16(%0, %%eax, 4) \n\t"
+                       MOVNTQ" %%mm6, 24(%0, %%eax, 4) \n\t"
+
+                       "addl $8, %%eax                 \n\t"
+                       "cmpl %4, %%eax                 \n\t"
+                       " jb 1b                         \n\t"
+                       ::"r"(dst), "r"(ysrc), "r"(usrc), "r"(vsrc), "g" (chromWidth)
+                       : "%eax"
+               );
+#else
+
+#if defined ARCH_ALPHA && defined HAVE_MVI
+#define pl2yuy2(n)                                     \
+       y1 = yc[n];                                     \
+       y2 = yc2[n];                                    \
+       u = uc[n];                                      \
+       v = vc[n];                                      \
+       asm("unpkbw %1, %0" : "=r"(y1) : "r"(y1));      \
+       asm("unpkbw %1, %0" : "=r"(y2) : "r"(y2));      \
+       asm("unpkbl %1, %0" : "=r"(u) : "r"(u));        \
+       asm("unpkbl %1, %0" : "=r"(v) : "r"(v));        \
+       yuv1 = (u << 8) + (v << 24);                    \
+       yuv2 = yuv1 + y2;                               \
+       yuv1 += y1;                                     \
+       qdst[n] = yuv1;                                 \
+       qdst2[n] = yuv2;
+
+               int i;
+               uint64_t *qdst = (uint64_t *) dst;
+               uint64_t *qdst2 = (uint64_t *) (dst + dstStride);
+               const uint32_t *yc = (uint32_t *) ysrc;
+               const uint32_t *yc2 = (uint32_t *) (ysrc + lumStride);
+               const uint16_t *uc = (uint16_t*) usrc, *vc = (uint16_t*) vsrc;
+               for(i = 0; i < chromWidth; i += 8){
+                       uint64_t y1, y2, yuv1, yuv2;
+                       uint64_t u, v;
+                       /* Prefetch */
+                       asm("ldq $31,64(%0)" :: "r"(yc));
+                       asm("ldq $31,64(%0)" :: "r"(yc2));
+                       asm("ldq $31,64(%0)" :: "r"(uc));
+                       asm("ldq $31,64(%0)" :: "r"(vc));
+
+                       pl2yuy2(0);
+                       pl2yuy2(1);
+                       pl2yuy2(2);
+                       pl2yuy2(3);
+
+                       yc += 4;
+                       yc2 += 4;
+                       uc += 4;
+                       vc += 4;
+                       qdst += 4;
+                       qdst2 += 4;
+               }
+               y++;
+               ysrc += lumStride;
+               dst += dstStride;
+
+#elif __WORDSIZE >= 64
+               int i;
+               uint64_t *ldst = (uint64_t *) dst;
+               const uint8_t *yc = ysrc, *uc = usrc, *vc = vsrc;
+               for(i = 0; i < chromWidth; i += 2){
+                       uint64_t k, l;
+                       k = yc[0] + (uc[0] << 8) +
+                           (yc[1] << 16) + (vc[0] << 24);
+                       l = yc[2] + (uc[1] << 8) +
+                           (yc[3] << 16) + (vc[1] << 24);
+                       *ldst++ = k + (l << 32);
+                       yc += 4;
+                       uc += 2;
+                       vc += 2;
+               }
+
+#else
+               int i, *idst = (int32_t *) dst;
+               const uint8_t *yc = ysrc, *uc = usrc, *vc = vsrc;
+               for(i = 0; i < chromWidth; i++){
+                       *idst++ = yc[0] + (uc[0] << 8) +
+                           (yc[1] << 16) + (vc[0] << 24);
+                       yc += 2;
+                       uc++;
+                       vc++;
+               }
+#endif
+#endif
+               if((y&(vertLumPerChroma-1))==(vertLumPerChroma-1) )
+               {
+                       usrc += chromStride;
+                       vsrc += chromStride;
+               }
+               ysrc += lumStride;
+               dst += dstStride;
+       }
+#ifdef HAVE_MMX
+asm(    EMMS" \n\t"
+        SFENCE" \n\t"
+        :::"memory");
+#endif
+}
 
 /**
- * Palette is assumed to contain bgr16, see rgb32to16 to convert the palette
+ *
+ * height should be a multiple of 2 and width should be a multiple of 16 (if this is a
+ * problem for anyone then tell me, and ill fix it)
  */
-void palette8torgb16(const uint8_t *src, uint8_t *dst, uint32_t num_pixels, const uint8_t *palette)
+static inline void RENAME(yv12toyuy2)(const uint8_t *ysrc, const uint8_t *usrc, const uint8_t *vsrc, uint8_t *dst,
+       unsigned int width, unsigned int height,
+       int lumStride, int chromStride, int dstStride)
 {
-       uint32_t i;
-       for(i=0; i<num_pixels; i++)
-               ((uint16_t *)dst)[i] = ((uint16_t *)palette)[ src[i] ];
+       //FIXME interpolate chroma
+       RENAME(yuvPlanartoyuy2)(ysrc, usrc, vsrc, dst, width, height, lumStride, chromStride, dstStride, 2);
 }
 
 /**
- * Pallete is assumed to contain bgr15, see rgb32to15 to convert the palette
+ *
+ * width should be a multiple of 16
  */
-void palette8torgb15(const uint8_t *src, uint8_t *dst, uint32_t num_pixels, const uint8_t *palette)
+static inline void RENAME(yuv422ptoyuy2)(const uint8_t *ysrc, const uint8_t *usrc, const uint8_t *vsrc, uint8_t *dst,
+       unsigned int width, unsigned int height,
+       int lumStride, int chromStride, int dstStride)
 {
-       uint32_t i;
-       for(i=0; i<num_pixels; i++)
-               ((uint16_t *)dst)[i] = ((uint16_t *)palette)[ src[i] ];
+       RENAME(yuvPlanartoyuy2)(ysrc, usrc, vsrc, dst, width, height, lumStride, chromStride, dstStride, 1);
 }
+
 /**
  *
- * num_pixels must be a multiple of 16 for the MMX version
+ * height should be a multiple of 2 and width should be a multiple of 16 (if this is a
+ * problem for anyone then tell me, and ill fix it)
  */
-void yv12toyuy2(const uint8_t *ysrc, const uint8_t *usrc, const uint8_t *vsrc, uint8_t *dst, uint32_t num_pixels)
+static inline void RENAME(yuy2toyv12)(const uint8_t *src, uint8_t *ydst, uint8_t *udst, uint8_t *vdst,
+       unsigned int width, unsigned int height,
+       int lumStride, int chromStride, int srcStride)
 {
+       unsigned y;
+       const unsigned chromWidth= width>>1;
+       for(y=0; y<height; y+=2)
+       {
 #ifdef HAVE_MMX
-       asm volatile(
-               "xorl %%eax, %%eax              \n\t"
-               "1:                             \n\t"
-               PREFETCH" 32(%1, %%eax, 2)      \n\t"
-               PREFETCH" 32(%2, %%eax)         \n\t"
-               PREFETCH" 32(%3, %%eax)         \n\t"
-               "movq (%2, %%eax), %%mm0        \n\t" // U(0)
-               "movq %%mm0, %%mm2              \n\t" // U(0)
-               "movq (%3, %%eax), %%mm1        \n\t" // V(0)
-               "punpcklbw %%mm1, %%mm0         \n\t" // UVUV UVUV(0)
-               "punpckhbw %%mm1, %%mm2         \n\t" // UVUV UVUV(8)
-
-               "movq (%1, %%eax,2), %%mm3      \n\t" // Y(0)
-               "movq 8(%1, %%eax,2), %%mm5     \n\t" // Y(8)
-               "movq %%mm3, %%mm4              \n\t" // Y(0)
-               "movq %%mm5, %%mm6              \n\t" // Y(8)
-               "punpcklbw %%mm0, %%mm3         \n\t" // YUYV YUYV(0)
-               "punpckhbw %%mm0, %%mm4         \n\t" // YUYV YUYV(4)
-               "punpcklbw %%mm2, %%mm5         \n\t" // YUYV YUYV(8)
-               "punpckhbw %%mm2, %%mm6         \n\t" // YUYV YUYV(12)
-
-               MOVNTQ" %%mm3, (%0, %%eax, 4)   \n\t"
-               MOVNTQ" %%mm4, 8(%0, %%eax, 4)  \n\t"
-               MOVNTQ" %%mm5, 16(%0, %%eax, 4) \n\t"
-               MOVNTQ" %%mm6, 24(%0, %%eax, 4) \n\t"
+               asm volatile(
+                       "xorl %%eax, %%eax              \n\t"
+                       "pcmpeqw %%mm7, %%mm7           \n\t"
+                       "psrlw $8, %%mm7                \n\t" // FF,00,FF,00...
+                       ".balign 16                     \n\t"
+                       "1:                             \n\t"
+                       PREFETCH" 64(%0, %%eax, 4)      \n\t"
+                       "movq (%0, %%eax, 4), %%mm0     \n\t" // YUYV YUYV(0)
+                       "movq 8(%0, %%eax, 4), %%mm1    \n\t" // YUYV YUYV(4)
+                       "movq %%mm0, %%mm2              \n\t" // YUYV YUYV(0)
+                       "movq %%mm1, %%mm3              \n\t" // YUYV YUYV(4)
+                       "psrlw $8, %%mm0                \n\t" // U0V0 U0V0(0)
+                       "psrlw $8, %%mm1                \n\t" // U0V0 U0V0(4)
+                       "pand %%mm7, %%mm2              \n\t" // Y0Y0 Y0Y0(0)
+                       "pand %%mm7, %%mm3              \n\t" // Y0Y0 Y0Y0(4)
+                       "packuswb %%mm1, %%mm0          \n\t" // UVUV UVUV(0)
+                       "packuswb %%mm3, %%mm2          \n\t" // YYYY YYYY(0)
 
-               "addl $8, %%eax                 \n\t"
-               "cmpl %4, %%eax                 \n\t"
-               " jb 1b                         \n\t"
-               EMMS" \n\t"
-               SFENCE
-               ::"r"(dst), "r"(ysrc), "r"(usrc), "r"(vsrc), "r" (num_pixels>>1)
-               : "memory", "%eax"
-       );
+                       MOVNTQ" %%mm2, (%1, %%eax, 2)   \n\t"
+
+                       "movq 16(%0, %%eax, 4), %%mm1   \n\t" // YUYV YUYV(8)
+                       "movq 24(%0, %%eax, 4), %%mm2   \n\t" // YUYV YUYV(12)
+                       "movq %%mm1, %%mm3              \n\t" // YUYV YUYV(8)
+                       "movq %%mm2, %%mm4              \n\t" // YUYV YUYV(12)
+                       "psrlw $8, %%mm1                \n\t" // U0V0 U0V0(8)
+                       "psrlw $8, %%mm2                \n\t" // U0V0 U0V0(12)
+                       "pand %%mm7, %%mm3              \n\t" // Y0Y0 Y0Y0(8)
+                       "pand %%mm7, %%mm4              \n\t" // Y0Y0 Y0Y0(12)
+                       "packuswb %%mm2, %%mm1          \n\t" // UVUV UVUV(8)
+                       "packuswb %%mm4, %%mm3          \n\t" // YYYY YYYY(8)
+
+                       MOVNTQ" %%mm3, 8(%1, %%eax, 2)  \n\t"
+
+                       "movq %%mm0, %%mm2              \n\t" // UVUV UVUV(0)
+                       "movq %%mm1, %%mm3              \n\t" // UVUV UVUV(8)
+                       "psrlw $8, %%mm0                \n\t" // V0V0 V0V0(0)
+                       "psrlw $8, %%mm1                \n\t" // V0V0 V0V0(8)
+                       "pand %%mm7, %%mm2              \n\t" // U0U0 U0U0(0)
+                       "pand %%mm7, %%mm3              \n\t" // U0U0 U0U0(8)
+                       "packuswb %%mm1, %%mm0          \n\t" // VVVV VVVV(0)
+                       "packuswb %%mm3, %%mm2          \n\t" // UUUU UUUU(0)
+
+                       MOVNTQ" %%mm0, (%3, %%eax)      \n\t"
+                       MOVNTQ" %%mm2, (%2, %%eax)      \n\t"
+
+                       "addl $8, %%eax                 \n\t"
+                       "cmpl %4, %%eax                 \n\t"
+                       " jb 1b                         \n\t"
+                       ::"r"(src), "r"(ydst), "r"(udst), "r"(vdst), "g" (chromWidth)
+                       : "memory", "%eax"
+               );
+
+               ydst += lumStride;
+               src  += srcStride;
+
+               asm volatile(
+                       "xorl %%eax, %%eax              \n\t"
+                       ".balign 16                     \n\t"
+                       "1:                             \n\t"
+                       PREFETCH" 64(%0, %%eax, 4)      \n\t"
+                       "movq (%0, %%eax, 4), %%mm0     \n\t" // YUYV YUYV(0)
+                       "movq 8(%0, %%eax, 4), %%mm1    \n\t" // YUYV YUYV(4)
+                       "movq 16(%0, %%eax, 4), %%mm2   \n\t" // YUYV YUYV(8)
+                       "movq 24(%0, %%eax, 4), %%mm3   \n\t" // YUYV YUYV(12)
+                       "pand %%mm7, %%mm0              \n\t" // Y0Y0 Y0Y0(0)
+                       "pand %%mm7, %%mm1              \n\t" // Y0Y0 Y0Y0(4)
+                       "pand %%mm7, %%mm2              \n\t" // Y0Y0 Y0Y0(8)
+                       "pand %%mm7, %%mm3              \n\t" // Y0Y0 Y0Y0(12)
+                       "packuswb %%mm1, %%mm0          \n\t" // YYYY YYYY(0)
+                       "packuswb %%mm3, %%mm2          \n\t" // YYYY YYYY(8)
 
+                       MOVNTQ" %%mm0, (%1, %%eax, 2)   \n\t"
+                       MOVNTQ" %%mm2, 8(%1, %%eax, 2)  \n\t"
+
+                       "addl $8, %%eax                 \n\t"
+                       "cmpl %4, %%eax                 \n\t"
+                       " jb 1b                         \n\t"
+
+                       ::"r"(src), "r"(ydst), "r"(udst), "r"(vdst), "g" (chromWidth)
+                       : "memory", "%eax"
+               );
 #else
-       int i;
-       num_pixels>>=1;
-       for(i=0; i<num_pixels; i++)
+               unsigned i;
+               for(i=0; i<chromWidth; i++)
+               {
+                       ydst[2*i+0]     = src[4*i+0];
+                       udst[i]         = src[4*i+1];
+                       ydst[2*i+1]     = src[4*i+2];
+                       vdst[i]         = src[4*i+3];
+               }
+               ydst += lumStride;
+               src  += srcStride;
+
+               for(i=0; i<chromWidth; i++)
+               {
+                       ydst[2*i+0]     = src[4*i+0];
+                       ydst[2*i+1]     = src[4*i+2];
+               }
+#endif
+               udst += chromStride;
+               vdst += chromStride;
+               ydst += lumStride;
+               src  += srcStride;
+       }
+#ifdef HAVE_MMX
+asm volatile(   EMMS" \n\t"
+               SFENCE" \n\t"
+               :::"memory");
+#endif
+}
+
+static inline void RENAME(yvu9toyv12)(const uint8_t *ysrc, const uint8_t *usrc, const uint8_t *vsrc,
+       uint8_t *ydst, uint8_t *udst, uint8_t *vdst,
+       unsigned int width, unsigned int height, int lumStride, int chromStride)
+{
+       /* Y Plane */
+       memcpy(ydst, ysrc, width*height);
+
+       /* XXX: implement upscaling for U,V */
+}
+
+static inline void RENAME(planar2x)(const uint8_t *src, uint8_t *dst, int srcWidth, int srcHeight, int srcStride, int dstStride)
+{
+       int x,y;
+       
+       dst[0]= src[0];
+        
+       // first line
+       for(x=0; x<srcWidth-1; x++){
+               dst[2*x+1]= (3*src[x] +   src[x+1])>>2;
+               dst[2*x+2]= (  src[x] + 3*src[x+1])>>2;
+       }
+       dst[2*srcWidth-1]= src[srcWidth-1];
+       
+        dst+= dstStride;
+
+       for(y=1; y<srcHeight; y++){
+#if defined (HAVE_MMX2) || defined (HAVE_3DNOW)
+               const int mmxSize= srcWidth&~15;
+               asm volatile(
+                       "movl %4, %%eax                 \n\t"
+                       "1:                             \n\t"
+                       "movq (%0, %%eax), %%mm0        \n\t"
+                       "movq (%1, %%eax), %%mm1        \n\t"
+                       "movq 1(%0, %%eax), %%mm2       \n\t"
+                       "movq 1(%1, %%eax), %%mm3       \n\t"
+                       "movq -1(%0, %%eax), %%mm4      \n\t"
+                       "movq -1(%1, %%eax), %%mm5      \n\t"
+                       PAVGB" %%mm0, %%mm5             \n\t"
+                       PAVGB" %%mm0, %%mm3             \n\t"
+                       PAVGB" %%mm0, %%mm5             \n\t"
+                       PAVGB" %%mm0, %%mm3             \n\t"
+                       PAVGB" %%mm1, %%mm4             \n\t"
+                       PAVGB" %%mm1, %%mm2             \n\t"
+                       PAVGB" %%mm1, %%mm4             \n\t"
+                       PAVGB" %%mm1, %%mm2             \n\t"
+                       "movq %%mm5, %%mm7              \n\t"
+                       "movq %%mm4, %%mm6              \n\t"
+                       "punpcklbw %%mm3, %%mm5         \n\t"
+                       "punpckhbw %%mm3, %%mm7         \n\t"
+                       "punpcklbw %%mm2, %%mm4         \n\t"
+                       "punpckhbw %%mm2, %%mm6         \n\t"
+#if 1
+                       MOVNTQ" %%mm5, (%2, %%eax, 2)   \n\t"
+                       MOVNTQ" %%mm7, 8(%2, %%eax, 2)  \n\t"
+                       MOVNTQ" %%mm4, (%3, %%eax, 2)   \n\t"
+                       MOVNTQ" %%mm6, 8(%3, %%eax, 2)  \n\t"
+#else
+                       "movq %%mm5, (%2, %%eax, 2)     \n\t"
+                       "movq %%mm7, 8(%2, %%eax, 2)    \n\t"
+                       "movq %%mm4, (%3, %%eax, 2)     \n\t"
+                       "movq %%mm6, 8(%3, %%eax, 2)    \n\t"
+#endif
+                       "addl $8, %%eax                 \n\t"
+                       " js 1b                         \n\t"
+                       :: "r" (src + mmxSize  ), "r" (src + srcStride + mmxSize  ),
+                          "r" (dst + mmxSize*2), "r" (dst + dstStride + mmxSize*2),
+                          "g" (-mmxSize)
+                       : "%eax"
+
+               );
+#else
+               const int mmxSize=1;
+#endif
+               dst[0        ]= (3*src[0] +   src[srcStride])>>2;
+               dst[dstStride]= (  src[0] + 3*src[srcStride])>>2;
+
+               for(x=mmxSize-1; x<srcWidth-1; x++){
+                       dst[2*x          +1]= (3*src[x+0] +   src[x+srcStride+1])>>2;
+                       dst[2*x+dstStride+2]= (  src[x+0] + 3*src[x+srcStride+1])>>2;
+                       dst[2*x+dstStride+1]= (  src[x+1] + 3*src[x+srcStride  ])>>2;
+                       dst[2*x          +2]= (3*src[x+1] +   src[x+srcStride  ])>>2;
+               }
+               dst[srcWidth*2 -1            ]= (3*src[srcWidth-1] +   src[srcWidth-1 + srcStride])>>2;
+               dst[srcWidth*2 -1 + dstStride]= (  src[srcWidth-1] + 3*src[srcWidth-1 + srcStride])>>2;
+
+               dst+=dstStride*2;
+               src+=srcStride;
+       }
+       
+       // last line
+#if 1
+       dst[0]= src[0];
+        
+       for(x=0; x<srcWidth-1; x++){
+               dst[2*x+1]= (3*src[x] +   src[x+1])>>2;
+               dst[2*x+2]= (  src[x] + 3*src[x+1])>>2;
+       }
+       dst[2*srcWidth-1]= src[srcWidth-1];
+#else
+       for(x=0; x<srcWidth; x++){
+               dst[2*x+0]=
+               dst[2*x+1]= src[x];
+       }
+#endif
+
+#ifdef HAVE_MMX
+asm volatile(   EMMS" \n\t"
+               SFENCE" \n\t"
+               :::"memory");
+#endif
+}
+
+/**
+ *
+ * height should be a multiple of 2 and width should be a multiple of 16 (if this is a
+ * problem for anyone then tell me, and ill fix it)
+ * chrominance data is only taken from every secound line others are ignored FIXME write HQ version
+ */
+static inline void RENAME(uyvytoyv12)(const uint8_t *src, uint8_t *ydst, uint8_t *udst, uint8_t *vdst,
+       unsigned int width, unsigned int height,
+       int lumStride, int chromStride, int srcStride)
+{
+       unsigned y;
+       const unsigned chromWidth= width>>1;
+       for(y=0; y<height; y+=2)
        {
-               dst[4*i+0] = ysrc[2*i+0];
-               dst[4*i+1] = usrc[i];
-               dst[4*i+2] = ysrc[2*i+1];
-               dst[4*i+3] = vsrc[i];
+#ifdef HAVE_MMX
+               asm volatile(
+                       "xorl %%eax, %%eax              \n\t"
+                       "pcmpeqw %%mm7, %%mm7           \n\t"
+                       "psrlw $8, %%mm7                \n\t" // FF,00,FF,00...
+                       ".balign 16                     \n\t"
+                       "1:                             \n\t"
+                       PREFETCH" 64(%0, %%eax, 4)      \n\t"
+                       "movq (%0, %%eax, 4), %%mm0     \n\t" // UYVY UYVY(0)
+                       "movq 8(%0, %%eax, 4), %%mm1    \n\t" // UYVY UYVY(4)
+                       "movq %%mm0, %%mm2              \n\t" // UYVY UYVY(0)
+                       "movq %%mm1, %%mm3              \n\t" // UYVY UYVY(4)
+                       "pand %%mm7, %%mm0              \n\t" // U0V0 U0V0(0)
+                       "pand %%mm7, %%mm1              \n\t" // U0V0 U0V0(4)
+                       "psrlw $8, %%mm2                \n\t" // Y0Y0 Y0Y0(0)
+                       "psrlw $8, %%mm3                \n\t" // Y0Y0 Y0Y0(4)
+                       "packuswb %%mm1, %%mm0          \n\t" // UVUV UVUV(0)
+                       "packuswb %%mm3, %%mm2          \n\t" // YYYY YYYY(0)
+
+                       MOVNTQ" %%mm2, (%1, %%eax, 2)   \n\t"
+
+                       "movq 16(%0, %%eax, 4), %%mm1   \n\t" // UYVY UYVY(8)
+                       "movq 24(%0, %%eax, 4), %%mm2   \n\t" // UYVY UYVY(12)
+                       "movq %%mm1, %%mm3              \n\t" // UYVY UYVY(8)
+                       "movq %%mm2, %%mm4              \n\t" // UYVY UYVY(12)
+                       "pand %%mm7, %%mm1              \n\t" // U0V0 U0V0(8)
+                       "pand %%mm7, %%mm2              \n\t" // U0V0 U0V0(12)
+                       "psrlw $8, %%mm3                \n\t" // Y0Y0 Y0Y0(8)
+                       "psrlw $8, %%mm4                \n\t" // Y0Y0 Y0Y0(12)
+                       "packuswb %%mm2, %%mm1          \n\t" // UVUV UVUV(8)
+                       "packuswb %%mm4, %%mm3          \n\t" // YYYY YYYY(8)
+
+                       MOVNTQ" %%mm3, 8(%1, %%eax, 2)  \n\t"
+
+                       "movq %%mm0, %%mm2              \n\t" // UVUV UVUV(0)
+                       "movq %%mm1, %%mm3              \n\t" // UVUV UVUV(8)
+                       "psrlw $8, %%mm0                \n\t" // V0V0 V0V0(0)
+                       "psrlw $8, %%mm1                \n\t" // V0V0 V0V0(8)
+                       "pand %%mm7, %%mm2              \n\t" // U0U0 U0U0(0)
+                       "pand %%mm7, %%mm3              \n\t" // U0U0 U0U0(8)
+                       "packuswb %%mm1, %%mm0          \n\t" // VVVV VVVV(0)
+                       "packuswb %%mm3, %%mm2          \n\t" // UUUU UUUU(0)
+
+                       MOVNTQ" %%mm0, (%3, %%eax)      \n\t"
+                       MOVNTQ" %%mm2, (%2, %%eax)      \n\t"
+
+                       "addl $8, %%eax                 \n\t"
+                       "cmpl %4, %%eax                 \n\t"
+                       " jb 1b                         \n\t"
+                       ::"r"(src), "r"(ydst), "r"(udst), "r"(vdst), "g" (chromWidth)
+                       : "memory", "%eax"
+               );
+
+               ydst += lumStride;
+               src  += srcStride;
+
+               asm volatile(
+                       "xorl %%eax, %%eax              \n\t"
+                       ".balign 16                     \n\t"
+                       "1:                             \n\t"
+                       PREFETCH" 64(%0, %%eax, 4)      \n\t"
+                       "movq (%0, %%eax, 4), %%mm0     \n\t" // YUYV YUYV(0)
+                       "movq 8(%0, %%eax, 4), %%mm1    \n\t" // YUYV YUYV(4)
+                       "movq 16(%0, %%eax, 4), %%mm2   \n\t" // YUYV YUYV(8)
+                       "movq 24(%0, %%eax, 4), %%mm3   \n\t" // YUYV YUYV(12)
+                       "psrlw $8, %%mm0                \n\t" // Y0Y0 Y0Y0(0)
+                       "psrlw $8, %%mm1                \n\t" // Y0Y0 Y0Y0(4)
+                       "psrlw $8, %%mm2                \n\t" // Y0Y0 Y0Y0(8)
+                       "psrlw $8, %%mm3                \n\t" // Y0Y0 Y0Y0(12)
+                       "packuswb %%mm1, %%mm0          \n\t" // YYYY YYYY(0)
+                       "packuswb %%mm3, %%mm2          \n\t" // YYYY YYYY(8)
+
+                       MOVNTQ" %%mm0, (%1, %%eax, 2)   \n\t"
+                       MOVNTQ" %%mm2, 8(%1, %%eax, 2)  \n\t"
+
+                       "addl $8, %%eax                 \n\t"
+                       "cmpl %4, %%eax                 \n\t"
+                       " jb 1b                         \n\t"
+
+                       ::"r"(src), "r"(ydst), "r"(udst), "r"(vdst), "g" (chromWidth)
+                       : "memory", "%eax"
+               );
+#else
+               unsigned i;
+               for(i=0; i<chromWidth; i++)
+               {
+                       udst[i]         = src[4*i+0];
+                       ydst[2*i+0]     = src[4*i+1];
+                       vdst[i]         = src[4*i+2];
+                       ydst[2*i+1]     = src[4*i+3];
+               }
+               ydst += lumStride;
+               src  += srcStride;
+
+               for(i=0; i<chromWidth; i++)
+               {
+                       ydst[2*i+0]     = src[4*i+1];
+                       ydst[2*i+1]     = src[4*i+3];
+               }
+#endif
+               udst += chromStride;
+               vdst += chromStride;
+               ydst += lumStride;
+               src  += srcStride;
        }
+#ifdef HAVE_MMX
+asm volatile(   EMMS" \n\t"
+               SFENCE" \n\t"
+               :::"memory");
 #endif
 }
 
-void yuy2toyv12(const uint8_t *src, uint8_t *ydst, uint8_t *udst, uint8_t *vdst, uint32_t num_pixels)
+/**
+ *
+ * height should be a multiple of 2 and width should be a multiple of 2 (if this is a
+ * problem for anyone then tell me, and ill fix it)
+ * chrominance data is only taken from every secound line others are ignored in the C version FIXME write HQ version
+ */
+static inline void RENAME(rgb24toyv12)(const uint8_t *src, uint8_t *ydst, uint8_t *udst, uint8_t *vdst,
+       unsigned int width, unsigned int height,
+       int lumStride, int chromStride, int srcStride)
 {
+       unsigned y;
+       const unsigned chromWidth= width>>1;
 #ifdef HAVE_MMX
-       asm volatile(
-               "xorl %%eax, %%eax              \n\t"
-               "pcmpeqw %%mm7, %%mm7           \n\t"
-               "psrlw $8, %%mm7                \n\t" // FF,00,FF,00...
-               "1:                             \n\t"
-               PREFETCH" 64(%0, %%eax, 4)      \n\t"
-               "movq (%0, %%eax, 4), %%mm0     \n\t" // YUYV YUYV(0)
-               "movq 8(%0, %%eax, 4), %%mm1    \n\t" // YUYV YUYV(4)
-               "movq %%mm0, %%mm2              \n\t" // YUYV YUYV(0)
-               "movq %%mm1, %%mm3              \n\t" // YUYV YUYV(4)
-               "psrlw $8, %%mm0                \n\t" // U0V0 U0V0(0)
-               "psrlw $8, %%mm1                \n\t" // U0V0 U0V0(4)
-               "pand %%mm7, %%mm2              \n\t" // Y0Y0 Y0Y0(0)
-               "pand %%mm7, %%mm3              \n\t" // Y0Y0 Y0Y0(4)
-               "packuswb %%mm1, %%mm0          \n\t" // UVUV UVUV(0)
-               "packuswb %%mm3, %%mm2          \n\t" // YYYY YYYY(0)
-
-               MOVNTQ" %%mm2, (%1, %%eax, 2)   \n\t"
-
-               "movq 16(%0, %%eax, 4), %%mm1   \n\t" // YUYV YUYV(8)
-               "movq 24(%0, %%eax, 4), %%mm2   \n\t" // YUYV YUYV(12)
-               "movq %%mm1, %%mm3              \n\t" // YUYV YUYV(8)
-               "movq %%mm2, %%mm4              \n\t" // YUYV YUYV(12)
-               "psrlw $8, %%mm1                \n\t" // U0V0 U0V0(8)
-               "psrlw $8, %%mm2                \n\t" // U0V0 U0V0(12)
-               "pand %%mm7, %%mm3              \n\t" // Y0Y0 Y0Y0(8)
-               "pand %%mm7, %%mm4              \n\t" // Y0Y0 Y0Y0(12)
-               "packuswb %%mm2, %%mm1          \n\t" // UVUV UVUV(8)
-               "packuswb %%mm4, %%mm3          \n\t" // YYYY YYYY(8)
-
-               MOVNTQ" %%mm3, 8(%1, %%eax, 2)  \n\t"
-
-               "movq %%mm0, %%mm2              \n\t" // UVUV UVUV(0)
-               "movq %%mm1, %%mm3              \n\t" // UVUV UVUV(8)
-               "psrlw $8, %%mm0                \n\t" // V0V0 V0V0(0)
-               "psrlw $8, %%mm1                \n\t" // V0V0 V0V0(8)
-               "pand %%mm7, %%mm2              \n\t" // U0U0 U0U0(0)
-               "pand %%mm7, %%mm3              \n\t" // U0U0 U0U0(8)
-               "packuswb %%mm1, %%mm0          \n\t" // VVVV VVVV(0)
-               "packuswb %%mm3, %%mm2          \n\t" // UUUU UUUU(0)
-
-               MOVNTQ" %%mm0, (%3, %%eax)      \n\t"
-               MOVNTQ" %%mm2, (%2, %%eax)      \n\t"
+       for(y=0; y<height-2; y+=2)
+       {
+               unsigned i;
+               for(i=0; i<2; i++)
+               {
+                       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"
 
-               "addl $8, %%eax                 \n\t"
-               "cmpl %4, %%eax                 \n\t"
-               " jb 1b                         \n\t"
-               EMMS" \n\t"
-               SFENCE
-               ::"r"(src), "r"(ydst), "r"(udst), "r"(vdst), "r" (num_pixels>>1)
-               : "memory", "%eax"
-       );
+                               "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"
+
+                               MOVNTQ" %%mm0, (%1, %%eax)      \n\t"
+                               "addl $8, %%eax                 \n\t"
+                               " js 1b                         \n\t"
+                               : : "r" (src+width*3), "r" (ydst+width), "g" (-width)
+                               : "%eax", "%ebx"
+                       );
+                       ydst += lumStride;
+                       src  += srcStride;
+               }
+               src -= srcStride*2;
+               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             \n\t"
+                       PAVGB" %%mm3, %%mm2             \n\t"
+                       "movq %%mm0, %%mm1              \n\t"
+                       "movq %%mm2, %%mm3              \n\t"
+                       "psrlq $24, %%mm0               \n\t"
+                       "psrlq $24, %%mm2               \n\t"
+                       PAVGB" %%mm1, %%mm0             \n\t"
+                       PAVGB" %%mm3, %%mm2             \n\t"
+                       "punpcklbw %%mm7, %%mm0         \n\t"
+                       "punpcklbw %%mm7, %%mm2         \n\t"
 #else
-       int i;
-       num_pixels>>=1;
-       for(i=0; i<num_pixels; i++)
+                       "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             \n\t"
+                       PAVGB" %%mm3, %%mm2             \n\t"
+                       "movq %%mm4, %%mm1              \n\t"
+                       "movq %%mm2, %%mm3              \n\t"
+                       "psrlq $24, %%mm4               \n\t"
+                       "psrlq $24, %%mm2               \n\t"
+                       PAVGB" %%mm1, %%mm4             \n\t"
+                       PAVGB" %%mm3, %%mm2             \n\t"
+                       "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" (src+width*6), "r" (src+srcStride+width*6), "r" (udst+width), "r" (vdst+width), "g" (-width)
+                       : "%eax", "%ebx"
+               );
+
+               udst += chromStride;
+               vdst += chromStride;
+               src  += srcStride*2;
+       }
+
+       asm volatile(   EMMS" \n\t"
+                       SFENCE" \n\t"
+                       :::"memory");
+#else
+       y=0;
+#endif
+       for(; y<height; y+=2)
+       {
+               unsigned i;
+               for(i=0; i<chromWidth; i++)
+               {
+                       unsigned int b= src[6*i+0];
+                       unsigned int g= src[6*i+1];
+                       unsigned int r= src[6*i+2];
+
+                       unsigned int Y  =  ((RY*r + GY*g + BY*b)>>RGB2YUV_SHIFT) + 16;
+                       unsigned int V  =  ((RV*r + GV*g + BV*b)>>RGB2YUV_SHIFT) + 128;
+                       unsigned int U  =  ((RU*r + GU*g + BU*b)>>RGB2YUV_SHIFT) + 128;
+
+                       udst[i]         = U;
+                       vdst[i]         = V;
+                       ydst[2*i]       = Y;
+
+                       b= src[6*i+3];
+                       g= src[6*i+4];
+                       r= src[6*i+5];
+
+                       Y  =  ((RY*r + GY*g + BY*b)>>RGB2YUV_SHIFT) + 16;
+                       ydst[2*i+1]     = Y;
+               }
+               ydst += lumStride;
+               src  += srcStride;
+
+               for(i=0; i<chromWidth; i++)
+               {
+                       unsigned int b= src[6*i+0];
+                       unsigned int g= src[6*i+1];
+                       unsigned int r= src[6*i+2];
+
+                       unsigned int Y  =  ((RY*r + GY*g + BY*b)>>RGB2YUV_SHIFT) + 16;
+
+                       ydst[2*i]       = Y;
+
+                       b= src[6*i+3];
+                       g= src[6*i+4];
+                       r= src[6*i+5];
+
+                       Y  =  ((RY*r + GY*g + BY*b)>>RGB2YUV_SHIFT) + 16;
+                       ydst[2*i+1]     = Y;
+               }
+               udst += chromStride;
+               vdst += chromStride;
+               ydst += lumStride;
+               src  += srcStride;
+       }
+}
+
+void RENAME(interleaveBytes)(uint8_t *src1, uint8_t *src2, uint8_t *dest,
+                           unsigned width, unsigned height, int src1Stride,
+                           int src2Stride, int dstStride){
+       unsigned h;
+
+       for(h=0; h < height; h++)
        {
-                ydst[2*i+0]    = src[4*i+0];
-                udst[i]        = src[4*i+1];
-                ydst[2*i+1]    = src[4*i+2];
-                vdst[i]        = src[4*i+3];
+               unsigned w;
+
+#ifdef HAVE_MMX
+#ifdef HAVE_SSE2
+               asm(
+                       "xorl %%eax, %%eax              \n\t"
+                       "1:                             \n\t"
+                       PREFETCH" 64(%1, %%eax)         \n\t"
+                       PREFETCH" 64(%2, %%eax)         \n\t"
+                       "movdqa (%1, %%eax), %%xmm0     \n\t"
+                       "movdqa (%1, %%eax), %%xmm1     \n\t"
+                       "movdqa (%2, %%eax), %%xmm2     \n\t"
+                       "punpcklbw %%xmm2, %%xmm0       \n\t"
+                       "punpckhbw %%xmm2, %%xmm1       \n\t"
+                       "movntdq %%xmm0, (%0, %%eax, 2) \n\t"
+                       "movntdq %%xmm1, 16(%0, %%eax, 2)\n\t"
+                       "addl $16, %%eax                        \n\t"
+                       "cmpl %3, %%eax                 \n\t"
+                       " jb 1b                         \n\t"
+                       ::"r"(dest), "r"(src1), "r"(src2), "r" (width-15)
+                       : "memory", "%eax"
+               );
+#else
+               asm(
+                       "xorl %%eax, %%eax              \n\t"
+                       "1:                             \n\t"
+                       PREFETCH" 64(%1, %%eax)         \n\t"
+                       PREFETCH" 64(%2, %%eax)         \n\t"
+                       "movq (%1, %%eax), %%mm0        \n\t"
+                       "movq 8(%1, %%eax), %%mm2       \n\t"
+                       "movq %%mm0, %%mm1              \n\t"
+                       "movq %%mm2, %%mm3              \n\t"
+                       "movq (%2, %%eax), %%mm4        \n\t"
+                       "movq 8(%2, %%eax), %%mm5       \n\t"
+                       "punpcklbw %%mm4, %%mm0         \n\t"
+                       "punpckhbw %%mm4, %%mm1         \n\t"
+                       "punpcklbw %%mm5, %%mm2         \n\t"
+                       "punpckhbw %%mm5, %%mm3         \n\t"
+                       MOVNTQ" %%mm0, (%0, %%eax, 2)   \n\t"
+                       MOVNTQ" %%mm1, 8(%0, %%eax, 2)  \n\t"
+                       MOVNTQ" %%mm2, 16(%0, %%eax, 2) \n\t"
+                       MOVNTQ" %%mm3, 24(%0, %%eax, 2) \n\t"
+                       "addl $16, %%eax                        \n\t"
+                       "cmpl %3, %%eax                 \n\t"
+                       " jb 1b                         \n\t"
+                       ::"r"(dest), "r"(src1), "r"(src2), "r" (width-15)
+                       : "memory", "%eax"
+               );
+#endif
+               for(w= (width&(~15)); w < width; w++)
+               {
+                       dest[2*w+0] = src1[w];
+                       dest[2*w+1] = src2[w];
+               }
+#else
+               for(w=0; w < width; w++)
+               {
+                       dest[2*w+0] = src1[w];
+                       dest[2*w+1] = src2[w];
+               }
+#endif
+               dest += dstStride;
+                src1 += src1Stride;
+                src2 += src2Stride;
        }
+#ifdef HAVE_MMX
+       asm(
+               EMMS" \n\t"
+               SFENCE" \n\t"
+               ::: "memory"
+               );
+#endif
+}
+
+static inline void RENAME(vu9_to_vu12)(const uint8_t *src1, const uint8_t *src2,
+                       uint8_t *dst1, uint8_t *dst2,
+                       unsigned width, unsigned height,
+                       int srcStride1, int srcStride2,
+                       int dstStride1, int dstStride2)
+{
+    unsigned int y,x,h;
+    int w;
+    w=width/2; h=height/2;
+#ifdef HAVE_MMX
+    asm volatile(
+       PREFETCH" %0\n\t"
+       PREFETCH" %1\n\t"
+       ::"m"(*(src1+srcStride1)),"m"(*(src2+srcStride2)):"memory");
 #endif
-}
\ No newline at end of file
+    for(y=0;y<h;y++){
+       const uint8_t* s1=src1+srcStride1*(y>>1);
+       uint8_t* d=dst1+dstStride1*y;
+       x=0;
+#ifdef HAVE_MMX
+       for(;x<w-31;x+=32)
+       {
+           asm volatile(
+               PREFETCH" 32%1\n\t"
+               "movq   %1, %%mm0\n\t"
+               "movq   8%1, %%mm2\n\t"
+               "movq   16%1, %%mm4\n\t"
+               "movq   24%1, %%mm6\n\t"
+               "movq   %%mm0, %%mm1\n\t"
+               "movq   %%mm2, %%mm3\n\t"
+               "movq   %%mm4, %%mm5\n\t"
+               "movq   %%mm6, %%mm7\n\t"
+               "punpcklbw %%mm0, %%mm0\n\t"
+               "punpckhbw %%mm1, %%mm1\n\t"
+               "punpcklbw %%mm2, %%mm2\n\t"
+               "punpckhbw %%mm3, %%mm3\n\t"
+               "punpcklbw %%mm4, %%mm4\n\t"
+               "punpckhbw %%mm5, %%mm5\n\t"
+               "punpcklbw %%mm6, %%mm6\n\t"
+               "punpckhbw %%mm7, %%mm7\n\t"
+               MOVNTQ" %%mm0, %0\n\t"
+               MOVNTQ" %%mm1, 8%0\n\t"
+               MOVNTQ" %%mm2, 16%0\n\t"
+               MOVNTQ" %%mm3, 24%0\n\t"
+               MOVNTQ" %%mm4, 32%0\n\t"
+               MOVNTQ" %%mm5, 40%0\n\t"
+               MOVNTQ" %%mm6, 48%0\n\t"
+               MOVNTQ" %%mm7, 56%0"
+               :"=m"(d[2*x])
+               :"m"(s1[x])
+               :"memory");
+       }
+#endif
+       for(;x<w;x++) d[2*x]=d[2*x+1]=s1[x];
+    }
+    for(y=0;y<h;y++){
+       const uint8_t* s2=src2+srcStride2*(y>>1);
+       uint8_t* d=dst2+dstStride2*y;
+       x=0;
+#ifdef HAVE_MMX
+       for(;x<w-31;x+=32)
+       {
+           asm volatile(
+               PREFETCH" 32%1\n\t"
+               "movq   %1, %%mm0\n\t"
+               "movq   8%1, %%mm2\n\t"
+               "movq   16%1, %%mm4\n\t"
+               "movq   24%1, %%mm6\n\t"
+               "movq   %%mm0, %%mm1\n\t"
+               "movq   %%mm2, %%mm3\n\t"
+               "movq   %%mm4, %%mm5\n\t"
+               "movq   %%mm6, %%mm7\n\t"
+               "punpcklbw %%mm0, %%mm0\n\t"
+               "punpckhbw %%mm1, %%mm1\n\t"
+               "punpcklbw %%mm2, %%mm2\n\t"
+               "punpckhbw %%mm3, %%mm3\n\t"
+               "punpcklbw %%mm4, %%mm4\n\t"
+               "punpckhbw %%mm5, %%mm5\n\t"
+               "punpcklbw %%mm6, %%mm6\n\t"
+               "punpckhbw %%mm7, %%mm7\n\t"
+               MOVNTQ" %%mm0, %0\n\t"
+               MOVNTQ" %%mm1, 8%0\n\t"
+               MOVNTQ" %%mm2, 16%0\n\t"
+               MOVNTQ" %%mm3, 24%0\n\t"
+               MOVNTQ" %%mm4, 32%0\n\t"
+               MOVNTQ" %%mm5, 40%0\n\t"
+               MOVNTQ" %%mm6, 48%0\n\t"
+               MOVNTQ" %%mm7, 56%0"
+               :"=m"(d[2*x])
+               :"m"(s2[x])
+               :"memory");
+       }
+#endif
+       for(;x<w;x++) d[2*x]=d[2*x+1]=s2[x];
+    }
+#ifdef HAVE_MMX
+       asm(
+               EMMS" \n\t"
+               SFENCE" \n\t"
+               ::: "memory"
+               );
+#endif
+}
+
+static inline void RENAME(yvu9_to_yuy2)(const uint8_t *src1, const uint8_t *src2, const uint8_t *src3,
+                       uint8_t *dst,
+                       unsigned width, unsigned height,
+                       int srcStride1, int srcStride2,
+                       int srcStride3, int dstStride)
+{
+    unsigned y,x,w,h;
+    w=width/2; h=height;
+    for(y=0;y<h;y++){
+       const uint8_t* yp=src1+srcStride1*y;
+       const uint8_t* up=src2+srcStride2*(y>>2);
+       const uint8_t* vp=src3+srcStride3*(y>>2);
+       uint8_t* d=dst+dstStride*y;
+       x=0;
+#ifdef HAVE_MMX
+       for(;x<w-7;x+=8)
+       {
+           asm volatile(
+               PREFETCH" 32(%1, %0)\n\t"
+               PREFETCH" 32(%2, %0)\n\t"
+               PREFETCH" 32(%3, %0)\n\t"
+               "movq   (%1, %0, 4), %%mm0\n\t"       /* Y0Y1Y2Y3Y4Y5Y6Y7 */
+               "movq   (%2, %0), %%mm1\n\t"       /* U0U1U2U3U4U5U6U7 */
+               "movq   (%3, %0), %%mm2\n\t"         /* V0V1V2V3V4V5V6V7 */
+               "movq   %%mm0, %%mm3\n\t"    /* Y0Y1Y2Y3Y4Y5Y6Y7 */
+               "movq   %%mm1, %%mm4\n\t"    /* U0U1U2U3U4U5U6U7 */
+               "movq   %%mm2, %%mm5\n\t"    /* V0V1V2V3V4V5V6V7 */
+               "punpcklbw %%mm1, %%mm1\n\t" /* U0U0 U1U1 U2U2 U3U3 */
+               "punpcklbw %%mm2, %%mm2\n\t" /* V0V0 V1V1 V2V2 V3V3 */
+               "punpckhbw %%mm4, %%mm4\n\t" /* U4U4 U5U5 U6U6 U7U7 */
+               "punpckhbw %%mm5, %%mm5\n\t" /* V4V4 V5V5 V6V6 V7V7 */
+
+               "movq   %%mm1, %%mm6\n\t"
+               "punpcklbw %%mm2, %%mm1\n\t" /* U0V0 U0V0 U1V1 U1V1*/
+               "punpcklbw %%mm1, %%mm0\n\t" /* Y0U0 Y1V0 Y2U0 Y3V0*/
+               "punpckhbw %%mm1, %%mm3\n\t" /* Y4U1 Y5V1 Y6U1 Y7V1*/
+               MOVNTQ" %%mm0, (%4, %0, 8)\n\t"
+               MOVNTQ" %%mm3, 8(%4, %0, 8)\n\t"
+               
+               "punpckhbw %%mm2, %%mm6\n\t" /* U2V2 U2V2 U3V3 U3V3*/
+               "movq   8(%1, %0, 4), %%mm0\n\t"
+               "movq   %%mm0, %%mm3\n\t"
+               "punpcklbw %%mm6, %%mm0\n\t" /* Y U2 Y V2 Y U2 Y V2*/
+               "punpckhbw %%mm6, %%mm3\n\t" /* Y U3 Y V3 Y U3 Y V3*/
+               MOVNTQ" %%mm0, 16(%4, %0, 8)\n\t"
+               MOVNTQ" %%mm3, 24(%4, %0, 8)\n\t"
+
+               "movq   %%mm4, %%mm6\n\t"
+               "movq   16(%1, %0, 4), %%mm0\n\t"
+               "movq   %%mm0, %%mm3\n\t"
+               "punpcklbw %%mm5, %%mm4\n\t"
+               "punpcklbw %%mm4, %%mm0\n\t" /* Y U4 Y V4 Y U4 Y V4*/
+               "punpckhbw %%mm4, %%mm3\n\t" /* Y U5 Y V5 Y U5 Y V5*/
+               MOVNTQ" %%mm0, 32(%4, %0, 8)\n\t"
+               MOVNTQ" %%mm3, 40(%4, %0, 8)\n\t"
+               
+               "punpckhbw %%mm5, %%mm6\n\t"
+               "movq   24(%1, %0, 4), %%mm0\n\t"
+               "movq   %%mm0, %%mm3\n\t"
+               "punpcklbw %%mm6, %%mm0\n\t" /* Y U6 Y V6 Y U6 Y V6*/
+               "punpckhbw %%mm6, %%mm3\n\t" /* Y U7 Y V7 Y U7 Y V7*/
+               MOVNTQ" %%mm0, 48(%4, %0, 8)\n\t"
+               MOVNTQ" %%mm3, 56(%4, %0, 8)\n\t"
+
+               : "+r" (x)
+                : "r"(yp), "r" (up), "r"(vp), "r"(d)
+               :"memory");
+       }
+#endif
+       for(; x<w; x++)
+       {
+           const int x2= x<<2;
+           d[8*x+0]=yp[x2];
+           d[8*x+1]=up[x];
+           d[8*x+2]=yp[x2+1];
+           d[8*x+3]=vp[x];
+           d[8*x+4]=yp[x2+2];
+           d[8*x+5]=up[x];
+           d[8*x+6]=yp[x2+3];
+           d[8*x+7]=vp[x];
+       }
+    }
+#ifdef HAVE_MMX
+       asm(
+               EMMS" \n\t"
+               SFENCE" \n\t"
+               ::: "memory"
+               );
+#endif
+}