mtd4.asm

Tiger sharec系列DSP的短点FFT代码(4,6,8,16点FFT),全汇编实现,将DSP的性能发挥到极致,在超高实时要求应用下有意义.

//15个时钟周期

#include "TS101_def.h"
.section M1DATA;
.align 4;
.var fft4Result[12];
.var win[8]="win8.dat";
.section M2DATA;
.align 4;
.var fft4input[8]="fft4input.dat";//order is re,im,re,im...

.section M0CODE;
.global _main;
_main: 
		k0=fft4input;;
		j0=win;;
		j1=fft4Result;;
		r0=l[k0+=2];r4=l[j0+=2];;//yr1=real[x(0),xr2=imag[x(0)]
		r1=l[k0+=2];r5=l[j0+=2];;//yr2=real[x(1)],xr2=imag[x(1)]
		r2=l[k0+=2];r6=l[j0+=2];fr0=r0*r4;;//yr3=real[x(2)],xr3=imag[x(2)]
		r3=l[k0+=0];r7=l[j0+=0];fr1=r1*r5;;//yr0=real[x(3)],xr0=imag[x(3)]
        fr2=r2*r6;r20=-1.0000;;
        fr3=r3*r7;;
        fr8=r0+r2,fr9=r0-r2;;
        fr10=r1+r3,fr11=r1-r3;;
        yr21=xr11;;
        fr14=r8+r10,fr16=r8-r10;yfr22=r20*r11;;
        xr21=yr22;;
        l[j1+=4]=r14;;
        fr15=r9+r21,fr17=r9-r21;l[j1+=-2]=r16;;
        l[j1+=4]=r15;;
        l[j1+=0]=r17;;
	    nop;;
		nop;;
		
		
		
		
		/*fr1=r1*r5;;
		
		fr8=r0+r2,fr9=r0-r2;;
		fr3=r3*r7;;
		fr10=r1+r3,fr11=r1-r3;;
		xr21=yr11;;
		fr4=r8+r10,fr6=r8-r10;yr20=xr11;;//R4放第一个X(0)的值,R6放第三个X(2)的值
		yfr21=r20*r22;l[j1+=4]=r4;;
		fr5=r9-r21;l[j1+=-2]=r6;;
        fr7=r9+r21;l[j1+=4]=r5;;
		l[j1+=0]=r7;;
		*/
_main.end: