addsub64.asm

《DSP原理及应用》源代码

		.title "ADDSUB64"
		.mmregs
		.def  	start,_c_int00 
		.bss 	x1 , 2 ,1,1
		.bss 	x3 , 2 ,1,1
		.bss 	y1 , 2 ,1,1
		.bss 	y3 , 1
		.bss 	y2 , 1
		.bss 	z1 , 2 ,1,1
		.bss 	z3 , 1
		.bss 	z2 , 1
		.bss 	w1 , 2 ,1,1
		.bss 	w3 , 2 ,1,1
table 	.long	12345678H		;x1x0
		.long	02468ACEH		;x3x2
		.long	22222222H		;y1y0
		.word  1357H,2468H		;y3,y2
		.long	44444444H		;z1z0
		.word  1020H,0B30AH		;z3,z2
        .text 
_c_int00
		b start 
		nop
		nop
start:  LD    	#x1 ,DP
		STM		#x1,AR1
		RPT	 	#11
		MVPD	table,*AR1+
		DLD 	x1,A 		; A = X1 X0
		DADD 	y1,A 		; A= X1 X0 + Y1 Y0,产生进位C

		DLD  	x3,B 		; B = X3 X2
		ADDC	y2,B		; B = X3 X2 + 00 Y2 + C
		ADD 	y3,16,B 	; B = X3 X2 + Y3 Y2 + C
		DSUB	z1,A		; A= X1 X0 + Y1 Y0-Z1Z0，产生借位C'
		DST 	A,w1		; W1W0= X1 X0 + Y1 Y0-Z1Z0
		SUBB	z2,B		; B = X3 X2 + Y3 Y2 + C - 00 Z2-C'
		SUB		z3,16,B		; B = X3 X2 + Y3 Y2 + C -Z3 Z2-C'
		DST 	B,w3		; W3 W2= X3 X2 + Y3 Y2 + C - Z3 Z2-C'
END:	B END
		.end