c**d54
cmx
cl
;
;no grad sort
;collected on glacier
;
;
def input fidc
ty ----------------------------------------------
	ty building matrix . . . 
	bld ir24_hbcbcgcdhd.mat 2 1024 128 0 y
	mat ir24_hbcbcgcdhd.mat w
;
;
;
;
;
esc out
if &out ne 0 quit
;
def phase0 0 
def phase1 0
def mult 1.0
for row 1 48
re &input
bc 
gv 1 rvalue
exc 
gv 1 ivalue
exc
eva real (&mult*&rvalue)
eva imag (&mult*&ivalue)
srv 1 1 &real &imag
;cnv 0 75
ss 512 90 
bc
zf 1024
ft
ph
;abq 
;pol 4
red 
sto 0 &row
ty row &row complete $
next row
;end
loa 0 1
mul 0.5 
sto 0 1
ty D1 transform complete
cl
def phase0 0.0
for col 1 1024
loa &col 0
def datype 1
def datsiz 48
;lpl 250 10 10 251 375
qsb 32 90
zf 128 
ft
ph
;abq 
;pol 4
rev
red
sto &col 0
ty column &col complete $
next col
ty
end