Convert DOS line endings to Unix line endings

Should fix issues some of us have with `misc/dist/uwp_template/AppxManifest.xml`
always showing up as modified.
Might cause issues on Windows due to the removal of BOMs or change of line endings
in some of the Mono, UWP or gradlew.bat files, we will test and adapt if need be.

1 files changed, 670 insertions, 670 deletions

diff --git a/thirdparty/libtheora/x86_vc/mmxfdct.c b/thirdparty/libtheora/x86_vc/mmxfdct.c
index dcf17c9fa7..d908ce2413 100644
--- a/thirdparty/libtheora/x86_vc/mmxfdct.c
+++ b/thirdparty/libtheora/x86_vc/mmxfdct.c
@@ -1,670 +1,670 @@
-/********************************************************************

- *                                                                  *

- * THIS FILE IS PART OF THE OggTheora SOFTWARE CODEC SOURCE CODE.   *

- * USE, DISTRIBUTION AND REPRODUCTION OF THIS LIBRARY SOURCE IS     *

- * GOVERNED BY A BSD-STYLE SOURCE LICENSE INCLUDED WITH THIS SOURCE *

- * IN 'COPYING'. PLEASE READ THESE TERMS BEFORE DISTRIBUTING.       *

- *                                                                  *

- * THE Theora SOURCE CODE IS COPYRIGHT (C) 1999-2006                *

- * by the Xiph.Org Foundation http://www.xiph.org/                  *

- *                                                                  *

- ********************************************************************/ 

- /*MMX fDCT implementation for x86_32*/

-/*$Id: fdct_ses2.c 14579 2008-03-12 06:42:40Z xiphmont $*/

-#include "x86enc.h"

-

-#if defined(OC_X86_ASM)

-

-#define OC_FDCT_STAGE1_8x4  __asm{ \

-  /*Stage 1:*/ \

-  /*mm0=t7'=t0-t7*/ \

-  __asm  psubw mm0,mm7 \

-  __asm  paddw mm7,mm7 \

-  /*mm1=t6'=t1-t6*/ \

-  __asm  psubw mm1, mm6 \

-  __asm  paddw mm6,mm6 \

-  /*mm2=t5'=t2-t5*/ \

-  __asm  psubw mm2,mm5 \

-  __asm  paddw mm5,mm5 \

-  /*mm3=t4'=t3-t4*/ \

-  __asm  psubw mm3,mm4 \

-  __asm  paddw mm4,mm4 \

-  /*mm7=t0'=t0+t7*/ \

-  __asm  paddw mm7,mm0 \

-  /*mm6=t1'=t1+t6*/  \

-  __asm  paddw mm6,mm1 \

-  /*mm5=t2'=t2+t5*/ \

-  __asm  paddw mm5,mm2 \

-  /*mm4=t3'=t3+t4*/ \

-  __asm  paddw mm4,mm3\

-}

-

-#define OC_FDCT8x4(_r0,_r1,_r2,_r3,_r4,_r5,_r6,_r7) __asm{ \

-  /*Stage 2:*/ \

-  /*mm7=t3''=t0'-t3'*/ \

-  __asm  psubw mm7,mm4 \

-  __asm  paddw mm4,mm4 \

-  /*mm6=t2''=t1'-t2'*/ \

-  __asm  psubw mm6,mm5 \

-  __asm  movq [Y+_r6],mm7 \

-  __asm  paddw mm5,mm5 \

-  /*mm1=t5''=t6'-t5'*/ \

-  __asm  psubw mm1,mm2 \

-  __asm  movq [Y+_r2],mm6 \

-  /*mm4=t0''=t0'+t3'*/ \

-  __asm  paddw mm4,mm7 \

-  __asm  paddw mm2,mm2 \

-  /*mm5=t1''=t1'+t2'*/ \

-  __asm  movq [Y+_r0],mm4 \

-  __asm  paddw mm5,mm6 \

-  /*mm2=t6''=t6'+t5'*/ \

-  __asm  paddw mm2,mm1 \

-  __asm  movq [Y+_r4],mm5 \

-  /*mm0=t7', mm1=t5'', mm2=t6'', mm3=t4'.*/ \

-  /*mm4, mm5, mm6, mm7 are free.*/ \

-  /*Stage 3:*/ \

-  /*mm6={2}x4, mm7={27146,0xB500>>1}x2*/ \

-  __asm  mov A,0x5A806A0A \

-  __asm  pcmpeqb mm6,mm6 \

-  __asm  movd mm7,A \

-  __asm  psrlw mm6,15 \

-  __asm  punpckldq mm7,mm7 \

-  __asm  paddw mm6,mm6 \

-  /*mm0=0, m2={-1}x4 \

-    mm5:mm4=t5''*27146+0xB500*/ \

-  __asm  movq mm4,mm1 \

-  __asm  movq mm5,mm1 \

-  __asm  punpcklwd mm4,mm6 \

-  __asm  movq [Y+_r3],mm2 \

-  __asm  pmaddwd mm4,mm7 \

-  __asm  movq [Y+_r7],mm0 \

-  __asm  punpckhwd mm5,mm6 \

-  __asm  pxor mm0,mm0 \

-  __asm  pmaddwd mm5,mm7 \

-  __asm  pcmpeqb mm2,mm2 \

-  /*mm2=t6'', mm1=t5''+(t5''!=0) \

-    mm4=(t5''*27146+0xB500>>16)*/ \

-  __asm  pcmpeqw mm0,mm1 \

-  __asm  psrad mm4,16 \

-  __asm  psubw mm0,mm2 \

-  __asm  movq mm2, [Y+_r3] \

-  __asm  psrad mm5,16 \

-  __asm  paddw mm1,mm0 \

-  __asm  packssdw mm4,mm5 \

-  /*mm4=s=(t5''*27146+0xB500>>16)+t5''+(t5''!=0)>>1*/ \

-  __asm  paddw mm4,mm1 \

-  __asm  movq mm0, [Y+_r7] \

-  __asm  psraw mm4,1 \

-  __asm  movq mm1,mm3 \

-  /*mm3=t4''=t4'+s*/ \

-  __asm  paddw mm3,mm4 \

-  /*mm1=t5'''=t4'-s*/ \

-  __asm  psubw mm1,mm4 \

-  /*mm1=0, mm3={-1}x4 \

-    mm5:mm4=t6''*27146+0xB500*/ \

-  __asm  movq mm4,mm2 \

-  __asm  movq mm5,mm2 \

-  __asm  punpcklwd mm4,mm6 \

-  __asm  movq [Y+_r5],mm1 \

-  __asm  pmaddwd mm4,mm7 \

-  __asm  movq [Y+_r1],mm3 \

-  __asm  punpckhwd mm5,mm6 \

-  __asm  pxor mm1,mm1 \

-  __asm  pmaddwd mm5,mm7 \

-  __asm  pcmpeqb mm3,mm3 \

-  /*mm2=t6''+(t6''!=0), mm4=(t6''*27146+0xB500>>16)*/ \

-  __asm  psrad mm4,16 \

-  __asm  pcmpeqw mm1,mm2 \

-  __asm  psrad mm5,16 \

-  __asm  psubw mm1,mm3 \

-  __asm  packssdw mm4,mm5 \

-  __asm  paddw mm2,mm1 \

-  /*mm1=t1'' \

-    mm4=s=(t6''*27146+0xB500>>16)+t6''+(t6''!=0)>>1*/ \

-  __asm  paddw mm4,mm2 \

-  __asm  movq mm1,[Y+_r4] \

-  __asm  psraw mm4,1 \

-  __asm  movq mm2,mm0 \

-  /*mm7={54491-0x7FFF,0x7FFF}x2 \

-    mm0=t7''=t7'+s*/ \

-  __asm  paddw mm0,mm4 \

-  /*mm2=t6'''=t7'-s*/ \

-  __asm  psubw mm2,mm4 \

-  /*Stage 4:*/ \

-  /*mm0=0, mm2=t0'' \

-    mm5:mm4=t1''*27146+0xB500*/ \

-  __asm  movq mm4,mm1 \

-  __asm  movq mm5,mm1 \

-  __asm  punpcklwd mm4,mm6 \

-  __asm  movq [Y+_r3],mm2 \

-  __asm  pmaddwd mm4,mm7 \

-  __asm  movq mm2,[Y+_r0] \

-  __asm  punpckhwd mm5,mm6 \

-  __asm  movq [Y+_r7],mm0 \

-  __asm  pmaddwd mm5,mm7 \

-  __asm  pxor mm0,mm0 \

-  /*mm7={27146,0x4000>>1}x2 \

-    mm0=s=(t1''*27146+0xB500>>16)+t1''+(t1''!=0)*/ \

-  __asm  psrad mm4,16 \

-  __asm  mov A,0x20006A0A \

-  __asm  pcmpeqw mm0,mm1 \

-  __asm  movd mm7,A \

-  __asm  psrad mm5,16 \

-  __asm  psubw mm0,mm3 \

-  __asm  packssdw mm4,mm5 \

-  __asm  paddw mm0,mm1 \

-  __asm  punpckldq mm7,mm7 \

-  __asm  paddw mm0,mm4 \

-  /*mm6={0x00000E3D}x2 \

-    mm1=-(t0''==0), mm5:mm4=t0''*27146+0x4000*/ \

-  __asm  movq mm4,mm2 \

-  __asm  movq mm5,mm2 \

-  __asm  punpcklwd mm4,mm6 \

-  __asm  mov A,0x0E3D \

-  __asm  pmaddwd mm4,mm7 \

-  __asm  punpckhwd mm5,mm6 \

-  __asm  movd mm6,A \

-  __asm  pmaddwd mm5,mm7 \

-  __asm  pxor mm1,mm1 \

-  __asm  punpckldq mm6,mm6 \

-  __asm  pcmpeqw mm1,mm2 \

-  /*mm4=r=(t0''*27146+0x4000>>16)+t0''+(t0''!=0)*/ \

-  __asm  psrad mm4,16 \

-  __asm  psubw mm1,mm3 \

-  __asm  psrad mm5,16 \

-  __asm  paddw mm2,mm1 \

-  __asm  packssdw mm4,mm5 \

-  __asm  movq mm1,[Y+_r5] \

-  __asm  paddw mm4,mm2 \

-  /*mm2=t6'', mm0=_y[0]=u=r+s>>1 \

-    The naive implementation could cause overflow, so we use \

-     u=(r&s)+((r^s)>>1).*/ \

-  __asm  movq mm2,[Y+_r3] \

-  __asm  movq mm7,mm0 \

-  __asm  pxor mm0,mm4 \

-  __asm  pand mm7,mm4 \

-  __asm  psraw mm0,1 \

-  __asm  mov A,0x7FFF54DC \

-  __asm  paddw mm0,mm7 \

-  __asm  movd mm7,A \

-  /*mm7={54491-0x7FFF,0x7FFF}x2 \

-    mm4=_y[4]=v=r-u*/ \

-  __asm  psubw mm4,mm0 \

-  __asm  punpckldq mm7,mm7 \

-  __asm  movq [Y+_r4],mm4 \

-  /*mm0=0, mm7={36410}x4 \

-    mm1=(t5'''!=0), mm5:mm4=54491*t5'''+0x0E3D*/ \

-  __asm  movq mm4,mm1 \

-  __asm  movq mm5,mm1 \

-  __asm  punpcklwd mm4,mm1 \

-  __asm  mov A,0x8E3A8E3A \

-  __asm  pmaddwd mm4,mm7 \

-  __asm  movq [Y+_r0],mm0 \

-  __asm  punpckhwd mm5,mm1 \

-  __asm  pxor mm0,mm0 \

-  __asm  pmaddwd mm5,mm7 \

-  __asm  pcmpeqw mm1,mm0 \

-  __asm  movd mm7,A \

-  __asm  psubw mm1,mm3 \

-  __asm  punpckldq mm7,mm7 \

-  __asm  paddd mm4,mm6 \

-  __asm  paddd mm5,mm6 \

-  /*mm0=0 \

-    mm3:mm1=36410*t6'''+((t5'''!=0)<<16)*/ \

-  __asm  movq mm6,mm2 \

-  __asm  movq mm3,mm2 \

-  __asm  pmulhw mm6,mm7 \

-  __asm  paddw mm1,mm2 \

-  __asm  pmullw mm3,mm7 \

-  __asm  pxor mm0,mm0 \

-  __asm  paddw mm6,mm1 \

-  __asm  movq mm1,mm3 \

-  __asm  punpckhwd mm3,mm6 \

-  __asm  punpcklwd mm1,mm6 \

-  /*mm3={-1}x4, mm6={1}x4 \

-    mm4=_y[5]=u=(54491*t5'''+36410*t6'''+0x0E3D>>16)+(t5'''!=0)*/ \

-  __asm  paddd mm5,mm3 \

-  __asm  paddd mm4,mm1 \

-  __asm  psrad mm5,16 \

-  __asm  pxor mm6,mm6 \

-  __asm  psrad mm4,16 \

-  __asm  pcmpeqb mm3,mm3 \

-  __asm  packssdw mm4,mm5 \

-  __asm  psubw mm6,mm3 \

-  /*mm1=t7'', mm7={26568,0x3400}x2 \

-    mm2=s=t6'''-(36410*u>>16)*/ \

-  __asm  movq mm1,mm4 \

-  __asm  mov A,0x340067C8 \

-  __asm  pmulhw mm4,mm7 \

-  __asm  movd mm7,A \

-  __asm  movq [Y+_r5],mm1 \

-  __asm  punpckldq mm7,mm7 \

-  __asm  paddw mm4,mm1 \

-  __asm  movq mm1,[Y+_r7] \

-  __asm  psubw mm2,mm4 \

-  /*mm6={0x00007B1B}x2 \

-    mm0=(s!=0), mm5:mm4=s*26568+0x3400*/ \

-  __asm  movq mm4,mm2 \

-  __asm  movq mm5,mm2 \

-  __asm  punpcklwd mm4,mm6 \

-  __asm  pcmpeqw mm0,mm2 \

-  __asm  pmaddwd mm4,mm7 \

-  __asm  mov A,0x7B1B \

-  __asm  punpckhwd mm5,mm6 \

-  __asm  movd mm6,A \

-  __asm  pmaddwd mm5,mm7 \

-  __asm  psubw mm0,mm3 \

-  __asm  punpckldq mm6,mm6 \

-  /*mm7={64277-0x7FFF,0x7FFF}x2 \

-    mm2=_y[3]=v=(s*26568+0x3400>>17)+s+(s!=0)*/ \

-  __asm  psrad mm4,17 \

-  __asm  paddw mm2,mm0 \

-  __asm  psrad mm5,17 \

-  __asm  mov A,0x7FFF7B16 \

-  __asm  packssdw mm4,mm5 \

-  __asm  movd mm7,A \

-  __asm  paddw mm2,mm4 \

-  __asm  punpckldq mm7,mm7 \

-  /*mm0=0, mm7={12785}x4 \

-    mm1=(t7''!=0), mm2=t4'', mm5:mm4=64277*t7''+0x7B1B*/ \

-  __asm  movq mm4,mm1 \

-  __asm  movq mm5,mm1 \

-  __asm  movq [Y+_r3],mm2 \

-  __asm  punpcklwd mm4,mm1 \

-  __asm  movq mm2,[Y+_r1] \

-  __asm  pmaddwd mm4,mm7 \

-  __asm  mov A,0x31F131F1 \

-  __asm  punpckhwd mm5,mm1 \

-  __asm  pxor mm0,mm0 \

-  __asm  pmaddwd mm5,mm7 \

-  __asm  pcmpeqw mm1,mm0 \

-  __asm  movd mm7,A \

-  __asm  psubw mm1,mm3 \

-  __asm  punpckldq mm7,mm7 \

-  __asm  paddd mm4,mm6 \

-  __asm  paddd mm5,mm6 \

-  /*mm3:mm1=12785*t4'''+((t7''!=0)<<16)*/ \

-  __asm  movq mm6,mm2 \

-  __asm  movq mm3,mm2 \

-  __asm  pmulhw mm6,mm7 \

-  __asm  pmullw mm3,mm7 \

-  __asm  paddw mm6,mm1 \

-  __asm  movq mm1,mm3 \

-  __asm  punpckhwd mm3,mm6 \

-  __asm  punpcklwd mm1,mm6 \

-  /*mm3={-1}x4, mm6={1}x4 \

-    mm4=_y[1]=u=(12785*t4'''+64277*t7''+0x7B1B>>16)+(t7''!=0)*/ \

-  __asm  paddd mm5,mm3 \

-  __asm  paddd mm4,mm1 \

-  __asm  psrad mm5,16 \

-  __asm  pxor mm6,mm6 \

-  __asm  psrad mm4,16 \

-  __asm  pcmpeqb mm3,mm3 \

-  __asm  packssdw mm4,mm5 \

-  __asm  psubw mm6,mm3 \

-  /*mm1=t3'', mm7={20539,0x3000}x2 \

-    mm4=s=(12785*u>>16)-t4''*/ \

-  __asm  movq [Y+_r1],mm4 \

-  __asm  pmulhw mm4,mm7 \

-  __asm  mov A,0x3000503B \

-  __asm  movq mm1,[Y+_r6] \

-  __asm  movd mm7,A \

-  __asm  psubw mm4,mm2 \

-  __asm  punpckldq mm7,mm7 \

-  /*mm6={0x00006CB7}x2 \

-    mm0=(s!=0), mm5:mm4=s*20539+0x3000*/ \

-  __asm  movq mm5,mm4 \

-  __asm  movq mm2,mm4 \

-  __asm  punpcklwd mm4,mm6 \

-  __asm  pcmpeqw mm0,mm2 \

-  __asm  pmaddwd mm4,mm7 \

-  __asm  mov A,0x6CB7 \

-  __asm  punpckhwd mm5,mm6 \

-  __asm  movd mm6,A \

-  __asm  pmaddwd mm5,mm7 \

-  __asm  psubw mm0,mm3 \

-  __asm  punpckldq mm6,mm6 \

-  /*mm7={60547-0x7FFF,0x7FFF}x2 \

-    mm2=_y[7]=v=(s*20539+0x3000>>20)+s+(s!=0)*/ \

-  __asm  psrad mm4,20 \

-  __asm  paddw mm2,mm0 \

-  __asm  psrad mm5,20 \

-  __asm  mov A,0x7FFF6C84 \

-  __asm  packssdw mm4,mm5 \

-  __asm  movd mm7,A \

-  __asm  paddw mm2,mm4 \

-  __asm  punpckldq mm7,mm7 \

-  /*mm0=0, mm7={25080}x4 \

-    mm2=t2'', mm5:mm4=60547*t3''+0x6CB7*/ \

-  __asm  movq mm4,mm1 \

-  __asm  movq mm5,mm1 \

-  __asm  movq [Y+_r7],mm2 \

-  __asm  punpcklwd mm4,mm1 \

-  __asm  movq mm2,[Y+_r2] \

-  __asm  pmaddwd mm4,mm7 \

-  __asm  mov A,0x61F861F8 \

-  __asm  punpckhwd mm5,mm1 \

-  __asm  pxor mm0,mm0 \

-  __asm  pmaddwd mm5,mm7 \

-  __asm  movd mm7,A \

-  __asm  pcmpeqw mm1,mm0 \

-  __asm  psubw mm1,mm3 \

-  __asm  punpckldq mm7,mm7 \

-  __asm  paddd mm4,mm6 \

-  __asm  paddd mm5,mm6 \

-  /*mm3:mm1=25080*t2''+((t3''!=0)<<16)*/ \

-  __asm  movq mm6,mm2 \

-  __asm  movq mm3,mm2 \

-  __asm  pmulhw mm6,mm7 \

-  __asm  pmullw mm3,mm7 \

-  __asm  paddw mm6,mm1 \

-  __asm  movq mm1,mm3 \

-  __asm  punpckhwd mm3,mm6 \

-  __asm  punpcklwd mm1,mm6 \

-  /*mm1={-1}x4 \

-    mm4=u=(25080*t2''+60547*t3''+0x6CB7>>16)+(t3''!=0)*/ \

-  __asm  paddd mm5,mm3 \

-  __asm  paddd mm4,mm1 \

-  __asm  psrad mm5,16 \

-  __asm  mov A,0x28005460 \

-  __asm  psrad mm4,16 \

-  __asm  pcmpeqb mm1,mm1 \

-  __asm  packssdw mm4,mm5 \

-  /*mm5={1}x4, mm6=_y[2]=u, mm7={21600,0x2800}x2 \

-    mm4=s=(25080*u>>16)-t2''*/ \

-  __asm  movq mm6,mm4 \

-  __asm  pmulhw mm4,mm7 \

-  __asm  pxor mm5,mm5 \

-  __asm  movd mm7,A \

-  __asm  psubw mm5,mm1 \

-  __asm  punpckldq mm7,mm7 \

-  __asm  psubw mm4,mm2 \

-  /*mm2=s+(s!=0) \

-    mm4:mm3=s*21600+0x2800*/ \

-  __asm  movq mm3,mm4 \

-  __asm  movq mm2,mm4 \

-  __asm  punpckhwd mm4,mm5 \

-  __asm  pcmpeqw mm0,mm2 \

-  __asm  pmaddwd mm4,mm7 \

-  __asm  psubw mm0,mm1 \

-  __asm  punpcklwd mm3,mm5 \

-  __asm  paddw mm2,mm0 \

-  __asm  pmaddwd mm3,mm7 \

-  /*mm0=_y[4], mm1=_y[7], mm4=_y[0], mm5=_y[5] \

-    mm3=_y[6]=v=(s*21600+0x2800>>18)+s+(s!=0)*/ \

-  __asm  movq mm0,[Y+_r4] \

-  __asm  psrad mm4,18 \

-  __asm  movq mm5,[Y+_r5] \

-  __asm  psrad mm3,18 \

-  __asm  movq mm1,[Y+_r7] \

-  __asm  packssdw mm3,mm4 \

-  __asm  movq mm4,[Y+_r0] \

-  __asm  paddw mm3,mm2 \

-}

-

-/*On input, mm4=_y[0], mm6=_y[2], mm0=_y[4], mm5=_y[5], mm3=_y[6], mm1=_y[7].

-  On output, {_y[4],mm1,mm2,mm3} contains the transpose of _y[4...7] and

-   {mm4,mm5,mm6,mm7} contains the transpose of _y[0...3].*/

-#define OC_TRANSPOSE8x4(_r0,_r1,_r2,_r3,_r4,_r5,_r6,_r7) __asm{ \

-  /*First 4x4 transpose:*/ \

-  /*mm0 = e3 e2 e1 e0 \

-    mm5 = f3 f2 f1 f0 \

-    mm3 = g3 g2 g1 g0 \

-    mm1 = h3 h2 h1 h0*/ \

-  __asm  movq mm2,mm0 \

-  __asm  punpcklwd mm0,mm5 \

-  __asm  punpckhwd mm2,mm5 \

-  __asm  movq mm5,mm3 \

-  __asm  punpcklwd mm3,mm1 \

-  __asm  punpckhwd mm5,mm1 \

-  /*mm0 = f1 e1 f0 e0 \

-    mm2 = f3 e3 f2 e2 \

-    mm3 = h1 g1 h0 g0 \

-    mm5 = h3 g3 h2 g2*/ \

-  __asm  movq mm1,mm0 \

-  __asm  punpckldq mm0,mm3 \

-  __asm  movq [Y+_r4],mm0 \

-  __asm  punpckhdq mm1,mm3 \

-  __asm  movq mm0,[Y+_r1] \

-  __asm  movq mm3,mm2 \

-  __asm  punpckldq mm2,mm5 \

-  __asm  punpckhdq mm3,mm5 \

-  __asm  movq mm5,[Y+_r3] \

-  /*_y[4] = h0 g0 f0 e0 \

-   mm1  = h1 g1 f1 e1 \

-   mm2  = h2 g2 f2 e2 \

-   mm3  = h3 g3 f3 e3*/ \

-  /*Second 4x4 transpose:*/ \

-  /*mm4 = a3 a2 a1 a0 \

-    mm0 = b3 b2 b1 b0 \

-    mm6 = c3 c2 c1 c0 \

-    mm5 = d3 d2 d1 d0*/ \

-  __asm  movq mm7,mm4 \

-  __asm  punpcklwd mm4,mm0 \

-  __asm  punpckhwd mm7,mm0 \

-  __asm  movq mm0,mm6 \

-  __asm  punpcklwd mm6,mm5 \

-  __asm  punpckhwd mm0,mm5 \

-  /*mm4 = b1 a1 b0 a0 \

-    mm7 = b3 a3 b2 a2 \

-    mm6 = d1 c1 d0 c0 \

-    mm0 = d3 c3 d2 c2*/ \

-  __asm  movq mm5,mm4 \

-  __asm  punpckldq mm4,mm6 \

-  __asm  punpckhdq mm5,mm6 \

-  __asm  movq mm6,mm7 \

-  __asm  punpckhdq mm7,mm0 \

-  __asm  punpckldq mm6,mm0 \

-  /*mm4 = d0 c0 b0 a0 \

-    mm5 = d1 c1 b1 a1 \

-    mm6 = d2 c2 b2 a2 \

-    mm7 = d3 c3 b3 a3*/ \

-}

-

-/*MMX implementation of the fDCT.*/

-void oc_enc_fdct8x8_mmx(ogg_int16_t _y[64],const ogg_int16_t _x[64]){

-  ptrdiff_t a;

-  __asm{

-#define Y eax

-#define A ecx

-#define X edx

-    /*Add two extra bits of working precision to improve accuracy; any more and

-       we could overflow.*/

-    /*We also add biases to correct for some systematic error that remains in

-       the full fDCT->iDCT round trip.*/

-    mov X, _x

-    mov Y, _y

-    movq mm0,[0x00+X]

-    movq mm1,[0x10+X]

-    movq mm2,[0x20+X]

-    movq mm3,[0x30+X]

-    pcmpeqb mm4,mm4

-    pxor mm7,mm7

-    movq mm5,mm0

-    psllw mm0,2

-    pcmpeqw mm5,mm7

-    movq mm7,[0x70+X]

-    psllw mm1,2

-    psubw mm5,mm4

-    psllw mm2,2

-    mov A,1

-    pslld mm5,16

-    movd mm6,A

-    psllq mm5,16

-    mov A,0x10001

-    psllw mm3,2

-    movd mm4,A

-    punpckhwd mm5,mm6

-    psubw mm1,mm6

-    movq mm6,[0x60+X]

-    paddw mm0,mm5

-    movq mm5,[0x50+X]

-    paddw mm0,mm4

-    movq mm4,[0x40+X]

-    /*We inline stage1 of the transform here so we can get better instruction

-       scheduling with the shifts.*/

-    /*mm0=t7'=t0-t7*/

-    psllw mm7,2

-    psubw mm0,mm7

-    psllw mm6,2

-    paddw mm7,mm7

-    /*mm1=t6'=t1-t6*/

-    psllw mm5,2

-    psubw mm1,mm6

-    psllw mm4,2

-    paddw mm6,mm6

-    /*mm2=t5'=t2-t5*/

-    psubw mm2,mm5

-    paddw mm5,mm5

-    /*mm3=t4'=t3-t4*/

-    psubw mm3,mm4

-    paddw mm4,mm4

-    /*mm7=t0'=t0+t7*/

-    paddw mm7,mm0

-    /*mm6=t1'=t1+t6*/

-    paddw mm6,mm1

-    /*mm5=t2'=t2+t5*/

-    paddw mm5,mm2

-    /*mm4=t3'=t3+t4*/

-    paddw mm4,mm3

-    OC_FDCT8x4(0x00,0x10,0x20,0x30,0x40,0x50,0x60,0x70)

-    OC_TRANSPOSE8x4(0x00,0x10,0x20,0x30,0x40,0x50,0x60,0x70)

-    /*Swap out this 8x4 block for the next one.*/

-    movq mm0,[0x08+X]

-    movq [0x30+Y],mm7

-    movq mm7,[0x78+X]

-    movq [0x50+Y],mm1

-    movq mm1,[0x18+X]

-    movq [0x20+Y],mm6

-    movq mm6,[0x68+X]

-    movq [0x60+Y],mm2

-    movq mm2,[0x28+X]

-    movq [0x10+Y],mm5

-    movq mm5,[0x58+X]

-    movq [0x70+Y],mm3

-    movq mm3,[0x38+X]

-    /*And increase its working precision, too.*/

-    psllw mm0,2

-    movq [0x00+Y],mm4

-    psllw mm7,2

-    movq mm4,[0x48+X]

-    /*We inline stage1 of the transform here so we can get better instruction

-       scheduling with the shifts.*/

-    /*mm0=t7'=t0-t7*/

-    psubw mm0,mm7

-    psllw mm1,2

-    paddw mm7,mm7

-    psllw mm6,2

-    /*mm1=t6'=t1-t6*/

-    psubw mm1,mm6

-    psllw mm2,2

-    paddw mm6,mm6

-    psllw mm5,2

-    /*mm2=t5'=t2-t5*/

-    psubw mm2,mm5

-    psllw mm3,2

-    paddw mm5,mm5

-    psllw mm4,2

-    /*mm3=t4'=t3-t4*/

-    psubw mm3,mm4

-    paddw mm4,mm4

-    /*mm7=t0'=t0+t7*/

-    paddw mm7,mm0

-    /*mm6=t1'=t1+t6*/

-    paddw mm6,mm1

-    /*mm5=t2'=t2+t5*/

-    paddw mm5,mm2

-    /*mm4=t3'=t3+t4*/

-    paddw mm4,mm3

-    OC_FDCT8x4(0x08,0x18,0x28,0x38,0x48,0x58,0x68,0x78)

-    OC_TRANSPOSE8x4(0x08,0x18,0x28,0x38,0x48,0x58,0x68,0x78)

-    /*Here the first 4x4 block of output from the last transpose is the second

-       4x4 block of input for the next transform.

-      We have cleverly arranged that it already be in the appropriate place,

-       so we only have to do half the stores and loads.*/

-    movq mm0,[0x00+Y]

-    movq [0x58+Y],mm1

-    movq mm1,[0x10+Y]

-    movq [0x68+Y],mm2

-    movq mm2,[0x20+Y]

-    movq [0x78+Y],mm3

-    movq mm3,[0x30+Y]

-    OC_FDCT_STAGE1_8x4

-    OC_FDCT8x4(0x00,0x10,0x20,0x30,0x08,0x18,0x28,0x38)

-    OC_TRANSPOSE8x4(0x00,0x10,0x20,0x30,0x08,0x18,0x28,0x38)

-    /*mm0={-2}x4*/

-    pcmpeqw mm0,mm0

-    paddw mm0,mm0

-    /*Round the results.*/

-    psubw mm1,mm0

-    psubw mm2,mm0

-    psraw mm1,2

-    psubw mm3,mm0

-    movq [0x18+Y],mm1

-    psraw mm2,2

-    psubw mm4,mm0

-    movq mm1,[0x08+Y]

-    psraw mm3,2

-    psubw mm5,mm0

-    psraw mm4,2

-    psubw mm6,mm0

-    psraw mm5,2

-    psubw mm7,mm0

-    psraw mm6,2

-    psubw mm1,mm0

-    psraw mm7,2

-    movq mm0,[0x40+Y]

-    psraw mm1,2

-    movq [0x30+Y],mm7

-    movq mm7,[0x78+Y]

-    movq [0x08+Y],mm1

-    movq mm1,[0x50+Y]

-    movq [0x20+Y],mm6

-    movq mm6,[0x68+Y]

-    movq [0x28+Y],mm2

-    movq mm2,[0x60+Y]

-    movq [0x10+Y],mm5

-    movq mm5,[0x58+Y]

-    movq [0x38+Y],mm3

-    movq mm3,[0x70+Y]

-    movq [0x00+Y],mm4

-    movq mm4,[0x48+Y]

-    OC_FDCT_STAGE1_8x4

-    OC_FDCT8x4(0x40,0x50,0x60,0x70,0x48,0x58,0x68,0x78)

-    OC_TRANSPOSE8x4(0x40,0x50,0x60,0x70,0x48,0x58,0x68,0x78)

-    /*mm0={-2}x4*/

-    pcmpeqw mm0,mm0

-    paddw mm0,mm0

-    /*Round the results.*/

-    psubw mm1,mm0

-    psubw mm2,mm0

-    psraw mm1,2

-    psubw mm3,mm0

-    movq [0x58+Y],mm1

-    psraw mm2,2

-    psubw mm4,mm0

-    movq mm1,[0x48+Y]

-    psraw mm3,2

-    psubw mm5,mm0

-    movq [0x68+Y],mm2

-    psraw mm4,2

-    psubw mm6,mm0

-    movq [0x78+Y],mm3

-    psraw mm5,2

-    psubw mm7,mm0

-    movq [0x40+Y],mm4

-    psraw mm6,2

-    psubw mm1,mm0

-    movq [0x50+Y],mm5

-    psraw mm7,2

-    movq [0x60+Y],mm6

-    psraw mm1,2

-    movq [0x70+Y],mm7

-    movq [0x48+Y],mm1

-#undef Y

-#undef A

-#undef X

-  }

-}

-

-#endif

+/********************************************************************
+ *                                                                  *
+ * THIS FILE IS PART OF THE OggTheora SOFTWARE CODEC SOURCE CODE.   *
+ * USE, DISTRIBUTION AND REPRODUCTION OF THIS LIBRARY SOURCE IS     *
+ * GOVERNED BY A BSD-STYLE SOURCE LICENSE INCLUDED WITH THIS SOURCE *
+ * IN 'COPYING'. PLEASE READ THESE TERMS BEFORE DISTRIBUTING.       *
+ *                                                                  *
+ * THE Theora SOURCE CODE IS COPYRIGHT (C) 1999-2006                *
+ * by the Xiph.Org Foundation http://www.xiph.org/                  *
+ *                                                                  *
+ ********************************************************************/ 
+ /*MMX fDCT implementation for x86_32*/
+/*$Id: fdct_ses2.c 14579 2008-03-12 06:42:40Z xiphmont $*/
+#include "x86enc.h"
+
+#if defined(OC_X86_ASM)
+
+#define OC_FDCT_STAGE1_8x4  __asm{ \
+  /*Stage 1:*/ \
+  /*mm0=t7'=t0-t7*/ \
+  __asm  psubw mm0,mm7 \
+  __asm  paddw mm7,mm7 \
+  /*mm1=t6'=t1-t6*/ \
+  __asm  psubw mm1, mm6 \
+  __asm  paddw mm6,mm6 \
+  /*mm2=t5'=t2-t5*/ \
+  __asm  psubw mm2,mm5 \
+  __asm  paddw mm5,mm5 \
+  /*mm3=t4'=t3-t4*/ \
+  __asm  psubw mm3,mm4 \
+  __asm  paddw mm4,mm4 \
+  /*mm7=t0'=t0+t7*/ \
+  __asm  paddw mm7,mm0 \
+  /*mm6=t1'=t1+t6*/  \
+  __asm  paddw mm6,mm1 \
+  /*mm5=t2'=t2+t5*/ \
+  __asm  paddw mm5,mm2 \
+  /*mm4=t3'=t3+t4*/ \
+  __asm  paddw mm4,mm3\
+}
+
+#define OC_FDCT8x4(_r0,_r1,_r2,_r3,_r4,_r5,_r6,_r7) __asm{ \
+  /*Stage 2:*/ \
+  /*mm7=t3''=t0'-t3'*/ \
+  __asm  psubw mm7,mm4 \
+  __asm  paddw mm4,mm4 \
+  /*mm6=t2''=t1'-t2'*/ \
+  __asm  psubw mm6,mm5 \
+  __asm  movq [Y+_r6],mm7 \
+  __asm  paddw mm5,mm5 \
+  /*mm1=t5''=t6'-t5'*/ \
+  __asm  psubw mm1,mm2 \
+  __asm  movq [Y+_r2],mm6 \
+  /*mm4=t0''=t0'+t3'*/ \
+  __asm  paddw mm4,mm7 \
+  __asm  paddw mm2,mm2 \
+  /*mm5=t1''=t1'+t2'*/ \
+  __asm  movq [Y+_r0],mm4 \
+  __asm  paddw mm5,mm6 \
+  /*mm2=t6''=t6'+t5'*/ \
+  __asm  paddw mm2,mm1 \
+  __asm  movq [Y+_r4],mm5 \
+  /*mm0=t7', mm1=t5'', mm2=t6'', mm3=t4'.*/ \
+  /*mm4, mm5, mm6, mm7 are free.*/ \
+  /*Stage 3:*/ \
+  /*mm6={2}x4, mm7={27146,0xB500>>1}x2*/ \
+  __asm  mov A,0x5A806A0A \
+  __asm  pcmpeqb mm6,mm6 \
+  __asm  movd mm7,A \
+  __asm  psrlw mm6,15 \
+  __asm  punpckldq mm7,mm7 \
+  __asm  paddw mm6,mm6 \
+  /*mm0=0, m2={-1}x4 \
+    mm5:mm4=t5''*27146+0xB500*/ \
+  __asm  movq mm4,mm1 \
+  __asm  movq mm5,mm1 \
+  __asm  punpcklwd mm4,mm6 \
+  __asm  movq [Y+_r3],mm2 \
+  __asm  pmaddwd mm4,mm7 \
+  __asm  movq [Y+_r7],mm0 \
+  __asm  punpckhwd mm5,mm6 \
+  __asm  pxor mm0,mm0 \
+  __asm  pmaddwd mm5,mm7 \
+  __asm  pcmpeqb mm2,mm2 \
+  /*mm2=t6'', mm1=t5''+(t5''!=0) \
+    mm4=(t5''*27146+0xB500>>16)*/ \
+  __asm  pcmpeqw mm0,mm1 \
+  __asm  psrad mm4,16 \
+  __asm  psubw mm0,mm2 \
+  __asm  movq mm2, [Y+_r3] \
+  __asm  psrad mm5,16 \
+  __asm  paddw mm1,mm0 \
+  __asm  packssdw mm4,mm5 \
+  /*mm4=s=(t5''*27146+0xB500>>16)+t5''+(t5''!=0)>>1*/ \
+  __asm  paddw mm4,mm1 \
+  __asm  movq mm0, [Y+_r7] \
+  __asm  psraw mm4,1 \
+  __asm  movq mm1,mm3 \
+  /*mm3=t4''=t4'+s*/ \
+  __asm  paddw mm3,mm4 \
+  /*mm1=t5'''=t4'-s*/ \
+  __asm  psubw mm1,mm4 \
+  /*mm1=0, mm3={-1}x4 \
+    mm5:mm4=t6''*27146+0xB500*/ \
+  __asm  movq mm4,mm2 \
+  __asm  movq mm5,mm2 \
+  __asm  punpcklwd mm4,mm6 \
+  __asm  movq [Y+_r5],mm1 \
+  __asm  pmaddwd mm4,mm7 \
+  __asm  movq [Y+_r1],mm3 \
+  __asm  punpckhwd mm5,mm6 \
+  __asm  pxor mm1,mm1 \
+  __asm  pmaddwd mm5,mm7 \
+  __asm  pcmpeqb mm3,mm3 \
+  /*mm2=t6''+(t6''!=0), mm4=(t6''*27146+0xB500>>16)*/ \
+  __asm  psrad mm4,16 \
+  __asm  pcmpeqw mm1,mm2 \
+  __asm  psrad mm5,16 \
+  __asm  psubw mm1,mm3 \
+  __asm  packssdw mm4,mm5 \
+  __asm  paddw mm2,mm1 \
+  /*mm1=t1'' \
+    mm4=s=(t6''*27146+0xB500>>16)+t6''+(t6''!=0)>>1*/ \
+  __asm  paddw mm4,mm2 \
+  __asm  movq mm1,[Y+_r4] \
+  __asm  psraw mm4,1 \
+  __asm  movq mm2,mm0 \
+  /*mm7={54491-0x7FFF,0x7FFF}x2 \
+    mm0=t7''=t7'+s*/ \
+  __asm  paddw mm0,mm4 \
+  /*mm2=t6'''=t7'-s*/ \
+  __asm  psubw mm2,mm4 \
+  /*Stage 4:*/ \
+  /*mm0=0, mm2=t0'' \
+    mm5:mm4=t1''*27146+0xB500*/ \
+  __asm  movq mm4,mm1 \
+  __asm  movq mm5,mm1 \
+  __asm  punpcklwd mm4,mm6 \
+  __asm  movq [Y+_r3],mm2 \
+  __asm  pmaddwd mm4,mm7 \
+  __asm  movq mm2,[Y+_r0] \
+  __asm  punpckhwd mm5,mm6 \
+  __asm  movq [Y+_r7],mm0 \
+  __asm  pmaddwd mm5,mm7 \
+  __asm  pxor mm0,mm0 \
+  /*mm7={27146,0x4000>>1}x2 \
+    mm0=s=(t1''*27146+0xB500>>16)+t1''+(t1''!=0)*/ \
+  __asm  psrad mm4,16 \
+  __asm  mov A,0x20006A0A \
+  __asm  pcmpeqw mm0,mm1 \
+  __asm  movd mm7,A \
+  __asm  psrad mm5,16 \
+  __asm  psubw mm0,mm3 \
+  __asm  packssdw mm4,mm5 \
+  __asm  paddw mm0,mm1 \
+  __asm  punpckldq mm7,mm7 \
+  __asm  paddw mm0,mm4 \
+  /*mm6={0x00000E3D}x2 \
+    mm1=-(t0''==0), mm5:mm4=t0''*27146+0x4000*/ \
+  __asm  movq mm4,mm2 \
+  __asm  movq mm5,mm2 \
+  __asm  punpcklwd mm4,mm6 \
+  __asm  mov A,0x0E3D \
+  __asm  pmaddwd mm4,mm7 \
+  __asm  punpckhwd mm5,mm6 \
+  __asm  movd mm6,A \
+  __asm  pmaddwd mm5,mm7 \
+  __asm  pxor mm1,mm1 \
+  __asm  punpckldq mm6,mm6 \
+  __asm  pcmpeqw mm1,mm2 \
+  /*mm4=r=(t0''*27146+0x4000>>16)+t0''+(t0''!=0)*/ \
+  __asm  psrad mm4,16 \
+  __asm  psubw mm1,mm3 \
+  __asm  psrad mm5,16 \
+  __asm  paddw mm2,mm1 \
+  __asm  packssdw mm4,mm5 \
+  __asm  movq mm1,[Y+_r5] \
+  __asm  paddw mm4,mm2 \
+  /*mm2=t6'', mm0=_y[0]=u=r+s>>1 \
+    The naive implementation could cause overflow, so we use \
+     u=(r&s)+((r^s)>>1).*/ \
+  __asm  movq mm2,[Y+_r3] \
+  __asm  movq mm7,mm0 \
+  __asm  pxor mm0,mm4 \
+  __asm  pand mm7,mm4 \
+  __asm  psraw mm0,1 \
+  __asm  mov A,0x7FFF54DC \
+  __asm  paddw mm0,mm7 \
+  __asm  movd mm7,A \
+  /*mm7={54491-0x7FFF,0x7FFF}x2 \
+    mm4=_y[4]=v=r-u*/ \
+  __asm  psubw mm4,mm0 \
+  __asm  punpckldq mm7,mm7 \
+  __asm  movq [Y+_r4],mm4 \
+  /*mm0=0, mm7={36410}x4 \
+    mm1=(t5'''!=0), mm5:mm4=54491*t5'''+0x0E3D*/ \
+  __asm  movq mm4,mm1 \
+  __asm  movq mm5,mm1 \
+  __asm  punpcklwd mm4,mm1 \
+  __asm  mov A,0x8E3A8E3A \
+  __asm  pmaddwd mm4,mm7 \
+  __asm  movq [Y+_r0],mm0 \
+  __asm  punpckhwd mm5,mm1 \
+  __asm  pxor mm0,mm0 \
+  __asm  pmaddwd mm5,mm7 \
+  __asm  pcmpeqw mm1,mm0 \
+  __asm  movd mm7,A \
+  __asm  psubw mm1,mm3 \
+  __asm  punpckldq mm7,mm7 \
+  __asm  paddd mm4,mm6 \
+  __asm  paddd mm5,mm6 \
+  /*mm0=0 \
+    mm3:mm1=36410*t6'''+((t5'''!=0)<<16)*/ \
+  __asm  movq mm6,mm2 \
+  __asm  movq mm3,mm2 \
+  __asm  pmulhw mm6,mm7 \
+  __asm  paddw mm1,mm2 \
+  __asm  pmullw mm3,mm7 \
+  __asm  pxor mm0,mm0 \
+  __asm  paddw mm6,mm1 \
+  __asm  movq mm1,mm3 \
+  __asm  punpckhwd mm3,mm6 \
+  __asm  punpcklwd mm1,mm6 \
+  /*mm3={-1}x4, mm6={1}x4 \
+    mm4=_y[5]=u=(54491*t5'''+36410*t6'''+0x0E3D>>16)+(t5'''!=0)*/ \
+  __asm  paddd mm5,mm3 \
+  __asm  paddd mm4,mm1 \
+  __asm  psrad mm5,16 \
+  __asm  pxor mm6,mm6 \
+  __asm  psrad mm4,16 \
+  __asm  pcmpeqb mm3,mm3 \
+  __asm  packssdw mm4,mm5 \
+  __asm  psubw mm6,mm3 \
+  /*mm1=t7'', mm7={26568,0x3400}x2 \
+    mm2=s=t6'''-(36410*u>>16)*/ \
+  __asm  movq mm1,mm4 \
+  __asm  mov A,0x340067C8 \
+  __asm  pmulhw mm4,mm7 \
+  __asm  movd mm7,A \
+  __asm  movq [Y+_r5],mm1 \
+  __asm  punpckldq mm7,mm7 \
+  __asm  paddw mm4,mm1 \
+  __asm  movq mm1,[Y+_r7] \
+  __asm  psubw mm2,mm4 \
+  /*mm6={0x00007B1B}x2 \
+    mm0=(s!=0), mm5:mm4=s*26568+0x3400*/ \
+  __asm  movq mm4,mm2 \
+  __asm  movq mm5,mm2 \
+  __asm  punpcklwd mm4,mm6 \
+  __asm  pcmpeqw mm0,mm2 \
+  __asm  pmaddwd mm4,mm7 \
+  __asm  mov A,0x7B1B \
+  __asm  punpckhwd mm5,mm6 \
+  __asm  movd mm6,A \
+  __asm  pmaddwd mm5,mm7 \
+  __asm  psubw mm0,mm3 \
+  __asm  punpckldq mm6,mm6 \
+  /*mm7={64277-0x7FFF,0x7FFF}x2 \
+    mm2=_y[3]=v=(s*26568+0x3400>>17)+s+(s!=0)*/ \
+  __asm  psrad mm4,17 \
+  __asm  paddw mm2,mm0 \
+  __asm  psrad mm5,17 \
+  __asm  mov A,0x7FFF7B16 \
+  __asm  packssdw mm4,mm5 \
+  __asm  movd mm7,A \
+  __asm  paddw mm2,mm4 \
+  __asm  punpckldq mm7,mm7 \
+  /*mm0=0, mm7={12785}x4 \
+    mm1=(t7''!=0), mm2=t4'', mm5:mm4=64277*t7''+0x7B1B*/ \
+  __asm  movq mm4,mm1 \
+  __asm  movq mm5,mm1 \
+  __asm  movq [Y+_r3],mm2 \
+  __asm  punpcklwd mm4,mm1 \
+  __asm  movq mm2,[Y+_r1] \
+  __asm  pmaddwd mm4,mm7 \
+  __asm  mov A,0x31F131F1 \
+  __asm  punpckhwd mm5,mm1 \
+  __asm  pxor mm0,mm0 \
+  __asm  pmaddwd mm5,mm7 \
+  __asm  pcmpeqw mm1,mm0 \
+  __asm  movd mm7,A \
+  __asm  psubw mm1,mm3 \
+  __asm  punpckldq mm7,mm7 \
+  __asm  paddd mm4,mm6 \
+  __asm  paddd mm5,mm6 \
+  /*mm3:mm1=12785*t4'''+((t7''!=0)<<16)*/ \
+  __asm  movq mm6,mm2 \
+  __asm  movq mm3,mm2 \
+  __asm  pmulhw mm6,mm7 \
+  __asm  pmullw mm3,mm7 \
+  __asm  paddw mm6,mm1 \
+  __asm  movq mm1,mm3 \
+  __asm  punpckhwd mm3,mm6 \
+  __asm  punpcklwd mm1,mm6 \
+  /*mm3={-1}x4, mm6={1}x4 \
+    mm4=_y[1]=u=(12785*t4'''+64277*t7''+0x7B1B>>16)+(t7''!=0)*/ \
+  __asm  paddd mm5,mm3 \
+  __asm  paddd mm4,mm1 \
+  __asm  psrad mm5,16 \
+  __asm  pxor mm6,mm6 \
+  __asm  psrad mm4,16 \
+  __asm  pcmpeqb mm3,mm3 \
+  __asm  packssdw mm4,mm5 \
+  __asm  psubw mm6,mm3 \
+  /*mm1=t3'', mm7={20539,0x3000}x2 \
+    mm4=s=(12785*u>>16)-t4''*/ \
+  __asm  movq [Y+_r1],mm4 \
+  __asm  pmulhw mm4,mm7 \
+  __asm  mov A,0x3000503B \
+  __asm  movq mm1,[Y+_r6] \
+  __asm  movd mm7,A \
+  __asm  psubw mm4,mm2 \
+  __asm  punpckldq mm7,mm7 \
+  /*mm6={0x00006CB7}x2 \
+    mm0=(s!=0), mm5:mm4=s*20539+0x3000*/ \
+  __asm  movq mm5,mm4 \
+  __asm  movq mm2,mm4 \
+  __asm  punpcklwd mm4,mm6 \
+  __asm  pcmpeqw mm0,mm2 \
+  __asm  pmaddwd mm4,mm7 \
+  __asm  mov A,0x6CB7 \
+  __asm  punpckhwd mm5,mm6 \
+  __asm  movd mm6,A \
+  __asm  pmaddwd mm5,mm7 \
+  __asm  psubw mm0,mm3 \
+  __asm  punpckldq mm6,mm6 \
+  /*mm7={60547-0x7FFF,0x7FFF}x2 \
+    mm2=_y[7]=v=(s*20539+0x3000>>20)+s+(s!=0)*/ \
+  __asm  psrad mm4,20 \
+  __asm  paddw mm2,mm0 \
+  __asm  psrad mm5,20 \
+  __asm  mov A,0x7FFF6C84 \
+  __asm  packssdw mm4,mm5 \
+  __asm  movd mm7,A \
+  __asm  paddw mm2,mm4 \
+  __asm  punpckldq mm7,mm7 \
+  /*mm0=0, mm7={25080}x4 \
+    mm2=t2'', mm5:mm4=60547*t3''+0x6CB7*/ \
+  __asm  movq mm4,mm1 \
+  __asm  movq mm5,mm1 \
+  __asm  movq [Y+_r7],mm2 \
+  __asm  punpcklwd mm4,mm1 \
+  __asm  movq mm2,[Y+_r2] \
+  __asm  pmaddwd mm4,mm7 \
+  __asm  mov A,0x61F861F8 \
+  __asm  punpckhwd mm5,mm1 \
+  __asm  pxor mm0,mm0 \
+  __asm  pmaddwd mm5,mm7 \
+  __asm  movd mm7,A \
+  __asm  pcmpeqw mm1,mm0 \
+  __asm  psubw mm1,mm3 \
+  __asm  punpckldq mm7,mm7 \
+  __asm  paddd mm4,mm6 \
+  __asm  paddd mm5,mm6 \
+  /*mm3:mm1=25080*t2''+((t3''!=0)<<16)*/ \
+  __asm  movq mm6,mm2 \
+  __asm  movq mm3,mm2 \
+  __asm  pmulhw mm6,mm7 \
+  __asm  pmullw mm3,mm7 \
+  __asm  paddw mm6,mm1 \
+  __asm  movq mm1,mm3 \
+  __asm  punpckhwd mm3,mm6 \
+  __asm  punpcklwd mm1,mm6 \
+  /*mm1={-1}x4 \
+    mm4=u=(25080*t2''+60547*t3''+0x6CB7>>16)+(t3''!=0)*/ \
+  __asm  paddd mm5,mm3 \
+  __asm  paddd mm4,mm1 \
+  __asm  psrad mm5,16 \
+  __asm  mov A,0x28005460 \
+  __asm  psrad mm4,16 \
+  __asm  pcmpeqb mm1,mm1 \
+  __asm  packssdw mm4,mm5 \
+  /*mm5={1}x4, mm6=_y[2]=u, mm7={21600,0x2800}x2 \
+    mm4=s=(25080*u>>16)-t2''*/ \
+  __asm  movq mm6,mm4 \
+  __asm  pmulhw mm4,mm7 \
+  __asm  pxor mm5,mm5 \
+  __asm  movd mm7,A \
+  __asm  psubw mm5,mm1 \
+  __asm  punpckldq mm7,mm7 \
+  __asm  psubw mm4,mm2 \
+  /*mm2=s+(s!=0) \
+    mm4:mm3=s*21600+0x2800*/ \
+  __asm  movq mm3,mm4 \
+  __asm  movq mm2,mm4 \
+  __asm  punpckhwd mm4,mm5 \
+  __asm  pcmpeqw mm0,mm2 \
+  __asm  pmaddwd mm4,mm7 \
+  __asm  psubw mm0,mm1 \
+  __asm  punpcklwd mm3,mm5 \
+  __asm  paddw mm2,mm0 \
+  __asm  pmaddwd mm3,mm7 \
+  /*mm0=_y[4], mm1=_y[7], mm4=_y[0], mm5=_y[5] \
+    mm3=_y[6]=v=(s*21600+0x2800>>18)+s+(s!=0)*/ \
+  __asm  movq mm0,[Y+_r4] \
+  __asm  psrad mm4,18 \
+  __asm  movq mm5,[Y+_r5] \
+  __asm  psrad mm3,18 \
+  __asm  movq mm1,[Y+_r7] \
+  __asm  packssdw mm3,mm4 \
+  __asm  movq mm4,[Y+_r0] \
+  __asm  paddw mm3,mm2 \
+}
+
+/*On input, mm4=_y[0], mm6=_y[2], mm0=_y[4], mm5=_y[5], mm3=_y[6], mm1=_y[7].
+  On output, {_y[4],mm1,mm2,mm3} contains the transpose of _y[4...7] and
+   {mm4,mm5,mm6,mm7} contains the transpose of _y[0...3].*/
+#define OC_TRANSPOSE8x4(_r0,_r1,_r2,_r3,_r4,_r5,_r6,_r7) __asm{ \
+  /*First 4x4 transpose:*/ \
+  /*mm0 = e3 e2 e1 e0 \
+    mm5 = f3 f2 f1 f0 \
+    mm3 = g3 g2 g1 g0 \
+    mm1 = h3 h2 h1 h0*/ \
+  __asm  movq mm2,mm0 \
+  __asm  punpcklwd mm0,mm5 \
+  __asm  punpckhwd mm2,mm5 \
+  __asm  movq mm5,mm3 \
+  __asm  punpcklwd mm3,mm1 \
+  __asm  punpckhwd mm5,mm1 \
+  /*mm0 = f1 e1 f0 e0 \
+    mm2 = f3 e3 f2 e2 \
+    mm3 = h1 g1 h0 g0 \
+    mm5 = h3 g3 h2 g2*/ \
+  __asm  movq mm1,mm0 \
+  __asm  punpckldq mm0,mm3 \
+  __asm  movq [Y+_r4],mm0 \
+  __asm  punpckhdq mm1,mm3 \
+  __asm  movq mm0,[Y+_r1] \
+  __asm  movq mm3,mm2 \
+  __asm  punpckldq mm2,mm5 \
+  __asm  punpckhdq mm3,mm5 \
+  __asm  movq mm5,[Y+_r3] \
+  /*_y[4] = h0 g0 f0 e0 \
+   mm1  = h1 g1 f1 e1 \
+   mm2  = h2 g2 f2 e2 \
+   mm3  = h3 g3 f3 e3*/ \
+  /*Second 4x4 transpose:*/ \
+  /*mm4 = a3 a2 a1 a0 \
+    mm0 = b3 b2 b1 b0 \
+    mm6 = c3 c2 c1 c0 \
+    mm5 = d3 d2 d1 d0*/ \
+  __asm  movq mm7,mm4 \
+  __asm  punpcklwd mm4,mm0 \
+  __asm  punpckhwd mm7,mm0 \
+  __asm  movq mm0,mm6 \
+  __asm  punpcklwd mm6,mm5 \
+  __asm  punpckhwd mm0,mm5 \
+  /*mm4 = b1 a1 b0 a0 \
+    mm7 = b3 a3 b2 a2 \
+    mm6 = d1 c1 d0 c0 \
+    mm0 = d3 c3 d2 c2*/ \
+  __asm  movq mm5,mm4 \
+  __asm  punpckldq mm4,mm6 \
+  __asm  punpckhdq mm5,mm6 \
+  __asm  movq mm6,mm7 \
+  __asm  punpckhdq mm7,mm0 \
+  __asm  punpckldq mm6,mm0 \
+  /*mm4 = d0 c0 b0 a0 \
+    mm5 = d1 c1 b1 a1 \
+    mm6 = d2 c2 b2 a2 \
+    mm7 = d3 c3 b3 a3*/ \
+}
+
+/*MMX implementation of the fDCT.*/
+void oc_enc_fdct8x8_mmx(ogg_int16_t _y[64],const ogg_int16_t _x[64]){
+  ptrdiff_t a;
+  __asm{
+#define Y eax
+#define A ecx
+#define X edx
+    /*Add two extra bits of working precision to improve accuracy; any more and
+       we could overflow.*/
+    /*We also add biases to correct for some systematic error that remains in
+       the full fDCT->iDCT round trip.*/
+    mov X, _x
+    mov Y, _y
+    movq mm0,[0x00+X]
+    movq mm1,[0x10+X]
+    movq mm2,[0x20+X]
+    movq mm3,[0x30+X]
+    pcmpeqb mm4,mm4
+    pxor mm7,mm7
+    movq mm5,mm0
+    psllw mm0,2
+    pcmpeqw mm5,mm7
+    movq mm7,[0x70+X]
+    psllw mm1,2
+    psubw mm5,mm4
+    psllw mm2,2
+    mov A,1
+    pslld mm5,16
+    movd mm6,A
+    psllq mm5,16
+    mov A,0x10001
+    psllw mm3,2
+    movd mm4,A
+    punpckhwd mm5,mm6
+    psubw mm1,mm6
+    movq mm6,[0x60+X]
+    paddw mm0,mm5
+    movq mm5,[0x50+X]
+    paddw mm0,mm4
+    movq mm4,[0x40+X]
+    /*We inline stage1 of the transform here so we can get better instruction
+       scheduling with the shifts.*/
+    /*mm0=t7'=t0-t7*/
+    psllw mm7,2
+    psubw mm0,mm7
+    psllw mm6,2
+    paddw mm7,mm7
+    /*mm1=t6'=t1-t6*/
+    psllw mm5,2
+    psubw mm1,mm6
+    psllw mm4,2
+    paddw mm6,mm6
+    /*mm2=t5'=t2-t5*/
+    psubw mm2,mm5
+    paddw mm5,mm5
+    /*mm3=t4'=t3-t4*/
+    psubw mm3,mm4
+    paddw mm4,mm4
+    /*mm7=t0'=t0+t7*/
+    paddw mm7,mm0
+    /*mm6=t1'=t1+t6*/
+    paddw mm6,mm1
+    /*mm5=t2'=t2+t5*/
+    paddw mm5,mm2
+    /*mm4=t3'=t3+t4*/
+    paddw mm4,mm3
+    OC_FDCT8x4(0x00,0x10,0x20,0x30,0x40,0x50,0x60,0x70)
+    OC_TRANSPOSE8x4(0x00,0x10,0x20,0x30,0x40,0x50,0x60,0x70)
+    /*Swap out this 8x4 block for the next one.*/
+    movq mm0,[0x08+X]
+    movq [0x30+Y],mm7
+    movq mm7,[0x78+X]
+    movq [0x50+Y],mm1
+    movq mm1,[0x18+X]
+    movq [0x20+Y],mm6
+    movq mm6,[0x68+X]
+    movq [0x60+Y],mm2
+    movq mm2,[0x28+X]
+    movq [0x10+Y],mm5
+    movq mm5,[0x58+X]
+    movq [0x70+Y],mm3
+    movq mm3,[0x38+X]
+    /*And increase its working precision, too.*/
+    psllw mm0,2
+    movq [0x00+Y],mm4
+    psllw mm7,2
+    movq mm4,[0x48+X]
+    /*We inline stage1 of the transform here so we can get better instruction
+       scheduling with the shifts.*/
+    /*mm0=t7'=t0-t7*/
+    psubw mm0,mm7
+    psllw mm1,2
+    paddw mm7,mm7
+    psllw mm6,2
+    /*mm1=t6'=t1-t6*/
+    psubw mm1,mm6
+    psllw mm2,2
+    paddw mm6,mm6
+    psllw mm5,2
+    /*mm2=t5'=t2-t5*/
+    psubw mm2,mm5
+    psllw mm3,2
+    paddw mm5,mm5
+    psllw mm4,2
+    /*mm3=t4'=t3-t4*/
+    psubw mm3,mm4
+    paddw mm4,mm4
+    /*mm7=t0'=t0+t7*/
+    paddw mm7,mm0
+    /*mm6=t1'=t1+t6*/
+    paddw mm6,mm1
+    /*mm5=t2'=t2+t5*/
+    paddw mm5,mm2
+    /*mm4=t3'=t3+t4*/
+    paddw mm4,mm3
+    OC_FDCT8x4(0x08,0x18,0x28,0x38,0x48,0x58,0x68,0x78)
+    OC_TRANSPOSE8x4(0x08,0x18,0x28,0x38,0x48,0x58,0x68,0x78)
+    /*Here the first 4x4 block of output from the last transpose is the second
+       4x4 block of input for the next transform.
+      We have cleverly arranged that it already be in the appropriate place,
+       so we only have to do half the stores and loads.*/
+    movq mm0,[0x00+Y]
+    movq [0x58+Y],mm1
+    movq mm1,[0x10+Y]
+    movq [0x68+Y],mm2
+    movq mm2,[0x20+Y]
+    movq [0x78+Y],mm3
+    movq mm3,[0x30+Y]
+    OC_FDCT_STAGE1_8x4
+    OC_FDCT8x4(0x00,0x10,0x20,0x30,0x08,0x18,0x28,0x38)
+    OC_TRANSPOSE8x4(0x00,0x10,0x20,0x30,0x08,0x18,0x28,0x38)
+    /*mm0={-2}x4*/
+    pcmpeqw mm0,mm0
+    paddw mm0,mm0
+    /*Round the results.*/
+    psubw mm1,mm0
+    psubw mm2,mm0
+    psraw mm1,2
+    psubw mm3,mm0
+    movq [0x18+Y],mm1
+    psraw mm2,2
+    psubw mm4,mm0
+    movq mm1,[0x08+Y]
+    psraw mm3,2
+    psubw mm5,mm0
+    psraw mm4,2
+    psubw mm6,mm0
+    psraw mm5,2
+    psubw mm7,mm0
+    psraw mm6,2
+    psubw mm1,mm0
+    psraw mm7,2
+    movq mm0,[0x40+Y]
+    psraw mm1,2
+    movq [0x30+Y],mm7
+    movq mm7,[0x78+Y]
+    movq [0x08+Y],mm1
+    movq mm1,[0x50+Y]
+    movq [0x20+Y],mm6
+    movq mm6,[0x68+Y]
+    movq [0x28+Y],mm2
+    movq mm2,[0x60+Y]
+    movq [0x10+Y],mm5
+    movq mm5,[0x58+Y]
+    movq [0x38+Y],mm3
+    movq mm3,[0x70+Y]
+    movq [0x00+Y],mm4
+    movq mm4,[0x48+Y]
+    OC_FDCT_STAGE1_8x4
+    OC_FDCT8x4(0x40,0x50,0x60,0x70,0x48,0x58,0x68,0x78)
+    OC_TRANSPOSE8x4(0x40,0x50,0x60,0x70,0x48,0x58,0x68,0x78)
+    /*mm0={-2}x4*/
+    pcmpeqw mm0,mm0
+    paddw mm0,mm0
+    /*Round the results.*/
+    psubw mm1,mm0
+    psubw mm2,mm0
+    psraw mm1,2
+    psubw mm3,mm0
+    movq [0x58+Y],mm1
+    psraw mm2,2
+    psubw mm4,mm0
+    movq mm1,[0x48+Y]
+    psraw mm3,2
+    psubw mm5,mm0
+    movq [0x68+Y],mm2
+    psraw mm4,2
+    psubw mm6,mm0
+    movq [0x78+Y],mm3
+    psraw mm5,2
+    psubw mm7,mm0
+    movq [0x40+Y],mm4
+    psraw mm6,2
+    psubw mm1,mm0
+    movq [0x50+Y],mm5
+    psraw mm7,2
+    movq [0x60+Y],mm6
+    psraw mm1,2
+    movq [0x70+Y],mm7
+    movq [0x48+Y],mm1
+#undef Y
+#undef A
+#undef X
+  }
+}
+
+#endif