Ad presented at Rutgers University on 2020-12-16HEADv1.0master

5 files changed, 89 insertions, 0 deletions

diff --git a/figs/relative_energy.fig/2020-08-27+2.dat b/figs/relative_energy.fig/2020-08-27+2.dat
new file mode 100644
index 0000000..ce0fd04
--- /dev/null
+++ b/figs/relative_energy.fig/2020-08-27+2.dat
@@ -0,0 +1,5 @@
+ 5e-3      0.09074   0.850    0.09086          0.851  
+ 1e-2      0.19640   0.836    0.19657           0.84  
+ 2e-2       0.4210    0.83     0.4214           0.83  
+ 1e-1       2.3477    0.84      2.348           0.86  
+ 1e+0        24.82    0.91      24.82           0.91  
diff --git a/figs/relative_energy.fig/Makefile b/figs/relative_energy.fig/Makefile
new file mode 100644
index 0000000..065a2d0
--- /dev/null
+++ b/figs/relative_energy.fig/Makefile
@@ -0,0 +1,43 @@
+PROJECTNAME=compare
+
+SIMPLEQ=~/Work/Research/2018+bose_gas/cmp/simpleq
+
+DATS=simpleq2.dat bigeq2.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, $@)
+
+bigeq2.dat:
+	for rho in 0.005 0.01 0.02 0.1 1; do \
+	  echo -n $$rho "  " >> $@-tmp ;\
+	  julia $(SIMPLEQ)/main.jl -p "eq=bigeq;N=8;P=9;J=10;maxiter=100;v_a=2;rho=$$rho;minlrho=-6;nlrho=50" -M medeq energy >> $@-tmp ;\
+	done
+	paste $@-tmp 2020-08-27+2.dat > $@
+	rm $@-tmp
+simpleq2.dat:
+	for lrho in -2.301029995 -2 -1.698970004 -1 0; do \
+	  julia $(SIMPLEQ)/main.jl -p "tolerance=1e-11;order=100;maxiter=100;v_a=2;minlrho=$$lrho;nlrho=1" -M simpleq energy_rho >> $@-tmp ;\
+	done
+	paste $@-tmp 2020-08-27+2.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)
+
+clean-tex:
+	rm -f $(PDFS)
+
+clean: clean-aux clean-tex
diff --git a/figs/relative_energy.fig/bigeq2.dat b/figs/relative_energy.fig/bigeq2.dat
new file mode 100644
index 0000000..d976640
--- /dev/null
+++ b/figs/relative_energy.fig/bigeq2.dat
@@ -0,0 +1,5 @@
+0.005    9.077274334771057e-02  4.210230174384609e-15	 5e-3      0.09074   0.850    0.09086          0.851  
+0.01    1.963333772412874e-01  5.103464326744255e-15	 1e-2      0.19640   0.836    0.19657           0.84  
+0.02    4.210128326609409e-01  4.686441858629907e-15	 2e-2       0.4210    0.83     0.4214           0.83  
+0.1    2.347849641433116e+00  1.470639915178939e-15	 1e-1       2.3477    0.84      2.348           0.86  
+1    2.481579120915402e+01  3.015163776163534e-15	 1e+0        24.82    0.91      24.82           0.91  
diff --git a/figs/relative_energy.fig/compare.gnuplot b/figs/relative_energy.fig/compare.gnuplot
new file mode 100644
index 0000000..7461923
--- /dev/null
+++ b/figs/relative_energy.fig/compare.gnuplot
@@ -0,0 +1,31 @@
+set ylabel "energy relative error"
+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]
+
+# default output canvas size: 12.5cm x 8.75cm
+set term lua tikz size 8,6 standalone
+
+set key off
+
+
+# set linestyle
+set style line 1 linetype rgbcolor "#4169E1" linewidth 2
+set style line 2 linetype rgbcolor "#DC143C" linewidth 2
+set style line 3 linetype rgbcolor "#32CD32" linewidth 2
+set style line 4 linetype rgbcolor "#4B0082" linewidth 2
+set style line 5 linetype rgbcolor "#DAA520" linewidth 2
+
+set pointsize 1
+
+set logscale xy
+
+plot \
+ "simpleq2.dat" using 1:($3<1e-5 ? abs($2-$5)/$5 : 1/0) with lp ls 1 ,\
+ "bigeq2.dat" using 1:($3<1e-5 ? abs($2-$5)/$5 : 1/0) with lp ls 4 ,\
+ "2020-08-27+2.dat" using 1:(abs($4-$2)/$2) with lp ls 2
+
diff --git a/figs/relative_energy.fig/simpleq2.dat b/figs/relative_energy.fig/simpleq2.dat
new file mode 100644
index 0000000..08240cd
--- /dev/null
+++ b/figs/relative_energy.fig/simpleq2.dat
@@ -0,0 +1,5 @@
+ 5.000000007644367e-03  9.849857484186206e-02  8.221496065916031e-13	 5e-3      0.09074   0.850    0.09086          0.851  
+ 1.000000000000000e-02  2.133766483043713e-01  2.425157716108297e-12	 1e-2      0.19640   0.836    0.19657           0.84  
+ 2.000000001547424e-02  4.541176566678342e-01  3.475741174353303e-12	 2e-2       0.4210    0.83     0.4214           0.83  
+ 1.000000000000000e-01  2.447792823379421e+00  1.580714513882226e-13	 1e-1       2.3477    0.84      2.348           0.86  
+ 1.000000000000000e+00  2.506019141134611e+01  5.730577322954954e-14	 1e+0        24.82    0.91      24.82           0.91