/*********************************************************************/
/*                                                                   */
/*             Optimized BLAS libraries                              */
/*                     By Kazushige Goto <kgoto@tacc.utexas.edu>     */
/*                                                                   */
/* Copyright (c) The University of Texas, 2009. All rights reserved. */
/* UNIVERSITY EXPRESSLY DISCLAIMS ANY AND ALL WARRANTIES CONCERNING  */
/* THIS SOFTWARE AND DOCUMENTATION, INCLUDING ANY WARRANTIES OF      */
/* MERCHANTABILITY, FITNESS FOR ANY PARTICULAR PURPOSE,              */
/* NON-INFRINGEMENT AND WARRANTIES OF PERFORMANCE, AND ANY WARRANTY  */
/* THAT MIGHT OTHERWISE ARISE FROM COURSE OF DEALING OR USAGE OF     */
/* TRADE. NO WARRANTY IS EITHER EXPRESS OR IMPLIED WITH RESPECT TO   */
/* THE USE OF THE SOFTWARE OR DOCUMENTATION.                         */
/* Under no circumstances shall University be liable for incidental, */
/* special, indirect, direct or consequential damages or loss of     */
/* profits, interruption of business, or related expenses which may  */
/* arise from use of Software or Documentation, including but not    */
/* limited to those resulting from defects in Software and/or        */
/* Documentation, or loss or inaccuracy of data of any kind.         */
/*********************************************************************/

#define ASSEMBLER
#include "common.h"
#include "l2param.h"

#if GEMV_UNROLL < 2
#undef  GEMV_UNROLL
#define GEMV_UNROLL 2
#endif

#ifndef WINDOWS_ABI

#define STACKSIZE	64
	
#define OLD_INCX	 8 + STACKSIZE(%rsp)
#define OLD_Y		16 + STACKSIZE(%rsp)
#define OLD_INCY	24 + STACKSIZE(%rsp)
#define OLD_BUFFER	32 + STACKSIZE(%rsp)
#define ALPHA		48	      (%rsp)
	
#define M	  %rdi
#define N	  %rsi
#define A	  %rcx
#define LDA	  %r8
#define X	  %r9
#define INCX	  %rdx
#define Y	  %rbp
#define INCY	  %r10

#else

#define STACKSIZE	256
	
#define OLD_ALPHA_I	 40 + STACKSIZE(%rsp)
#define OLD_A		 48 + STACKSIZE(%rsp)
#define OLD_LDA		 56 + STACKSIZE(%rsp)
#define OLD_X		 64 + STACKSIZE(%rsp)
#define OLD_INCX	 72 + STACKSIZE(%rsp)
#define OLD_Y		 80 + STACKSIZE(%rsp)
#define OLD_INCY	 88 + STACKSIZE(%rsp)
#define OLD_BUFFER	 96 + STACKSIZE(%rsp)
#define ALPHA		224	       (%rsp)

#define M	  %rcx
#define N	  %rdx
#define A	  %r8
#define LDA	  %r9
#define X	  %rdi
#define INCX	  %rsi
#define Y	  %rbp
#define INCY	  %r10

#endif

#define I	%rax
#define A1	%r11
#define A2	%r12

#define X1	%rbx
#define Y1	%r13
#define BUFFER	%r14

#ifdef ALIGNED_ACCESS
#define MM	%r15
#else
#define MM	M
#endif

#undef SUBPS

#ifndef CONJ
#define SUBPS	   addps
#else
#define SUBPS	   subps
#endif

	PROLOGUE
	PROFCODE

	subq	$STACKSIZE, %rsp
	movq	%rbx,  0(%rsp)
	movq	%rbp,  8(%rsp)
	movq	%r12, 16(%rsp)
	movq	%r13, 24(%rsp)
	movq	%r14, 32(%rsp)
	movq	%r15, 40(%rsp)

#ifdef WINDOWS_ABI
	movq	%rdi,    48(%rsp)
	movq	%rsi,    56(%rsp)
	movups	%xmm6,   64(%rsp)
	movups	%xmm7,   80(%rsp)
	movups	%xmm8,   96(%rsp)
	movups	%xmm9,  112(%rsp)
	movups	%xmm10, 128(%rsp)
	movups	%xmm11, 144(%rsp)
	movups	%xmm12, 160(%rsp)
	movups	%xmm13, 176(%rsp)
	movups	%xmm14, 192(%rsp)
	movups	%xmm15, 208(%rsp)

	movq	OLD_A,     A
	movq	OLD_LDA,   LDA
	movq	OLD_X,     X

	movaps	%xmm3,       %xmm0
	movss	OLD_ALPHA_I, %xmm1
#endif

	movq	OLD_INCX,  INCX
	movq	OLD_Y,     Y
	movq	OLD_INCY,  INCY
	movq	OLD_BUFFER, BUFFER

	salq	$ZBASE_SHIFT,   LDA
	salq	$ZBASE_SHIFT,   INCX
	salq	$ZBASE_SHIFT,   INCY

	unpcklps %xmm1, %xmm0

	movlps	 %xmm0, ALPHA

	testq	M, M
	jle	.L999
	testq	N, N
	jle	.L999

	subq	$-32 * SIZE, A
	
	movq	BUFFER, X1

#ifdef ALIGNED_ACCESS
	movq	M, MM
	movq	A, %rax
	andq	$4 * SIZE - 1, %rax
	cmpq	$2 * SIZE, %rax

	jl	.L0X

	movsd	(X), %xmm0
	addq	INCX, X
	movlps	%xmm0, 2 * SIZE(X1)

	addq	$2 * SIZE, BUFFER
 	addq	$4 * SIZE, X1
	decq	MM

.L0X:
#endif

	movq	MM, I
	sarq	$3, I
	jle	.L05
	ALIGN_4

.L02:
	movsd	(X), %xmm0
	addq	INCX, X
	movhps	(X), %xmm0
	addq	INCX, X

	movsd	(X), %xmm1
	addq	INCX, X
	movhps	(X), %xmm1
	addq	INCX, X

	movsd	(X), %xmm2
	addq	INCX, X
	movhps	(X), %xmm2
	addq	INCX, X

	movsd	(X), %xmm3
	addq	INCX, X
	movhps	(X), %xmm3
	addq	INCX, X

	movaps	%xmm0,  0 * SIZE(X1)
	movaps	%xmm1,  4 * SIZE(X1)
	movaps	%xmm2,  8 * SIZE(X1)
	movaps	%xmm3, 12 * SIZE(X1)

	addq	$16 * SIZE, X1
	decq	I
	jg	.L02
	ALIGN_4

.L05:
	movq	MM, I
	andq	$7, I
	jle	.L10
	ALIGN_2

.L06:
	movsd	(X), %xmm0
	addq	INCX, X
	movlps	%xmm0, 0 * SIZE(X1)
	addq	$2 * SIZE, X1
	decq	I
	jg	.L06
	ALIGN_4

.L10:
	movq	Y, Y1

#ifdef ALIGNED_ACCESS
	testq	$SIZE, A
	jne	.L200

	testq	$2 * SIZE, LDA
	jne	.L100
#endif

#if GEMV_UNROLL >= 4
	cmpq	$4, N
	jl	.L20
	ALIGN_3

.L11:
	subq	$4, N

	leaq	32 * SIZE(BUFFER), X1

	movq	A, A1
	leaq	(A1, LDA, 2), A2
	leaq	(A1, LDA, 4), A

	xorps	%xmm0, %xmm0
	xorps	%xmm1, %xmm1
	xorps	%xmm2, %xmm2
	xorps	%xmm3, %xmm3
	xorps	%xmm4, %xmm4
	xorps	%xmm5, %xmm5
	xorps	%xmm6, %xmm6
	xorps	%xmm7, %xmm7

#ifdef ALIGNED_ACCESS
	cmpq	M, MM
	je	.L1X

#ifdef movsd
	xorps	%xmm8, %xmm8
#endif
	movsd	-32 * SIZE(A1), %xmm8
#ifdef movsd
	xorps	%xmm9, %xmm9
#endif
	movsd	-32 * SIZE(A1, LDA), %xmm9
#ifdef movsd
	xorps	%xmm10, %xmm10
#endif
	movsd	-32 * SIZE(A2), %xmm10
#ifdef movsd
	xorps	%xmm11, %xmm11
#endif
	movsd	-32 * SIZE(A2, LDA), %xmm11

#ifdef movsd
	xorps	%xmm12, %xmm12
#endif
	movsd	-32 * SIZE(X1), %xmm12

	pshufd	$0xb1,  %xmm8, %xmm14
	mulps	%xmm12, %xmm8
	addps	%xmm8,  %xmm0
	pshufd	$0xb1,  %xmm9, %xmm15
	mulps	%xmm12, %xmm9
	addps	%xmm9,  %xmm2

	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm1
	mulps	%xmm12, %xmm15
	SUBPS	%xmm15, %xmm3

	pshufd	$0xb1,  %xmm10, %xmm14
	mulps	%xmm12, %xmm10
	addps	%xmm10, %xmm4
	pshufd	$0xb1,  %xmm11, %xmm15
	mulps	%xmm12, %xmm11
	addps	%xmm11, %xmm6

	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm5
	mulps	%xmm12, %xmm15
	SUBPS	%xmm15, %xmm7

	addq	 $2 * SIZE, A1
	addq	 $2 * SIZE, A2
	addq	 $2 * SIZE, X1
	ALIGN_3
.L1X:
#endif

	movaps	-32 * SIZE(X1), %xmm12
	movaps	-28 * SIZE(X1), %xmm13

#ifdef PREFETCHW
	PREFETCHW	7 * SIZE(Y1)
#endif

	movq	MM,  I
	sarq	$3,  I
	jle	.L15

	MOVUPS_A1(-32 * SIZE, A1, %xmm8)
	MOVUPS_A2(-32 * SIZE, A1, LDA, 1, %xmm9)
	MOVUPS_A1(-32 * SIZE, A2, %xmm10)
	MOVUPS_A2(-32 * SIZE, A2, LDA, 1, %xmm11)

	decq	 I
	jle	 .L14
	ALIGN_3

.L13:
#ifdef PREFETCH
	PREFETCH	(PREFETCHSIZE) - 128 + PREOFFSET(A1)
#endif

	pshufd	$0xb1,  %xmm8, %xmm14
	mulps	%xmm12, %xmm8
	addps	%xmm8,  %xmm0
	MOVUPS_A1(-28 * SIZE, A1, %xmm8)
	pshufd	$0xb1,  %xmm9, %xmm15
	mulps	%xmm12, %xmm9
	addps	%xmm9,  %xmm2
	MOVUPS_A2(-28 * SIZE, A1, LDA, 1, %xmm9)

	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm1
	mulps	%xmm12, %xmm15
	SUBPS	%xmm15, %xmm3

	pshufd	$0xb1,  %xmm10, %xmm14
	mulps	%xmm12, %xmm10
	addps	%xmm10, %xmm4
	MOVUPS_A1(-28 * SIZE, A2, %xmm10)
	pshufd	$0xb1,  %xmm11, %xmm15
	mulps	%xmm12, %xmm11
	addps	%xmm11, %xmm6
	MOVUPS_A2(-28 * SIZE, A2, LDA, 1, %xmm11)

	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm5
	mulps	%xmm12, %xmm15
	movaps	-24 * SIZE(X1), %xmm12
	SUBPS	%xmm15, %xmm7

#ifdef PREFETCH
	PREFETCH	(PREFETCHSIZE) - 128 + PREOFFSET(A1, LDA)
#endif

	pshufd	$0xb1,  %xmm8, %xmm14
	mulps	%xmm13, %xmm8
	addps	%xmm8,  %xmm0
	MOVUPS_A1(-24 * SIZE, A1, %xmm8)
	pshufd	$0xb1,  %xmm9, %xmm15
	mulps	%xmm13, %xmm9
	addps	%xmm9,  %xmm2
	MOVUPS_A2(-24 * SIZE, A1, LDA, 1, %xmm9)

	mulps	%xmm13, %xmm14
	SUBPS	%xmm14, %xmm1
	mulps	%xmm13, %xmm15
	SUBPS	%xmm15, %xmm3

	pshufd	$0xb1,  %xmm10, %xmm14
	mulps	%xmm13, %xmm10
	addps	%xmm10, %xmm4
	MOVUPS_A1(-24 * SIZE, A2, %xmm10)
	pshufd	$0xb1,  %xmm11, %xmm15
	mulps	%xmm13, %xmm11
	addps	%xmm11, %xmm6
	MOVUPS_A2(-24 * SIZE, A2, LDA, 1, %xmm11)

	mulps	%xmm13, %xmm14
	SUBPS	%xmm14, %xmm5
	mulps	%xmm13, %xmm15
	movaps	-20 * SIZE(X1), %xmm13
	SUBPS	%xmm15, %xmm7

#ifdef PREFETCH
	PREFETCH	(PREFETCHSIZE) - 128 + PREOFFSET(A2)
#endif

	pshufd	$0xb1,  %xmm8, %xmm14
	mulps	%xmm12, %xmm8
	addps	%xmm8,  %xmm0
	MOVUPS_A1(-20 * SIZE, A1, %xmm8)
	pshufd	$0xb1,  %xmm9, %xmm15
	mulps	%xmm12, %xmm9
	addps	%xmm9,  %xmm2
	MOVUPS_A2(-20 * SIZE, A1, LDA, 1, %xmm9)

	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm1
	mulps	%xmm12, %xmm15
	SUBPS	%xmm15, %xmm3

	pshufd	$0xb1,  %xmm10, %xmm14
	mulps	%xmm12, %xmm10
	addps	%xmm10, %xmm4
	MOVUPS_A1(-20 * SIZE, A2, %xmm10)
	pshufd	$0xb1,  %xmm11, %xmm15
	mulps	%xmm12, %xmm11
	addps	%xmm11, %xmm6
	MOVUPS_A2(-20 * SIZE, A2, LDA, 1, %xmm11)

	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm5
	mulps	%xmm12, %xmm15
	movaps	-16 * SIZE(X1), %xmm12
	SUBPS	%xmm15, %xmm7

#ifdef PREFETCH
	PREFETCH	(PREFETCHSIZE) - 128 + PREOFFSET(A2, LDA)
#endif

	pshufd	$0xb1,  %xmm8, %xmm14
	mulps	%xmm13, %xmm8
	addps	%xmm8,  %xmm0
	MOVUPS_A1(-16 * SIZE, A1, %xmm8)
	pshufd	$0xb1,  %xmm9, %xmm15
	mulps	%xmm13, %xmm9
	addps	%xmm9,  %xmm2
	MOVUPS_A2(-16 * SIZE, A1, LDA, 1, %xmm9)

	mulps	%xmm13, %xmm14
	SUBPS	%xmm14, %xmm1
	mulps	%xmm13, %xmm15
	SUBPS	%xmm15, %xmm3

#ifdef PREFETCHW
	PREFETCH	(PREFETCHSIZE) - 128 + PREOFFSET(X1)
#endif

	pshufd	$0xb1,  %xmm10, %xmm14
	mulps	%xmm13, %xmm10
	addps	%xmm10, %xmm4
	MOVUPS_A1(-16 * SIZE, A2, %xmm10)
	pshufd	$0xb1,  %xmm11, %xmm15
	mulps	%xmm13, %xmm11
	addps	%xmm11, %xmm6
	MOVUPS_A2(-16 * SIZE, A2, LDA, 1, %xmm11)

	mulps	%xmm13, %xmm14
	SUBPS	%xmm14, %xmm5
	mulps	%xmm13, %xmm15
	movaps	-12 * SIZE(X1), %xmm13
	SUBPS	%xmm15, %xmm7

	subq	 $-16 * SIZE, A1
	subq	 $-16 * SIZE, A2
	subq	 $-16 * SIZE, X1

	subq	 $1, I
	BRANCH
	jg	.L13
	ALIGN_3

.L14:
	pshufd	$0xb1,  %xmm8, %xmm14
	mulps	%xmm12, %xmm8
	addps	%xmm8,  %xmm0
	MOVUPS_A1(-28 * SIZE, A1, %xmm8)
	pshufd	$0xb1,  %xmm9, %xmm15
	mulps	%xmm12, %xmm9
	addps	%xmm9,  %xmm2
	MOVUPS_A2(-28 * SIZE, A1, LDA, 1, %xmm9)

	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm1
	mulps	%xmm12, %xmm15
	SUBPS	%xmm15, %xmm3

	pshufd	$0xb1,  %xmm10, %xmm14
	mulps	%xmm12, %xmm10
	addps	%xmm10, %xmm4
	MOVUPS_A1(-28 * SIZE, A2, %xmm10)
	pshufd	$0xb1,  %xmm11, %xmm15
	mulps	%xmm12, %xmm11
	addps	%xmm11, %xmm6
	MOVUPS_A2(-28 * SIZE, A2, LDA, 1, %xmm11)

	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm5
	mulps	%xmm12, %xmm15
	movaps	-24 * SIZE(X1), %xmm12
	SUBPS	%xmm15, %xmm7

	pshufd	$0xb1,  %xmm8, %xmm14
	mulps	%xmm13, %xmm8
	addps	%xmm8,  %xmm0
	MOVUPS_A1(-24 * SIZE, A1, %xmm8)
	pshufd	$0xb1,  %xmm9, %xmm15
	mulps	%xmm13, %xmm9
	addps	%xmm9,  %xmm2
	MOVUPS_A2(-24 * SIZE, A1, LDA, 1, %xmm9)

	mulps	%xmm13, %xmm14
	SUBPS	%xmm14, %xmm1
	mulps	%xmm13, %xmm15
	SUBPS	%xmm15, %xmm3

	pshufd	$0xb1,  %xmm10, %xmm14
	mulps	%xmm13, %xmm10
	addps	%xmm10, %xmm4
	MOVUPS_A1(-24 * SIZE, A2, %xmm10)
	pshufd	$0xb1,  %xmm11, %xmm15
	mulps	%xmm13, %xmm11
	addps	%xmm11, %xmm6
	MOVUPS_A2(-24 * SIZE, A2, LDA, 1, %xmm11)

	mulps	%xmm13, %xmm14
	SUBPS	%xmm14, %xmm5
	mulps	%xmm13, %xmm15
	movaps	-20 * SIZE(X1), %xmm13
	SUBPS	%xmm15, %xmm7

	pshufd	$0xb1,  %xmm8, %xmm14
	mulps	%xmm12, %xmm8
	addps	%xmm8,  %xmm0
	MOVUPS_A1(-20 * SIZE, A1, %xmm8)
	pshufd	$0xb1,  %xmm9, %xmm15
	mulps	%xmm12, %xmm9
	addps	%xmm9,  %xmm2
	MOVUPS_A2(-20 * SIZE, A1, LDA, 1, %xmm9)

	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm1
	mulps	%xmm12, %xmm15
	SUBPS	%xmm15, %xmm3

	pshufd	$0xb1,  %xmm10, %xmm14
	mulps	%xmm12, %xmm10
	addps	%xmm10, %xmm4
	MOVUPS_A1(-20 * SIZE, A2, %xmm10)
	pshufd	$0xb1,  %xmm11, %xmm15
	mulps	%xmm12, %xmm11
	addps	%xmm11, %xmm6
	MOVUPS_A2(-20 * SIZE, A2, LDA, 1, %xmm11)

	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm5
	mulps	%xmm12, %xmm15
	movaps	-16 * SIZE(X1), %xmm12
	SUBPS	%xmm15, %xmm7

	pshufd	$0xb1,  %xmm8, %xmm14
	mulps	%xmm13, %xmm8
	addps	%xmm8,  %xmm0
	pshufd	$0xb1,  %xmm9, %xmm15
	mulps	%xmm13, %xmm9
	addps	%xmm9,  %xmm2

	mulps	%xmm13, %xmm14
	SUBPS	%xmm14, %xmm1
	mulps	%xmm13, %xmm15
	SUBPS	%xmm15, %xmm3

	pshufd	$0xb1,  %xmm10, %xmm14
	mulps	%xmm13, %xmm10
	addps	%xmm10, %xmm4
	pshufd	$0xb1,  %xmm11, %xmm15
	mulps	%xmm13, %xmm11
	addps	%xmm11, %xmm6

	mulps	%xmm13, %xmm14
	SUBPS	%xmm14, %xmm5
	mulps	%xmm13, %xmm15
	movaps	-12 * SIZE(X1), %xmm13
	SUBPS	%xmm15, %xmm7

	subq	 $-16 * SIZE, A1
	subq	 $-16 * SIZE, A2
	subq	 $-16 * SIZE, X1
	ALIGN_3

.L15:
	testq	$4, MM
	je	.L17

	MOVUPS_A1(-32 * SIZE, A1, %xmm8)
	MOVUPS_A2(-32 * SIZE, A1, LDA, 1, %xmm9)
	MOVUPS_A1(-32 * SIZE, A2, %xmm10)
	MOVUPS_A2(-32 * SIZE, A2, LDA, 1, %xmm11)

	pshufd	$0xb1,  %xmm8, %xmm14
	mulps	%xmm12, %xmm8
	addps	%xmm8,  %xmm0
	MOVUPS_A1(-28 * SIZE, A1, %xmm8)
	pshufd	$0xb1,  %xmm9, %xmm15
	mulps	%xmm12, %xmm9
	addps	%xmm9,  %xmm2
	MOVUPS_A2(-28 * SIZE, A1, LDA, 1, %xmm9)

	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm1
	mulps	%xmm12, %xmm15
	SUBPS	%xmm15, %xmm3

	pshufd	$0xb1,  %xmm10, %xmm14
	mulps	%xmm12, %xmm10
	addps	%xmm10, %xmm4
	MOVUPS_A1(-28 * SIZE, A2, %xmm10)
	pshufd	$0xb1,  %xmm11, %xmm15
	mulps	%xmm12, %xmm11
	addps	%xmm11, %xmm6
	MOVUPS_A2(-28 * SIZE, A2, LDA, 1, %xmm11)

	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm5
	mulps	%xmm12, %xmm15
	movaps	-24 * SIZE(X1), %xmm12
	SUBPS	%xmm15, %xmm7

	pshufd	$0xb1,  %xmm8, %xmm14
	mulps	%xmm13, %xmm8
	addps	%xmm8,  %xmm0
	pshufd	$0xb1,  %xmm9, %xmm15
	mulps	%xmm13, %xmm9
	addps	%xmm9,  %xmm2

	mulps	%xmm13, %xmm14
	SUBPS	%xmm14, %xmm1
	mulps	%xmm13, %xmm15
	SUBPS	%xmm15, %xmm3

	pshufd	$0xb1,  %xmm10, %xmm14
	mulps	%xmm13, %xmm10
	addps	%xmm10, %xmm4
	pshufd	$0xb1,  %xmm11, %xmm15
	mulps	%xmm13, %xmm11
	addps	%xmm11, %xmm6

	mulps	%xmm13, %xmm14
	SUBPS	%xmm14, %xmm5
	mulps	%xmm13, %xmm15
	movaps	-20 * SIZE(X1), %xmm13
	SUBPS	%xmm15, %xmm7

	addq	 $8 * SIZE, A1
	addq	 $8 * SIZE, A2
	addq	 $8 * SIZE, X1
	ALIGN_3

.L17:
	testq	$2, MM
	je	.L18

	MOVUPS_A1(-32 * SIZE, A1, %xmm8)
	MOVUPS_A2(-32 * SIZE, A1, LDA, 1, %xmm9)
	MOVUPS_A1(-32 * SIZE, A2, %xmm10)
	MOVUPS_A2(-32 * SIZE, A2, LDA, 1, %xmm11)

	pshufd	$0xb1,  %xmm8, %xmm14
	mulps	%xmm12, %xmm8
	addps	%xmm8,  %xmm0
	pshufd	$0xb1,  %xmm9, %xmm15
	mulps	%xmm12, %xmm9
	addps	%xmm9,  %xmm2

	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm1
	mulps	%xmm12, %xmm15
	SUBPS	%xmm15, %xmm3

	pshufd	$0xb1,  %xmm10, %xmm14
	mulps	%xmm12, %xmm10
	addps	%xmm10, %xmm4
	pshufd	$0xb1,  %xmm11, %xmm15
	mulps	%xmm12, %xmm11
	addps	%xmm11, %xmm6

	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm5
	mulps	%xmm12, %xmm15
	SUBPS	%xmm15, %xmm7

	movaps	%xmm13, %xmm12

	addq	 $4 * SIZE, A1
	addq	 $4 * SIZE, A2
	ALIGN_3

.L18:
	testq	$1, MM
	je	.L19

#ifdef movsd
	xorps	%xmm8, %xmm8
#endif
	movsd	-32 * SIZE(A1), %xmm8
#ifdef movsd
	xorps	%xmm9, %xmm9
#endif
	movsd	-32 * SIZE(A1, LDA), %xmm9
#ifdef movsd
	xorps	%xmm10, %xmm10
#endif
	movsd	-32 * SIZE(A2), %xmm10
#ifdef movsd
	xorps	%xmm11, %xmm11
#endif
	movsd	-32 * SIZE(A2, LDA), %xmm11

	pshufd	$0xb1,  %xmm8, %xmm14
	mulps	%xmm12, %xmm8
	addps	%xmm8,  %xmm0
	pshufd	$0xb1,  %xmm9, %xmm15
	mulps	%xmm12, %xmm9
	addps	%xmm9,  %xmm2

	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm1
	mulps	%xmm12, %xmm15
	SUBPS	%xmm15, %xmm3

	pshufd	$0xb1,  %xmm10, %xmm14
	mulps	%xmm12, %xmm10
	addps	%xmm10, %xmm4
	pshufd	$0xb1,  %xmm11, %xmm15
	mulps	%xmm12, %xmm11
	addps	%xmm11, %xmm6

	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm5
	mulps	%xmm12, %xmm15
	SUBPS	%xmm15, %xmm7
	ALIGN_3

.L19:
	pcmpeqb	%xmm11, %xmm11
	psllq	$63,    %xmm11

#if (!defined(CONJ) && !defined(XCONJ)) || (defined(CONJ) && defined(XCONJ))
	xorps	%xmm11, %xmm0
	xorps	%xmm11, %xmm2
	xorps	%xmm11, %xmm4
	xorps	%xmm11, %xmm6
#else
	xorps	%xmm11, %xmm1
	xorps	%xmm11, %xmm3
	xorps	%xmm11, %xmm5
	xorps	%xmm11, %xmm7
#endif

#ifdef HAVE_SSE3
	haddps	%xmm1, %xmm0
	haddps	%xmm3, %xmm2
	haddps	%xmm2, %xmm0

	haddps	%xmm5, %xmm4
	haddps	%xmm7, %xmm6
	haddps	%xmm6, %xmm4
#else

	movaps	%xmm0, %xmm8
	unpcklps %xmm1, %xmm0
	unpckhps %xmm1, %xmm8

	movaps	%xmm2, %xmm9
	unpcklps %xmm3, %xmm2
	unpckhps %xmm3, %xmm9

	movaps	%xmm4, %xmm10
	unpcklps %xmm5, %xmm4
	unpckhps %xmm5, %xmm10

	movaps	%xmm6, %xmm12
	unpcklps %xmm7, %xmm6
	unpckhps %xmm7, %xmm12

	addps	%xmm8,  %xmm0
	addps	%xmm9,  %xmm2
	addps	%xmm10, %xmm4
	addps	%xmm12, %xmm6

	movhlps	%xmm0, %xmm1
	movhlps	%xmm2, %xmm3
	movhlps	%xmm4, %xmm5
	movhlps	%xmm6, %xmm7

	addps	%xmm1, %xmm0
	addps	%xmm3, %xmm2
	addps	%xmm5, %xmm4
	addps	%xmm7, %xmm6

	movlhps %xmm2, %xmm0
	movlhps %xmm6, %xmm4
#endif

	pshufd	$0xb1, %xmm0, %xmm1
	pshufd	$0xb1, %xmm4, %xmm5

#ifdef HAVE_SSE3
	movddup	ALPHA, %xmm15
#else
	movsd	ALPHA, %xmm15
	pshufd	$0x44, %xmm15, %xmm15
#endif

	mulps	%xmm15, %xmm0
	mulps	%xmm15, %xmm1
	mulps	%xmm15, %xmm4
	mulps	%xmm15, %xmm5

	xorps	%xmm11, %xmm0
	xorps	%xmm11, %xmm4

#ifdef HAVE_SSE3
	haddps	%xmm1, %xmm0
	haddps	%xmm5, %xmm4
#else
	movaps	%xmm0, %xmm2
	shufps	$0x88, %xmm1, %xmm0
	shufps	$0xdd, %xmm1, %xmm2

	movaps	%xmm4, %xmm6
	shufps	$0x88, %xmm5, %xmm4
	shufps	$0xdd, %xmm5, %xmm6
	
	addps	%xmm2, %xmm0
	addps	%xmm6, %xmm4
#endif

	movsd	(Y), %xmm2
	addq	INCY, Y
	movhps	(Y), %xmm2
	addq	INCY, Y
	movsd	(Y), %xmm6
	addq	INCY, Y
	movhps	(Y), %xmm6
	addq	INCY, Y

	shufps	$0xd8, %xmm0, %xmm0
	shufps	$0xd8, %xmm4, %xmm4

	addps	%xmm2, %xmm0
	addps	%xmm6, %xmm4

	movlps	%xmm0, (Y1)
	addq	INCY, Y1
	movhps	%xmm0, (Y1)
	addq	INCY, Y1
	movlps	%xmm4, (Y1)
	addq	INCY, Y1
	movhps	%xmm4, (Y1)
	addq	INCY, Y1

	cmpq	$4, N
	jge	.L11
	ALIGN_3

.L20:
#endif

	cmpq	$2, N
	jl	.L30
#if GEMV_UNROLL == 2
	ALIGN_3

.L21:
#endif
	subq	$2, N

	leaq	32 * SIZE(BUFFER), X1

	movq	A, A1
	leaq	(A1, LDA), A2
	leaq	(A1, LDA, 2), A

	xorps	%xmm0, %xmm0
	xorps	%xmm1, %xmm1
	xorps	%xmm2, %xmm2
	xorps	%xmm3, %xmm3

#ifdef ALIGNED_ACCESS
	cmpq	M, MM
	je	.L2X

#ifdef movsd
	xorps	%xmm8, %xmm8
#endif
	movsd	-32 * SIZE(A1), %xmm8
#ifdef movsd
	xorps	%xmm9, %xmm9
#endif
	movsd	-32 * SIZE(A2), %xmm9

#ifdef movsd
	xorps	%xmm12, %xmm12
#endif
	movsd	-32 * SIZE(X1), %xmm12

	pshufd	$0xb1, %xmm8,  %xmm4
	mulps	%xmm12,  %xmm8
	addps	%xmm8,  %xmm0
	mulps	%xmm12,  %xmm4
	SUBPS	%xmm4,  %xmm1

	pshufd	$0xb1, %xmm9, %xmm5
	mulps	%xmm12,  %xmm9
	addps	%xmm9, %xmm2
	mulps	%xmm12,  %xmm5
	SUBPS	%xmm5, %xmm3

	addq	 $2 * SIZE, A1
	addq	 $2 * SIZE, A2
	addq	 $2 * SIZE, X1
	ALIGN_3
.L2X:
#endif

	movaps	-32 * SIZE(X1), %xmm12
	movaps	-28 * SIZE(X1), %xmm13

#if (GEMV_UNROLL == 2) && defined(PREFETCHW)
	PREFETCHW	3 * SIZE(Y1)
#endif

	movq	MM,  I
	sarq	$3,  I
	jle	.L25

	MOVUPS_A1(-32 * SIZE, A1, %xmm8)
	MOVUPS_A1(-32 * SIZE, A2, %xmm9)
	MOVUPS_A1(-28 * SIZE, A1, %xmm10)
	MOVUPS_A1(-28 * SIZE, A2, %xmm11)

	decq	 I
	jle	 .L24
	ALIGN_3

.L23:
#ifdef PREFETCH
	PREFETCH	(PREFETCHSIZE) * 2 - 128 + PREOFFSET(A1)
#endif

	pshufd	$0xb1, %xmm8,  %xmm4
	mulps	%xmm12,  %xmm8
	addps	%xmm8,  %xmm0
	MOVUPS_A1(-24 * SIZE, A1, %xmm8)
	mulps	%xmm12,  %xmm4
	SUBPS	%xmm4,  %xmm1

	pshufd	$0xb1, %xmm9, %xmm5
	mulps	%xmm12,  %xmm9
	addps	%xmm9, %xmm2
	MOVUPS_A1(-24 * SIZE, A2, %xmm9)
	mulps	%xmm12,  %xmm5
	SUBPS	%xmm5, %xmm3

	movaps	-24 * SIZE(X1), %xmm12

	pshufd	$0xb1, %xmm10,  %xmm6
	mulps	%xmm13,  %xmm10
	addps	%xmm10,  %xmm0
	MOVUPS_A1(-20 * SIZE, A1, %xmm10)
	mulps	%xmm13,  %xmm6
	SUBPS	%xmm6,  %xmm1

	pshufd	$0xb1, %xmm11, %xmm7
	mulps	%xmm13,  %xmm11
	addps	%xmm11, %xmm2
	MOVUPS_A1(-20 * SIZE, A2, %xmm11)
	mulps	%xmm13,  %xmm7
	SUBPS	%xmm7, %xmm3

	movaps	-20 * SIZE(X1), %xmm13

#ifdef PREFETCH
	PREFETCH	(PREFETCHSIZE) * 2 - 128 + PREOFFSET(A2)
#endif

	pshufd	$0xb1, %xmm8,  %xmm4
	mulps	%xmm12,  %xmm8
	addps	%xmm8,  %xmm0
	MOVUPS_A1(-16 * SIZE, A1, %xmm8)
	mulps	%xmm12,  %xmm4
	SUBPS	%xmm4,  %xmm1

	pshufd	$0xb1, %xmm9, %xmm5
	mulps	%xmm12,  %xmm9
	addps	%xmm9, %xmm2
	MOVUPS_A1(-16 * SIZE, A2, %xmm9)
	mulps	%xmm12,  %xmm5
	SUBPS	%xmm5, %xmm3

	movaps	-16 * SIZE(X1), %xmm12

#ifdef PREFETCHW
	PREFETCH	(PREFETCHSIZE) * 2 - 128 + PREOFFSET(X1)
#endif

	pshufd	$0xb1, %xmm10,  %xmm6
	mulps	%xmm13,  %xmm10
	addps	%xmm10,  %xmm0
	MOVUPS_A1(-12 * SIZE, A1, %xmm10)
	mulps	%xmm13,  %xmm6
	SUBPS	%xmm6,  %xmm1

	pshufd	$0xb1, %xmm11, %xmm7
	mulps	%xmm13,  %xmm11
	addps	%xmm11, %xmm2
	MOVUPS_A1(-12 * SIZE, A2, %xmm11)
	mulps	%xmm13,  %xmm7
	SUBPS	%xmm7, %xmm3

	movaps	-12 * SIZE(X1), %xmm13

	subq	 $-16 * SIZE, A1
	subq	 $-16 * SIZE, A2
	subq	 $-16 * SIZE, X1

	subq	 $1, I
	BRANCH
	jg	.L23
	ALIGN_3

.L24:
	pshufd	$0xb1, %xmm8,  %xmm4
	mulps	%xmm12,  %xmm8
	addps	%xmm8,  %xmm0
	MOVUPS_A1(-24 * SIZE, A1, %xmm8)
	mulps	%xmm12,  %xmm4
	SUBPS	%xmm4,  %xmm1

	pshufd	$0xb1, %xmm9, %xmm5
	mulps	%xmm12,  %xmm9
	addps	%xmm9, %xmm2
	MOVUPS_A1(-24 * SIZE, A2, %xmm9)
	mulps	%xmm12,  %xmm5
	SUBPS	%xmm5, %xmm3

	movaps	-24 * SIZE(X1), %xmm12

	pshufd	$0xb1, %xmm10,  %xmm6
	mulps	%xmm13,  %xmm10
	addps	%xmm10,  %xmm0
	MOVUPS_A1(-20 * SIZE, A1, %xmm10)
	mulps	%xmm13,  %xmm6
	SUBPS	%xmm6,  %xmm1

	pshufd	$0xb1, %xmm11, %xmm7
	mulps	%xmm13,  %xmm11
	addps	%xmm11, %xmm2
	MOVUPS_A1(-20 * SIZE, A2, %xmm11)
	mulps	%xmm13,  %xmm7
	SUBPS	%xmm7, %xmm3

	movaps	-20 * SIZE(X1), %xmm13

	pshufd	$0xb1, %xmm8,  %xmm4
	mulps	%xmm12,  %xmm8
	addps	%xmm8,  %xmm0
	mulps	%xmm12,  %xmm4
	SUBPS	%xmm4,  %xmm1

	pshufd	$0xb1, %xmm9, %xmm5
	mulps	%xmm12,  %xmm9
	addps	%xmm9, %xmm2
	mulps	%xmm12,  %xmm5
	SUBPS	%xmm5, %xmm3

	movaps	-16 * SIZE(X1), %xmm12

	pshufd	$0xb1, %xmm10,  %xmm6
	mulps	%xmm13,  %xmm10
	addps	%xmm10,  %xmm0
	mulps	%xmm13,  %xmm6
	SUBPS	%xmm6,  %xmm1

	pshufd	$0xb1, %xmm11, %xmm7
	mulps	%xmm13,  %xmm11
	addps	%xmm11, %xmm2
	mulps	%xmm13,  %xmm7
	SUBPS	%xmm7, %xmm3

	movaps	-12 * SIZE(X1), %xmm13

	subq	 $-16 * SIZE, A1
	subq	 $-16 * SIZE, A2
	subq	 $-16 * SIZE, X1
	ALIGN_3

.L25:
	testq	$4, MM
	je	.L27

	MOVUPS_A1(-32 * SIZE, A1, %xmm8)
	MOVUPS_A1(-32 * SIZE, A2, %xmm9)
	MOVUPS_A1(-28 * SIZE, A1, %xmm10)
	MOVUPS_A1(-28 * SIZE, A2, %xmm11)

	pshufd	$0xb1, %xmm8,  %xmm4
	mulps	%xmm12,  %xmm8
	addps	%xmm8,  %xmm0
	mulps	%xmm12,  %xmm4
	SUBPS	%xmm4,  %xmm1

	pshufd	$0xb1, %xmm9, %xmm5
	mulps	%xmm12,  %xmm9
	addps	%xmm9, %xmm2
	mulps	%xmm12,  %xmm5
	SUBPS	%xmm5, %xmm3

	movaps	-24 * SIZE(X1), %xmm12

	pshufd	$0xb1, %xmm10,  %xmm6
	mulps	%xmm13,  %xmm10
	addps	%xmm10,  %xmm0
	mulps	%xmm13,  %xmm6
	SUBPS	%xmm6,  %xmm1

	pshufd	$0xb1, %xmm11, %xmm7
	mulps	%xmm13,  %xmm11
	addps	%xmm11, %xmm2
	mulps	%xmm13,  %xmm7
	SUBPS	%xmm7, %xmm3

	movaps	-20 * SIZE(X1), %xmm13

	addq	 $8 * SIZE, A1
	addq	 $8 * SIZE, A2
	ALIGN_3

.L27:
	testq	$2, MM
	je	.L28

	MOVUPS_A1(-32 * SIZE, A1, %xmm8)
	MOVUPS_A1(-32 * SIZE, A2, %xmm9)

	pshufd	$0xb1, %xmm8,  %xmm4
	mulps	%xmm12,  %xmm8
	addps	%xmm8,  %xmm0
	mulps	%xmm12,  %xmm4
	SUBPS	%xmm4,  %xmm1

	pshufd	$0xb1, %xmm9, %xmm5
	mulps	%xmm12,  %xmm9
	addps	%xmm9, %xmm2
	mulps	%xmm12,  %xmm5
	SUBPS	%xmm5, %xmm3

	movaps	%xmm13, %xmm12

	addq	 $4 * SIZE, A1
	addq	 $4 * SIZE, A2
	ALIGN_3

.L28:
	testq	$1, MM
	je	.L29

#ifdef movsd
	xorps	%xmm8, %xmm8
#endif
	movsd	-32 * SIZE(A1), %xmm8
#ifdef movsd
	xorps	%xmm9, %xmm9
#endif
	movsd	-32 * SIZE(A2), %xmm9

	pshufd	$0xb1, %xmm8,  %xmm4
	mulps	%xmm12,  %xmm8
	addps	%xmm8,  %xmm0
	mulps	%xmm12,  %xmm4
	SUBPS	%xmm4,  %xmm1

	pshufd	$0xb1, %xmm9, %xmm5
	mulps	%xmm12,  %xmm9
	addps	%xmm9, %xmm2
	mulps	%xmm12,  %xmm5
	SUBPS	%xmm5, %xmm3
	ALIGN_3

.L29:
	pcmpeqb	%xmm5, %xmm5
	psllq	$63,    %xmm5

#if (!defined(CONJ) && !defined(XCONJ)) || (defined(CONJ) && defined(XCONJ))
	xorps	%xmm5, %xmm0
	xorps	%xmm5, %xmm2
#else
	xorps	%xmm5, %xmm1
	xorps	%xmm5, %xmm3
#endif

#ifdef HAVE_SSE3
	haddps	%xmm1, %xmm0
	haddps	%xmm3, %xmm2
	haddps	%xmm2, %xmm0
#else
	movaps	%xmm0, %xmm8
	unpcklps %xmm1, %xmm0
	unpckhps %xmm1, %xmm8

	movaps	%xmm2, %xmm4
	unpcklps %xmm3, %xmm2
	unpckhps %xmm3, %xmm4

	addps	%xmm8, %xmm0
	addps	%xmm4, %xmm2

	movhlps	%xmm0, %xmm1
	movhlps	%xmm2, %xmm3

	addps	%xmm1, %xmm0
	addps	%xmm3, %xmm2
	movlhps %xmm2, %xmm0
#endif

	pshufd	$0xb1, %xmm0, %xmm1

#ifdef HAVE_SSE3
	movddup	ALPHA, %xmm15
#else
	movsd	ALPHA, %xmm15
	pshufd	$0x44, %xmm15, %xmm15
#endif

	mulps	%xmm15, %xmm0
	mulps	%xmm15, %xmm1

	xorps	%xmm5, %xmm0

#ifdef HAVE_SSE3
	haddps	%xmm1, %xmm0
#else
	movaps	%xmm0, %xmm2
	shufps	$0x88, %xmm1, %xmm0
	shufps	$0xdd, %xmm1, %xmm2
	
	addps	%xmm2, %xmm0
#endif

	movsd	(Y), %xmm12
	addq	INCY, Y
	movhps	(Y), %xmm12
	addq	INCY, Y

	shufps	$0xd8, %xmm0, %xmm0

	addps	%xmm12, %xmm0

	movlps	%xmm0, (Y1)
	addq	INCY, Y1
	movhps	%xmm0, (Y1)
	addq	INCY, Y1

#if GEMV_UNROLL == 2
	cmpq	$2, N
	jge	.L21
#endif
	ALIGN_3

.L30:
	cmpq	$1, N
	jl	.L999

	leaq	32 * SIZE(BUFFER), X1

	movq	A, A1

	xorps	%xmm0, %xmm0
	xorps	%xmm1, %xmm1

#ifdef ALIGNED_ACCESS
	cmpq	M, MM
	je	.L3X

#ifdef movsd
	xorps	%xmm8, %xmm8
#endif
	movsd	-32 * SIZE(A1), %xmm8
#ifdef movsd
	xorps	%xmm12, %xmm12
#endif
	movsd	-32 * SIZE(X1), %xmm12

	pshufd	$0xb1, %xmm8,  %xmm4
	mulps	%xmm12,  %xmm8
	addps	%xmm8,  %xmm0
	mulps	%xmm12,  %xmm4
	SUBPS	%xmm4,  %xmm1

	addq	 $2 * SIZE, A1
	addq	 $2 * SIZE, X1
	ALIGN_3
.L3X:
#endif

	movaps	-32 * SIZE(X1), %xmm12
	movaps	-28 * SIZE(X1), %xmm13

	movq	MM,  I
	sarq	$3,  I
	jle	.L35

	MOVUPS_A1(-32 * SIZE, A1, %xmm8)
	MOVUPS_A1(-28 * SIZE, A1, %xmm10)

	decq	 I
	jle	 .L34
	ALIGN_3

.L33:
#ifdef PREFETCH
	PREFETCH	(PREFETCHSIZE) * 4 - 128 + PREOFFSET(A1)
#endif

	pshufd	$0xb1, %xmm8,  %xmm4
	mulps	%xmm12,  %xmm8
	addps	%xmm8,  %xmm0
	MOVUPS_A1(-24 * SIZE, A1, %xmm8)
	mulps	%xmm12,  %xmm4
	SUBPS	%xmm4,  %xmm1

	movaps	-24 * SIZE(X1), %xmm12

	pshufd	$0xb1, %xmm10,  %xmm6
	mulps	%xmm13,  %xmm10
	addps	%xmm10,  %xmm0
	MOVUPS_A1(-20 * SIZE, A1, %xmm10)
	mulps	%xmm13,  %xmm6
	SUBPS	%xmm6,  %xmm1

	movaps	-20 * SIZE(X1), %xmm13

	pshufd	$0xb1, %xmm8,  %xmm4
	mulps	%xmm12,  %xmm8
	addps	%xmm8,  %xmm0
	MOVUPS_A1(-16 * SIZE, A1, %xmm8)
	mulps	%xmm12,  %xmm4
	SUBPS	%xmm4,  %xmm1

	movaps	-16 * SIZE(X1), %xmm12

#ifdef PREFETCHW
	PREFETCH	(PREFETCHSIZE) * 4 - 128 + PREOFFSET(X1)
#endif

	pshufd	$0xb1, %xmm10,  %xmm6
	mulps	%xmm13,  %xmm10
	addps	%xmm10,  %xmm0
	MOVUPS_A1(-12 * SIZE, A1, %xmm10)
	mulps	%xmm13,  %xmm6
	SUBPS	%xmm6,  %xmm1

	movaps	-12 * SIZE(X1), %xmm13

	subq	 $-16 * SIZE, A1
	subq	 $-16 * SIZE, X1

	subq	 $1, I
	BRANCH
	jg	.L33
	ALIGN_3

.L34:
	pshufd	$0xb1, %xmm8,  %xmm4
	mulps	%xmm12,  %xmm8
	addps	%xmm8,  %xmm0
	MOVUPS_A1(-24 * SIZE, A1, %xmm8)
	mulps	%xmm12,  %xmm4
	SUBPS	%xmm4,  %xmm1

	movaps	-24 * SIZE(X1), %xmm12

	pshufd	$0xb1, %xmm10,  %xmm6
	mulps	%xmm13,  %xmm10
	addps	%xmm10,  %xmm0
	MOVUPS_A1(-20 * SIZE, A1, %xmm10)
	mulps	%xmm13,  %xmm6
	SUBPS	%xmm6,  %xmm1

	movaps	-20 * SIZE(X1), %xmm13

	pshufd	$0xb1, %xmm8,  %xmm4
	mulps	%xmm12,  %xmm8
	addps	%xmm8,  %xmm0
	mulps	%xmm12,  %xmm4
	SUBPS	%xmm4,  %xmm1

	movaps	-16 * SIZE(X1), %xmm12

	pshufd	$0xb1, %xmm10,  %xmm6
	mulps	%xmm13,  %xmm10
	addps	%xmm10,  %xmm0
	mulps	%xmm13,  %xmm6
	SUBPS	%xmm6,  %xmm1

	movaps	-12 * SIZE(X1), %xmm13

	subq	 $-16 * SIZE, A1
	subq	 $-16 * SIZE, X1
	ALIGN_3

.L35:
	testq	$4, MM
	je	.L37

	MOVUPS_A1(-32 * SIZE, A1, %xmm8)
	MOVUPS_A1(-28 * SIZE, A1, %xmm10)

	pshufd	$0xb1, %xmm8,  %xmm4
	mulps	%xmm12,  %xmm8
	addps	%xmm8,  %xmm0
	mulps	%xmm12,  %xmm4
	SUBPS	%xmm4,  %xmm1

	movaps	-24 * SIZE(X1), %xmm12

	pshufd	$0xb1, %xmm10,  %xmm6
	mulps	%xmm13,  %xmm10
	addps	%xmm10,  %xmm0
	mulps	%xmm13,  %xmm6
	SUBPS	%xmm6,  %xmm1

	movaps	-20 * SIZE(X1), %xmm13

	addq	 $8 * SIZE, A1
	ALIGN_3

.L37:
	testq	$2, MM
	je	.L38

	MOVUPS_A1(-32 * SIZE, A1, %xmm8)

	pshufd	$0xb1, %xmm8,  %xmm4
	mulps	%xmm12,  %xmm8
	addps	%xmm8,  %xmm0
	mulps	%xmm12,  %xmm4
	SUBPS	%xmm4,  %xmm1

	movaps	%xmm13, %xmm12

	addq	 $4 * SIZE, A1
	ALIGN_3

.L38:
	testq	$1, MM
	je	.L39

#ifdef movsd
	xorps	%xmm8, %xmm8
#endif
	movsd	-32 * SIZE(A1), %xmm8

	pshufd	$0xb1, %xmm8,  %xmm4
	mulps	%xmm12,  %xmm8
	addps	%xmm8,  %xmm0
	mulps	%xmm12,  %xmm4
	SUBPS	%xmm4,  %xmm1
	ALIGN_3

.L39:
	pcmpeqb	%xmm5, %xmm5
	psllq	$63,    %xmm5

#if (!defined(CONJ) && !defined(XCONJ)) || (defined(CONJ) && defined(XCONJ))
	xorps	%xmm5, %xmm0
#else
	xorps	%xmm5, %xmm1
#endif

#ifdef HAVE_SSE3
	haddps	%xmm1, %xmm0
	haddps	%xmm0, %xmm0
#else
	movaps	%xmm0, %xmm8
	unpcklps %xmm1, %xmm0
	unpckhps %xmm1, %xmm8

	addps	%xmm8, %xmm0

	movhlps	%xmm0, %xmm1

	addps	%xmm1, %xmm0
#endif

	pshufd	$0xb1, %xmm0, %xmm1

#ifdef HAVE_SSE3
	movddup	ALPHA, %xmm15
#else
	movsd	ALPHA, %xmm15
	pshufd	$0x44, %xmm15, %xmm15
#endif

	mulps	%xmm15, %xmm0
	mulps	%xmm15, %xmm1

	xorps	%xmm5, %xmm0

#ifdef HAVE_SSE3
	haddps	%xmm1, %xmm0
#else
	movaps	%xmm0, %xmm2
	shufps	$0x88, %xmm1, %xmm0
	shufps	$0xdd, %xmm1, %xmm2
	
	addps	%xmm2, %xmm0
#endif

	movsd	(Y), %xmm12
	addq	INCY, Y

	shufps	$0xd8, %xmm0, %xmm0

	addps	%xmm12, %xmm0

	movlps	%xmm0, (Y1)
	addq	INCY, Y1
#ifdef ALIGNED_ACCESS
	jmp	.L999
	ALIGN_3

.L100:

#if GEMV_UNROLL >= 4
	cmpq	$4, N
	jl	.L110
	ALIGN_3

.L101:
	subq	$4, N

	leaq	32 * SIZE(BUFFER), X1

	movq	A, A1
	leaq	(A1, LDA, 2), A2
	leaq	(A1, LDA, 4), A

	xorps	%xmm0, %xmm0
	xorps	%xmm1, %xmm1
	xorps	%xmm2, %xmm2
	xorps	%xmm3, %xmm3
	xorps	%xmm4, %xmm4
	xorps	%xmm5, %xmm5
	xorps	%xmm6, %xmm6
	xorps	%xmm7, %xmm7

#ifdef ALIGNED_ACCESS
	cmpq	M, MM
	je	.L10X

#ifdef movsd
	xorps	%xmm8, %xmm8
#endif
	movsd	-32 * SIZE(A1), %xmm8
#ifdef movsd
	xorps	%xmm9, %xmm9
#endif
	movsd	-32 * SIZE(A1, LDA), %xmm9
#ifdef movsd
	xorps	%xmm10, %xmm10
#endif
	movsd	-32 * SIZE(A2), %xmm10
#ifdef movsd
	xorps	%xmm11, %xmm11
#endif
	movsd	-32 * SIZE(A2, LDA), %xmm11

#ifdef movsd
	xorps	%xmm12, %xmm12
#endif
	movsd	-32 * SIZE(X1), %xmm12

	pshufd	$0xb1,  %xmm8, %xmm14
	mulps	%xmm12, %xmm8
	addps	%xmm8,  %xmm0
	pshufd	$0xb1,  %xmm9, %xmm15
	mulps	%xmm12, %xmm9
	addps	%xmm9,  %xmm2

	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm1
	mulps	%xmm12, %xmm15
	SUBPS	%xmm15, %xmm3

	pshufd	$0xb1,  %xmm10, %xmm14
	mulps	%xmm12, %xmm10
	addps	%xmm10, %xmm4
	pshufd	$0xb1,  %xmm11, %xmm15
	mulps	%xmm12, %xmm11
	addps	%xmm11, %xmm6

	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm5
	mulps	%xmm12, %xmm15
	SUBPS	%xmm15, %xmm7

	addq	 $2 * SIZE, A1
	addq	 $2 * SIZE, A2
	addq	 $2 * SIZE, X1
	ALIGN_3
.L10X:
#endif

	movaps	-32 * SIZE(X1), %xmm12
	movaps	-28 * SIZE(X1), %xmm13

#ifdef PREFETCHW
	PREFETCHW	7 * SIZE(Y1)
#endif

	movq	MM,  I
	sarq	$3,  I
	jle	.L105

	movaps	-32 * SIZE(A1), %xmm8
	movsd	-32 * SIZE(A1, LDA), %xmm9
	movhps	-30 * SIZE(A1, LDA), %xmm9

	movaps	-32 * SIZE(A2), %xmm10
	movsd	-32 * SIZE(A2, LDA), %xmm11
	movhps	-30 * SIZE(A2, LDA), %xmm11

	decq	 I
	jle	 .L104
	ALIGN_3

.L103:
#ifdef PREFETCH
	PREFETCH	(PREFETCHSIZE) - 128 + PREOFFSET(A1)
#endif

	pshufd	$0xb1,  %xmm8, %xmm14
	mulps	%xmm12, %xmm8
	addps	%xmm8,  %xmm0
	movaps	-28 * SIZE(A1), %xmm8
	pshufd	$0xb1,  %xmm9, %xmm15
	mulps	%xmm12, %xmm9
	addps	%xmm9,  %xmm2
	movsd	-28 * SIZE(A1, LDA), %xmm9
	movhps	-26 * SIZE(A1, LDA), %xmm9

	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm1
	mulps	%xmm12, %xmm15
	SUBPS	%xmm15, %xmm3

	pshufd	$0xb1,  %xmm10, %xmm14
	mulps	%xmm12, %xmm10
	addps	%xmm10, %xmm4
	movaps	-28 * SIZE(A2), %xmm10
	pshufd	$0xb1,  %xmm11, %xmm15
	mulps	%xmm12, %xmm11
	addps	%xmm11, %xmm6
	movsd	-28 * SIZE(A2, LDA), %xmm11
	movhps	-26 * SIZE(A2, LDA), %xmm11

	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm5
	mulps	%xmm12, %xmm15
	movaps	-24 * SIZE(X1), %xmm12
	SUBPS	%xmm15, %xmm7

#ifdef PREFETCH
	PREFETCH	(PREFETCHSIZE) - 128 + PREOFFSET(A1, LDA)
#endif

	pshufd	$0xb1,  %xmm8, %xmm14
	mulps	%xmm13, %xmm8
	addps	%xmm8,  %xmm0
	movaps	-24 * SIZE(A1), %xmm8
	pshufd	$0xb1,  %xmm9, %xmm15
	mulps	%xmm13, %xmm9
	addps	%xmm9,  %xmm2
	movsd	-24 * SIZE(A1, LDA), %xmm9
	movhps	-22 * SIZE(A1, LDA), %xmm9

	mulps	%xmm13, %xmm14
	SUBPS	%xmm14, %xmm1
	mulps	%xmm13, %xmm15
	SUBPS	%xmm15, %xmm3

	pshufd	$0xb1,  %xmm10, %xmm14
	mulps	%xmm13, %xmm10
	addps	%xmm10, %xmm4
	movaps	-24 * SIZE(A2), %xmm10
	pshufd	$0xb1,  %xmm11, %xmm15
	mulps	%xmm13, %xmm11
	addps	%xmm11, %xmm6
	movsd	-24 * SIZE(A2, LDA), %xmm11
	movhps	-22 * SIZE(A2, LDA), %xmm11

	mulps	%xmm13, %xmm14
	SUBPS	%xmm14, %xmm5
	mulps	%xmm13, %xmm15
	movaps	-20 * SIZE(X1), %xmm13
	SUBPS	%xmm15, %xmm7

#ifdef PREFETCH
	PREFETCH	(PREFETCHSIZE) - 128 + PREOFFSET(A2)
#endif

	pshufd	$0xb1,  %xmm8, %xmm14
	mulps	%xmm12, %xmm8
	addps	%xmm8,  %xmm0
	movaps	-20 * SIZE(A1), %xmm8
	pshufd	$0xb1,  %xmm9, %xmm15
	mulps	%xmm12, %xmm9
	addps	%xmm9,  %xmm2
	movsd	-20 * SIZE(A1, LDA), %xmm9
	movhps	-18 * SIZE(A1, LDA), %xmm9

	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm1
	mulps	%xmm12, %xmm15
	SUBPS	%xmm15, %xmm3

	pshufd	$0xb1,  %xmm10, %xmm14
	mulps	%xmm12, %xmm10
	addps	%xmm10, %xmm4
	movaps	-20 * SIZE(A2), %xmm10
	pshufd	$0xb1,  %xmm11, %xmm15
	mulps	%xmm12, %xmm11
	addps	%xmm11, %xmm6
	movsd	-20 * SIZE(A2, LDA), %xmm11
	movhps	-18 * SIZE(A2, LDA), %xmm11

	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm5
	mulps	%xmm12, %xmm15
	movaps	-16 * SIZE(X1), %xmm12
	SUBPS	%xmm15, %xmm7

#ifdef PREFETCH
	PREFETCH	(PREFETCHSIZE) - 128 + PREOFFSET(A2, LDA)
#endif

	pshufd	$0xb1,  %xmm8, %xmm14
	mulps	%xmm13, %xmm8
	addps	%xmm8,  %xmm0
	movaps	-16 * SIZE(A1), %xmm8
	pshufd	$0xb1,  %xmm9, %xmm15
	mulps	%xmm13, %xmm9
	addps	%xmm9,  %xmm2
	movsd	-16 * SIZE(A1, LDA), %xmm9
	movhps	-14 * SIZE(A1, LDA), %xmm9

	mulps	%xmm13, %xmm14
	SUBPS	%xmm14, %xmm1
	mulps	%xmm13, %xmm15
	SUBPS	%xmm15, %xmm3

#ifdef PREFETCHW
	PREFETCH	(PREFETCHSIZE) - 128 + PREOFFSET(X1)
#endif

	pshufd	$0xb1,  %xmm10, %xmm14
	mulps	%xmm13, %xmm10
	addps	%xmm10, %xmm4
	movaps	-16 * SIZE(A2), %xmm10
	pshufd	$0xb1,  %xmm11, %xmm15
	mulps	%xmm13, %xmm11
	addps	%xmm11, %xmm6
	movsd	-16 * SIZE(A2, LDA), %xmm11
	movhps	-14 * SIZE(A2, LDA), %xmm11

	mulps	%xmm13, %xmm14
	SUBPS	%xmm14, %xmm5
	mulps	%xmm13, %xmm15
	movaps	-12 * SIZE(X1), %xmm13
	SUBPS	%xmm15, %xmm7

	subq	 $-16 * SIZE, A1
	subq	 $-16 * SIZE, A2
	subq	 $-16 * SIZE, X1

	subq	 $1, I
	BRANCH
	jg	.L103
	ALIGN_3

.L104:
	pshufd	$0xb1,  %xmm8, %xmm14
	mulps	%xmm12, %xmm8
	addps	%xmm8,  %xmm0
	movaps	-28 * SIZE(A1), %xmm8
	pshufd	$0xb1,  %xmm9, %xmm15
	mulps	%xmm12, %xmm9
	addps	%xmm9,  %xmm2
	movsd	-28 * SIZE(A1, LDA), %xmm9
	movhps	-26 * SIZE(A1, LDA), %xmm9

	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm1
	mulps	%xmm12, %xmm15
	SUBPS	%xmm15, %xmm3

	pshufd	$0xb1,  %xmm10, %xmm14
	mulps	%xmm12, %xmm10
	addps	%xmm10, %xmm4
	movaps	-28 * SIZE(A2), %xmm10
	pshufd	$0xb1,  %xmm11, %xmm15
	mulps	%xmm12, %xmm11
	addps	%xmm11, %xmm6
	movsd	-28 * SIZE(A2, LDA), %xmm11
	movhps	-26 * SIZE(A2, LDA), %xmm11

	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm5
	mulps	%xmm12, %xmm15
	movaps	-24 * SIZE(X1), %xmm12
	SUBPS	%xmm15, %xmm7

	pshufd	$0xb1,  %xmm8, %xmm14
	mulps	%xmm13, %xmm8
	addps	%xmm8,  %xmm0
	movaps	-24 * SIZE(A1), %xmm8
	pshufd	$0xb1,  %xmm9, %xmm15
	mulps	%xmm13, %xmm9
	addps	%xmm9,  %xmm2
	movsd	-24 * SIZE(A1, LDA), %xmm9
	movhps	-22 * SIZE(A1, LDA), %xmm9

	mulps	%xmm13, %xmm14
	SUBPS	%xmm14, %xmm1
	mulps	%xmm13, %xmm15
	SUBPS	%xmm15, %xmm3

	pshufd	$0xb1,  %xmm10, %xmm14
	mulps	%xmm13, %xmm10
	addps	%xmm10, %xmm4
	movaps	-24 * SIZE(A2), %xmm10
	pshufd	$0xb1,  %xmm11, %xmm15
	mulps	%xmm13, %xmm11
	addps	%xmm11, %xmm6
	movsd	-24 * SIZE(A2, LDA), %xmm11
	movhps	-22 * SIZE(A2, LDA), %xmm11

	mulps	%xmm13, %xmm14
	SUBPS	%xmm14, %xmm5
	mulps	%xmm13, %xmm15
	movaps	-20 * SIZE(X1), %xmm13
	SUBPS	%xmm15, %xmm7

	pshufd	$0xb1,  %xmm8, %xmm14
	mulps	%xmm12, %xmm8
	addps	%xmm8,  %xmm0
	movaps	-20 * SIZE(A1), %xmm8
	pshufd	$0xb1,  %xmm9, %xmm15
	mulps	%xmm12, %xmm9
	addps	%xmm9,  %xmm2
	movsd	-20 * SIZE(A1, LDA), %xmm9
	movhps	-18 * SIZE(A1, LDA), %xmm9

	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm1
	mulps	%xmm12, %xmm15
	SUBPS	%xmm15, %xmm3

	pshufd	$0xb1,  %xmm10, %xmm14
	mulps	%xmm12, %xmm10
	addps	%xmm10, %xmm4
	movaps	-20 * SIZE(A2), %xmm10
	pshufd	$0xb1,  %xmm11, %xmm15
	mulps	%xmm12, %xmm11
	addps	%xmm11, %xmm6
	movsd	-20 * SIZE(A2, LDA), %xmm11
	movhps	-18 * SIZE(A2, LDA), %xmm11

	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm5
	mulps	%xmm12, %xmm15
	movaps	-16 * SIZE(X1), %xmm12
	SUBPS	%xmm15, %xmm7

	pshufd	$0xb1,  %xmm8, %xmm14
	mulps	%xmm13, %xmm8
	addps	%xmm8,  %xmm0
	pshufd	$0xb1,  %xmm9, %xmm15
	mulps	%xmm13, %xmm9
	addps	%xmm9,  %xmm2

	mulps	%xmm13, %xmm14
	SUBPS	%xmm14, %xmm1
	mulps	%xmm13, %xmm15
	SUBPS	%xmm15, %xmm3

	pshufd	$0xb1,  %xmm10, %xmm14
	mulps	%xmm13, %xmm10
	addps	%xmm10, %xmm4
	pshufd	$0xb1,  %xmm11, %xmm15
	mulps	%xmm13, %xmm11
	addps	%xmm11, %xmm6

	mulps	%xmm13, %xmm14
	SUBPS	%xmm14, %xmm5
	mulps	%xmm13, %xmm15
	movaps	-12 * SIZE(X1), %xmm13
	SUBPS	%xmm15, %xmm7

	subq	 $-16 * SIZE, A1
	subq	 $-16 * SIZE, A2
	subq	 $-16 * SIZE, X1
	ALIGN_3

.L105:
	testq	$4, MM
	je	.L107

	movaps	-32 * SIZE(A1), %xmm8
	movsd	-32 * SIZE(A1, LDA), %xmm9
	movhps	-30 * SIZE(A1, LDA), %xmm9

	movaps	-32 * SIZE(A2), %xmm10
	movsd	-32 * SIZE(A2, LDA), %xmm11
	movhps	-30 * SIZE(A2, LDA), %xmm11

	pshufd	$0xb1,  %xmm8, %xmm14
	mulps	%xmm12, %xmm8
	addps	%xmm8,  %xmm0
	movaps	-28 * SIZE(A1), %xmm8
	pshufd	$0xb1,  %xmm9, %xmm15
	mulps	%xmm12, %xmm9
	addps	%xmm9,  %xmm2
	movsd	-28 * SIZE(A1, LDA), %xmm9
	movhps	-26 * SIZE(A1, LDA), %xmm9

	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm1
	mulps	%xmm12, %xmm15
	SUBPS	%xmm15, %xmm3

	pshufd	$0xb1,  %xmm10, %xmm14
	mulps	%xmm12, %xmm10
	addps	%xmm10, %xmm4
	movaps	-28 * SIZE(A2), %xmm10
	pshufd	$0xb1,  %xmm11, %xmm15
	mulps	%xmm12, %xmm11
	addps	%xmm11, %xmm6
	movsd	-28 * SIZE(A2, LDA), %xmm11
	movhps	-26 * SIZE(A2, LDA), %xmm11

	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm5
	mulps	%xmm12, %xmm15
	movaps	-24 * SIZE(X1), %xmm12
	SUBPS	%xmm15, %xmm7

	pshufd	$0xb1,  %xmm8, %xmm14
	mulps	%xmm13, %xmm8
	addps	%xmm8,  %xmm0
	pshufd	$0xb1,  %xmm9, %xmm15
	mulps	%xmm13, %xmm9
	addps	%xmm9,  %xmm2

	mulps	%xmm13, %xmm14
	SUBPS	%xmm14, %xmm1
	mulps	%xmm13, %xmm15
	SUBPS	%xmm15, %xmm3

	pshufd	$0xb1,  %xmm10, %xmm14
	mulps	%xmm13, %xmm10
	addps	%xmm10, %xmm4
	pshufd	$0xb1,  %xmm11, %xmm15
	mulps	%xmm13, %xmm11
	addps	%xmm11, %xmm6

	mulps	%xmm13, %xmm14
	SUBPS	%xmm14, %xmm5
	mulps	%xmm13, %xmm15
	movaps	-20 * SIZE(X1), %xmm13
	SUBPS	%xmm15, %xmm7

	addq	 $8 * SIZE, A1
	addq	 $8 * SIZE, A2
	addq	 $8 * SIZE, X1
	ALIGN_3

.L107:
	testq	$2, MM
	je	.L108

	movaps	-32 * SIZE(A1), %xmm8
	movsd	-32 * SIZE(A1, LDA), %xmm9
	movhps	-30 * SIZE(A1, LDA), %xmm9

	movaps	-32 * SIZE(A2), %xmm10
	movsd	-32 * SIZE(A2, LDA), %xmm11
	movhps	-30 * SIZE(A2, LDA), %xmm11

	pshufd	$0xb1,  %xmm8, %xmm14
	mulps	%xmm12, %xmm8
	addps	%xmm8,  %xmm0
	pshufd	$0xb1,  %xmm9, %xmm15
	mulps	%xmm12, %xmm9
	addps	%xmm9,  %xmm2

	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm1
	mulps	%xmm12, %xmm15
	SUBPS	%xmm15, %xmm3

	pshufd	$0xb1,  %xmm10, %xmm14
	mulps	%xmm12, %xmm10
	addps	%xmm10, %xmm4
	pshufd	$0xb1,  %xmm11, %xmm15
	mulps	%xmm12, %xmm11
	addps	%xmm11, %xmm6

	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm5
	mulps	%xmm12, %xmm15
	SUBPS	%xmm15, %xmm7

	movaps	%xmm13, %xmm12

	addq	 $4 * SIZE, A1
	addq	 $4 * SIZE, A2
	ALIGN_3

.L108:
	testq	$1, MM
	je	.L109

#ifdef movsd
	xorps	%xmm8, %xmm8
#endif
	movsd	-32 * SIZE(A1), %xmm8
#ifdef movsd
	xorps	%xmm9, %xmm9
#endif
	movsd	-32 * SIZE(A1, LDA), %xmm9
#ifdef movsd
	xorps	%xmm10, %xmm10
#endif
	movsd	-32 * SIZE(A2), %xmm10
#ifdef movsd
	xorps	%xmm11, %xmm11
#endif
	movsd	-32 * SIZE(A2, LDA), %xmm11

	pshufd	$0xb1,  %xmm8, %xmm14
	mulps	%xmm12, %xmm8
	addps	%xmm8,  %xmm0
	pshufd	$0xb1,  %xmm9, %xmm15
	mulps	%xmm12, %xmm9
	addps	%xmm9,  %xmm2

	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm1
	mulps	%xmm12, %xmm15
	SUBPS	%xmm15, %xmm3

	pshufd	$0xb1,  %xmm10, %xmm14
	mulps	%xmm12, %xmm10
	addps	%xmm10, %xmm4
	pshufd	$0xb1,  %xmm11, %xmm15
	mulps	%xmm12, %xmm11
	addps	%xmm11, %xmm6

	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm5
	mulps	%xmm12, %xmm15
	SUBPS	%xmm15, %xmm7
	ALIGN_3

.L109:
	pcmpeqb	%xmm11, %xmm11
	psllq	$63,    %xmm11

#if (!defined(CONJ) && !defined(XCONJ)) || (defined(CONJ) && defined(XCONJ))
	xorps	%xmm11, %xmm0
	xorps	%xmm11, %xmm2
	xorps	%xmm11, %xmm4
	xorps	%xmm11, %xmm6
#else
	xorps	%xmm11, %xmm1
	xorps	%xmm11, %xmm3
	xorps	%xmm11, %xmm5
	xorps	%xmm11, %xmm7
#endif

#ifdef HAVE_SSE3
	haddps	%xmm1, %xmm0
	haddps	%xmm3, %xmm2
	haddps	%xmm2, %xmm0

	haddps	%xmm5, %xmm4
	haddps	%xmm7, %xmm6
	haddps	%xmm6, %xmm4
#else

	movaps	%xmm0, %xmm8
	unpcklps %xmm1, %xmm0
	unpckhps %xmm1, %xmm8

	movaps	%xmm2, %xmm9
	unpcklps %xmm3, %xmm2
	unpckhps %xmm3, %xmm9

	movaps	%xmm4, %xmm10
	unpcklps %xmm5, %xmm4
	unpckhps %xmm5, %xmm10

	movaps	%xmm6, %xmm11
	unpcklps %xmm7, %xmm6
	unpckhps %xmm7, %xmm11

	addps	%xmm8,  %xmm0
	addps	%xmm9,  %xmm2
	addps	%xmm10, %xmm4
	addps	%xmm11, %xmm6

	movhlps	%xmm0, %xmm1
	movhlps	%xmm2, %xmm3
	movhlps	%xmm4, %xmm5
	movhlps	%xmm6, %xmm7

	addps	%xmm1, %xmm0
	addps	%xmm3, %xmm2
	addps	%xmm5, %xmm4
	addps	%xmm7, %xmm6

	movlhps %xmm2, %xmm0
	movlhps %xmm6, %xmm4
#endif

	pshufd	$0xb1, %xmm0, %xmm1
	pshufd	$0xb1, %xmm4, %xmm5

#ifdef HAVE_SSE3
	movddup	ALPHA, %xmm15
#else
	movsd	ALPHA, %xmm15
	pshufd	$0x44, %xmm15, %xmm15
#endif

	mulps	%xmm15, %xmm0
	mulps	%xmm15, %xmm1
	mulps	%xmm15, %xmm4
	mulps	%xmm15, %xmm5

	xorps	%xmm11, %xmm0
	xorps	%xmm11, %xmm4

#ifdef HAVE_SSE3
	haddps	%xmm1, %xmm0
	haddps	%xmm5, %xmm4
#else
	movaps	%xmm0, %xmm2
	shufps	$0x88, %xmm1, %xmm0
	shufps	$0xdd, %xmm1, %xmm2

	movaps	%xmm4, %xmm6
	shufps	$0x88, %xmm5, %xmm4
	shufps	$0xdd, %xmm5, %xmm6
	
	addps	%xmm2, %xmm0
	addps	%xmm6, %xmm4
#endif

	movsd	(Y), %xmm2
	addq	INCY, Y
	movhps	(Y), %xmm2
	addq	INCY, Y
	movsd	(Y), %xmm6
	addq	INCY, Y
	movhps	(Y), %xmm6
	addq	INCY, Y

	shufps	$0xd8, %xmm0, %xmm0
	shufps	$0xd8, %xmm4, %xmm4

	addps	%xmm2, %xmm0
	addps	%xmm6, %xmm4

	movlps	%xmm0, (Y1)
	addq	INCY, Y1
	movhps	%xmm0, (Y1)
	addq	INCY, Y1
	movlps	%xmm4, (Y1)
	addq	INCY, Y1
	movhps	%xmm4, (Y1)
	addq	INCY, Y1

	cmpq	$4, N
	jge	.L101
	ALIGN_3

.L110:
#endif

	cmpq	$2, N
	jl	.L120
#if GEMV_UNROLL == 2
	ALIGN_3

.L111:
#endif
	subq	$2, N

	leaq	32 * SIZE(BUFFER), X1

	movq	A, A1
	leaq	(A1, LDA), A2
	leaq	(A1, LDA, 2), A

	xorps	%xmm0, %xmm0
	xorps	%xmm1, %xmm1
	xorps	%xmm2, %xmm2
	xorps	%xmm3, %xmm3

#ifdef ALIGNED_ACCESS
	cmpq	M, MM
	je	.L11X

#ifdef movsd
	xorps	%xmm8, %xmm8
#endif
	movsd	-32 * SIZE(A1), %xmm8
#ifdef movsd
	xorps	%xmm9, %xmm9
#endif
	movsd	-32 * SIZE(A2), %xmm9

#ifdef movsd
	xorps	%xmm12, %xmm12
#endif
	movsd	-32 * SIZE(X1), %xmm12

	pshufd	$0xb1, %xmm8,  %xmm4
	mulps	%xmm12,  %xmm8
	addps	%xmm8,  %xmm0
	mulps	%xmm12,  %xmm4
	SUBPS	%xmm4,  %xmm1

	pshufd	$0xb1, %xmm9, %xmm5
	mulps	%xmm12,  %xmm9
	addps	%xmm9, %xmm2
	mulps	%xmm12,  %xmm5
	SUBPS	%xmm5, %xmm3

	addq	 $2 * SIZE, A1
	addq	 $2 * SIZE, A2
	addq	 $2 * SIZE, X1
	ALIGN_3
.L11X:
#endif

	movaps	-32 * SIZE(X1), %xmm12
	movaps	-28 * SIZE(X1), %xmm13

#if (GEMV_UNROLL == 2) && defined(PREFETCHW)
	PREFETCHW	3 * SIZE(Y1)
#endif

	movq	MM,  I
	sarq	$3,  I
	jle	.L115

	movaps	-32 * SIZE(A1), %xmm8
	movsd	-32 * SIZE(A2), %xmm9
	movhps	-30 * SIZE(A2), %xmm9

	movaps	-28 * SIZE(A1), %xmm10
	movsd	-28 * SIZE(A2), %xmm11
	movhps	-26 * SIZE(A2), %xmm11

	decq	 I
	jle	 .L114
	ALIGN_3

.L113:
#ifdef PREFETCH
	PREFETCH	(PREFETCHSIZE) * 2 - 128 + PREOFFSET(A1)
#endif

	pshufd	$0xb1, %xmm8,  %xmm4
	mulps	%xmm12,  %xmm8
	addps	%xmm8,  %xmm0
	movaps	-24 * SIZE(A1), %xmm8
	mulps	%xmm12,  %xmm4
	SUBPS	%xmm4,  %xmm1

	pshufd	$0xb1, %xmm9, %xmm5
	mulps	%xmm12,  %xmm9
	addps	%xmm9, %xmm2
	movsd	-24 * SIZE(A2), %xmm9
	movhps	-22 * SIZE(A2), %xmm9
	mulps	%xmm12,  %xmm5
	SUBPS	%xmm5, %xmm3

	movaps	-24 * SIZE(X1), %xmm12

	pshufd	$0xb1, %xmm10,  %xmm6
	mulps	%xmm13,  %xmm10
	addps	%xmm10,  %xmm0
	movaps	-20 * SIZE(A1), %xmm10
	mulps	%xmm13,  %xmm6
	SUBPS	%xmm6,  %xmm1

	pshufd	$0xb1, %xmm11, %xmm7
	mulps	%xmm13,  %xmm11
	addps	%xmm11, %xmm2
	movsd	-20 * SIZE(A2), %xmm11
	movhps	-18 * SIZE(A2), %xmm11
	mulps	%xmm13,  %xmm7
	SUBPS	%xmm7, %xmm3

	movaps	-20 * SIZE(X1), %xmm13

#ifdef PREFETCH
	PREFETCH	(PREFETCHSIZE) * 2 - 128 + PREOFFSET(A2)
#endif

	pshufd	$0xb1, %xmm8,  %xmm4
	mulps	%xmm12,  %xmm8
	addps	%xmm8,  %xmm0
	movaps	-16 * SIZE(A1), %xmm8
	mulps	%xmm12,  %xmm4
	SUBPS	%xmm4,  %xmm1

	pshufd	$0xb1, %xmm9, %xmm5
	mulps	%xmm12,  %xmm9
	addps	%xmm9, %xmm2
	movsd	-16 * SIZE(A2), %xmm9
	movhps	-14 * SIZE(A2), %xmm9
	mulps	%xmm12,  %xmm5
	SUBPS	%xmm5, %xmm3

	movaps	-16 * SIZE(X1), %xmm12

#ifdef PREFETCHW
	PREFETCH	(PREFETCHSIZE) * 2 - 128 + PREOFFSET(X1)
#endif

	pshufd	$0xb1, %xmm10,  %xmm6
	mulps	%xmm13,  %xmm10
	addps	%xmm10,  %xmm0
	movaps	-12 * SIZE(A1), %xmm10
	mulps	%xmm13,  %xmm6
	SUBPS	%xmm6,  %xmm1

	pshufd	$0xb1, %xmm11, %xmm7
	mulps	%xmm13,  %xmm11
	addps	%xmm11, %xmm2
	movsd	-12 * SIZE(A2), %xmm11
	movhps	-10 * SIZE(A2), %xmm11
	mulps	%xmm13,  %xmm7
	SUBPS	%xmm7, %xmm3

	movaps	-12 * SIZE(X1), %xmm13

	subq	 $-16 * SIZE, A1
	subq	 $-16 * SIZE, A2
	subq	 $-16 * SIZE, X1

	subq	 $1, I
	BRANCH
	jg	.L113
	ALIGN_3

.L114:
	pshufd	$0xb1, %xmm8,  %xmm4
	mulps	%xmm12,  %xmm8
	addps	%xmm8,  %xmm0
	movaps	-24 * SIZE(A1), %xmm8
	mulps	%xmm12,  %xmm4
	SUBPS	%xmm4,  %xmm1

	pshufd	$0xb1, %xmm9, %xmm5
	mulps	%xmm12,  %xmm9
	addps	%xmm9, %xmm2
	movsd	-24 * SIZE(A2), %xmm9
	movhps	-22 * SIZE(A2), %xmm9
	mulps	%xmm12,  %xmm5
	SUBPS	%xmm5, %xmm3

	movaps	-24 * SIZE(X1), %xmm12

	pshufd	$0xb1, %xmm10,  %xmm6
	mulps	%xmm13,  %xmm10
	addps	%xmm10,  %xmm0
	movaps	-20 * SIZE(A1), %xmm10
	mulps	%xmm13,  %xmm6
	SUBPS	%xmm6,  %xmm1

	pshufd	$0xb1, %xmm11, %xmm7
	mulps	%xmm13,  %xmm11
	addps	%xmm11, %xmm2
	movsd	-20 * SIZE(A2), %xmm11
	movhps	-18 * SIZE(A2), %xmm11
	mulps	%xmm13,  %xmm7
	SUBPS	%xmm7, %xmm3

	movaps	-20 * SIZE(X1), %xmm13

	pshufd	$0xb1, %xmm8,  %xmm4
	mulps	%xmm12,  %xmm8
	addps	%xmm8,  %xmm0
	mulps	%xmm12,  %xmm4
	SUBPS	%xmm4,  %xmm1

	pshufd	$0xb1, %xmm9, %xmm5
	mulps	%xmm12,  %xmm9
	addps	%xmm9, %xmm2
	mulps	%xmm12,  %xmm5
	SUBPS	%xmm5, %xmm3

	movaps	-16 * SIZE(X1), %xmm12

	pshufd	$0xb1, %xmm10,  %xmm6
	mulps	%xmm13,  %xmm10
	addps	%xmm10,  %xmm0
	mulps	%xmm13,  %xmm6
	SUBPS	%xmm6,  %xmm1

	pshufd	$0xb1, %xmm11, %xmm7
	mulps	%xmm13,  %xmm11
	addps	%xmm11, %xmm2
	mulps	%xmm13,  %xmm7
	SUBPS	%xmm7, %xmm3

	movaps	-12 * SIZE(X1), %xmm13

	subq	 $-16 * SIZE, A1
	subq	 $-16 * SIZE, A2
	subq	 $-16 * SIZE, X1
	ALIGN_3

.L115:
	testq	$4, MM
	je	.L117

	movaps	-32 * SIZE(A1), %xmm8
	movsd	-32 * SIZE(A2), %xmm9
	movhps	-30 * SIZE(A2), %xmm9

	movaps	-28 * SIZE(A1), %xmm10
	movsd	-28 * SIZE(A2), %xmm11
	movhps	-26 * SIZE(A2), %xmm11

	pshufd	$0xb1, %xmm8,  %xmm4
	mulps	%xmm12,  %xmm8
	addps	%xmm8,  %xmm0
	mulps	%xmm12,  %xmm4
	SUBPS	%xmm4,  %xmm1

	pshufd	$0xb1, %xmm9, %xmm5
	mulps	%xmm12,  %xmm9
	addps	%xmm9, %xmm2
	mulps	%xmm12,  %xmm5
	SUBPS	%xmm5, %xmm3

	movaps	-24 * SIZE(X1), %xmm12

	pshufd	$0xb1, %xmm10,  %xmm6
	mulps	%xmm13,  %xmm10
	addps	%xmm10,  %xmm0
	mulps	%xmm13,  %xmm6
	SUBPS	%xmm6,  %xmm1

	pshufd	$0xb1, %xmm11, %xmm7
	mulps	%xmm13,  %xmm11
	addps	%xmm11, %xmm2
	mulps	%xmm13,  %xmm7
	SUBPS	%xmm7, %xmm3

	movaps	-20 * SIZE(X1), %xmm13

	addq	 $8 * SIZE, A1
	addq	 $8 * SIZE, A2
	ALIGN_3

.L117:
	testq	$2, MM
	je	.L118

	movaps	-32 * SIZE(A1), %xmm8
	movsd	-32 * SIZE(A2), %xmm9
	movhps	-30 * SIZE(A2), %xmm9

	pshufd	$0xb1, %xmm8,  %xmm4
	mulps	%xmm12,  %xmm8
	addps	%xmm8,  %xmm0
	mulps	%xmm12,  %xmm4
	SUBPS	%xmm4,  %xmm1

	pshufd	$0xb1, %xmm9, %xmm5
	mulps	%xmm12,  %xmm9
	addps	%xmm9, %xmm2
	mulps	%xmm12,  %xmm5
	SUBPS	%xmm5, %xmm3

	movaps	%xmm13, %xmm12

	addq	 $4 * SIZE, A1
	addq	 $4 * SIZE, A2
	ALIGN_3

.L118:
	testq	$1, MM
	je	.L119

#ifdef movsd
	xorps	%xmm8, %xmm8
#endif
	movsd	-32 * SIZE(A1), %xmm8
#ifdef movsd
	xorps	%xmm9, %xmm9
#endif
	movsd	-32 * SIZE(A2), %xmm9

	pshufd	$0xb1, %xmm8,  %xmm4
	mulps	%xmm12,  %xmm8
	addps	%xmm8,  %xmm0
	mulps	%xmm12,  %xmm4
	SUBPS	%xmm4,  %xmm1

	pshufd	$0xb1, %xmm9, %xmm5
	mulps	%xmm12,  %xmm9
	addps	%xmm9, %xmm2
	mulps	%xmm12,  %xmm5
	SUBPS	%xmm5, %xmm3
	ALIGN_3

.L119:
	pcmpeqb	%xmm5, %xmm5
	psllq	$63,    %xmm5

#if (!defined(CONJ) && !defined(XCONJ)) || (defined(CONJ) && defined(XCONJ))
	xorps	%xmm5, %xmm0
	xorps	%xmm5, %xmm2
#else
	xorps	%xmm5, %xmm1
	xorps	%xmm5, %xmm3
#endif

#ifdef HAVE_SSE3
	haddps	%xmm1, %xmm0
	haddps	%xmm3, %xmm2
	haddps	%xmm2, %xmm0
#else
	movaps	%xmm0, %xmm8
	unpcklps %xmm1, %xmm0
	unpckhps %xmm1, %xmm8

	movaps	%xmm2, %xmm4
	unpcklps %xmm3, %xmm2
	unpckhps %xmm3, %xmm4

	addps	%xmm8, %xmm0
	addps	%xmm4, %xmm2

	movhlps	%xmm0, %xmm1
	movhlps	%xmm2, %xmm3

	addps	%xmm1, %xmm0
	addps	%xmm3, %xmm2
	movlhps %xmm2, %xmm0
#endif

	pshufd	$0xb1, %xmm0, %xmm1

#ifdef HAVE_SSE3
	movddup	ALPHA, %xmm15
#else
	movsd	ALPHA, %xmm15
	pshufd	$0x44, %xmm15, %xmm15
#endif

	mulps	%xmm15, %xmm0
	mulps	%xmm15, %xmm1

	xorps	%xmm5, %xmm0

#ifdef HAVE_SSE3
	haddps	%xmm1, %xmm0
#else
	movaps	%xmm0, %xmm2
	shufps	$0x88, %xmm1, %xmm0
	shufps	$0xdd, %xmm1, %xmm2
	
	addps	%xmm2, %xmm0
#endif

	movsd	(Y), %xmm12
	addq	INCY, Y
	movhps	(Y), %xmm12
	addq	INCY, Y

	shufps	$0xd8, %xmm0, %xmm0

	addps	%xmm12, %xmm0

	movlps	%xmm0, (Y1)
	addq	INCY, Y1
	movhps	%xmm0, (Y1)
	addq	INCY, Y1

#if GEMV_UNROLL == 2
	cmpq	$2, N
	jge	.L111
#endif
	ALIGN_3

.L120:
	cmpq	$1, N
	jl	.L999

	leaq	32 * SIZE(BUFFER), X1

	movq	A, A1

	xorps	%xmm0, %xmm0
	xorps	%xmm1, %xmm1

#ifdef ALIGNED_ACCESS
	cmpq	M, MM
	je	.L12X

#ifdef movsd
	xorps	%xmm8, %xmm8
#endif
	movsd	-32 * SIZE(A1), %xmm8
#ifdef movsd
	xorps	%xmm12, %xmm12
#endif
	movsd	-32 * SIZE(X1), %xmm12

	pshufd	$0xb1, %xmm8,  %xmm4
	mulps	%xmm12,  %xmm8
	addps	%xmm8,  %xmm0
	mulps	%xmm12,  %xmm4
	SUBPS	%xmm4,  %xmm1

	addq	 $2 * SIZE, A1
	addq	 $2 * SIZE, X1
	ALIGN_3
.L12X:
#endif

	movaps	-32 * SIZE(X1), %xmm12
	movaps	-28 * SIZE(X1), %xmm13

	movq	MM,  I
	sarq	$3,  I
	jle	.L125

	MOVUPS_A1(-32 * SIZE, A1, %xmm8)
	MOVUPS_A1(-28 * SIZE, A1, %xmm10)

	decq	 I
	jle	 .L124
	ALIGN_3

.L123:
#ifdef PREFETCH
	PREFETCH	(PREFETCHSIZE) * 4 - 128 + PREOFFSET(A1)
#endif

	pshufd	$0xb1, %xmm8,  %xmm4
	mulps	%xmm12,  %xmm8
	addps	%xmm8,  %xmm0
	MOVUPS_A1(-24 * SIZE, A1, %xmm8)
	mulps	%xmm12,  %xmm4
	SUBPS	%xmm4,  %xmm1

	movaps	-24 * SIZE(X1), %xmm12

	pshufd	$0xb1, %xmm10,  %xmm6
	mulps	%xmm13,  %xmm10
	addps	%xmm10,  %xmm0
	MOVUPS_A1(-20 * SIZE, A1, %xmm10)
	mulps	%xmm13,  %xmm6
	SUBPS	%xmm6,  %xmm1

	movaps	-20 * SIZE(X1), %xmm13

	pshufd	$0xb1, %xmm8,  %xmm4
	mulps	%xmm12,  %xmm8
	addps	%xmm8,  %xmm0
	MOVUPS_A1(-16 * SIZE, A1, %xmm8)
	mulps	%xmm12,  %xmm4
	SUBPS	%xmm4,  %xmm1

	movaps	-16 * SIZE(X1), %xmm12

#ifdef PREFETCHW
	PREFETCH	(PREFETCHSIZE) * 4 - 128 + PREOFFSET(X1)
#endif

	pshufd	$0xb1, %xmm10,  %xmm6
	mulps	%xmm13,  %xmm10
	addps	%xmm10,  %xmm0
	MOVUPS_A1(-12 * SIZE, A1, %xmm10)
	mulps	%xmm13,  %xmm6
	SUBPS	%xmm6,  %xmm1

	movaps	-12 * SIZE(X1), %xmm13

	subq	 $-16 * SIZE, A1
	subq	 $-16 * SIZE, X1

	subq	 $1, I
	BRANCH
	jg	.L123
	ALIGN_3

.L124:
	pshufd	$0xb1, %xmm8,  %xmm4
	mulps	%xmm12,  %xmm8
	addps	%xmm8,  %xmm0
	MOVUPS_A1(-24 * SIZE, A1, %xmm8)
	mulps	%xmm12,  %xmm4
	SUBPS	%xmm4,  %xmm1

	movaps	-24 * SIZE(X1), %xmm12

	pshufd	$0xb1, %xmm10,  %xmm6
	mulps	%xmm13,  %xmm10
	addps	%xmm10,  %xmm0
	MOVUPS_A1(-20 * SIZE, A1, %xmm10)
	mulps	%xmm13,  %xmm6
	SUBPS	%xmm6,  %xmm1

	movaps	-20 * SIZE(X1), %xmm13

	pshufd	$0xb1, %xmm8,  %xmm4
	mulps	%xmm12,  %xmm8
	addps	%xmm8,  %xmm0
	mulps	%xmm12,  %xmm4
	SUBPS	%xmm4,  %xmm1

	movaps	-16 * SIZE(X1), %xmm12

	pshufd	$0xb1, %xmm10,  %xmm6
	mulps	%xmm13,  %xmm10
	addps	%xmm10,  %xmm0
	mulps	%xmm13,  %xmm6
	SUBPS	%xmm6,  %xmm1

	movaps	-12 * SIZE(X1), %xmm13

	subq	 $-16 * SIZE, A1
	subq	 $-16 * SIZE, X1
	ALIGN_3

.L125:
	testq	$4, MM
	je	.L127

	MOVUPS_A1(-32 * SIZE, A1, %xmm8)
	MOVUPS_A1(-28 * SIZE, A1, %xmm10)

	pshufd	$0xb1, %xmm8,  %xmm4
	mulps	%xmm12,  %xmm8
	addps	%xmm8,  %xmm0
	mulps	%xmm12,  %xmm4
	SUBPS	%xmm4,  %xmm1

	movaps	-24 * SIZE(X1), %xmm12

	pshufd	$0xb1, %xmm10,  %xmm6
	mulps	%xmm13,  %xmm10
	addps	%xmm10,  %xmm0
	mulps	%xmm13,  %xmm6
	SUBPS	%xmm6,  %xmm1

	movaps	-20 * SIZE(X1), %xmm13

	addq	 $8 * SIZE, A1
	ALIGN_3

.L127:
	testq	$2, MM
	je	.L128

	MOVUPS_A1(-32 * SIZE, A1, %xmm8)

	pshufd	$0xb1, %xmm8,  %xmm4
	mulps	%xmm12,  %xmm8
	addps	%xmm8,  %xmm0
	mulps	%xmm12,  %xmm4
	SUBPS	%xmm4,  %xmm1

	movaps	%xmm13, %xmm12

	addq	 $4 * SIZE, A1
	ALIGN_3

.L128:
	testq	$1, MM
	je	.L129

#ifdef movsd
	xorps	%xmm8, %xmm8
#endif
	movsd	-32 * SIZE(A1), %xmm8

	pshufd	$0xb1, %xmm8,  %xmm4
	mulps	%xmm12,  %xmm8
	addps	%xmm8,  %xmm0
	mulps	%xmm12,  %xmm4
	SUBPS	%xmm4,  %xmm1
	ALIGN_3

.L129:
	pcmpeqb	%xmm5, %xmm5
	psllq	$63,    %xmm5

#if (!defined(CONJ) && !defined(XCONJ)) || (defined(CONJ) && defined(XCONJ))
	xorps	%xmm5, %xmm0
#else
	xorps	%xmm5, %xmm1
#endif

#ifdef HAVE_SSE3
	haddps	%xmm1, %xmm0
	haddps	%xmm0, %xmm0
#else
	movaps	%xmm0, %xmm8
	unpcklps %xmm1, %xmm0
	unpckhps %xmm1, %xmm8

	addps	%xmm8, %xmm0

	movhlps	%xmm0, %xmm1

	addps	%xmm1, %xmm0
#endif

	pshufd	$0xb1, %xmm0, %xmm1

#ifdef HAVE_SSE3
	movddup	ALPHA, %xmm15
#else
	movsd	ALPHA, %xmm15
	pshufd	$0x44, %xmm15, %xmm15
#endif

	mulps	%xmm15, %xmm0
	mulps	%xmm15, %xmm1

	xorps	%xmm5, %xmm0

#ifdef HAVE_SSE3
	haddps	%xmm1, %xmm0
#else
	movaps	%xmm0, %xmm2
	shufps	$0x88, %xmm1, %xmm0
	shufps	$0xdd, %xmm1, %xmm2
	
	addps	%xmm2, %xmm0
#endif

	movsd	(Y), %xmm12
	addq	INCY, Y

	shufps	$0xd8, %xmm0, %xmm0

	addps	%xmm12, %xmm0

	movlps	%xmm0, (Y1)
	addq	INCY, Y1
	jmp	.L999
	ALIGN_3


.L200:
	testq	$2 * SIZE, LDA
	jne	.L300

	cmpq	$2, N
	jl	.L210
	ALIGN_3

.L201:
	subq	$2, N

	leaq	32 * SIZE(BUFFER), X1

	movq	A, A1
	leaq	(A1, LDA), A2
	leaq	(A1, LDA, 2), A

	xorps	%xmm0, %xmm0
	xorps	%xmm1, %xmm1
	xorps	%xmm2, %xmm2
	xorps	%xmm3, %xmm3

#ifdef ALIGNED_ACCESS
	cmpq	M, MM
	je	.L20X

#ifdef movsd
	xorps	%xmm8, %xmm8
#endif
	movsd	-32 * SIZE(A1), %xmm8
#ifdef movsd
	xorps	%xmm9, %xmm9
#endif
	movsd	-32 * SIZE(A2), %xmm9

#ifdef movsd
	xorps	%xmm12, %xmm12
#endif
	movsd	-32 * SIZE(X1), %xmm12

	pshufd	$0xb1, %xmm8,  %xmm4
	mulps	%xmm12,  %xmm8
	addps	%xmm8,  %xmm0
	mulps	%xmm12,  %xmm4
	SUBPS	%xmm4,  %xmm1

	pshufd	$0xb1, %xmm9, %xmm5
	mulps	%xmm12,  %xmm9
	addps	%xmm9, %xmm2
	mulps	%xmm12,  %xmm5
	SUBPS	%xmm5, %xmm3

	addq	 $2 * SIZE, A1
	addq	 $2 * SIZE, A2
	addq	 $2 * SIZE, X1
	ALIGN_3
.L20X:
#endif

	movaps	-33 * SIZE(A1), %xmm4
	movaps	-33 * SIZE(A2), %xmm5

	movaps	-32 * SIZE(X1), %xmm12
	movaps	-28 * SIZE(X1), %xmm13

#ifdef PREFETCHW
	PREFETCHW	3 * SIZE(Y1)
#endif

	movq	MM,  I
	sarq	$3,  I
	jle	.L205

	movaps	 -29 * SIZE(A1), %xmm6
	movaps	 -29 * SIZE(A2), %xmm7

	decq	 I
	jle	 .L204
	ALIGN_3

.L203:
#ifdef PREFETCH
	PREFETCH	(PREFETCHSIZE) * 2 - 128 + PREOFFSET(A1)
#endif

	movss	 %xmm6, %xmm4
	shufps	 $0x39, %xmm4, %xmm4
	pshufd	$0xb1,  %xmm4,  %xmm14
	mulps	%xmm12, %xmm4
	addps	%xmm4,  %xmm0
	movaps	 -25 * SIZE(A1), %xmm4
	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm1

	movss	 %xmm7, %xmm5
	shufps	 $0x39, %xmm5, %xmm5
	pshufd	$0xb1,  %xmm5,  %xmm14
	mulps	%xmm12, %xmm5
	addps	%xmm5,  %xmm2
	movaps	 -25 * SIZE(A2), %xmm5
	mulps	%xmm12, %xmm14
	movaps	-24 * SIZE(X1), %xmm12
	SUBPS	%xmm14, %xmm3

	movss	 %xmm4, %xmm6
	shufps	 $0x39, %xmm6, %xmm6
	pshufd	$0xb1,  %xmm6,  %xmm14
	mulps	%xmm13, %xmm6
	addps	%xmm6,  %xmm0
	movaps	 -21 * SIZE(A1), %xmm6
	mulps	%xmm13, %xmm14
	SUBPS	%xmm14, %xmm1

	movss	 %xmm5, %xmm7
	shufps	 $0x39, %xmm7, %xmm7
	pshufd	$0xb1,  %xmm7,  %xmm14
	mulps	%xmm13, %xmm7
	addps	%xmm7,  %xmm2
	movaps	 -21 * SIZE(A2), %xmm7
	mulps	%xmm13, %xmm14
	movaps	-20 * SIZE(X1), %xmm13
	SUBPS	%xmm14, %xmm3

#ifdef PREFETCH
	PREFETCH	(PREFETCHSIZE) * 2 - 128 + PREOFFSET(A2)
#endif

	movss	 %xmm6, %xmm4
	shufps	 $0x39, %xmm4, %xmm4
	pshufd	$0xb1,  %xmm4,  %xmm14
	mulps	%xmm12, %xmm4
	addps	%xmm4,  %xmm0
	movaps	 -17 * SIZE(A1), %xmm4
	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm1

	movss	 %xmm7, %xmm5
	shufps	 $0x39, %xmm5, %xmm5
	pshufd	$0xb1,  %xmm5,  %xmm14
	mulps	%xmm12, %xmm5
	addps	%xmm5,  %xmm2
	movaps	 -17 * SIZE(A2), %xmm5
	mulps	%xmm12, %xmm14
	movaps	-16 * SIZE(X1), %xmm12
	SUBPS	%xmm14, %xmm3

#ifdef PREFETCHW
	PREFETCH	(PREFETCHSIZE) * 2 - 128 + PREOFFSET(X1)
#endif

	movss	 %xmm4, %xmm6
	shufps	 $0x39, %xmm6, %xmm6
	pshufd	$0xb1,  %xmm6,  %xmm14
	mulps	%xmm13, %xmm6
	addps	%xmm6,  %xmm0
	movaps	 -13 * SIZE(A1), %xmm6
	mulps	%xmm13, %xmm14
	SUBPS	%xmm14, %xmm1

	movss	 %xmm5, %xmm7
	shufps	 $0x39, %xmm7, %xmm7
	pshufd	$0xb1,  %xmm7,  %xmm14
	mulps	%xmm13, %xmm7
	addps	%xmm7,  %xmm2
	movaps	 -13 * SIZE(A2), %xmm7
	mulps	%xmm13, %xmm14
	movaps	-12 * SIZE(X1), %xmm13
	SUBPS	%xmm14, %xmm3

	subq	 $-16 * SIZE, A1
	subq	 $-16 * SIZE, A2
	subq	 $-16 * SIZE, X1

	subq	 $1, I
	BRANCH
	jg	.L203
	ALIGN_3

.L204:
	movss	 %xmm6, %xmm4
	shufps	 $0x39, %xmm4, %xmm4
	pshufd	$0xb1,  %xmm4,  %xmm14
	mulps	%xmm12, %xmm4
	addps	%xmm4,  %xmm0
	movaps	 -25 * SIZE(A1), %xmm4
	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm1

	movss	 %xmm7, %xmm5
	shufps	 $0x39, %xmm5, %xmm5
	pshufd	$0xb1,  %xmm5,  %xmm14
	mulps	%xmm12, %xmm5
	addps	%xmm5,  %xmm2
	movaps	 -25 * SIZE(A2), %xmm5
	mulps	%xmm12, %xmm14
	movaps	-24 * SIZE(X1), %xmm12
	SUBPS	%xmm14, %xmm3

	movss	 %xmm4, %xmm6
	shufps	 $0x39, %xmm6, %xmm6
	pshufd	$0xb1,  %xmm6,  %xmm14
	mulps	%xmm13, %xmm6
	addps	%xmm6,  %xmm0
	movaps	 -21 * SIZE(A1), %xmm6
	mulps	%xmm13, %xmm14
	SUBPS	%xmm14, %xmm1

	movss	 %xmm5, %xmm7
	shufps	 $0x39, %xmm7, %xmm7
	pshufd	$0xb1,  %xmm7,  %xmm14
	mulps	%xmm13, %xmm7
	addps	%xmm7,  %xmm2
	movaps	 -21 * SIZE(A2), %xmm7
	mulps	%xmm13, %xmm14
	movaps	-20 * SIZE(X1), %xmm13
	SUBPS	%xmm14, %xmm3

	movss	 %xmm6, %xmm4
	shufps	 $0x39, %xmm4, %xmm4
	pshufd	$0xb1,  %xmm4,  %xmm14
	mulps	%xmm12, %xmm4
	addps	%xmm4,  %xmm0
	movaps	 -17 * SIZE(A1), %xmm4
	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm1

	movss	 %xmm7, %xmm5
	shufps	 $0x39, %xmm5, %xmm5
	pshufd	$0xb1,  %xmm5,  %xmm14
	mulps	%xmm12, %xmm5
	addps	%xmm5,  %xmm2
	movaps	 -17 * SIZE(A2), %xmm5
	mulps	%xmm12, %xmm14
	movaps	-16 * SIZE(X1), %xmm12
	SUBPS	%xmm14, %xmm3

	movss	 %xmm4, %xmm6
	shufps	 $0x39, %xmm6, %xmm6
	pshufd	$0xb1,  %xmm6,  %xmm14
	mulps	%xmm13, %xmm6
	addps	%xmm6,  %xmm0
	mulps	%xmm13, %xmm14
	SUBPS	%xmm14, %xmm1

	movss	 %xmm5, %xmm7
	shufps	 $0x39, %xmm7, %xmm7
	pshufd	$0xb1,  %xmm7,  %xmm14
	mulps	%xmm13, %xmm7
	addps	%xmm7,  %xmm2
	mulps	%xmm13, %xmm14
	movaps	-12 * SIZE(X1), %xmm13
	SUBPS	%xmm14, %xmm3

	subq	 $-16 * SIZE, A1
	subq	 $-16 * SIZE, A2
	subq	 $-16 * SIZE, X1
	ALIGN_3

.L205:
	testq	$4, MM
	je	.L207

	movaps	 -29 * SIZE(A1), %xmm6

	movss	 %xmm6, %xmm4
	shufps	 $0x39, %xmm4, %xmm4
	pshufd	$0xb1,  %xmm4,  %xmm14
	mulps	%xmm12, %xmm4
	addps	%xmm4,  %xmm0
	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm1

	movaps	 -29 * SIZE(A2), %xmm7

	movss	 %xmm7, %xmm5
	shufps	 $0x39, %xmm5, %xmm5
	pshufd	$0xb1,  %xmm5,  %xmm14
	mulps	%xmm12, %xmm5
	addps	%xmm5,  %xmm2
	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm3

	movaps	 -25 * SIZE(A1), %xmm8

	movss	 %xmm8, %xmm6
	shufps	 $0x39, %xmm6, %xmm6
	pshufd	$0xb1,  %xmm6,  %xmm14
	mulps	%xmm13, %xmm6
	addps	%xmm6,  %xmm0
	mulps	%xmm13, %xmm14
	SUBPS	%xmm14, %xmm1

	movaps	 -25 * SIZE(A2), %xmm9

	movss	 %xmm9, %xmm7
	shufps	 $0x39, %xmm7, %xmm7
	pshufd	$0xb1,  %xmm7,  %xmm14
	mulps	%xmm13, %xmm7
	addps	%xmm7,  %xmm2
	mulps	%xmm13, %xmm14
	SUBPS	%xmm14, %xmm3

	movaps	%xmm8, %xmm4
	movaps	%xmm9, %xmm5

	movaps	-24 * SIZE(X1), %xmm12
	movaps	-20 * SIZE(X1), %xmm13

	addq	 $8 * SIZE, A1
	addq	 $8 * SIZE, A2
	ALIGN_3

.L207:
	testq	$2, MM
	je	.L208

	movaps	 -29 * SIZE(A1), %xmm6
	movaps	 -29 * SIZE(A2), %xmm7

	movss	 %xmm6, %xmm4
	shufps	 $0x39, %xmm4, %xmm4
	pshufd	$0xb1,  %xmm4,  %xmm14
	mulps	%xmm12, %xmm4
	addps	%xmm4,  %xmm0
	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm1

	movss	 %xmm7, %xmm5
	shufps	 $0x39, %xmm5, %xmm5
	pshufd	$0xb1,  %xmm5,  %xmm14
	mulps	%xmm12, %xmm5
	addps	%xmm5,  %xmm2
	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm3

	movaps	 %xmm6, %xmm4
	movaps	 %xmm7, %xmm5
	movaps	 %xmm13, %xmm12

	addq	 $4 * SIZE, A1
	addq	 $4 * SIZE, A2
	ALIGN_3

.L208:
	testq	$1, MM
	je	.L209

#ifdef movsd
	xorps	%xmm8, %xmm8
#endif
	movsd	-32 * SIZE(A1), %xmm8
#ifdef movsd
	xorps	%xmm9, %xmm9
#endif
	movsd	-32 * SIZE(A2), %xmm9

	pshufd	$0xb1, %xmm8,  %xmm4
	mulps	%xmm12,  %xmm8
	addps	%xmm8,  %xmm0
	mulps	%xmm12,  %xmm4
	SUBPS	%xmm4,  %xmm1

	pshufd	$0xb1, %xmm9, %xmm5
	mulps	%xmm12,  %xmm9
	addps	%xmm9, %xmm2
	mulps	%xmm12,  %xmm5
	SUBPS	%xmm5, %xmm3
	ALIGN_3

.L209:
	pcmpeqb	%xmm5, %xmm5
	psllq	$63,    %xmm5

#if (!defined(CONJ) && !defined(XCONJ)) || (defined(CONJ) && defined(XCONJ))
	xorps	%xmm5, %xmm0
	xorps	%xmm5, %xmm2
#else
	xorps	%xmm5, %xmm1
	xorps	%xmm5, %xmm3
#endif

#ifdef HAVE_SSE3
	haddps	%xmm1, %xmm0
	haddps	%xmm3, %xmm2
	haddps	%xmm2, %xmm0
#else
	movaps	%xmm0, %xmm8
	unpcklps %xmm1, %xmm0
	unpckhps %xmm1, %xmm8

	movaps	%xmm2, %xmm4
	unpcklps %xmm3, %xmm2
	unpckhps %xmm3, %xmm4

	addps	%xmm8, %xmm0
	addps	%xmm4, %xmm2

	movhlps	%xmm0, %xmm1
	movhlps	%xmm2, %xmm3

	addps	%xmm1, %xmm0
	addps	%xmm3, %xmm2
	movlhps %xmm2, %xmm0
#endif

	pshufd	$0xb1, %xmm0, %xmm1

#ifdef HAVE_SSE3
	movddup	ALPHA, %xmm15
#else
	movsd	ALPHA, %xmm15
	pshufd	$0x44, %xmm15, %xmm15
#endif

	mulps	%xmm15, %xmm0
	mulps	%xmm15, %xmm1

	xorps	%xmm5, %xmm0

#ifdef HAVE_SSE3
	haddps	%xmm1, %xmm0
#else
	movaps	%xmm0, %xmm2
	shufps	$0x88, %xmm1, %xmm0
	shufps	$0xdd, %xmm1, %xmm2
	
	addps	%xmm2, %xmm0
#endif

	movsd	(Y), %xmm12
	addq	INCY, Y
	movhps	(Y), %xmm12
	addq	INCY, Y

	shufps	$0xd8, %xmm0, %xmm0

	addps	%xmm12, %xmm0

	movlps	%xmm0, (Y1)
	addq	INCY, Y1
	movhps	%xmm0, (Y1)
	addq	INCY, Y1

	cmpq	$2, N
	jge	.L201
	ALIGN_3

.L210:
	cmpq	$1, N
	jl	.L999

	leaq	32 * SIZE(BUFFER), X1

	movq	A, A1

	xorps	%xmm0, %xmm0
	xorps	%xmm1, %xmm1

#ifdef ALIGNED_ACCESS
	cmpq	M, MM
	je	.L21X

#ifdef movsd
	xorps	%xmm8, %xmm8
#endif
	movsd	-32 * SIZE(A1), %xmm8
#ifdef movsd
	xorps	%xmm12, %xmm12
#endif
	movsd	-32 * SIZE(X1), %xmm12

	pshufd	$0xb1, %xmm8,  %xmm4
	mulps	%xmm12,  %xmm8
	addps	%xmm8,  %xmm0
	mulps	%xmm12,  %xmm4
	SUBPS	%xmm4,  %xmm1

	addq	 $2 * SIZE, A1
	addq	 $2 * SIZE, X1
	ALIGN_3
.L21X:
#endif

	movaps	-33 * SIZE(A1), %xmm4

	movaps	-32 * SIZE(X1), %xmm12
	movaps	-28 * SIZE(X1), %xmm13

	movq	MM,  I
	sarq	$3,  I
	jle	.L215

	movaps	 -29 * SIZE(A1), %xmm5
	movaps	 -25 * SIZE(A1), %xmm6
	movaps	 -21 * SIZE(A1), %xmm7

	decq	 I
	jle	 .L214
	ALIGN_3

.L213:
#ifdef PREFETCH
	PREFETCH	(PREFETCHSIZE) * 4 - 128 + PREOFFSET(A1)
#endif

	movss	 %xmm5, %xmm4
	shufps	 $0x39, %xmm4, %xmm4
	pshufd	$0xb1,  %xmm4,  %xmm14
	mulps	%xmm12, %xmm4
	addps	%xmm4,  %xmm0
	movaps	 -17 * SIZE(A1), %xmm4
	mulps	%xmm12, %xmm14
	movaps	-24 * SIZE(X1), %xmm12
	SUBPS	%xmm14, %xmm1

	movss	 %xmm6, %xmm5
	shufps	 $0x39, %xmm5, %xmm5
	pshufd	$0xb1,  %xmm5,  %xmm15
	mulps	%xmm13, %xmm5
	addps	%xmm5,  %xmm0
	movaps	 -13 * SIZE(A1), %xmm5
	mulps	%xmm13, %xmm15
	movaps	-20 * SIZE(X1), %xmm13
	SUBPS	%xmm15, %xmm1

#ifdef PREFETCHW
	PREFETCH	(PREFETCHSIZE) * 4 - 128 + PREOFFSET(X1)
#endif

	movss	 %xmm7, %xmm6
	shufps	 $0x39, %xmm6, %xmm6
	pshufd	$0xb1,  %xmm6,  %xmm14
	mulps	%xmm12, %xmm6
	addps	%xmm6,  %xmm0
	movaps	  -9 * SIZE(A1), %xmm6
	mulps	%xmm12, %xmm14
	movaps	-16 * SIZE(X1), %xmm12
	SUBPS	%xmm14, %xmm1

	movss	 %xmm4, %xmm7
	shufps	 $0x39, %xmm7, %xmm7
	pshufd	$0xb1,  %xmm7,  %xmm15
	mulps	%xmm13, %xmm7
	addps	%xmm7,  %xmm0
	movaps	  -5 * SIZE(A1), %xmm7
	mulps	%xmm13, %xmm15
	movaps	-12 * SIZE(X1), %xmm13
	SUBPS	%xmm15, %xmm1

	subq	 $-16 * SIZE, A1
	subq	 $-16 * SIZE, X1

	subq	 $1, I
	BRANCH
	jg	.L213
	ALIGN_3

.L214:
	movss	 %xmm5, %xmm4
	shufps	 $0x39, %xmm4, %xmm4
	pshufd	$0xb1,  %xmm4,  %xmm14
	mulps	%xmm12, %xmm4
	addps	%xmm4,  %xmm0
	movaps	 -17 * SIZE(A1), %xmm4
	mulps	%xmm12, %xmm14
	movaps	-24 * SIZE(X1), %xmm12
	SUBPS	%xmm14, %xmm1

	movss	 %xmm6, %xmm5
	shufps	 $0x39, %xmm5, %xmm5
	pshufd	$0xb1,  %xmm5,  %xmm15
	mulps	%xmm13, %xmm5
	addps	%xmm5,  %xmm0
	mulps	%xmm13, %xmm15
	movaps	-20 * SIZE(X1), %xmm13
	SUBPS	%xmm15, %xmm1

	movss	 %xmm7, %xmm6
	shufps	 $0x39, %xmm6, %xmm6
	pshufd	$0xb1,  %xmm6,  %xmm14
	mulps	%xmm12, %xmm6
	addps	%xmm6,  %xmm0
	mulps	%xmm12, %xmm14
	movaps	-16 * SIZE(X1), %xmm12
	SUBPS	%xmm14, %xmm1

	movss	 %xmm4, %xmm7
	shufps	 $0x39, %xmm7, %xmm7
	pshufd	$0xb1,  %xmm7,  %xmm15
	mulps	%xmm13, %xmm7
	addps	%xmm7,  %xmm0
	mulps	%xmm13, %xmm15
	movaps	-12 * SIZE(X1), %xmm13
	SUBPS	%xmm15, %xmm1

	subq	 $-16 * SIZE, A1
	subq	 $-16 * SIZE, X1
	ALIGN_3

.L215:
	testq	$4, MM
	je	.L217

	movaps	 -29 * SIZE(A1), %xmm5
	movaps	 -25 * SIZE(A1), %xmm6

	movss	 %xmm5, %xmm4
	shufps	 $0x39, %xmm4, %xmm4
	pshufd	$0xb1,  %xmm4,  %xmm14
	mulps	%xmm12, %xmm4
	addps	%xmm4,  %xmm0
	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm1

	movss	 %xmm6, %xmm5
	shufps	 $0x39, %xmm5, %xmm5
	pshufd	$0xb1,  %xmm5,  %xmm15
	mulps	%xmm13, %xmm5
	addps	%xmm5,  %xmm0
	mulps	%xmm13, %xmm15
	SUBPS	%xmm15, %xmm1

	movaps	-24 * SIZE(X1), %xmm12
	movaps	-20 * SIZE(X1), %xmm13
	movaps	%xmm6, %xmm4

	addq	 $8 * SIZE, A1
	ALIGN_3

.L217:
	testq	$2, MM
	je	.L218

	movaps	 -29 * SIZE(A1), %xmm5

	movss	 %xmm5, %xmm4
	shufps	 $0x39, %xmm4, %xmm4

	pshufd	$0xb1,  %xmm4,  %xmm14
	mulps	%xmm12, %xmm4
	addps	%xmm4,  %xmm0
	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm1

	movaps	%xmm13, %xmm12

	addq	 $4 * SIZE, A1
	ALIGN_3

.L218:
	testq	$1, MM
	je	.L219

#ifdef movsd
	xorps	%xmm8, %xmm8
#endif
	movsd	-32 * SIZE(A1), %xmm8

	pshufd	$0xb1, %xmm8,  %xmm4
	mulps	%xmm12,  %xmm8
	addps	%xmm8,  %xmm0
	mulps	%xmm12,  %xmm4
	SUBPS	%xmm4,  %xmm1
	ALIGN_3

.L219:
	pcmpeqb	%xmm5, %xmm5
	psllq	$63,    %xmm5

#if (!defined(CONJ) && !defined(XCONJ)) || (defined(CONJ) && defined(XCONJ))
	xorps	%xmm5, %xmm0
#else
	xorps	%xmm5, %xmm1
#endif

#ifdef HAVE_SSE3
	haddps	%xmm1, %xmm0
	haddps	%xmm0, %xmm0
#else
	movaps	%xmm0, %xmm8
	unpcklps %xmm1, %xmm0
	unpckhps %xmm1, %xmm8

	addps	%xmm8, %xmm0

	movhlps	%xmm0, %xmm1

	addps	%xmm1, %xmm0
#endif

	pshufd	$0xb1, %xmm0, %xmm1

#ifdef HAVE_SSE3
	movddup	ALPHA, %xmm15
#else
	movsd	ALPHA, %xmm15
	pshufd	$0x44, %xmm15, %xmm15
#endif

	mulps	%xmm15, %xmm0
	mulps	%xmm15, %xmm1

	xorps	%xmm5, %xmm0

#ifdef HAVE_SSE3
	haddps	%xmm1, %xmm0
#else
	movaps	%xmm0, %xmm2
	shufps	$0x88, %xmm1, %xmm0
	shufps	$0xdd, %xmm1, %xmm2
	
	addps	%xmm2, %xmm0
#endif

	movsd	(Y), %xmm12
	addq	INCY, Y

	shufps	$0xd8, %xmm0, %xmm0

	addps	%xmm12, %xmm0

	movlps	%xmm0, (Y1)
	addq	INCY, Y1
	jmp	.L999

.L300:
	cmpq	$2, N
	jl	.L310
	ALIGN_3

.L301:
	subq	$2, N

	leaq	32 * SIZE(BUFFER), X1

	movq	A, A1
	leaq	(A1, LDA), A2
	leaq	(A1, LDA, 2), A

	xorps	%xmm0, %xmm0
	xorps	%xmm1, %xmm1
	xorps	%xmm2, %xmm2
	xorps	%xmm3, %xmm3

#ifdef ALIGNED_ACCESS
	cmpq	M, MM
	je	.L30X

#ifdef movsd
	xorps	%xmm8, %xmm8
#endif
	movsd	-32 * SIZE(A1), %xmm8
#ifdef movsd
	xorps	%xmm9, %xmm9
#endif
	movsd	-32 * SIZE(A2), %xmm9

#ifdef movsd
	xorps	%xmm12, %xmm12
#endif
	movsd	-32 * SIZE(X1), %xmm12

	pshufd	$0xb1, %xmm8,  %xmm4
	mulps	%xmm12,  %xmm8
	addps	%xmm8,  %xmm0
	mulps	%xmm12,  %xmm4
	SUBPS	%xmm4,  %xmm1

	pshufd	$0xb1, %xmm9, %xmm5
	mulps	%xmm12,  %xmm9
	addps	%xmm9, %xmm2
	mulps	%xmm12,  %xmm5
	SUBPS	%xmm5, %xmm3

	addq	 $2 * SIZE, A1
	addq	 $2 * SIZE, A2
	addq	 $2 * SIZE, X1
	ALIGN_3
.L30X:
#endif

	movaps	-33 * SIZE(A1), %xmm4
	movaps	-35 * SIZE(A2), %xmm5

	movaps	-32 * SIZE(X1), %xmm12
	movaps	-28 * SIZE(X1), %xmm13

#ifdef PREFETCHW
	PREFETCHW	3 * SIZE(Y1)
#endif

	movq	MM,  I
	sarq	$3,  I
	jle	.L305

	movaps	 -29 * SIZE(A1), %xmm6
	movaps	 -31 * SIZE(A2), %xmm7

	decq	 I
	jle	 .L304
	ALIGN_3

.L303:
#ifdef PREFETCH
	PREFETCH	(PREFETCHSIZE) * 2 - 128 + PREOFFSET(A1)
#endif

	movss	 %xmm6, %xmm4
	shufps	 $0x39, %xmm4, %xmm4
	pshufd	$0xb1,  %xmm4,  %xmm14
	mulps	%xmm12, %xmm4
	addps	%xmm4,  %xmm0
	movaps	 -25 * SIZE(A1), %xmm4
	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm1

	movss	 %xmm7, %xmm5
	shufps	 $0x93, %xmm7, %xmm5
	pshufd	$0xb1,  %xmm5,  %xmm14
	mulps	%xmm12, %xmm5
	addps	%xmm5,  %xmm2
	movaps	 -27 * SIZE(A2), %xmm5
	mulps	%xmm12, %xmm14
	movaps	-24 * SIZE(X1), %xmm12
	SUBPS	%xmm14, %xmm3

	movss	 %xmm4, %xmm6
	shufps	 $0x39, %xmm6, %xmm6
	pshufd	$0xb1,  %xmm6,  %xmm14
	mulps	%xmm13, %xmm6
	addps	%xmm6,  %xmm0
	movaps	 -21 * SIZE(A1), %xmm6
	mulps	%xmm13, %xmm14
	SUBPS	%xmm14, %xmm1

	movss	 %xmm5, %xmm7
	shufps	 $0x93, %xmm5, %xmm7
	pshufd	$0xb1,  %xmm7,  %xmm14
	mulps	%xmm13, %xmm7
	addps	%xmm7,  %xmm2
	movaps	 -23 * SIZE(A2), %xmm7
	mulps	%xmm13, %xmm14
	movaps	-20 * SIZE(X1), %xmm13
	SUBPS	%xmm14, %xmm3

#ifdef PREFETCH
	PREFETCH	(PREFETCHSIZE) * 2 - 128 + PREOFFSET(A2)
#endif

	movss	 %xmm6, %xmm4
	shufps	 $0x39, %xmm4, %xmm4
	pshufd	$0xb1,  %xmm4,  %xmm14
	mulps	%xmm12, %xmm4
	addps	%xmm4,  %xmm0
	movaps	 -17 * SIZE(A1), %xmm4
	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm1

	movss	 %xmm7, %xmm5
	shufps	 $0x93, %xmm7, %xmm5
	pshufd	$0xb1,  %xmm5,  %xmm14
	mulps	%xmm12, %xmm5
	addps	%xmm5,  %xmm2
	movaps	 -19 * SIZE(A2), %xmm5
	mulps	%xmm12, %xmm14
	movaps	-16 * SIZE(X1), %xmm12
	SUBPS	%xmm14, %xmm3

#ifdef PREFETCHW
	PREFETCH	(PREFETCHSIZE) * 2 - 128 + PREOFFSET(X1)
#endif

	movss	 %xmm4, %xmm6
	shufps	 $0x39, %xmm6, %xmm6
	pshufd	$0xb1,  %xmm6,  %xmm14
	mulps	%xmm13, %xmm6
	addps	%xmm6,  %xmm0
	movaps	 -13 * SIZE(A1), %xmm6
	mulps	%xmm13, %xmm14
	SUBPS	%xmm14, %xmm1

	movss	 %xmm5, %xmm7
	shufps	 $0x93, %xmm5, %xmm7
	pshufd	$0xb1,  %xmm7,  %xmm14
	mulps	%xmm13, %xmm7
	addps	%xmm7,  %xmm2
	movaps	 -15 * SIZE(A2), %xmm7
	mulps	%xmm13, %xmm14
	movaps	-12 * SIZE(X1), %xmm13
	SUBPS	%xmm14, %xmm3

	subq	 $-16 * SIZE, A1
	subq	 $-16 * SIZE, A2
	subq	 $-16 * SIZE, X1

	subq	 $1, I
	BRANCH
	jg	.L303
	ALIGN_3

.L304:
	movss	 %xmm6, %xmm4
	shufps	 $0x39, %xmm4, %xmm4
	pshufd	$0xb1,  %xmm4,  %xmm14
	mulps	%xmm12, %xmm4
	addps	%xmm4,  %xmm0
	movaps	 -25 * SIZE(A1), %xmm4
	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm1

	movss	 %xmm7, %xmm5
	shufps	 $0x93, %xmm7, %xmm5
	pshufd	$0xb1,  %xmm5,  %xmm14
	mulps	%xmm12, %xmm5
	addps	%xmm5,  %xmm2
	movaps	 -27 * SIZE(A2), %xmm5
	mulps	%xmm12, %xmm14
	movaps	-24 * SIZE(X1), %xmm12
	SUBPS	%xmm14, %xmm3

	movss	 %xmm4, %xmm6
	shufps	 $0x39, %xmm6, %xmm6
	pshufd	$0xb1,  %xmm6,  %xmm14
	mulps	%xmm13, %xmm6
	addps	%xmm6,  %xmm0
	movaps	 -21 * SIZE(A1), %xmm6
	mulps	%xmm13, %xmm14
	SUBPS	%xmm14, %xmm1

	movss	 %xmm5, %xmm7
	shufps	 $0x93, %xmm5, %xmm7
	pshufd	$0xb1,  %xmm7,  %xmm14
	mulps	%xmm13, %xmm7
	addps	%xmm7,  %xmm2
	movaps	 -23 * SIZE(A2), %xmm7
	mulps	%xmm13, %xmm14
	movaps	-20 * SIZE(X1), %xmm13
	SUBPS	%xmm14, %xmm3

	movss	 %xmm6, %xmm4
	shufps	 $0x39, %xmm4, %xmm4
	pshufd	$0xb1,  %xmm4,  %xmm14
	mulps	%xmm12, %xmm4
	addps	%xmm4,  %xmm0
	movaps	 -17 * SIZE(A1), %xmm4
	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm1

	movss	 %xmm7, %xmm5
	shufps	 $0x93, %xmm7, %xmm5
	pshufd	$0xb1,  %xmm5,  %xmm14
	mulps	%xmm12, %xmm5
	addps	%xmm5,  %xmm2
	movaps	 -19 * SIZE(A2), %xmm5
	mulps	%xmm12, %xmm14
	movaps	-16 * SIZE(X1), %xmm12
	SUBPS	%xmm14, %xmm3

	movss	 %xmm4, %xmm6
	shufps	 $0x39, %xmm6, %xmm6
	pshufd	$0xb1,  %xmm6,  %xmm14
	mulps	%xmm13, %xmm6
	addps	%xmm6,  %xmm0
	mulps	%xmm13, %xmm14
	SUBPS	%xmm14, %xmm1

	movss	 %xmm5, %xmm7
	shufps	 $0x93, %xmm5, %xmm7
	pshufd	$0xb1,  %xmm7,  %xmm14
	mulps	%xmm13, %xmm7
	addps	%xmm7,  %xmm2
	mulps	%xmm13, %xmm14
	movaps	-12 * SIZE(X1), %xmm13
	SUBPS	%xmm14, %xmm3

	subq	 $-16 * SIZE, A1
	subq	 $-16 * SIZE, A2
	subq	 $-16 * SIZE, X1
	ALIGN_3

.L305:
	testq	$4, MM
	je	.L307

	movaps	 -29 * SIZE(A1), %xmm6

	movss	 %xmm6, %xmm4
	shufps	 $0x39, %xmm4, %xmm4
	pshufd	$0xb1,  %xmm4,  %xmm14
	mulps	%xmm12, %xmm4
	addps	%xmm4,  %xmm0
	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm1

	movaps	 -31 * SIZE(A2), %xmm7

	movss	 %xmm7, %xmm5
	shufps	 $0x93, %xmm7, %xmm5
	pshufd	$0xb1,  %xmm5,  %xmm14
	mulps	%xmm12, %xmm5
	addps	%xmm5,  %xmm2
	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm3

	movaps	 -25 * SIZE(A1), %xmm8

	movss	 %xmm8, %xmm6
	shufps	 $0x39, %xmm6, %xmm6
	pshufd	$0xb1,  %xmm6,  %xmm14
	mulps	%xmm13, %xmm6
	addps	%xmm6,  %xmm0
	mulps	%xmm13, %xmm14
	SUBPS	%xmm14, %xmm1

	movaps	 -27 * SIZE(A2), %xmm9

	movss	 %xmm9, %xmm7
	shufps	 $0x93, %xmm9, %xmm7
	pshufd	$0xb1,  %xmm7,  %xmm14
	mulps	%xmm13, %xmm7
	addps	%xmm7,  %xmm2
	mulps	%xmm13, %xmm14
	SUBPS	%xmm14, %xmm3

	movaps	%xmm8, %xmm4
	movaps	%xmm9, %xmm5

	movaps	-24 * SIZE(X1), %xmm12
	movaps	-20 * SIZE(X1), %xmm13

	addq	 $8 * SIZE, A1
	addq	 $8 * SIZE, A2
	ALIGN_3

.L307:
	testq	$2, MM
	je	.L308

	movaps	 -29 * SIZE(A1), %xmm6
	movaps	 -31 * SIZE(A2), %xmm7

	movss	 %xmm6, %xmm4
	shufps	 $0x39, %xmm4, %xmm4
	pshufd	$0xb1,  %xmm4,  %xmm14
	mulps	%xmm12, %xmm4
	addps	%xmm4,  %xmm0
	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm1

	movss	 %xmm7, %xmm5
	shufps	 $0x93, %xmm7, %xmm5
	pshufd	$0xb1,  %xmm5,  %xmm14
	mulps	%xmm12, %xmm5
	addps	%xmm5,  %xmm2
	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm3

	movaps	 %xmm6, %xmm4
	movaps	 %xmm7, %xmm5
	movaps	 %xmm13, %xmm12

	addq	 $4 * SIZE, A1
	addq	 $4 * SIZE, A2
	ALIGN_3

.L308:
	testq	$1, MM
	je	.L309

#ifdef movsd
	xorps	%xmm8, %xmm8
#endif
	movsd	-32 * SIZE(A1), %xmm8
#ifdef movsd
	xorps	%xmm9, %xmm9
#endif
	movsd	-32 * SIZE(A2), %xmm9

	pshufd	$0xb1, %xmm8,  %xmm4
	mulps	%xmm12,  %xmm8
	addps	%xmm8,  %xmm0
	mulps	%xmm12,  %xmm4
	SUBPS	%xmm4,  %xmm1

	pshufd	$0xb1, %xmm9, %xmm5
	mulps	%xmm12,  %xmm9
	addps	%xmm9, %xmm2
	mulps	%xmm12,  %xmm5
	SUBPS	%xmm5, %xmm3
	ALIGN_3

.L309:
	pcmpeqb	%xmm5, %xmm5
	psllq	$63,    %xmm5

#if (!defined(CONJ) && !defined(XCONJ)) || (defined(CONJ) && defined(XCONJ))
	xorps	%xmm5, %xmm0
	xorps	%xmm5, %xmm2
#else
	xorps	%xmm5, %xmm1
	xorps	%xmm5, %xmm3
#endif

#ifdef HAVE_SSE3
	haddps	%xmm1, %xmm0
	haddps	%xmm3, %xmm2
	haddps	%xmm2, %xmm0
#else
	movaps	%xmm0, %xmm8
	unpcklps %xmm1, %xmm0
	unpckhps %xmm1, %xmm8

	movaps	%xmm2, %xmm4
	unpcklps %xmm3, %xmm2
	unpckhps %xmm3, %xmm4

	addps	%xmm8, %xmm0
	addps	%xmm4, %xmm2

	movhlps	%xmm0, %xmm1
	movhlps	%xmm2, %xmm3

	addps	%xmm1, %xmm0
	addps	%xmm3, %xmm2
	movlhps %xmm2, %xmm0
#endif

	pshufd	$0xb1, %xmm0, %xmm1

#ifdef HAVE_SSE3
	movddup	ALPHA, %xmm15
#else
	movsd	ALPHA, %xmm15
	pshufd	$0x44, %xmm15, %xmm15
#endif

	mulps	%xmm15, %xmm0
	mulps	%xmm15, %xmm1

	xorps	%xmm5, %xmm0

#ifdef HAVE_SSE3
	haddps	%xmm1, %xmm0
#else
	movaps	%xmm0, %xmm2
	shufps	$0x88, %xmm1, %xmm0
	shufps	$0xdd, %xmm1, %xmm2
	
	addps	%xmm2, %xmm0
#endif

	movsd	(Y), %xmm12
	addq	INCY, Y
	movhps	(Y), %xmm12
	addq	INCY, Y

	shufps	$0xd8, %xmm0, %xmm0

	addps	%xmm12, %xmm0

	movlps	%xmm0, (Y1)
	addq	INCY, Y1
	movhps	%xmm0, (Y1)
	addq	INCY, Y1

	cmpq	$2, N
	jge	.L301
	ALIGN_3

.L310:
	cmpq	$1, N
	jl	.L999

	leaq	32 * SIZE(BUFFER), X1

	movq	A, A1

	xorps	%xmm0, %xmm0
	xorps	%xmm1, %xmm1

#ifdef ALIGNED_ACCESS
	cmpq	M, MM
	je	.L31X

#ifdef movsd
	xorps	%xmm8, %xmm8
#endif
	movsd	-32 * SIZE(A1), %xmm8
#ifdef movsd
	xorps	%xmm12, %xmm12
#endif
	movsd	-32 * SIZE(X1), %xmm12

	pshufd	$0xb1, %xmm8,  %xmm4
	mulps	%xmm12,  %xmm8
	addps	%xmm8,  %xmm0
	mulps	%xmm12,  %xmm4
	SUBPS	%xmm4,  %xmm1

	addq	 $2 * SIZE, A1
	addq	 $2 * SIZE, X1
	ALIGN_3
.L31X:
#endif

	movaps	-33 * SIZE(A1), %xmm4

	movaps	-32 * SIZE(X1), %xmm12
	movaps	-28 * SIZE(X1), %xmm13

	movq	MM,  I
	sarq	$3,  I
	jle	.L315

	movaps	 -29 * SIZE(A1), %xmm5
	movaps	 -25 * SIZE(A1), %xmm6
	movaps	 -21 * SIZE(A1), %xmm7

	decq	 I
	jle	 .L314
	ALIGN_3

.L313:
#ifdef PREFETCH
	PREFETCH	(PREFETCHSIZE) * 4 - 128 + PREOFFSET(A1)
#endif

	movss	 %xmm5, %xmm4
	shufps	 $0x39, %xmm4, %xmm4
	pshufd	$0xb1,  %xmm4,  %xmm14
	mulps	%xmm12, %xmm4
	addps	%xmm4,  %xmm0
	movaps	 -17 * SIZE(A1), %xmm4
	mulps	%xmm12, %xmm14
	movaps	-24 * SIZE(X1), %xmm12
	SUBPS	%xmm14, %xmm1

	movss	 %xmm6, %xmm5
	shufps	 $0x39, %xmm5, %xmm5
	pshufd	$0xb1,  %xmm5,  %xmm15
	mulps	%xmm13, %xmm5
	addps	%xmm5,  %xmm0
	movaps	 -13 * SIZE(A1), %xmm5
	mulps	%xmm13, %xmm15
	movaps	-20 * SIZE(X1), %xmm13
	SUBPS	%xmm15, %xmm1

#ifdef PREFETCHW
	PREFETCH	(PREFETCHSIZE) * 4 - 128 + PREOFFSET(X1)
#endif

	movss	 %xmm7, %xmm6
	shufps	 $0x39, %xmm6, %xmm6
	pshufd	$0xb1,  %xmm6,  %xmm14
	mulps	%xmm12, %xmm6
	addps	%xmm6,  %xmm0
	movaps	  -9 * SIZE(A1), %xmm6
	mulps	%xmm12, %xmm14
	movaps	-16 * SIZE(X1), %xmm12
	SUBPS	%xmm14, %xmm1

	movss	 %xmm4, %xmm7
	shufps	 $0x39, %xmm7, %xmm7
	pshufd	$0xb1,  %xmm7,  %xmm15
	mulps	%xmm13, %xmm7
	addps	%xmm7,  %xmm0
	movaps	  -5 * SIZE(A1), %xmm7
	mulps	%xmm13, %xmm15
	movaps	-12 * SIZE(X1), %xmm13
	SUBPS	%xmm15, %xmm1

	subq	 $-16 * SIZE, A1
	subq	 $-16 * SIZE, X1

	subq	 $1, I
	BRANCH
	jg	.L313
	ALIGN_3

.L314:
	movss	 %xmm5, %xmm4
	shufps	 $0x39, %xmm4, %xmm4
	pshufd	$0xb1,  %xmm4,  %xmm14
	mulps	%xmm12, %xmm4
	addps	%xmm4,  %xmm0
	movaps	 -17 * SIZE(A1), %xmm4
	mulps	%xmm12, %xmm14
	movaps	-24 * SIZE(X1), %xmm12
	SUBPS	%xmm14, %xmm1

	movss	 %xmm6, %xmm5
	shufps	 $0x39, %xmm5, %xmm5
	pshufd	$0xb1,  %xmm5,  %xmm15
	mulps	%xmm13, %xmm5
	addps	%xmm5,  %xmm0
	mulps	%xmm13, %xmm15
	movaps	-20 * SIZE(X1), %xmm13
	SUBPS	%xmm15, %xmm1

	movss	 %xmm7, %xmm6
	shufps	 $0x39, %xmm6, %xmm6
	pshufd	$0xb1,  %xmm6,  %xmm14
	mulps	%xmm12, %xmm6
	addps	%xmm6,  %xmm0
	mulps	%xmm12, %xmm14
	movaps	-16 * SIZE(X1), %xmm12
	SUBPS	%xmm14, %xmm1

	movss	 %xmm4, %xmm7
	shufps	 $0x39, %xmm7, %xmm7
	pshufd	$0xb1,  %xmm7,  %xmm15
	mulps	%xmm13, %xmm7
	addps	%xmm7,  %xmm0
	mulps	%xmm13, %xmm15
	movaps	-12 * SIZE(X1), %xmm13
	SUBPS	%xmm15, %xmm1

	subq	 $-16 * SIZE, A1
	subq	 $-16 * SIZE, X1
	ALIGN_3

.L315:
	testq	$4, MM
	je	.L317

	movaps	 -29 * SIZE(A1), %xmm5
	movaps	 -25 * SIZE(A1), %xmm6

	movss	 %xmm5, %xmm4
	shufps	 $0x39, %xmm4, %xmm4
	pshufd	$0xb1,  %xmm4,  %xmm14
	mulps	%xmm12, %xmm4
	addps	%xmm4,  %xmm0
	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm1

	movss	 %xmm6, %xmm5
	shufps	 $0x39, %xmm5, %xmm5
	pshufd	$0xb1,  %xmm5,  %xmm15
	mulps	%xmm13, %xmm5
	addps	%xmm5,  %xmm0
	mulps	%xmm13, %xmm15
	SUBPS	%xmm15, %xmm1

	movaps	-24 * SIZE(X1), %xmm12
	movaps	-20 * SIZE(X1), %xmm13
	movaps	%xmm6, %xmm4

	addq	 $8 * SIZE, A1
	ALIGN_3

.L317:
	testq	$2, MM
	je	.L318

	movaps	 -29 * SIZE(A1), %xmm5

	movss	 %xmm5, %xmm4
	shufps	 $0x39, %xmm4, %xmm4

	pshufd	$0xb1,  %xmm4,  %xmm14
	mulps	%xmm12, %xmm4
	addps	%xmm4,  %xmm0
	mulps	%xmm12, %xmm14
	SUBPS	%xmm14, %xmm1

	movaps	%xmm13, %xmm12

	addq	 $4 * SIZE, A1
	ALIGN_3

.L318:
	testq	$1, MM
	je	.L319

#ifdef movsd
	xorps	%xmm8, %xmm8
#endif
	movsd	-32 * SIZE(A1), %xmm8

	pshufd	$0xb1, %xmm8,  %xmm4
	mulps	%xmm12,  %xmm8
	addps	%xmm8,  %xmm0
	mulps	%xmm12,  %xmm4
	SUBPS	%xmm4,  %xmm1
	ALIGN_3

.L319:
	pcmpeqb	%xmm5, %xmm5
	psllq	$63,    %xmm5

#if (!defined(CONJ) && !defined(XCONJ)) || (defined(CONJ) && defined(XCONJ))
	xorps	%xmm5, %xmm0
#else
	xorps	%xmm5, %xmm1
#endif

#ifdef HAVE_SSE3
	haddps	%xmm1, %xmm0
	haddps	%xmm0, %xmm0
#else
	movaps	%xmm0, %xmm8
	unpcklps %xmm1, %xmm0
	unpckhps %xmm1, %xmm8

	addps	%xmm8, %xmm0

	movhlps	%xmm0, %xmm1

	addps	%xmm1, %xmm0
#endif

	pshufd	$0xb1, %xmm0, %xmm1

#ifdef HAVE_SSE3
	movddup	ALPHA, %xmm15
#else
	movsd	ALPHA, %xmm15
	pshufd	$0x44, %xmm15, %xmm15
#endif

	mulps	%xmm15, %xmm0
	mulps	%xmm15, %xmm1

	xorps	%xmm5, %xmm0

#ifdef HAVE_SSE3
	haddps	%xmm1, %xmm0
#else
	movaps	%xmm0, %xmm2
	shufps	$0x88, %xmm1, %xmm0
	shufps	$0xdd, %xmm1, %xmm2
	
	addps	%xmm2, %xmm0
#endif

	movsd	(Y), %xmm12
	addq	INCY, Y

	shufps	$0xd8, %xmm0, %xmm0

	addps	%xmm12, %xmm0

	movlps	%xmm0, (Y1)
	addq	INCY, Y1
#endif
	ALIGN_3

.L999:
	movq	  0(%rsp), %rbx
	movq	  8(%rsp), %rbp
	movq	 16(%rsp), %r12
	movq	 24(%rsp), %r13
	movq	 32(%rsp), %r14
	movq	 40(%rsp), %r15

#ifdef WINDOWS_ABI
	movq	 48(%rsp), %rdi
	movq	 56(%rsp), %rsi
	movups	 64(%rsp), %xmm6
	movups	 80(%rsp), %xmm7
	movups	 96(%rsp), %xmm8
	movups	112(%rsp), %xmm9
	movups	128(%rsp), %xmm10
	movups	144(%rsp), %xmm11
	movups	160(%rsp), %xmm12
	movups	176(%rsp), %xmm13
	movups	192(%rsp), %xmm14
	movups	208(%rsp), %xmm15
#endif

	addq	$STACKSIZE, %rsp
	ret

	EPILOGUE
