#CFLAGS=-fullwarn -O -g3
#CFLAGS=-Wall -O3 -g -DDEBUG -DVERBOSE -pg
CFLAGS=-Wall -O3 -g -DDEBUG -DVERBOSE

LDFLAGS=-lm

SHELL=/bin/sh


vortex: vortex.o runge_kutta_4.o my_malloc.o
	$(CC) $(CFLAGS) -o $@ vortex.o runge_kutta_4.o my_malloc.o $(LDFLAGS)

always:

02-co-different: always vortex
	cd $@ ; ../vortex -c -o 0.1 -t 26.3 -f vinit.dat -v -p
	cd $@ ; ../traj-02.gnuplot

02-co-equal: always vortex
	cd $@ ; ../vortex -c -o 0.1 -t 19.7 -f vinit.dat -v -p
	cd $@ ; ../traj-02.gnuplot

02-counter-different: always vortex
	cd $@ ; ../vortex -c -o 0.10 -t 8.8 -f vinit.dat -v -p
	cd $@ ; ../traj-02.gnuplot

02-counter-equal: always vortex
	cd $@ ; ../vortex -c -o 0.1 -t 10.0 -f vinit.dat -v -p
	cd $@ ; ../traj-02.gnuplot

two-vortices: 02-co-different 02-co-equal 02-counter-different 02-counter-equal

03-co-equal-line: always vortex
	cd $@ ; ../vortex -t 100.0 -f vinit.dat -p -v
	cd $@ ; ../traj-03.gnuplot

03-triangle-equal: always vortex
	cd $@ ; ../vortex -p -o 0.1 -t 13.2 -f vinit.dat -S -v -r -c
	cd $@ ; ../traj-03.gnuplot

04-colinear: always vortex
	cp vortex $@
	#cd $@ ; vortex -t 10.0    -f vinit.dat -l -v; mv vortex.out vortex-0010.out
	#cd $@ ; vortex -t 20.0    -f vinit.dat -l -v; mv vortex.out vortex-0020.out
	#cd $@ ; vortex -t 50.0    -f vinit.dat -l -v; mv vortex.out vortex-0050.out
	#cd $@ ; vortex -t 100.0   -f vinit.dat -l ; mv vortex.out vortex-0100.out
	#cd $@ ; vortex -t 200.0   -f vinit.dat -l ; mv vortex.out vortex-0200.out
	#cd $@ ; vortex -t 500.0   -f vinit.dat -l ; mv vortex.out vortex-0500.out
	#cd $@ ; vortex -t 1000.0  -f vinit.dat -l ; mv vortex.out vortex-1000.out
	#cd $@ ; vortex -t 2000.0  -f vinit.dat -l ; mv vortex.out vortex-2000.out
	#cd $@ ; vortex -t 5000.0  -f vinit.dat -l ; mv vortex.out vortex-5000.out
	#cd $@ ; vortex -t 10000.0 -f vinit.dat -l ; mv vortex.out vortex-10000.out
	#cd $@ ; vortex -t 100.0 -f vinit.dat -S
	#cd $@ ; vortex -t 100000.0 -f vinit.dat -s 0.48 -S -p -o 1000.0
	#cd $@ ; vortex -t 100.0 -f vinit.dat -o 0.2 -S -s 0.48
	cd $@ ; vortex -t 100.0 -f vinit.dat -o 0.2 -p
	#cd $@ ; ../sep.gnuplot

04-square: always vortex
	cd $@ ; ../vortex -t 10000.0 -S -s 2.0 -f v4.dat
	cd $@ ; ../sep.gnuplot

04-rect-1.5-00: always vortex
	#cd $@ ; ../vortex -t 1000.0 -f vinit.dat -S -s 0.85 -p -o 0.2
	cd $@ ; ../vortex -t 500.0 -f vinit.dat -p

04-rect-1.5-01: always vortex
	cd $@ ; ../vortex -t 10000.0 -S -s 2.2 -f v4.dat

04-rect-1.5-02: always vortex
	cd $@ ; ../vortex -t 10000.0 -S -s 2.2 -f v4.dat

04-rect-1.5-03: always vortex
	cd $@ ; ../vortex -t 10000.0 -S -s 2.2 -f v4.dat

04-rect-1.5-04: always vortex
	cd $@ ; ../vortex -t 10000.0 -S -s 2.2 -f v4.dat

04-rect-1.5-05: always vortex
	cd $@ ; ../vortex -t 10000.0 -S -s 2.2 -f v4.dat

04-rect-1.5-06: always vortex
	cd $@ ; ../vortex -t 50000.0 -S -s 2.2 -f v4.dat

04-rect-1.5-07: always vortex
	cd $@ ; ../vortex -t 1000.0 -S -s 2.2 -f v4.dat -l -v

04-rect-1.5-08: always vortex
	cd $@ ; ../vortex -t 20000.0 -S -s 2.2 -f v4.dat

04-rect-1.5-09: always vortex
	cd $@ ; ../vortex -t 10000.0 -S -s 2.2 -f v4.dat

04-rect-1.5-10: always vortex
	cd $@ ; ../vortex -t 1000000.0 -S -s 2.2 -f v4.dat

04-rect-1.5: always vortex
	make 04-rect-1.5-00
	make 04-rect-1.5-01
	make 04-rect-1.5-02
	make 04-rect-1.5-03
	make 04-rect-1.5-04
	make 04-rect-1.5-05
	make 04-rect-1.5-06
	make 04-rect-1.5-07
	make 04-rect-1.5-08
	make 04-rect-1.5-09
	make 04-rect-1.5-10
	04-rect-1.5-sep.gnuplot

04-rect-1.75-00: always vortex
	#cd $@ ; ../vortex -t 10.0 -f vinit.dat -S
	#cd $@ ; ../vortex -t 1000.0 -f vinit.dat -S -s 0.84 -p -o 0.2
	#cd $@ ; ../vortex -t   10.0 -f vinit.dat -l ; mv vortex.out vortex-0010.out
	#cd $@ ; ../vortex -t   20.0 -f vinit.dat -l ; mv vortex.out vortex-0020.out
	#cd $@ ; ../vortex -t   50.0 -f vinit.dat -l ; mv vortex.out vortex-0050.out
	ecd $@ ; ../vortex -t  100.0 -f vinit.dat -l ; mv vortex.out vortex-0100.out
	#cd $@ ; ../vortex -t  200.0 -f vinit.dat -l ; mv vortex.out vortex-0200.out
	#cd $@ ; ../vortex -t  500.0 -f vinit.dat -l ; mv vortex.out vortex-0500.out
	#cd $@ ; ../vortex -t 1000.0 -f vinit.dat -l ; mv vortex.out vortex-1000.out
	cd $@ ; ../vortex -t   100.0 -f vinit.dat -p -o 0.1

04-rect-1.75-01: always vortex
	#cd $@ ; ../vortex -t 2000.0 -f vinit.dat -S -s 0.84 -r -c -C 2.54962 -o 0.2 -p
	cd $@ ; ../vortex -t 100.0 -f vinit.dat -r -c -C 2.54962 -o 0.2 -p

04-rect-1.75-02: always vortex
	#cd $@ ; ../vortex -t 2000.0 -f vinit.dat -S -s 0.84 -r -c -C 2.54962 -o 0.2 -p
	cd $@ ; ../vortex -t 100.0 -f vinit.dat -r -c -C 2.54962 -o 0.2 -p

04-rect-1.75-03: always vortex
	cd $@ ; ../vortex -t 2000.0 -f vinit.dat -S -s 0.84 -r -c -C 2.54962 -o 0.2 -p

04-rect-1.75-04: always vortex
	#cd $@ ; ../vortex -t 100000.0 -f vinit.dat -S -s 0.84 -r -c -C 2.54962 -o 0.2 -p
	cd $@ ; ../vortex -t 100.0 -f vinit.dat -r -c -C 2.54962 -o 0.2 -p

04-rect-1.75: 04-rect-1.75-01 04-rect-1.75-02 04-rect-1.75-03 04-rect-1.75-04

04-rect-2.0-00: always vortex
	#cd $@ ; ../vortex -t 100000.0 -S -s 0.95 -f vinit.dat -v
	#cd $@ ; ../vortex -t   10.0 -f vinit.dat -v -l ; mv vortex.out vortex-0010.out
	#cd $@ ; ../vortex -t   20.0 -f vinit.dat -v -l ; mv vortex.out vortex-0020.out
	#cd $@ ; ../vortex -t   50.0 -f vinit.dat -v -l ; mv vortex.out vortex-0050.out
	#cd $@ ; ../vortex -t  100.0 -f vinit.dat -v -l ; mv vortex.out vortex-0100.out
	#cd $@ ; ../vortex -t  200.0 -f vinit.dat -v -l ; mv vortex.out vortex-0200.out
	#cd $@ ; ../vortex -t  500.0 -f vinit.dat -v -l ; mv vortex.out vortex-0500.out
	#cd $@ ; ../vortex -t 1000.0 -f vinit.dat -v -l ; mv vortex.out vortex-1000.out
	#cd $@ ; ../vortex -t 2000.0 -f vinit.dat -v -l ; mv vortex.out vortex-2000.out
	#cd $@ ; ../vortex -t 5000.0 -f vinit.dat -v -l ; mv vortex.out vortex-5000.out
	cd $@ ; ../vortex -t   200.0 -f vinit.dat -p -o 0.2

04-triangle-a0.001-01: always vortex
	#cd $@ ; ../vortex -t 100.0 -T 0.001 -p -v -x 0 -y 0
	#cd $@ ; ../traj-04.gnuplot
	#cd $@ ; ../tri.gnuplot
	#cd $@ ; ../vortex -t 100.0 -T 0.001 -p -v -x 0 -y 0 -l
	cd $@ ; ../vortex -t 1000.0 -T 0.001 -p -v -x 0 -y 0 -l

04-triangle-a0.001-02: always vortex
	#cd $@ ; ../vortex -t 10000.0 -T 0.001 -v -x 0 -y 0.1
	#cd $@ ; ../tri.gnuplot
	cd $@ ; ../vortex -t 100.0 -T 0.001 -v -x 0 -y 0.1

04-triangle-a0.001-03: always vortex
	cd $@ ; ../vortex -t 5000.0 -T 0.001 -v -x 0 -y 0.2
	cd $@ ; ../tri.gnuplot

04-triangle-a0.001-04: always vortex
	cd $@ ; ../vortex -t 5000.0 -T 0.001 -v -x 0.2 -y -0.3
	cd $@ ; ../tri.gnuplot

04-triangle-a0.001-05: always vortex
	cd $@ ; ../vortex -t 5000.0 -T 0.001 -v -x 0.0 -y -0.5
	cd $@ ; ../tri.gnuplot

04-triangle-a0.001-06: always vortex
	cd $@ ; ../vortex -t 5000.0 -T 0.001 -v -x 1.0 -y 0.0
	cd $@ ; ../tri.gnuplot

04-triangle-a0.001-07: always vortex
	cd $@ ; ../vortex -t 5000.0 -T 0.001 -v -x 0.0 -y 1.2
	cd $@ ; ../tri.gnuplot

04-triangle-a0.001: 04-triangle-a0.001-01 04-triangle-a0.001-02
04-triangle-a0.001: 04-triangle-a0.001-03 04-triangle-a0.001-04
04-triangle-a0.001: 04-triangle-a0.001-05 04-triangle-a0.001-06
04-triangle-a0.001: 04-triangle-a0.001-07

04-triangle-a0.200: always vortex
	cd $@ ; ../vortex -t 20.0 -T 0.200 -p -v
	cd $@ ; ../traj-04.gnuplot
	cd $@ ; ../tri.gnuplot

04-triangle-a0.400: always vortex
	cd $@ ; ../vortex -t 20.0 -T 0.400 -p -v
	cd $@ ; ../traj-04.gnuplot
	cd $@ ; ../tri.gnuplot

04-triangle-a0.400-00: always vortex
	cd $@ ; ../vortex -t 500000.0 -T 0.400 -v -p -o 1000.0
	cd $@ ; ../traj-04.gnuplot
	cd $@ ; ../tri.gnuplot

04-triangle-a0.600: always vortex
	cd $@ ; ../vortex -t 20.0 -T 0.600 -p -v
	cd $@ ; ../traj-04.gnuplot
	cd $@ ; ../tri.gnuplot

04-triangle-a0.600-00: always vortex
	cd $@ ; ../vortex -t 2000.0 -T 0.600 -v
	#cd $@ ; ../traj-04.gnuplot
	cd $@ ; ../tri.gnuplot

04-triangle-a0.700: always vortex
	cd $@ ; ../vortex -t 20.0 -T 0.700 -p -v
	cd $@ ; ../traj-04.gnuplot
	cd $@ ; ../tri.gnuplot

04-triangle-a0.700-00: always vortex
	cp vortex $@
	#cd $@ ; vortex -t 20000.0 -T 0.700 -v -o 1000.0 -p
	#cd $@ ; vortex -t 21001.0 -T 0.700 -f vinit-19000.dat -v -o 1000.0 -p
	#cd $@ ; vortex -t 40001.0 -T 0.700 -f vinit-29000.dat -v -o 1000.0 -p
	#cd $@ ; ../traj-04.gnuplot
	#cd $@ ; ../tri.gnuplot
	cd $@ ; vortex -t   10.0 -T 0.700 -v -l;mv vortex.out vortex-0010.out
	cd $@ ; vortex -t   20.0 -T 0.700 -v -l;mv vortex.out vortex-0020.out
	cd $@ ; vortex -t   50.0 -T 0.700 -v -l;mv vortex.out vortex-0050.out
	cd $@ ; vortex -t  100.0 -T 0.700 -v -l;mv vortex.out vortex-0100.out
	cd $@ ; vortex -t  200.0 -T 0.700 -v -l;mv vortex.out vortex-0200.out
	cd $@ ; vortex -t  500.0 -T 0.700 -v -l;mv vortex.out vortex-0500.out
	cd $@ ; vortex -t 1000.0 -T 0.700 -v -l;mv vortex.out vortex-1000.out
	cd $@ ; vortex -t 2000.0 -T 0.700 -v -l;mv vortex.out vortex-2000.out

04-triangle-a0.700-01: always vortex
	cp vortex $@
	cd $@ ; vortex -t 2001.0 -T 0.700 -v -o 1000.0 -p -x 0 -y 1.0
	cd $@ ; ../tri.gnuplot

04-triangle-a0.700-02: always vortex
	cp vortex $@
	cd $@ ; vortex -t 2001.0 -T 0.700 -v -o 1000.0 -p -x 1.2 -y 0.0
	cd $@ ; ../tri.gnuplot

04-triangle-a0.700-03: always vortex
	cp vortex $@
	cd $@ ; vortex -t 2001.0 -T 0.700 -v -o 1000.0 -p -x 0.9 -y -0.7
	cd $@ ; ../tri.gnuplot

04-triangle-a0.700-04: always vortex
	cp vortex $@
	cd $@ ; vortex -t 5001.0 -T 0.700 -v -o 1000.0 -p -x 0.8 -y -0.91
	cd $@ ; ../tri.gnuplot

04-triangle-a0.700-05: always vortex
	cp vortex $@
	cd $@ ; vortex -t 5001.0 -T 0.700 -v -o 1000.0 -p -x 0.9 -y -0.7
	cd $@ ; ../tri.gnuplot

04-triangle-a0.700: 04-triangle-a0.700-00 04-triangle-a0.700-01
04-triangle-a0.700: 04-triangle-a0.700-02 04-triangle-a0.700-03
04-triangle-a0.700: 04-triangle-a0.700-04 04-triangle-a0.700-05

04-triangle-a0.800: always vortex
	cd $@ ; ../vortex -t 20.0 -T 0.800 -p -v
	cd $@ ; ../traj-04.gnuplot
	cd $@ ; ../tri.gnuplot

04-triangle-a1.000: always vortex
	cd $@ ; ../vortex -t 20.0 -T 1.000 -p -v
	cd $@ ; ../traj-04.gnuplot
	cd $@ ; ../tri.gnuplot

04-triangles: 04-triangle-a0.200 04-triangle-a0.400 04-triangle-a0.600
04-triangles: 04-triangle-a0.800 04-triangle-a1.000

04-line-passive: always vortex
	cd $@ ; ../vortex -t 1000.0 -S -s 1.5 -f v4.dat
	cd $@ ; ../traj-04.gnuplot
	cd $@ ; ../traj-04-v4.gnuplot
	cd $@ ; ../sep.gnuplot

05-line: always vortex
	cd $@ ; ../vortex -t   10.0 -f vinit.dat -v -l; mv vortex.out vortex-0010.out
	#cd $@ ; ../vortex -t   20.0 -f vinit.dat -v -l; mv vortex.out vortex-0020.out
	#cd $@ ; ../vortex -t   50.0 -f vinit.dat -v -l; mv vortex.out vortex-0050.out
	#cd $@ ; ../vortex -t  100.0 -f vinit.dat -v -l; mv vortex.out vortex-0100.out
	#cd $@ ; ../vortex -t  200.0 -f vinit.dat -v -l; mv vortex.out vortex-0200.out
	#cd $@ ; ../vortex -t  500.0 -f vinit.dat -v -l; mv vortex.out vortex-0500.out
	#cd $@ ; ../vortex -t 1000.0 -f vinit.dat -v -l; mv vortex.out vortex-1000.out
	#cd $@ ; ../vortex -t 2000.0 -f vinit.dat -v -l; mv vortex.out vortex-2000.out
	cd $@ ; ../vortex -t  100.0 -f vinit.dat -v -p -o 0.2

10-line: always vortex
	#cd $@ ; ../vortex -t   5.0 -f vinit.dat -v -l; mv vortex.out vortex-0005.out
	#cd $@ ; ../vortex -t  10.0 -f vinit.dat -v -l; mv vortex.out vortex-0010.out
	#cd $@ ; ../vortex -t  20.0 -f vinit.dat -v -l; mv vortex.out vortex-0020.out
	#cd $@ ; ../vortex -t  50.0 -f vinit.dat -v -l; mv vortex.out vortex-0050.out
	#cd $@ ; ../vortex -t 100.0 -f vinit.dat -v -l; mv vortex.out vortex-0100.out
	#cd $@ ; ../vortex -t 200.0 -f vinit.dat -v -l; mv vortex.out vortex-0200.out
	#cd $@ ; ../vortex -t 500.0 -f vinit.dat -v -l; mv vortex.out vortex-0500.out
	#cd $@ ; ../vortex -t 1000.0 -f vinit.dat -v -l; mv vortex.out vortex-1000.out
	cd $@ ; ../vortex -t  10.0 -f vinit.dat -p

SUBDIRS=

DIST=*.[ch] *.gnuplot trajAnim Makefile paper/paper.tex

BACKDIR=`date +%y%b%d | tr '[A-Z]' '[a-z]'`

DISTDIR=pv-$(BACKDIR)

DISTFILE=$(DISTDIR).tar

backup:
	mkdir $(BACKDIR)
	tar cvfh $(BACKDIR)/$(DISTFILE) $(DIST) $(SUBDIRS)
	#tar cvf $(BACKDIR)/$(DISTFILE) $(DIST) $(SUBDIRS)
	( cd $(BACKDIR) ; tar xvf $(DISTFILE) && rm $(DISTFILE) )

dist: backup
	@if [ -d $(DISTDIR) ]; then echo "$(DISTDIR) already exists";exit 1;fi
	mv $(BACKDIR) $(DISTDIR)
	tar cvf $(DISTFILE) $(DISTDIR)
	gzip $(DISTFILE)

ctags:
	ctags *.[ch]

clean: ctags
	-rm -f *.o

veryclean: clean
	-rm -f vortex