atl_smm14x1x84_sse.c
来自「基于Blas CLapck的.用过的人知道是干啥的」· C语言 代码 · 共 2,566 行 · 第 1/5 页
C
2,566 行
mulps rB0, rA0 addps rA0, rC12 mulps 64-120(pA5,ldab,8), rB0 addps rB0, rC13 prefB(128-176(pB,ldab)) prefB(192-176(pB,ldab)) movaps 80-120(pA10,mldab5,2), rA0 movaps 80-120(pB0), rB0 mulps rB0, rA0 addps rA0, rC0 movaps 80-120(pA5, mldab,4), rA0 mulps rB0, rA0 addps rA0, rC1 movaps 80-120(pA10, mldab,8), rA0 mulps rB0, rA0 addps rA0, rC2 movaps 80-120(pA5, mldab,2), rA0 mulps rB0, rA0 addps rA0, rC3 movaps 80-120(pA5, mldab), rA0 mulps rB0, rA0 addps rA0, rC4 movaps 80-120(pA5), rA0 mulps rB0, rA0 addps rA0, rC5 movaps 80-120(pA5, ldab), rA0 mulps rB0, rA0 addps rA0, rC6 movaps 80-120(pA5, ldab,2), rA0 mulps rB0, rA0 addps rA0, rC7 movaps 80-120(pA10, mldab,2), rA0 mulps rB0, rA0 addps rA0, rC8 movaps 80-120(pA5,ldab,4), rA0 mulps rB0, rA0 addps rA0, rC9 movaps 80-120(pA10), rA0 mulps rB0, rA0 addps rA0, rC10 movaps 80-120(pA10,ldab), rA0 mulps rB0, rA0 addps rA0, rC11 movaps 80-120(pA10,ldab,2), rA0 mulps rB0, rA0 addps rA0, rC12 mulps 80-120(pA5,ldab,8), rB0 addps rB0, rC13 movaps 96-120(pA10,mldab5,2), rA0 movaps 96-120(pB0), rB0 mulps rB0, rA0 addps rA0, rC0 movaps 96-120(pA5, mldab,4), rA0 mulps rB0, rA0 addps rA0, rC1 movaps 96-120(pA10, mldab,8), rA0 mulps rB0, rA0 addps rA0, rC2 movaps 96-120(pA5, mldab,2), rA0 mulps rB0, rA0 addps rA0, rC3 movaps 96-120(pA5, mldab), rA0 mulps rB0, rA0 addps rA0, rC4 movaps 96-120(pA5), rA0 mulps rB0, rA0 addps rA0, rC5 movaps 96-120(pA5, ldab), rA0 mulps rB0, rA0 addps rA0, rC6 movaps 96-120(pA5, ldab,2), rA0 mulps rB0, rA0 addps rA0, rC7 movaps 96-120(pA10, mldab,2), rA0 mulps rB0, rA0 addps rA0, rC8 movaps 96-120(pA5,ldab,4), rA0 mulps rB0, rA0 addps rA0, rC9 movaps 96-120(pA10), rA0 mulps rB0, rA0 addps rA0, rC10 movaps 96-120(pA10,ldab), rA0 mulps rB0, rA0 addps rA0, rC11 movaps 96-120(pA10,ldab,2), rA0 mulps rB0, rA0 addps rA0, rC12 mulps 96-120(pA5,ldab,8), rB0 addps rB0, rC13 movaps 112-120(pA10,mldab5,2), rA0 movaps 112-120(pB0), rB0 mulps rB0, rA0 addps rA0, rC0 movaps 112-120(pA5, mldab,4), rA0 mulps rB0, rA0 addps rA0, rC1 movaps 112-120(pA10, mldab,8), rA0 mulps rB0, rA0 addps rA0, rC2 movaps 112-120(pA5, mldab,2), rA0 mulps rB0, rA0 addps rA0, rC3 movaps 112-120(pA5, mldab), rA0 mulps rB0, rA0 addps rA0, rC4 movaps 112-120(pA5), rA0 mulps rB0, rA0 addps rA0, rC5 movaps 112-120(pA5, ldab), rA0 mulps rB0, rA0 addps rA0, rC6 movaps 112-120(pA5, ldab,2), rA0 mulps rB0, rA0 addps rA0, rC7 movaps 112-120(pA10, mldab,2), rA0 mulps rB0, rA0 addps rA0, rC8 movaps 112-120(pA5,ldab,4), rA0 mulps rB0, rA0 addps rA0, rC9 movaps 112-120(pA10), rA0 mulps rB0, rA0 addps rA0, rC10 movaps 112-120(pA10,ldab), rA0 mulps rB0, rA0 addps rA0, rC11 movaps 112-120(pA10,ldab,2), rA0 mulps rB0, rA0 addps rA0, rC12 mulps 112-120(pA5,ldab,8), rB0 prefC((pC)) prefC((pC,incCn)) addps rB0, rC13 movaps 128-120(pA10,mldab5,2), rA0 movaps 128-120(pB0), rB0 mulps rB0, rA0 addps rA0, rC0 movaps 128-120(pA5, mldab,4), rA0 mulps rB0, rA0 addps rA0, rC1 movaps 128-120(pA10, mldab,8), rA0 mulps rB0, rA0 addps rA0, rC2 movaps 128-120(pA5, mldab,2), rA0 mulps rB0, rA0 addps rA0, rC3 movaps 128-120(pA5, mldab), rA0 mulps rB0, rA0 addps rA0, rC4 movaps 128-120(pA5), rA0 mulps rB0, rA0 addps rA0, rC5 movaps 128-120(pA5, ldab), rA0 mulps rB0, rA0 addps rA0, rC6 movaps 128-120(pA5, ldab,2), rA0 mulps rB0, rA0 addps rA0, rC7 movaps 128-120(pA10, mldab,2), rA0 mulps rB0, rA0 addps rA0, rC8 movaps 128-120(pA5,ldab,4), rA0 mulps rB0, rA0 addps rA0, rC9 movaps 128-120(pA10), rA0 mulps rB0, rA0 addps rA0, rC10 movaps 128-120(pA10,ldab), rA0 mulps rB0, rA0 addps rA0, rC11 movaps 128-120(pA10,ldab,2), rA0 mulps rB0, rA0 addps rA0, rC12 mulps 128-120(pA5,ldab,8), rB0 addps rB0, rC13 movaps 144-120(pA10,mldab5,2), rA0 movaps 144-120(pB0), rB0 mulps rB0, rA0 addps rA0, rC0 movaps 144-120(pA5, mldab,4), rA0#ifdef BETA1 addss (pC), rC0#endif mulps rB0, rA0 addps rA0, rC1 movaps 144-120(pA10, mldab,8), rA0#ifdef BETA1 addss CMUL(4)(pC), rC1#endif mulps rB0, rA0 addps rA0, rC2 movaps 144-120(pA5, mldab,2), rA0#ifdef BETA1 addss CMUL(8)(pC), rC2#endif mulps rB0, rA0 addps rA0, rC3 movaps 144-120(pA5, mldab), rA0#ifdef BETA1 addss CMUL(12)(pC), rC3#endif mulps rB0, rA0 addps rA0, rC4 movaps 144-120(pA5), rA0#ifdef BETA1 addss CMUL(16)(pC), rC4#endif mulps rB0, rA0 addps rA0, rC5 movaps 144-120(pA5, ldab), rA0#ifdef BETA1 addss CMUL(20)(pC), rC5#endif mulps rB0, rA0 addps rA0, rC6 movaps 144-120(pA5, ldab,2), rA0#ifdef BETA1 addss CMUL(24)(pC), rC6#endif mulps rB0, rA0 addps rA0, rC7 movaps 144-120(pA10, mldab,2), rA0#ifdef BETA1 addss CMUL(28)(pC), rC7#endif mulps rB0, rA0 addps rA0, rC8 movaps 144-120(pA5,ldab,4), rA0#ifdef BETA1 addss CMUL(32)(pC), rC8#endif mulps rB0, rA0 addps rA0, rC9 movaps 144-120(pA10), rA0#ifdef BETA1 addss CMUL(36)(pC), rC9#endif mulps rB0, rA0 addps rA0, rC10 movaps 144-120(pA10,ldab), rA0#ifdef BETA1 addss CMUL(40)(pC), rC10#endif mulps rB0, rA0 addps rA0, rC11 movaps 144-120(pA10,ldab,2), rA0#ifdef BETA1 addss CMUL(44)(pC), rC11#endif mulps rB0, rA0 addps rA0, rC12 mulps 144-120(pA5,ldab,8), rB0 addps rB0, rC13/*UKLOOP *//* * Get these bastard things summed up correctly */ /* rC0 = c0a c0b c0c c0d */ /* rC1 = c1a c1b c1c c1d */ /* rC2 = c2a c2b c2c c2d */ /* rC3 = c3a c3b c3c c3d *//* */ movaps rC2, rB0 /* rB0 = c2a c2b c2c c2d */ prefC(64(pC,incCn)) prefB(256-176(pB,ldab)) movaps rC0, rA0 /* rA0 = c0a c0b c0c c0d */ unpckhps rC3, rB0 /* rB0 = c2c c3c c2d c3d */#ifdef BETA1 addss CMUL(52)(pC), rC13#endif unpckhps rC1, rA0 /* rA0 = c0c c1c c0d c1d */#ifdef BETA1 addss CMUL(48)(pC), rC12#endif unpcklps rC3, rC2 /* rC2 = c2a c3a c2b c3b */ movlhps rB0, rC3 /* rC3 = c3a c3b c2c c3c */ unpcklps rC1, rC0 /* rC0 = c0a c1a c0b c1b */ movhlps rA0, rC3 /* rC3 = c0d c1d c2c c3c */ movlhps rC2, rA0 /* rA0 = c0c c1c c2a c3a */ movhlps rC0, rB0 /* rB0 = c0b c1b c2d c3d */ addps rA0, rC3 /* rC3 = c0cd c1cd c2ac c3ac */ movlhps rC0, rC1 /* rC1 = c1a c1b c0a c1a */ movhlps rC1, rC2 /* rC2 = c0a c1a c2b c3b */ movaps rC4, rA0 /* rA0 = c4a c4b c4c c4d */ addps rB0, rC2 /* rC2 = c0ab c1ab c2bd c3bd */ movaps rC6, rB0 /* rB0 = c6a c6b c6c c6d */ addps rC2, rC3 /* rC3 = c0abcd c1abcd c2bdac c3bdac */ /* rC4 = c4a c4b c4c c4d */ /* rC5 = c5a c5b c5c c5d */ /* rC6 = c6a c6b c6c c6d */ /* rC7 = c7a c7b c7c c7d */ /* rC8 = c08a c08b c08c c08d */ /* rC9 = c09a c09b c09c c09d */ /* rC10 = c10a c10b c10c c10d */ /* rC11 = c11a c11b c11c c11d */ /* rC12 = c12a c12b c12c c12d */ /* rC13 = c13a c13b c13c c13d *//* */ movaps rC10, rC0 /* rC0 = c10a c10b c10c c10d */ movaps rC8 , rC1 /* rC1 = c08a c08b c08c c08d */ movaps rC12, rC2 /* rC2 = c12a c12b c12c c12d */ unpckhps rC7, rB0 /* rB0 = c6c c7c c6d c7d */ unpckhps rC5, rA0 /* rA0 = c4c c5c c4d c5d */ unpcklps rC7, rC6 /* rC6 = c6a c7a c6b c7b */ unpckhps rC11, rC0 /* rC0 = c10c c11c c10d c11d */ unpckhps rC9 , rC1 /* rC1 = c08c c09c c08d c09d */ movlhps rB0, rC7 /* rC7 = c7a c7b c6c c7c */ unpcklps rC5, rC4 /* rC4 = c4a c5a c4b c5b */ movhlps rA0, rC7 /* rC7 = c4d c5d c6c c7c */ movlhps rC6, rA0 /* rA0 = c4c c5c c6a c7a */ unpcklps rC11, rC10 /* rC10 = c10a c11a c10b c11b */ movhlps rC4, rB0 /* rB0 = c4b c5b c6d c7d */ movlhps rC0, rC11 /* rC11 = c11a c11b c10c c11c */ addps rA0, rC7 /* rC7 = c4cd c5cd c6ac c7ac */#ifdef BETAX #ifdef SREAL movups (pC), rA0 movlhps rC4, rC5 /* rC5 = c5a c5b c4a c5a */ movups 16(pC), rC4 unpcklps rC9 , rC8 /* rC8 = c08a c09a c08b c09b */ movhlps rC1, rC11 /* rC11 = c08d c09d c10c c11c */ movlhps rC10, rC1 /* rC1 = c08c c09c c10a c11a */ movhlps rC5, rC6 /* rC6 = c4a c5a c6b c7b */ movups 32(pC), rC5 movhlps rC8 , rC0 /* rC0 = c08b c09b c10d c11d */ unpcklps rC13, rC2 /* rC2 = c12a c13a c12b c13b */ addps rC1, rC11 /* rC11 = c08cd c09cd c10ac c11ac */ movlps 48(pC), rC1 addps rB0, rC6 /* rC6 = c4ab c5ab c6bd c7bd */ movlhps rC8 , rC9 /* rC9 = c09a c09b c08a c09a */ unpckhps rC13, rC12 /* rC12 = c12c c13c c12d c13d */ movhlps rC9 , rC10 /* rC10 = c08a c09a c10b c11b */ addps rC6, rC7 /* rC7 = c4abcd c5abcd c6bdac c7bdac */ mulps BOF(%rsp), rA0 addps rC0, rC10 /* rC10 = c08ab c09ab c10bd c11bd */ mulps BOF(%rsp), rC4 addps rC2, rC12 /* rC12 = c12ac c13ac c12bd c13bd */ mulps BOF(%rsp), rC5 addps rC10, rC11 /* rC11 = c08abcd c09abcd c10bdac c11bdac */ mulps BOF(%rsp), rC1/* */ movhlps rC12, rC13 /* rC13 = c12bd c13bd X X */ addps rA0, rC3 addps rC13, rC12 /* rC12 = c12abcd c13abcd X X */ addps rC4, rC7 addps rC5, rC11 prefB(320-176(pB,ldab)) addps rC1, rC12 #else /* BETA = X, complex type */ movups (pC), rA0 movlhps rC4, rC5 /* rC5 = c5a c5b c4a c5a */ movups 16(pC), rC4 unpcklps rC9 , rC8 /* rC8 = c08a c09a c08b c09b */ shufps $0x88, rC4, rA0 /* rA0 = c0 c1 c2 c3 */ movups 32(pC), rC4 /* rC4 = c4 X c5 X */ movhlps rC1, rC11 /* rC11 = c08d c09d c10c c11c */ movlhps rC10, rC1 /* rC1 = c08c c09c c10a c11a */ movhlps rC5, rC6 /* rC6 = c4a c5a c6b c7b */ movups 48(pC), rC5 /* rC5 = c6 X c7 X */ movhlps rC8 , rC0 /* rC0 = c08b c09b c10d c11d */ unpcklps rC13, rC2 /* rC2 = c12a c13a c12b c13b */ addps rC1, rC11 /* rC11 = c08cd c09cd c10ac c11ac */ shufps $0x88, rC5, rC4 /* rC4 = c4 c5 c6 c7 */ movups 64(pC), rC5 /* rC5 = c8 X c9 X */ movups 80(pC), rC1 /* rC1 = c10 X c11 X */ addps rB0, rC6 /* rC6 = c4ab c5ab c6bd c7bd */ shufps $0x88, rC1, rC5 /* rC5 = c8 c9 c10 c11 */ movlhps rC8 , rC9 /* rC9 = c09a c09b c08a c09a */ movss 96(pC), rC1 unpckhps rC13, rC12 /* rC12 = c12c c13c c12d c13d */ movss 104(pC), rB0 movhlps rC9 , rC10 /* rC10 = c08a c09a c10b c11b */ unpcklps rB0, rC1 addps rC6, rC7 /* rC7 = c4abcd c5abcd c6bdac c7bdac */ mulps BOF(%rsp), rA0 addps rC0, rC10 /* rC10 = c08ab c09ab c10bd c11bd */ mulps BOF(%rsp), rC4 addps rC2, rC12 /* rC12 = c12ac c13ac c12bd c13bd */ mulps BOF(%rsp), rC5 addps rC10, rC11 /* rC11 = c08abcd c09abcd c10bdac c11bdac */ mulps BOF(%rsp), rC1/* */ movhlps rC12, rC13 /* rC13 = c12bd c13bd X X */ addps rA0, rC3 addps rC13, rC12 /* rC12 = c12abcd c13abcd X X */ addps rC4, rC7 addps rC5, rC11 prefB(320-176(pB,ldab)) addps rC1, rC12 #endif#else movlhps rC4, rC5 /* rC5 = c5a c5b c4a c5a */ unpcklps rC9 , rC8 /* rC8 = c08a c09a c08b c09b */ movhlps rC1, rC11 /* rC11 = c08d c09d c10c c11c */ movlhps rC10, rC1 /* rC1 = c08c c09c c10a c11a */ movhlps rC5, rC6 /* rC6 = c4a c5a c6b c7b */ movhlps rC8 , rC0 /* rC0 = c08b c09b c10d c11d */ unpcklps rC13, rC2 /* rC2 = c12a c13a c12b c13b */ addps rC1, rC11 /* rC11 = c08cd c09cd c10ac c11ac */ addps rB0, rC6 /* rC6 = c4ab c5ab c6bd c7bd */ movlhps rC8 , rC9 /* rC9 = c09a c09b c08a c09a */ unpckhps rC13, rC12 /* rC12 = c12c c13c c12d c13d */ movhlps rC9 , rC10 /* rC10 = c08a c09a c10b c11b */ addps rC6, rC7 /* rC7 = c4abcd c5abcd c6bdac c7bdac */ addps rC0, rC10 /* rC10 = c08ab c09ab c10bd c11bd */ addps rC2, rC12 /* rC12 = c12ac c13ac c12bd c13bd */ addps rC10, rC11 /* rC11 = c08abcd c09abcd c10bdac c11bdac *//* */ movhlps rC12, rC13 /* rC13 = c12bd c13bd X X */ prefB(320-176(pB,ldab)) addps rC13, rC12 /* rC12 = c12abcd c13abcd X X */#endif/* * Write results back to C; pC += 14; */#ifdef SREAL movups rC3, (pC) movups rC7, 16(pC) movups rC11, 32(pC) movlps rC12, 48(pC)/* addq $56, pC */#else movss rC3, (pC) movss rC7, 32(pC) movhlps rC3, rC0 movhlps rC7, rC6 movss rC0, 16(pC) movss rC6, 48(pC) shufps $0x55, rC3, rC3 shufps $0x55, rC7, rC7 movss rC3, 8(pC) movss rC7, 40(pC) shufps $0x55, rC0, rC0 shufps $0x55, rC6, rC6 movss rC0, 24(pC) movss rC6, 56(pC) movss rC11, 64(pC) movhlps rC11, rC2 movss rC12, 96(pC) movss rC2, 80(pC) shufps $0x55, rC11, rC11 shufps $0x55, rC12, rC12 movss rC11, 72(pC) shufps $0x55, rC2, rC2 movss rC12, 104(pC) movss rC2, 88(pC)/* addq $112, pC */#endif/* * Write results back to C *//* * pC += 14; pA += 14*NB; pB -= NB; *//* * while (pA != stM); *//* subq $1, stM *//* jne UMLOOP *//* * pC += 14; pA += 14*NB; pB -= NB; */ subq $NBNBso-NB14so+176, pA5 subq $NBNBso-NB14so+176, pA10 addq $NBso-176, pB0/* * while (pA != stM); *//* subq $1, stM *//* jne UMLOOP *//* * pC += incCn; pA -= NBNB; pB += NB; */ addq incCn, pC/* * while (pB != stN); */ sub $1, stN jne UNLOOP/* * Restore callee-saved iregs */DONE: movq -8(%rsp), %rbp movq -16(%rsp), %rbx/* movq -32(%rsp), %r12 *//* movq -40(%rsp), %r13 */ ret⌨️ 快捷键说明
复制代码Ctrl + C
搜索代码Ctrl + F
全屏模式F11
增大字号Ctrl + =
减小字号Ctrl + -
显示快捷键?