Initial commitv0.0

3 files changed, 146 insertions, 0 deletions

diff --git a/figs/energy_relative_error.fig/Makefile b/figs/energy_relative_error.fig/Makefile
new file mode 100644
index 0000000..8c75f0f
--- /dev/null
+++ b/figs/energy_relative_error.fig/Makefile
@@ -0,0 +1,71 @@
+PROJECTNAME=cmp_energy16 cmp_energy1
+
+SIMPLEQ=simpleq
+
+DATS= simpleq16.dat mueq16.dat bigeq16.dat simpleq1.dat mueq1.dat bigeq1.dat
+PDFS=$(addsuffix .pdf, $(PROJECTNAME))
+TEXS=$(addsuffix .tikz.tex, $(PROJECTNAME))
+
+all: $(PDFS)
+
+$(PDFS): $(DATS)
+	gnuplot $(patsubst %.pdf, %.gnuplot, $@) > $(patsubst %.pdf, %.tikz.tex, $@)
+	pdflatex -jobname $(basename $@) -file-line-error $(patsubst %.pdf, %.tikz.tex, $@)
+
+bigeq16.dat:
+	for rho in 1e-3 5e-3 1e-2 2e-2 1e-1 1e+0; do \
+	  echo -n $$rho "  " >> $@-tmp ;\
+	  julia $(SIMPLEQ)/main.jl -p "eq=bigeq;N=12;P=8;J=10;maxiter=100;v_a=16;rho=$$rho;minlrho=-6;nlrho=50" -M medeq energy >> $@-tmp ;\
+	done
+	paste $@-tmp 2020-08-27+16.dat > $@
+	rm $@-tmp
+simpleq16.dat:
+	for lrho in -3 -2.3010299956639813 -2 -1.6989700043360187 -1 0; do \
+	  julia $(SIMPLEQ)/main.jl -p "tolerance=1e-11;order=100;maxiter=100;v_a=16;minlrho=$$lrho;nlrho=1" -M simpleq energy_rho >> $@-tmp ;\
+	done
+	paste $@-tmp 2020-08-27+16.dat > $@
+	rm $@-tmp
+mueq16.dat:
+	for lrho in -3 -2.3010299956639813 -2 -1.6989700043360187 -1 0; do \
+	  julia $(SIMPLEQ)/main.jl -p "tolerance=1e-11;order=100;maxiter=100;v_a=16;minlrho=$$lrho;nlrho=1" -M mueq energy_rho >> $@-tmp ;\
+	done
+	paste $@-tmp 2020-08-27+16.dat > $@
+	rm $@-tmp
+
+bigeq1.dat:
+	for rho in 1e-6 1e-4 1e-3 5e-3 1e-2 1.5e-2 2e-2 5e-2 1e-1 1.0 1e1 5e1; do \
+	  echo -n $$rho "  " >> $@-tmp ;\
+	  julia $(SIMPLEQ)/main.jl -p "N=12;P=8;J=10;maxiter=100;v_a=1;rho=$$rho;minlrho=-6;nlrho=50;eq=bigeq" -M medeq energy >> $@-tmp ;\
+	done
+	paste $@-tmp 2020-10-15+energy.dat > $@
+	rm $@-tmp
+simpleq1.dat:
+	for lrho in -6 -4 -3 -2.3010299956639813 -2 -1.8239087409443189 -1.6989700043360187 -1.3010299956639813 -1 0 1 1.6989700043360187; do\
+	  julia $(SIMPLEQ)/main.jl -p "tolerance=1e-11;order=100;maxiter=100;v_a=1;minlrho=$$lrho;nlrho=1" -M simpleq energy_rho >> $@-tmp ;\
+	done
+	paste $@-tmp 2020-10-15+energy.dat > $@
+	rm $@-tmp
+mueq1.dat:
+	for lrho in -6 -4 -3 -2.3010299956639813 -2 -1.8239087409443189 -1.6989700043360187 -1.3010299956639813 -1 0 1 1.6989700043360187; do\
+	  julia $(SIMPLEQ)/main.jl -p "tolerance=1e-11;order=100;maxiter=100;v_a=1;minlrho=$$lrho;nlrho=1" -M mueq energy_rho >> $@-tmp ;\
+	done
+	paste $@-tmp 2020-10-15+energy.dat > $@
+	rm $@-tmp
+
+
+install: $(PDFS)
+	cp $^ $(INSTALLDIR)/
+
+clean-aux:
+	rm -f $(addsuffix .tikz.tex, $(PROJECTNAME))
+	rm -f $(addsuffix .aux, $(PROJECTNAME))
+	rm -f $(addsuffix .log, $(PROJECTNAME))
+
+clean-dat:
+	rm -f $(DATS)
+	rm -f Abar.dat
+
+clean-tex:
+	rm -f $(PDFS)
+
+clean: clean-aux clean-tex
diff --git a/figs/energy_relative_error.fig/cmp_energy1.gnuplot b/figs/energy_relative_error.fig/cmp_energy1.gnuplot
new file mode 100644
index 0000000..977d249
--- /dev/null
+++ b/figs/energy_relative_error.fig/cmp_energy1.gnuplot
@@ -0,0 +1,38 @@
+set title "$v(\\mathbf x)=e^{-|\\mathbf x|}$"
+set ylabel "energy relative error" offset +3
+set xlabel "$\\rho$"
+
+set xtics 1e-6, 100, 100
+set xtics add ("$10^{-6}$" 0.000001, "$10^{-4}$" 0.0001, "$10^{-2}$" 0.01, "$1$" 1.0, "$10^2$" 100)
+unset mxtics
+set xrange [1e-7:1e2]
+
+set ytics 1e-5, 10, 1e-1
+set ytics add ("$10^{-5}$" 1e-5, "$10^{-4}$" 1e-4, "$10^{-3}$" 1e-3, "$10^{-2}$" 1e-2, "$10^{-1}$" 1e-1)
+set yrange [1e-6:1e-1]
+
+# default output canvas size: 12.5cm x 8.75cm
+set term lua tikz size 8,6 standalone
+
+set key bottom left box linetype rgbcolor"#999999" width 0.3 height 0.3 spacing 1.3
+
+
+# set linestyle
+set style line 1 linetype rgbcolor "#4169E1" linewidth 2 dashtype "."
+set style line 2 linetype rgbcolor "#DC143C" linewidth 2
+set style line 3 linetype rgbcolor "#32CD32" linewidth 2 dashtype "-"
+set style line 4 linetype rgbcolor "#4B0082" linewidth 2
+set style line 5 linetype rgbcolor "#DAA520" linewidth 2
+
+set pointsize 1
+
+max(x,y)=(x>y ? x : y)
+
+set logscale xy
+
+plot \
+ "simpleq1.dat" using 1:($3<1e-5 ? max(abs($2-$5)/$5,abs($2-$6)/$6) : 1/0) every ::1 with lp ls 1 title "simple" ,\
+ "mueq1.dat" using 1:($3<1e-5 ? max(abs($2-$5)/$5,abs($2-$6)/$6) : 1/0) every ::1 with lp ls 3 title "med" ,\
+ "bigeq1.dat" using 1:($3<1e-5 ? max(abs($2-$5)/$5,abs($2-$6)/$6) : 1/0) every ::1 with lp ls 5 title "big" ,\
+ "2020-10-15+energy.dat" using 1:(max(max(abs($6-$2)/$2,abs($6-$3)/$3),max(abs($7-$2)/$2,abs($7-$3)/$3))) every ::1 with lp ls 2 title "BDJ"
+
diff --git a/figs/energy_relative_error.fig/cmp_energy16.gnuplot b/figs/energy_relative_error.fig/cmp_energy16.gnuplot
new file mode 100644
index 0000000..b2ffb4c
--- /dev/null
+++ b/figs/energy_relative_error.fig/cmp_energy16.gnuplot
@@ -0,0 +1,37 @@
+set title "$v(\\mathbf x)=16e^{-|\\mathbf x|}$"
+set ylabel "energy relative error" offset +3
+set xlabel "$\\rho$"
+
+set xtics 1e-6, 100, 100
+set xtics add ("$10^{-6}$" 0.000001, "$10^{-4}$" 0.0001, "$10^{-2}$" 0.01, "$1$" 1.0, "$10^2$" 100)
+unset mxtics
+set xrange [0.000001:100]
+
+set ytics 1e-5, 10, 1
+set ytics add ("$10^{-5}$" 1e-5, "$10^{-4}$" 1e-4, "$10^{-3}$" 1e-3, "$10^{-2}$" 1e-2, "$10^{-1}$" 1e-1, "$1$" 1)
+set yrange [1e-5:1]
+
+# default output canvas size: 12.5cm x 8.75cm
+set term lua tikz size 8,6 standalone
+
+set key bottom left box linetype rgbcolor"#999999" width 0.3 height 0.3 spacing 1.3
+
+
+# set linestyle
+set style line 1 linetype rgbcolor "#4169E1" linewidth 2 dashtype "."
+set style line 2 linetype rgbcolor "#DC143C" linewidth 2
+set style line 3 linetype rgbcolor "#32CD32" linewidth 2 dashtype "-"
+set style line 4 linetype rgbcolor "#4B0082" linewidth 2
+set style line 5 linetype rgbcolor "#DAA520" linewidth 2
+
+set pointsize 1
+
+set logscale xy
+
+max(x,y)=(x>y ? x : y)
+
+plot \
+ "simpleq16.dat" using 1:($3<1e-5 ? max(abs($2-$5)/$5,abs($2-$6)/$6) : 1/0) with lp ls 1 title "simple" ,\
+ "mueq16.dat" using 1:($3<1e-5 ? max(abs($2-$5)/$5,abs($2-$6)/$6) : 1/0) with lp ls 3 title "med" ,\
+ "bigeq16.dat" using 1:($3<1e-5 ? max(abs($2-$5)/$5,abs($2-$6)/$6) : 1/0) with lp ls 5 title "big" ,\
+ "2020-08-27+16.dat" using 1:(max(max(abs($6-$2)/$2,abs($6-$3)/$3),max(abs($7-$2)/$2,abs($7-$3)/$3))) with lp ls 2 title "BDJ"