⭐ 欢迎来到虫虫下载站! | 📦 资源下载 📁 资源专辑 ℹ️ 关于我们
⭐ 虫虫下载站

📄 ztrsm_kernel_lt.s

📁 Optimized GotoBLAS libraries
💻 S
📖 第 1 页 / 共 5 页
字号:
/*********************************************************************//*                                                                   *//*             Optimized BLAS libraries                              *//*                     By Kazushige Goto <kgoto@tacc.utexas.edu>     *//*                                                                   *//* Copyright (c) The University of Texas, 2005. All rights reserved. *//* UNIVERSITY EXPRESSLY DISCLAIMS ANY AND ALL WARRANTIES CONCERNING  *//* THIS SOFTWARE AND DOCUMENTATION, INCLUDING ANY WARRANTIES OF      *//* MERCHANTABILITY, FITNESS FOR ANY PARTICULAR PURPOSE,              *//* NON-INFRINGEMENT AND WARRANTIES OF PERFORMANCE, AND ANY WARRANTY  *//* THAT MIGHT OTHERWISE ARISE FROM COURSE OF DEALING OR USAGE OF     *//* TRADE. NO WARRANTY IS EITHER EXPRESS OR IMPLIED WITH RESPECT TO   *//* THE USE OF THE SOFTWARE OR DOCUMENTATION.                         *//* Under no circumstances shall University be liable for incidental, *//* special, indirect, direct or consequential damages or loss of     *//* profits, interruption of business, or related expenses which may  *//* arise from use of Software or Documentation, including but not    *//* limited to those resulting from defects in Software and/or        *//* Documentation, or loss or inaccuracy of data of any kind.         *//*********************************************************************/#define ASSEMBLER#include "common.h"#ifdef DOUBLE#define PREFETCHSIZE  (16 *  8)#else#define PREFETCHSIZE  (32 *  8)#endif#ifndef LN#define CPREFETCHSIZE  7#else#define CPREFETCHSIZE -8#endif#define CPREFETCH     lfetch.excl.nt1#define M	r32#define N	r33#define K	r34#define A	r37#define B	r38#define C	r39#define LDC	r35#define I	r15#define J	r16#define AOFFSET	r17#define BOFFSET	r18#define TEMP	r19#define L	r20#define C1	r21#define C2	r22#define C3	r23#define C4	r24#define C5	r25#define C6	r26#define C7	r27#define C8	r28#define PREA	r8#define PREB	r9#define PREC	r10#define SP	r12#define ARLC	r29#define PR	r30#define ARPFS	r31#define ALPHA_R	f8#define ALPHA_I	f9#define AORIG	loc0#define KK	loc1#define KK8	loc2#define OFFSET	loc3#define AOFFSET2 loc4#define BOFFSET2 loc5#ifndef CONJ#define FCALC_A	  FSUB#define FCALC_B	  FADD#define FMA_A	  FNMA#define FMA_B	  FMA#else#define FCALC_A	  FADD#define FCALC_B	  FSUB#define FMA_A	  FMA#define FMA_B	  FNMA#endif#ifndef CONJ#define FCALC_C	  FMA#define FCALC_D	  FNMA#else#define FCALC_C	  FNMA#define FCALC_D	  FMA#endif#ifndef CONJ#define FMA_C	FNMA#define FMA_D	FMA#define FSUB_A	FSUB#else#define FMA_C	FMA#define FMA_D	FMS#define FSUB_A	FADD#endif	PROLOGUE	.prologue	PROFCODE	{ .mfi	.save	ar.pfs, ARPFS	alloc	ARPFS = ar.pfs, 8, 8, 0, 0	mov	f64  = f0	adds	r14 = 16, SP	}	{ .mfi	nop	__LINE__	mov	f65  = f0	adds	r15 = 24, SP	}	;;	{ .mfi	ld8	LDC    = [r14]	mov	f81  = f0	mov	PR = pr	}	{ .mfi	ld8	OFFSET = [r15]	mov	f96  = f0	shr	J = N, 2	}	;;	{ .mfi	shladd	LDC = LDC, ZBASE_SHIFT, r0	mov	f97  = f0	}	{ .mfi	nop	__LINE__	mov	f113 = f0	}	;;#ifdef LN       { .mmi	setf.sig f32 = M	setf.sig f33 = K	shladd	C = M,  ZBASE_SHIFT, C	}	;;	{.mmf	nop	__LINE__	nop	__LINE__	xmpy.l	f32  = f32, f33	}	;;	{ .mmi	getf.sig r2 = f32	;;	nop	__LINE__	shladd	A = r2, ZBASE_SHIFT, A	}	;;#endif#ifdef RN	sub	KK  = r0, OFFSET#endif#ifdef RT       { .mmi	setf.sig f32 = N	setf.sig f33 = K	nop	 __LINE__	}	;;	{ .mmi	setf.sig f34 = LDC	nop	 __LINE__	nop	 __LINE__	}	;;	{ .mmf	nop	 __LINE__	nop	 __LINE__	xmpy.l	f33  = f32, f33	}	{ .mmf	nop	 __LINE__	sub	KK = N, OFFSET	xmpy.l	f34  = f32, f34	}	;;	{ .mmi	getf.sig r2 = f33	getf.sig r3 = f34	}	;;	shladd	B = r2, ZBASE_SHIFT, B	add	C = r3, C#endif	;;	.body	{ .mfi	nop	__LINE__	mov	f80  = f0	mov	ARLC  = ar.lc	}	{ .mfb	cmp.ge	p6, p0  = 0, J	mov	f112 = f0	(p6)	br.cond.dpnt .L050	}	;;	.align 16.L010:#ifdef RT       { .mmi	shladd	r3 = LDC, 2, r0	nop	__LINE__	shl	r2 = K, 2 + ZBASE_SHIFT	}	;;	{ .mmi	sub	B = B, r2	sub	C = C, r3	nop	__LINE__	}	;;#endif	{ .mmi	mov	C1 = C			// coffset1 = c + 0 * ldc	add	C2 = LDC, C		// coffset2 = c + 1 * ldc	shr	I  = M, 2	}	{ .mmi	adds	J = -1, J#ifdef LN	add	KK = M, OFFSET#elif defined LT	mov	KK = OFFSET#else	nop	__LINE__#endif#if defined(LN) || defined(RT)	mov	AORIG = A#else	mov	AOFFSET = A#endif	}	;;	;;	{ .mmi	shladd	C3 = LDC, 1, C		// coffset3 = c + 2 * ldc	shladd	C4 = LDC, 1, C2		// coffset4 = c + 3 * ldc#if defined(LT) || defined(RN)	mov	L = KK#else	sub	L = K, KK#endif	}	{ .mib	cmp.eq	p6, p7 = 0, I#ifndef RT	shladd	C = LDC, 2, C		// coffset += 8 * ldc#else	nop	__LINE__#endif	(p6)	br.cond.dpnt .L020	}	;;	.align 16.L011:	{ .mmi	cmp.ne	p7, p0 = r0, L	adds	BOFFSET = 0 * SIZE, B	shl	r2 = K, 2 + ZBASE_SHIFT	}	{ .mfi	shladd	r3 = KK, ZBASE_SHIFT, r0	mov	f118 = f0	nop	__LINE__	}	;;#if defined(LT) || defined(RN)	{ .mfb	(p7) LDFPD	f48, f49 = [BOFFSET], 2 * SIZE	mov	f66  = f0	nop	__LINE__	}	{ .mmf	nop	__LINE__	nop	__LINE__	mov	f67  = f0	}	;;#else	{ .mfi	shladd	BOFFSET = r3, 2, B	mov	f66  = f0#ifdef LN	sub	AORIG = AORIG, r2#else	nop	__LINE__#endif	}	;;	{ .mfi	(p7) LDFPD	f48, f49 = [BOFFSET], 2 * SIZE	mov	f67  = f0	shladd	AOFFSET = r3, 2, AORIG	}	;;#endif	;;	{ .mfi	(p7) LDFPD	f32, f33 = [AOFFSET], 2 * SIZE	mov	f82  = f0	nop	__LINE__	}	{ .mfi	(p7) LDFPD	f50, f51 = [BOFFSET], 2 * SIZE	mov	f83  = f0	adds	PREC = CPREFETCHSIZE * SIZE, C1	}	;;	{ .mfi	(p7) LDFPD	f34, f35  = [AOFFSET], 2 * SIZE	mov	f98  = f0	adds	L =  1, L	}	{ .mfi	(p7) LDFPD	f52, f53 = [BOFFSET], 2 * SIZE	mov	f99  = f0	adds	C5 = 4 * SIZE, C1	}	;;	{ .mfi	(p7) LDFPD	f36, f37  = [AOFFSET], 2 * SIZE	mov	f114 = f0	tbit.z	p12, p0 = L, 0	}	{ .mfi	(p7) LDFPD	f54, f55 = [BOFFSET], 2 * SIZE	mov	f115 = f0	adds	C6 = 4 * SIZE, C2	}	;;	{ .mfi	(p7) LDFPD	f38, f39  = [AOFFSET], 2 * SIZE	mov	f68  = f0	shr	L = L, 1	}	{ .mfi	setf.d	f86  = r0	mov	f69  = f0	adds	C7 = 4 * SIZE, C3	}	;;	{ .mfi	CPREFETCH [PREC], LDC	mov	f84  = f0	adds	L =  -1, L	}	{ .mfi	setf.d	f87  = r0	mov	f85  = f0	adds	C8 = 4 * SIZE, C4	}	;;	{ .mfi	CPREFETCH [PREC], LDC	mov	f100 = f0	mov	ar.lc = L	}	{ .mfi	setf.d	f102 = r0	mov	f101 = f0	cmp.eq	p3, p0 = r0, r0	}	;;	{ .mfi	CPREFETCH [PREC], LDC	mov	f116 = f0	adds	PREA = (PREFETCHSIZE + 0) * SIZE, AOFFSET	}	{ .mfi	setf.d	f103 = r0	mov	f117 = f0	adds	PREB = (PREFETCHSIZE + 0) * SIZE, BOFFSET	}	;;	{ .mfi	CPREFETCH [PREC]	mov	f70  = f0	cmp.eq  p6, p0 = -1, L	}	{ .mfb	setf.d	f119 = r0 	mov	f71  = f0	(p6) br.cond.dpnt   .L018	}	;;	.align 16.L012:/*  1 */	{ .mfi	lfetch.nt1	[PREA],  16 * SIZE	FMA	f64   = f32, f48, f64	// A1 * B1	nop	__LINE__	}	{ .mfb	(p12) cmp.ne p3, p0 =  0, L	FMA_B	f65   = f32, f49, f65	// A1 * B2	nop	__LINE__	}	;;/*  2 */	{ .mfi	lfetch.nt1	[PREB],  16 * SIZE	FMA	f80   = f32, f50, f80	// A1 * B3	nop	__LINE__	}	{ .mfb	cmp.ne	p4, p5 =  0, L	FMA_B	f81   = f32, f51, f81	// A1 * B4	nop	__LINE__	}	;;/*  3 */	{ .mfb	(p3) LDFPD	f56, f57 = [BOFFSET],   2 * SIZE	FMA	f96   = f32, f52, f96	// A1 * B5	nop	__LINE__	}	{ .mfb	FMA_B	f97  = f32, f53, f97	// A1 * B6	nop	__LINE__	}	;;/*  4 */	{ .mfb	(p3) LDFPD	f40, f41 = [AOFFSET], 2 * SIZE	FMA	f112  = f32, f54, f112	// A1 * B7	nop	__LINE__	}	{ .mfb	FMA_B	f113  = f32, f55, f113	// A1 * B8	nop	__LINE__	}	;;/*  5 */	{ .mfb	(p3) LDFPD	f58, f59 = [BOFFSET],  2 * SIZE	FMA	f65   = f33, f48, f65	// A2 * B1	nop	__LINE__	}	{ .mfb	FMA_A	f64   = f33, f49, f64	// A2 * B2	nop	__LINE__	}	;;/*  6 */	{ .mfb	(p3) LDFPD	f60, f61 = [BOFFSET], 2 * SIZE	FMA	f81   = f33, f50, f81	// A2 * B3	nop	__LINE__	}	{ .mfb	FMA_A	f80   = f33, f51, f80	// A2 * B4	nop	__LINE__	}	;;/*  7 */	{ .mfb	(p3) LDFPD	f62, f63 = [BOFFSET], 2 * SIZE	FMA	f97   = f33, f52, f97	// A2 * B5	nop	__LINE__	}	{ .mfb	FMA_A	f96  = f33, f53, f96	// A2 * B6	nop	__LINE__	}	;;/*  8 */	{ .mfb	(p3) LDFPD	f42, f43 = [AOFFSET], 2 * SIZE	FMA	f113  = f33, f54, f113	// A2 * B7	nop	__LINE__	}	{ .mfb	FMA_A	f112  = f33, f55, f112	// A2 * B8	nop	__LINE__	}	;;/*  9 */	{ .mfb	(p3) LDFPD	f44, f45 = [AOFFSET], 2 * SIZE	FMA	f66   = f34, f48, f66	// A3 * B1	nop	__LINE__	}	{ .mfb	FMA_B	f67   = f34, f49, f67	// A3 * B2	nop	__LINE__	}	;;/* 10 */	{ .mfb	(p3) LDFPD	f46, f47 = [AOFFSET], 2 * SIZE	FMA	f82   = f34, f50, f82	// A3 * B3	nop	__LINE__	}	{ .mfb	FMA_B	f83   = f34, f51, f83	// A3 * B4	nop	__LINE__	}	;;/* 11 */	{ .mfb	FMA	f98   = f34, f52, f98	// A3 * B5	nop	__LINE__	}	{ .mfb	nop	__LINE__	FMA_B	f99  = f34, f53, f99	// A3 * B6	nop	__LINE__	}	;; /* 12 */	{ .mfb	FMA	f114  = f34, f54, f114	// A3 * B7	nop	__LINE__	}	{ .mfb	nop	__LINE__	FMA_B	f115  = f34, f55, f115	// A3 * B8	nop	__LINE__	}	;;/* 13 */	{ .mfb	nop	__LINE__	FMA	f67   = f35, f48, f67	// A4 * B1	}	{ .mfb	nop	__LINE__	FMA_A	f66   = f35, f49, f66	// A4 * B2	nop	__LINE__	}	;;/* 14 */	{ .mfb	FMA	f83   = f35, f50, f83	// A4 * B3	nop	__LINE__	}	{ .mfb	nop	__LINE__	FMA_A	f82   = f35, f51, f82	// A4 * B4	nop	__LINE__	}	;;/* 15 */	{ .mfb	FMA	f99   = f35, f52, f99	// A4 * B5	nop	__LINE__	}	{ .mfb	nop	__LINE__	FMA_A	f98  = f35, f53, f98	// A4 * B6	nop	__LINE__	}	;;/* 16 */	{ .mfb	FMA	f115  = f35, f54, f115	// A4 * B7	nop	__LINE__	}	{ .mfb	nop	__LINE__	FMA_A	f114  = f35, f55, f114	// A4 * B8	nop	__LINE__	}	;;/* 17 */	{ .mfb	nop	__LINE__	FMA	f68   = f36, f48, f68	// A5 * B1	nop	__LINE__	}	{ .mfb	nop	__LINE__	FMA_B	f69   = f36, f49, f69	// A5 * B2	nop	__LINE__	}	;;/* 18 */	{ .mfb	nop	__LINE__	FMA	f84   = f36, f50, f84	// A5 * B3	nop	__LINE__	}	{ .mfb	nop	__LINE__	FMA_B	f85   = f36, f51, f85	// A5 * B4	nop	__LINE__	}	;;/* 19 */	{ .mfb	nop	__LINE__	FMA	f100  = f36, f52, f100	// A5 * B5	nop	__LINE__	}	{ .mfb	nop	__LINE__	FMA_B	f101  = f36, f53, f101	// A5 * B6	nop	__LINE__	}	;;/* 20 */	{ .mfb	nop	__LINE__	FMA	f116  = f36, f54, f116	// A5 * B7	nop	__LINE__	}	{ .mfb	nop	__LINE__	FMA_B	f117  = f36, f55, f117	// A5 * B8	nop	__LINE__	}	;;/* 21 */	{ .mfb	nop	__LINE__	FMA	f69   = f37, f48, f69	// A6 * B1	nop	__LINE__	}	{ .mfb	nop	__LINE__	FMA_A	f68   = f37, f49, f68	// A6 * B2	nop	__LINE__	}	;;/* 22 */	{ .mfb	nop	__LINE__	FMA	f85   = f37, f50, f85	// A6 * B3	nop	__LINE__	}	{ .mfb	nop	__LINE__	FMA_A	f84   = f37, f51, f84	// A6 * B4	nop	__LINE__	}	;;/* 23 */	{ .mfb	nop	__LINE__	FMA	f101  = f37, f52, f101	// A6 * B5	nop	__LINE__	}	{ .mfb	nop	__LINE__	FMA_A	f100  = f37, f53, f100	// A6 * B6	nop	__LINE__	}	;;/* 24 */	{ .mfb	nop	__LINE__	FMA	f117  = f37, f54, f117	// A6 * B7	nop	__LINE__	}	{ .mfb	nop	__LINE__	FMA_A	f116  = f37, f55, f116	// A6 * B8	nop	__LINE__	}	;;/* 25 */	{ .mfb	nop	__LINE__	FMA	f70   = f38, f48, f70	// A7 * B1	nop	__LINE__	}	{ .mfb	nop	__LINE__	FMA_B	f71   = f38, f49, f71	// A7 * B2	nop	__LINE__	}	;;/* 26 */	{ .mfb	nop	__LINE__	FMA	f86   = f38, f50, f86	// A7 * B3	nop	__LINE__	}	{ .mfb	nop	__LINE__	FMA_B	f87   = f38, f51, f87	// A7 * B4	nop	__LINE__	}	;;/* 27 */	{ .mfb	nop	__LINE__	FMA	f102  = f38, f52, f102	// A7 * B5	nop	__LINE__	}	{ .mfb	nop	__LINE__	FMA_B	f103  = f38, f53, f103	// A7 * B6	nop	__LINE__	}	;;/* 28 */	{ .mfb	nop	__LINE__	FMA	f118  = f38, f54, f118	// A7 * B7	nop	__LINE__	}	{ .mfb	nop	__LINE__	FMA_B	f119  = f38, f55, f119	// A7 * B8	nop	__LINE__	}	;;/* 29 */	{ .mfb	nop	__LINE__	FMA	f71   = f39, f48, f71	// A8 * B1	nop	__LINE__	}	{ .mfb	nop	__LINE__	FMA_A	f70   = f39, f49, f70	// A8 * B2	nop	__LINE__	}	;;/* 30 */	{ .mfb	(p4) LDFPD	f32, f33 = [AOFFSET],   2 * SIZE	FMA	f87   = f39, f50, f87	// A8 * B3	nop	__LINE__	}	{ .mfb	nop	__LINE__	FMA_A	f86   = f39, f51, f86	// A8 * B4	nop	__LINE__	}	;;/* 31 */	{ .mfb	(p4) LDFPD	f48, f49 = [BOFFSET],   2 * SIZE	FMA	f103  = f39, f52, f103	// A8 * B5	nop	__LINE__	}	{ .mfb	nop	__LINE__	FMA_A	f102  = f39, f53, f102	// A8 * B6	nop	__LINE__	}	;;/* 32 */	{ .mfb	nop	__LINE__	FMA	f119  = f39, f54, f119	// A8 * B7	nop	__LINE__	}	{ .mfb	nop	__LINE__	FMA_A	f118  = f39, f55, f118	// A8 * B8	nop	__LINE__	}	;;/* 33 */	{ .mfb	nop	__LINE__	(p3) FMA	f64   = f40, f56, f64	// A1 * B1	nop	__LINE__	}	{ .mfb	nop	__LINE__	(p3) FMA_B	f65   = f40, f57, f65	// A1 * B2	nop	__LINE__	}	;;/* 34 */	{ .mfb	(p4) LDFPD	f50, f51 = [BOFFSET],  2 * SIZE	(p3) FMA	f80   = f40, f58, f80	// A1 * B3	nop	__LINE__	}	{ .mfb	nop	__LINE__	(p3) FMA_B	f81   = f40, f59, f81	// A1 * B4	nop	__LINE__	}	;;/* 35 */	{ .mfb	(p4) LDFPD	f52, f53 = [BOFFSET], 2 * SIZE	(p3) FMA	f96   = f40, f60, f96	// A1 * B5	nop	__LINE__	}	{ .mfb	nop	__LINE__	(p3) FMA_B	f97  = f40, f61, f97	// A1 * B6	nop	__LINE__	}	;;/* 36 */	{ .mfb	(p4) LDFPD	f54, f55 = [BOFFSET], 2 * SIZE	(p3) FMA	f112  = f40, f62, f112	// A1 * B7	nop	__LINE__	}	{ .mfb	nop	__LINE__	(p3) FMA_B	f113  = f40, f63, f113	// A1 * B8	nop	__LINE__	}	;;/* 37 */	{ .mfb	(p4) LDFPD	f34, f35 = [AOFFSET], 2 * SIZE	(p3) FMA	f65   = f41, f56, f65	// A2 * B1	nop	__LINE__	}	{ .mfb	nop	__LINE__	(p3) FMA_A	f64   = f41, f57, f64	// A2 * B2	nop	__LINE__	}	;;/* 38 */	{ .mfb	(p4) LDFPD	f36, f37 = [AOFFSET], 2 * SIZE	(p3) FMA	f81   = f41, f58, f81	// A2 * B3	nop	__LINE__	}	{ .mfb	nop	__LINE__	(p3) FMA_A	f80   = f41, f59, f80	// A2 * B4	nop	__LINE__	}	;;/* 39 */	{ .mfb	(p4) LDFPD	f38, f39 = [AOFFSET], 2 * SIZE	(p3) FMA	f97   = f41, f60, f97	// A2 * B5	nop	__LINE__	}	{ .mfb	nop	__LINE__	(p3) FMA_A	f96  = f41, f61, f96	// A2 * B6	nop	__LINE__	}	;;/* 40 */	{ .mfb	nop	__LINE__	(p3) FMA	f113  = f41, f62, f113	// A2 * B7	nop	__LINE__	}	{ .mfb	nop	__LINE__	(p3) FMA_A	f112  = f41, f63, f112	// A2 * B8	nop	__LINE__	}	;; /* 41 */	{ .mfb	nop	__LINE__	(p3) FMA	f66   = f42, f56, f66	// A3 * B1	nop	__LINE__	}	{ .mfb	nop	__LINE__	(p3) FMA_B	f67   = f42, f57, f67	// A3 * B2	nop	__LINE__	}	;;/* 42 */	{ .mfb	nop	__LINE__	(p3) FMA	f82   = f42, f58, f82	// A3 * B3	nop	__LINE__	}	{ .mfb	nop	__LINE__	(p3) FMA_B	f83   = f42, f59, f83	// A3 * B4	nop	__LINE__	}	;;/* 43 */	{ .mfb	nop	__LINE__	(p3) FMA	f98   = f42, f60, f98	// A3 * B5	nop	__LINE__	}	{ .mfb	nop	__LINE__	(p3) FMA_B	f99  = f42, f61, f99	// A3 * B6	nop	__LINE__	}	;;/* 44 */	{ .mfb	nop	__LINE__	(p3) FMA	f114  = f42, f62, f114	// A3 * B7	nop	__LINE__	}

⌨️ 快捷键说明

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