From 19042107cb35092c674517c8bfe476f1a66c5874 Mon Sep 17 00:00:00 2001 From: Ian Jauslin Date: Wed, 25 Nov 2020 01:10:59 -0500 Subject: Initial commit --- figs/condensate.fig/Makefile | 43 +++++++++++++ figs/condensate.fig/condensate05.gnuplot | 37 +++++++++++ figs/condensate.fig/condensate16.gnuplot | 36 +++++++++++ figs/correlation.fig/Makefile | 43 +++++++++++++ figs/correlation.fig/correlation0001.gnuplot | 28 +++++++++ figs/correlation.fig/correlation02.gnuplot | 28 +++++++++ figs/energy.fig/Makefile | 43 +++++++++++++ figs/energy.fig/energy1.gnuplot | 37 +++++++++++ figs/energy.fig/energy16.gnuplot | 37 +++++++++++ figs/energy_relative_error.fig/Makefile | 71 ++++++++++++++++++++++ figs/energy_relative_error.fig/cmp_energy1.gnuplot | 38 ++++++++++++ .../energy_relative_error.fig/cmp_energy16.gnuplot | 37 +++++++++++ figs/hardcore.fig/Makefile | 47 ++++++++++++++ figs/hardcore.fig/hardcore_condensate.gnuplot | 39 ++++++++++++ figs/hardcore.fig/hardcore_energy.gnuplot | 38 ++++++++++++ figs/tan.fig/Makefile | 35 +++++++++++ figs/tan.fig/tan4.gnuplot | 48 +++++++++++++++ figs/tan.fig/tan7.gnuplot | 48 +++++++++++++++ figs/ux.fig/Makefile | 43 +++++++++++++ figs/ux.fig/ux-0001.gnuplot | 26 ++++++++ figs/ux.fig/ux-02.gnuplot | 26 ++++++++ 21 files changed, 828 insertions(+) create mode 100644 figs/condensate.fig/Makefile create mode 100644 figs/condensate.fig/condensate05.gnuplot create mode 100644 figs/condensate.fig/condensate16.gnuplot create mode 100644 figs/correlation.fig/Makefile create mode 100644 figs/correlation.fig/correlation0001.gnuplot create mode 100644 figs/correlation.fig/correlation02.gnuplot create mode 100644 figs/energy.fig/Makefile create mode 100644 figs/energy.fig/energy1.gnuplot create mode 100644 figs/energy.fig/energy16.gnuplot create mode 100644 figs/energy_relative_error.fig/Makefile create mode 100644 figs/energy_relative_error.fig/cmp_energy1.gnuplot create mode 100644 figs/energy_relative_error.fig/cmp_energy16.gnuplot create mode 100644 figs/hardcore.fig/Makefile create mode 100644 figs/hardcore.fig/hardcore_condensate.gnuplot create mode 100644 figs/hardcore.fig/hardcore_energy.gnuplot create mode 100644 figs/tan.fig/Makefile create mode 100644 figs/tan.fig/tan4.gnuplot create mode 100644 figs/tan.fig/tan7.gnuplot create mode 100644 figs/ux.fig/Makefile create mode 100644 figs/ux.fig/ux-0001.gnuplot create mode 100644 figs/ux.fig/ux-02.gnuplot (limited to 'figs') diff --git a/figs/condensate.fig/Makefile b/figs/condensate.fig/Makefile new file mode 100644 index 0000000..626032e --- /dev/null +++ b/figs/condensate.fig/Makefile @@ -0,0 +1,43 @@ +PROJECTNAME=condensate16 condensate05 + +SIMPLEQ=simpleq + +DATS=simpleq16.dat mueq16.dat bigeq16.dat simpleq0.5.dat mueq0.5.dat bigeq0.5.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, $@) + +simpleq16.dat: + julia $(SIMPLEQ)/main.jl -p "tolerance=1e-11;order=100;maxiter=100;v_a=16" -M simpleq condensate_fraction_rho > $@ +mueq16.dat: + julia $(SIMPLEQ)/main.jl -p "tolerance=1e-11;order=100;maxiter=100;v_a=16" -M mueq condensate_fraction_rho > $@ +bigeq16.dat: + julia -p 8 $(SIMPLEQ)/main.jl -p "eq=bigeq;N=12;P=8;v_a=16" -M medeq condensate_fraction_rho > $@ +simpleq0.5.dat: + julia $(SIMPLEQ)/main.jl -p "tolerance=1e-11;order=100;maxiter=100;v_a=0.5" -M simpleq condensate_fraction_rho > $@ +mueq0.5.dat: + julia $(SIMPLEQ)/main.jl -p "tolerance=1e-11;order=100;maxiter=100;v_a=0.5" -M mueq condensate_fraction_rho > $@ +bigeq0.5.dat: + julia -p 8 $(SIMPLEQ)/main.jl -p "eq=bigeq;N=12;P=8;v_a=0.5" -M medeq condensate_fraction_rho > $@ + + +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/condensate.fig/condensate05.gnuplot b/figs/condensate.fig/condensate05.gnuplot new file mode 100644 index 0000000..66a642e --- /dev/null +++ b/figs/condensate.fig/condensate05.gnuplot @@ -0,0 +1,37 @@ +set ylabel "$\\tilde\\eta$" norotate offset -1,0 +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 ytics 0.6, 0.1 +set mytics 2 + +set xrange [4e-7:100] + +set yrange [:0.06] + +# default output canvas size: 12.5cm x 8.75cm +set term lua tikz size 8,6 standalone + +set key top 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 x + +plot \ + "simpleq0.5.dat" using 1:($3<1e-5 ? $2 : 1/0) with lines ls 1 title "simple",\ + "mueq0.5.dat" using 1:($3<1e-5 ? $2 : 1/0) with lines ls 3 title "med",\ + "bigeq0.5.dat" using 1:($3<1e-5 ? $2 : 1/0) with lines ls 5 title "big",\ + "2020-08-27+0.5.dat" using 1:(1-$3) with points ls 2 title "QMC" + diff --git a/figs/condensate.fig/condensate16.gnuplot b/figs/condensate.fig/condensate16.gnuplot new file mode 100644 index 0000000..b581781 --- /dev/null +++ b/figs/condensate.fig/condensate16.gnuplot @@ -0,0 +1,36 @@ +set ylabel "$\\tilde\\eta$" norotate offset -1,0 +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 ytics 0.6, 0.1 +set mytics 2 + +set xrange [1e-7:100] + + +# default output canvas size: 12.5cm x 8.75cm +set term lua tikz size 8,6 standalone + +set key top 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 x + +plot \ + "simpleq16.dat" using 1:($3<1e-5 ? $2 : 1/0) with lines ls 1 title "simple",\ + "mueq16.dat" using 1:($3<1e-5 ? $2 : 1/0) with lines ls 3 title "med",\ + "bigeq16.dat" using 1:($3<1e-5 ? $2 : 1/0) with lines ls 5 title "big",\ + "2020-08-27+16.dat" using 1:(1-$3) with points ls 2 title "QMC" + diff --git a/figs/correlation.fig/Makefile b/figs/correlation.fig/Makefile new file mode 100644 index 0000000..4a87b58 --- /dev/null +++ b/figs/correlation.fig/Makefile @@ -0,0 +1,43 @@ +PROJECTNAME=correlation0001 correlation02 + +SIMPLEQ=simpleq + +DATS=bigeq02.dat mueq02.dat simpleq02.dat bigeq0001.dat simpleq0001.dat mueq0001.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, $@) + +simpleq0001.dat: + julia $(SIMPLEQ)/main.jl -p "rho=1e-4;minlrho=-6;nlrho=20;eq=simpleq;N=12;P=8;xmax=55" -M medeq 2pt > $@ +mueq0001.dat: + julia $(SIMPLEQ)/main.jl -p "rho=1e-4;minlrho=-6;nlrho=20;eq=mueq;N=12;P=8;xmax=55" -M medeq 2pt > $@ +bigeq0001.dat: + julia $(SIMPLEQ)/main.jl -p "rho=1e-4;minlrho=-6;nlrho=20;v_a=16;eq=bigeq;N=12;P=8;xmax=55" -M medeq 2pt > $@ + +simpleq02.dat: + julia $(SIMPLEQ)/main.jl -p "rho=2e-2;minlrho=-6;nlrho=60;eq=simpleq;N=12;P=8;xmax=10" -M medeq 2pt > $@ +mueq02.dat: + julia $(SIMPLEQ)/main.jl -p "rho=2e-2;minlrho=-6;nlrho=60;eq=mueq;N=12;P=8;xmax=10" -M medeq 2pt > $@ +bigeq02.dat: + julia $(SIMPLEQ)/main.jl -p "rho=2e-2;minlrho=-6;nlrho=60;v_a=16;eq=bigeq;N=12;P=8;xmax=10" -M medeq 2pt > $@ + +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/correlation.fig/correlation0001.gnuplot b/figs/correlation.fig/correlation0001.gnuplot new file mode 100644 index 0000000..24dd900 --- /dev/null +++ b/figs/correlation.fig/correlation0001.gnuplot @@ -0,0 +1,28 @@ +set title "$\\rho=0.0001$" +set ylabel "$\\frac{\\tilde C_2}{\\rho^2}$" norotate offset -1,0 +set xlabel "$|\\mathbf x|$" + +# default output canvas size: 12.5cm x 8.75cm +set term lua tikz size 8,6 standalone + +set key bottom right 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 1 +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 style line 6 linetype rgbcolor "#555500" linewidth 2 + +set pointsize 0.5 + +rho=1e-4 + +plot \ + "simpleq0001.dat" using 1:($2/rho**2) with lines ls 1 title "simple",\ + "mueq0001.dat" using 1:($2/rho**2) with lines ls 3 title "med",\ + "bigeq0001.dat" using 1:($2/rho**2) with lines ls 5 title "big",\ + "2020-10-01+correlation16-0001.dat" using 1:2 ls 2 title "QMC" + diff --git a/figs/correlation.fig/correlation02.gnuplot b/figs/correlation.fig/correlation02.gnuplot new file mode 100644 index 0000000..eefa865 --- /dev/null +++ b/figs/correlation.fig/correlation02.gnuplot @@ -0,0 +1,28 @@ +set title "$\\rho=0.02$" +set ylabel "$\\frac{\\tilde C_2}{\\rho^2}$" norotate offset -1,0 +set xlabel "$|\\mathbf x|$" + +# default output canvas size: 12.5cm x 8.75cm +set term lua tikz size 8,6 standalone + +set key bottom right 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 1 +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 style line 6 linetype rgbcolor "#555500" linewidth 2 + +set pointsize 0.5 + +rho=2e-2 + +plot \ + "simpleq02.dat" using 1:($2/rho**2) with lines ls 1 title "simple",\ + "mueq02.dat" using 1:($2/rho**2) with lines ls 3 title "med",\ + "bigeq02.dat" using 1:($2/rho**2) with lines ls 5 title "big",\ + "2020-08-27+correlation.dat" using 1:2 ls 2 title "QMC" + diff --git a/figs/energy.fig/Makefile b/figs/energy.fig/Makefile new file mode 100644 index 0000000..6ffa316 --- /dev/null +++ b/figs/energy.fig/Makefile @@ -0,0 +1,43 @@ +PROJECTNAME=energy16 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, $@) + +simpleq16.dat: + julia $(SIMPLEQ)/main.jl -p "tolerance=1e-11;order=100;maxiter=100;v_a=16" -M simpleq energy_rho > $@ +mueq16.dat: + julia $(SIMPLEQ)/main.jl -p "tolerance=1e-11;order=100;maxiter=100;v_a=16" -M mueq energy_rho > $@ +bigeq16.dat: + julia -p 8 $(SIMPLEQ)/main.jl -p "eq=bigeq;N=12;P=8;v_a=16" -M medeq energy_rho > $@ +simpleq1.dat: + julia $(SIMPLEQ)/main.jl -p "tolerance=1e-11;order=100;maxiter=100;v_a=1" -M simpleq energy_rho > $@ +mueq1.dat: + julia $(SIMPLEQ)/main.jl -p "tolerance=1e-11;order=100;maxiter=100;v_a=1" -M mueq energy_rho > $@ +bigeq1.dat: + julia -p 8 $(SIMPLEQ)/main.jl -p "eq=bigeq;N=12;P=8;v_a=1" -M medeq energy_rho > $@ + + +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/energy.fig/energy1.gnuplot b/figs/energy.fig/energy1.gnuplot new file mode 100644 index 0000000..9667dba --- /dev/null +++ b/figs/energy.fig/energy1.gnuplot @@ -0,0 +1,37 @@ +set title "$v(\\mathbf x)=e^{-|\\mathbf x|}$" +set ylabel "$\\displaystyle\\frac{\\tilde e}{\\rho}$" norotate offset -1,0 +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 8, 1, 13 +set mytics 2 +set yrange [7.5:13] + + +# default output canvas size: 12.5cm x 8.75cm +set term lua tikz size 8,6 standalone + +set key bottom right 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 x + +plot \ + "simpleq1.dat" using 1:($3<1e-5 ? ($2/$1) : 1/0) with lines ls 1 title "simple" ,\ + "mueq1.dat" using 1:($3<1e-5 ? ($2/$1) : 1/0) with lines ls 3 title "med",\ + "bigeq1.dat" using 1:($3<1e-5 ? ($2/$1) : 1/0) with lines ls 5 title "big",\ + "2020-10-15+energy.dat" using 1:($2/$1) with points ls 2 title "QMC" + diff --git a/figs/energy.fig/energy16.gnuplot b/figs/energy.fig/energy16.gnuplot new file mode 100644 index 0000000..4fd8de6 --- /dev/null +++ b/figs/energy.fig/energy16.gnuplot @@ -0,0 +1,37 @@ +set title "$v(\\mathbf x)=16e^{-|\\mathbf x|}$" +set ylabel "$\\displaystyle\\frac{\\tilde e}{\\rho}$" norotate offset -1,0 +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 20, 40, 220 +set mytics 2 +set yrange [20:220] + + +# default output canvas size: 12.5cm x 8.75cm +set term lua tikz size 8,6 standalone + +set key bottom right 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 x + +plot \ + "simpleq16.dat" using 1:($3<1e-5 ? ($2/$1) : 1/0) with lines ls 1 title "simple",\ + "mueq16.dat" using 1:($3<1e-5 ? ($2/$1) : 1/0) with lines ls 3 title "med",\ + "bigeq16.dat" using 1:($3<1e-5 ? ($2/$1) : 1/0) with lines ls 5 title "big",\ + "2020-08-27+16.dat" using 1:($2/$1) with points ls 2 title "QMC" + 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" diff --git a/figs/hardcore.fig/Makefile b/figs/hardcore.fig/Makefile new file mode 100644 index 0000000..1cf8737 --- /dev/null +++ b/figs/hardcore.fig/Makefile @@ -0,0 +1,47 @@ +PROJECTNAME=hardcore_energy hardcore_condensate +SIMPLEQ=simpleq + +DATS=simpleq-energy.dat simpleq_soft-energy.dat mueq_soft-energy.dat bigeq_soft-energy.dat simpleq-condensate.dat simpleq_soft-condensate.dat mueq_soft-condensate.dat bigeq_soft-condensate.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, $@) + +simpleq-energy.dat: + julia -p 8 $(SIMPLEQ)/main.jl -p "N=12;P=8;J=5;nlrho=24;maxlrho=0.15" -M simpleq-hardcore energy_rho > $@ +simpleq_soft-energy.dat: + julia $(SIMPLEQ)/main.jl -p "maxiter=100;maxlrho=0.15;v_a=512;v_b=1.2494804800401416" -U tent -M simpleq energy_rho > $@ +mueq_soft-energy.dat: + julia $(SIMPLEQ)/main.jl -p "maxiter=100;maxlrho=0.15;v_a=512;v_b=1.2494804800401416" -U tent -M mueq energy_rho > $@ +bigeq_soft-energy.dat: + julia -p 8 $(SIMPLEQ)/main.jl -p "N=12;P=8;eq=bigeq;maxiter=100;maxlrho=0.15;nlrho=64;v_a=512;v_b=1.2494804800401416" -U tent -M medeq energy_rho_init_prevrho > $@ + +simpleq-condensate.dat: + julia -p 8 $(SIMPLEQ)/main.jl -p "N=12;P=8;J=5;maxlrho=0.15;nlrho=24" -M simpleq-hardcore condensate_fraction_rho > $@ +simpleq_soft-condensate.dat: + julia $(SIMPLEQ)/main.jl -p "maxiter=100;maxlrho=0.15;v_a=512;v_b=1.2494804800401416" -U tent -M simpleq condensate_fraction_rho > $@ +mueq_soft-condensate.dat: + julia $(SIMPLEQ)/main.jl -p "maxiter=100;maxlrho=0.15;v_a=512;v_b=1.2494804800401416" -U tent -M mueq condensate_fraction_rho > $@ +bigeq_soft-condensate.dat: + julia -p 8 $(SIMPLEQ)/main.jl -p "N=12;P=8;eq=bigeq;maxiter=100;maxlrho=0.15;nlrho=64;v_a=512;v_b=1.2494804800401416" -U tent -M medeq condensate_fraction_rho > $@ + +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/hardcore.fig/hardcore_condensate.gnuplot b/figs/hardcore.fig/hardcore_condensate.gnuplot new file mode 100644 index 0000000..974c9fd --- /dev/null +++ b/figs/hardcore.fig/hardcore_condensate.gnuplot @@ -0,0 +1,39 @@ +set ylabel "$\\tilde\\eta$" norotate +set xlabel "$\\rho$" + +set xtics 1e-6, 100, 1 +set xtics add ("$10^{-6}$" 0.000001, "$10^{-4}$" 0.0001, "$10^{-2}$" 0.01, "$1$" 1.0) +unset mxtics + +set xrange [:sqrt(2)] + +#set ytics 1e-6, 10, 1 +#set ytics add ("$10^{-6}$" 0.000001, "$10^{-5}$" 1e-5, "$10^{-4}$" 0.0001, "$10^{-3}$" 1e-3, "$10^{-2}$" 0.01, "$10^{-1}$" 1e-1, "$1$" 1.0) +#unset mytics +set yrange [0:1] + +# default output canvas size: 12.5cm x 8.75cm +set term lua tikz size 8,6 standalone + +set key top 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 style line 6 linetype rgbcolor "#999999" linewidth 2 + +set pointsize 1 + +set logscale x + +plot \ + "simpleq-condensate.dat" using 1:($3<1e-2 ? $2 : 1/0) ls 6 title "simple hc",\ + "simpleq_soft-condensate.dat" using 1:($3<1e-2 ? $2 : 1/0) with lines ls 1 title "simple" ,\ + "mueq_soft-condensate.dat" using 1:($3<1e-2 ? $2 : 1/0) with lines ls 3 title "med",\ + "bigeq_soft-condensate.dat" using 1:($3<1e-2 ? $2 : 1/0) with lines ls 5 title "big",\ + "GBC99-condensate.dat" using 1:(1-$2) ls 2 title "QMC [15]" + diff --git a/figs/hardcore.fig/hardcore_energy.gnuplot b/figs/hardcore.fig/hardcore_energy.gnuplot new file mode 100644 index 0000000..9967232 --- /dev/null +++ b/figs/hardcore.fig/hardcore_energy.gnuplot @@ -0,0 +1,38 @@ +set ylabel "$\\frac{\\tilde e}\\rho$" norotate +set xlabel "$\\rho$" + +set xtics 1e-6, 100, 1 +set xtics add ("$10^{-6}$" 0.000001, "$10^{-4}$" 0.0001, "$10^{-2}$" 0.01, "$1$" 1.0) +unset mxtics + +set xrange [:sqrt(2)] + +#set ytics 1e-6, 10, 1 +#set ytics add ("$10^{-6}$" 0.000001, "$10^{-5}$" 1e-5, "$10^{-4}$" 0.0001, "$10^{-3}$" 1e-3, "$10^{-2}$" 0.01, "$10^{-1}$" 1e-1, "$1$" 1.0) +#unset mytics + +# default output canvas size: 12.5cm x 8.75cm +set term lua tikz size 8,6 standalone + +set key top 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 style line 6 linetype rgbcolor "#999999" linewidth 2 + +set pointsize 1 + +set logscale xy + +plot \ + "simpleq-energy.dat" using 1:($3<1e-2 ? $2/$1 : 1/0) ls 6 title "simple hc",\ + "simpleq_soft-energy.dat" using 1:($3<1e-2 ? $2/$1 : 1/0) with lines ls 1 title "simple",\ + "mueq_soft-energy.dat" using 1:($3<1e-2 ? $2/$1 : 1/0) every ::::95 with lines ls 3 title "med",\ + "bigeq_soft-energy.dat" using 1:($3<1e-2 ? $2/$1 : 1/0) with lines ls 5 title "big",\ + "GBC99-energy.dat" using 1:($2/(2*$1)) ls 2 title "QMC [15]" + diff --git a/figs/tan.fig/Makefile b/figs/tan.fig/Makefile new file mode 100644 index 0000000..02ea899 --- /dev/null +++ b/figs/tan.fig/Makefile @@ -0,0 +1,35 @@ +PROJECTNAME=tan7 tan4 + +SIMPLEQ=simpleq + +DATS=bigeq7.dat bigeq4.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, $@) + +bigeq7.dat: + julia $(SIMPLEQ)/main.jl -p "eq=bigeq;N=12;P=8;J=10;maxiter=100;v_a=1;rho=1e-7;minlrho=-7;nlrho=1" -M medeq momentum_distribution > $@ +bigeq4.dat: + julia $(SIMPLEQ)/main.jl -p "eq=bigeq;N=12;P=8;J=10;maxiter=100;v_a=1;rho=1e-4;minlrho=-4;nlrho=1" -M medeq momentum_distribution > $@ + +install: $(PDFS) + cp $^ $(INSTALLDIR)/ + +clean-aux: + rm -f $(addsuffix .tikz.tex, $(PROJECTNAME)) + rm -f $(addsuffix .aux, $(PROJECTNAME)) + rm -f $(addsuffix .log, $(PROJECTNAME)) + rm -f fit.log + +clean-dat: + rm -f $(DATS) + +clean-tex: + rm -f $(PDFS) + +clean: clean-aux clean-tex diff --git a/figs/tan.fig/tan4.gnuplot b/figs/tan.fig/tan4.gnuplot new file mode 100644 index 0000000..f1424b6 --- /dev/null +++ b/figs/tan.fig/tan4.gnuplot @@ -0,0 +1,48 @@ +set title "$\\rho=10^{-4}$" +set ylabel norotate "$\\tilde{\\mathcal M}$" +set xlabel "$|\\mathbf k|$" + +set xtics 1e-3, 100, 1e3 +set xtics add ("$10^{-3}$" 0.001, "$10^{-1}$" 1e-1, "$10^{1}$" 10, "$10^{3}$" 1000) +unset mxtics +set xrange [1e-3:2e3] + +set ytics 1e-40, 1e10, 1 +set ytics add ("$10^{-40}$" 1e-40, "$10^{-30}$" 1e-30, "$10^{-20}$" 1e-20, "$10^{-10}$" 1e-10) +unset mytics +set yrange [1e-41:1e10] + +# 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 3 +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 3 +set style line 5 linetype rgbcolor "#DAA520" linewidth 3 +set style line 6 linetype rgbcolor "#660000" linewidth 3 +set style line 7 linetype rgbcolor "#003333" linewidth 3 + +set pointsize 1 + +set logscale xy + +f(x)=a-4*x +g(x)=c-12*x + +fit [-2:-1] f(x) "bigeq4.dat" u (log10($1)):(log10(abs($2))) via a +fit [1:] g(x) "bigeq4.dat" u (log10($1)):(log10(abs($2))) via c + +set label "$|\\mathbf k|^{-4}$" at 1e-2,10**(a) textcolor "#660000" +set label "$|\\mathbf k|^{-12}$" at 1e2,10**(c-22) textcolor "#003333" + +plot \ + "bigeq4.dat" using 1:(abs($2)) every::2::118 with lines ls 5 ,\ + 10**(f(log10(x))) ls 6 dashtype "." ,\ + 10**(g(log10(x))) ls 7 dashtype "." ,\ + + diff --git a/figs/tan.fig/tan7.gnuplot b/figs/tan.fig/tan7.gnuplot new file mode 100644 index 0000000..6ccda8b --- /dev/null +++ b/figs/tan.fig/tan7.gnuplot @@ -0,0 +1,48 @@ +set title "$\\rho=10^{-7}$" +set ylabel norotate "$\\tilde{\\mathcal M}$" +set xlabel "$|\\mathbf k|$" + +set xtics 1e-3, 1000, 1e6 +set xtics add ("$10^{-3}$" 0.001, "$10^{0}$" 1, "$10^{3}$" 1000, "$10^6$" 1000000) +unset mxtics +set xrange [1e-3:1000000] + +set ytics 1e-60, 1e20, 1 +set ytics add ("$10^{-60}$" 1e-60, "$10^{-40}$" 1e-40, "$10^{-20}$" 1e-20) +unset mytics +set yrange [1e-61:1e10] + +# 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 3 +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 3 +set style line 5 linetype rgbcolor "#DAA520" linewidth 3 +set style line 6 linetype rgbcolor "#660000" linewidth 3 +set style line 7 linetype rgbcolor "#003333" linewidth 3 + +set pointsize 1 + +set logscale xy + +f(x)=a-4*x +g(x)=c-12*x + +fit [-2:-1] f(x) "bigeq7.dat" u (log10($1)):(log10(abs($2))) via a +fit [1:] g(x) "bigeq7.dat" u (log10($1)):(log10(abs($2))) via c + +set label "$|\\mathbf k|^{-4}$" at 1e-2,10**(a-4) textcolor "#660000" +set label "$|\\mathbf k|^{-12}$" at 1e3,10**(c-34) textcolor "#003333" + +plot \ + "bigeq7.dat" using 1:(abs($2)) with lines ls 5 ,\ + 10**(f(log10(x))) ls 6 dashtype "." ,\ + 10**(g(log10(x))) ls 7 dashtype "." ,\ + + diff --git a/figs/ux.fig/Makefile b/figs/ux.fig/Makefile new file mode 100644 index 0000000..bd20f97 --- /dev/null +++ b/figs/ux.fig/Makefile @@ -0,0 +1,43 @@ +PROJECTNAME=ux-02 ux-0001 + +SIMPLEQ=simpleq + +DATS=simpleq16-0.02.dat mueq16-0.02.dat bigeq16-0.02.dat simpleq16-0.0001.dat mueq16-0.0001.dat bigeq16-0.0001.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, $@) + +simpleq16-0.02.dat: + julia $(SIMPLEQ)/main.jl -p "v_a=16;xmax=10;nx=200;rho=2e-2" -M simpleq ux >$@ +mueq16-0.02.dat: + julia $(SIMPLEQ)/main.jl -p "v_a=16;xmax=10;nx=200;rho=2e-2" -M mueq ux > $@ +bigeq16-0.02.dat: + julia $(SIMPLEQ)/main.jl -p "rho=2e-2;minlrho=-6;nlrho=60;v_a=16;eq=bigeq;N=16;P=8;xmax=10;nx=200;N=16;P=8" -M medeq ux > $@ + +simpleq16-0.0001.dat: + julia $(SIMPLEQ)/main.jl -p "v_a=16;xmax=40;nx=200;rho=1e-4" -M simpleq ux >$@ +mueq16-0.0001.dat: + julia $(SIMPLEQ)/main.jl -p "v_a=16;xmax=40;nx=200;rho=1e-4" -M mueq ux > $@ +bigeq16-0.0001.dat: + julia $(SIMPLEQ)/main.jl -p "rho=1e-4;minlrho=-6;nlrho=20;v_a=16;eq=bigeq;N=16;P=8;xmax=40;nx=200;N=16;P=8" -M medeq ux > $@ + +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/ux.fig/ux-0001.gnuplot b/figs/ux.fig/ux-0001.gnuplot new file mode 100644 index 0000000..4cc6fd7 --- /dev/null +++ b/figs/ux.fig/ux-0001.gnuplot @@ -0,0 +1,26 @@ +set title "$\\rho=0.0001$" +set ylabel "$\\tilde g_2$" norotate offset -1,0 +set xlabel "$|\\mathbf x|$" + +# default output canvas size: 12.5cm x 8.75cm +set term lua tikz size 8,6 standalone + +set key bottom right 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 style line 6 linetype rgbcolor "#555500" linewidth 2 + +set pointsize 1 + +plot \ + "simpleq16-0.0001.dat" using 1:(1-$2) with lines ls 1 title "simple",\ + "mueq16-0.0001.dat" using 1:(1-$2) with lines ls 3 title "med",\ + "bigeq16-0.0001.dat" using 1:(1-$2) with lines ls 5 title "big",\ + "2020-10-01+u16-0001.dat" using 1:2 ls 2 title "QMC" + diff --git a/figs/ux.fig/ux-02.gnuplot b/figs/ux.fig/ux-02.gnuplot new file mode 100644 index 0000000..6379449 --- /dev/null +++ b/figs/ux.fig/ux-02.gnuplot @@ -0,0 +1,26 @@ +set title "$\\rho=0.02$" +set ylabel "$\\tilde g_2$" norotate offset -1,0 +set xlabel "$|\\mathbf x|$" + +# default output canvas size: 12.5cm x 8.75cm +set term lua tikz size 8,6 standalone + +set key bottom right 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 style line 6 linetype rgbcolor "#555500" linewidth 2 + +set pointsize 0.5 + +plot \ + "simpleq16-0.02.dat" using 1:(1-$2) with lines ls 1 title "simple",\ + "mueq16-0.02.dat" using 1:(1-$2) with lines ls 3 title "med",\ + "bigeq16-0.02.dat" using 1:(1-$2) with lines ls 5 title "big",\ + "2020-08-27+u.dat" using 1:2 ls 2 title "QMC" + -- cgit v1.2.3-54-g00ecf