#define USRLIBS TDFLIB
#define TUBEARGS see_progress=y
#include <SEP.defs>
$#define FIGLIST im1-1b    ir1-1b    fm1-1b    fr1-1b    rr1-1b  mattschwab1-1b
$#define FIGLIST im1-1c    ir1-1c    fm1-1c    fr1-1c    rr1-1c  mattschwab1-1c
$#define FIGLIST im1-2+1   ir1-2+1   fm1-2+1   fr1-2+1     mattschwab1-2+1
$#define FIGLIST im1-.5-.5 ir1-.5-.5 fm1-.5-.5 fr1-.5-.5   mattschwab1-.5-.5
$#define FIGLIST im1-1a ir1-1a fm1-1a fr1-1a rr1-1a mattschwab1-1a
$#define FIGLIST im1-2+1 im1-1a mattschwab1-1a rr1-1a
$#define FIGLIST im1-2+1 im1-1a fm1-1a tm1-1a mattschwab1-1a rr1-1a
#define FIGLIST im1-2+1 im1-1a fm1-1a mattschwab1-1a rr1-1a
$#define FIGLIST im10001

default& : tube

frob&: Frobtest.x
	Frobtest.x n=13

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

FIGDIR/im%.v! FIGDIR/im%.v3 FIGDIR/ir%.v! FIGDIR/ir%.v3 im%.H ir%.H: Alias.x Invint.x %.H 
$FIGDIR/%.v!            : 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=FIGDIR/im%.v3 fastplot=20
	<ir%.H Transp plane=23     > transp.H
	<transp.H GRAPH out=FIGDIR/ir%.v3 fastplot=20
	<transp.H Window f3=200 > junk.H
	<  junk.H DOTY connect=1 out=FIGDIR/ir%.v
	<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 > FIGDIR/im%.v

FIGDIR/fm%.v! FIGDIR/fm%.v3 FIGDIR/fr%.v! FIGDIR/fr%.v3 fm%.H fr%.H: Alias.x Invfrob.x %.H 
$FIGDIR/%.v!            : Alias.x Invfrob.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 Invfrob.x n1=200 O1D1 rr=fr%.H mm=fm%.H > levint.H filtin=%.H
	<fm%.H Transp plane=23     > transp.H
	<transp.H GRAPH out=FIGDIR/fm%.v3 fastplot=20
	<fr%.H Transp plane=23     > transp.H
	<transp.H GRAPH out=FIGDIR/fr%.v3 fastplot=20
	<transp.H Window f3=200 > junk.H
	<  junk.H DOTY connect=1 out=FIGDIR/fr%.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 > FIGDIR/fm%.v

FIGDIR/tm%.v! FIGDIR/tm%.v3 FIGDIR/tr%.v! FIGDIR/tr%.v3 tm%.H tr%.H: Alias.x Invtrob.x %.H 
$FIGDIR/%.v!            : Alias.x Invtrob.x %.H 
	Zero tm%.H
	Zero tr%.H
	Alias.x tmax=80. pairs=1   n1=200 > alias.P
	Alias.x tmax=80. pairs=0   n1=200 > alias.H
	<alias.P Invtrob.x n1=200 O1D1 rr=tr%.H mm=tm%.H > levint.H filtin=%.H
	<tm%.H Transp plane=23     > transp.H
	<transp.H GRAPH out=FIGDIR/tm%.v3 fastplot=20
	<tr%.H Transp plane=23     > transp.H
	<transp.H GRAPH out=FIGDIR/tr%.v3 fastplot=20
	<transp.H Window f3=200 > junk.H
	<  junk.H DOTY connect=1 out=FIGDIR/tr%.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 > FIGDIR/tm%.v


FIGDIR/mattschwab%.v : im%.H fm%.H
$FIGDIR/mattschwab%.v : im%.H fm%.H tm%.H
	< im%.H Window n2=1 f2=5   > iwind.H
	< fm%.H Window n2=1 f2=5   > fwind.H
$	< tm%.H Window n2=1 f2=5   > twind.H
	< im%.H Window n2=1 f2=200 > 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_frob_5:frob_5:trob_5:niter_200 out=FIGDIR/mattschwab%.v
	< junk.H DOTY connect=1 overlap=2.0 title=" " labels=no_frob_5:frob_5:niter_200 out=FIGDIR/mattschwab%.v

FIGDIR/rr%.v :  Vecnorm.x  ir%.H fr%.H 
$FIGDIR/rr%.v :  Vecnorm.x  ir%.H fr%.H tr%.H
	<ir%.H Vecnorm.x > junk1.H
	<fr%.H Vecnorm.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 n1=20 > wind.H
$	< wind.H Graph LBL  symbol='nft' symbolsz=10 symbolsz=7 title="Total Residual" out=FIGDIR/rr%.v > junk.V 
	< wind.H Graph LBL  symbol='nf' symbolsz=10 symbolsz=7 title="Total Residual" out=FIGDIR/rr%.v > junk.V 

Invtrob.x : Invtrob.o invtrob1.o trob1.o
        FLINK (Invtrob.o invtrob1.o trob1.o, Invtrob.x)

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

clean&: jclean
	@-/bin/rm \*.P

$       hardcopy to make the mean lines darker.
#define HARDCOPY
#include <SEP.idoc.rules>
#include <SEP.obj.rules>
#include <SEP.prog.rules>

