#define USRLIBS TDFLIB
#include <SEP.defs>
#define FIGLIST mshole
#define FIGLIST msiter mshole

#define TA2VPLOT Ta2vplot titlesz=15 > /dev/null title=" " wantaxis=0
#define BYTE Byte Xpclip=100

default& : tube

#define SITER 1
#define NITER 30
FIGDIR/msiter.v : make.x hole.x Msmiss.x j1.H j1-4.H
	Zero  aa.H
	        make.x  n2=20 > make.H
	<make.H hole.x        > hole.H
	<hole.H Msmiss.x msc=j1-4.H niter=SITER > mbeg.H; <aa.H Disfil
	<hole.H Msmiss.x msc=j1-4.H niter=NITER > mend.H; grep norm mend.H
	<hole.H Msmiss.x msc=j1.H   niter=SITER > sbeg.H
	<hole.H Msmiss.x msc=j1.H   niter=NITER > send.H; grep norm send.H
	<sbeg.H BYTE|TA2VPLOT title=single_scale:_iter=SITER out=sbeg.v
	<send.H BYTE|TA2VPLOT title=single_scale:_iter=NITER out=send.v
	<mbeg.H BYTE|TA2VPLOT title=multiscale:___iter=SITER out=mbeg.v
	<mend.H BYTE|TA2VPLOT title=multiscale:___iter=NITER out=mend.v
	vp_SideBySideAniso sbeg.v send.v > single.v
	vp_SideBySideAniso mbeg.v mend.v > multip.v
	vp_OverUnderAniso  single.v multip.v > junk.v
	<junk.v  vppen vpstyle=n txscale=2. >FIGDIR/msiter.v

#define SKIP 4
#define SKIP 3
FIGDIR/mshole.v! : make.x hole.x strp.x Msmiss.x j1-SKIP.H
	Zero aa.H
	        make.x  n2=20     > make.H
	<make.H hole.x            > hole.H
	<hole.H strp.x  skip=SKIP > strp.H
	<hole.H Msmiss.x msc=j1-SKIP.H > gap.H; < aa.H Disfil; grep norm gap.H
	        Zero aa.H
	<strp.H Msmiss.x msc=j1-SKIP.H > gap.H; < aa.H Disfil; grep norm gap.H
	<make.H BYTE | TA2VPLOT title=original out=orig.v
	<strp.H BYTE | TA2VPLOT title=gapped   out=strp.v
	<gap.H  BYTE | TA2VPLOT title=restored out=gap.v
	vp_SideBySideAniso orig.v strp.v gap.v > junk.v
	<junk.v  vppen vpstyle=n txscale=2. >FIGDIR/mshole.v

test2&: testtcai2.x testicaf2.x
	testtcai2.x hello
	testicaf2.x hello

test1&: testicaf1.x testtcai1.x
	testicaf1.x hello
	testtcai1.x hello

j%.H : j%.HH
	< j%.HH Dd esize=4 > j%.H 

clean&: jclean texclean

#include <SEP.idoc.rules>
#include <SEP.prog.rules>
#include <SEP.obj.rules>
