MJ = cd ~/research.ln/Source/Oneway/; make $(bin)/joint_encode.x ; cd 

norne-dmult-flat-negpos-06-%:
	$(MJ)
	##+ve
	#Window3d squeeze=n min2=0  < norne-dmult-negpos-$*-b4mult-06.H 	       			> junk$@.H ; < junk$@.H Transp reshape=2,3,4 plane=13 | Transp   > $@-pos.H 
	#Window3d squeeze=n min2=0  < norne-dmult-pos-pri-$*-06.H 				> junk$@.H ; < junk$@.H Transp reshape=2,3,4 plane=13 | Transp   > $@-pos-rm.H 
	#$B/joint_encode.x report=$@.log modl=$@-pos.H  \
		modl1=$@-pos.H rm1=$@-pos-rm.H  smodl1=$@-pos-s.H ts1=$@-pos-ts.H sd1=mod-flat-tsd1.H cc1=mod-flat-cc1.H rr1=mod-flat-rr1.H \
        	nmix=1 par=flat-ang.p nmix=5 nodes=1 shift_nodes=1 lflat=y &
	##-ve
	#Window3d squeeze=n n2=61   < norne-dmult-negpos-$*-b4mult-06.H | Reverse which=2       	> junk$@.H ; < junk$@.H Transp reshape=2,3,4 plane=13 | Transp   > $@-neg.H
	#Window3d squeeze=n n2=61   < norne-dmult-neg-pri-$*-06.H       | Reverse which=2       	> junk$@.H ; < junk$@.H Transp reshape=2,3,4 plane=13 | Transp   > $@-neg-rm.H 
	#$B/joint_encode.x report=$@.log modl=$@-neg.H  \
		modl1=$@-neg.H rm1=$@-neg-rm.H  smodl1=$@-neg-s.H ts1=$@-neg-ts.H sd1=mod-flat-tsd1.H cc1=mod-flat-cc1.H rr1=mod-flat-rr1.H \
        	nmix=1 par=flat-ang.p nmix=5 nodes=1 shift_nodes=1 lflat=y 
	#touch $@
	Window3d n4=1 n6=61 	 < $@-neg-ts.H | Transp reshape=1,3,4 | Reverse which=2 	  > $@-tmp1-ts.H 
	Window3d n4=1 n6=60 f6=1 < $@-pos-ts.H | Transp reshape=1,3,4 			  > $@-tmp2-ts.H 
	Cat3d $@-tmp1-ts.H $@-tmp2-ts.H axis=2 > $@-ts.H
	Transp plane=24 < $@-neg-s.H | Window3d n2=61      | Transp plane=13 | Reverse which=2 > $@-tmp1.H
	Transp plane=24 < $@-pos-s.H | Window3d n2=60 f2=1 | Transp plane=13 	   	       > $@-tmp2.H
	Cat3d $@-tmp1.H $@-tmp2.H axis=2 > $@-s.H
# Aplitude balance
rr = 1.10228805034138824083 #ratio of the rms wnergy btw base and moni
norne-dmult-amp-08-%:
	Math exp="file1*$(rr)" file1=norne-dmult-negpos-$*-b4mult-08.H           >  norne-dmult-negpos-$*-b4mult-amp-08.H 
norne-hole-sxgx-dmult-amp-08-%:
	Math exp="file1*$(rr)" file1=norne-hole-sxgx-dmult-negpos-$*-b4mult-08.H >  norne-hole-sxgx-dmult-negpos-$*-b4mult-amp-08.H 
norne-fill-hole-sxgx-dmult-amp-08-%:
	Math exp="file1*$(rr)" file1=norne-fill-hole-sxgx-dmult-negpos-$*-b4mult-08.H >  norne-fill-hole-sxgx-dmult-negpos-$*-b4mult-amp-08.H 

# ADCIG Warp after Flatten
norne-dmult-flat-warp-%:
	$(MJ)
	Window3d squeeze= n f2= 0 n2= 11 j1= 2 n1= 320 n4= 100 j2= 2 j3= 2 < norne-dmult-flat-negpos-06-$*-s.H     > junk$@.H ; < junk$@.H Transp plane=13  > $@-06.H 
	Window3d squeeze= n f2= 0 n2= 11 j1= 2 n1= 320 n4= 100 j2= 2 j3= 2 < norne-dmult-negpos-$*-b4mult-amp-08.H > junk$@.H ; < junk$@.H Transp plane=13  > $@-08.H 
	$B/joint_encode.x report=$@.log modl=$@-06.H  \
		modl1=$@-06.H rm1=$@-06.H  smodl1=$@-06-s.H .H ts1=$@-06-ts.H 2.H sd1=$@-06-tsd.H .H cc1=$@-06-cc.H H rr1=$@-06-rr.H  \
		modl2=$@-08.H rm2=$@-08.H  smodl2=$@-08-s.H .H ts2=$@-08-ts.H 2.H sd2=$@-08-tsd.H .H cc2=$@-08-cc.H H rr2=$@-08-rr.H  \
        	nmix=1 par=warp-ang.P nmix=7 nodes=1 shift_nodes=1 verb=8 &
	touch $@
norne-hole-sxgx-dmult-flat-warp-%:
	$(MJ)
	Window3d squeeze= n f2= 0 n2= 11 j1= 2 n1= 320 n4= 100 j2= 2 j3= 2 < norne-dmult-flat-negpos-06-$*-s.H    	      > junk$@.H ; < junk$@.H Transp plane=13  > $@-06.H 
	Window3d squeeze= n f2= 0 n2= 11 j1= 2 n1= 320 n4= 100 j2= 2 j3= 2 < norne-fill-hole-sxgx-dmult-negpos-$*-b4mult-08.H > junk$@.H ; < junk$@.H Transp plane=13  > $@-08-rm.H 
	Window3d squeeze= n f2= 0 n2= 11 j1= 2 n1= 320 n4= 100 j2= 2 j3= 2 < norne-hole-sxgx-dmult-negpos-$*-b4mult-amp-08.H  > junk$@.H ; < junk$@.H Transp plane=13  > $@-08.H 
	$B/joint_encode.x report=$@.log modl=$@-06.H  \
		modl1=$@-06.H rm1=$@-06.H    smodl1=$@-06-s.H .H ts1=$@-06-ts.H 2.H sd1=$@-06-tsd.H .H cc1=$@-06-cc.H H rr1=$@-06-rr.H  \
		modl2=$@-08.H rm2=$@-08-rm.H smodl2=$@-08-s.H .H ts2=$@-08-ts.H 2.H sd2=$@-08-tsd.H .H cc2=$@-08-cc.H H rr2=$@-08-rr.H  \
        	nmix=1 par=warp-ang.P nmix=7 nodes=1 shift_nodes=1 verb=8 &
norne-hole-hole-sxgx-dmult-flat-warp-%:
	$(MJ)
	Window3d squeeze= n f2= 0 n2= 11 j1= 2 n1= 320 n4= 100 j2= 2 j3= 2 < norne-dmult-flat-negpos-06-$*-s.H           > junk$@.H ; < junk$@.H Transp plane=13  > $@-06.H 
	Window3d squeeze= n f2= 0 n2= 11 j1= 2 n1= 320 n4= 100 j2= 2 j3= 2 < norne-dmult-negpos-$*-b4mult-06.H           > junk$@.H ; < junk$@.H Transp plane=13  > $@-06-rm.H 
	Window3d squeeze= n f2= 0 n2= 11 j1= 2 n1= 320 n4= 100 j2= 2 j3= 2 < norne-hole-sxgx-dmult-negpos-$*-b4mult-06.H > junk$@.H ; < junk$@.H Transp plane=13  > $@-066.H 
	$B/joint_encode.x report=$@.log modl=$@-06.H  \
		modl1=$@-06.H  rm1=$@-06.H    smodl1=$@-00-s.H .H ts1=$@-00-ts.H 2.H sd1=$@-00-tsd.H .H cc1=$@-00-cc.H H rr1=$@-00-rr.H  \
		modl2=$@-066.H rm2=$@-06-rm.H smodl2=$@-06-s.H .H ts2=$@-06-ts.H 2.H sd2=$@-06-tsd.H .H cc2=$@-06-cc.H H rr2=$@-06-rr.H  \
        	nmix=1 par=warp-ang.P nmix=7 nodes=1 shift_nodes=1 verb=8 &
	touch $@
norne-fill-hole-sxgx-dmult-flat-warp-%:
	$(MJ)
	Window3d squeeze= n f2= 0 n2= 11 j1= 2 n1= 320 n4= 100 j2= 2 j3= 2 < norne-dmult-flat-negpos-06-$*-s.H                     > junk$@.H ; < junk$@.H Transp plane=13  > $@-06.H 
	Window3d squeeze= n f2= 0 n2= 11 j1= 2 n1= 320 n4= 100 j2= 2 j3= 2 < norne-fill-hole-sxgx-dmult-negpos-$*-b4mult-amp-08.H  > junk$@.H ; < junk$@.H Transp plane=13  > $@-08.H 
	$B/joint_encode.x report=$@.log modl=$@-06.H  \
		modl1=$@-06.H rm1=$@-06.H smodl1=$@-06-s.H .H ts1=$@-06-ts.H 2.H sd1=$@-06-tsd.H .H cc1=$@-06-cc.H H rr1=$@-06-rr.H  \
		modl2=$@-08.H rm2=$@-08.H smodl2=$@-08-s.H .H ts2=$@-08-ts.H 2.H sd2=$@-08-tsd.H .H cc2=$@-08-cc.H H rr2=$@-08-rr.H  \
        	nmix=1 par=warp-ang.P nmix=7 nodes=1 shift_nodes=1 verb=8 &
	touch $@
warp-dmult-negpos-%:
	make -n norne-dmult-flat-warp-$*
	make -n norne-hole-dmult-flat-warp-$* &
	make -n norne-hole-sxgx-dmult-flat-warp-$*
	make -n norne-fill-sxgx-dmult-flat-warp-$*

# stacked images
mig-norne-dmult-flat-warp-%:
	Window3d min3=1500 n3=501 < norne-dmult-flat-warp-$*-06.H   | Stack3d | Transp plane=23 > $@-06.H; echo o2=0 d2=25 >> $@-06.H & 
	Window3d min3=1500 n3=501 < norne-dmult-flat-warp-$*-08-s.H | Stack3d | Transp plane=23 > $@-08.H; echo o2=0 d2=25 >> $@-08.H 
mig-norne-hole-sxgx-dmult-flat-warp-%:
	Window3d min3=1500 n3=501 < norne-hole-sxgx-dmult-flat-warp-$*-06.H   | Stack3d | Transp plane=23 > $@-06.H; echo o2=0 d2=25 >> $@-06.H  &
	Window3d min3=1500 n3=501 < norne-hole-sxgx-dmult-flat-warp-$*-08-s.H | Stack3d | Transp plane=23 > $@-08.H; echo o2=0 d2=25 >> $@-08.H 
mig-norne-fill-hole-sxgx-dmult-flat-warp-%:
	Window3d min3=1500 n3=501 < norne-fill-hole-sxgx-dmult-flat-warp-$*-06.H   | Stack3d | Transp plane=23 > $@-06.H; echo o2=0 d2=25 >> $@-06.H  %
	Window3d min3=1500 n3=501 < norne-fill-hole-sxgx-dmult-flat-warp-$*-08-s.H | Stack3d | Transp plane=23 > $@-08.H; echo o2=0 d2=25 >> $@-08.H 
mig-norne-hole-hole-sxgx-dmult-flat-warp-%:
	Window3d min3=1500 n3=501 < norne-hole-hole-sxgx-dmult-flat-warp-$*-06.H   | Stack3d | Transp plane=23 > $@-06.H; echo o2=0 d2=25 >> $@-06.H  
#warp stacked images before inversion
mig-shift-norne-dmult-flat-warp-%: mig-norne-dmult-flat-warp-%
	$(MJ)
	$B/joint_encode.x report=$@.log modl=mig-norne-dmult-flat-warp-$*-06.H   \
		modl1=mig-norne-dmult-flat-warp-$*-06.H rm1=mig-norne-dmult-flat-warp-$*-06.H  smodl1=$@-06-s.H .H ts1=$@-06-ts.H 2.H sd1=$@-06-tsd.H .H cc1=$@-06-cc.H H rr1=$@-06-rr.H  \
		modl2=mig-norne-dmult-flat-warp-$*-08.H rm2=mig-norne-dmult-flat-warp-$*-08.H  smodl2=$@-08-s.H .H ts2=$@-08-ts.H 2.H sd2=$@-08-tsd.H .H cc2=$@-08-cc.H H rr2=$@-08-rr.H  \
        	nmix=1 par=warp-ang.P nmix=7 nodes=1 shift_nodes=1 verb=8 &
	touch $@
mig-shift-norne-hole-sxgx-dmult-flat-warp-%: mig-hole-norne-dmult-flat-warp-%
	$(MJ)
	$B/joint_encode.x report=$@.log modl=mig-norne-hole-sxgx-dmult-flat-warp-$*-06.H   \
		modl1=mig-norne-hole-sxgx-dmult-flat-warp-$*-06.H rm1=mig-norne-hole-sxgx-dmult-flat-warp-$*-06.H       smodl1=$@-06-s.H .H ts1=$@-06-ts.H 2.H sd1=$@-06-tsd.H .H cc1=$@-06-cc.H H rr1=$@-06-rr.H  \
		modl2=mig-norne-hole-sxgx-dmult-flat-warp-$*-08.H rm2=mig-norne-fill-hole-sxgx-dmult-flat-warp-$*-08.H  smodl2=$@-08-s.H .H ts2=$@-08-ts.H 2.H sd2=$@-08-tsd.H .H cc2=$@-08-cc.H H rr2=$@-08-rr.H  \
        	nmix=1 par=warp-ang.P nmix=7 nodes=1 shift_nodes=1 verb=8 &
	touch $@
mig-shift-norne-fill-hole-sxgx-dmult-flat-warp-%: mig-norne-fill-hole-dmult-flat-warp-%
	$(MJ)
	$B/joint_encode.x report=$@.log modl=mig-norne-fill-hole-sxgx-dmult-flat-warp-$*-06.H   \
		modl1=mig-norne-fill-hole-sxgx-dmult-flat-warp-$*-06.H rm1=mig-norne-fill-hole-sxgx-dmult-flat-warp-$*-06.H  smodl1=$@-06-s.H .H ts1=$@-06-ts.H 2.H sd1=$@-06-tsd.H .H cc1=$@-06-cc.H H rr1=$@-06-rr.H  \
		modl2=mig-norne-fill-hole-sxgx-dmult-flat-warp-$*-08.H rm2=mig-norne-fill-hole-sxgx-dmult-flat-warp-$*-08.H  smodl2=$@-08-s.H .H ts2=$@-08-ts.H 2.H sd2=$@-08-tsd.H .H cc2=$@-08-cc.H H rr2=$@-08-rr.H  \
        	nmix=1 par=warp-ang.P nmix=7 nodes=1 shift_nodes=1 verb=8 &
	touch $@
rwarp-dmult-%:
	make -n mig-shift-nornes-sxgx-dmult-flat-warp-$*
	make -n mig-shift-norne-hole-sxgx-dmult-flat-warp-$*
	make -n mig-shift-norne-fill-hole-sxgx-dmult-flat-warp-$*

warp-all:
	make norne-dmult-flat-negpos-06-2759 	#rmo
	make warp-dmult-negpos-2759		#warp
	make rwarp-dmult-2759			#warp
