Motorola DSP56000 Manual page 584

sub
move
move
add
move
nop
move
add
move
move
jmp
_loop move
_endl
move
move
do #points,_add_offset
move
add
move
_add_offset
endm
;convert normal order to berglang order
norm2berg macro points,bergtable,twiddle
;points is actual size of table to be converting
move
move
move
do
move
move
move
move
move
data_temp
endm
; Real-Valued FFT for MOTOROLA DSP56000/1/2
;
; based on Glenn Bergland's algorithm
;
; ______________________________
rifft macro points,binlogsz,idata,odata,twiddle,bergtable
move
move
B - 22
BENCHMARK PROGRAMS
y0,a
b,x1
r0,y0
y0,b
b,r0
a,x:(r0+n0)
x1,b
b,x0
x1,b
_star
y1,a
#>offset,a
#bergtable,r0
x:(r0),B
A,B
B,x:(r0)+
;end of sincos macro
#bergtable,r0
#twiddle,r2
r2,r6
#points,data_temp
x:(r0)+,r3
r3,r7
x:(r3),a
y:(r7),b
a,x:(r2)+ b,y:(r6)+
#idata,r0
#points/4,n0
BENCHMARK PROGRAMS
;k-bergtabl[j]
;save b, x1=i
;y0=j=i+i
;b=j+i
;r0=j+i
;store bergtabl[j+i]
;b=j+i+i
;save b
;recover b=i
;recover a
;offset is the location of output data or twiddle
;r0=pointer of bergland table
;r2=twiddle pointer for X
;r6=twiddle pointer for Y
;get index
;get value
;write back
;r0 = ptr to a
;bflys in ea group, half at ea pass
MOTOROLA