#define BINDIR ../../Bin/MTYPE
#define FIGDIR ../Fig
#define FIGLIST chaos
#define ALLPLOTS FIGDIR/undamped.v FIGDIR/damped.v FIGDIR/topharm.v FIGDIR/subharm.v FIGDIR/nonlin.v FIGDIR/nonlin2.v FIGDIR/nonlin3.v 
#include <SEP.defs>
#define FATMULT 1
#define PLOTSIZE xsize=2.5 ysize=1.5
#define GRIDSIZE grid=-1 gridnum=2,1 gridsize=2.6,1.6
#define TOTALGRID  grid=-1 gridnum=1,6 gridsize=6,1.6

default&: tube

BINDIR/Logistic Logistic: Logistic.o
	FC FCOPTS Logistic.o SEPLIB SYSLIB -o BINDIR/Logistic

#define DUFFOBJS  driver.o dlsode.o duffing.o jac.o
#define DUFFOBJ2  driver.o dlsode.o duffing2.o jac.o

BINDIR/duffing duffing:        DUFFOBJS
        FC LDOPTS DUFFOBJS -o BINDIR/duffing USRLIBS  SEPLIB SYSLIB

BINDIR/duffing2 duffing2:       DUFFOBJ2
        FC LDOPTS DUFFOBJ2  -o BINDIR/duffing2 USRLIBS  SEPLIB SYSLIB

$ generate figures

undamped.ascii: BINDIR/duffing duffing
	BINDIR/duffing omega=1.0 epsilon=1.0 delta=0.0 gamma=0.0 y1=0.0 y2=1.0 pts=100 tmax=20.00
	cp phase undamped.ascii

FIGDIR/undamped.v: undamped.H 
	Graph <undamped.H title='Phase Space t=20' label1='u' label2='v' out=j1.v >/dev/null
	Real <undamped.H |Ftplot out=j2.v >/dev/null
	vppen align=lb <j1.v |      vppen  PLOTSIZE  >j11.v
	vppen big=n align=lb <j2.v |vppen  PLOTSIZE  >j22.v
	vppen GRIDSIZE j11.v j22.v |vppen align=lb >FIGDIR/undamped.v
	
damped.ascii: BINDIR/duffing duffing
	BINDIR/duffing omega=1.0 epsilon=1.0 delta=2.0 gamma=0.0 y1=0.0 y2=1.0 pts=100 tmax=20.00
	cp phase damped.ascii

FIGDIR/damped.v: damped.H
	Graph <damped.H title='Phase Space t=20' label1='u' label2='v' out=j1.v >/dev/null
	Real <damped.H |Ftplot out=j2.v >/dev/null
	vppen align=lb <j1.v |      vppen  PLOTSIZE  >j11.v
	vppen big=n align=lb <j2.v |vppen  PLOTSIZE  >j22.v
	vppen GRIDSIZE j11.v j22.v |vppen align=lb >FIGDIR/damped.v


topharm.ascii: BINDIR/duffing duffing
	BINDIR/duffing omega=4.67 epsilon=40.0 delta=0.0 gamma=10.0 y1=0.0 y2=1.0 pts=100 tmax=100.00
	cp phase topharm.ascii

FIGDIR/topharm.v: topharm.H 
	Graph <topharm.H title='Phase Space t=100' label1='u' label2='v' out=j1.v >/dev/null
	Real <topharm.H |Ftplot out=j2.v >/dev/null
	vppen align=lb <j1.v |      vppen  PLOTSIZE  >j11.v
	vppen big=n align=lb <j2.v |vppen  PLOTSIZE  >j22.v
	vppen GRIDSIZE j11.v j22.v |vppen align=lb >FIGDIR/topharm.v

subharm.ascii:  BINDIR/duffing duffing
	BINDIR/duffing omega=4.67 epsilon=40.0 delta=0.0 gamma=0.67 y1=0.0 y2=1.0 pts=100 tmax=100.00
	cp phase subharm.ascii

FIGDIR/subharm.v: subharm.H 
	Graph <subharm.H title='Phase Space t=100' label1='u' label2='v' out=j1.v >/dev/null
	Real <subharm.H |Ftplot out=j2.v >/dev/null
	vppen align=lb <j1.v |      vppen  PLOTSIZE  >j11.v
	vppen big=n align=lb <j2.v |vppen  PLOTSIZE  >j22.v
	vppen GRIDSIZE j11.v j22.v |vppen align=lb >FIGDIR/subharm.v
	
nonlin.ascii:  BINDIR/duffing duffing
	BINDIR/duffing omega=1 epsilon=10.0 delta=2.0 gamma=10 y1=0.0 y2=1.0 pts=100 tmax=100.00
	cp phase nonlin.ascii

FIGDIR/nonlin.v: nonlin.H 
	Graph <nonlin.H title='Phase Space' label1='u' label2='v' out=j1.v >/dev/null
	Real <nonlin.H |Ftplot out=j2.v >/dev/null
	vppen align=lb <j1.v |      vppen  PLOTSIZE  >j11.v
	vppen big=n align=lb <j2.v |vppen  PLOTSIZE  >j22.v
	vppen GRIDSIZE j11.v j22.v |vppen align=lb >FIGDIR/nonlin.v
	
nonlin2.ascii:  BINDIR/duffing duffing
	BINDIR/duffing omega=10 epsilon=10.0 delta=2.0 gamma=10 y1=0.0 y2=1.0 pts=100 tmax=100.00
	cp phase nonlin2.ascii

FIGDIR/nonlin2.v: nonlin2.H 
	Graph <nonlin2.H title='Phase Space' label1='u' label2='v' out=j1.v >/dev/null
	Real <nonlin2.H |Ftplot out=j2.v >/dev/null
	vppen align=lb <j1.v |      vppen  PLOTSIZE  >j11.v
	vppen big=n align=lb <j2.v |vppen  PLOTSIZE  >j22.v
	vppen GRIDSIZE j11.v j22.v |vppen align=lb >FIGDIR/nonlin2.v
	
nonlin3.ascii:  BINDIR/duffing duffing
	BINDIR/duffing omega=9.3 epsilon=40.0 delta=2.0 gamma=10 y1=0.0 y2=1.0 pts=100 tmax=100.00
	cp phase nonlin3.ascii

FIGDIR/nonlin3.v: nonlin3.H 
	Graph <nonlin3.H title='Phase Space' label1='u' label2='v' out=j1.v >/dev/null
	Real <nonlin3.H |Ftplot out=j2.v >/dev/null
	vppen align=lb <j1.v |      vppen  PLOTSIZE  >j11.v
	vppen big=n align=lb <j2.v |vppen  PLOTSIZE  >j22.v
	vppen GRIDSIZE j11.v j22.v |vppen align=lb >FIGDIR/nonlin3.v
	


FIGDIR/chaos.v:  ALLPLOTS
	vppen TOTALGRID  ALLPLOTS |vppen scale=0.4 |vppen align=lb >FIGDIR/chaos.v
	vppen ALLPLOTS scale=2.5 >FIGDIR/chaos.v3

	

%.H: %.ascii
	atoF <%.ascii >%.H
	SETn1.csh %.H


clean&: 
	-/bin/rm phase duffing duffing2
	cake jclean
	-RM_CMD -f *.v *.ascii 


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