OSchip
/
OpenBLAS
Not watched
1
0
Fork 0
Code
Releases 66 Wiki evaluate Activity
Issues 0 Pull Requests 0 Datasets Model Cloudbrain HPC
Browse Source

updated dgemv_n kernel for nehalem

tags/v0.2.15^2
Werner Saar 10 years ago
parent
30f52d53df
commit
133c11a156
1 changed files with 72 additions and 165 deletions
Unified View
Diff Options
Show Stats Download Patch File Download Diff File
  1. +72
    -165
      kernel/x86_64/dgemv_n_microk_nehalem-4.c

+ 72
- 165
kernel/x86_64/dgemv_n_microk_nehalem-4.c View File

@@ -27,10 +27,10 @@ USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.






#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; BLASLONG register i = 0;
@@ -46,200 +46,107 @@ static void dgemv_kernel_4x8( BLASLONG n, FLOAT **ap, FLOAT *x, FLOAT *y, BLASLO
"shufpd $0, %%xmm14, %%xmm14\n\t" "shufpd $0, %%xmm14, %%xmm14\n\t"
"shufpd $0, %%xmm15, %%xmm15\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" "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 (%4,%0,8), %%xmm8 \n\t"
"movups 16(%4,%0,8), %%xmm0 \n\t"
"movups (%5,%0,8), %%xmm9 \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 (%6,%0,8), %%xmm10 \n\t"
"movups 16(%6,%0,8), %%xmm2 \n\t"
"movups (%7,%0,8), %%xmm11 \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 %%xmm4 , %%xmm4 \n\t"
"xorpd %%xmm5 , %%xmm5 \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 %%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 %%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 %%xmm4 , %%xmm4 \n\t"
"xorpd %%xmm5 , %%xmm5 \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 %%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 %%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 %%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 %%xmm9 , %%xmm4 \n\t"
"mulpd %%xmm14, %%xmm10 \n\t"
"addpd %%xmm10 , %%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, 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


: :
: :
@@ -253,8 +160,8 @@ static void dgemv_kernel_4x4( BLASLONG n, FLOAT **ap, FLOAT *x, FLOAT *y, FLOAT
"r" (ap[3]), // 7 "r" (ap[3]), // 7
"r" (alpha) // 8 "r" (alpha) // 8
: "cc", : "cc",
"%xmm4", "%xmm5",
"%xmm6", "%xmm7",
"%xmm0", "%xmm1", "%xmm2", "%xmm3",
"%xmm4", "%xmm5", "%xmm6", "%xmm7",
"%xmm8", "%xmm9", "%xmm10", "%xmm11", "%xmm8", "%xmm9", "%xmm10", "%xmm11",
"%xmm12", "%xmm13", "%xmm14", "%xmm15", "%xmm12", "%xmm13", "%xmm14", "%xmm15",
"memory" "memory"


Write Preview
Loading…
Cancel
Save