-#define HAVE_KERNEL_4x8 1
-static void dgemv_kernel_4x8( BLASLONG n, FLOAT **ap, FLOAT *x, FLOAT *y, BLASLONG lda4, FLOAT *alpha) __attribute__ ((noinline));
+#define HAVE_KERNEL_4x4 1
+static void dgemv_kernel_4x4( BLASLONG n, FLOAT **ap, FLOAT *x, FLOAT *y, FLOAT *alpha) __attribute__ ((noinline));
-static void dgemv_kernel_4x8( BLASLONG n, FLOAT **ap, FLOAT *x, FLOAT *y, BLASLONG lda4, FLOAT *alpha)
+static void dgemv_kernel_4x4( BLASLONG n, FLOAT **ap, FLOAT *x, FLOAT *y, FLOAT *alpha)
{
BLASLONG register i = 0;
"shufpd $0, %%xmm14, %%xmm14\n\t"
"shufpd $0, %%xmm15, %%xmm15\n\t"
- "movsd 32(%2), %%xmm0 \n\t" // x4
- "movsd 40(%2), %%xmm1 \n\t" // x5
- "movsd 48(%2), %%xmm2 \n\t" // x6
- "movsd 56(%2), %%xmm3 \n\t" // x7
- "shufpd $0, %%xmm0 , %%xmm0 \n\t"
- "shufpd $0, %%xmm1 , %%xmm1 \n\t"
- "shufpd $0, %%xmm2 , %%xmm2 \n\t"
- "shufpd $0, %%xmm3 , %%xmm3 \n\t"
-
- "movsd (%9), %%xmm6 \n\t" // alpha
+ "movsd (%8), %%xmm6 \n\t" // alpha
"shufpd $0, %%xmm6 , %%xmm6 \n\t"
-
- ".align 16 \n\t"
- "1: \n\t"
- "xorpd %%xmm4 , %%xmm4 \n\t"
- "xorpd %%xmm5 , %%xmm5 \n\t"
- "movups (%3,%0,8), %%xmm7 \n\t" // 2 * y
-
- ".align 2 \n\t"
"movups (%4,%0,8), %%xmm8 \n\t"
+ "movups 16(%4,%0,8), %%xmm0 \n\t"
"movups (%5,%0,8), %%xmm9 \n\t"
+ "movups 16(%5,%0,8), %%xmm1 \n\t"
"movups (%6,%0,8), %%xmm10 \n\t"
+ "movups 16(%6,%0,8), %%xmm2 \n\t"
"movups (%7,%0,8), %%xmm11 \n\t"
- ".align 2 \n\t"
- "mulpd %%xmm12, %%xmm8 \n\t"
- "mulpd %%xmm13, %%xmm9 \n\t"
- "mulpd %%xmm14, %%xmm10 \n\t"
- "mulpd %%xmm15, %%xmm11 \n\t"
- "addpd %%xmm8 , %%xmm4 \n\t"
- "addpd %%xmm9 , %%xmm5 \n\t"
- "addpd %%xmm10, %%xmm4 \n\t"
- "addpd %%xmm11, %%xmm5 \n\t"
-
- "movups (%4,%8,8), %%xmm8 \n\t"
- "movups (%5,%8,8), %%xmm9 \n\t"
- "movups (%6,%8,8), %%xmm10 \n\t"
- "movups (%7,%8,8), %%xmm11 \n\t"
- ".align 2 \n\t"
- "mulpd %%xmm0 , %%xmm8 \n\t"
- "mulpd %%xmm1 , %%xmm9 \n\t"
- "mulpd %%xmm2 , %%xmm10 \n\t"
- "mulpd %%xmm3 , %%xmm11 \n\t"
- "addpd %%xmm8 , %%xmm4 \n\t"
- "addpd %%xmm9 , %%xmm5 \n\t"
- "addpd %%xmm10, %%xmm4 \n\t"
- "addpd %%xmm11, %%xmm5 \n\t"
+ "movups 16(%7,%0,8), %%xmm3 \n\t"
- "addpd %%xmm5 , %%xmm4 \n\t"
- "mulpd %%xmm6 , %%xmm4 \n\t"
- "addpd %%xmm4 , %%xmm7 \n\t"
+ "addq $4 , %0 \n\t"
+ "subq $4 , %1 \n\t"
+ "jz 2f \n\t"
- "movups %%xmm7 , (%3,%0,8) \n\t" // 2 * y
+ ".align 16 \n\t"
+ "1: \n\t"
"xorpd %%xmm4 , %%xmm4 \n\t"
"xorpd %%xmm5 , %%xmm5 \n\t"
- "movups 16(%3,%0,8), %%xmm7 \n\t" // 2 * y
-
- ".align 2 \n\t"
- "movups 16(%4,%0,8), %%xmm8 \n\t"
- "movups 16(%5,%0,8), %%xmm9 \n\t"
- "movups 16(%6,%0,8), %%xmm10 \n\t"
- "movups 16(%7,%0,8), %%xmm11 \n\t"
- ".align 2 \n\t"
+ "movups -32(%3,%0,8), %%xmm7 \n\t" // 2 * y
+
"mulpd %%xmm12, %%xmm8 \n\t"
- "mulpd %%xmm13, %%xmm9 \n\t"
- "mulpd %%xmm14, %%xmm10 \n\t"
- "mulpd %%xmm15, %%xmm11 \n\t"
- "addpd %%xmm8 , %%xmm4 \n\t"
- "addpd %%xmm9 , %%xmm5 \n\t"
- "addpd %%xmm10, %%xmm4 \n\t"
- "addpd %%xmm11, %%xmm5 \n\t"
-
- "movups 16(%4,%8,8), %%xmm8 \n\t"
- "movups 16(%5,%8,8), %%xmm9 \n\t"
- "movups 16(%6,%8,8), %%xmm10 \n\t"
- "movups 16(%7,%8,8), %%xmm11 \n\t"
- ".align 2 \n\t"
- "mulpd %%xmm0 , %%xmm8 \n\t"
- "mulpd %%xmm1 , %%xmm9 \n\t"
- "mulpd %%xmm2 , %%xmm10 \n\t"
- "mulpd %%xmm3 , %%xmm11 \n\t"
+ "mulpd %%xmm12, %%xmm0 \n\t"
"addpd %%xmm8 , %%xmm4 \n\t"
- "addpd %%xmm9 , %%xmm5 \n\t"
- "addpd %%xmm10, %%xmm4 \n\t"
- "addpd %%xmm11, %%xmm5 \n\t"
+ "addpd %%xmm0 , %%xmm5 \n\t"
- "addq $4 , %8 \n\t"
- "addpd %%xmm5 , %%xmm4 \n\t"
- "mulpd %%xmm6 , %%xmm4 \n\t"
- "addpd %%xmm4 , %%xmm7 \n\t"
-
- "movups %%xmm7 , 16(%3,%0,8) \n\t" // 2 * y
+ "movups (%4,%0,8), %%xmm8 \n\t"
+ "movups 16(%4,%0,8), %%xmm0 \n\t"
- "addq $4 , %0 \n\t"
- "subq $4 , %1 \n\t"
- "jnz 1b \n\t"
+ "mulpd %%xmm13, %%xmm9 \n\t"
+ "mulpd %%xmm13, %%xmm1 \n\t"
+ "addpd %%xmm9 , %%xmm4 \n\t"
+ "addpd %%xmm1 , %%xmm5 \n\t"
- :
- :
- "r" (i), // 0
- "r" (n), // 1
- "r" (x), // 2
- "r" (y), // 3
- "r" (ap[0]), // 4
- "r" (ap[1]), // 5
- "r" (ap[2]), // 6
- "r" (ap[3]), // 7
- "r" (lda4), // 8
- "r" (alpha) // 9
- : "cc",
- "%xmm0", "%xmm1",
- "%xmm2", "%xmm3",
- "%xmm4", "%xmm5",
- "%xmm6", "%xmm7",
- "%xmm8", "%xmm9", "%xmm10", "%xmm11",
- "%xmm12", "%xmm13", "%xmm14", "%xmm15",
- "memory"
- );
+ "movups (%5,%0,8), %%xmm9 \n\t"
+ "movups 16(%5,%0,8), %%xmm1 \n\t"
-}
+ "mulpd %%xmm14, %%xmm10 \n\t"
+ "mulpd %%xmm14, %%xmm2 \n\t"
+ "addpd %%xmm10 , %%xmm4 \n\t"
+ "addpd %%xmm2 , %%xmm5 \n\t"
+ "movups (%6,%0,8), %%xmm10 \n\t"
+ "movups 16(%6,%0,8), %%xmm2 \n\t"
+ "mulpd %%xmm15, %%xmm11 \n\t"
+ "mulpd %%xmm15, %%xmm3 \n\t"
+ "addpd %%xmm11 , %%xmm4 \n\t"
+ "addpd %%xmm3 , %%xmm5 \n\t"
+ "movups (%7,%0,8), %%xmm11 \n\t"
+ "movups 16(%7,%0,8), %%xmm3 \n\t"
-#define HAVE_KERNEL_4x4 1
-static void dgemv_kernel_4x4( BLASLONG n, FLOAT **ap, FLOAT *x, FLOAT *y, FLOAT *alpha) __attribute__ ((noinline));
-static void dgemv_kernel_4x4( BLASLONG n, FLOAT **ap, FLOAT *x, FLOAT *y, FLOAT *alpha)
-{
+ "mulpd %%xmm6 , %%xmm4 \n\t"
+ "addpd %%xmm7 , %%xmm4 \n\t"
+ "movups -16(%3,%0,8), %%xmm7 \n\t" // 2 * y
+ "movups %%xmm4 , -32(%3,%0,8) \n\t" // 2 * y
- BLASLONG register i = 0;
+ "mulpd %%xmm6 , %%xmm5 \n\t"
+ "addpd %%xmm7 , %%xmm5 \n\t"
+ "movups %%xmm5 , -16(%3,%0,8) \n\t" // 2 * y
- __asm__ __volatile__
- (
- "movsd (%2), %%xmm12 \n\t" // x0
- "movsd 8(%2), %%xmm13 \n\t" // x1
- "movsd 16(%2), %%xmm14 \n\t" // x2
- "movsd 24(%2), %%xmm15 \n\t" // x3
- "shufpd $0, %%xmm12, %%xmm12\n\t"
- "shufpd $0, %%xmm13, %%xmm13\n\t"
- "shufpd $0, %%xmm14, %%xmm14\n\t"
- "shufpd $0, %%xmm15, %%xmm15\n\t"
+ "addq $4 , %0 \n\t"
+ "subq $4 , %1 \n\t"
+ "jnz 1b \n\t"
- "movsd (%8), %%xmm6 \n\t" // alpha
- "shufpd $0, %%xmm6 , %%xmm6 \n\t"
+ "2: \n\t"
- ".align 16 \n\t"
- "1: \n\t"
"xorpd %%xmm4 , %%xmm4 \n\t"
"xorpd %%xmm5 , %%xmm5 \n\t"
- "movups (%3,%0,8), %%xmm7 \n\t" // 2 * y
- "movups (%4,%0,8), %%xmm8 \n\t"
- "movups (%5,%0,8), %%xmm9 \n\t"
- "movups (%6,%0,8), %%xmm10 \n\t"
- "movups (%7,%0,8), %%xmm11 \n\t"
"mulpd %%xmm12, %%xmm8 \n\t"
- "mulpd %%xmm13, %%xmm9 \n\t"
- "mulpd %%xmm14, %%xmm10 \n\t"
- "mulpd %%xmm15, %%xmm11 \n\t"
"addpd %%xmm8 , %%xmm4 \n\t"
- "addpd %%xmm9 , %%xmm4 \n\t"
- "addpd %%xmm10 , %%xmm4 \n\t"
- "addpd %%xmm4 , %%xmm11 \n\t"
-
- "mulpd %%xmm6 , %%xmm11 \n\t"
- "addpd %%xmm7 , %%xmm11 \n\t"
- "movups %%xmm11, (%3,%0,8) \n\t" // 2 * y
-
- "xorpd %%xmm4 , %%xmm4 \n\t"
- "xorpd %%xmm5 , %%xmm5 \n\t"
- "movups 16(%3,%0,8), %%xmm7 \n\t" // 2 * y
-
- "movups 16(%4,%0,8), %%xmm8 \n\t"
- "movups 16(%5,%0,8), %%xmm9 \n\t"
- "movups 16(%6,%0,8), %%xmm10 \n\t"
- "movups 16(%7,%0,8), %%xmm11 \n\t"
- "mulpd %%xmm12, %%xmm8 \n\t"
"mulpd %%xmm13, %%xmm9 \n\t"
- "mulpd %%xmm14, %%xmm10 \n\t"
- "mulpd %%xmm15, %%xmm11 \n\t"
- "addpd %%xmm8 , %%xmm4 \n\t"
"addpd %%xmm9 , %%xmm4 \n\t"
+ "mulpd %%xmm14, %%xmm10 \n\t"
"addpd %%xmm10 , %%xmm4 \n\t"
- "addpd %%xmm4 , %%xmm11 \n\t"
-
- "mulpd %%xmm6 , %%xmm11 \n\t"
- "addpd %%xmm7 , %%xmm11 \n\t"
- "movups %%xmm11, 16(%3,%0,8) \n\t" // 2 * y
-
- "addq $4 , %0 \n\t"
- "subq $4 , %1 \n\t"
- "jnz 1b \n\t"
+ "mulpd %%xmm15, %%xmm11 \n\t"
+ "addpd %%xmm11 , %%xmm4 \n\t"
+
+ "mulpd %%xmm12, %%xmm0 \n\t"
+ "addpd %%xmm0 , %%xmm5 \n\t"
+ "mulpd %%xmm13, %%xmm1 \n\t"
+ "addpd %%xmm1 , %%xmm5 \n\t"
+ "mulpd %%xmm14, %%xmm2 \n\t"
+ "addpd %%xmm2 , %%xmm5 \n\t"
+ "mulpd %%xmm15, %%xmm3 \n\t"
+ "addpd %%xmm3 , %%xmm5 \n\t"
+
+ "movups -32(%3,%0,8), %%xmm7 \n\t" // 2 * y
+ "mulpd %%xmm6 , %%xmm4 \n\t"
+ "addpd %%xmm7 , %%xmm4 \n\t"
+ "movups %%xmm4 , -32(%3,%0,8) \n\t" // 2 * y
+
+ "movups -16(%3,%0,8), %%xmm7 \n\t" // 2 * y
+ "mulpd %%xmm6 , %%xmm5 \n\t"
+ "addpd %%xmm7 , %%xmm5 \n\t"
+ "movups %%xmm5 , -16(%3,%0,8) \n\t" // 2 * y
:
:
"r" (ap[3]), // 7
"r" (alpha) // 8
: "cc",
- "%xmm4", "%xmm5",
- "%xmm6", "%xmm7",
+ "%xmm0", "%xmm1", "%xmm2", "%xmm3",
+ "%xmm4", "%xmm5", "%xmm6", "%xmm7",
"%xmm8", "%xmm9", "%xmm10", "%xmm11",
"%xmm12", "%xmm13", "%xmm14", "%xmm15",
"memory"