Motorola DSP96002 User Manual page 631

page
;
; Implements the Walsh-Hadamard Transform
;
iord equ
n equ
org
data
dc
dc
dc
dc
dc
dc
dc
dc
dc
dc
dc
dc
dc
dc
dc
dc
org
start
move
move
do
move
move
faddsub.s d0,d1
move
faddsub.s d2,d3
move
move
B-112
132,60,1,1
4 ;order of transform=log2(npoints)
1<<iord ;length of transform
x:$1000
0.0000000E+00
2.000000
3.000000
8.000000
9.000000
12.00000
15.00000
19.00000
20.00000
22.00000
23.00000
24.00000
25.00000
26.00000
27.00000
28.00000
p:$100
#data,r0 ;point to upper leg
#data+n/2,r4 ;point to lower leg
#n/4,_firststage ;do first stage. split into X and Y
x:(r0)+,d0.s
x:(r4)+,d1.s
x:(r0)-,d2.s
x:(r4)+,d3.s
d1.s,x:(r0)
d0.s,y:(r0)+
d3.s,x:(r0)
DSP96002 USER'S MANUAL
;get upper leg of bfly 1
;get lower leg of bfly 1
;get upper leg of bfly 2
;get lower leg of bfly 2
;save sum 1
;save dif 1
;save sum 2
MOTOROLA