###  -*- make -*-
### This makefile can be used to build and run the XMDS examples


XSIL_FILES = Nlevels_no_dopler_with_z_4wm_with_perturbations.xsil Nlevels_no_dopler_with_z_4wm.xsil
M_FILES = $(patsubst %.xsil,%.m,$(XSIL_FILES))

XSIL2GRAPHICS = xsil2graphics

# fast light
#PARAMS =  --delta1=0 --delta2=0 --delta3=0 --E1o=1.9e7 --E2o=3.1e5 --E3o=3.8e7 --E4o=6.3e4
# slow light EIT
PARAMS =  --delta1=0 --delta2=0 --delta3=0 --E1o=1.9e7 --E2o=3.1e5 --E3o=0 --E4o=0

all: $(XSIL_FILES) Nlevels_no_dopler_with_z_4wm.xsil $(M_FILES) plot

Nlevels_no_dopler_with_z_4wm_with_perturbations.xsil: ../Nlevels_no_dopler_with_z_4wm_with_perturbations/Nlevels_no_dopler_with_z_4wm_with_perturbations.run
	$< $(PARAMS) | grep "Time elapsed for simulation is:" > perturbative_analysis_execution_time.txt

Nlevels_no_dopler_with_z_4wm.xsil: ../Nlevels_no_dopler_with_z_4wm/Nlevels_no_dopler_with_z_4wm.run
	$< $(PARAMS) | grep "Time elapsed for simulation is:" > exact_analysis_execution_time.txt

%.m: %.xsil
	$(XSIL2GRAPHICS) $<

plot: $(M_FILES)
	octave pp.m

clean:
	rm -f $(CC_FILES) $(RUN_FILES) $(M_FILES) $(XSIL_FILES) *.wisdom.fftw3 *.dat octave-core *.wisdom *.pdf
	rm -f $(png_targets)

real_clean: clean
	rm -f *.eps

eps_targets  = $(wildcard *.eps)
pdf_targets  = $(eps_targets:%.eps=%.pdf)	
png_targets  = $(pdf_targets:%.pdf=%.png)

pdf: $(pdf_targets) 

$(pdf_targets): %.pdf : %.eps
	cat $< | ps2eps -B > __tt.eps
	epspdf __tt.eps $@
	rm -f __tt.eps
	#ps2eps -B $< | epspdf $< $@

png: pdf $(png_targets) 

$(png_targets): %.png : %.pdf
	convert -density 300 $< $@

.PRECIOUS: %.run %.xsil %.m
.PHONY: all clean