--- /dev/null
+
+#include <stdio.h>
+#include <stdlib.h>
+
+
+/* DO NOT COMPILE WITH -O/-O2/-O3 ! GENERATES INVALID ASSEMBLY. */
+
+
+/* mmx.h
+
+ MultiMedia eXtensions GCC interface library for IA32.
+
+ To use this library, simply include this header file
+ and compile with GCC. You MUST have inlining enabled
+ in order for mmx_ok() to work; this can be done by
+ simply using -O on the GCC command line.
+
+ Compiling with -DMMX_TRACE will cause detailed trace
+ output to be sent to stderr for each mmx operation.
+ This adds lots of code, and obviously slows execution to
+ a crawl, but can be very useful for debugging.
+
+ THIS SOFTWARE IS PROVIDED ``AS IS'' AND WITHOUT ANY
+ EXPRESS OR IMPLIED WARRANTIES, INCLUDING, WITHOUT
+ LIMITATION, THE IMPLIED WARRANTIES OF MERCHANTABILITY
+ AND FITNESS FOR ANY PARTICULAR PURPOSE.
+
+ June 11, 1998 by H. Dietz and R. Fisher
+*/
+
+
+/* The type of an value that fits in an MMX register
+ (note that long long constant values MUST be suffixed
+ by LL and unsigned long long values by ULL, lest
+ they be truncated by the compiler)
+*/
+typedef union {
+ long long q; /* Quadword (64-bit) value */
+ unsigned long long uq; /* Unsigned Quadword */
+ int d[2]; /* 2 Doubleword (32-bit) values */
+ unsigned int ud[2]; /* 2 Unsigned Doubleword */
+ short w[4]; /* 4 Word (16-bit) values */
+ unsigned short uw[4]; /* 4 Unsigned Word */
+ char b[8]; /* 8 Byte (8-bit) values */
+ unsigned char ub[8]; /* 8 Unsigned Byte */
+} mmx_t;
+
+
+/* Function to test if mmx instructions are supported...
+*/
+inline extern int
+mmx_ok(void)
+{
+ /* Returns 1 if mmx instructions are ok,
+ 0 if hardware does not support mmx
+ */
+ register int ok = 0;
+
+ __asm__ __volatile__ (
+ /* Get CPU version information */
+ "movl $1, %%eax\n\t"
+ "cpuid\n\t"
+ "movl %%edx, %0"
+ : "=a" (ok)
+ : /* no input */
+ );
+ return((ok & 0x800000) == 0x800000);
+}
+
+
+/* Helper functions for the instruction macros that follow...
+ (note that memory-to-register, m2r, instructions are nearly
+ as efficient as register-to-register, r2r, instructions;
+ however, memory-to-memory instructions are really simulated
+ as a convenience, and are only 1/3 as efficient)
+*/
+#ifdef MMX_TRACE
+
+/* Include the stuff for printing a trace to stderr...
+*/
+
+#include <stdio.h>
+
+#define mmx_m2r(op, mem, reg) \
+ { \
+ mmx_t mmx_trace; \
+ mmx_trace = (mem); \
+ fprintf(stderr, #op "_m2r(" #mem "=0x%016llx, ", mmx_trace.q); \
+ __asm__ __volatile__ ("movq %%" #reg ", %0" \
+ : "=X" (mmx_trace) \
+ : /* nothing */ ); \
+ fprintf(stderr, #reg "=0x%016llx) => ", mmx_trace.q); \
+ __asm__ __volatile__ (#op " %0, %%" #reg \
+ : /* nothing */ \
+ : "X" (mem)); \
+ __asm__ __volatile__ ("movq %%" #reg ", %0" \
+ : "=X" (mmx_trace) \
+ : /* nothing */ ); \
+ fprintf(stderr, #reg "=0x%016llx\n", mmx_trace.q); \
+ }
+
+#define mmx_r2m(op, reg, mem) \
+ { \
+ mmx_t mmx_trace; \
+ __asm__ __volatile__ ("movq %%" #reg ", %0" \
+ : "=X" (mmx_trace) \
+ : /* nothing */ ); \
+ fprintf(stderr, #op "_r2m(" #reg "=0x%016llx, ", mmx_trace.q); \
+ mmx_trace = (mem); \
+ fprintf(stderr, #mem "=0x%016llx) => ", mmx_trace.q); \
+ __asm__ __volatile__ (#op " %%" #reg ", %0" \
+ : "=X" (mem) \
+ : /* nothing */ ); \
+ mmx_trace = (mem); \
+ fprintf(stderr, #mem "=0x%016llx\n", mmx_trace.q); \
+ }
+
+#define mmx_r2r(op, regs, regd) \
+ { \
+ mmx_t mmx_trace; \
+ __asm__ __volatile__ ("movq %%" #regs ", %0" \
+ : "=X" (mmx_trace) \
+ : /* nothing */ ); \
+ fprintf(stderr, #op "_r2r(" #regs "=0x%016llx, ", mmx_trace.q); \
+ __asm__ __volatile__ ("movq %%" #regd ", %0" \
+ : "=X" (mmx_trace) \
+ : /* nothing */ ); \
+ fprintf(stderr, #regd "=0x%016llx) => ", mmx_trace.q); \
+ __asm__ __volatile__ (#op " %" #regs ", %" #regd); \
+ __asm__ __volatile__ ("movq %%" #regd ", %0" \
+ : "=X" (mmx_trace) \
+ : /* nothing */ ); \
+ fprintf(stderr, #regd "=0x%016llx\n", mmx_trace.q); \
+ }
+
+#define mmx_m2m(op, mems, memd) \
+ { \
+ mmx_t mmx_trace; \
+ mmx_trace = (mems); \
+ fprintf(stderr, #op "_m2m(" #mems "=0x%016llx, ", mmx_trace.q); \
+ mmx_trace = (memd); \
+ fprintf(stderr, #memd "=0x%016llx) => ", mmx_trace.q); \
+ __asm__ __volatile__ ("movq %0, %%mm0\n\t" \
+ #op " %1, %%mm0\n\t" \
+ "movq %%mm0, %0" \
+ : "=X" (memd) \
+ : "X" (mems)); \
+ mmx_trace = (memd); \
+ fprintf(stderr, #memd "=0x%016llx\n", mmx_trace.q); \
+ }
+
+#else
+
+/* These macros are a lot simpler without the tracing...
+*/
+
+#define mmx_m2r(op, mem, reg) \
+ __asm__ __volatile__ (#op " %0, %%" #reg \
+ : /* nothing */ \
+ : "X" (mem))
+
+#define mmx_r2m(op, reg, mem) \
+ __asm__ __volatile__ (#op " %%" #reg ", %0" \
+ : "=X" (mem) \
+ : /* nothing */ )
+
+#define mmx_r2r(op, regs, regd) \
+ __asm__ __volatile__ (#op " %" #regs ", %" #regd)
+
+#define mmx_m2m(op, mems, memd) \
+ __asm__ __volatile__ ("movq %0, %%mm0\n\t" \
+ #op " %1, %%mm0\n\t" \
+ "movq %%mm0, %0" \
+ : "=X" (memd) \
+ : "X" (mems))
+
+#endif
+
+
+/* 1x64 MOVe Quadword
+ (this is both a load and a store...
+ in fact, it is the only way to store)
+*/
+#define movq_m2r(var, reg) mmx_m2r(movq, var, reg)
+#define movq_r2m(reg, var) mmx_r2m(movq, reg, var)
+#define movq_r2r(regs, regd) mmx_r2r(movq, regs, regd)
+#define movq(vars, vard) \
+ __asm__ __volatile__ ("movq %1, %%mm0\n\t" \
+ "movq %%mm0, %0" \
+ : "=X" (vard) \
+ : "X" (vars))
+
+
+/* 1x64 MOVe Doubleword
+ (like movq, this is both load and store...
+ but is most useful for moving things between
+ mmx registers and ordinary registers)
+*/
+#define movd_m2r(var, reg) mmx_m2r(movd, var, reg)
+#define movd_r2m(reg, var) mmx_r2m(movd, reg, var)
+#define movd_r2r(regs, regd) mmx_r2r(movd, regs, regd)
+#define movd(vars, vard) \
+ __asm__ __volatile__ ("movd %1, %%mm0\n\t" \
+ "movd %%mm0, %0" \
+ : "=X" (vard) \
+ : "X" (vars))
+
+
+/* 2x32, 4x16, and 8x8 Parallel ADDs
+*/
+#define paddd_m2r(var, reg) mmx_m2r(paddd, var, reg)
+#define paddd_r2r(regs, regd) mmx_r2r(paddd, regs, regd)
+#define paddd(vars, vard) mmx_m2m(paddd, vars, vard)
+
+#define paddw_m2r(var, reg) mmx_m2r(paddw, var, reg)
+#define paddw_r2r(regs, regd) mmx_r2r(paddw, regs, regd)
+#define paddw(vars, vard) mmx_m2m(paddw, vars, vard)
+
+#define paddb_m2r(var, reg) mmx_m2r(paddb, var, reg)
+#define paddb_r2r(regs, regd) mmx_r2r(paddb, regs, regd)
+#define paddb(vars, vard) mmx_m2m(paddb, vars, vard)
+
+
+/* 4x16 and 8x8 Parallel ADDs using Saturation arithmetic
+*/
+#define paddsw_m2r(var, reg) mmx_m2r(paddsw, var, reg)
+#define paddsw_r2r(regs, regd) mmx_r2r(paddsw, regs, regd)
+#define paddsw(vars, vard) mmx_m2m(paddsw, vars, vard)
+
+#define paddsb_m2r(var, reg) mmx_m2r(paddsb, var, reg)
+#define paddsb_r2r(regs, regd) mmx_r2r(paddsb, regs, regd)
+#define paddsb(vars, vard) mmx_m2m(paddsb, vars, vard)
+
+
+/* 4x16 and 8x8 Parallel ADDs using Unsigned Saturation arithmetic
+*/
+#define paddusw_m2r(var, reg) mmx_m2r(paddusw, var, reg)
+#define paddusw_r2r(regs, regd) mmx_r2r(paddusw, regs, regd)
+#define paddusw(vars, vard) mmx_m2m(paddusw, vars, vard)
+
+#define paddusb_m2r(var, reg) mmx_m2r(paddusb, var, reg)
+#define paddusb_r2r(regs, regd) mmx_r2r(paddusb, regs, regd)
+#define paddusb(vars, vard) mmx_m2m(paddusb, vars, vard)
+
+
+/* 2x32, 4x16, and 8x8 Parallel SUBs
+*/
+#define psubd_m2r(var, reg) mmx_m2r(psubd, var, reg)
+#define psubd_r2r(regs, regd) mmx_r2r(psubd, regs, regd)
+#define psubd(vars, vard) mmx_m2m(psubd, vars, vard)
+
+#define psubw_m2r(var, reg) mmx_m2r(psubw, var, reg)
+#define psubw_r2r(regs, regd) mmx_r2r(psubw, regs, regd)
+#define psubw(vars, vard) mmx_m2m(psubw, vars, vard)
+
+#define psubb_m2r(var, reg) mmx_m2r(psubb, var, reg)
+#define psubb_r2r(regs, regd) mmx_r2r(psubb, regs, regd)
+#define psubb(vars, vard) mmx_m2m(psubb, vars, vard)
+
+
+/* 4x16 and 8x8 Parallel SUBs using Saturation arithmetic
+*/
+#define psubsw_m2r(var, reg) mmx_m2r(psubsw, var, reg)
+#define psubsw_r2r(regs, regd) mmx_r2r(psubsw, regs, regd)
+#define psubsw(vars, vard) mmx_m2m(psubsw, vars, vard)
+
+#define psubsb_m2r(var, reg) mmx_m2r(psubsb, var, reg)
+#define psubsb_r2r(regs, regd) mmx_r2r(psubsb, regs, regd)
+#define psubsb(vars, vard) mmx_m2m(psubsb, vars, vard)
+
+
+/* 4x16 and 8x8 Parallel SUBs using Unsigned Saturation arithmetic
+*/
+#define psubusw_m2r(var, reg) mmx_m2r(psubusw, var, reg)
+#define psubusw_r2r(regs, regd) mmx_r2r(psubusw, regs, regd)
+#define psubusw(vars, vard) mmx_m2m(psubusw, vars, vard)
+
+#define psubusb_m2r(var, reg) mmx_m2r(psubusb, var, reg)
+#define psubusb_r2r(regs, regd) mmx_r2r(psubusb, regs, regd)
+#define psubusb(vars, vard) mmx_m2m(psubusb, vars, vard)
+
+
+/* 4x16 Parallel MULs giving Low 4x16 portions of results
+*/
+#define pmullw_m2r(var, reg) mmx_m2r(pmullw, var, reg)
+#define pmullw_r2r(regs, regd) mmx_r2r(pmullw, regs, regd)
+#define pmullw(vars, vard) mmx_m2m(pmullw, vars, vard)
+
+
+/* 4x16 Parallel MULs giving High 4x16 portions of results
+*/
+#define pmulhw_m2r(var, reg) mmx_m2r(pmulhw, var, reg)
+#define pmulhw_r2r(regs, regd) mmx_r2r(pmulhw, regs, regd)
+#define pmulhw(vars, vard) mmx_m2m(pmulhw, vars, vard)
+
+
+/* 4x16->2x32 Parallel Mul-ADD
+ (muls like pmullw, then adds adjacent 16-bit fields
+ in the multiply result to make the final 2x32 result)
+*/
+#define pmaddwd_m2r(var, reg) mmx_m2r(pmaddwd, var, reg)
+#define pmaddwd_r2r(regs, regd) mmx_r2r(pmaddwd, regs, regd)
+#define pmaddwd(vars, vard) mmx_m2m(pmaddwd, vars, vard)
+
+
+/* 1x64 bitwise AND
+*/
+#ifdef BROKEN_PAND
+#define pand_m2r(var, reg) \
+ { \
+ mmx_m2r(pandn, (mmx_t) -1LL, reg); \
+ mmx_m2r(pandn, var, reg); \
+ }
+#define pand_r2r(regs, regd) \
+ { \
+ mmx_m2r(pandn, (mmx_t) -1LL, regd); \
+ mmx_r2r(pandn, regs, regd) \
+ }
+#define pand(vars, vard) \
+ { \
+ movq_m2r(vard, mm0); \
+ mmx_m2r(pandn, (mmx_t) -1LL, mm0); \
+ mmx_m2r(pandn, vars, mm0); \
+ movq_r2m(mm0, vard); \
+ }
+#else
+#define pand_m2r(var, reg) mmx_m2r(pand, var, reg)
+#define pand_r2r(regs, regd) mmx_r2r(pand, regs, regd)
+#define pand(vars, vard) mmx_m2m(pand, vars, vard)
+#endif
+
+
+/* 1x64 bitwise AND with Not the destination
+*/
+#define pandn_m2r(var, reg) mmx_m2r(pandn, var, reg)
+#define pandn_r2r(regs, regd) mmx_r2r(pandn, regs, regd)
+#define pandn(vars, vard) mmx_m2m(pandn, vars, vard)
+
+
+/* 1x64 bitwise OR
+*/
+#define por_m2r(var, reg) mmx_m2r(por, var, reg)
+#define por_r2r(regs, regd) mmx_r2r(por, regs, regd)
+#define por(vars, vard) mmx_m2m(por, vars, vard)
+
+
+/* 1x64 bitwise eXclusive OR
+*/
+#define pxor_m2r(var, reg) mmx_m2r(pxor, var, reg)
+#define pxor_r2r(regs, regd) mmx_r2r(pxor, regs, regd)
+#define pxor(vars, vard) mmx_m2m(pxor, vars, vard)
+
+
+/* 2x32, 4x16, and 8x8 Parallel CoMPare for EQuality
+ (resulting fields are either 0 or -1)
+*/
+#define pcmpeqd_m2r(var, reg) mmx_m2r(pcmpeqd, var, reg)
+#define pcmpeqd_r2r(regs, regd) mmx_r2r(pcmpeqd, regs, regd)
+#define pcmpeqd(vars, vard) mmx_m2m(pcmpeqd, vars, vard)
+
+#define pcmpeqw_m2r(var, reg) mmx_m2r(pcmpeqw, var, reg)
+#define pcmpeqw_r2r(regs, regd) mmx_r2r(pcmpeqw, regs, regd)
+#define pcmpeqw(vars, vard) mmx_m2m(pcmpeqw, vars, vard)
+
+#define pcmpeqb_m2r(var, reg) mmx_m2r(pcmpeqb, var, reg)
+#define pcmpeqb_r2r(regs, regd) mmx_r2r(pcmpeqb, regs, regd)
+#define pcmpeqb(vars, vard) mmx_m2m(pcmpeqb, vars, vard)
+
+
+/* 2x32, 4x16, and 8x8 Parallel CoMPare for Greater Than
+ (resulting fields are either 0 or -1)
+*/
+#define pcmpgtd_m2r(var, reg) mmx_m2r(pcmpgtd, var, reg)
+#define pcmpgtd_r2r(regs, regd) mmx_r2r(pcmpgtd, regs, regd)
+#define pcmpgtd(vars, vard) mmx_m2m(pcmpgtd, vars, vard)
+
+#define pcmpgtw_m2r(var, reg) mmx_m2r(pcmpgtw, var, reg)
+#define pcmpgtw_r2r(regs, regd) mmx_r2r(pcmpgtw, regs, regd)
+#define pcmpgtw(vars, vard) mmx_m2m(pcmpgtw, vars, vard)
+
+#define pcmpgtb_m2r(var, reg) mmx_m2r(pcmpgtb, var, reg)
+#define pcmpgtb_r2r(regs, regd) mmx_r2r(pcmpgtb, regs, regd)
+#define pcmpgtb(vars, vard) mmx_m2m(pcmpgtb, vars, vard)
+
+
+/* 1x64, 2x32, and 4x16 Parallel Shift Left Logical
+*/
+#define psllq_m2r(var, reg) mmx_m2r(psllq, var, reg)
+#define psllq_r2r(regs, regd) mmx_r2r(psllq, regs, regd)
+#define psllq(vars, vard) mmx_m2m(psllq, vars, vard)
+
+#define pslld_m2r(var, reg) mmx_m2r(pslld, var, reg)
+#define pslld_r2r(regs, regd) mmx_r2r(pslld, regs, regd)
+#define pslld(vars, vard) mmx_m2m(pslld, vars, vard)
+
+#define psllw_m2r(var, reg) mmx_m2r(psllw, var, reg)
+#define psllw_r2r(regs, regd) mmx_r2r(psllw, regs, regd)
+#define psllw(vars, vard) mmx_m2m(psllw, vars, vard)
+
+
+/* 1x64, 2x32, and 4x16 Parallel Shift Right Logical
+*/
+#define psrlq_m2r(var, reg) mmx_m2r(psrlq, var, reg)
+#define psrlq_r2r(regs, regd) mmx_r2r(psrlq, regs, regd)
+#define psrlq(vars, vard) mmx_m2m(psrlq, vars, vard)
+
+#define psrld_m2r(var, reg) mmx_m2r(psrld, var, reg)
+#define psrld_r2r(regs, regd) mmx_r2r(psrld, regs, regd)
+#define psrld(vars, vard) mmx_m2m(psrld, vars, vard)
+
+#define psrlw_m2r(var, reg) mmx_m2r(psrlw, var, reg)
+#define psrlw_r2r(regs, regd) mmx_r2r(psrlw, regs, regd)
+#define psrlw(vars, vard) mmx_m2m(psrlw, vars, vard)
+
+
+/* 2x32 and 4x16 Parallel Shift Right Arithmetic
+*/
+#define psrad_m2r(var, reg) mmx_m2r(psrad, var, reg)
+#define psrad_r2r(regs, regd) mmx_r2r(psrad, regs, regd)
+#define psrad(vars, vard) mmx_m2m(psrad, vars, vard)
+
+#define psraw_m2r(var, reg) mmx_m2r(psraw, var, reg)
+#define psraw_r2r(regs, regd) mmx_r2r(psraw, regs, regd)
+#define psraw(vars, vard) mmx_m2m(psraw, vars, vard)
+
+
+/* 2x32->4x16 and 4x16->8x8 PACK and Signed Saturate
+ (packs source and dest fields into dest in that order)
+*/
+#define packssdw_m2r(var, reg) mmx_m2r(packssdw, var, reg)
+#define packssdw_r2r(regs, regd) mmx_r2r(packssdw, regs, regd)
+#define packssdw(vars, vard) mmx_m2m(packssdw, vars, vard)
+
+#define packsswb_m2r(var, reg) mmx_m2r(packsswb, var, reg)
+#define packsswb_r2r(regs, regd) mmx_r2r(packsswb, regs, regd)
+#define packsswb(vars, vard) mmx_m2m(packsswb, vars, vard)
+
+
+/* 4x16->8x8 PACK and Unsigned Saturate
+ (packs source and dest fields into dest in that order)
+*/
+#define packuswb_m2r(var, reg) mmx_m2r(packuswb, var, reg)
+#define packuswb_r2r(regs, regd) mmx_r2r(packuswb, regs, regd)
+#define packuswb(vars, vard) mmx_m2m(packuswb, vars, vard)
+
+
+/* 2x32->1x64, 4x16->2x32, and 8x8->4x16 UNPaCK Low
+ (interleaves low half of dest with low half of source
+ as padding in each result field)
+*/
+#define punpckldq_m2r(var, reg) mmx_m2r(punpckldq, var, reg)
+#define punpckldq_r2r(regs, regd) mmx_r2r(punpckldq, regs, regd)
+#define punpckldq(vars, vard) mmx_m2m(punpckldq, vars, vard)
+
+#define punpcklwd_m2r(var, reg) mmx_m2r(punpcklwd, var, reg)
+#define punpcklwd_r2r(regs, regd) mmx_r2r(punpcklwd, regs, regd)
+#define punpcklwd(vars, vard) mmx_m2m(punpcklwd, vars, vard)
+
+#define punpcklbw_m2r(var, reg) mmx_m2r(punpcklbw, var, reg)
+#define punpcklbw_r2r(regs, regd) mmx_r2r(punpcklbw, regs, regd)
+#define punpcklbw(vars, vard) mmx_m2m(punpcklbw, vars, vard)
+
+
+/* 2x32->1x64, 4x16->2x32, and 8x8->4x16 UNPaCK High
+ (interleaves high half of dest with high half of source
+ as padding in each result field)
+*/
+#define punpckhdq_m2r(var, reg) mmx_m2r(punpckhdq, var, reg)
+#define punpckhdq_r2r(regs, regd) mmx_r2r(punpckhdq, regs, regd)
+#define punpckhdq(vars, vard) mmx_m2m(punpckhdq, vars, vard)
+
+#define punpckhwd_m2r(var, reg) mmx_m2r(punpckhwd, var, reg)
+#define punpckhwd_r2r(regs, regd) mmx_r2r(punpckhwd, regs, regd)
+#define punpckhwd(vars, vard) mmx_m2m(punpckhwd, vars, vard)
+
+#define punpckhbw_m2r(var, reg) mmx_m2r(punpckhbw, var, reg)
+#define punpckhbw_r2r(regs, regd) mmx_r2r(punpckhbw, regs, regd)
+#define punpckhbw(vars, vard) mmx_m2m(punpckhbw, vars, vard)
+
+
+/* Empty MMx State
+ (used to clean-up when going from mmx to float use
+ of the registers that are shared by both; note that
+ there is no float-to-mmx operation needed, because
+ only the float tag word info is corruptible)
+*/
+#ifdef MMX_TRACE
+
+#define emms() \
+ { \
+ fprintf(stderr, "emms()\n"); \
+ __asm__ __volatile__ ("emms"); \
+ }
+
+#else
+
+#define emms() __asm__ __volatile__ ("emms")
+
+#endif
+
+
+
+int main()
+{
+ // int rval;
+ mmx_t ma;
+ mmx_t mb;
+
+ movq_r2r(mm0, mm1);
+
+// rval = mmx_ok();
+
+ /* Announce return value of mmx_ok() */
+// printf("Value returned from init was %x.", rval);
+// printf(" (Indicates MMX %s available)\n\n",(rval)? "is" : "not");
+// fflush(stdout); fflush(stdout);
+
+// if(rval)
+ {
+ /* PADD *****************************************************/
+ ma.q = 0x1111111180000000LL;
+ mb.q = 0x7fffffff00000001LL;
+ paddd(ma, mb);
+ fprintf(stdout, "paddd: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "paddd: mb.q is 9111111080000001\n");
+ fflush(stdout);
+
+ ma.q = 0x0001000100010001LL;
+ mb.q = 0x80007fffffff0001LL;
+ paddw(ma, mb);
+ fprintf(stdout, "paddw: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "paddw: mb.q is 8001800000000002\n");
+ fflush(stdout);
+
+ ma.q = 0x80007fffffff0001LL;
+ mb.q = 0x0001000100010000LL;
+ paddw(ma, mb);
+ fprintf(stdout, "paddw: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "paddw: mb.q is 8001800000000001\n");
+ fflush(stdout);
+
+ ma.q = 0x01010101807fff01LL;
+ mb.q = 0x807fff0101010101LL;
+ paddb(ma, mb);
+ fprintf(stdout, "paddb: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "paddb: mb.q is 8180000281800002\n");
+ fflush(stdout);
+
+
+ /* PADDS ****************************************************/
+ ma.q = 0x0001000100010001LL;
+ mb.q = 0x80007fffffff0001LL;
+ paddsw(ma, mb);
+ fprintf(stdout, "paddsw: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "paddsw: mb.q is 80017fff00000002\n");
+
+ ma.q = 0x80007fffffff0001LL;
+ mb.q = 0x0001000100010000LL;
+ paddsw(ma, mb);
+ fprintf(stdout, "paddsw: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "paddsw: mb.q is 80017fff00000001\n");
+
+ ma.q = 0x01010101807fff01LL;
+ mb.q = 0x807fff0101010101LL;
+ paddsb(ma, mb);
+ fprintf(stdout, "paddsb: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "paddsb: mb.q is 817f0002817f0002\n");
+ fflush(stdout);
+
+
+ /* PADDUS ***************************************************/
+ ma.q = 0x0001000100010001LL;
+ mb.q = 0x80007fffffff0001LL;
+ paddusw(ma, mb);
+ fprintf(stdout, "paddusw: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "paddusw: mb.q is 80018000ffff0002\n");
+ fflush(stdout);
+
+ ma.q = 0x80007fffffff0001LL;
+ mb.q = 0x0001000100010000LL;
+ paddusw(ma, mb);
+ fprintf(stdout, "paddusw: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "paddusw: mb.q is 80018000ffff0001\n");
+ fflush(stdout);
+
+ ma.q = 0x01010101807fff01LL;
+ mb.q = 0x807fff0101010101LL;
+ paddusb(ma, mb);
+ fprintf(stdout, "paddusb: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "paddusb: mb.q is 8180ff028180ff02\n");
+ fflush(stdout);
+
+
+ /* PSUB *****************************************************/
+ ma.q = 0x7fffffff00000001LL;
+ mb.q = 0x1111111180000000LL;
+ psubd(ma, mb);
+ fprintf(stdout, "psubd: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "psubd: mb.q is 911111127fffffff\n");
+ fflush(stdout);
+
+ ma.q = 0x80007fffffff0001LL;
+ mb.q = 0x0001000100010001LL;
+ psubw(ma, mb);
+ fprintf(stdout, "psubw: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "psubw: mb.q is 8001800200020000\n");
+ fflush(stdout);
+
+ ma.q = 0x0001000100010000LL;
+ mb.q = 0x80007fffffff0001LL;
+ psubw(ma, mb);
+ fprintf(stdout, "psubw: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "psubw: mb.q is 7fff7ffefffe0001\n");
+ fflush(stdout);
+
+ ma.q = 0x807fff0101010101LL;
+ mb.q = 0x01010101807fff01LL;
+ psubb(ma, mb);
+ fprintf(stdout, "psubb: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "psubb: mb.q is 818202007f7efe00\n");
+ fflush(stdout);
+
+
+ /* PSUBS ****************************************************/
+ ma.q = 0x80007fffffff0001LL;
+ mb.q = 0x0001000100010001LL;
+ psubsw(ma, mb);
+ fprintf(stdout, "psubsw: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "psubsw: mb.q is 7fff800200020000\n");
+ fflush(stdout);
+
+ ma.q = 0x0001000100010000LL;
+ mb.q = 0x80007fffffff0001LL;
+ psubsw(ma, mb);
+ fprintf(stdout, "psubsw: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "psubsw: mb.q is 80007ffefffe0001\n");
+ fflush(stdout);
+
+ ma.q = 0x807fff0101010101LL;
+ mb.q = 0x01010101807fff01LL;
+ psubsb(ma, mb);
+ fprintf(stdout, "psubsb: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "psubsb: mb.q is 7f820200807efe00\n");
+ fflush(stdout);
+
+
+ /* PSUBUS ***************************************************/
+ ma.q = 0x80007fffffff0001LL;
+ mb.q = 0x0001000100010001LL;
+ psubusw(ma, mb);
+ fprintf(stdout, "psubusw: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "psubusw: mb.q is 0000000000000000\n");
+ fflush(stdout);
+
+ ma.q = 0x0001000100010000LL;
+ mb.q = 0x80007fffffff0001LL;
+ psubusw(ma, mb);
+ fprintf(stdout, "psubusw: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "psubusw: mb.q is 7fff7ffefffe0001\n");
+ fflush(stdout);
+
+ ma.q = 0x807fff0101010101LL;
+ mb.q = 0x01010101807fff01LL;
+ psubusb(ma, mb);
+ fprintf(stdout, "psubusb: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "psubusb: mb.q is 000000007f7efe00\n");
+ fflush(stdout);
+
+
+ /* PMUL *****************************************************/
+ ma.q = 0x8000ffff00ff0000LL;
+ mb.q = 0x0200ffff00ffffffLL;
+ pmulhw(ma, mb);
+ fprintf(stdout, "pmulhw: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "pmulhw: mb.q is ff00000000000000\n");
+ fflush(stdout);
+
+ mb.q = 0x0200ffff00ffffffLL;
+ pmullw(ma, mb);
+ fprintf(stdout, "pmullw: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "pmullw: mb.q is 00000001fe010000\n");
+ fflush(stdout);
+
+
+ /* PMADD ****************************************************/
+ ma.q = 0x8000345680007f34LL;
+ mb.q = 0x93234a27ffff1707LL;
+
+ pmaddwd(ma, mb);
+ fprintf(stdout, "pmaddwd: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "pmaddwd: mb.q is 4597551a0b71a66c\n");
+ fflush(stdout);
+
+ /* PMADD ****************************************************/
+ ma.q = 0x8000800080008000LL;
+ mb.q = 0x8000800080008000LL;
+
+ pmaddwd(ma, mb);
+ fprintf(stdout, "pmaddwd: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "pmaddwd: mb.q is 8000000080000000\n");
+ fflush(stdout);
+
+
+ /* PCMPEQ ***************************************************/
+ ma.q = 0x800034568f237f34LL;
+ mb.q = 0x93009a568f237f34LL;
+
+ pcmpeqd(ma, mb);
+ fprintf(stdout, "pcmpeqd: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "pcmpeqd: mb.q is 00000000ffffffff\n");
+ fflush(stdout);
+
+ mb.q = 0x93009a568f237f34LL;
+ pcmpeqw(ma, mb);
+ fprintf(stdout, "pcmpeqw: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "pcmpeqw: mb.q is 00000000ffffffff\n");
+ fflush(stdout);
+
+ mb.q = 0x93009a568f237f34LL;
+ pcmpeqb(ma, mb);
+ fprintf(stdout, "pcmpeqb: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "pcmpeqb: mb.q is 00ff00ffffffffff\n");
+ fflush(stdout);
+
+
+
+ /* PCMPGT ***************************************************/
+ ma.q = 0x666688884477aaffLL;
+ mb.q = 0x1234567890abcdefLL;
+
+ pcmpgtd(ma, mb);
+ fprintf(stdout, "pcmpgtd: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "pcmpgtd: mb.q is 0000000000000000\n");
+ fflush(stdout);
+
+ mb.q = 0x1234567890abcdefLL;
+ pcmpgtw(ma, mb);
+ fprintf(stdout, "pcmpgtw: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "pcmpgtw: mb.q is 0000ffff0000ffff\n");
+ fflush(stdout);
+
+ mb.q = 0x1234567890abcdefLL;
+ pcmpgtb(ma, mb);
+ fprintf(stdout, "pcmpgtb: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "pcmpgtb: mb.q is 0000ffff0000ff00\n");
+ fflush(stdout);
+
+
+ /* PACKSS ***************************************************/
+ ma.q = 0x00012222000abbbbLL;
+ mb.q = 0x0000888800003333LL;
+
+ packssdw(ma, mb);
+ fprintf(stdout, "packssdw: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "packssdw: mb.q is 7fff7fff7fff3333\n");
+ fflush(stdout);
+
+ ma.q = 0x00aa00dd01009999LL;
+ mb.q = 0x0011002200330044LL;
+
+ packsswb(ma, mb);
+ fprintf(stdout, "packsswb: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "packsswb: mb.q is 7f7f7f8011223344\n");
+ fflush(stdout);
+
+
+ /* PACKUS ***************************************************/
+ ma.q = 0x00aa00dd01009999LL;
+ mb.q = 0x0011002200330044LL;
+
+ packuswb(ma, mb);
+ fprintf(stdout, "packuswb: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "packuswb: mb.q is aaddff0011223344\n");
+ fflush(stdout);
+
+
+ /* PUNPCKH **************************************************/
+ ma.q = 0x090a0b0c0d0e0f00LL;
+ mb.q = 0x0102030405060708LL;
+
+ punpckhdq(ma, mb);
+ fprintf(stdout, "punpckhdq: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "punpckhdq: mb.q is 090a0b0c01020304\n");
+ fflush(stdout);
+
+ mb.q = 0x0102030405060708LL;
+ punpckhwd(ma, mb);
+ fprintf(stdout, "punpckhwd: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "punpckhwd: mb.q is 090a01020b0c0304\n");
+ fflush(stdout);
+
+ mb.q = 0x0102030405060708LL;
+ punpckhbw(ma, mb);
+ fprintf(stdout, "punpckhbw: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "punpckhbw: mb.q is 09010a020b030c04\n");
+ fflush(stdout);
+
+
+ /* PUNPCKL **************************************************/
+ ma.q = 0x090a0b0c0d0e0f00LL;
+ mb.q = 0x0102030405060708LL;
+
+ punpckldq(ma, mb);
+ fprintf(stdout, "punpckldq: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "punpckldq: mb.q is 0d0e0f0005060708\n");
+ fflush(stdout);
+
+ mb.q = 0x0102030405060708LL;
+ punpcklwd(ma, mb);
+ fprintf(stdout, "punpcklwd: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "punpcklwd: mb.q is 0d0e05060f000708\n");
+ fflush(stdout);
+
+ mb.q = 0x0102030405060708LL;
+ punpcklbw(ma, mb);
+ fprintf(stdout, "punpcklbw: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "punpcklbw: mb.q is 0d050e060f070008\n");
+ fflush(stdout);
+
+
+
+ /* PAND, PANDN, POR, PXOR ***********************************/
+ ma.q = 0x5555555555555555LL;
+ mb.q = 0x3333333333333333LL;
+
+ pand(ma, mb);
+ fprintf(stdout, "pand: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "pand: mb.q is 1111111111111111\n");
+ fflush(stdout);
+
+ mb.q = 0x3333333333333333LL;
+ pandn(ma, mb);
+ fprintf(stdout, "pandn: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "pandn: mb.q is 4444444444444444\n");
+ fflush(stdout);
+
+ mb.q = 0x3333333333333333LL;
+ por(ma, mb);
+ fprintf(stdout, "por: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "por: mb.q is 7777777777777777\n");
+ fflush(stdout);
+
+ mb.q = 0x3333333333333333LL;
+ pxor(ma, mb);
+ fprintf(stdout, "pxor: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "pxor: mb.q is 6666666666666666\n");
+ fflush(stdout);
+
+
+
+ /* PSLL *****************************************************/
+ ma.q = 0x0000000000000018LL;
+ mb.q = 0x0123456789abcdefLL;
+
+ psllq(ma, mb);
+ fprintf(stdout, "psllq: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "psllq: mb.q is 6789abcdef000000\n");
+ fflush(stdout);
+
+ mb.q = 0x0123456789abcdefLL;
+ pslld(ma, mb);
+ fprintf(stdout, "pslld: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "pslld: mb.q is 67000000ef000000\n");
+ fflush(stdout);
+
+ mb.q = 0x0123456789abcdefLL;
+ psllw(ma, mb);
+ fprintf(stdout, "psllw: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "psllw: mb.q is 0000000000000000\n");
+ fflush(stdout);
+
+
+
+ /* PSRL *****************************************************/
+ ma.q = 0x0000000000000018LL;
+ mb.q = 0x0123456789abcdefLL;
+
+ psrlq(ma, mb);
+ fprintf(stdout, "psrlq: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "psrlq: mb.q is 0000000123456789\n");
+ fflush(stdout);
+
+ mb.q = 0x0123456789abcdefLL;
+ psrld(ma, mb);
+ fprintf(stdout, "psrld: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "psrld: mb.q is 0000000100000089\n");
+ fflush(stdout);
+
+ mb.q = 0x0123456789abcdefLL;
+ psrlw(ma, mb);
+ fprintf(stdout, "psrlw: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "psrlw: mb.q is 0000000000000000\n");
+ fflush(stdout);
+
+
+
+ /* PSRA *****************************************************/
+ ma.q = 0x0000000000000018LL;
+ mb.q = 0x0123456789abcdefLL;
+
+ psrad(ma, mb);
+ fprintf(stdout, "psrad: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "psrad: mb.q is 00000001ffffff89\n");
+ fflush(stdout);
+
+ mb.q = 0x0123456789abcdefLL;
+ psraw(ma, mb);
+ fprintf(stdout, "psraw: mb.q is %016llx\n", mb.q);
+ fprintf(stdout, "psraw: mb.q is 00000000ffffffff\n");
+ fflush(stdout);
+
+ /* Exit MXX *************************************************/
+ emms();
+ }
+
+ /* Clean-up and exit nicely */
+ exit(0);
+}