UF77LIBS = ${TDFLIB}
UTUBEFLAGS = see_progress=y
include ${SEPINC}/SEP.top
RESULTSER = # im1-2+1 im1-1a fm1-1a early1-1a usual1-1a schwab1-1a whitvar whititer
RESULTSER = im1-2+1 im1-1a early1-1a usual1-1a schwab1-1a whitvar whititer

default : tube

GRAPH = Graph clip=1 max2=1 min2=-1 crowd=.99 > junk.H
DOTY = Dots dots=0 connect=0 gaineach=0 constsep=1 strings=0 clip=1.2 > /dev/null
O1D1 = `<alias.P Get o1 d1`
LBL = label1="Iterations"

${RESDIR}/im%.v ${RESDIR}/im%.v3 im%.H ir%.H: Alias.x Invint.x %.H 
	Zero im$*.H
	Zero ir$*.H
	Alias.x tmax=80. pairs=1   n1=200 > alias.P
	Alias.x tmax=80. pairs=0   n1=200 > alias.H
	<alias.P Invint.x n1=200 ${O1D1} rr=ir$*.H mm=im$*.H > levint.H filtin=$*.H
	<im$*.H Transp plane=23     > transp.H
	<transp.H ${GRAPH} out=${RESDIR}/im$*.v3 fastplot=20
	<levint.H ${GRAPH} out=levint.v
	< alias.H ${GRAPH} out=alias.v
	< alias.H ${DOTY}  strings=1 out=alias.v
	<levint.H ${DOTY}  connect=1 out=levint.v
	vp_Overlay alias.v levint.v > ${RESDIR}/im$*.v

# ${RESDIR}/fm%.v ${RESDIR}/fm%.v3 fm%.H fr%.H: Alias.x Invintp.x %.H
fm%.H                               fr%.H: Alias.x Invintp.x %.H
	Zero fm$*.H
	Zero fr$*.H
	Alias.x tmax=80. pairs=1   n1=200 > alias.P
	Alias.x tmax=80. pairs=0   n1=200 > alias.H
	<alias.P Invintp.x n1=200 ${O1D1} rr=fr$*.H mm=fm$*.H rho=.7 > levint.H filtin=$*.H
	<fm$*.H Transp plane=23     > transp.H
#	<transp.H ${GRAPH} out=${RESDIR}/fm%.v3 fastplot=20
	<levint.H ${GRAPH} out=levint.v
	< alias.H ${GRAPH} out=alias.v
	< alias.H ${DOTY}  strings=1 out=alias.v
	<levint.H ${DOTY}  connect=1 out=levint.v
#	vp_Overlay alias.v levint.v > ${RESDIR}/fm%.v

NITER = 250
ITERATION = 5
${RESDIR}/early%.v : im%.H fm%.H    
	< im$*.H Window n2=1 f2=${ITERATION}   > iwind.H
	< fm$*.H Window n2=1 f2=${ITERATION}   > fwind.H
	< fm$*.H Window n2=1 f2=${NITER} > icorr.H
	Merge <iwind.H fwind.H         icorr.H space=no axis=2 > junk.H
	< junk.H ${DOTY} connect=1 overlap=2.0 title=" " labels=no_prec.${ITERATION}:precond.${ITERATION}:niter=${NITER} labelsz=10 out=${RESDIR}/early$*.v

${RESDIR}/usual%.v :  Resnorm.x  ir%.H fr%.H 
	<ir$*.H Resnorm.x > junk1.H
	<fr$*.H Resnorm.x > junk2.H
#	Merge  junk1.H junk2.H         space=n axis=1 > junk.H
	Merge <junk1.H junk2.H         space=n axis=1 > junk.H
	< junk.H Transp plane=12  > transp.H
	<transp.H Window f1=0 n1=20 > wind.H
	< wind.H Graph ${LBL}  symbol='np' symbolsz=10 symbolsz=7 title="Data Residual Vector Length" min2=0 min1=-1 out=${RESDIR}/usual$*.v > junk.V 

${RESDIR}/schwab%.v :  Modnorm.x  im%.H fm%.H 
	<im$*.H Modnorm.x > junk1.H
	<fm$*.H Modnorm.x > junk2.H
#	Merge  junk1.H junk2.H         space=n axis=1 > junk.H
	Merge <junk1.H junk2.H         space=n axis=1 > junk.H
	< junk.H Transp plane=12  > transp.H
	<transp.H Window f1=0 n1=100 > wind.H
	< wind.H Graph ${LBL}  symbol='np' symbolsz=10 symbolsz=7 title="Model Residual Vector Length" min2=0 min1=-1 out=${RESDIR}/schwab$*.v > junk.V 

${RESDIR}/whitvar.v ${RESDIR}/whitvar.v3 whitvar.H: Alias.x Invintp2.x 1-1a.H precond.p
	Zero whitvar.H
	Alias.x tmax=80. pairs=1   n1=200 > alias.P
	Alias.x tmax=80. pairs=0   n1=200 > alias.H
	<alias.P Invintp2.x n1=200 ${O1D1} mm=whitvar.H> invintp.H  par=precond.p
	<whitvar.H Transp plane=23 > transp.H;
	<transp.H  ${GRAPH} fastplot=20	    	out=${RESDIR}/whitvar.v3
	<invintp.H ${GRAPH} out=invintp.v
	< alias.H  ${GRAPH} out=alias.v
	< alias.H  ${DOTY}  strings=1 out=alias.v
	<invintp.H ${DOTY}  connect=1 out=invintp.v
	vp_Overlay alias.v invintp.v >              ${RESDIR}/whitvar.v

LABITER = labels=iter=200:iter=5:iter=2 labelsz=10 
${RESDIR}/whititer.v : whitvar.H
	< whitvar.H Window n2=1 f2=2   > it.2.H
	< whitvar.H Window n2=1 f2=5   > it.5.H
	< whitvar.H Window n2=1 f2=200 > it.200.H
	Merge < it.200.H it.5.H it.2.H space=no axis=2 > junk.H
	< junk.H ${DOTY} connect=1 overlap=2.0 title=" " ${LABITER} out=${RESDIR}/whititer.v

1%.H : 1%.HH
	<1$*.HH Dd esize=4 > 1$*.H

clean: jclean
	@-touch junk.P
	@-/bin/rm *.P

#       hardcopy to make the mean lines darker.
HARDCOPY =

include ${SEPINC}/SEP.bottom
