F ?= 10
E ?= 100
geom ?= elastic

TMPFILE := $(shell mktemp)

all: meshing 
	$(eval TGT := $(shell ls *_$(geom)_$(F)g$(E)kPa.unv))
	$(eval TGT := $(TGT:%.unv=%.xml))
	make $(TGT)
	mkdir -p $(geom)/$(F)g/$(E)kPa
	mv *_$(geom)_$(F)g$(E)kPa.xml $(geom)/$(F)g/$(E)kPa

meshing: 
	salome start -t --ns-port-log=$(TMPFILE)
	salome shell -p `cat $(TMPFILE)` -- python -tt geometries_$(geom).py [$(F),$(E)]
	salome kill `cat $(TMPFILE)`
	rm -f $(TMPFILE)

clean: 
	rm -f *.msh *.unv

distclean: clean
	rm -f *.xml

%.msh: %.unv
	gmsh -3 $<

%.xml: %.msh
	meshio-convert $< $(<:.msh=.xml)
	rm -f $< $(<:.msh=_gmsh:geometrical.xml)
	rm -f $< $(<:.msh=_gmsh:physical.xml)

print: 
	@echo 'geom = elastic, 0 < F < 100, E>=50'

figures:
	for j in *.unv; do gmsh -3 $j -o "$(basename "$j" .unv).vtk"; done;
	mv *.vtk ../Comp_domain_figures/
