From d4b3b0cc4c5b53c5ba4f7d0dd5070d7457f964b8 Mon Sep 17 00:00:00 2001 From: Ian Jauslin Date: Sun, 26 Feb 2023 19:10:42 -0500 Subject: Initial commit --- figs/condensate.fig/Makefile | 43 ++++++++++++ figs/condensate.fig/bigeq.dat | 100 +++++++++++++++++++++++++++ figs/condensate.fig/condensate.gnuplot | 40 +++++++++++ figs/condensate.fig/condensate_medeq.gnuplot | 43 ++++++++++++ figs/condensate.fig/medeq.dat | 100 +++++++++++++++++++++++++++ 5 files changed, 326 insertions(+) create mode 100644 figs/condensate.fig/Makefile create mode 100644 figs/condensate.fig/bigeq.dat create mode 100644 figs/condensate.fig/condensate.gnuplot create mode 100644 figs/condensate.fig/condensate_medeq.gnuplot create mode 100644 figs/condensate.fig/medeq.dat (limited to 'figs/condensate.fig') diff --git a/figs/condensate.fig/Makefile b/figs/condensate.fig/Makefile new file mode 100644 index 0000000..cac5f40 --- /dev/null +++ b/figs/condensate.fig/Makefile @@ -0,0 +1,43 @@ +PROJECT=$(basename $(wildcard *.gnuplot)) + +SIMPLESOLV=simplesolv + +all: $(addsuffix .pdf, $(PROJECT)) + +define gnuplot2pdf + gnuplot $(1).gnuplot > $(1).tikz.tex + pdflatex -jobname $(1) -file-line-error $(1).tikz.tex +endef + +define medeq_run + julia $(SIMPLESOLV)/main.jl -p "eq=medeq;order=200;maxiter=100;v_a=8;minlrho=-6;maxlrho=2;nlrho=100" -M easyeq condensate_fraction_rho > medeq.dat +endef +define bigeq_run + julia -p 16 $(SIMPLESOLV)/main.jl -p "eq=bigeq;N=12;P=8;J=10;maxiter=100;v_a=8;minlrho=-6;nlrho=100;maxlrho=2" -M anyeq condensate_fraction_rho > bigeq.dat +endef + +medeq.dat: + $(call medeq_run) + +bigeq.dat: + $(call bigeq_run) + +DATS=medeq.dat bigeq.dat + +condensate.pdf: $(DATS) + $(call gnuplot2pdf,$(basename $@)) +condensate_medeq.pdf: $(DATS) + $(call gnuplot2pdf,$(basename $@)) + + +clean-aux: + rm -f $(addsuffix .tikz.tex, $(PROJECT)) + rm -f *.aux *.log + +clean-dat: + rm -f $(DATS) + +clean-pdf: + rm -f $(addsuffix .pdf, $(PROJECT)) + +clean: clean-aux clean-pdf diff --git a/figs/condensate.fig/bigeq.dat b/figs/condensate.fig/bigeq.dat new file mode 100644 index 0000000..7beede0 --- /dev/null +++ b/figs/condensate.fig/bigeq.dat @@ -0,0 +1,100 @@ + 1.000000000000000e-06 8.524653565633833e-03 7.195282669866215e-15 + 1.204503540258781e-06 9.482931278294420e-03 1.219933599614962e-14 + 1.450828778495940e-06 1.052073432648738e-02 1.476802021057201e-14 + 1.747528400007683e-06 1.164867071649953e-02 2.986683032499419e-15 + 2.104904144512022e-06 1.288059904788134e-02 1.880210643928359e-14 + 2.535364493970111e-06 1.423165106499452e-02 1.524596171157415e-14 + 3.053855508833412e-06 1.571746025308502e-02 1.021384261300024e-13 + 3.678379771828634e-06 1.735409651757319e-02 8.142001190132888e-15 + 4.430621457583877e-06 1.915832851889812e-02 1.213702793846234e-14 + 5.336699231206313e-06 2.114801579568841e-02 1.737691905063488e-14 + 6.428073117284319e-06 2.334252989270571e-02 8.865535490910355e-12 + 7.742636826811277e-06 2.576315226604090e-02 2.686577129397137e-12 + 9.326033468832200e-06 2.843343064254577e-02 1.438084696630850e-12 + 1.123324032978027e-05 3.137950314614005e-02 2.937439611729668e-12 + 1.353047774579808e-05 3.463041485755209e-02 5.399796836950549e-12 + 1.629750834620643e-05 3.821844838166732e-02 7.534307445347654e-12 + 1.963040650040273e-05 4.217947063976608e-02 8.540778248911911e-12 + 2.364489412645407e-05 4.655327365355863e-02 8.130847144221019e-12 + 2.848035868435805e-05 5.138387252824026e-02 6.515425217203499e-12 + 3.430469286314919e-05 5.671972969979729e-02 4.911369751416716e-12 + 4.132012400115334e-05 6.261389833596930e-02 8.912668688516455e-15 + 4.977023564332114e-05 6.912410160811572e-02 3.766883359036800e-14 + 5.994842503189409e-05 7.631276517641099e-02 6.820252575321234e-15 + 7.220809018385471e-05 8.424698806468546e-02 2.022663500277492e-14 + 8.697490026177834e-05 9.299838732648825e-02 2.090481005707888e-14 + 1.047615752789665e-04 1.026427162611128e-01 1.226684774874274e-14 + 1.261856883066021e-04 1.132591594617568e-01 1.006143610935187e-14 + 1.519911082952935e-04 1.249292448150214e-01 1.020576205644618e-14 + 1.830738280295370e-04 1.377353406641975e-01 1.570903430629419e-14 + 2.205130739903046e-04 1.517586789773243e-01 2.619798908910617e-14 + 2.656087782946686e-04 1.670767645022099e-01 7.901048564229978e-15 + 3.199267137797384e-04 1.837599621474728e-01 3.746068672715705e-15 + 3.853528593710532e-04 2.018670716203985e-01 3.001662303433755e-15 + 4.641588833612782e-04 2.214397964264479e-01 1.103630516809694e-14 + 5.590810182512228e-04 2.424961195590747e-01 1.366652070574403e-14 + 6.734150657750829e-04 2.650226615959725e-01 4.701570959947225e-14 + 8.111308307896872e-04 2.889661589017276e-01 1.102455991777986e-13 + 9.770099572992256e-04 3.142243325472197e-01 2.194837449083438e-13 + 1.176811952434999e-03 3.406366355397510e-01 4.291634987063047e-13 + 1.417474162926806e-03 3.679756223891710e-01 7.877925955154667e-13 + 1.707352647470692e-03 3.959399384752566e-01 1.337225989099042e-12 + 2.056512308348653e-03 4.241501502647564e-01 2.087839965182089e-12 + 2.477076355991711e-03 4.521487893796861e-01 3.006779535875992e-12 + 2.983647240283340e-03 4.794059963287035e-01 3.981564512935711e-12 + 3.593813663804629e-03 5.053319487611297e-01 4.827221170083811e-12 + 4.328761281083062e-03 5.292967869119043e-01 5.317794467285372e-12 + 5.214008287999690e-03 5.506579873265355e-01 5.310929880704551e-12 + 6.280291441834260e-03 5.687941257240574e-01 4.789376962806014e-12 + 7.564633275546291e-03 5.831428216724064e-01 3.897122325557593e-12 + 9.111627561154896e-03 5.932395516139335e-01 2.854055547985561e-12 + 1.097498765493057e-02 5.987531905030805e-01 1.868866977359501e-12 + 1.321941148466031e-02 5.995138571635353e-01 1.082542791147572e-12 + 1.592282793341094e-02 5.955291021078862e-01 5.656821139464731e-13 + 1.917910261672489e-02 5.869857426377714e-01 2.660387742778624e-13 + 2.310129700083163e-02 5.742365544932845e-01 1.116641921715526e-13 + 2.782559402207126e-02 5.577732070731861e-01 4.265686091722418e-14 + 3.351602650938848e-02 5.381888149470879e-01 1.499258345731993e-14 + 4.037017258596558e-02 5.161348489423120e-01 1.586748589205255e-14 + 4.862601580065353e-02 4.922776396647696e-01 1.074536787081187e-14 + 5.857020818056673e-02 4.672592573491898e-01 4.516275882902432e-15 + 7.054802310718646e-02 4.416663252899820e-01 4.462190018920522e-15 + 8.497534359086456e-02 4.160086630420100e-01 2.795856670789435e-15 + 1.023531021899027e-01 3.907079735362658e-01 3.070112609272139e-15 + 1.232846739442068e-01 3.660954314530985e-01 3.711572850823813e-15 + 1.484968262254467e-01 3.424161914304082e-01 5.624460331260991e-15 + 1.788649529057435e-01 3.198385361615758e-01 3.949985464132756e-15 + 2.154434690031887e-01 2.984655183858066e-01 5.239547190388315e-15 + 2.595024211399737e-01 2.783473502776277e-01 3.356133932711423e-15 + 3.125715849688241e-01 2.594932958267652e-01 1.981023711310264e-15 + 3.764935806792472e-01 2.418823046553480e-01 2.359177607668277e-15 + 4.534878508128591e-01 2.254720208381897e-01 1.402803994597511e-15 + 5.462277217684348e-01 2.102060850751125e-01 2.553306469267983e-15 + 6.579332246575682e-01 1.960198295921745e-01 1.235408236856520e-15 + 7.924828983539186e-01 1.828445625979321e-01 2.839059014773400e-15 + 9.545484566618347e-01 1.706106766192757e-01 1.234537422134605e-15 + 1.149756995397738e+00 1.592498142496109e-01 2.193439580869223e-15 + 1.384886371393875e+00 1.486963029452997e-01 1.648286202055480e-15 + 1.668100537200059e+00 1.388880394362550e-01 3.124563475483895e-15 + 2.009233002565050e+00 1.297669712794426e-01 2.212240236489514e-15 + 2.420128264794383e+00 1.212792919665664e-01 1.822578297170560e-15 + 2.915053062825182e+00 1.133754387934347e-01 1.293726888514874e-15 + 3.511191734215135e+00 1.060099602551520e-01 1.638195775953276e-15 + 4.229242874389508e+00 9.914130199191040e-02 1.752453175614045e-15 + 5.094138014816386e+00 9.273154661320844e-02 2.383213588406398e-15 + 6.135907273413176e+00 8.674613226064723e-02 4.673666845465575e-15 + 7.390722033525790e+00 8.115356686357338e-02 3.843207367033268e-15 + 8.902150854450392e+00 7.592514919615562e-02 2.889463903956984e-15 + 1.072267222010325e+01 7.103470363590553e-02 2.577031650678448e-15 + 1.291549665014885e+01 6.645833257136877e-02 3.567371947044331e-15 + 1.555676143930475e+01 6.217418839914374e-02 2.678037011024236e-15 + 1.873817422860387e+01 5.816226574695751e-02 3.465354011432907e-15 + 2.257019719633922e+01 5.440421374532894e-02 2.975998245290574e-15 + 2.718588242732946e+01 5.088316761973691e-02 3.428201805340647e-15 + 3.274549162877732e+01 4.758359840543907e-02 3.879983823861438e-15 + 3.944206059437664e+01 4.449117916775601e-02 3.988638949817691e-15 + 4.750810162102803e+01 4.159266589368779e-02 4.204239894057915e-15 + 5.722367659350220e+01 3.887579136661706e-02 4.017134010366431e-15 + 6.892612104349709e+01 3.632917078016813e-02 4.381717563239069e-15 + 8.302175681319753e+01 3.394221826626880e-02 9.103579900515236e-12 + 1.000000000000000e+02 3.170507358779736e-02 6.838220457800429e-12 diff --git a/figs/condensate.fig/condensate.gnuplot b/figs/condensate.fig/condensate.gnuplot new file mode 100644 index 0000000..51b2e98 --- /dev/null +++ b/figs/condensate.fig/condensate.gnuplot @@ -0,0 +1,40 @@ +set ylabel norotate "$\\eta$" +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, 1e-5, 2e-5 +#set ytics add ("$10^{-5}$" 1e-5) +#set mytics +set yrange [0:1] + +# default output canvas size: 12.5cm x 8.75cm +set term lua tikz size 8,6 standalone + +set key off + + +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 +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 dashtype (2,2) +set style line 11 linetype rgbcolor "#4169E1" linewidth 1 dashtype (1,2) +set style line 12 linetype rgbcolor "#DC143C" linewidth 1 dashtype (1,2) + +a=3.233946179033422 +bog(rho)=1-8./(3*sqrt(pi))*sqrt(rho*a**3) + +set logscale x + +set arrow from 0.9e-3,1 to 0.9e-3,0 nohead ls 12 +set arrow from 0.2,1 to 0.2,0 nohead ls 12 + +plot \ + "bigeq.dat" using 1:(1-$2) with lines ls 2 + + #bog(x) with lines ls 6 title "Bog" ,\ + diff --git a/figs/condensate.fig/condensate_medeq.gnuplot b/figs/condensate.fig/condensate_medeq.gnuplot new file mode 100644 index 0000000..71d84fa --- /dev/null +++ b/figs/condensate.fig/condensate_medeq.gnuplot @@ -0,0 +1,43 @@ +set ylabel norotate "$\\eta$" +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, 1e-5, 2e-5 +#set ytics add ("$10^{-5}$" 1e-5) +#set mytics +set yrange [0:1] + +# 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 2 height 0.3 spacing 1.3 + + +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 +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 dashtype (2,2) +set style line 11 linetype rgbcolor "#4169E1" linewidth 1 dashtype (1,2) +set style line 12 linetype rgbcolor "#DC143C" linewidth 1 dashtype (1,2) + +a=3.233946179033422 +bog(rho)=1-8./(3*sqrt(pi))*sqrt(rho*a**3) + +set logscale x + +set arrow from 1.9e-3,1 to 1.9e-3,0 nohead ls 11 +set arrow from 0.05,1 to 0.05,0 nohead ls 11 +set arrow from 0.9e-3,1 to 0.9e-3,0 nohead ls 12 +set arrow from 0.2,1 to 0.2,0 nohead ls 12 + +plot \ + "medeq.dat" using 1:(1-$2) with lines ls 1 title "med" ,\ + "bigeq.dat" using 1:(1-$2) with lines ls 2 title "big" + + #bog(x) with lines ls 6 title "Bog" ,\ + diff --git a/figs/condensate.fig/medeq.dat b/figs/condensate.fig/medeq.dat new file mode 100644 index 0000000..bd24c39 --- /dev/null +++ b/figs/condensate.fig/medeq.dat @@ -0,0 +1,100 @@ + 1.000000000000000e-06 8.839225668273828e-03 1.112362506292058e-13 + 1.204503540258781e-06 9.710202722409071e-03 2.799081721322997e-14 + 1.450828778495940e-06 1.066787984589104e-02 9.646668617782145e-14 + 1.747528400007683e-06 1.172104457893248e-02 2.837740288150831e-14 + 2.104904144512022e-06 1.287940010882542e-02 1.940422291133097e-14 + 2.535364493970111e-06 1.415366354830729e-02 2.160144539549974e-14 + 3.053855508833412e-06 1.555567470706874e-02 6.913033576413322e-14 + 3.678379771828634e-06 1.709851622208548e-02 3.529301449695497e-14 + 4.430621457583877e-06 1.879664595680302e-02 3.708718434666157e-14 + 5.336699231206313e-06 2.066604245630733e-02 3.257075643408975e-14 + 6.428073117284319e-06 2.272436410952421e-02 7.614257224517049e-14 + 7.742636826811277e-06 2.499112240989452e-02 1.345107953771322e-13 + 9.326033468832200e-06 2.748786930284112e-02 1.564109988333924e-13 + 1.123324032978027e-05 3.023839801257560e-02 3.524401616970253e-13 + 1.353047774579808e-05 3.326895588691060e-02 5.915310424583823e-13 + 1.629750834620643e-05 3.660846660285050e-02 9.654172456076661e-13 + 1.963040650040273e-05 4.028875744290850e-02 1.567378423156396e-12 + 2.364489412645407e-05 4.434478515682510e-02 2.565768885039563e-12 + 2.848035868435805e-05 4.881485103206280e-02 4.198938828273798e-12 + 3.430469286314919e-05 5.374079204785320e-02 6.693326701361960e-12 + 4.132012400115334e-05 5.916813022465738e-02 1.282891351842603e-14 + 4.977023564332114e-05 6.514615634788637e-02 2.685971612008854e-14 + 5.994842503189409e-05 7.172791703687389e-02 1.510292982835255e-14 + 7.220809018385471e-05 7.897006562610288e-02 8.225023204630623e-15 + 8.697490026177834e-05 8.693252767306676e-02 1.181995900107375e-14 + 1.047615752789665e-04 9.567792151389125e-02 1.579760391473529e-14 + 1.261856883066021e-04 1.052706639611355e-01 1.323647673296498e-14 + 1.519911082952935e-04 1.157756823392874e-01 1.399160931186238e-14 + 1.830738280295370e-04 1.272566486606280e-01 1.455628208441605e-14 + 2.205130739903046e-04 1.397736527937312e-01 2.276953826701700e-14 + 2.656087782946686e-04 1.533802427823905e-01 2.927894244062454e-14 + 3.199267137797384e-04 1.681197865588261e-01 3.798057512829229e-14 + 3.853528593710532e-04 1.840211548677170e-01 1.376762463408949e-13 + 4.641588833612782e-04 2.010937941680861e-01 4.971157121525018e-14 + 5.590810182512228e-04 2.193223521156883e-01 6.022869581305702e-14 + 6.734150657750829e-04 2.386611340274900e-01 7.000429865660912e-14 + 8.111308307896872e-04 2.590287969632276e-01 1.761575689941084e-13 + 9.770099572992256e-04 2.803038110155694e-01 9.195321217574250e-14 + 1.176811952434999e-03 3.023213092811143e-01 9.942888124370229e-14 + 1.417474162926806e-03 3.248719778026897e-01 9.771829753425081e-14 + 1.707352647470692e-03 3.477035747274209e-01 9.676129230098743e-14 + 2.056512308348653e-03 3.705254946950684e-01 9.808951587632371e-14 + 2.477076355991711e-03 3.930165111923665e-01 9.469446613818251e-14 + 2.983647240283340e-03 4.148354647231526e-01 9.135639239279189e-14 + 3.593813663804629e-03 4.356342735753405e-01 9.452103047087137e-14 + 4.328761281083062e-03 4.550722993419716e-01 6.340576466252691e-14 + 5.214008287999690e-03 4.728308731866177e-01 5.381048929967230e-14 + 6.280291441834260e-03 4.886267320820731e-01 7.810882165096904e-14 + 7.564633275546291e-03 5.022232412549393e-01 4.156292768592150e-14 + 9.111627561154896e-03 5.134385631744048e-01 4.964091113337942e-14 + 1.097498765493057e-02 5.221503153414601e-01 1.136043314532288e-13 + 1.321941148466031e-02 5.282966646005823e-01 1.305427879388923e-13 + 1.592282793341094e-02 5.318741654433237e-01 1.749382689682514e-13 + 1.917910261672489e-02 5.329329149385416e-01 6.101521560074490e-14 + 2.310129700083163e-02 5.315697459184463e-01 5.254468894478471e-14 + 2.782559402207126e-02 5.279202168223454e-01 4.941416876084902e-14 + 3.351602650938848e-02 5.221501031127233e-01 6.210669281637460e-14 + 4.037017258596558e-02 5.144469819063703e-01 6.389758657089891e-14 + 4.862601580065353e-02 5.050123590590333e-01 3.681418963040387e-14 + 5.857020818056673e-02 4.940546417356843e-01 3.407697372845994e-14 + 7.054802310718646e-02 4.817831270606888e-01 3.615272082845640e-14 + 8.497534359086456e-02 4.684030685709907e-01 2.912494195584581e-14 + 1.023531021899027e-01 4.541118003091507e-01 3.144980910790648e-14 + 1.232846739442068e-01 4.390958424788206e-01 5.581706351674976e-14 + 1.484968262254467e-01 4.235288790611684e-01 1.927874089184264e-14 + 1.788649529057435e-01 4.075704819878456e-01 3.863818027108164e-14 + 2.154434690031887e-01 3.913654536469757e-01 1.804039318856111e-14 + 2.595024211399737e-01 3.750436654344996e-01 1.638478753877823e-14 + 3.125715849688241e-01 3.587202813108329e-01 1.144229025072781e-14 + 3.764935806792472e-01 3.424962692608498e-01 1.098907948271736e-14 + 4.534878508128591e-01 3.264591183275271e-01 8.628384257483112e-15 + 5.462277217684348e-01 3.106836932620877e-01 1.370550022963008e-14 + 6.579332246575682e-01 2.952331720691905e-01 1.577348752061684e-14 + 7.924828983539186e-01 2.801600234349851e-01 8.730687111231310e-15 + 9.545484566618347e-01 2.655069910756730e-01 4.704026023152253e-15 + 1.149756995397738e+00 2.513080604538897e-01 5.330853540670894e-15 + 1.384886371393875e+00 2.375893901967694e-01 5.400481732136556e-15 + 1.668100537200059e+00 2.243701960821875e-01 5.723368239287297e-15 + 2.009233002565050e+00 2.116635798206379e-01 3.194676092401876e-15 + 2.420128264794383e+00 1.994772982363891e-01 3.663016915613689e-15 + 2.915053062825182e+00 1.878144710122765e-01 3.211327223263231e-15 + 3.511191734215135e+00 1.766742270661003e-01 2.945694384790395e-15 + 4.229242874389508e+00 1.660522909969821e-01 8.267137164157695e-12 + 5.094138014816386e+00 1.559415120020650e-01 6.102930450395757e-12 + 6.135907273413176e+00 1.463323382951112e-01 4.483869821456010e-12 + 7.390722033525790e+00 1.372132404523200e-01 3.300054740858767e-12 + 8.902150854450392e+00 1.285710873070426e-01 2.602366007129917e-12 + 1.072267222010325e+01 1.203914780853577e-01 2.687248603050329e-12 + 1.291549665014885e+01 1.126590344367703e-01 1.285249290116455e-12 + 1.555676143930475e+01 1.053576559071436e-01 1.028238133255642e-12 + 1.873817422860387e+01 9.847074225014571e-02 6.851396690422008e-13 + 2.257019719633922e+01 9.198138578711640e-02 5.016676390508252e-13 + 2.718588242732946e+01 8.587253682555671e-02 3.667481033275500e-13 + 3.274549162877732e+01 8.012714492708864e-02 3.505963119876864e-13 + 3.944206059437664e+01 7.472827862038722e-02 2.114274892193983e-13 + 4.750810162102803e+01 6.965922592333049e-02 3.149668716367259e-13 + 5.722367659350220e+01 6.490357785626068e-02 1.843708685782719e-13 + 6.892612104349709e+01 6.044529692383219e-02 1.471169005195864e-13 + 8.302175681319753e+01 5.626877236516024e-02 1.667671360512704e-13 + 1.000000000000000e+02 5.235886379952866e-02 5.945909772949099e-14 -- cgit v1.2.3-70-g09d2