PRJ=ODE_sensitivity-0.1
CPP=g++
OPT=-O3
NEWMAT=/home/kbriggs/newmat11/ # set path to your newmat here (http://www.robertnz.net/nm_intro.htm)

bruss: bruss.o dopri5a.o
	$(CPP) $(OPT) bruss.o -L ${NEWMAT} dopri5a.o -lnewmat -o bruss

bruss.o: bruss.cc ode.h ad.h
	$(CPP)  $(OPT) -c bruss.cc -I ${NEWMAT}

dopri5a.o: dopri5a.cc dopri5a.h
	$(CPP)  $(OPT) -c dopri5a.cc 

plot: bruss
	./bruss >| bruss.plt
	colex 1 2 < bruss.plt >| bruss1.plt
	colex 1 3 < bruss.plt >| bruss2.plt
	graph -Tps -C -g3 --font-name="times-italic" --font-size="0.04" -X"B" -Y"d(solution)/dB" -L"Brusselator parameter sensitivity" --line-width=0.005 -m1 bruss1.plt -m3 bruss2.plt | epstopdf --filter >| bruss.pdf
	#graph -Tpng -C -g3  --font-size="0.05" -X"B" -Y"d(solution)/dB" -L"Brusselator parameter sensitivity" bruss1.plt bruss2.plt >| bruss.png
	convert bruss.pdf bruss.png
	@echo
	@echo "acroread bruss.pdf"

tar:
	(cd ..; tar zcvf tarfiles/${PRJ}.tgz ${PRJ}/README ${PRJ}/*h ${PRJ}/*c ${PRJ}/Makefile)

www:
	make tar
	cp ODE_sensitivity.html ../website-1.1/ODE_sensitivity.txt
	(cd ../website-1.1; make && ./upload www/ODE_sensitivity.html)
