From 6ee0adb808e0d268ccfedbde2456e447c54d70bb Mon Sep 17 00:00:00 2001 From: Ian Jauslin Date: Mon, 30 Nov 2020 13:46:20 -0500 Subject: As presented at Texas A&M on 2020-11-30 --- figs/relative_energy.fig/2020-08-27+2.dat | 5 ++++ figs/relative_energy.fig/Makefile | 43 +++++++++++++++++++++++++++++++ figs/relative_energy.fig/bigeq2.dat | 5 ++++ figs/relative_energy.fig/compare.gnuplot | 31 ++++++++++++++++++++++ figs/relative_energy.fig/simpleq2.dat | 5 ++++ 5 files changed, 89 insertions(+) create mode 100644 figs/relative_energy.fig/2020-08-27+2.dat create mode 100644 figs/relative_energy.fig/Makefile create mode 100644 figs/relative_energy.fig/bigeq2.dat create mode 100644 figs/relative_energy.fig/compare.gnuplot create mode 100644 figs/relative_energy.fig/simpleq2.dat (limited to 'figs/relative_energy.fig') 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 -- cgit v1.2.3-70-g09d2