📄 ztrsm_kernel_lt_4x2_sse.s

📁 Optimized GotoBLAS libraries
💻 S
📖 第 1 页 / 共 5 页
字号:
	movaps	 48 * SIZE(AO), %xmm10	addq   $32  * SIZE, AO	addq   $128 * SIZE, BO	decq   %rax	jne    .L22	ALIGN_4.L25:#if defined(LT) || defined(RN)	movq	KK, %rax#else	movq	K, %rax	subq	KK, %rax#endif	movaps	POSINV,  %xmm15	andq	$7, %rax		# if (k & 1)	BRANCH	je .L28	ALIGN_4.L26:	mulps	%xmm8, %xmm9	addps	%xmm9, %xmm0	movaps	 4 * SIZE(BO), %xmm9	mulps	%xmm8, %xmm9	addps	%xmm9, %xmm1	movaps	 8 * SIZE(BO), %xmm9	mulps	%xmm8, %xmm9	mulps	12 * SIZE(BO), %xmm8	addps	%xmm9, %xmm2	movaps	 16 * SIZE(BO), %xmm9	addps	%xmm8, %xmm3	movaps	  4 * SIZE(AO), %xmm8	addq	$ 4 * SIZE, AO		# aoffset  += 4	addq	$16 * SIZE, BO		# boffset1 += 8	decq	%rax	jg	.L26	ALIGN_4.L28:	shufps	$0xb1, %xmm1, %xmm1	shufps	$0xb1, %xmm3, %xmm3#if defined(LN) || defined(LT)#ifndef CONJ	xorps	 %xmm15, %xmm1	xorps	 %xmm15, %xmm3#else	xorps	 %xmm15, %xmm0	xorps	 %xmm15, %xmm2#endif#else	xorps	 %xmm15, %xmm1	xorps	 %xmm15, %xmm3#endif	addps	%xmm1, %xmm0	addps	%xmm3, %xmm2#if defined(LN) || defined(RT)	movq	KK, %rax#ifdef LN	subq	$2, %rax#else	subq	$2, %rax#endif	movq	AORIG, AO	movq	BORIG, B	leaq	BUFFER, BO	salq	$ZBASE_SHIFT, %rax	leaq	(AO, %rax, 2), AO	leaq	(B,  %rax, 2), B	leaq	(BO, %rax, 8), BO#endif#if defined(LN) || defined(LT)	movaps	%xmm0, %xmm1	unpcklpd %xmm2, %xmm0	unpckhpd %xmm2, %xmm1	movaps	 0 * SIZE(B), %xmm2	movaps	 4 * SIZE(B), %xmm3	subps	%xmm0,  %xmm2	subps	%xmm1,  %xmm3#else	movaps	 0 * SIZE(AO), %xmm1	movaps	 4 * SIZE(AO), %xmm5	subps	%xmm0,  %xmm1	subps	%xmm2,  %xmm5#endif#ifdef LN	movaps	 4 * SIZE(AO), %xmm8	pshufd	 $0xee, %xmm8, %xmm9	pshufd	 $0xbb, %xmm8, %xmm10	pshufd	 $0xa0, %xmm3, %xmm0	pshufd	 $0xf5, %xmm3, %xmm3#ifndef CONJ	xorps	 %xmm15, %xmm3#else	xorps	 %xmm15, %xmm0#endif	mulps	%xmm9,  %xmm0	mulps	%xmm10, %xmm3	addps	%xmm0,  %xmm3	pshufd	 $0x44, %xmm8, %xmm9	pshufd	 $0x11, %xmm8, %xmm10	pshufd	 $0xa0, %xmm3, %xmm0	pshufd	 $0xf5, %xmm3, %xmm1#ifndef CONJ	xorps	 %xmm15, %xmm1#else	xorps	 %xmm15, %xmm0#endif	mulps	%xmm9,  %xmm0	mulps	%xmm10, %xmm1	subps	%xmm0,  %xmm2	subps	%xmm1,  %xmm2	movaps	 0 * SIZE(AO), %xmm8	pshufd	 $0x44, %xmm8, %xmm9	pshufd	 $0x11, %xmm8, %xmm10	pshufd	 $0xa0, %xmm2, %xmm0	pshufd	 $0xf5, %xmm2, %xmm2#ifndef CONJ	xorps	 %xmm15, %xmm2#else	xorps	 %xmm15, %xmm0#endif	mulps	%xmm9,  %xmm0	mulps	%xmm10, %xmm2	addps	%xmm0,  %xmm2#endif#ifdef LT	movaps	 0 * SIZE(AO), %xmm8	pshufd	 $0x44, %xmm8, %xmm9	pshufd	 $0x11, %xmm8, %xmm10	pshufd	 $0xa0, %xmm2, %xmm0	pshufd	 $0xf5, %xmm2, %xmm2#ifndef CONJ	xorps	 %xmm15, %xmm2#else	xorps	 %xmm15, %xmm0#endif	mulps	%xmm9,  %xmm0	mulps	%xmm10, %xmm2	addps	%xmm0,  %xmm2	pshufd	 $0xee, %xmm8, %xmm9	pshufd	 $0xbb, %xmm8, %xmm10	pshufd	 $0xa0, %xmm2, %xmm0	pshufd	 $0xf5, %xmm2, %xmm1#ifndef CONJ	xorps	 %xmm15, %xmm1#else	xorps	 %xmm15, %xmm0#endif	mulps	%xmm9,  %xmm0	mulps	%xmm10, %xmm1	subps	%xmm0,  %xmm3	subps	%xmm1,  %xmm3	movaps	 4 * SIZE(AO), %xmm8	pshufd	 $0xee, %xmm8, %xmm9	pshufd	 $0xbb, %xmm8, %xmm10	pshufd	 $0xa0, %xmm3, %xmm0	pshufd	 $0xf5, %xmm3, %xmm3#ifndef CONJ	xorps	 %xmm15, %xmm3#else	xorps	 %xmm15, %xmm0#endif	mulps	%xmm9,  %xmm0	mulps	%xmm10, %xmm3	addps	%xmm0,  %xmm3#endif#ifdef RN	movaps	 0 * SIZE(B), %xmm8	pshufd	 $0x44, %xmm8, %xmm9	pshufd	 $0x11, %xmm8, %xmm10	pshufd	 $0xa0, %xmm1, %xmm0	pshufd	 $0xf5, %xmm1, %xmm1#ifndef CONJ	xorps	 %xmm15, %xmm1#else	xorps	 %xmm15, %xmm0#endif	mulps	%xmm9,  %xmm0	mulps	%xmm10, %xmm1	addps	%xmm0,  %xmm1	pshufd	 $0xee, %xmm8, %xmm9	pshufd	 $0xbb, %xmm8, %xmm10	pshufd	 $0xa0, %xmm1, %xmm0	pshufd	 $0xf5, %xmm1, %xmm2#ifndef CONJ	xorps	 %xmm15, %xmm2#else	xorps	 %xmm15, %xmm0#endif	mulps	%xmm9,  %xmm0	mulps	%xmm10, %xmm2	subps	%xmm0,  %xmm5	subps	%xmm2,  %xmm5	movaps	 4 * SIZE(B), %xmm8	pshufd	 $0xee, %xmm8, %xmm9	pshufd	 $0xbb, %xmm8, %xmm10	pshufd	 $0xa0, %xmm5, %xmm4	pshufd	 $0xf5, %xmm5, %xmm5#ifndef CONJ	xorps	 %xmm15, %xmm5#else	xorps	 %xmm15, %xmm4#endif	mulps	%xmm9,  %xmm4	mulps	%xmm10, %xmm5	addps	%xmm4,  %xmm5#endif#ifdef RT	movaps	 4 * SIZE(B), %xmm8	pshufd	 $0xee, %xmm8, %xmm9	pshufd	 $0xbb, %xmm8, %xmm10	pshufd	 $0xa0, %xmm5, %xmm0	pshufd	 $0xf5, %xmm5, %xmm5#ifndef CONJ	xorps	 %xmm15, %xmm5#else	xorps	 %xmm15, %xmm0#endif	mulps	%xmm9,  %xmm0	mulps	%xmm10, %xmm5	addps	%xmm0,  %xmm5	pshufd	 $0x44, %xmm8, %xmm9	pshufd	 $0x11, %xmm8, %xmm10	pshufd	 $0xa0, %xmm5, %xmm0	pshufd	 $0xf5, %xmm5, %xmm2#ifndef CONJ	xorps	 %xmm15, %xmm2#else	xorps	 %xmm15, %xmm0#endif	mulps	%xmm9,  %xmm0	mulps	%xmm10, %xmm2	subps	%xmm0,  %xmm1	subps	%xmm2,  %xmm1	movaps	 0 * SIZE(B), %xmm8	pshufd	 $0x44, %xmm8, %xmm9	pshufd	 $0x11, %xmm8, %xmm10	pshufd	 $0xa0, %xmm1, %xmm4	pshufd	 $0xf5, %xmm1, %xmm1#ifndef CONJ	xorps	 %xmm15, %xmm1#else	xorps	 %xmm15, %xmm4#endif	mulps	%xmm9,  %xmm4	mulps	%xmm10, %xmm1	addps	%xmm4,  %xmm1#endif#ifdef LN	subq	$4 * SIZE, CO1	subq	$4 * SIZE, CO2#endif#if defined(LN) || defined(LT)	movaps	%xmm2,   0 * SIZE(B)	movaps	%xmm3,   4 * SIZE(B)	pshufd	$0x00, %xmm2, %xmm0	pshufd	$0x55, %xmm2, %xmm1	pshufd	$0xaa, %xmm2, %xmm4	pshufd	$0xff, %xmm2, %xmm5	movaps	%xmm0,   0 * SIZE(BO)	movaps	%xmm1,   4 * SIZE(BO)	movaps	%xmm4,   8 * SIZE(BO)	movaps	%xmm5,  12 * SIZE(BO)	pshufd	$0x00, %xmm3, %xmm0	pshufd	$0x55, %xmm3, %xmm1	pshufd	$0xaa, %xmm3, %xmm4	pshufd	$0xff, %xmm3, %xmm5	movaps	%xmm0,  16 * SIZE(BO)	movaps	%xmm1,  20 * SIZE(BO)	movaps	%xmm4,  24 * SIZE(BO)	movaps	%xmm5,  28 * SIZE(BO)	movlps	%xmm2,   0 * SIZE(CO1)	movlps	%xmm3,   2 * SIZE(CO1)	movhps	%xmm2,   0 * SIZE(CO2)	movhps	%xmm3,   2 * SIZE(CO2)#else	movaps	%xmm1,   0 * SIZE(AO)	movaps	%xmm5,   4 * SIZE(AO)	movlps	%xmm1,   0 * SIZE(CO1)	movhps	%xmm1,   2 * SIZE(CO1)	movlps	%xmm5,   0 * SIZE(CO2)	movhps	%xmm5,   2 * SIZE(CO2)#endif#ifndef LN	addq	$4 * SIZE, CO1	addq	$4 * SIZE, CO2#endif#if defined(LT) || defined(RN)	movq	K,  %rax	subq	KK, %rax	salq	$ZBASE_SHIFT, %rax	leaq	(AO, %rax, 2), AO#ifdef LT	addq	$8 * SIZE, B#endif#endif#ifdef LN	subq	$2, KK	movq	BORIG, B#endif#ifdef LT	addq	$2, KK#endif#ifdef RT	movq	K, %rax	movq	BORIG, B	salq	$1 + ZBASE_SHIFT, %rax	addq	%rax, AORIG#endif	ALIGN_4		.L30:	testq	$1, M	je	.L39#ifdef LN       movq	K, %rax       salq	$ZBASE_SHIFT, %rax       subq	%rax, AORIG#endif#if defined(LN) || defined(RT)	movq	KK, %rax	movq	AORIG, AO	salq	$ZBASE_SHIFT, %rax	leaq	(AO, %rax, 1), AO#endif	leaq	BUFFER, BO#if defined(LN) || defined(RT)	movq	KK, %rax	salq	$1 + ZBASE_SHIFT, %rax	leaq	(BO, %rax, 4), BO#endif		movlps	 0 * SIZE(AO), %xmm8	movhps	 2 * SIZE(AO), %xmm8	movlps	 8 * SIZE(AO), %xmm10	movhps	10 * SIZE(AO), %xmm10	movaps	 0 * SIZE(BO), %xmm9	movaps	16 * SIZE(BO), %xmm11	movaps	32 * SIZE(BO), %xmm13	movaps	48 * SIZE(BO), %xmm15	pxor	%xmm0, %xmm0	pxor	%xmm1, %xmm1	pxor	%xmm2, %xmm2	pxor	%xmm3, %xmm3#if defined(LT) || defined(RN)	movq	KK, %rax#else	movq	K, %rax	subq	KK, %rax#endif	sarq	$3, %rax	je	.L35	ALIGN_4.L32:	mulps	%xmm8, %xmm9	addps	%xmm9, %xmm0#if defined(OPTERON) && defined(HAVE_PREFETCH)	PREFETCH	(PREFETCHSIZE +  0) * SIZE(AO)#endif	movsd	  4 * SIZE(BO), %xmm9	mulps	%xmm8, %xmm9	addps	%xmm9, %xmm1	movsd	  8 * SIZE(BO), %xmm9	mulps	%xmm8, %xmm9	addps	%xmm9, %xmm2	movsd	 12 * SIZE(BO), %xmm9	mulps	%xmm8, %xmm9	movsd	  2 * SIZE(AO), %xmm8	addps	%xmm9, %xmm3	movsd	 64 * SIZE(BO), %xmm9	mulps	%xmm8, %xmm11	addps	%xmm11, %xmm0	movsd	 20 * SIZE(BO), %xmm11	mulps	%xmm8, %xmm11	addps	%xmm11, %xmm1	movsd	 24 * SIZE(BO), %xmm11	mulps	%xmm8, %xmm11	addps	%xmm11, %xmm2	movsd	 28 * SIZE(BO), %xmm11	mulps	%xmm8, %xmm11	movsd	  4 * SIZE(AO), %xmm8	addps	%xmm11, %xmm3	movsd	 80 * SIZE(BO), %xmm11	mulps	%xmm8, %xmm13	addps	%xmm13, %xmm0	movsd	 36 * SIZE(BO), %xmm13	mulps	%xmm8, %xmm13	addps	%xmm13, %xmm1	movsd	 40 * SIZE(BO), %xmm13	mulps	%xmm8, %xmm13	addps	%xmm13, %xmm2	movsd	 44 * SIZE(BO), %xmm13	mulps	%xmm8, %xmm13	movsd	  6 * SIZE(AO), %xmm8	addps	%xmm13, %xmm3	movsd	 96 * SIZE(BO), %xmm13	mulps	%xmm8, %xmm15	addps	%xmm15, %xmm0	movsd	 52 * SIZE(BO), %xmm15	mulps	%xmm8, %xmm15	addps	%xmm15, %xmm1	movsd	 56 * SIZE(BO), %xmm15	mulps	%xmm8, %xmm15	addps	%xmm15, %xmm2	movsd	 60 * SIZE(BO), %xmm15	mulps	%xmm8, %xmm15	movsd	 16 * SIZE(AO), %xmm8	addps	%xmm15, %xmm3	movsd	112 * SIZE(BO), %xmm15	mulps	%xmm10, %xmm9	addps	%xmm9, %xmm0	movsd	 68 * SIZE(BO), %xmm9	mulps	%xmm10, %xmm9	addps	%xmm9, %xmm1	movsd	 72 * SIZE(BO), %xmm9	mulps	%xmm10, %xmm9	addps	%xmm9, %xmm2	movsd	 76 * SIZE(BO), %xmm9	mulps	%xmm10, %xmm9	movsd	 10 * SIZE(AO), %xmm10	addps	%xmm9, %xmm3	movsd	128 * SIZE(BO), %xmm9	mulps	%xmm10, %xmm11	addps	%xmm11, %xmm0	movsd	 84 * SIZE(BO), %xmm11	mulps	%xmm10, %xmm11	addps	%xmm11, %xmm1	movsd	 88 * SIZE(BO), %xmm11	mulps	%xmm10, %xmm11	addps	%xmm11, %xmm2	movsd	 92 * SIZE(BO), %xmm11	mulps	%xmm10, %xmm11	movsd	 12 * SIZE(AO), %xmm10	addps	%xmm11, %xmm3	movsd	144 * SIZE(BO), %xmm11	mulps	%xmm10, %xmm13	addps	%xmm13, %xmm0	movsd	100 * SIZE(BO), %xmm13	mulps	%xmm10, %xmm13	addps	%xmm13, %xmm1	movsd	104 * SIZE(BO), %xmm13	mulps	%xmm10, %xmm13	addps	%xmm13, %xmm2	movsd	108 * SIZE(BO), %xmm13	mulps	%xmm10, %xmm13	movsd	 14 * SIZE(AO), %xmm10	addps	%xmm13, %xmm3	movsd	160 * SIZE(BO), %xmm13	mulps	%xmm10, %xmm15	addps	%xmm15, %xmm0	movsd	116 * SIZE(BO), %xmm15	mulps	%xmm10, %xmm15	addps	%xmm15, %xmm1	movsd	120 * SIZE(BO), %xmm15	mulps	%xmm10, %xmm15	addps	%xmm15, %xmm2	movsd	124 * SIZE(BO), %xmm15	mulps	%xmm10, %xmm15	movsd	 24 * SIZE(AO), %xmm10	addps	%xmm15, %xmm3	movsd	176 * SIZE(BO), %xmm15	addq   $16  * SIZE, AO	addq   $128 * SIZE, BO	decq   %rax	jne    .L32	ALIGN_4.L35:#if defined(LT) || defined(RN)	movq	KK, %rax#else	movq	K, %rax	subq	KK, %rax#endif	movaps	POSINV,  %xmm15	andq	$7, %rax		# if (k & 1)	BRANCH	je .L38	ALIGN_4.L36:	mulps	%xmm8, %xmm9	addps	%xmm9, %xmm0	movsd	 4 * SIZE(BO), %xmm9	mulps	%xmm8, %xmm9	addps	%xmm9, %xmm1	movsd	 8 * SIZE(BO), %xmm9	mulps	%xmm8, %xmm9	addps	%xmm9, %xmm2	movsd	12 * SIZE(BO), %xmm9	mulps	%xmm8, %xmm9	movsd	 2 * SIZE(AO), %xmm8	addps	%xmm9, %xmm3	movsd	16 * SIZE(BO), %xmm9	addq	$ 2 * SIZE, AO		# aoffset  += 4	addq	$16 * SIZE, BO		# boffset1 += 8	decq	%rax	jg	.L36	ALIGN_4.L38:	shufps	$0xb1, %xmm1, %xmm1	shufps	$0xb1, %xmm3, %xmm3#if defined(LN) || defined(LT)#ifndef CONJ	xorps	 %xmm15, %xmm1	xorps	 %xmm15, %xmm3#else	xorps	 %xmm15, %xmm0	xorps	 %xmm15, %xmm2#endif#else	xorps	 %xmm15, %xmm1	xorps	 %xmm15, %xmm3#endif	addps	%xmm1, %xmm0	addps	%xmm3, %xmm2#if defined(LN) || defined(RT)	movq	KK, %rax#ifdef LN	subq	$1, %rax#else	subq	$2, %rax#endif	movq	AORIG, AO	movq	BORIG, B	leaq	BUFFER, BO	salq	$ZBASE_SHIFT, %rax	leaq	(AO, %rax, 1), AO	leaq	(B,  %rax, 2), B	leaq	(BO, %rax, 8), BO#endif#if defined(LN) || defined(LT)	unpcklpd %xmm2, %xmm0	movaps	 0 * SIZE(B), %xmm2	subps	%xmm0,  %xmm2#else	movlps	 0 * SIZE(AO), %xmm1	movlps	 2 * SIZE(AO), %xmm5	subps	%xmm0,  %xmm1	subps	%xmm2,  %xmm5#endif#ifdef LN	movaps	 0 * SIZE(AO), %xmm8	pshufd	 $0x44, %xmm8, %xmm9	pshufd	 $0x11, %xmm8, %xmm10	pshufd	 $0xa0, %xmm2, %xmm0	pshufd	 $0xf5, %xmm2, %xmm2#ifndef CONJ	xorps	 %xmm15, %xmm2#else	xorps	 %xmm15, %xmm0#endif	mulps	%xmm9,  %xmm0	mulps	%xmm10, %xmm2	addps	%xmm0,  %xmm2#endif#ifdef LT	movaps	 0 * SIZE(AO), %xmm8	pshufd	 $0x44, %xmm8, %xmm9	pshufd	 $0x11, %xmm8, %xmm10	pshufd	 $0xa0, %xmm2, %xmm0	pshufd	 $0xf5, %xmm2, %xmm2#ifndef CONJ	xorps	 %xmm15, %xmm2#else	xorps	 %xmm15, %xmm0#endif	mulps	%xmm9,  %xmm0	mulps	%xmm10, %xmm2	addps	%xmm0,  %xmm2#endif#ifdef RN	movaps	 0 * SIZE(B), %xmm8	pshufd	 $0x44, %xmm8, %xmm9	pshufd	 $0x11, %xmm8, %xmm10	pshufd	 $0xa0, %xmm1, %xmm0	pshufd	 $0xf5, %xmm1, %xmm1#ifndef CONJ	xorps	 %xmm15, %xmm1#else	xorps	 %xmm15, %xmm0#endif	mulps	%xmm9,  %xmm0	mulps	%xmm10, %xmm1	addps	%xmm0,  %xmm1	pshufd	 $0xee, %xmm8, %xmm9	pshufd	 $0xbb, %xmm8, %xmm10	pshufd	 $0xa0, %xmm1, %xmm0	pshufd	 $0xf5, %xmm1, %xmm2#ifndef CONJ	xorps	 %xmm15, %xmm2#else	xorps	 %xmm15, %xmm0#endif	mulps	%xmm9,  %xmm0	mulps	%xmm10, %xmm2	subps	%xmm0,  %xmm5	subps	%xmm2,  %xmm5
⌨️ 快捷键说明

复制代码 Ctrl + C
搜索代码 Ctrl + F
全屏模式 F11
切换主题 Ctrl + Shift + D
显示快捷键 ?
增大字号 Ctrl + =
减小字号 Ctrl + -