c**d54    process first 1H-1H planes (real and imag in 15n)
cmx
cl
;
;
;
;
;
def input fidc
ty -----------------------------------------
   ty Building matrix...
   bld ec283_n15noe_wnoe  2 512 512 0 y
   mat ec283_n15noe_wnoe w
   ty Matrix opened and write-enabled...
   ty
;
; 				D1 processing
;
;
  esc out
  if &out ne 0 quit
	def phase0 -16
	def phase1 0
def swidth 6000
 for tier 1 32
     eva itier (2*&tier)
     eva rtier (&itier-1)
     ty  Processing tiers &rtier and &itier $
     esc out
;     if &out ne 0 scram
       ; Process phase=1, phase2=1
       re &input
       cnv 1 10
       qss 512 90 1
       ;gm -7 0.1
       zf 1024
       ft
       ph
       def datsiz 512
       ;pol 4
       ;abl 4 40
       red
sto 0 &rtier 
       ;
       ;process phase=2, phase2=1
       re &input
       cnv 1 10
       qss 512 90 1
       ;gm -7 0.1
       zf 1024
       ft
       ph
       def datsiz 512
       ;pol 4
       ;abl 4 40
       red
sto 0 &itier
   nex
   tim r t3
;end
bun 0
   bun 2
   ty There are &vector vectors to process.
   tim z
def phase0 0
def phase1 0
def swidth 1500
   for ivec 1 &vector
     lwb
     def datype 1 
;
; parameters for mirror image LP
   def datsiz 48 ;=27
   eva datso (&datsiz) ;=27
   eva shift (&datsiz-1) ;=26
   eva fpnt (&datsiz+1) ;=28
   eva lpnt (&shift+&datsiz) ;=53
   eva lpf (&lpnt+1) ;=54
   eva lpl (&lpnt+&datso) ;=80
   eva hafpnt (nint(&lpnt/2)) ;=27
;
; conduct mirror image LP
     def datsiz &lpnt
     srv &fpnt &lpnt 0 0
     shr &shift
     stb 1
     shl &shift
     def datsiz &datso
     cnj
     rev
     def datsiz &lpnt
     srv &datso &lpnt 0 0
     adb 1
     ldb 1
     ;
     ;The following LP setup looks good but defies std setup...
     lpl &lpnt 6 4 &lpf &lpl
     shl &shift
     eva datsiz (&datso+&datso)
     sb &datsiz 90
     zf &d2size
gv 1 reel imagin
eva reel (&reel*0.5)
eva imagin (&imagin*0.5)
pv 1 &reel &imagin
     ft
     ph
     ;rev
     red
     bas zero
     ;abq
     ;pol 4
     swb
     ty  D2-vector=&ivec $
   nex
   tim r t2
   bun 0
quit:
end