SDKDIR=$(HOME)/.dotnet/sdk/6.0.102# or whatever your sdk is
CSCPATH=$(SDKDIR)/Roslyn/bincore/csc.dll
NETSTANDARD=$(SDKDIR)/ref/netstandard.dll
CONFIG=$(SDKDIR)/vstest.console.runtimeconfig.json
DOTNET=DOTNET_CLI_TELEMETRY_OPTOUT=1; dotnet
CSC=$(DOTNET) $(CSCPATH) -reference:$(NETSTANDARD)
RUN=$(DOTNET) exec --runtimeconfig $(CONFIG)
# export DOTNET_CLI_TELEMETRY_OPTOUT=1
CSC = mcs
RUN = mono

all:Psi.svg E0.rmax.svg E1.rmax.svg E0.dr.svg

E0.rmax.svg:out.E0.rmax.txt Makefile
	echo '\
	set term svg;\
	set out "$@";\
	set title "Hydrogen, matrix, E_0 convergence with r_{max}";\
	set xlabel "r_{max}, Bohr radius";\
	set ylabel "E_{0}, Hartree";\
	plot [][-0.6:]\
	 "$<" using 1:3 with linespoints pointtype 5 title "calculated"\
	,-0.5 with lines title "exact"\
	'|gnuplot

E1.rmax.svg:out.E1.rmax.txt Makefile
	echo '\
	set term svg;\
	set out "$@";\
	set title "Hydrogen, matrix, E_1 convergence with r_{max}";\
	set xlabel "r_{max}, Bohr radius";\
	set ylabel "E_{1}, Hartree";\
	plot []\
	 "$<" using 1:4 with linespoints pointtype 5 title "calculated"\
	,-0.125 with lines title "exact"\
	'|gnuplot

E0.dr.svg:out.E0.dr.txt Makefile
	echo '\
	set term svg;\
	set out "$@";\
	set title "Hydrogen, matrix, E_0 convergence with dr";\
	set xlabel "dr, Bohr radius";\
	set ylabel "E_{0}, Hartree";\
	set key left;\
	plot []\
	 "$<" using 2:3 with linespoints pointtype 5 title "calculated"\
	,-0.5 with lines title "exact (-0.5 Hartree)"\
	'|gnuplot

out.E0.dr.txt:main.exe
	>$@
	for dr in 0.02 0.05 0.1 0.15 0.2 0.25 0.3; do \
	echo $$dr; mono main.exe 6 $$dr 1>/dev/null 2>>$@ \
	; done

out.E1.rmax.txt:main.exe
	>$@
	for r in 7 8 9 10 11 12 13; do \
	echo $$r; mono main.exe $$r 0.1 1>/dev/null 2>>$@ \
	; done

out.E0.rmax.txt:main.exe
	>$@
	for r in 2 3 4 5 6; do \
	mono main.exe $$r 0.1 1>/dev/null 2>>$@ \
	; done

Psi.svg:out.data.txt Makefile
	echo '\
	set term svg;\
	set out "$@";\
	set xlabel "r, Bohr radius";\
	set ylabel "f(r), arb. u.";\
	set title "Hydrogen atom ground state radial wave-function";\
	plot\
	 "$<" index 1 using 1:2 with points pointtype 7 pointsize 0.5 title "calculated"\
	,"$<" index 1 using 1:3 with lines title "exact"\
	'| gnuplot

out.data.txt:main.exe
	$(RUN) $< > $@

main.exe:main.cs matlib.dll
	$(CSC) -reference:matlib.dll $<

matlib.dll: \
	../../matlib/matrix/vector.cs \
	../../matlib/matrix/matrix.cs \
	../../matlib/jacobi-eigenvalue/EVD.cs \
	../../matlib/jacobi-eigenvalue/jacobi.cs
	$(CSC) -target:library -out:./$@ $^

clean:
	$(RM) *.dll *.exe *.svg out*

#%.exe: %.cs
#	$(CSC) -t:exe $(addprefix -r:,$(filter %.dll,$^)) $(filter %.cs,$^)
#
#%.dll: %.cs
#	$(CSC) -t:library -out:./%.dll $(addprefix -r:,$(filter %.dll,$^)) $(filter %.cs,$^)