/*	ia32.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
*/



/*	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) \
	{ \
		p64_t p64_trace; \
		p64_trace = (mem); \
		fprintf(stderr, #op "_m2r(" #mem "=0x%016llx, ", p64_trace.q); \
		__asm__ __volatile__ ("movq %%" #reg ", %0" \
				      : "=X" (p64_trace) \
				      : /* nothing */ ); \
		fprintf(stderr, #reg "=0x%016llx) => ", p64_trace.q); \
		__asm__ __volatile__ (#op " %0, %%" #reg \
				      : /* nothing */ \
				      : "X" (mem)); \
		__asm__ __volatile__ ("movq %%" #reg ", %0" \
				      : "=X" (p64_trace) \
				      : /* nothing */ ); \
		fprintf(stderr, #reg "=0x%016llx\n", p64_trace.q); \
	}

#define	mmx_r2m(op, reg, mem) \
	{ \
		p64_t p64_trace; \
		__asm__ __volatile__ ("movq %%" #reg ", %0" \
				      : "=X" (p64_trace) \
				      : /* nothing */ ); \
		fprintf(stderr, #op "_r2m(" #reg "=0x%016llx, ", p64_trace.q); \
		p64_trace = (mem); \
		fprintf(stderr, #mem "=0x%016llx) => ", p64_trace.q); \
		__asm__ __volatile__ (#op " %%" #reg ", %0" \
				      : "=X" (mem) \
				      : /* nothing */ ); \
		p64_trace = (mem); \
		fprintf(stderr, #mem "=0x%016llx\n", p64_trace.q); \
	}

#define	mmx_i2r(op, immed, regd) \
	{ \
		p64_t p64_trace; \
		fprintf(stderr, #op "_i2r(" #immed ", "); \
		__asm__ __volatile__ ("movq %%" #regd ", %0" \
				      : "=X" (p64_trace) \
				      : /* nothing */ ); \
		fprintf(stderr, #regd "=0x%016llx) => ", p64_trace.q); \
		__asm__ __volatile__ (#op " $" #immed ", %" #regd) \
		__asm__ __volatile__ ("movq %%" #regd ", %0" \
				      : "=X" (p64_trace) \
				      : /* nothing */ ); \
		fprintf(stderr, #regd "=0x%016llx\n", p64_trace.q); \
	}


#define	mmx_r2r(op, regs, regd) \
	{ \
		p64_t p64_trace; \
		__asm__ __volatile__ ("movq %%" #regs ", %0" \
				      : "=X" (p64_trace) \
				      : /* nothing */ ); \
		fprintf(stderr, #op "_r2r(" #regs "=0x%016llx, ", p64_trace.q); \
		__asm__ __volatile__ ("movq %%" #regd ", %0" \
				      : "=X" (p64_trace) \
				      : /* nothing */ ); \
		fprintf(stderr, #regd "=0x%016llx) => ", p64_trace.q); \
		__asm__ __volatile__ (#op " %" #regs ", %" #regd); \
		__asm__ __volatile__ ("movq %%" #regd ", %0" \
				      : "=X" (p64_trace) \
				      : /* nothing */ ); \
		fprintf(stderr, #regd "=0x%016llx\n", p64_trace.q); \
	}

#define	mmx_m2m(op, mems, memd) \
	{ \
		p64_t p64_trace; \
		p64_trace = (mems); \
		fprintf(stderr, #op "_m2m(" #mems "=0x%016llx, ", p64_trace.q); \
		p64_trace = (memd); \
		fprintf(stderr, #memd "=0x%016llx) => ", p64_trace.q); \
		__asm__ __volatile__ ("movq %0, %%mm0\n\t" \
				      #op " %1, %%mm0\n\t" \
				      "movq %%mm0, %0" \
				      : "=X" (memd) \
				      : "X" (mems)); \
		p64_trace = (memd); \
		fprintf(stderr, #memd "=0x%016llx\n", p64_trace.q); \
	}

#else

/*	These macros are a lot simpler without the tracing...
*/

#define	ia32_m2r(op, mem, reg) \
	__asm__ __volatile__ (#op " %0, %%" #reg \
			      : /* nothing */ \
			      : "X" (mem))

#define	ia32_r2m(op, reg, mem) \
	__asm__ __volatile__ (#op " %%" #reg ", %0" \
			      : "=X" (mem) \
			      : /* nothing */ )

#define	ia32_i2r(op, immed, regd) \
	__asm__ __volatile__ (#op " $" #immed ", %" #regd)

#define	ia32_r2r(op, regs, regd) \
	__asm__ __volatile__ (#op " %" #regs ", %" #regd)

#define	ia32_m2m(op, mems, memd) \
	__asm__ __volatile__ ("movl %0, %%mm0\n\t" \
			      #op " %1, %%mm0\n\t" \
			      "movl %%mm0, %0" \
			      : "=X" (memd) \
			      : "X" (mems))

#endif


/*	1x32 MOVe Long
*/
#define	movl_m2r(var, reg)	ia32_m2r(movl, var, reg)
#define	movl_r2m(reg, var)	ia32_r2m(movl, reg, var)
#define	movl_r2r(regs, regd)	ia32_r2r(movl, regs, regd)


/*	1x32, 2x16, and 4x8 Parallel ADDs
*/
#define	addl_m2r(var, reg)	ia32_m2r(addl, var, reg)
#define	addl_r2r(regs, regd)	ia32_r2r(addl, regs, regd)
#define	addl(vars, vard)	ia32_m2m(addl, vars, vard)

#define	paddw_m2r(var, reg)	ia32_m2r(paddw, var, reg)
#define	paddw_r2r(regs, regd)	ia32_r2r(paddw, regs, regd)
#define	paddw(vars, vard)	ia32_m2m(paddw, vars, vard)

#define	paddb_m2r(var, reg)	ia32_m2r(paddb, var, reg)
#define	paddb_r2r(regs, regd)	ia32_r2r(paddb, regs, regd)
#define	paddb(vars, vard)	ia32_m2m(paddb, vars, vard)


/*	4x16 and 8x8 Parallel ADDs using Saturation arithmetic
*/
#define	paddsw_m2r(var, reg)	ia32_m2r(paddsw, var, reg)
#define	paddsw_r2r(regs, regd)	ia32_r2r(paddsw, regs, regd)
#define	paddsw(vars, vard)	ia32_m2m(paddsw, vars, vard)

#define	paddsb_m2r(var, reg)	ia32_m2r(paddsb, var, reg)
#define	paddsb_r2r(regs, regd)	ia32_r2r(paddsb, regs, regd)
#define	paddsb(vars, vard)	ia32_m2m(paddsb, vars, vard)


/*	4x16 and 8x8 Parallel ADDs using Unsigned Saturation arithmetic
*/
#define	paddusw_m2r(var, reg)	ia32_m2r(paddusw, var, reg)
#define	paddusw_r2r(regs, regd)	ia32_r2r(paddusw, regs, regd)
#define	paddusw(vars, vard)	ia32_m2m(paddusw, vars, vard)

#define	paddusb_m2r(var, reg)	ia32_m2r(paddusb, var, reg)
#define	paddusb_r2r(regs, regd)	ia32_r2r(paddusb, regs, regd)
#define	paddusb(vars, vard)	ia32_m2m(paddusb, vars, vard)


/*	2x32, 4x16, and 8x8 Parallel SUBs
*/
#define	psubd_m2r(var, reg)	ia32_m2r(psubd, var, reg)
#define	psubd_r2r(regs, regd)	ia32_r2r(psubd, regs, regd)
#define	psubd(vars, vard)	ia32_m2m(psubd, vars, vard)

#define	psubw_m2r(var, reg)	ia32_m2r(psubw, var, reg)
#define	psubw_r2r(regs, regd)	ia32_r2r(psubw, regs, regd)
#define	psubw(vars, vard)	ia32_m2m(psubw, vars, vard)

#define	psubb_m2r(var, reg)	ia32_m2r(psubb, var, reg)
#define	psubb_r2r(regs, regd)	ia32_r2r(psubb, regs, regd)
#define	psubb(vars, vard)	ia32_m2m(psubb, vars, vard)


/*	4x16 and 8x8 Parallel SUBs using Saturation arithmetic
*/
#define	psubsw_m2r(var, reg)	ia32_m2r(psubsw, var, reg)
#define	psubsw_r2r(regs, regd)	ia32_r2r(psubsw, regs, regd)
#define	psubsw(vars, vard)	ia32_m2m(psubsw, vars, vard)

#define	psubsb_m2r(var, reg)	ia32_m2r(psubsb, var, reg)
#define	psubsb_r2r(regs, regd)	ia32_r2r(psubsb, regs, regd)
#define	psubsb(vars, vard)	ia32_m2m(psubsb, vars, vard)


/*	4x16 and 8x8 Parallel SUBs using Unsigned Saturation arithmetic
*/
#define	psubusw_m2r(var, reg)	ia32_m2r(psubusw, var, reg)
#define	psubusw_r2r(regs, regd)	ia32_r2r(psubusw, regs, regd)
#define	psubusw(vars, vard)	ia32_m2m(psubusw, vars, vard)

#define	psubusb_m2r(var, reg)	ia32_m2r(psubusb, var, reg)
#define	psubusb_r2r(regs, regd)	ia32_r2r(psubusb, regs, regd)
#define	psubusb(vars, vard)	ia32_m2m(psubusb, vars, vard)


/*	4x16 Parallel MULs giving Low 4x16 portions of results
*/
#define	pmullw_m2r(var, reg)	ia32_m2r(pmullw, var, reg)
#define	pmullw_r2r(regs, regd)	ia32_r2r(pmullw, regs, regd)
#define	pmullw(vars, vard)	ia32_m2m(pmullw, vars, vard)


/*	4x16 Parallel MULs giving High 4x16 portions of results
*/
#define	pmulhw_m2r(var, reg)	ia32_m2r(pmulhw, var, reg)
#define	pmulhw_r2r(regs, regd)	ia32_r2r(pmulhw, regs, regd)
#define	pmulhw(vars, vard)	ia32_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)	ia32_m2r(pmaddwd, var, reg)
#define	pmaddwd_r2r(regs, regd)	ia32_r2r(pmaddwd, regs, regd)
#define	pmaddwd(vars, vard)	ia32_m2m(pmaddwd, vars, vard)


/*	1x32 bitwise AND
*/
#define	andl_m2r(var, reg)	ia32_m2r(andl, var, reg)
#define	andl_r2r(regs, regd)	ia32_r2r(andl, regs, regd)
#define	andl(vars, vard)	ia32_m2m(andl, vars, vard)


/*	1x32 bitwise AND with Not the destination
	This doesn't exist in IA32?
*/
#define	andnl_m2r(var, reg)	ia32_m2r(andnl, var, reg)
#define	andnl_r2r(regs, regd)	ia32_r2r(andnl, regs, regd)
#define	andnl(vars, vard)	ia32_m2m(andnl, vars, vard)


/*	1x32 bitwise OR
*/
#define	orl_m2r(var, reg)	ia32_m2r(orl, var, reg)
#define	orl_r2r(regs, regd)	ia32_r2r(orl, regs, regd)
#define	orl(vars, vard)		ia32_m2m(orl, vars, vard)


/*	1x64 bitwise eXclusive OR
*/
#define	xorl_m2r(var, reg)	ia32_m2r(xorl, var, reg)
#define	xorl_r2r(regs, regd)	ia32_r2r(xorl, regs, regd)
#define	xorl(vars, vard)	ia32_m2m(xorl, vars, vard)


/*	2x32, 4x16, and 8x8 Parallel CoMPare for EQuality
	(resulting fields are either 0 or -1)
*/
#define	pcmpeqd_m2r(var, reg)	ia32_m2r(pcmpeqd, var, reg)
#define	pcmpeqd_r2r(regs, regd)	ia32_r2r(pcmpeqd, regs, regd)
#define	pcmpeqd(vars, vard)	ia32_m2m(pcmpeqd, vars, vard)

#define	pcmpeqw_m2r(var, reg)	ia32_m2r(pcmpeqw, var, reg)
#define	pcmpeqw_r2r(regs, regd)	ia32_r2r(pcmpeqw, regs, regd)
#define	pcmpeqw(vars, vard)	ia32_m2m(pcmpeqw, vars, vard)

#define	pcmpeqb_m2r(var, reg)	ia32_m2r(pcmpeqb, var, reg)
#define	pcmpeqb_r2r(regs, regd)	ia32_r2r(pcmpeqb, regs, regd)
#define	pcmpeqb(vars, vard)	ia32_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)	ia32_m2r(pcmpgtd, var, reg)
#define	pcmpgtd_r2r(regs, regd)	ia32_r2r(pcmpgtd, regs, regd)
#define	pcmpgtd(vars, vard)	ia32_m2m(pcmpgtd, vars, vard)

#define	pcmpgtw_m2r(var, reg)	ia32_m2r(pcmpgtw, var, reg)
#define	pcmpgtw_r2r(regs, regd)	ia32_r2r(pcmpgtw, regs, regd)
#define	pcmpgtw(vars, vard)	ia32_m2m(pcmpgtw, vars, vard)

#define	pcmpgtb_m2r(var, reg)	ia32_m2r(pcmpgtb, var, reg)
#define	pcmpgtb_r2r(regs, regd)	ia32_r2r(pcmpgtb, regs, regd)
#define	pcmpgtb(vars, vard)	ia32_m2m(pcmpgtb, vars, vard)


/*	1x64, 2x32, and 4x16 Parallel Shift Left Logical
*/
#define	psllq_m2r(var, reg)	ia32_m2r(psllq, var, reg)
#define	psllq_i2r(immed, regd)	ia32_i2r(psllq, immed, regd)
#define	psllq_r2r(regs, regd)	ia32_r2r(psllq, regs, regd)
#define	psllq(vars, vard)	ia32_m2m(psllq, vars, vard)

#define	pslld_m2r(var, reg)	ia32_m2r(pslld, var, reg)
#define	pslld_i2r(immed, regd)	ia32_i2r(pslld, immed, regd)
#define	pslld_r2r(regs, regd)	ia32_r2r(pslld, regs, regd)
#define	pslld(vars, vard)	ia32_m2m(pslld, vars, vard)

#define	psllw_m2r(var, reg)	ia32_m2r(psllw, var, reg)
#define	psllw_i2r(immed, regd)	ia32_i2r(psllw, immed, regd)
#define	psllw_r2r(regs, regd)	ia32_r2r(psllw, regs, regd)
#define	psllw(vars, vard)	ia32_m2m(psllw, vars, vard)


/*	1x64, 2x32, and 4x16 Parallel Shift Right Logical
*/
#define	psrlq_m2r(var, reg)	ia32_m2r(psrlq, var, reg)
#define	psrlq_i2r(immed, regd)	ia32_i2r(psrlq, immed, regd)
#define	psrlq_r2r(regs, regd)	ia32_r2r(psrlq, regs, regd)
#define	psrlq(vars, vard)	ia32_m2m(psrlq, vars, vard)

#define	psrld_m2r(var, reg)	ia32_m2r(psrld, var, reg)
#define	psrld_i2r(immed, regd)	ia32_i2r(psrld, immed, regd)
#define	psrld_r2r(regs, regd)	ia32_r2r(psrld, regs, regd)
#define	psrld(vars, vard)	ia32_m2m(psrld, vars, vard)

#define	psrlw_m2r(var, reg)	ia32_m2r(psrlw, var, reg)
#define	psrlw_i2r(immed, regd)	ia32_i2r(psrlw, immed, regd)
#define	psrlw_r2r(regs, regd)	ia32_r2r(psrlw, regs, regd)
#define	psrlw(vars, vard)	ia32_m2m(psrlw, vars, vard)


/*	2x32 and 4x16 Parallel Shift Right Arithmetic
*/
#define	psrad_m2r(var, reg)	ia32_m2r(psrad, var, reg)
#define	psrad_i2r(immed, regd)	ia32_i2r(psrad, immed, regd)
#define	psrad_r2r(regs, regd)	ia32_r2r(psrad, regs, regd)
#define	psrad(vars, vard)	ia32_m2m(psrad, vars, vard)

#define	psraw_m2r(var, reg)	ia32_m2r(psraw, var, reg)
#define	psraw_i2r(immed, regd)	ia32_i2r(psraw, immed, regd)
#define	psraw_r2r(regs, regd)	ia32_r2r(psraw, regs, regd)
#define	psraw(vars, vard)	ia32_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)	ia32_m2r(packssdw, var, reg)
#define	packssdw_r2r(regs, regd) ia32_r2r(packssdw, regs, regd)
#define	packssdw(vars, vard)	ia32_m2m(packssdw, vars, vard)

#define	packsswb_m2r(var, reg)	ia32_m2r(packsswb, var, reg)
#define	packsswb_r2r(regs, regd) ia32_r2r(packsswb, regs, regd)
#define	packsswb(vars, vard)	ia32_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)	ia32_m2r(packuswb, var, reg)
#define	packuswb_r2r(regs, regd) ia32_r2r(packuswb, regs, regd)
#define	packuswb(vars, vard)	ia32_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)	ia32_m2r(punpckldq, var, reg)
#define	punpckldq_r2r(regs, regd) ia32_r2r(punpckldq, regs, regd)
#define	punpckldq(vars, vard)	ia32_m2m(punpckldq, vars, vard)

#define	punpcklwd_m2r(var, reg)	ia32_m2r(punpcklwd, var, reg)
#define	punpcklwd_r2r(regs, regd) ia32_r2r(punpcklwd, regs, regd)
#define	punpcklwd(vars, vard)	ia32_m2m(punpcklwd, vars, vard)

#define	punpcklbw_m2r(var, reg)	ia32_m2r(punpcklbw, var, reg)
#define	punpcklbw_r2r(regs, regd) ia32_r2r(punpcklbw, regs, regd)
#define	punpcklbw(vars, vard)	ia32_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)	ia32_m2r(punpckhdq, var, reg)
#define	punpckhdq_r2r(regs, regd) ia32_r2r(punpckhdq, regs, regd)
#define	punpckhdq(vars, vard)	ia32_m2m(punpckhdq, vars, vard)

#define	punpckhwd_m2r(var, reg)	ia32_m2r(punpckhwd, var, reg)
#define	punpckhwd_r2r(regs, regd) ia32_r2r(punpckhwd, regs, regd)
#define	punpckhwd(vars, vard)	ia32_m2m(punpckhwd, vars, vard)

#define	punpckhbw_m2r(var, reg)	ia32_m2r(punpckhbw, var, reg)
#define	punpckhbw_r2r(regs, regd) ia32_r2r(punpckhbw, regs, regd)
#define	punpckhbw(vars, vard)	ia32_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


/*	derived functions...
*/

#define	pdivb_m2r(mem, reg) \
	{ \
		p64_t _reg; \
		p64_t _mem = mem; \
		movq_r2m(reg, _reg); \
		if (_mem.b[0]) _reg.b[0] /= _mem.b[0]; \
		if (_mem.b[1]) _reg.b[1] /= _mem.b[1]; \
		if (_mem.b[2]) _reg.b[2] /= _mem.b[2]; \
		if (_mem.b[3]) _reg.b[3] /= _mem.b[3]; \
		if (_mem.b[4]) _reg.b[4] /= _mem.b[4]; \
		if (_mem.b[5]) _reg.b[5] /= _mem.b[5]; \
		if (_mem.b[6]) _reg.b[6] /= _mem.b[6]; \
		if (_mem.b[7]) _reg.b[7] /= _mem.b[7]; \
		movq_m2r(_reg, reg); \
	}

#define	pdivb_r2r(regs, regd) \
	{ \
		p64_t _regs; \
		movq_r2m(regs, _regs); \
		pdivb_m2r(_regs, regd); \
	}

#define	pdivw_m2r(mem, reg) \
	{ \
		p64_t _reg; \
		p64_t _mem = mem; \
		movq_r2m(reg, _reg); \
		if (_mem.w[0]) _reg.w[0] /= _mem.w[0]; \
		if (_mem.w[1]) _reg.w[1] /= _mem.w[1]; \
		if (_mem.w[2]) _reg.w[2] /= _mem.w[2]; \
		if (_mem.w[3]) _reg.w[3] /= _mem.w[3]; \
		movq_m2r(_reg, reg); \
	}

#define	pdivw_r2r(regs, regd) \
	{ \
		p64_t _regs; \
		movq_r2m(regs, _regs); \
		pdivw_m2r(_regs, regd); \
	}

#define	pdivd_m2r(mem, reg) \
	{ \
		p64_t _reg; \
		p64_t _mem = mem; \
		movq_r2m(reg, _reg); \
		if (_mem.d[0]) _reg.d[0] /= _mem.d[0]; \
		if (_mem.d[1]) _reg.d[1] /= _mem.d[1]; \
		movq_m2r(_reg, reg); \
	}

#define	pdivd_r2r(regs, regd) \
	{ \
		p64_t _regs; \
		movq_r2m(regs, _regs); \
		pdivd_m2r(_regs, regd); \
	}


#define	pmodb_m2r(mem, reg) \
	{ \
		p64_t _reg; \
		p64_t _mem = mem; \
		movq_r2m(reg, _reg); \
		if (_mem.b[0]) _reg.b[0] %= _mem.b[0]; \
		if (_mem.b[1]) _reg.b[1] %= _mem.b[1]; \
		if (_mem.b[2]) _reg.b[2] %= _mem.b[2]; \
		if (_mem.b[3]) _reg.b[3] %= _mem.b[3]; \
		if (_mem.b[4]) _reg.b[4] %= _mem.b[4]; \
		if (_mem.b[5]) _reg.b[5] %= _mem.b[5]; \
		if (_mem.b[6]) _reg.b[6] %= _mem.b[6]; \
		if (_mem.b[7]) _reg.b[7] %= _mem.b[7]; \
		movq_m2r(_reg, reg); \
	}

#define	pmodb_r2r(regs, regd) \
	{ \
		p64_t _regs; \
		movq_r2m(regs, _regs); \
		pmodb_m2r(_regs, regd); \
	}

#define	pmodw_m2r(mem, reg) \
	{ \
		p64_t _reg; \
		p64_t _mem = mem; \
		movq_r2m(reg, _reg); \
		if (_mem.w[0]) _reg.w[0] %= _mem.w[0]; \
		if (_mem.w[1]) _reg.w[1] %= _mem.w[1]; \
		if (_mem.w[2]) _reg.w[2] %= _mem.w[2]; \
		if (_mem.w[3]) _reg.w[3] %= _mem.w[3]; \
		movq_m2r(_reg, reg); \
	}

#define	pmodw_r2r(regs, regd) \
	{ \
		p64_t _regs; \
		movq_r2m(regs, _regs); \
		pmodw_m2r(_regs, regd); \
	}

#define	pmodd_m2r(mem, reg) \
	{ \
		p64_t _reg; \
		p64_t _mem = mem; \
		movq_r2m(reg, _reg); \
		if (_mem.d[0]) _reg.d[0] %= _mem.d[0]; \
		if (_mem.d[1]) _reg.d[1] %= _mem.d[1]; \
		movq_m2r(_reg, reg); \
	}

#define	pmodd_r2r(regs, regd) \
	{ \
		p64_t _regs; \
		movq_r2m(regs, _regs); \
		pmodd_m2r(_regs, regd); \
	}

#define	pmulld_m2r(mem, reg) \
	{ \
		p64_t _reg; \
		p64_t _mem = mem; \
		movq_r2m(reg, _reg); \
		_reg.d[0] *= _mem.d[0]; \
		_reg.d[1] *= _mem.d[1]; \
		movq_m2r(_reg, reg); \
	}

#define	pmulld_r2r(regs, regd) \
	{ \
		p64_t _regs; \
		movq_r2m(regs, _regs); \
		pmulld_m2r(_regs, regd); \
	}
