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


XSIL_FILES = Nlevels_no_dopler_with_z_4wm.xsil
M_FILES = $(patsubst %.xsil,%.m,$(XSIL_FILES))
GNUPLOT_FILES = $(wildcard *.gp)

SCRIPTS_DIR = .

XSIL2GRAPHICS = xsil2graphics

# fast light
#PARAMS =  --delta1=0 --delta2=0 --delta3=0 --E1o=1.9e7 --E2o=3.1e5 --E3o=3.8e7 --E4o=1e1
#PARAMS =  --delta1=0 --delta2=0 --delta3=0 --E1o=0.1e7 --E2o=1e4 --E3o=.3e7 --E4o=0 --Lcell=1.5e-2 --Temperature=.0001
PARAMS =  \
	  --Ndens=1e15 \
	  --Lcell=10.0e-2 \
	  --Temperature=5 \
	  --Pwidth=0.4e-6 \
	  --delta1=0 --delta2=0 --delta3=0 \
	  --E1o=2e7 --E2o=1e2 --E3o=4e7 --E4o=1e0


include $(PARAMS_FILE)


# 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 fig

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

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

fig: pp_I2.stamp
	$(MAKE) -f $(SCRIPTS_DIR)/Makefile.fig $@


pretty_plots: pp_I2.stamp $(GNUPLOT_FILES)
	gnuplot $(SCRIPTS_DIR)/plot_fields_propagation_I2.gp
	gnuplot $(SCRIPTS_DIR)/plot_fields_propagation_I4.gp

plot: pp_I2.stamp

pp_I2.stamp: $(XSIL_FILES) $(M_FILES)
	octave -q $(SCRIPTS_DIR)/pp_I2.m

clean:
	rm -f $(M_FILES) $(XSIL_FILES) *.dat octave-core 
	rm -f pp_I2.stamp
	$(MAKE) -f $(SCRIPTS_DIR)/Makefile.fig $@

real_clean: clean
	$(MAKE) -f $(SCRIPTS_DIR)/Makefile.fig $@
	rm -f exact_analysis_execution_time.txt

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