PROJECT=$(basename $(wildcard *.gnuplot))

SIMPLESOLV=simplesolv~/Work/Research/2018+bose_gas/cmp/simplesolv/source/src

all: $(addsuffix .pdf, $(PROJECT))

define gnuplot2pdf
	gnuplot $(1).gnuplot > $(1).tikz.tex
	pdflatex -jobname $(1) -file-line-error $(1).tikz.tex
endef

define bigeq_run
	julia -p 16 $(SIMPLESOLV)/main.jl -p "eq=bigeq;N=12;P=8;J=10;maxiter=100;v_a=8;rho=1e-$(1);minlrho_init=-6;nlrho_init=100;window_L=100;kmin=0;kmax=8;nk=100" -M anyeq 2pt_fourier > 2pt_fourier_$(1).dat
endef
define bigeq_run_fine
	julia -p 16 $(SIMPLESOLV)/main.jl -p "eq=bigeq;N=12;P=8;J=10;maxiter=100;v_a=8;rho=1e-$(1);minlrho_init=-6;nlrho_init=100;window_L=100;kmin=0;kmax=8;nk=150" -M anyeq 2pt_fourier > 2pt_fourier_$(1).dat
endef
define bigeq_run_plus
	julia -p 16 $(SIMPLESOLV)/main.jl -p "eq=bigeq;N=18;P=12;J=15;maxiter=100;v_a=8;rho=1e-$(1);minlrho_init=-6;nlrho_init=100;window_L=50;kmin=0;kmax=8;nk=100" -M anyeq 2pt_fourier > 2pt_fourier_$(1).dat
endef

2pt_fourier_5.dat:
	$(call bigeq_run,5)
2pt_fourier_4.dat:
	$(call bigeq_run,4)
2pt_fourier_3.dat:
	$(call bigeq_run,3)
2pt_fourier_2.dat:
	$(call bigeq_run_fine,2)
2pt_fourier_1.dat:
	$(call bigeq_run_plus,1)
2pt_fourier_0.3.dat:
	julia -p 16 $(SIMPLESOLV)/main.jl -p "eq=bigeq;N=18;P=12;J=15;maxiter=100;v_a=8;rho=0.3;minlrho_init=-6;nlrho_init=100;window_L=25;kmin=0;kmax=8;nk=100" -M anyeq 2pt_fourier > 2pt_fourier_0.3.dat

DATS=2pt_fourier_5.dat 2pt_fourier_4.dat 2pt_fourier_3.dat 2pt_fourier_2.dat 2pt_fourier_1.dat 2pt_fourier_0.3.dat

2pt_fourier_peak.pdf: $(DATS)
	$(call gnuplot2pdf,$(basename $@))
2pt_fourier_full.pdf: $(DATS)
	$(call gnuplot2pdf,$(basename $@))


clean-aux:
	rm -f $(addsuffix .tikz.tex, $(PROJECT))
	rm -f *.aux *.log

clean-dat:
	rm -f $(DATS)

clean-pdf:
	rm -f $(addsuffix .pdf, $(PROJECT))

clean: clean-aux clean-pdf