📄 trsm_kernel_ln_4x2_sse2.s
字号:
#if defined(LN) || defined(LT) movapd %xmm4, %xmm0 unpcklpd %xmm5, %xmm4 unpckhpd %xmm5, %xmm0 movapd 0 * SIZE(B), %xmm2 movapd 2 * SIZE(B), %xmm3 subpd %xmm4, %xmm2 subpd %xmm0, %xmm3#else movapd 0 * SIZE(AA), %xmm0 movapd 2 * SIZE(AA), %xmm1 subpd %xmm4, %xmm0 subpd %xmm5, %xmm1#endif#ifdef LN movsd 3 * SIZE(AA), %xmm0 movhpd 3 * SIZE(AA), %xmm0 mulpd %xmm0, %xmm3 movsd 2 * SIZE(AA), %xmm0 movhpd 2 * SIZE(AA), %xmm0 mulpd %xmm3, %xmm0 subpd %xmm0, %xmm2 movsd 0 * SIZE(AA), %xmm0 movhpd 0 * SIZE(AA), %xmm0 mulpd %xmm0, %xmm2#endif#ifdef LT movsd 0 * SIZE(AA), %xmm0 movhpd 0 * SIZE(AA), %xmm0 mulpd %xmm0, %xmm2 movsd 1 * SIZE(AA), %xmm0 movhpd 1 * SIZE(AA), %xmm0 mulpd %xmm2, %xmm0 subpd %xmm0, %xmm3 movsd 3 * SIZE(AA), %xmm0 movhpd 3 * SIZE(AA), %xmm0 mulpd %xmm0, %xmm3#endif#ifdef RN movsd 0 * SIZE(B), %xmm4 movhpd 0 * SIZE(B), %xmm4 mulpd %xmm4, %xmm0 movsd 1 * SIZE(B), %xmm4 movhpd 1 * SIZE(B), %xmm4 mulpd %xmm0, %xmm4 subpd %xmm4, %xmm1 movsd 3 * SIZE(B), %xmm4 movhpd 3 * SIZE(B), %xmm4 mulpd %xmm4, %xmm1#endif#ifdef RT movsd 3 * SIZE(B), %xmm4 movhpd 3 * SIZE(B), %xmm4 mulpd %xmm4, %xmm1 movsd 2 * SIZE(B), %xmm4 movhpd 2 * SIZE(B), %xmm4 mulpd %xmm1, %xmm4 subpd %xmm4, %xmm0 movsd 0 * SIZE(B), %xmm4 movhpd 0 * SIZE(B), %xmm4 mulpd %xmm4, %xmm0#endif#if defined(LN) || defined(LT) movapd %xmm2, 0 * SIZE(B) movapd %xmm3, 2 * SIZE(B) movsd %xmm2, 0 * SIZE(BB) movsd %xmm2, 1 * SIZE(BB) movhpd %xmm2, 2 * SIZE(BB) movhpd %xmm2, 3 * SIZE(BB) movsd %xmm3, 4 * SIZE(BB) movsd %xmm3, 5 * SIZE(BB) movhpd %xmm3, 6 * SIZE(BB) movhpd %xmm3, 7 * SIZE(BB)#else movapd %xmm0, 0 * SIZE(AA) movapd %xmm1, 2 * SIZE(AA)#endif#ifdef LN subl $2 * SIZE, %esi#endif#if defined(LN) || defined(LT) movsd %xmm2, 0 * SIZE(%esi) movsd %xmm3, 1 * SIZE(%esi) movhpd %xmm2, 0 * SIZE(%esi, LDC) movhpd %xmm3, 1 * SIZE(%esi, LDC)#else movsd %xmm0, 0 * SIZE(%esi) movhpd %xmm0, 1 * SIZE(%esi) movsd %xmm1, 0 * SIZE(%esi, LDC) movhpd %xmm1, 1 * SIZE(%esi, LDC)#endif#ifndef LN addl $2 * SIZE, %esi#endif#if defined(LT) || defined(RN) movl K, %eax subl KK, %eax leal (,%eax, SIZE), %eax leal (AA, %eax, 2), AA#ifdef LT addl $4 * SIZE, B#endif#endif#ifdef LN subl $2, KK movl BORIG, B#endif#ifdef LT addl $2, KK#endif#ifdef RT movl K, %eax movl BORIG, B sall $1 + BASE_SHIFT, %eax addl %eax, AORIG#endif ALIGN_2.L50: movl M, %ebx sarl $2, %ebx # i = (m >> 2) jle .L99 ALIGN_4.L10:#ifdef LN movl K, %eax sall $2 + BASE_SHIFT, %eax subl %eax, AORIG#endif#if defined(LN) || defined(RT) movl KK, %eax movl AORIG, AA leal (, %eax, SIZE), %eax leal (AA, %eax, 4), AA#endif leal BUFFER, BB#if defined(LN) || defined(RT) movl KK, %eax sall $1 + BASE_SHIFT, %eax leal (BB, %eax, 2), BB#endif movapd 0 * SIZE(BB), %xmm2 pxor %xmm4, %xmm4 movapd 0 * SIZE(AA), %xmm0 pxor %xmm5, %xmm5 movapd 8 * SIZE(BB), %xmm3 pxor %xmm6, %xmm6 movapd 8 * SIZE(AA), %xmm1 pxor %xmm7, %xmm7#ifdef LN prefetcht2 -4 * SIZE(%esi) prefetcht2 -4 * SIZE(%esi, LDC)#else prefetcht2 4 * SIZE(%esi) prefetcht2 4 * SIZE(%esi, LDC)#endif #if defined(LT) || defined(RN) movl KK, %eax#else movl K, %eax subl KK, %eax#endif#ifdef PENTIUM4 andl $-8, %eax je,pn .L12 sall $3, %eax .L1X: KERNEL1(32 * 0) KERNEL2(32 * 0) KERNEL3(32 * 0) KERNEL4(32 * 0) KERNEL5(32 * 0) KERNEL6(32 * 0) KERNEL7(32 * 0) KERNEL8(32 * 0) cmpl $64 * 1, %eax jle,pn .L11 KERNEL1(32 * 1) KERNEL2(32 * 1) KERNEL3(32 * 1) KERNEL4(32 * 1) KERNEL5(32 * 1) KERNEL6(32 * 1) KERNEL7(32 * 1) KERNEL8(32 * 1) cmpl $64 * 2, %eax jle,pn .L11 KERNEL1(32 * 2) KERNEL2(32 * 2) KERNEL3(32 * 2) KERNEL4(32 * 2) KERNEL5(32 * 2) KERNEL6(32 * 2) KERNEL7(32 * 2) KERNEL8(32 * 2) cmpl $64 * 3, %eax jle,pn .L11 KERNEL1(32 * 3) KERNEL2(32 * 3) KERNEL3(32 * 3) KERNEL4(32 * 3) KERNEL5(32 * 3) KERNEL6(32 * 3) KERNEL7(32 * 3) KERNEL8(32 * 3) cmpl $64 * 4, %eax jle,pn .L11 KERNEL1(32 * 4) KERNEL2(32 * 4) KERNEL3(32 * 4) KERNEL4(32 * 4) KERNEL5(32 * 4) KERNEL6(32 * 4) KERNEL7(32 * 4) KERNEL8(32 * 4) cmpl $64 * 5, %eax jle,pn .L11 KERNEL1(32 * 5) KERNEL2(32 * 5) KERNEL3(32 * 5) KERNEL4(32 * 5) KERNEL5(32 * 5) KERNEL6(32 * 5) KERNEL7(32 * 5) KERNEL8(32 * 5) cmpl $64 * 6, %eax jle,pn .L11 KERNEL1(32 * 6) KERNEL2(32 * 6) KERNEL3(32 * 6) KERNEL4(32 * 6) KERNEL5(32 * 6) KERNEL6(32 * 6) KERNEL7(32 * 6) KERNEL8(32 * 6) cmpl $64 * 7, %eax jle,pn .L11 KERNEL1(32 * 7) KERNEL2(32 * 7) KERNEL3(32 * 7) KERNEL4(32 * 7) KERNEL5(32 * 7) KERNEL6(32 * 7) KERNEL7(32 * 7) KERNEL8(32 * 7) addl $64 * 4 * SIZE, AA addl $64 * 4 * SIZE, BB subl $64 * 8, %eax jg,pt .L1X .L11: leal (AA, %eax, 4), AA leal (BB, %eax, 4), BB#else sarl $3, %eax je .L12.L11: KERNEL1(32 * 0) KERNEL2(32 * 0) KERNEL3(32 * 0) KERNEL4(32 * 0) KERNEL5(32 * 0) KERNEL6(32 * 0) KERNEL7(32 * 0) KERNEL8(32 * 0) addl $32 * SIZE, %ecx addl $32 * SIZE, %edx decl %eax jne .L11#endif.L12:#if defined(LT) || defined(RN) movl KK, %eax#else movl K, %eax subl KK, %eax#endif andl $7, %eax # if (k & 1) BRANCH je .L14.L13: mulpd %xmm0, %xmm2 mulpd 2 * SIZE(BB), %xmm0 addpd %xmm2, %xmm4 movapd 0 * SIZE(BB), %xmm2 addpd %xmm0, %xmm5 movapd 2 * SIZE(AA), %xmm0 mulpd %xmm0, %xmm2 mulpd 2 * SIZE(BB), %xmm0 addpd %xmm2, %xmm6 movapd 4 * SIZE(BB), %xmm2 addpd %xmm0, %xmm7 movapd 4 * SIZE(AA), %xmm0 addl $4 * SIZE, AA # aoffset += 8 addl $4 * SIZE, BB # boffset1 += 8 subl $1, %eax jg .L13 ALIGN_4 .L14:#if defined(LN) || defined(RT) movl KK, %eax#ifdef LN subl $4, %eax#else subl $2, %eax#endif movl AORIG, AA movl BORIG, B leal BUFFER, BB leal (, %eax, SIZE), %eax leal (AA, %eax, 4), AA leal (B, %eax, 2), B leal (BB, %eax, 4), BB#endif#if defined(LN) || defined(LT) movapd %xmm4, %xmm0 unpcklpd %xmm5, %xmm4 unpckhpd %xmm5, %xmm0 movapd %xmm6, %xmm1 unpcklpd %xmm7, %xmm6 unpckhpd %xmm7, %xmm1 movapd 0 * SIZE(B), %xmm2 movapd 2 * SIZE(B), %xmm3 movapd 4 * SIZE(B), %xmm5 movapd 6 * SIZE(B), %xmm7 subpd %xmm4, %xmm2 subpd %xmm0, %xmm3 subpd %xmm6, %xmm5 subpd %xmm1, %xmm7#else movapd 0 * SIZE(AA), %xmm0 movapd 2 * SIZE(AA), %xmm1 movapd 4 * SIZE(AA), %xmm2 movapd 6 * SIZE(AA), %xmm3 subpd %xmm4, %xmm0 subpd %xmm6, %xmm1 subpd %xmm5, %xmm2 subpd %xmm7, %xmm3#endif#ifdef LN movsd 15 * SIZE(AA), %xmm0 movhpd 15 * SIZE(AA), %xmm0 mulpd %xmm0, %xmm7 movsd 14 * SIZE(AA), %xmm0 movhpd 14 * SIZE(AA), %xmm0 mulpd %xmm7, %xmm0 subpd %xmm0, %xmm5 movsd 13 * SIZE(AA), %xmm0 movhpd 13 * SIZE(AA), %xmm0 mulpd %xmm7, %xmm0 subpd %xmm0, %xmm3 movsd 12 * SIZE(AA), %xmm0 movhpd 12 * SIZE(AA), %xmm0 mulpd %xmm7, %xmm0 subpd %xmm0, %xmm2 movsd 10 * SIZE(AA), %xmm0 movhpd 10 * SIZE(AA), %xmm0 mulpd %xmm0, %xmm5 movsd 9 * SIZE(AA), %xmm0 movhpd 9 * SIZE(AA), %xmm0 mulpd %xmm5, %xmm0 subpd %xmm0, %xmm3 movsd 8 * SIZE(AA), %xmm0 movhpd 8 * SIZE(AA), %xmm0 mulpd %xmm5, %xmm0 subpd %xmm0, %xmm2 movsd 5 * SIZE(AA), %xmm0 movhpd 5 * SIZE(AA), %xmm0 mulpd %xmm0, %xmm3 movsd 4 * SIZE(AA), %xmm0 movhpd 4 * SIZE(AA), %xmm0 mulpd %xmm3, %xmm0 subpd %xmm0, %xmm2 movsd 0 * SIZE(AA), %xmm0 movhpd 0 * SIZE(AA), %xmm0 mulpd %xmm0, %xmm2#endif#ifdef LT movsd 0 * SIZE(AA), %xmm0 movhpd 0 * SIZE(AA), %xmm0 mulpd %xmm0, %xmm2 movsd 1 * SIZE(AA), %xmm0 movhpd 1 * SIZE(AA), %xmm0 mulpd %xmm2, %xmm0 subpd %xmm0, %xmm3 movsd 2 * SIZE(AA), %xmm0 movhpd 2 * SIZE(AA), %xmm0 mulpd %xmm2, %xmm0 subpd %xmm0, %xmm5 movsd 3 * SIZE(AA), %xmm0 movhpd 3 * SIZE(AA), %xmm0 mulpd %xmm2, %xmm0 subpd %xmm0, %xmm7 movsd 5 * SIZE(AA), %xmm0 movhpd 5 * SIZE(AA), %xmm0 mulpd %xmm0, %xmm3 movsd 6 * SIZE(AA), %xmm0 movhpd 6 * SIZE(AA), %xmm0 mulpd %xmm3, %xmm0 subpd %xmm0, %xmm5 movsd 7 * SIZE(AA), %xmm0 movhpd 7 * SIZE(AA), %xmm0 mulpd %xmm3, %xmm0 subpd %xmm0, %xmm7 movsd 10 * SIZE(AA), %xmm0 movhpd 10 * SIZE(AA), %xmm0 mulpd %xmm0, %xmm5 movsd 11 * SIZE(AA), %xmm0 movhpd 11 * SIZE(AA), %xmm0 mulpd %xmm5, %xmm0 subpd %xmm0, %xmm7 movsd 15 * SIZE(AA), %xmm0 movhpd 15 * SIZE(AA), %xmm0 mulpd %xmm0, %xmm7#endif#ifdef RN movsd 0 * SIZE(B), %xmm4 movhpd 0 * SIZE(B), %xmm4 mulpd %xmm4, %xmm0 mulpd %xmm4, %xmm1 movsd 1 * SIZE(B), %xmm4 movhpd 1 * SIZE(B), %xmm4 mulpd %xmm0, %xmm4 subpd %xmm4, %xmm2 movsd 1 * SIZE(B), %xmm4 movhpd 1 * SIZE(B), %xmm4 mulpd %xmm1, %xmm4 subpd %xmm4, %xmm3 movsd 3 * SIZE(B), %xmm4 movhpd 3 * SIZE(B), %xmm4 mulpd %xmm4, %xmm2 mulpd %xmm4, %xmm3#endif#ifdef RT movsd 3 * SIZE(B), %xmm4 movhpd 3 * SIZE(B), %xmm4 mulpd %xmm4, %xmm2 mulpd %xmm4, %xmm3 movsd 2 * SIZE(B), %xmm4 movhpd 2 * SIZE(B), %xmm4 mulpd %xmm2, %xmm4 subpd %xmm4, %xmm0 movsd 2 * SIZE(B), %xmm4 movhpd 2 * SIZE(B), %xmm4 mulpd %xmm3, %xmm4 subpd %xmm4, %xmm1 movsd 0 * SIZE(B), %xmm4 movhpd 0 * SIZE(B), %xmm4 mulpd %xmm4, %xmm0 mulpd %xmm4, %xmm1#endif#if defined(LN) || defined(LT) movapd %xmm2, 0 * SIZE(B) movapd %xmm3, 2 * SIZE(B) movapd %xmm5, 4 * SIZE(B) movapd %xmm7, 6 * SIZE(B) movsd %xmm2, 0 * SIZE(BB) movsd %xmm2, 1 * SIZE(BB) movhpd %xmm2, 2 * SIZE(BB) movhpd %xmm2, 3 * SIZE(BB) movsd %xmm3, 4 * SIZE(BB) movsd %xmm3, 5 * SIZE(BB) movhpd %xmm3, 6 * SIZE(BB) movhpd %xmm3, 7 * SIZE(BB) movsd %xmm5, 8 * SIZE(BB) movsd %xmm5, 9 * SIZE(BB) movhpd %xmm5, 10 * SIZE(BB) movhpd %xmm5, 11 * SIZE(BB) movsd %xmm7, 12 * SIZE(BB) movsd %xmm7, 13 * SIZE(BB) movhpd %xmm7, 14 * SIZE(BB) movhpd %xmm7, 15 * SIZE(BB)#else movapd %xmm0, 0 * SIZE(AA) movapd %xmm1, 2 * SIZE(AA) movapd %xmm2, 4 * SIZE(AA) movapd %xmm3, 6 * SIZE(AA)#endif#ifdef LN subl $4 * SIZE, %esi#endif#if defined(LN) || defined(LT) movsd %xmm2, 0 * SIZE(%esi) movsd %xmm3, 1 * SIZE(%esi) movsd %xmm5, 2 * SIZE(%esi) movsd %xmm7, 3 * SIZE(%esi) movhpd %xmm2, 0 * SIZE(%esi, LDC) movhpd %xmm3, 1 * SIZE(%esi, LDC) movhpd %xmm5, 2 * SIZE(%esi, LDC) movhpd %xmm7, 3 * SIZE(%esi, LDC)#else movsd %xmm0, 0 * SIZE(%esi) movhpd %xmm0, 1 * SIZE(%esi) movsd %xmm1, 2 * SIZE(%esi) movhpd %xmm1, 3 * SIZE(%esi) movsd %xmm2, 0 * SIZE(%esi, LDC) movhpd %xmm2, 1 * SIZE(%esi, LDC) movsd %xmm3, 2 * SIZE(%esi, LDC) movhpd %xmm3, 3 * SIZE(%esi, LDC)#endif#ifndef LN addl $4 * SIZE, %esi#endif#if defined(LT) || defined(RN) movl K, %eax subl KK, %eax leal (,%eax, SIZE), %eax leal (AA, %eax, 4), AA#ifdef LT addl $8 * SIZE, B#endif#endif#ifdef LN subl $4, KK movl BORIG, B#endif#ifdef LT addl $4, KK#endif#ifdef RT movl K, %eax movl BORIG, B sall $2 + BASE_SHIFT, %eax addl %eax, AORIG#endif decl %ebx # i -- jg .L10 ALIGN_2 .L99:#ifdef LN movl K, %eax leal (, %eax, SIZE), %eax leal (B, %eax, 2), B#endif#if defined(LT) || defined(RN) movl K, %eax subl KK, %eax leal (,%eax, SIZE), %eax leal (B, %eax, 2), B#endif#ifdef RN addl $2, KK#endif#ifdef RT subl $2, KK#endif decl J # j -- jg .L01 ALIGN_2.L100: movl N, %eax testl $1, %eax jle .L999 ALIGN_2 .L101:/* Copying to Sub Buffer */#ifdef LN movl OFFSET, %eax addl M, %eax movl %eax, KK#endif leal BUFFER, %ecx#ifdef RT movl K, %eax sall $0 + BASE_SHIFT, %eax subl %eax, B#endif#if defined(LN) || defined(RT) movl KK, %eax movl B, BORIG leal (, %eax, SIZE), %eax leal (B, %eax, 1), B leal (BB, %eax, 2), BB#endif #ifdef LT movl OFFSET, %eax movl %eax, KK#endif#if defined(LT) || defined(RN) movl KK, %eax#else movl K, %eax subl KK, %eax#endif sarl $3, %eax jle .L103 ALIGN_4 .L102: movsd 0 * SIZE(B), %xmm0 movsd 1 * SIZE(B), %xmm1 movsd 2 * SIZE(B), %xmm2 movsd 3 * SIZE(B), %xmm3 movsd 4 * SIZE(B), %xmm4 movsd 5 * SIZE(B), %xmm5 movsd 6 * SIZE(B), %xmm6 movsd 7 * SIZE(B), %xmm7 unpcklpd %xmm0, %xmm0 unpcklpd %xmm1, %xmm1 unpcklpd %xmm2, %xmm2 unpcklpd %xmm3, %xmm3 unpcklpd %xmm4, %xmm4 unpcklpd %xmm5, %xmm5 unpcklpd %xmm6, %xmm6 unpcklpd %xmm7, %xmm7 movapd %xmm0, 0 * SIZE(%ecx) movapd %xmm1, 2 * SIZE(%ecx) movapd %xmm2, 4 * SIZE(%ecx) movapd %xmm3, 6 * SIZE(%ecx) movapd %xmm4, 8 * SIZE(%ecx) movapd %xmm5, 10 * SIZE(%ecx) movapd %xmm6, 12 * SIZE(%ecx) movapd %xmm7, 14 * SIZE(%ecx) prefetcht0 104 * SIZE(B) addl $ 8 * SIZE, B addl $16 * SIZE, %ecx decl %eax BRANCH jne .L102 ALIGN_2.L103:#if defined(LT) || defined(RN) movl KK, %eax#else movl K, %eax subl KK, %eax#endif andl $7, %eax BRANCH jle .L105 ALIGN_2.L104: movsd 0 * SIZE(B), %xmm0 unpcklpd %xmm0, %xmm0 movapd %xmm0, 0 * SIZE(%ecx) addl $1 * SIZE, B addl $2 * SIZE, %ecx decl %eax jne .L104 ALIGN_4 .L105:#if defined(LT) || defined(RN) movl A, AA#else movl A, %eax movl %eax, AORIG#endif#ifdef RT subl LDC, C#endif movl C, %esi # coffset = c#ifndef RT addl LDC, C#endif movl M, %ebx testl $1, %ebx jle .L130
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -