BIN=../Bin
DATA=../../Data

# Match simulated primaries to the data with the standard goal
# Mf=d to get an estimative of multiplaes. The kinematics of 
#the estmated multiples is perfect

#!!!!!!!!!!!!!!!!!!!!!!!!!!!
#
# Returning to the problem
#      (03/23/07)
#
#!!!!!!!!!!!!!!!!!!!!!!!!!!!


# do adaptive subtraction of the primaries to get the multiples
#
primest01.H:
	Window < prim_data2_h.H f3=110 j3=400 n3=2> target2cig.H
	Window < prim2_prc.H f3=110 j3=400 n3=2> junk2cig.H
	time ${BIN}/Ultimatch-new-all.x < target2cig.H\
		multiples_in=junk2cig.H mu=1. \
		mode=2 nsmode=2 wmode=1 nfilt=11,11 npatch=50,1 \
		op=1 niter=50 2D=1 eps=30.0 > $@

multest03.H:
	Window < prim_data2_h.H f3=30 n3=1> targetcig.H
	Window < multest01.H f3=30 n3=1> junk2cig.H
	time ${BIN}/Ultimatch-new-all.x < target2cig.H\
                multiples_in=junk2cig.H mu=1. \
                mode=2 nsmode=2 wmode=1 nfilt=11,11 npatch=50,1 \
                op=1 niter=50 2D=1 eps=10.0 > $@


mulest02_mu1_nf2141_p501_ni60_ep5.H:
	Window < prim_data2_h.H f3=430 n3=1> targetcig.H
	Window < primest01_tiny.H f3=430 n3=1> junkcig.H
	Window < multest01_tiny.H f3=430 n3=1> junk2cig.H
	time ${BIN}/Ultimatch-new-all.x < targetcig.H\
		primaries_in=junkcig.H multiples_in=junk2cig.H mu=1. \
                primaries_out=primest02_mu1_nf2111_np5011_ni60_ep30.H\
		mode=2 nsmode=2 wmode=1 nfilt=21,41 npatch=50,1 \
		op=2 niter=60 2D=1 eps=5.0 > $@

# Not implemented
mulest02_mode3.H:
	time ${BIN}/Ultimatch-new-all.x < prim_data1_h_samp.H \
                primaries_in=primest01_e80_samp.H multiples_in=mulest01_e80_samp.H mu=.6 \
                primaries_out=primest02_mode3.H\
                mode=3 nsmode=2 wmode=1 nfilt=31,31 npatch=10,1 \
                op=2 niter=60 2D=1 eps1=10.0 eps2=5. esp3=5. > $@
	
mulest02_mode2.H:
	time ${BIN}/Ultimatch-new-all.x < prim_data1_h_samp.H \
                primaries_in=primest01_e80_samp.H multiples_in=mulest01_e80_samp.H mu=.6 \
                primaries_out=primest02_mode3.H\
                mode=2 nsmode=2 wmode=1 nfilt=11,31 npatch=10,1 \
                op=2 niter=60 2D=1 eps=60.0 > $@

mul02_tst.H:
	time ${BIN}/Ultimatch-new-all.x < prim_data1_h_samp.H\
                multiples_in=mul_samp.H \
		mu=1. mode=2 nsmode=2 wmode=1 nfilt=31,31 npatch=5,1 \
		op=1 niter=25 2D=1 eps=30.0 > $@

###############################################################
# TEST THE PATCH SIZE
#
# do adaptive subtraction of the primaries to get the multiples
# in common angle panels for the whole line
mult_patch_tst.H:
	time ${BIN}/Ultimatch-new-all.x < target_CAP.H\
                primaries_in=junk_CAP.H multiples_in=diff01_prc_CAP.H mu=0.8 \
                primaries_out=prim_patch_tst.H\
                mode=2 nsmode=2 wmode=1 nfilt=21,11 npatch=35,35 \
                op=2 niter=40 2D=1 eps=15.0 > $@

################################################################
#
# Now after getting the primairies usinf inverse theory 
#
#
mult01.H: 
	time ${BIN}/Ultimatch-new-all.x < ${DATA}/prim_data1_h.H\
		primaries_in=${DATA}/prim_prc.H multiples_in=${DATA}/diff01_prc.H mu=0.8 \
		primaries_out=${DATA}/prim01.H\
		mode=2 nsmode=2 wmode=1 nfilt=21,21 npatch=30,5 \

####### Parallel test
prim_par01.H:
	time Parallel command=${BIN}/Ultimatch-new-all.x \
		files="prim_data1_h_samp.H,prim_prc_samp.H,$@" \
		tags="stdin,multiples_in,stdout" \
		mach_mfile=machinefile_200 axis=3,3,3 \
		usage="INPUT,INPUT,OUTPUT" \
                file_type="DISTRIBUTE,DISTRIBUTE,DISTRIBUTE" \
	        pverb=2 nblock=3 \
                mu=1. mode=2 nsmode=2 wmode=1 nfilt=11,11 npatch=50,50 \
                op=1 niter=30 2D=1 eps=30.0 > prim_par01.log

mul_par02_tst.H:
	time Parallel command=${BIN}/Ultimatch-new-all.x \
                files="prim_data1_h_samp.H,mul_samp.H,prim_prc_samp.H,$@" \
                tags="stdin,multiples_in,primaries_in,stdout" \
                mach_mfile=machinefile_200 axis=3,3,3,3 \
                usage="INPUT,INPUT,INPUT,OUTPUT" \
                file_type="DISTRIBUTE,DISTRIBUTE,DISTRIBUTE,DISTRIBUTE" \
                pverb=2 nblock=3 \
                mu=1. mode=2 nsmode=2 wmode=1 nfilt=11,11 npatch=50,50 \
                op=2 niter=30 2D=1 eps=30.0 > mul_par02.log

mult01_par_sing.H: 
	time Parallel command=${BIN}/Ultimatch-new-all.x \
		files="prim_data1_h.H,mulest01.H,$@" \
		tags="stdin,multiples_in,stdout" \
		mach_mfile=machinefile_200 axis=3,3,3 \
		usage="INPUT,INPUT,OUTPUT" \
		file_type="DISTRIBUTE,DISTRIBUTE,DISTRIBUTE" \
		pverb=2 nblock=90 \
		mu=1. mode=2 nsmode=2 wmode=1 nfilt=21,21 npatch=25,5 \
		op=1 niter=60 2D=1 eps=60. > mult01_par_sing.log 

mult02_par_sing.H: 
	time Parallel command=${BIN}/Ultimatch-new-all.x \
		files="prim_data1_h.H,mulest01.H,$@" \
		tags="stdin,multiples_in,stdout" \
		mach_mfile=machinefile_200 axis=3,3,3 \
		usage="INPUT,INPUT,OUTPUT" \
		file_type="DISTRIBUTE,DISTRIBUTE,DISTRIBUTE" \
		pverb=2 nblock=90 \
		mu=1. mode=2 nsmode=2 wmode=1 nfilt=21,31 npatch=50,1 \
		op=1 niter=50 2D=1 eps=5.0 > mult02_par_sing.log

mult01_tst.H:
	time ${BIN}/Ultimatch-new-all.x <prim_data1_h_samp.H \
		primaries_in=prim_prc_samp.H multiples_in=mul_samp.H mu=0.8 \
		primaries_out=primout_samp.H \
		mode=2 nsmode=2 wmode=1 nfilt=11,11 npatch=50,50 \
		op=2 niter=50 2D=1 eps=15. > $@

mult01_par.H:
	time Parallel command=${BIN}/Ultimatch-new-all.x \
		files="${DATA}/prim_data1_h.H,${DATA}/prim_prc.H,${DATA}/mulest01_e80.H,$@" \
		tags="stdin,primaries_in,multiples_in,stdout" \
		mach_mfile=machinefile_200 axis=3,3,3,3 \
		usage="INPUT,INPUT,INPUT,OUTPUT" \
		file_type="DISTRIBUTE,DISTRIBUTE,DISTRIBUTE,DISTRIBUTE" \
		pverb=2 nblock=20\
		mu=0.8 mode=2 nsmode=2 wmode=1 nfilt=11,11 npatch=50,50 \
		op=2 niter=50 2D=1 eps=15. > mult01_par.log