OpenBLAS/kernel/mips64/gemm_kernel_loongson3a.S

#define gsLQC1(base,fq,ft,offset) .word(0x32<<26|base<<21|ft<<16|0x1<<15|offset<<6|0x1<<5|fq)
#define gsSQC1(base,fq,ft,offset) .word(0x3A<<26|base<<21|ft<<16|0x1<<15|offset<<6|0x1<<5|fq)
#define FETCH   ld

#define REALNAME ASMNAME
#define ASSEMBLER
#include "common.h"

#define M	$4
#define	N	$5
#define	K	$6
#define A	$8
#define B	$9
#define C	$10
#define LDC	$11

#define AO	$12
#define BO	$13

#define I	$2
#define J	$3
#define L	$7

#define CO1	$14
#define CO2	$15
#define CO3	$16
#define CO4	$17

#define KCO	$18
#define MCO	$19
#define NCO	$20

#define SPANB	$21
#define SPANC	$22
#define PREB	$23
#define PREA	$24
#define SPANA	$25

#define ALPHA	$f15

#define R8	8
#define	R9	9
#define R14	14
#define R15	15
#define R16	16
#define R17 17

#define	t11	$f30
#define	t21	$f31
#define	t31	$f28
#define	t41	$f29

#define	t12	$f26
#define	t22	$f27
#define	t32	$f24
#define	t42	$f25

#define	t13	$f22
#define	t23	$f23
#define	t33	$f20
#define	t43	$f21

#define	t14	$f18
#define	t24	$f19
#define	t34	$f16
#define	t44	$f17

#define	c11	$f0
#define	c21	$f1
#define	c31	$f2
#define	c41	$f3

#define	c12	$f4
#define	c22	$f5
#define	c32	$f6
#define	c42	$f7

#define	c13	$f8
#define	c23	$f9
#define	c33	$f10
#define c43	$f11

#define	c14	$f12
#define	c24	$f13
#define	c34	$f14
#define	c44	$f0

#define	a0	$f0
#define	a1	$f1
#define	a2	$f2
#define	a3	$f3
#define	a4	$f4
#define	a5	$f5
#define	a6	$f6
#define	a7	$f7
#define	b0	$f8
#define	b1	$f9
#define	b2	$f10
#define b3	$f11
#define	b4	$f12
#define	b5	$f13
#define	b6	$f14
#define	b7	$f15

#define F31 31
#define F30 30
#define F29 29
#define F28 28
#define F27 27
#define F26 26
#define F25 25
#define F24 24 
#define F23 23
#define F22 22
#define F21 21
#define F20 20
#define F19 19
#define F18 18
#define F17 17
#define F16 16 
#define F15 15
#define F14 14
#define F13 13
#define F12 12
#define F11 11
#define F10 10
#define F9 9
#define F8 8
#define F7 7
#define F6 6
#define F5 5
#define F4 4 
#define F3 3 
#define F2 2 
#define F1 1 
#define F0 0

	PROLOGUE
	
	daddiu	$sp, $sp, -160
	sd	$16,   0($sp)
	sd	$17,   8($sp)
	sd	$18,  16($sp)
	sd	$19,  24($sp)
	sd	$20,  32($sp)
	sd	$21,  40($sp)
	sd	$22,  48($sp)
	ST	$f24, 56($sp)
	ST	$f25, 64($sp)
	ST	$f26, 72($sp)
	ST	$f27, 80($sp)
	ST	$f28, 88($sp)
	sd	$23,  96($sp)
	sd	$24, 104($sp)
	sd	$25, 112($sp)
	ST	$f20,120($sp)
	ST	$f21,128($sp)
	ST	$f22,136($sp)
	ST	$f23,144($sp)


	.align	5					#  BACKUP
.L0_N4:							#  Loop N
	ST	ALPHA,152($sp)		#  Backup	ALPHA
	move	MCO,M				#  Backup	M

	move	NCO,N				#  Backup	N
	move	KCO,K				#  Backup	K

	move	AO,A				#  Backup	A_addr
	move	BO,B				#  Backup	B_addr
	
	dsll	LDC,LDC,BASE_SHIFT			#  LDC*8Byte
	dsll	SPANB,KCO,2+BASE_SHIFT			#  SPANB=KC*NR(4)*8Byte=KC*2^5
	
	dsll	SPANA,KCO,1+BASE_SHIFT			#  SPANA = KCO*4mr*8Byte
	dsra	N,NCO,2				#  N=NCO/2

	beq		N,$0,.L0_N2			#  N=0,NCO<4
	dsll	SPANC,LDC,2			#  SPANC=LDC*4

.L0_N4_Lb:
	move	CO1,C				#  Set C	
	dsra	M,MCO,2				#  M=MCO/2
	
	move	A,AO				#  Reset A
	daddu	CO2,CO1,LDC

	daddu	CO3,CO2,LDC
	daddu	PREB,BO,SPANB		#  PreB point next panelB

	daddu	CO4,CO3,LDC
	beqz	M,.L14_M2
	daddu	PREA,AO,SPANA

.L10:
	MTC		$0,t11
	MOV	t21,t11
	gsLQC1(R8,F1,F0,0)			#a0,a1

	MOV	t31,t11
	MOV	t41,t11
	gsLQC1(R9,F9,F8,0)			#b0,b1

	MOV	t12,t11
	MOV	t22,t11
	gsLQC1(R8,F3,F2,1)			#a2,a3
	
	MOV	t32,t11
	MOV	t42,t11
	gsLQC1(R9,F11,F10,1)		#b2,b3

	dsra	K,KCO,2				#  K=KCO/2
	MOV	t13,t11
	
	MOV	t23,t11
	MOV	t33,t11

	MOV	t43,t11
	MOV	t14,t11
	
	MOV	t24,t11
	MOV	t34,t11
	
	MOV	t44,t11
	beqz	K,.L15
	nop

.L11:							#  N=M=K=4
	gsLQC1(R8,F5,F4,2)			#  R8=A
	MADD	t11,t11,a0,b0
	MADD	t21,t21,a1,b0

	gsLQC1(R9,F13,F12,2)		#  R9=B
	MADD	t12,t12,a0,b1
	MADD	t22,t22,a1,b1

	gsLQC1(R8,F7,F6,3)
	MADD	t31,t31,a2,b0
	MADD	t41,t41,a3,b0
	
	gsLQC1(R9,F15,F14,3)
	MADD	t32,t32,a2,b1
	MADD	t42,t42,a3,b1

	FETCH		$0,(PREB)
	MADD	t13,t13,a0,b2
	MADD	t23,t23,a1,b2
	
	MADD	t14,t14,a0,b3
	MADD	t24,t24,a1,b3
	
	FETCH		$0,(PREA)
	MADD	t33,t33,a2,b2
	MADD	t43,t43,a3,b2

	MADD	t34,t34,a2,b3
	MADD	t44,t44,a3,b3
								#load2 comp1
.L12:
	gsLQC1(R8,F1,F0,4)
	MADD	t11,t11,a4,b4
	MADD	t21,t21,a5,b4

	gsLQC1(R9,F9,F8,4)
	MADD	t12,t12,a4,b5
	MADD	t22,t22,a5,b5

	gsLQC1(R8,F3,F2,5)
	MADD	t31,t31,a6,b4
	MADD	t41,t41,a7,b4

	gsLQC1(R9,F11,F10,5)
	MADD	t32,t32,a6,b5
	MADD	t42,t42,a7,b5

	FETCH		$0,4*SIZE(PREB)
	MADD	t13,t13,a4,b6
	MADD	t23,t23,a5,b6

	MADD	t14,t14,a4,b7
	MADD	t24,t24,a5,b7

	FETCH		$0,4*SIZE(PREA)
	MADD	t33,t33,a6,b6
	MADD	t43,t43,a7,b6
	
	MADD	t34,t34,a6,b7
	MADD	t44,t44,a7,b7

.L13:
	gsLQC1(R8,F5,F4,6)	
	MADD	t11,t11,a0,b0
	MADD	t21,t21,a1,b0

	gsLQC1(R9,F13,F12,6)
	MADD	t12,t12,a0,b1
	MADD	t22,t22,a1,b1

	gsLQC1(R8,F7,F6,7)
	MADD	t31,t31,a2,b0
	MADD	t41,t41,a3,b0

	gsLQC1(R9,F15,F14,7)
	MADD	t32,t32,a2,b1
	MADD	t42,t42,a3,b1
	daddu	A,A,16*SIZE				#  A+=4(mr)*4(kr)*8Byte=16*SIZE

	FETCH		$0,8*SIZE(PREB)
	MADD	t13,t13,a0,b2
	MADD	t23,t23,a1,b2
	daddu	B,B,16*SIZE

	MADD	t14,t14,a0,b3
	MADD	t24,t24,a1,b3

	FETCH		$0,8*SIZE(PREA)
	MADD	t33,t33,a2,b2
	MADD	t43,t43,a3,b2

	MADD	t34,t34,a2,b3
	MADD	t44,t44,a3,b3
	
.L14:
	gsLQC1(R8,F1,F0,0)
	MADD	t11,t11,a4,b4
	MADD	t21,t21,a5,b4

	gsLQC1(R9,F9,F8,0)
	MADD	t12,t12,a4,b5
	MADD	t22,t22,a5,b5

	gsLQC1(R8,F3,F2,1)
	MADD	t31,t31,a6,b4
	MADD	t41,t41,a7,b4
	daddiu	K,K,-1

	gsLQC1(R9,F11,F10,1)
	MADD	t32,t32,a6,b5
	MADD	t42,t42,a7,b5

	FETCH		$0,12*SIZE(PREB)
	MADD	t13,t13,a4,b6
	MADD	t23,t23,a5,b6

	MADD	t14,t14,a4,b7
	MADD	t24,t24,a5,b7

	FETCH		$0,12*SIZE(PREA)
	MADD	t33,t33,a6,b6
	MADD	t43,t43,a7,b6
	daddu	PREB,PREB,16*SIZE
	
	MADD	t34,t34,a6,b7
	daddu	PREA,PREA,16*SIZE
	bnez 	K,.L11
	MADD	t44,t44,a7,b7

.L15:							#  N=4 M=4 K=2
	and		K,KCO,2				#  k = KCO&2
	beqz	K,.L18
	nop

.L16:			
	gsLQC1(R8,F5,F4,2)			#  R8=A
	MADD	t11,t11,a0,b0
	MADD	t21,t21,a1,b0

	gsLQC1(R9,F13,F12,2)		#  R9=B
	MADD	t12,t12,a0,b1
	MADD	t22,t22,a1,b1

	gsLQC1(R8,F7,F6,3)
	MADD	t31,t31,a2,b0
	MADD	t41,t41,a3,b0

	gsLQC1(R9,F15,F14,3)
	MADD	t32,t32,a2,b1
	MADD	t42,t42,a3,b1
	daddu	A,A,8*SIZE				#  A+=4(mr)*2(kr)*8Byte=8*SIZE

	FETCH		$0,0(PREB)
	MADD	t13,t13,a0,b2
	MADD	t23,t23,a1,b2
	daddu	B,B,8*SIZE

	MADD	t14,t14,a0,b3
	MADD	t24,t24,a1,b3

	FETCH		$0,0(PREA)
	MADD	t33,t33,a2,b2
	MADD	t43,t43,a3,b2

	MADD	t34,t34,a2,b3
	MADD	t44,t44,a3,b3
								
.L17:
	gsLQC1(R8,F1,F0,0)
	MADD	t11,t11,a4,b4
	MADD	t21,t21,a5,b4

	gsLQC1(R9,F9,F8,0)
	MADD	t12,t12,a4,b5
	MADD	t22,t22,a5,b5

	gsLQC1(R8,F3,F2,1)
	MADD	t31,t31,a6,b4
	MADD	t41,t41,a7,b4

	gsLQC1(R9,F11,F10,1)
	MADD	t32,t32,a6,b5
	MADD	t42,t42,a7,b5

	FETCH		$0,4*SIZE(PREB)
	MADD	t13,t13,a4,b6
	MADD	t23,t23,a5,b6

	MADD	t14,t14,a4,b7
	MADD	t24,t24,a5,b7

	FETCH		$0,4*SIZE(PREA)
	MADD	t33,t33,a6,b6
	MADD	t43,t43,a7,b6
	daddu	PREB,PREB,8*SIZE
	
	MADD	t34,t34,a6,b7
	MADD	t44,t44,a7,b7
	daddu	PREA,PREA,8*SIZE
	
.L18:							#  N=4, M=4, K=1
	and		K,KCO,1
	beqz	K,.L19				#  
	LD	ALPHA,152($sp)		#  Get ALPHA
	
	FETCH		$0,0(PREB)
	MADD	t11,t11,a0,b0
	MADD	t21,t21,a1,b0
	daddu	A,A,4*SIZE				#  A+=4(mr)*1(kr)*8Byte=32

	MADD	t12,t12,a0,b1
	MADD	t22,t22,a1,b1
	daddu	B,B,4*SIZE

	FETCH		$0,0(PREA)
	MADD	t31,t31,a2,b0
	MADD	t41,t41,a3,b0
	daddu	PREB,PREB,4*SIZE

	MADD	t32,t32,a2,b1
	MADD	t42,t42,a3,b1
	daddu	PREA,PREA,4*SIZE

	MADD	t13,t13,a0,b2
	MADD	t23,t23,a1,b2

	MADD	t14,t14,a0,b3
	MADD	t24,t24,a1,b3

	MADD	t33,t33,a2,b2
	MADD	t43,t43,a3,b2

	MADD	t34,t34,a2,b3
	MADD	t44,t44,a3,b3

.L19:							#  Write Back
	LD	c11,0(CO1)			#  Fetch 16 C
	LD	c21,1*SIZE(CO1)			
	LD	c31,2*SIZE(CO1)
	LD	c41,3*SIZE(CO1)

	LD	c12,0(CO2)
	MADD	t11,c11,t11,ALPHA
	LD	c22,1*SIZE(CO2)
	MADD	t21,c21,t21,ALPHA
	LD	c32,2*SIZE(CO2)
	MADD	t31,c31,t31,ALPHA
	LD	c42,3*SIZE(CO2)
	MADD	t41,c41,t41,ALPHA

	LD	c13,0(CO3)
	MADD	t12,c12,t12,ALPHA
	LD	c23,1*SIZE(CO3)
	MADD	t22,c22,t22,ALPHA
	LD	c33,2*SIZE(CO3)
	MADD	t32,c32,t32,ALPHA
	LD	c43,3*SIZE(CO3)
	MADD	t42,c42,t42,ALPHA

	LD	c14,0(CO4)
	MADD	t13,c13,t13,ALPHA
	LD	c24,1*SIZE(CO4)
	MADD	t23,c23,t23,ALPHA
	LD	c34,2*SIZE(CO4)
	MADD	t33,c33,t33,ALPHA
	LD	c44,3*SIZE(CO4)
	MADD	t43,c43,t43,ALPHA

	ST	t11,0(CO1)
	MADD	t14,c14,t14,ALPHA
	ST	t21,1*SIZE(CO1)
	MADD	t24,c24,t24,ALPHA
	ST	t31,2*SIZE(CO1)
	MADD	t34,c34,t34,ALPHA
	ST	t41,3*SIZE(CO1)
	MADD	t44,c44,t44,ALPHA
	daddiu	M,M,-1				#  M--

	ST	t12,0(CO2)
	ST	t22,1*SIZE(CO2)
	ST	t32,2*SIZE(CO2)
	ST	t42,3*SIZE(CO2)

	ST	t13,0(CO3)
	ST	t23,1*SIZE(CO3)
	ST	t33,2*SIZE(CO3)
	ST	t43,3*SIZE(CO3)

	FETCH	$0,4*SIZE(CO1)
	FETCH	$0,4*SIZE(CO2)
	FETCH	$0,4*SIZE(CO3)
	FETCH	$0,4*SIZE(CO4)

	FETCH	$0,8*SIZE(CO1)
	FETCH	$0,8*SIZE(CO2)
	FETCH	$0,8*SIZE(CO3)
	FETCH	$0,8*SIZE(CO4)

	ST	t14,0(CO4)
	daddu	CO1,CO1,4*SIZE			#  COx += 4*8Byte
	ST	t24,1*SIZE(CO4)
	daddu	CO2,CO2,4*SIZE
	ST	t34,2*SIZE(CO4)
	daddu	CO3,CO3,4*SIZE
	ST	t44,3*SIZE(CO4)
	move	B,BO				#  Reset B
	daddu	PREB,BO,SPANB
	bnez	M,.L10				#  M!=0
	daddu	CO4,CO4,4*SIZE



.L14_M2:
	and		M,MCO,2				#  Remainder M = 2
	beqz	M,.L14_M1			
	nop

.L20:
	MTC		$0,t11
	MOV	t21,t11
	gsLQC1(R8,F1,F0,0)			#a0,a1

	MOV	t12,t11
	MOV	t22,t11
	gsLQC1(R9,F9,F8,0)			#b0,b1
	
	dsra	K,KCO,2				#  K=KCO/2
	MOV	t13,t11
	gsLQC1(R9,F11,F10,1)		#b2,b3
	
	MOV	t23,t11
	MOV	t14,t11
	
	MOV	t24,t11
	beqz	K,.L25
	nop

.L21:							#  N=4 m=2,=K=4
	gsLQC1(R8,F5,F4,1)			#  R8=A
	MADD	t11,t11,a0,b0
	MADD	t21,t21,a1,b0

	gsLQC1(R9,F13,F12,2)		#  R9=B
	MADD	t12,t12,a0,b1
	MADD	t22,t22,a1,b1

	gsLQC1(R9,F15,F14,3)
	MADD	t13,t13,a0,b2
	MADD	t23,t23,a1,b2
	
	gsLQC1(R8,F3,F2,2)
	MADD	t14,t14,a0,b3
	MADD	t24,t24,a1,b3
	
	gsLQC1(R9,F9,F8,4)
	MADD	t11,t11,a4,b4
	MADD	t21,t21,a5,b4

	gsLQC1(R9,F11,F10,5)
	MADD	t12,t12,a4,b5
	MADD	t22,t22,a5,b5

	gsLQC1(R8,F7,F6,3)	
	MADD	t13,t13,a4,b6
	MADD	t23,t23,a5,b6

	MADD	t14,t14,a4,b7
	MADD	t24,t24,a5,b7

	gsLQC1(R9,F13,F12,6)
	MADD	t11,t11,a2,b0
	MADD	t21,t21,a3,b0
	daddu	A,A,8*SIZE				#  A+=2(mr)*4(kr)*8Byte=8*SIZE

	gsLQC1(R9,F15,F14,7)
	MADD	t12,t12,a2,b1
	MADD	t22,t22,a3,b1
	daddiu	K,K,-1

	gsLQC1(R8,F1,F0,0)
	MADD	t13,t13,a2,b2
	MADD	t23,t23,a3,b2
	daddu	B,B,16*SIZE				#  B+=4(nr)*4(kr)*8Byte=16*SIZE

	MADD	t14,t14,a2,b3
	MADD	t24,t24,a3,b3

	gsLQC1(R9,F9,F8,0)
	MADD	t11,t11,a6,b4
	MADD	t21,t21,a7,b4

	gsLQC1(R9,F11,F10,1)
	MADD	t12,t12,a6,b5
	MADD	t22,t22,a7,b5

	MADD	t13,t13,a6,b6
	MADD	t23,t23,a7,b6

	MADD	t14,t14,a6,b7
	bnez 	K,.L21
	MADD	t24,t24,a7,b7

.L25:							#  N=4 M=2 K=2
	and		K,KCO,2				#  k = KCO&2
	beqz	K,.L28
	nop

.L26:			
	gsLQC1(R8,F5,F4,1)			#  R8=A
	MADD	t11,t11,a0,b0
	MADD	t21,t21,a1,b0

	gsLQC1(R9,F13,F12,2)		#  R9=B
	MADD	t12,t12,a0,b1
	MADD	t22,t22,a1,b1
	daddu	A,A,4*SIZE				#  A+=2(mr)*2(kr)*8Byte=32

	gsLQC1(R9,F15,F14,3)
	MADD	t13,t13,a0,b2
	MADD	t23,t23,a1,b2
	daddu	B,B,8*SIZE
	
	MADD	t14,t14,a0,b3
	MADD	t24,t24,a1,b3

.L27:
	gsLQC1(R8,F1,F0,0)
	MADD	t11,t11,a4,b4
	MADD	t21,t21,a5,b4

	gsLQC1(R9,F9,F8,0)
	MADD	t12,t12,a4,b5
	MADD	t22,t22,a5,b5

	gsLQC1(R9,F11,F10,1)
	MADD	t13,t13,a4,b6
	MADD	t23,t23,a5,b6

	MADD	t14,t14,a4,b7
	MADD	t24,t24,a5,b7
	
.L28:							#  N=4, M=2, K=1
	and		K,KCO,1
	beqz	K,.L29				#  
	LD	ALPHA,152($sp)		#  Get ALPHA
	
	MADD	t11,t11,a0,b0
	MADD	t21,t21,a1,b0
	daddu	A,A,2*SIZE				#  A+=2(mr)*1(kr)*8Byte=16
	daddu	B,B,4*SIZE

	MADD	t12,t12,a0,b1
	MADD	t22,t22,a1,b1

	MADD	t13,t13,a0,b2
	MADD	t23,t23,a1,b2

	MADD	t14,t14,a0,b3
	MADD	t24,t24,a1,b3

.L29:							#  Write Back
	LD	c11,0(CO1)			#  Fetch 16 C
	LD	c21,1*SIZE(CO1)			

	LD	c12,0(CO2)
	LD	c22,1*SIZE(CO2)
	
	LD	c13,0(CO3)
	MADD	t11,c11,t11,ALPHA
	LD	c23,1*SIZE(CO3)
	MADD	t21,c21,t21,ALPHA

	LD	c14,0(CO4)
	MADD	t12,c12,t12,ALPHA
	LD	c24,1*SIZE(CO4)
	MADD	t22,c22,t22,ALPHA

	ST	t11,0(CO1)
	MADD	t13,c13,t13,ALPHA
	ST	t21,1*SIZE(CO1)
	MADD	t23,c23,t23,ALPHA

	ST	t12,0(CO2)
	MADD	t14,c14,t14,ALPHA
	ST	t22,1*SIZE(CO2)
	MADD	t24,c24,t24,ALPHA

	ST	t13,0(CO3)
	move	B,BO				#  Reset B
	ST	t23,1*SIZE(CO3)
	daddu	CO1,CO1,2*SIZE			#  COx += 2*8Byte

	FETCH	$0,0(CO1)
	FETCH	$0,2*SIZE(CO2)
	FETCH	$0,2*SIZE(CO3)
	FETCH	$0,2*SIZE(CO4)

	ST	t14,0(CO4)
	daddu	CO2,CO2,2*SIZE
	ST	t24,1*SIZE(CO4)
	daddu	CO3,CO3,2*SIZE
	daddu	CO4,CO4,2*SIZE



.L14_M1:
	and		M,MCO,1				#  Remainder M = 1
	beqz	M,.L0_N4_Loop		#  M = 0, finishing one panel B
	nop

.L30:
	gsLQC1(R8,F1,F0,0)
	dsra	K,KCO,2				#  K=KCO/2
	gsLQC1(R9,F9,F8,0)			#b0,b1
	MTC		$0,t11
	gsLQC1(R9,F11,F10,1)		#b2,b3
	MOV	t12,t11
	MOV	t13,t11
	beqz	K,.L35
	MOV	t14,t11

.L31:							#  N=4 m=1,=K=4
	gsLQC1(R8,F3,F2,1)	
	gsLQC1(R9,F13,F12,2)		#  R9=B
	MADD	t11,t11,a0,b0
	MADD	t12,t12,a0,b1
	
	gsLQC1(R9,F15,F14,3)
	MADD	t13,t13,a0,b2
	MADD	t14,t14,a0,b3
	
	gsLQC1(R9,F9,F8,4)
	MADD	t11,t11,a1,b4
	MADD	t12,t12,a1,b5
	
	gsLQC1(R9,F11,F10,5)
	MADD	t13,t13,a1,b6
	MADD	t14,t14,a1,b7
	daddiu	K,K,-1
	
	gsLQC1(R9,F13,F12,6)
	MADD	t11,t11,a2,b0
	MADD	t12,t12,a2,b1
	daddu	A,A,4*SIZE				#  A+=1(mr)*4(kr)*8Byte=8*SIZE
	
	gsLQC1(R9,F15,F14,7)
	MADD	t13,t13,a2,b2
	MADD	t14,t14,a2,b3
	daddu	B,B,16*SIZE				#  B+=4(nr)*4(kr)*8Byte=16*SIZE

	gsLQC1(R8,F1,F0,0)	
	gsLQC1(R9,F9,F8,0)
	MADD	t11,t11,a3,b4
	MADD	t12,t12,a3,b5
	
	gsLQC1(R9,F11,F10,1)
	MADD	t13,t13,a3,b6
	bnez 	K,.L31
	MADD	t14,t14,a3,b7

.L35:							#  N=4 M=1 K=2
	and		K,KCO,2				#  k = KCO&2
	beqz	K,.L38
	nop

.L36:			
	gsLQC1(R9,F13,F12,2)		#  R9=B
	MADD	t11,t11,a0,b0
	MADD	t12,t12,a0,b1
	daddu	A,A,2*SIZE				#  A+=1(mr)*2(kr)*8Byte=32
	
	gsLQC1(R9,F15,F14,3)
	MADD	t13,t13,a0,b2
	MADD	t14,t14,a0,b3
	daddu	B,B,8*SIZE


.L37:
	LD	a0,0(A)
	
	gsLQC1(R9,F9,F8,0)
	MADD	t11,t11,a1,b4
	MADD	t12,t12,a1,b5
	
	gsLQC1(R9,F11,F10,1)
	MADD	t13,t13,a1,b6
	MADD	t14,t14,a1,b7
	
.L38:							#  N=4, M=1, K=1
	and		K,KCO,1
	beqz	K,.L39				#  
	LD	ALPHA,152($sp)		#  Get ALPHA
	
	MADD	t11,t11,a0,b0
	MADD	t12,t12,a0,b1
	daddu	A,A,1*SIZE				#  A+=1(mr)*1(kr)*8Byte=16
	daddu	B,B,4*SIZE
	
	MADD	t13,t13,a0,b2
	MADD	t14,t14,a0,b3

.L39:							#  Write Back
	LD	c11,0(CO1)			#  Fetch 16 C
	LD	c12,0(CO2)
	LD	c13,0(CO3)
	LD	c14,0(CO4)
	
	MADD	t11,c11,t11,ALPHA
	MADD	t12,c12,t12,ALPHA
	MADD	t13,c13,t13,ALPHA
	MADD	t14,c14,t14,ALPHA

	ST	t11,0(CO1)
	ST	t12,0(CO2)
	ST	t13,0(CO3)
	ST	t14,0(CO4)


.L0_N4_Loop:
	daddu	BO,BO,SPANB			#  BO point to next panel B
	daddiu	N,N,-1				#  N--
	daddu	C,C,SPANC			#  C pointe to next panel C
	bnez	N,.L0_N4_Lb			#  N!=0
	move	B,BO				#  Set B



	.align	5					
.L0_N2:
	and		N,NCO,2				#  Remainder N = 2
	beqz	N,.L0_N1			#  N=0,NCO<2
	dsll	SPANC,LDC,1			#  SPANC=LDC*2

.L0_N2_Lb:
	move	CO1,C				#  Set C	
	dsra	M,MCO,2				#  M=MCO/2
	
	dsll	SPANB,KCO,1+BASE_SHIFT			#  SPANB=KC*NR(2)*8Byte=KC*16=KC*2^4
	move	A,AO				#  Reset A

	daddu	CO2,CO1,LDC
	beqz	M,.L12_M2
	daddu	PREA,AO,SPANA

.L40:						
	MTC		$0,t11
	MOV	t21,t11
	gsLQC1(R8,F1,F0,0)			#a0,a1

	MOV	t31,t11
	MOV	t41,t11
	gsLQC1(R9,F9,F8,0)			#b0,b1

	dsra	K,KCO,2				#  K=KCO/2
	MOV	t12,t11
	gsLQC1(R8,F3,F2,1)			#a2,a3
	
	MOV	t22,t11
	MOV	t32,t11
	
	MOV	t42,t11
	beqz	K,.L45
	nop

.L41:							#  N=2,M=K=4
	gsLQC1(R8,F5,F4,2)			#  R8=A
	MADD	t11,t11,a0,b0
	MADD	t21,t21,a1,b0

	gsLQC1(R9,F13,F12,1)		#  R9=B
	MADD	t12,t12,a0,b1
	MADD	t22,t22,a1,b1

	gsLQC1(R8,F7,F6,3)
	MADD	t31,t31,a2,b0
	MADD	t41,t41,a3,b0
	
	FETCH		$0,(PREA)
	MADD	t32,t32,a2,b1
	MADD	t42,t42,a3,b1

.L42:
	gsLQC1(R8,F1,F0,4)
	MADD	t11,t11,a4,b4
	MADD	t21,t21,a5,b4

	gsLQC1(R9,F11,F10,2)
	MADD	t12,t12,a4,b5
	MADD	t22,t22,a5,b5

	gsLQC1(R8,F3,F2,5)
	MADD	t31,t31,a6,b4
	MADD	t41,t41,a7,b4

	FETCH		$0,4*SIZE(PREA)
	MADD	t32,t32,a6,b5
	MADD	t42,t42,a7,b5

.L43:
	gsLQC1(R8,F5,F4,6)	
	MADD	t11,t11,a0,b2
	MADD	t21,t21,a1,b2

	gsLQC1(R9,F15,F14,3)
	MADD	t12,t12,a0,b3
	MADD	t22,t22,a1,b3

	gsLQC1(R8,F7,F6,7)
	MADD	t31,t31,a2,b2
	MADD	t41,t41,a3,b2
	daddu	B,B,8*SIZE				#  B+=2(nr)*4(kr)*8Byte=8*SIZE

	FETCH		$0,8*SIZE(PREA)
	MADD	t32,t32,a2,b3
	MADD	t42,t42,a3,b3
	daddu	A,A,16*SIZE				#  A+=4(mr)*4(kr)*8Byte=16*SIZE

.L44:
	gsLQC1(R8,F1,F0,0)
	MADD	t11,t11,a4,b6
	MADD	t21,t21,a5,b6
	daddiu	K,K,-1

	gsLQC1(R9,F9,F8,0)
	MADD	t12,t12,a4,b7
	MADD	t22,t22,a5,b7
	daddu	PREA,PREA,16*SIZE

	gsLQC1(R8,F3,F2,1)
	MADD	t31,t31,a6,b6
	MADD	t41,t41,a7,b6

	FETCH		$0,-4*SIZE(PREA)
	MADD	t32,t32,a6,b7
	bnez 	K,.L41
	MADD	t42,t42,a7,b7


.L45:							#  N=2 M=4 K=2
	and		K,KCO,2				#  k = KCO&2
	beqz	K,.L48
	nop

.L46:			
	gsLQC1(R8,F5,F4,2)			#  R8=A
	MADD	t11,t11,a0,b0
	MADD	t21,t21,a1,b0

	gsLQC1(R9,F13,F12,1)		#  R9=B
	MADD	t12,t12,a0,b1
	MADD	t22,t22,a1,b1

	gsLQC1(R8,F7,F6,3)
	MADD	t31,t31,a2,b0
	MADD	t41,t41,a3,b0
	daddu	B,B,4*SIZE				#  B+=2(nr)*2(kr)*8Byte=32

	FETCH		$0,0(PREA)
	MADD	t32,t32,a2,b1
	MADD	t42,t42,a3,b1
	daddu	A,A,8*SIZE				#  A+=4(mr)*2(kr)*8Byte=8*SIZE

.L47:
	gsLQC1(R8,F1,F0,0)
	MADD	t11,t11,a4,b4
	MADD	t21,t21,a5,b4

	gsLQC1(R9,F9,F8,0)
	MADD	t12,t12,a4,b5
	MADD	t22,t22,a5,b5

	gsLQC1(R8,F3,F2,1)
	MADD	t31,t31,a6,b4
	MADD	t41,t41,a7,b4

	FETCH		$0,4*SIZE(PREA)
	MADD	t32,t32,a6,b5
	MADD	t42,t42,a7,b5
	daddu	PREA,PREA,8*SIZE

	
.L48:							#  N=2, M=4, K=1
	and		K,KCO,1
	beqz	K,.L49				#  
	LD	ALPHA,152($sp)		#  Get ALPHA
	
	FETCH		$0,0(PREA)
	MADD	t11,t11,a0,b0
	MADD	t21,t21,a1,b0
	daddu	A,A,4*SIZE				#  A+=4(mr)*1(kr)*8Byte=32

	MADD	t12,t12,a0,b1
	MADD	t22,t22,a1,b1
	daddu	B,B,2*SIZE
	daddu	PREA,PREA,4*SIZE

	MADD	t31,t31,a2,b0
	MADD	t41,t41,a3,b0

	MADD	t32,t32,a2,b1
	MADD	t42,t42,a3,b1

.L49:							#  Write Back
	LD	c11,0(CO1)			#  Fetch 16 C
	LD	c21,1*SIZE(CO1)			
	LD	c31,2*SIZE(CO1)
	LD	c41,3*SIZE(CO1)

	LD	c12,0(CO2)
	MADD	t11,c11,t11,ALPHA
	LD	c22,1*SIZE(CO2)
	MADD	t21,c21,t21,ALPHA
	LD	c32,2*SIZE(CO2)
	MADD	t31,c31,t31,ALPHA
	LD	c42,3*SIZE(CO2)
	MADD	t41,c41,t41,ALPHA

	ST	t11,0(CO1)
	MADD	t12,c12,t12,ALPHA
	ST	t21,1*SIZE(CO1)
	MADD	t22,c22,t22,ALPHA
	ST	t31,2*SIZE(CO1)
	MADD	t32,c32,t32,ALPHA
	ST	t41,3*SIZE(CO1)
	MADD	t42,c42,t42,ALPHA
	daddiu	M,M,-1				#  M--

	ST	t12,0(CO2)
	ST	t22,1*SIZE(CO2)
	ST	t32,2*SIZE(CO2)
	ST	t42,3*SIZE(CO2)

	FETCH	$0,4*SIZE(CO1)
	FETCH	$0,4*SIZE(CO2)

	FETCH	$0,8*SIZE(CO1)
	FETCH	$0,8*SIZE(CO2)

	daddu	CO1,CO1,4*SIZE			#  COx += 4*8Byte
	daddu	CO2,CO2,4*SIZE
	bnez	M,.L40				#  M!=0
	move	B,BO				#  Reset B


.L12_M2:
	and		M,MCO,2				#  Remainder M = 2
	beqz	M,.L12_M1			
	nop

.L50:
	dsra	K,KCO,2				#  K=KCO/2
	MTC		$0,t11
	gsLQC1(R8,F1,F0,0)			#a0,a1

	MOV	t21,t11
	MOV	t12,t11
	gsLQC1(R9,F9,F8,0)			#b0,b1
	
	MOV	t22,t11
	beqz	K,.L55
	nop

.L51:							#  N=2 m=2,=K=4
	gsLQC1(R8,F5,F4,1)			#  R8=A
	MADD	t11,t11,a0,b0
	MADD	t21,t21,a1,b0

	gsLQC1(R9,F13,F12,1)		#  R9=B
	MADD	t12,t12,a0,b1
	MADD	t22,t22,a1,b1

	gsLQC1(R8,F3,F2,2)
	MADD	t11,t11,a4,b4
	MADD	t21,t21,a5,b4

	gsLQC1(R9,F11,F10,2)
	MADD	t12,t12,a4,b5
	MADD	t22,t22,a5,b5
	daddiu	K,K,-1

	gsLQC1(R8,F7,F6,3)	
	MADD	t11,t11,a2,b2
	MADD	t21,t21,a3,b2
	daddu	A,A,8*SIZE				#  A+=2(mr)*4(kr)*8Byte=8*SIZE

	gsLQC1(R9,F15,F14,3)
	MADD	t12,t12,a2,b3
	MADD	t22,t22,a3,b3
	daddu	B,B,8*SIZE				#  B+=2(nr)*4(kr)*8Byte=16*SIZE

	gsLQC1(R8,F1,F0,0)
	MADD	t11,t11,a6,b6
	MADD	t21,t21,a7,b6

	gsLQC1(R9,F9,F8,0)
	MADD	t12,t12,a6,b7
	bnez 	K,.L51
	MADD	t22,t22,a7,b7

.L55:							#  N=2 M=2 K=2
	and		K,KCO,2				#  k = KCO&2
	beqz	K,.L58
	nop

.L56:			
	gsLQC1(R8,F5,F4,1)			#  R8=A
	MADD	t11,t11,a0,b0
	MADD	t21,t21,a1,b0
	daddu	A,A,4*SIZE				#  A+=2(mr)*2(kr)*8Byte=32

	gsLQC1(R9,F13,F12,1)		#  R9=B
	MADD	t12,t12,a0,b1
	MADD	t22,t22,a1,b1
	daddu	B,B,4*SIZE

.L57:
	gsLQC1(R8,F1,F0,0)
	MADD	t11,t11,a4,b4
	MADD	t21,t21,a5,b4

	gsLQC1(R9,F9,F8,0)
	MADD	t12,t12,a4,b5
	MADD	t22,t22,a5,b5

	
.L58:							#  N=2, M=2, K=1
	and		K,KCO,1
	beqz	K,.L59				#  
	LD	ALPHA,152($sp)		#  Get ALPHA
	
	MADD	t11,t11,a0,b0
	MADD	t21,t21,a1,b0
	daddu	A,A,2*SIZE				#  A+=2(mr)*1(kr)*8Byte=16
	daddu	B,B,2*SIZE

	MADD	t12,t12,a0,b1
	MADD	t22,t22,a1,b1


.L59:							#  Write Back
	LD	c11,0(CO1)			#  Fetch 16 C
	LD	c21,1*SIZE(CO1)			
	LD	c12,0(CO2)
	LD	c22,1*SIZE(CO2)
	
	MADD	t11,c11,t11,ALPHA
	MADD	t21,c21,t21,ALPHA
	MADD	t12,c12,t12,ALPHA
	MADD	t22,c22,t22,ALPHA

	ST	t11,0(CO1)
	ST	t21,1*SIZE(CO1)
	ST	t12,0(CO2)
	move	B,BO				#  Reset B
	ST	t22,1*SIZE(CO2)

	daddu	CO1,CO1,2*SIZE			#  COx += 2*8Byte
	daddu	CO2,CO2,2*SIZE

	FETCH	$0,0(CO1)
	FETCH	$0,0(CO2)


.L12_M1:
	and		M,MCO,1				#  Remainder M = 1
	beqz	M,.L0_N2_Loop		#  M = 0, finishing one panel B
	nop

.L60:
	dsra	K,KCO,2				#  K=KCO/2
	MTC		$0,t11
	gsLQC1(R8,F4,F0,0)

	MOV	t21,t11
	MOV	t12,t11
	gsLQC1(R9,F9,F8,0)			#b0,b1
	
	MOV	t22,t11
	beqz	K,.L65
	nop

.L61:							#  N=2 m=1,=K=4
	gsLQC1(R9,F13,F12,1)		#  R9=B
	MADD	t11,t11,a0,b0
	MADD	t12,t12,a0,b1

	gsLQC1(R9,F11,F10,2)
	MADD	t11,t11,a4,b4
	MADD	t12,t12,a4,b5
	daddiu	K,K,-1

	gsLQC1(R8,F6,F2,1)
	MADD	t11,t11,a2,b2
	
	gsLQC1(R9,F15,F14,3)
	MADD	t12,t12,a2,b3
	daddu	A,A,4*SIZE				#  A+=1(mr)*4(kr)*8Byte=32

	gsLQC1(R8,F4,F0,0)
	daddu	B,B,8*SIZE				#  B+=2(nr)*4(kr)*8Byte=8*SIZE
	
	gsLQC1(R9,F9,F8,0)
	MADD	t11,t11,a6,b6
	bnez 	K,.L61
	MADD	t12,t12,a6,b7

.L65:							#  N=2 M=1 K=2
	and		K,KCO,2				#  k = KCO&2
	beqz	K,.L68
	nop

.L66:			
	gsLQC1(R9,F13,F12,1)		#  R9=B
	MADD	t11,t11,a0,b0
	MADD	t12,t12,a0,b1
	daddu	A,A,2*SIZE				#  A+=1(mr)*2(kr)*8Byte=16
	daddu	B,B,4*SIZE

.L67:
	LD	a0,0(A)
	gsLQC1(R9,F9,F8,0)
	MADD	t11,t11,a4,b4
	MADD	t12,t12,a4,b5

	
.L68:							#  N=2, M=1, K=1
	and		K,KCO,1
	beqz	K,.L69				#  
	LD	ALPHA,152($sp)		#  Get ALPHA
	
	MADD	t11,t11,a0,b0
	MADD	t12,t12,a0,b1
	daddu	A,A,1*SIZE				#  A+=1(mr)*1(kr)*8Byte=16
	daddu	B,B,2*SIZE


.L69:							#  Write Back
	LD	c11,0(CO1)			#  Fetch 16 C
	LD	c12,0(CO2)
	
	MADD	t11,c11,t11,ALPHA
	MADD	t12,c12,t12,ALPHA

	ST	t11,0(CO1)
	ST	t12,0(CO2)
	move	B,BO				#  Reset B

	daddu	CO1,CO1,1*SIZE			#  COx += 2*8Byte
	daddu	CO2,CO2,1*SIZE

	FETCH	$0,0(CO1)
	FETCH	$0,0(CO2)


.L0_N2_Loop:
	daddu	BO,BO,SPANB			#  BO+=KC*2N
	move	B,BO				#  Set B
	daddu	C,C,SPANC			#  C+=LDC*2



	.align	5					
.L0_N1:
	and		N,NCO,1				#  Remainder N = 1
	beqz	N,.L999				#  N=0,NCO<1
	nop

	move	CO1,C				#  Set C	
	dsra	M,MCO,2				#  M=MCO/2
	
	move	A,AO				#  Reset A
	beqz	M,.L11_M2
	daddu	PREA,AO,SPANA


.L70:						
	dsra	K,KCO,2				#  K=KCO/2
	gsLQC1(R9,F12,F8,0)
	MTC		$0,t11
	gsLQC1(R8,F1,F0,0)			#a0,a1
	MOV	t21,t11
	gsLQC1(R8,F3,F2,1)			#a2,a3
	MOV	t31,t11
	beqz	K,.L75
	MOV	t41,t11

.L71:							#  N=1,M=K=4
	gsLQC1(R8,F5,F4,2)			#  R8=A
	gsLQC1(R8,F7,F6,3)
	MADD	t11,t11,a0,b0
	MADD	t21,t21,a1,b0

	FETCH		$0,(PREA)
	MADD	t31,t31,a2,b0
	MADD	t41,t41,a3,b0

.L72:
	gsLQC1(R9,F14,F10,1)
	gsLQC1(R8,F1,F0,4)
	gsLQC1(R8,F3,F2,5)
	MADD	t11,t11,a4,b4
	MADD	t21,t21,a5,b4

	FETCH		$0,4*SIZE(PREA)
	MADD	t31,t31,a6,b4
	MADD	t41,t41,a7,b4

.L73:
	gsLQC1(R8,F5,F4,6)	
	gsLQC1(R8,F7,F6,7)
	MADD	t11,t11,a0,b2
	MADD	t21,t21,a1,b2
	daddu	B,B,4*SIZE				#  B+=1(nr)*4(kr)*8Byte=32

	FETCH		$0,8*SIZE(PREA)
	MADD	t31,t31,a2,b2
	MADD	t41,t41,a3,b2
	daddu	A,A,16*SIZE				#  A+=4(mr)*4(kr)*8Byte=16*SIZE

.L74:
	gsLQC1(R9,F12,F8,0)
	gsLQC1(R8,F1,F0,0)
	daddu	PREA,PREA,16*SIZE
	gsLQC1(R8,F3,F2,1)
	MADD	t11,t11,a4,b6
	MADD	t21,t21,a5,b6
	daddiu	K,K,-1

	FETCH		$0,-32(PREA)
	MADD	t31,t31,a6,b6
	bnez 	K,.L71
	MADD	t41,t41,a7,b6


.L75:							#  N=2 M=4 K=2
	and		K,KCO,2				#  k = KCO&2
	beqz	K,.L78
	nop

.L76:			
	gsLQC1(R8,F5,F4,2)			#  R8=A
	gsLQC1(R8,F7,F6,3)
	MADD	t11,t11,a0,b0
	MADD	t21,t21,a1,b0
	daddu	B,B,2*SIZE				#  B+=1(nr)*2(kr)*8Byte=32

	FETCH		$0,0(PREA)
	MADD	t31,t31,a2,b0
	MADD	t41,t41,a3,b0
	daddu	A,A,8*SIZE				#  A+=4(mr)*2(kr)*8Byte=8*SIZE

.L77:
	LD	b0,0(B)
	gsLQC1(R8,F1,F0,0)
	gsLQC1(R8,F3,F2,1)
	MADD	t11,t11,a4,b4
	MADD	t21,t21,a5,b4

	FETCH		$0,4*SIZE(PREA)
	MADD	t31,t31,a6,b4
	MADD	t41,t41,a7,b4
	daddu	PREA,PREA,8*SIZE


	
.L78:							#  N=2, M=4, K=1
	and		K,KCO,1
	beqz	K,.L79				#  
	LD	ALPHA,152($sp)		#  Get ALPHA
	
	FETCH		$0,0(PREA)
	MADD	t11,t11,a0,b0
	MADD	t21,t21,a1,b0
	daddu	A,A,4*SIZE				#  A+=4(mr)*1(kr)*8Byte=32

	MADD	t31,t31,a2,b0
	MADD	t41,t41,a3,b0
	daddu	B,B,1*SIZE
	daddu	PREA,PREA,4*SIZE


.L79:							#  Write Back
	LD	c11,0(CO1)			#  Fetch 16 C
	LD	c21,1*SIZE(CO1)			
	LD	c31,2*SIZE(CO1)
	LD	c41,3*SIZE(CO1)

	MADD	t11,c11,t11,ALPHA
	MADD	t21,c21,t21,ALPHA
	MADD	t31,c31,t31,ALPHA
	MADD	t41,c41,t41,ALPHA

	ST	t11,0(CO1)
	ST	t21,1*SIZE(CO1)
	ST	t31,2*SIZE(CO1)
	ST	t41,3*SIZE(CO1)
	daddiu	M,M,-1				#  M--

	FETCH	$0,4*SIZE(CO1)

	daddu	CO1,CO1,4*SIZE			#  COx += 4*8Byte
	bnez	M,.L70				#  M!=0
	move	B,BO				#  Reset B



.L11_M2:
	and		M,MCO,2				#  Remainder M = 2
	beqz	M,.L11_M1			
	nop

.L80:						
	dsra	K,KCO,2				#  K=KCO/2
	gsLQC1(R9,F12,F8,0)
	MTC		$0,t11
	gsLQC1(R8,F1,F0,0)			#a0,a1
	MOV		t21,t11
	beqz	K,.L85
	nop

.L81:							#  N=1,M=2,K=4
	gsLQC1(R8,F5,F4,1)			#  R8=A
	MADD	t11,t11,a0,b0
	MADD	t21,t21,a1,b0

	gsLQC1(R8,F3,F2,2)
	MADD	t11,t11,a4,b4
	MADD	t21,t21,a5,b4

	gsLQC1(R9,F14,F10,1)
	daddu	B,B,4*SIZE				#  B+=1(nr)*4(kr)*8Byte=32
	
	gsLQC1(R8,F7,F6,3)
	MADD	t11,t11,a2,b2
	MADD	t21,t21,a3,b2
	daddu	A,A,8*SIZE				#  A+=2(mr)*4(kr)*8Byte=8*SIZE

	gsLQC1(R9,F12,F8,0)
	daddiu	K,K,-1
	
	gsLQC1(R8,F1,F0,0)
	MADD	t11,t11,a6,b6
	bnez 	K,.L81
	MADD	t21,t21,a7,b6


.L85:							#  N=2 M=4 K=2
	and		K,KCO,2				#  k = KCO&2
	beqz	K,.L88
	nop

.L86:			
	gsLQC1(R8,F5,F4,1)			#  R8=A
	MADD	t11,t11,a0,b0
	MADD	t21,t21,a1,b0
	daddu	B,B,2*SIZE				#  B+=1(nr)*2(kr)*8Byte=16

	LD	b0,0(B)
	daddu	A,A,4*SIZE				#  A+=2(mr)*2(kr)*8Byte=32
	
	gsLQC1(R8,F1,F0,0)
	MADD	t11,t11,a4,b4
	MADD	t21,t21,a5,b4

	
.L88:							#  N=2, M=4, K=1
	and		K,KCO,1
	beqz	K,.L89				#  
	LD	ALPHA,152($sp)		#  Get ALPHA
	
	MADD	t11,t11,a0,b0
	MADD	t21,t21,a1,b0
	daddu	A,A,2*SIZE				#  A+=2(mr)*1(kr)*8Byte=16
	daddu	B,B,1*SIZE


.L89:							#  Write Back
	LD	c11,0(CO1)			#  Fetch 16 C
	LD	c21,1*SIZE(CO1)			

	MADD	t11,c11,t11,ALPHA
	MADD	t21,c21,t21,ALPHA

	ST	t11,0(CO1)
	ST	t21,1*SIZE(CO1)

	FETCH	$0,2*SIZE(CO1)
	
	daddu	CO1,CO1,2*SIZE			#  COx += 2*8Byte
	move	B,BO				#  Reset B


.L11_M1:
	and		M,MCO,1				#  Remainder M = 1
	beqz	M,.L999				#  M = 0, End
	nop

.L90:						
	dsra	K,KCO,2				#  K=KCO/2
	gsLQC1(R8,F4,F0,0)
	gsLQC1(R9,F12,F8,0)
	beqz	K,.L95
	MTC		$0,t11

.L91:							#  N=1,M=1,K=4
	gsLQC1(R8,F6,F2,1)
	MADD	t11,t11,a0,b0
	gsLQC1(R9,F14,F10,1)
	MADD	t11,t11,a4,b4
	daddu	A,A,4*SIZE				#  A+=1(mr)*4(kr)*8Byte=32


	gsLQC1(R8,F4,F0,0)
	MADD	t11,t11,a2,b2
	daddu	B,B,4*SIZE				#  B+=1(nr)*4(kr)*8Byte=32
	
	gsLQC1(R9,F12,F8,0)
	MADD	t11,t11,a6,b6
	daddiu	K,K,-1
	bnez 	K,.L91
	nop

.L95:							#  N=2 M=4 K=2
	and		K,KCO,2				#  k = KCO&2
	beqz	K,.L98
	nop

.L96:			
	MADD	t11,t11,a0,b0
	MADD	t11,t11,a4,b4
	daddu	B,B,2*SIZE				#  B+=1(nr)*2(kr)*8Byte=16
	daddu	A,A,2*SIZE				#  A+=1(mr)*2(kr)*8Byte=32

	LD	b0,0(B)
	LD	a0,0(A)

	
.L98:							#  N=2, M=4, K=1
	and		K,KCO,1
	beqz	K,.L99				#  
	LD	ALPHA,152($sp)		#  Get ALPHA
	MADD	t11,t11,a0,b0


.L99:							#  Write Back
	LD	c11,0(CO1)			#  Fetch 16 C
	MADD	t11,c11,t11,ALPHA
	ST	t11,0(CO1)




.L999:							#  End
	ld	$16,   0($sp)
	ld	$17,   8($sp)
	ld	$18,  16($sp)
	ld	$19,  24($sp)
	ld	$20,  32($sp)
	ld	$21,  40($sp)
	ld	$22,  48($sp)
	LD	$f24, 56($sp)
	LD	$f25, 64($sp)
	LD	$f26, 72($sp)
	LD	$f27, 80($sp)
	LD	$f28, 88($sp)
	ld	$23,  96($sp)
	ld	$24, 104($sp)
	ld	$25, 112($sp)
	LD	$f20,120($sp)
	LD	$f21,128($sp)
	LD	$f22,136($sp)
	LD	$f23,144($sp)

	j	$31
	daddiu	$sp, $sp, 160

	EPILOGUE