atl_dmm1x6x72_sse2.c
来自「基于Blas CLapck的.用过的人知道是干啥的」· C语言 代码 · 共 1,724 行 · 第 1/3 页
C
1,724 行
addpd rA0, rC4 mulpd 144-120(pA3,ldab,2), rB0 addpd rB0, rC5 movapd 160-120(pB0), rB0 movapd 160-120(pA0), rA0 mulpd rB0, rA0 addpd rA0, rC0 movapd 160-120(pA0,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC1 movapd 160-120(pA0,ldab,2), rA0 mulpd rB0, rA0 addpd rA0, rC2 movapd 160-120(pA3), rA0 mulpd rB0, rA0 addpd rA0, rC3 movapd 160-120(pA3,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC4 mulpd 160-120(pA3,ldab,2), rB0 addpd rB0, rC5 movapd 176-120(pB0), rB0 movapd 176-120(pA0), rA0 mulpd rB0, rA0 addpd rA0, rC0 movapd 176-120(pA0,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC1 movapd 176-120(pA0,ldab,2), rA0 mulpd rB0, rA0 addpd rA0, rC2 movapd 176-120(pA3), rA0 mulpd rB0, rA0 addpd rA0, rC3 movapd 176-120(pA3,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC4 mulpd 176-120(pA3,ldab,2), rB0 addpd rB0, rC5 movapd 192-120(pB0), rB0 movapd 192-120(pA0), rA0 mulpd rB0, rA0 addpd rA0, rC0 movapd 192-120(pA0,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC1 movapd 192-120(pA0,ldab,2), rA0 mulpd rB0, rA0 addpd rA0, rC2 movapd 192-120(pA3), rA0 mulpd rB0, rA0 addpd rA0, rC3 movapd 192-120(pA3,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC4 mulpd 192-120(pA3,ldab,2), rB0 addpd rB0, rC5 movapd 208-120(pB0), rB0 movapd 208-120(pA0), rA0 mulpd rB0, rA0 addpd rA0, rC0 movapd 208-120(pA0,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC1 movapd 208-120(pA0,ldab,2), rA0 mulpd rB0, rA0 addpd rA0, rC2 movapd 208-120(pA3), rA0 mulpd rB0, rA0 addpd rA0, rC3 movapd 208-120(pA3,ldab), rA0 mulpd rB0, rA0 prefB(448(pB0,ldab)) addpd rA0, rC4 mulpd 208-120(pA3,ldab,2), rB0 addpd rB0, rC5 movapd 224-120(pB0), rB0 movapd 224-120(pA0), rA0 mulpd rB0, rA0 addpd rA0, rC0 movapd 224-120(pA0,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC1 movapd 224-120(pA0,ldab,2), rA0 mulpd rB0, rA0 addpd rA0, rC2 movapd 224-120(pA3), rA0 mulpd rB0, rA0 addpd rA0, rC3 movapd 224-120(pA3,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC4 mulpd 224-120(pA3,ldab,2), rB0 addpd rB0, rC5 movapd 240-120(pB0), rB0 movapd 240-120(pA0), rA0 mulpd rB0, rA0 addpd rA0, rC0 movapd 240-120(pA0,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC1 movapd 240-120(pA0,ldab,2), rA0 mulpd rB0, rA0 addpd rA0, rC2 movapd 240-120(pA3), rA0 mulpd rB0, rA0 addpd rA0, rC3 movapd 240-120(pA3,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC4 mulpd 240-120(pA3,ldab,2), rB0 addpd rB0, rC5 movapd 256-120(pB0), rB0 movapd 256-120(pA0), rA0 mulpd rB0, rA0 addpd rA0, rC0 movapd 256-120(pA0,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC1 movapd 256-120(pA0,ldab,2), rA0 mulpd rB0, rA0 addpd rA0, rC2 movapd 256-120(pA3), rA0 mulpd rB0, rA0 addpd rA0, rC3 movapd 256-120(pA3,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC4 mulpd 256-120(pA3,ldab,2), rB0 addpd rB0, rC5 movapd 272-120(pB0), rB0 movapd 272-120(pA0), rA0 mulpd rB0, rA0 addpd rA0, rC0 movapd 272-120(pA0,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC1 movapd 272-120(pA0,ldab,2), rA0 mulpd rB0, rA0 addpd rA0, rC2 movapd 272-120(pA3), rA0 mulpd rB0, rA0 addpd rA0, rC3 movapd 272-120(pA3,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC4 mulpd 272-120(pA3,ldab,2), rB0 addpd rB0, rC5 movapd 288-120(pB0), rB0 movapd 288-120(pA0), rA0 mulpd rB0, rA0 addpd rA0, rC0 movapd 288-120(pA0,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC1 movapd 288-120(pA0,ldab,2), rA0 mulpd rB0, rA0 addpd rA0, rC2 movapd 288-120(pA3), rA0 mulpd rB0, rA0 addpd rA0, rC3 movapd 288-120(pA3,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC4 mulpd 288-120(pA3,ldab,2), rB0 addpd rB0, rC5 movapd 304-120(pB0), rB0 movapd 304-120(pA0), rA0 mulpd rB0, rA0 addpd rA0, rC0 movapd 304-120(pA0,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC1 movapd 304-120(pA0,ldab,2), rA0 mulpd rB0, rA0 addpd rA0, rC2 movapd 304-120(pA3), rA0 mulpd rB0, rA0 addpd rA0, rC3 movapd 304-120(pA3,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC4 mulpd 304-120(pA3,ldab,2), rB0 addpd rB0, rC5 movapd 320-120(pB0), rB0 movapd 320-120(pA0), rA0 mulpd rB0, rA0 addpd rA0, rC0 movapd 320-120(pA0,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC1 movapd 320-120(pA0,ldab,2), rA0 mulpd rB0, rA0 addpd rA0, rC2 movapd 320-120(pA3), rA0 mulpd rB0, rA0 addpd rA0, rC3 movapd 320-120(pA3,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC4 mulpd 320-120(pA3,ldab,2), rB0 addpd rB0, rC5 movapd 336-120(pB0), rB0 movapd 336-120(pA0), rA0 mulpd rB0, rA0 addpd rA0, rC0 movapd 336-120(pA0,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC1 movapd 336-120(pA0,ldab,2), rA0 mulpd rB0, rA0 addpd rA0, rC2 movapd 336-120(pA3), rA0 mulpd rB0, rA0 addpd rA0, rC3 movapd 336-120(pA3,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC4 mulpd 336-120(pA3,ldab,2), rB0 addpd rB0, rC5 movapd 352-120(pB0), rB0 movapd 352-120(pA0), rA0 mulpd rB0, rA0 addpd rA0, rC0 movapd 352-120(pA0,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC1 movapd 352-120(pA0,ldab,2), rA0 mulpd rB0, rA0 addpd rA0, rC2 movapd 352-120(pA3), rA0 mulpd rB0, rA0 addpd rA0, rC3 movapd 352-120(pA3,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC4 mulpd 352-120(pA3,ldab,2), rB0 addpd rB0, rC5 movapd 368-120(pB0), rB0 movapd 368-120(pA0), rA0 mulpd rB0, rA0 addpd rA0, rC0 movapd 368-120(pA0,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC1 movapd 368-120(pA0,ldab,2), rA0 mulpd rB0, rA0 addpd rA0, rC2 movapd 368-120(pA3), rA0 mulpd rB0, rA0 addpd rA0, rC3 movapd 368-120(pA3,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC4 mulpd 368-120(pA3,ldab,2), rB0 addpd rB0, rC5 movapd 384-120(pB0), rB0 movapd 384-120(pA0), rA0 mulpd rB0, rA0 addpd rA0, rC0 movapd 384-120(pA0,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC1 movapd 384-120(pA0,ldab,2), rA0 mulpd rB0, rA0 addpd rA0, rC2 movapd 384-120(pA3), rA0 mulpd rB0, rA0 addpd rA0, rC3 movapd 384-120(pA3,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC4 mulpd 384-120(pA3,ldab,2), rB0 addpd rB0, rC5 movapd 400-120(pB0), rB0 movapd 400-120(pA0), rA0 mulpd rB0, rA0 addpd rA0, rC0 movapd 400-120(pA0,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC1 movapd 400-120(pA0,ldab,2), rA0 mulpd rB0, rA0 addpd rA0, rC2 movapd 400-120(pA3), rA0 mulpd rB0, rA0 addpd rA0, rC3 movapd 400-120(pA3,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC4 mulpd 400-120(pA3,ldab,2), rB0 addpd rB0, rC5 movapd 416-120(pB0), rB0 movapd 416-120(pA0), rA0 mulpd rB0, rA0 addpd rA0, rC0 movapd 416-120(pA0,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC1 movapd 416-120(pA0,ldab,2), rA0 mulpd rB0, rA0 addpd rA0, rC2 movapd 416-120(pA3), rA0 mulpd rB0, rA0 addpd rA0, rC3 movapd 416-120(pA3,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC4 mulpd 416-120(pA3,ldab,2), rB0 addpd rB0, rC5 movapd 432-120(pB0), rB0 movapd 432-120(pA0), rA0 mulpd rB0, rA0 addpd rA0, rC0 movapd 432-120(pA0,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC1 movapd 432-120(pA0,ldab,2), rA0 mulpd rB0, rA0 addpd rA0, rC2 movapd 432-120(pA3), rA0 mulpd rB0, rA0 addpd rA0, rC3 movapd 432-120(pA3,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC4 mulpd 432-120(pA3,ldab,2), rB0 addpd rB0, rC5 movapd 448-120(pB0), rB0 movapd 448-120(pA0), rA0 mulpd rB0, rA0 addpd rA0, rC0 movapd 448-120(pA0,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC1 movapd 448-120(pA0,ldab,2), rA0 mulpd rB0, rA0 addpd rA0, rC2 movapd 448-120(pA3), rA0 mulpd rB0, rA0 addpd rA0, rC3 movapd 448-120(pA3,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC4 mulpd 448-120(pA3,ldab,2), rB0 addpd rB0, rC5 movapd 464-120(pB0), rB0 movapd 464-120(pA0), rA0 mulpd rB0, rA0 addpd rA0, rC0 movapd 464-120(pA0,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC1 movapd 464-120(pA0,ldab,2), rA0 mulpd rB0, rA0 addpd rA0, rC2 movapd 464-120(pA3), rA0 mulpd rB0, rA0 addpd rA0, rC3 movapd 464-120(pA3,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC4 mulpd 464-120(pA3,ldab,2), rB0 addpd rB0, rC5 movapd 480-120(pB0), rB0 movapd 480-120(pA0), rA0 mulpd rB0, rA0 addpd rA0, rC0 movapd 480-120(pA0,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC1 movapd 480-120(pA0,ldab,2), rA0 mulpd rB0, rA0 addpd rA0, rC2 movapd 480-120(pA3), rA0 mulpd rB0, rA0 addpd rA0, rC3 movapd 480-120(pA3,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC4 mulpd 480-120(pA3,ldab,2), rB0 addpd rB0, rC5 movapd 496-120(pB0), rB0 movapd 496-120(pA0), rA0 mulpd rB0, rA0 addpd rA0, rC0 movapd 496-120(pA0,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC1 movapd 496-120(pA0,ldab,2), rA0 mulpd rB0, rA0 addpd rA0, rC2 movapd 496-120(pA3), rA0 mulpd rB0, rA0 addpd rA0, rC3 movapd 496-120(pA3,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC4 mulpd 496-120(pA3,ldab,2), rB0 addpd rB0, rC5 movapd 512-120(pB0), rB0 movapd 512-120(pA0), rA0 mulpd rB0, rA0 addpd rA0, rC0 movapd 512-120(pA0,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC1 movapd 512-120(pA0,ldab,2), rA0 mulpd rB0, rA0 addpd rA0, rC2 movapd 512-120(pA3), rA0 mulpd rB0, rA0 addpd rA0, rC3 movapd 512-120(pA3,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC4 mulpd 512-120(pA3,ldab,2), rB0 addpd rB0, rC5 movapd 528-120(pB0), rB0 movapd 528-120(pA0), rA0 mulpd rB0, rA0 addpd rA0, rC0 movapd 528-120(pA0,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC1 movapd 528-120(pA0,ldab,2), rA0 mulpd rB0, rA0 addpd rA0, rC2 movapd 528-120(pA3), rA0 mulpd rB0, rA0 addpd rA0, rC3 movapd 528-120(pA3,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC4 mulpd 528-120(pA3,ldab,2), rB0 addpd rB0, rC5/* pref2((pfA)) *//* addl $48, pfA */ movapd 544-120(pB0), rB0 movapd 544-120(pA0), rA0 mulpd rB0, rA0 addpd rA0, rC0 movapd 544-120(pA0,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC1 movapd 544-120(pA0,ldab,2), rA0 mulpd rB0, rA0 addpd rA0, rC2 movapd 544-120(pA3), rA0 mulpd rB0, rA0 addpd rA0, rC3 movapd 544-120(pA3,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC4 mulpd 544-120(pA3,ldab,2), rB0 addpd rB0, rC5 movapd 560-120(pB0), rB0 movapd 560-120(pA0), rA0 mulpd rB0, rA0 addpd rA0, rC0 movapd 560-120(pA0,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC1 movapd 560-120(pA0,ldab,2), rA0 mulpd rB0, rA0 addpd rA0, rC2 movapd 560-120(pA3), rA0 mulpd rB0, rA0 addpd rA0, rC3 movapd 560-120(pA3,ldab), rA0 mulpd rB0, rA0 addpd rA0, rC4 mulpd 560-120(pA3,ldab,2), rB0 addpd rB0, rC5/* * Get these bastard things summed up correctly */ /* rC0 = c0a c0b */ /* rC1 = c1a c1b */ /* rC2 = c2a c2b */ /* rC3 = c3a c3b *//* */ /* rC4 = c4a c4b */ /* rC5 = c5a c5b */ movapd rC0, rA0 addl $NBso, pB0 movapd rC4, rB0 unpcklpd rC1, rC0 /* rC0 = c0a c1a */ subl $NBNBso-NB6so, pA0 unpcklpd rC5, rC4 /* rC4 = c4a c5a */ unpckhpd rC1, rA0 /* rA0 = c0b c1b */ unpckhpd rC5, rB0 /* rB0 = c4b c5b */ addpd rA0, rC0 /* rC0 = c0ab c1ab */ addpd rB0, rC4 /* rC4 = c4ab c5ab */ movapd rC2, rA0 unpcklpd rC3, rC2 /* rC2 = c2a c3a */ unpckhpd rC3, rA0 /* rA0 = c2b c3b */ subl $NBNBso-NB6so, pA3 addpd rA0, rC2 /* rC2 = c2ab c3ab *//* * Write results back to C */ movl COFF(%esp), ldc movlpd rC0, (pC) movhpd rC0, (pC,ldc) movlpd rC2, (pC,ldc,2) movlpd rC4, (pC,ldc,4) addl ldc, pC movhpd rC2, (pC,ldc,2) movhpd rC4, (pC,ldc,4)/* * pC += 6; pA += 6*NB *//* * pC += ldc; pA -= NBNB; pB += NB; */ movl COFF+4(%esp), pC#ifdef DCPLX addl $16, pC#else addl $8, pC#endif movl pC, COFF+4(%esp)/* * while (pB != stN); */ sub $1, stN jnz NLOOP/* * Restore callee-saved iregs */ #ifdef BETAX movl 32(%esp), %ebp movl 28(%esp), %ebx movl 24(%esp), %esi movl 20(%esp), %edi movl 16(%esp), %esp #else movl 12(%esp), %ebp movl 8(%esp), %ebx movl 4(%esp), %esi movl (%esp), %edi addl $24, %esp #endif ret
⌨️ 快捷键说明
复制代码Ctrl + C
搜索代码Ctrl + F
全屏模式F11
增大字号Ctrl + =
减小字号Ctrl + -
显示快捷键?