; file firlib\src\fir_sp_asm_manual.s
; brief FIR filtering function in C6600 assembler with typical VLIW optimization
; (single precision IEEE754)
; author 
; date 

	.global fir_sp_asm_vliw

fir_sp_asm_vliw:
				;int i, j;
				MV				A8,A0		;i<nyk
				SUB				B6,1,B6
												
L1:				; input array loop			
				;for (i=0; i<nyk; i++) {
				;	tmp = 0;
				ZERO			A5			
	||			MV				B6,A1		;j<na	
	||			MV				A4,A19
	||			MV				B4,B19					
											
L2:				; FIR filter algorithm - dot product	
				;for (j=0; j<na; j++){
				;	tmp += a[j]*xk[i+j];
				LDW				*A19++,A3			
	 ||			LDW				*B19++,B9
				NOP				4				
				MPYSP			A3,B9,A17
		[A1] 	BDEC			L2,A1
				NOP				2
				FADDSP			A17,A5,A5					
				NOP				2			
				; end of FIR filter algorithm - dot product	
				
		[A0]	BDEC			L1,A0				
				STW				A5,*A6++	;yk(i) = tmp;
				ADD				A4,4,A4  	;i+j (float -> 4b)										
				NOP				3
				; end of input array loop	
						
				B				B3
				NOP				5
				; end of fir_sp_asm_vliw procedure