29 files changed, 6837 insertions, 34 deletions
diff --git a/content.tex b/content.tex
index 0111c9d..3a9152d 100644
--- a/content.tex
+++ b/content.tex
@@ -850,9 +850,8 @@ Mit \(u_\text{max}:=1\) vereinfacht sich damit die analytische Lösung der \(x\)
 Das Geschwindigkeitsprofil des Poiseuille-Flusses ist also parabelförmig.
 \end{Definition}
 
-Wir wollen in einem \(1 \times 4\) Längeneinheiten bemessenden Rohr einen Poiseuille-Fluss simulieren. Als Auflösung einer Längeneinheit sei dabei \(N=10\) gewählt, was in der Gitterdiskretisierung durch \(11 \times 21\) grobe und \(21 \times 43\) feine Knoten abgebildet wird.
-
-In Abbildung~\ref{fig:PoiseuilleGridSetup} sehen wir das resultierende Gitter zusammen mit den zugewiesenen Materialzahlen. Wand- und Einflusszellen werden nach dieser Vorlage mit lokalen Geschwindigkeitsrandbedingungen umgesetzt. Während für den Einfluss dabei das Geschwindigkeitsprofil der analytischen Poiseuille-Lösung vorausgesetzt wird, erhält der Ausfluss eine Druckrandbedingung. Die noch verbleibenden gewöhnlichen Fluidzellen erfahren, abgesehen von den üblichen Kollisions- und Strömungsschritten, keine besondere Behandlung.
+Wir wollen in einem \(1 \times 4\) Meter bemessenden Rohr einen Poiseuille-Fluss simulieren. Als Auflösung einer Längeneinheit sei dabei \(N=10\) gewählt, was in der Diskretisierung durch \(11 \times 21\) grobe und \(21 \times 43\) feine Knoten abgebildet wird.
+In Abbildung~\ref{fig:PoiseuilleGridSetup} sehen wir das resultierende Gitter zusammen mit den zugewiesenen Materialzahlen. Wand- und Einflusszellen werden nach dieser Vorlage mit lokalen Geschwindigkeitsrandbedingungen umgesetzt. Während für den Einfluss dabei das Geschwindigkeitsprofil der analytischen Poiseuille-Lösung vorausgesetzt wird, erhält der Ausfluss eine Druckrandbedingung.
 
 \begin{figure}[H]
 \begin{adjustbox}{center}
@@ -862,7 +861,7 @@ In Abbildung~\ref{fig:PoiseuilleGridSetup} sehen wir das resultierende Gitter zu
 \label{fig:PoiseuilleGridSetup}
 \end{figure}
 
-Neben diesen knotenspezifischen Eigenschaften sei \(u=0.01\) die charateristische Geschwindigkeit in Lattice-Einheiten und \(\text{Re}=10\) die modellierte Reynolds-Zahl. Erstellen wir unsere grobe \class{Grid2D} Instanz mit diesen, die Relaxionszeit \(\tau\) fixierenden, Werten und führen die Simulation bis zur Konvergenz aus, erblicken wir bei geeigneter Aufbereitung in ParaView~\cite{paraview05} schließlich das in Abbildung~\ref{fig:PoiseuilleVelocityGrid} ersichtliche Strömungsbild. Konvergenz bedeutet in diesem Kontext, dass die durchschnittliche Energie des feinen Gitters unter einen Residuumswert, hier \(1\mathrm{e}{-5}\), gefallen ist.
+Neben diesen knotenspezifischen Eigenschaften sei \(u=\num[round-mode=off]{0.01}\) die charateristische Geschwindigkeit in Lattice-Einheiten und \(\text{Re}=10\) die modellierte Reynolds-Zahl. Erstellen wir unsere grobe \class{Grid2D} Instanz mit diesen, die Relaxionszeit \(\tau\) fixierenden, Werten und führen die Simulation bis zur Konvergenz aus, erblicken wir bei geeigneter Aufbereitung in ParaView~\cite{paraview05} schließlich das in Abbildung~\ref{fig:PoiseuilleVelocityGrid} ersichtliche Strömungsbild. Konvergenz bedeutet in diesem Kontext, dass die durchschnittliche Energie des feinen Gitters unter einen Residuumswert, hier \num{1e-5}, gefallen ist.
 
 \begin{figure}[h]
 \begin{adjustbox}{center}
@@ -872,7 +871,7 @@ Neben diesen knotenspezifischen Eigenschaften sei \(u=0.01\) die charateristisch
 \label{fig:PoiseuilleVelocityGrid}
 \end{figure}
 
-Bei erster Betrachtung lässt sich erkennen, dass die Strömung den Gitterübergang subjektiv ideal bestritten hat. Es treten also keine ungewöhnlichen Artefakte im Geschwindigkeitsbild auf und dieses setzt sich nach dem Übergang in, bis auf die neuen Zwischenwerte, unveränderter Weise fort. Tatsächlich ist bei Interpolation der Knotenzwischenbereiche zur Bildung einer geschlossenen Fläche kein Gitterübergang erkennbar.
+Bei erster Betrachtung lässt sich erkennen, dass die Strömung den Gitterübergang subjektiv ideal bestritten hat. Es treten also keine ungewöhnlichen Artefakte im Geschwindigkeitsbild auf und dieses setzt sich nach dem Übergang in, bis auf die neuen Zwischenwerte, unveränderter Weise fort. Tatsächlich ist bei Bildung einer geschlossenen Fläche durch Interpolation der Zwischenbereiche kein Gitterübergang erkennbar.
 
 \newpage
 \subsubsection{Vergleich mit der analytischen Lösung}
@@ -1143,14 +1142,14 @@ Auch im Vergleich mit diesem 145314 Knoten umfassenden Gitter bewähren sich die
 & Uniform & Verfeinert & Referenzintervall \cite{SchaeferTurek96} \\
 \hline
 \hline
-\(\widehat{c_w}\) & 3.28433 & 3.21927 & \([3.22,\ 3.24]\) \\
-\(|\widehat{c_w}-3.23|\) & 0.05433 & 0.01073 & [0, 0.01] \\
+\(\widehat{c_w}\) & \num{3.28433} & \num{3.21927} & \([\num{3.22},\ \num{3.24}]\) \\
+\(|\widehat{c_w}-\num{3.23}|\) & \num{0.05433} & \num{0.01073} & [\num{0}, \num{0.01}] \\
 \hline
-\(\widehat{c_a}\) & 1.07046 & 1.10359 & \([0.99,\ 1.01]\) \\
-\(|\widehat{c_a}-1.0|\) & 0.07046 & 0.10359 & [0, 0.01]  \\
+\(\widehat{c_a}\) & \num{1.07046} & \num{1.10359} & \([\num{0.99},\ \num{1.01}]\) \\
+\(|\widehat{c_a}-1.0|\) & \num{0.07046} & \num{0.10359} & [\num{0}, \num{0.01}]  \\
 \hline
-\(\Delta P\) & 2.5793 & 2.44285 & \([2.46,\ 2.5]\)  \\
-\(|\Delta P-2.48|\) & 0.0993 & 0.03715 & [0, 0.02] \\
+\(\Delta P\) & \num{2.5793} & \num{2.44285} & \([\num{2.46},\ \num{2.5}]\)  \\
+\(|\Delta P-2.48|\) & \num{0.0993} & \num{0.03715} & [\num{0}, \num{0.02}] \\
 \hline
 \hline
 Knotenanzahl & 145314 & 13454 & \\
@@ -1175,5 +1174,10 @@ Dazu sehen wir in Abbildung~\ref{fig:142000nodes} die charakteristischen Messwer
 \label{fig:142000nodes}
 \end{figure}
 
+\begin{figure}[H]
+\centering
+\input{img/cylinder2d_generously_refined_comparison.tikz}
+\end{figure}
+
 \newpage
 \section{Fazit}
diff --git a/img/common/velocity_fluid.gnuplot b/img/common/velocity_fluid.gnuplot
index 5535111..c0db2a9 100644
--- a/img/common/velocity_fluid.gnuplot
+++ b/img/common/velocity_fluid.gnuplot
@@ -7,6 +7,7 @@ set output 'tmp/'.plotname.'.png'
 load 'common/moreland.pal'
 
 set datafile separator ','
+set decimalsign ','
 
 set size 1,1
 set margin 0,0,0,0
diff --git a/img/cylinder2d_deltap_comparison.tikz b/img/cylinder2d_deltap_comparison.tikz
index 97ecd37..2d6fd99 100644
--- a/img/cylinder2d_deltap_comparison.tikz
+++ b/img/cylinder2d_deltap_comparison.tikz
@@ -15,7 +15,8 @@
 	xlabel={Simulierte physikalische Zeit},
 	ylabel={Druckdifferenz},
 	x unit=s,
-	y unit=N/m^2
+	y unit=N/m^2,
+	y tick label style={/pgf/number format/.cd, use comma}
 ]
 
 \addplot[
@@ -35,7 +36,7 @@
 \addlegendentry {Problembezogen verfeinertes \(N=5\) Gitter};
 
 \addplot[color=black]{2.48};
-\addlegendentry {Gemittelte Referenzlösung \(\Delta P := 2.48\)};
+\addlegendentry {Gemittelte Referenzlösung \(\Delta P := \num{2.48}\)};
 
 \end{axis}
 \end{tikzpicture}
diff --git a/img/cylinder2d_drag_comparison.tikz b/img/cylinder2d_drag_comparison.tikz
index 98336ec..ee0c783 100644
--- a/img/cylinder2d_drag_comparison.tikz
+++ b/img/cylinder2d_drag_comparison.tikz
@@ -34,7 +34,7 @@
 \addlegendentry {Problembezogen verfeinertes \(N=5\) Gitter};
 
 \addplot[color=black]{3.23};
-\addlegendentry {Gemittelte Referenzlösung \(c_\text{Dmax} := 3.23\)};
+\addlegendentry {Gemittelte Referenzlösung \(c_\text{Dmax} := \num{3.23}\)};
 
 \end{axis}
 \end{tikzpicture}
diff --git a/img/cylinder2d_high_res_comparsion.tikz b/img/cylinder2d_high_res_comparsion.tikz
index aa8b2a7..9bde49b 100644
--- a/img/cylinder2d_high_res_comparsion.tikz
+++ b/img/cylinder2d_high_res_comparsion.tikz
@@ -1,6 +1,7 @@
 \begin{tikzpicture}
 \pgfplotstableread[col sep=comma]{img/data/cylinder2d_optimized_refinement_n40_re100_drag_lift_deltap.csv}\refined
-\pgfplotstableread[col sep=comma]{img/data/cylinder2d_unrefined_n40_re100_drag_lift_deltap.csv}\uniformVeryHigh
+\pgfplotstableread[col sep=comma]{img/data/cylinder2d_unrefined_n40_re100_drag_lift_deltap.csv}\uniformHigh
+\pgfplotstableread[col sep=comma]{img/data/cylinder2d_unrefined_n80_re100_drag_lift_deltap.csv}\uniformVeryHigh
 \pgfplotstableread[col sep=comma]{img/data/cylinder2d_unrefined_n160_re100_drag_lift_deltap.csv}\uniformVeryVeryHigh
 
 \begin{axis}[
@@ -15,11 +16,13 @@
 	ylabel={Widerstandskoeffizient},
 	ylabel absolute, every axis y label/.append style={yshift=0.4cm}
 ]
+
 \addplot[thick,color=black]{3.23};
-\addlegendentry {\(\widehat{c_w} := 3.23\)};
+\addlegendentry {\(\widehat{c_w} := \num[round-mode=off]{3.23}\)};
 \addplot[color=blue!50!white,thin ] table [x expr=\thisrow{time}, y=drag] {\uniformVeryVeryHigh};
 \addplot[color=red!50!white,thin  ] table [x expr=\thisrow{time}, y=drag] {\uniformVeryHigh};
 \addplot[color=green!70!black,thin] table [x expr=8*\thisrow{time}, y=drag] {\refined};
+
 \end{axis}
 
 \begin{axis}[
@@ -35,11 +38,13 @@
 	ylabel={Auftriebskoeffizient},
 	ylabel absolute, every axis y label/.append style={yshift=0.4cm}
 ]
+
 \addplot[thick,color=black]{1.0};
 \addlegendentry {\(\widehat{c_a} := 1\)};
 \addplot[color=blue!50!white,thin ] table [x expr=\thisrow{time}, y=lift] {\uniformVeryVeryHigh};
 \addplot[color=red!50!white,thin  ] table [x expr=\thisrow{time}, y=lift] {\uniformVeryHigh};
 \addplot[color=green!70!black,thin] table [x expr=8*\thisrow{time}, y=lift] {\refined};
+
 \end{axis}
 
 \begin{axis}[
@@ -58,14 +63,16 @@
 	y unit=N/m^2,
 	ylabel absolute, every axis y label/.append style={yshift=0.4cm}
 ]
+
+\addplot[thick,color=black]{2.48};
+\addlegendentry {\(\Delta P := \num[round-mode=off]{2.48}\)};
 \addplot[color=blue!50!white,thin  ] table [x expr=\thisrow{time}, y=deltap] {\uniformVeryVeryHigh};
-\addlegendentry {Uniformes \(N=160\) Gitter mit 2298014 Knoten};
+\addlegendentry {Uniformes \(N=160\) Gitter mit \(\sim 2300000\) Knoten};
 \addplot[color=red!50!white,thin   ] table [x expr=\thisrow{time}, y=deltap] {\uniformVeryHigh};
-\addlegendentry {Uniformes \(N=40\) Gitter mit 576758 Knoten};
+\addlegendentry {Uniformes \(N=80\) Gitter mit \(\sim 577000\) Knoten};
 \addplot[color=green!70!black,thin] table [x expr=8*\thisrow{time}, y=deltap] {\refined};
-\addlegendentry {Verfeinertes \(N=40\) Gitter mit 787592 Knoten};
-\addplot[thick,color=black]{2.48};
-\addlegendentry {\(\Delta P := 2.48\)};
+\addlegendentry {Verfeinertes \(N=40\) Gitter mit \(\sim 790000\) Knoten};
+
 \end{axis}
 
 \end{tikzpicture}
diff --git a/img/cylinder2d_overview.tikz b/img/cylinder2d_overview.tikz
index 3fa0da8..d77c1c1 100644
--- a/img/cylinder2d_overview.tikz
+++ b/img/cylinder2d_overview.tikz
@@ -10,8 +10,8 @@
 \draw[very thick,fill=gray!20!white] (2,2) circle (0.5);
 
 \draw[Latex-Latex,thin] (3,1.5) -- (3,2.5) node[pos=0.5,right] {\(D\)};
-\draw[Latex-Latex,thin] (3,0) -- (3,1.5) node[pos=0.5,right] {\(1.5 D\)};
-\draw[Latex-Latex,thin] (3,2.5) -- (3,4.1) node[pos=0.5,right] {\(1.6 D\)};
+\draw[Latex-Latex,thin] (3,0) -- (3,1.5) node[pos=0.5,right] {\(\num{1.5} D\)};
+\draw[Latex-Latex,thin] (3,2.5) -- (3,4.1) node[pos=0.5,right] {\(\num{1.6} D\)};
 
 \draw[Latex-Latex,thin] (0,-0.5) -- (2,-0.5) node[pos=0.5,below] {\(2 D\)};
 \draw[Latex-Latex,thin] (2,-0.5) -- (22,-0.5) node[pos=0.5,below] {\(20 D\)};
diff --git a/img/data/cylinder2d_generously_refined_n15_re100_drag_lift_deltap.csv b/img/data/cylinder2d_generously_refined_n15_re100_drag_lift_deltap.csv
new file mode 100644
index 0000000..df3ffdf
--- /dev/null
+++ b/img/data/cylinder2d_generously_refined_n15_re100_drag_lift_deltap.csv
@@ -0,0 +1,1613 @@
+time,drag,lift,deltap
+0,0,0,0
+0.00124074,2.16707e-07,-1.01615e-09,1.41561e-07
+0.00248148,5.66404e-06,-2.32384e-08,3.48643e-06
+0.00372222,2.07988e-05,-4.13144e-08,1.26355e-05
+0.00496296,4.51912e-05,-7.25249e-08,2.64518e-05
+0.0062037,8.05954e-05,-1.30149e-07,4.81813e-05
+0.00744444,0.000125882,-1.51336e-07,7.29354e-05
+0.00868519,0.000179798,-2.59974e-07,0.000104258
+0.00992593,0.000249092,-2.85972e-07,0.000144944
+0.0111667,0.000322546,-3.82066e-07,0.000182332
+0.0124074,0.000413101,-5.36216e-07,0.000238647
+0.0136481,0.000513467,-4.832e-07,0.000289239
+0.0148889,0.000618901,-8.09399e-07,0.000349125
+0.0161296,0.000750258,-7.48793e-07,0.000425646
+0.0173704,0.000874775,-9.12487e-07,0.000483313
+0.0186111,0.00102628,-1.23767e-06,0.000579434
+0.0198519,0.00118623,-9.88473e-07,0.000654675
+0.0210926,0.00134556,-1.62793e-06,0.000743902
+0.0223333,0.00155243,-1.41798e-06,0.000864309
+0.0235741,0.00174183,-1.65263e-06,0.000946437
+0.0248148,0.00197321,-2.19395e-06,0.00109526
+0.0260556,0.00221027,-1.6465e-06,0.00120166
+0.0272963,0.00243669,-2.70767e-06,0.00132701
+0.028537,0.002725,-2.26102e-06,0.00149418
+0.0297778,0.00296786,-2.56064e-06,0.00158973
+0.0310185,0.00326876,-3.37514e-06,0.00178677
+0.0322593,0.00357276,-2.42635e-06,0.00191507
+0.0335,0.00385458,-3.98675e-06,0.00206945
+0.0347407,0.00422257,-3.24489e-06,0.00228287
+0.0359815,0.00451762,-3.62553e-06,0.00238808
+0.0372222,0.00489117,-4.7446e-06,0.00263737
+0.038463,0.00526355,-3.3024e-06,0.00278564
+0.0397037,0.00559963,-5.4718e-06,0.00296833
+0.0409444,0.00605,-4.35728e-06,0.00323001
+0.0421852,0.00639465,-4.81179e-06,0.00334038
+0.0434259,0.00683807,-6.30127e-06,0.00364226
+0.0446667,0.00726866,-4.2744e-06,0.00380186
+0.0459074,0.00764363,-7.12343e-06,0.00400396
+0.0471481,0.00816632,-5.58634e-06,0.00430873
+0.0483889,0.00854746,-6.1322e-06,0.00441428
+0.0496296,0.00905768,-8.02073e-06,0.00476959
+0.0508704,0.0095552,-5.3359e-06,0.00494393
+0.0521111,0.00998171,-8.96114e-06,0.00517342
+0.0533519,0.0105919,-6.93372e-06,0.00553075
+0.0545926,0.0110188,-7.57015e-06,0.00563438
+0.0558333,0.0116049,-9.90253e-06,0.00605014
+0.0570741,0.0121718,-6.4789e-06,0.00623808
+0.0583148,0.0126479,-1.09423e-05,0.00649377
+0.0595556,0.0133479,-8.36819e-06,0.00690521
+0.0607963,0.0138195,-9.11157e-06,0.00700391
+0.062037,0.014484,-1.19074e-05,0.00748331
+0.0632778,0.0151227,-7.66954e-06,0.00768405
+0.0645185,0.0156511,-1.30517e-05,0.00796762
+0.0657593,0.0164518,-9.86902e-06,0.00843995
+0.067,0.01698,-1.07163e-05,0.0085384
+0.0682407,0.0177372,-1.40013e-05,0.00909125
+0.0694815,0.0184582,-8.87881e-06,0.00930843
+0.0707222,0.0190443,-1.52312e-05,0.00962269
+0.071963,0.019944,-1.13646e-05,0.0101544
+0.0732037,0.0205145,-1.2334e-05,0.0102417
+0.0744444,0.0213522,-1.61241e-05,0.0108622
+0.0756852,0.0221451,-1.00421e-05,0.0110877
+0.0769259,0.0227795,-1.74376e-05,0.0114277
+0.0781667,0.0237761,-1.28308e-05,0.0120179
+0.0794074,0.0243875,-1.39253e-05,0.0120918
+0.0806481,0.0253069,-1.82353e-05,0.0127823
+0.0818889,0.0261708,-1.11377e-05,0.013014
+0.0831296,0.0268508,-1.96262e-05,0.0133787
+0.0843704,0.027943,-1.42167e-05,0.0140268
+0.0856111,0.0285895,-1.54563e-05,0.0140825
+0.0868519,0.0295842,-2.02904e-05,0.0148408
+0.0880926,0.0305066,-1.21208e-05,0.0150697
+0.0893333,0.0312174,-2.17588e-05,0.0154511
+0.0905741,0.032396,-1.54919e-05,0.0161524
+0.0918148,0.0330688,-1.68894e-05,0.0161839
+0.0930556,0.0341394,-2.2253e-05,0.0170124
+0.0942963,0.0351289,-1.2973e-05,0.0172427
+0.095537,0.0358826,-2.38004e-05,0.017649
+0.0967778,0.0371586,-1.66095e-05,0.01841
+0.0980185,0.0378634,-1.81946e-05,0.0184195
+0.0992593,0.0390147,-2.40723e-05,0.0193224
+0.1005,0.0400709,-1.36239e-05,0.0195521
+0.101741,0.0408651,-2.56872e-05,0.0199825
+0.102981,0.0422392,-1.75037e-05,0.0208036
+0.104222,0.042975,-1.92885e-05,0.0207897
+0.105463,0.0442098,-2.56592e-05,0.0217697
+0.106704,0.0453352,-1.40047e-05,0.0219989
+0.107944,0.0461746,-2.73359e-05,0.0224571
+0.109185,0.0476567,-1.80808e-05,0.0233438
+0.110426,0.0484344,-2.00899e-05,0.0233122
+0.111667,0.0497636,-2.69241e-05,0.0243766
+0.112907,0.0509635,-1.4011e-05,0.0246067
+0.114148,0.0518488,-2.86271e-05,0.0250936
+0.115389,0.0534327,-1.82255e-05,0.0260415
+0.11663,0.0542379,-2.04757e-05,0.0259829
+0.11787,0.0556505,-2.77208e-05,0.0271263
+0.119111,0.0569158,-1.35184e-05,0.0273509
+0.120352,0.0578412,-2.94426e-05,0.0278645
+0.121593,0.0595255,-1.78016e-05,0.0288728
+0.122833,0.0603573,-2.03241e-05,0.0287867
+0.124074,0.0618542,-2.79341e-05,0.030011
+0.125315,0.0631824,-1.23953e-05,0.0302278
+0.126556,0.0641452,-2.96405e-05,0.0307678
+0.127796,0.0659263,-1.66796e-05,0.0318342
+0.129037,0.0667778,-1.94912e-05,0.0317167
+0.130278,0.0683517,-2.73988e-05,0.0330189
+0.131519,0.0697303,-1.04963e-05,0.0332198
+0.132759,0.0707192,-2.90649e-05,0.0337816
+0.134,0.0725904,-1.46902e-05,0.0349022
+0.135241,0.0734578,-1.78219e-05,0.0347523
+0.136481,0.0751134,-2.59559e-05,0.0361368
+0.137722,0.0765515,-7.65761e-06,0.0363268
+0.138963,0.0775787,-2.75388e-05,0.0369195
+0.140204,0.0795487,-1.16557e-05,0.0380992
+0.141444,0.0804362,-1.51265e-05,0.0379199
+0.142685,0.0821772,-2.3385e-05,0.0393898
+0.143926,0.0836735,-3.66457e-06,0.0395677
+0.145167,0.084739,-2.48355e-05,0.0401931
+0.146407,0.0868092,-7.32944e-06,0.0414321
+0.147648,0.0877182,-1.11659e-05,0.0412254
+0.148889,0.0895499,-1.94374e-05,0.0427851
+0.15013,0.091109,1.73986e-06,0.0429529
+0.15137,0.0922212,-2.06875e-05,0.0436182
+0.152611,0.0944018,-1.43766e-06,0.0449221
+0.153852,0.0953435,-5.65762e-06,0.0446956
+0.155093,0.0972761,-1.3807e-05,0.0463517
+0.156333,0.0989004,8.86168e-06,0.0465103
+0.157574,0.100059,-1.47759e-05,0.0472176
+0.158815,0.102344,6.35132e-06,0.0485818
+0.160056,0.103308,1.73404e-06,0.0483308
+0.161296,0.105335,-6.15241e-06,0.0500808
+0.162537,0.107019,1.80359e-05,0.0502274
+0.163778,0.108225,-6.74843e-06,0.0509797
+0.165019,0.110615,1.63867e-05,0.0524048
+0.166259,0.111603,1.13511e-05,0.052133
+0.1675,0.113728,3.89304e-06,0.0539805
+0.168741,0.115472,2.96285e-05,0.0541149
+0.169981,0.116726,3.7627e-06,0.0549158
+0.171222,0.11922,2.90447e-05,0.0564004
+0.172463,0.120228,2.35806e-05,0.0561081
+0.173704,0.122448,1.67192e-05,0.0580527
+0.174944,0.124242,4.40192e-05,0.0581709
+0.176185,0.125541,1.71588e-05,0.0590224
+0.177426,0.128137,4.47214e-05,0.0605665
+0.178667,0.12917,3.88034e-05,0.0602601
+0.179907,0.131496,3.2729e-05,0.0623107
+0.181148,0.133354,6.16013e-05,0.0624213
+0.182389,0.134714,3.38364e-05,0.0633368
+0.18363,0.137424,6.3826e-05,0.0649464
+0.18487,0.138489,5.74355e-05,0.0646338
+0.186111,0.140929,5.23454e-05,0.0667966
+0.187352,0.142854,8.27933e-05,0.0669027
+0.188593,0.144283,5.42296e-05,0.0678908
+0.189833,0.147112,8.67805e-05,0.0695699
+0.191074,0.148219,7.98869e-05,0.0692606
+0.192315,0.150786,7.59937e-05,0.0715448
+0.193556,0.152788,0.000108,0.0716543
+0.194796,0.154302,7.87408e-05,0.072729
+0.196037,0.157267,0.000113987,0.0744867
+0.197278,0.158431,0.000106555,0.0741943
+0.198519,0.16114,0.000104067,0.0766099
+0.199759,0.163225,0.000137599,0.0767278
+0.201,0.16483,0.000107754,0.0778974
+0.202241,0.167929,0.000145808,0.0797333
+0.203481,0.169151,0.000137778,0.0794623
+0.204722,0.172004,0.000136895,0.0820138
+0.205963,0.174176,0.000171886,0.0821448
+0.207204,0.175884,0.000141544,0.0834213
+0.208444,0.179126,0.000182487,0.0853412
+0.209685,0.180415,0.000173775,0.0851037
+0.210926,0.183425,0.000174679,0.0877999
+0.212167,0.18569,0.000211018,0.0879505
+0.213407,0.18751,0.000180254,0.089345
+0.214648,0.1909,0.000224131,0.0913523
+0.215889,0.192261,0.000214617,0.0911574
+0.21713,0.195432,0.000217456,0.0940042
+0.21837,0.197792,0.000254983,0.0941787
+0.219611,0.199734,0.000223841,0.0957028
+0.220852,0.203281,0.000270635,0.0978044
+0.222093,0.204731,0.000260156,0.0976687
+0.223333,0.208086,0.000265042,0.100681
+0.224574,0.210559,0.000303528,0.100894
+0.225815,0.212648,0.000272025,0.102568
+0.227056,0.216368,0.000321668,0.104774
+0.228296,0.217923,0.000310009,0.104713
+0.229537,0.221476,0.000317009,0.107902
+0.230778,0.224074,0.000356169,0.108164
+0.232019,0.226327,0.00032427,0.110003
+0.233259,0.230235,0.000376622,0.112322
+0.2345,0.231913,0.000363516,0.112353
+0.235741,0.235686,0.000372638,0.115732
+0.236981,0.238427,0.000412111,0.116057
+0.238222,0.240868,0.000379736,0.11808
+0.239463,0.244983,0.000434588,0.120524
+0.240704,0.246806,0.00041971,0.120666
+0.241944,0.250817,0.000430902,0.124245
+0.243185,0.253711,0.000470261,0.124642
+0.244426,0.256352,0.000437279,0.126861
+0.245667,0.260677,0.000494338,0.12943
+0.246907,0.262652,0.000477314,0.129692
+0.248148,0.266909,0.00049045,0.133479
+0.249389,0.269965,0.000529189,0.133955
+0.25063,0.272823,0.000495418,0.136383
+0.25187,0.27737,0.000554312,0.139085
+0.253111,0.27951,0.000534713,0.139481
+0.254352,0.284025,0.000549607,0.14348
+0.255593,0.287249,0.00058716,0.144042
+0.256833,0.290334,0.000552374,0.146689
+0.258074,0.295105,0.000612661,0.149522
+0.259315,0.297417,0.000590017,0.150061
+0.260556,0.302194,0.000606422,0.154275
+0.261796,0.305588,0.000642169,0.154927
+0.263037,0.30891,0.000606098,0.1578
+0.264278,0.313912,0.00066727,0.16077
+0.265519,0.316411,0.000641083,0.161465
+0.266759,0.321467,0.000658693,0.165902
+0.268,0.325044,0.000691973,0.166654
+0.269241,0.328619,0.000654324,0.169767
+0.270481,0.333859,0.000715818,0.172875
+0.271722,0.336553,0.000685565,0.173737
+0.272963,0.341892,0.000704033,0.178397
+0.274204,0.345656,0.00073414,0.179256
+0.275444,0.349494,0.000694581,0.182615
+0.276685,0.354975,0.000755775,0.185863
+0.277926,0.357875,0.000720902,0.186903
+0.279167,0.363506,0.000739816,0.191793
+0.280407,0.367463,0.000765998,0.192765
+0.281648,0.371575,0.000724154,0.196381
+0.282889,0.377303,0.000784353,0.199771
+0.28413,0.380415,0.000744265,0.201
+0.28537,0.38634,0.000763142,0.206118
+0.286611,0.390485,0.000784576,0.207206
+0.287852,0.394869,0.00074001,0.211081
+0.289093,0.400832,0.000798427,0.214609
+0.290333,0.404153,0.000752458,0.216031
+0.291574,0.410366,0.000770716,0.221376
+0.292815,0.414697,0.000786494,0.222584
+0.294056,0.419355,0.000738687,0.226723
+0.295296,0.425548,0.000794425,0.230388
+0.296537,0.429078,0.000741845,0.232012
+0.297778,0.435575,0.000758805,0.237582
+0.299019,0.440085,0.000767937,0.238912
+0.300259,0.445013,0.000716304,0.243318
+0.3015,0.451426,0.000768375,0.247116
+0.302741,0.455161,0.000708402,0.248946
+0.303981,0.461933,0.000723312,0.254737
+0.305222,0.466615,0.000724756,0.256193
+0.306463,0.471815,0.000668679,0.26087
+0.307704,0.478445,0.000716055,0.264803
+0.308944,0.482392,0.000647914,0.266852
+0.310185,0.489444,0.000660016,0.272867
+0.311426,0.494302,0.000652768,0.274458
+0.312667,0.499779,0.000591684,0.279414
+0.313907,0.506624,0.00063339,0.283482
+0.315148,0.510785,0.000556418,0.285758
+0.316389,0.518115,0.000565046,0.291996
+0.31763,0.523149,0.000548242,0.293728
+0.31887,0.528908,0.000481743,0.298967
+0.320111,0.535967,0.000516971,0.303173
+0.321352,0.540351,0.000430726,0.305687
+0.322593,0.547962,0.000435427,0.31215
+0.323833,0.553178,0.000408473,0.314033
+0.325074,0.559226,0.000336431,0.319562
+0.326315,0.5665,0.000364667,0.323907
+0.327556,0.571111,0.00026906,0.326669
+0.328796,0.579,0.000269706,0.333353
+0.330037,0.584392,0.000232391,0.33539
+0.331278,0.590726,0.000155067,0.341208
+0.332519,0.598205,0.000176191,0.345689
+0.333759,0.603043,7.157e-05,0.348702
+0.335,0.611204,6.84533e-05,0.355605
+0.336241,0.616773,2.10385e-05,0.3578
+0.337481,0.623394,-6.08444e-05,0.363909
+0.338722,0.631072,-4.6476e-05,0.368525
+0.339963,0.636137,-0.000159245,0.371794
+0.341204,0.644564,-0.000165353,0.37891
+0.342444,0.650302,-0.000222078,0.381266
+0.343685,0.657205,-0.000307288,0.387663
+0.344926,0.665072,-0.000298809,0.392409
+0.346167,0.670359,-0.000418328,0.395938
+0.347407,0.679042,-0.00042616,0.403259
+0.348648,0.684946,-0.000490873,0.405779
+0.349889,0.692129,-0.00057768,0.412463
+0.35113,0.700182,-0.000573721,0.41734
+0.35237,0.705697,-0.000698079,0.421135
+0.353611,0.714635,-0.000705903,0.428659
+0.354852,0.720706,-0.000776786,0.431349
+0.356093,0.72817,-0.000863008,0.43832
+0.357333,0.736401,-0.000861764,0.443324
+0.358574,0.742144,-0.000988612,0.447389
+0.359815,0.751329,-0.000994324,0.455109
+0.361056,0.757564,-0.00106917,0.457974
+0.362296,0.765307,-0.00115229,0.465229
+0.363537,0.773713,-0.00115165,0.47036
+0.364778,0.779687,-0.00127835,0.474698
+0.366019,0.789116,-0.00127962,0.482611
+0.367259,0.795516,-0.00135599,0.485654
+0.3685,0.803538,-0.00143337,0.493191
+0.369741,0.812112,-0.00143108,0.498448
+0.370981,0.818314,-0.00155492,0.503061
+0.372222,0.827974,-0.00154943,0.511156
+0.373463,0.834531,-0.00162489,0.514377
+0.374704,0.84282,-0.00169398,0.522186
+0.375944,0.851548,-0.00168796,0.52756
+0.377185,0.857972,-0.00180641,0.532446
+0.378426,0.867852,-0.00179213,0.540714
+0.379667,0.87456,-0.00186457,0.544113
+0.380907,0.883108,-0.00192323,0.552188
+0.382148,0.891979,-0.00191183,0.557675
+0.383389,0.898619,-0.00202289,0.56283
+0.38463,0.908703,-0.00199839,0.571258
+0.38587,0.915551,-0.00206633,0.574833
+0.387111,0.924346,-0.00211317,0.583162
+0.388352,0.933345,-0.00209553,0.588754
+0.389593,0.94019,-0.00219807,0.594172
+0.390833,0.950464,-0.00216285,0.60275
+0.392074,0.957446,-0.00222577,0.606499
+0.393315,0.966478,-0.00226047,0.615073
+0.394556,0.975598,-0.00223681,0.620768
+0.395796,0.98265,-0.0023308,0.626449
+0.397037,0.993104,-0.00228559,0.635167
+0.398278,1.00022,-0.00234416,0.639091
+0.399519,1.00948,-0.00236767,0.647901
+0.400759,1.01871,-0.00233948,0.653692
+0.402,1.02596,-0.00242622,0.659631
+0.403241,1.03658,-0.00237305,0.668478
+0.404481,1.04382,-0.00242924,0.672577
+0.405722,1.0533,-0.00244384,0.681613
+0.406963,1.06263,-0.0024139,0.687497
+0.408204,1.07008,-0.00249595,0.693692
+0.409444,1.08086,-0.00243813,0.702659
+0.410685,1.08822,-0.00249508,0.706932
+0.411926,1.09792,-0.00250424,0.716184
+0.413167,1.10734,-0.00247637,0.722156
+0.414407,1.11498,-0.00255732,0.728598
+0.415648,1.1259,-0.00249905,0.737672
+0.416889,1.13337,-0.00256071,0.742113
+0.41813,1.14326,-0.00256862,0.751565
+0.41937,1.15276,-0.00254723,0.757616
+0.420611,1.16057,-0.00263113,0.764297
+0.421852,1.17162,-0.00257699,0.773464
+0.423093,1.17919,-0.00264748,0.77807
+0.424333,1.18927,-0.00265838,0.787708
+0.425574,1.19883,-0.00264774,0.793831
+0.426815,1.20681,-0.00273835,0.800739
+0.428056,1.21796,-0.00269243,0.809985
+0.429296,1.22562,-0.00277519,0.814748
+0.430537,1.23586,-0.00279253,0.824554
+0.431778,1.24547,-0.00279587,0.83074
+0.433019,1.25361,-0.00289579,0.837864
+0.434259,1.26484,-0.00286082,0.847174
+0.4355,1.27259,-0.00295774,0.852089
+0.436741,1.28298,-0.00298328,0.86205
+0.437981,1.29263,-0.00300194,0.868294
+0.439222,1.30091,-0.00311174,0.875625
+0.440463,1.31222,-0.00308828,0.884988
+0.441704,1.32004,-0.00319896,0.890051
+0.442944,1.33056,-0.00323207,0.900151
+0.444185,1.34024,-0.00326484,0.906446
+0.445426,1.34865,-0.00338255,0.913972
+0.446667,1.36003,-0.00336853,0.923376
+0.447907,1.36791,-0.00348986,0.928582
+0.449148,1.37856,-0.00352727,0.938806
+0.450389,1.38826,-0.0035703,0.945147
+0.45163,1.3968,-0.00369135,0.95286
+0.45287,1.40822,-0.00368217,0.962292
+0.454111,1.41617,-0.00380864,0.967636
+0.455352,1.42693,-0.00384482,0.97797
+0.456593,1.43663,-0.0038921,0.984349
+0.457833,1.44528,-0.00401001,0.992233
+0.459074,1.45673,-0.00399937,1.00168
+0.460315,1.46472,-0.004124,1.00715
+0.461556,1.47557,-0.00415225,1.01757
+0.462796,1.48526,-0.0041969,1.02398
+0.464037,1.49402,-0.00430468,1.03202
+0.465278,1.50548,-0.00428607,1.04147
+0.466519,1.51351,-0.00440206,1.04706
+0.467759,1.52443,-0.00441626,1.05755
+0.469,1.5341,-0.0044523,1.06398
+0.470241,1.54294,-0.0045443,1.07216
+0.471481,1.55438,-0.00451292,1.08159
+0.472722,1.56243,-0.00461558,1.08729
+0.473963,1.5734,-0.00461213,1.09783
+0.475204,1.58303,-0.00463642,1.10427
+0.476444,1.59194,-0.00471026,1.11257
+0.477685,1.60335,-0.00466487,1.12198
+0.478926,1.61142,-0.00475335,1.12777
+0.480167,1.62242,-0.00473279,1.13835
+0.481407,1.632,-0.00474653,1.14479
+0.482648,1.64096,-0.00480435,1.1532
+0.483889,1.65234,-0.00474833,1.16256
+0.48513,1.6604,-0.00482652,1.16844
+0.48637,1.67142,-0.00479411,1.17904
+0.487611,1.68094,-0.00480316,1.18547
+0.488852,1.68994,-0.00485172,1.19398
+0.490093,1.70125,-0.00479285,1.20328
+0.491333,1.70932,-0.0048688,1.20924
+0.492574,1.72033,-0.00483369,1.21984
+0.493815,1.72978,-0.00484737,1.22626
+0.495056,1.73881,-0.00489648,1.23484
+0.496296,1.75006,-0.00484506,1.24408
+0.497537,1.7581,-0.00492881,1.25012
+0.498778,1.7691,-0.00490151,1.26069
+0.500019,1.77847,-0.00492987,1.26709
+0.501259,1.78751,-0.00498939,1.27573
+0.5025,1.79866,-0.004955,1.28489
+0.503741,1.80668,-0.00505516,1.29098
+0.504981,1.81764,-0.00504402,1.30151
+0.506222,1.8269,-0.00509411,1.30789
+0.507463,1.83595,-0.00517024,1.31657
+0.508704,1.84699,-0.00515803,1.32563
+0.509944,1.85496,-0.00527813,1.33176
+0.511185,1.86586,-0.00528573,1.34224
+0.512426,1.87502,-0.00535828,1.34857
+0.513667,1.88403,-0.00545037,1.35727
+0.514907,1.89495,-0.00545822,1.36623
+0.516148,1.90286,-0.00559418,1.37239
+0.517389,1.91368,-0.0056153,1.38278
+0.51863,1.92271,-0.0057031,1.38907
+0.51987,1.93168,-0.00580265,1.39776
+0.521111,1.94245,-0.00582064,1.4066
+0.522352,1.9503,-0.00596092,1.41278
+0.523593,1.96101,-0.00598325,1.42308
+0.524833,1.96991,-0.00607252,1.4293
+0.526074,1.97883,-0.00616516,1.43798
+0.527315,1.98943,-0.0061783,1.44668
+0.528556,1.9972,-0.00630722,1.45286
+0.529796,2.0078,-0.00631535,1.46305
+0.531037,2.01655,-0.00639036,1.4692
+0.532278,2.0254,-0.00646099,1.47785
+0.533519,2.03583,-0.00645481,1.48641
+0.534759,2.04351,-0.00655858,1.49258
+0.536,2.05398,-0.00654033,1.50264
+0.537241,2.06257,-0.00658993,1.50872
+0.538481,2.07133,-0.00662941,1.51731
+0.539722,2.08158,-0.00659665,1.52572
+0.540963,2.08916,-0.00666997,1.53187
+0.542204,2.09948,-0.00662273,1.54179
+0.543444,2.10791,-0.00664636,1.54778
+0.544685,2.11657,-0.00665699,1.5563
+0.545926,2.12662,-0.00660249,1.56455
+0.547167,2.13408,-0.00665267,1.57066
+0.548407,2.14423,-0.00658671,1.58041
+0.549648,2.15248,-0.00659671,1.58631
+0.550889,2.16102,-0.00659353,1.59474
+0.55213,2.17085,-0.0065344,1.60281
+0.55337,2.17818,-0.0065803,1.60887
+0.554611,2.18815,-0.00651647,1.61844
+0.555852,2.19621,-0.00653458,1.62424
+0.557093,2.2046,-0.00654057,1.63256
+0.558333,2.2142,-0.00650013,1.64044
+0.559574,2.22139,-0.00656505,1.64643
+0.560815,2.23115,-0.00652654,1.6558
+0.562056,2.23901,-0.00657464,1.66148
+0.563296,2.24725,-0.00661068,1.66968
+0.564537,2.2566,-0.00660781,1.67736
+0.565778,2.26364,-0.00670823,1.68327
+0.567019,2.27317,-0.00670906,1.69243
+0.568259,2.28082,-0.00679767,1.69799
+0.5695,2.28889,-0.00687111,1.70604
+0.570741,2.29799,-0.0069092,1.71351
+0.571981,2.30486,-0.00704446,1.71934
+0.573222,2.31416,-0.00708008,1.72827
+0.574463,2.32159,-0.00720001,1.7337
+0.575704,2.32948,-0.00729782,1.74159
+0.576944,2.33831,-0.00735953,1.74885
+0.578185,2.34501,-0.00750832,1.75457
+0.579426,2.35406,-0.00755406,1.76327
+0.580667,2.36127,-0.00767705,1.76855
+0.581907,2.36896,-0.00776857,1.77627
+0.583148,2.37752,-0.0078208,1.78331
+0.584389,2.38404,-0.0079485,1.78892
+0.58563,2.39283,-0.00796905,1.79736
+0.58687,2.3998,-0.00805933,1.8025
+0.588111,2.40728,-0.00810979,1.81003
+0.589352,2.41556,-0.00811905,1.81684
+0.590593,2.42188,-0.00819445,1.82232
+0.591833,2.43039,-0.00816189,1.8305
+0.593074,2.43712,-0.00819507,1.83549
+0.594315,2.44439,-0.00818501,1.84281
+0.595556,2.45237,-0.0081368,1.84938
+0.596796,2.45848,-0.00815122,1.85472
+0.598037,2.46671,-0.00806325,1.86263
+0.599278,2.47318,-0.00804328,1.86745
+0.600519,2.48021,-0.00798368,1.87455
+0.601759,2.48788,-0.00789552,1.88087
+0.603,2.49378,-0.00787319,1.88605
+0.604241,2.5017,-0.00776054,1.89368
+0.605481,2.50791,-0.00772376,1.89833
+0.606722,2.51469,-0.00765662,1.9052
+0.607963,2.52205,-0.00757478,1.91127
+0.609204,2.52772,-0.00756568,1.91629
+0.610444,2.53533,-0.00748101,1.92362
+0.611685,2.54129,-0.00748126,1.92809
+0.612926,2.5478,-0.0074609,1.93472
+0.614167,2.55484,-0.00743833,1.94053
+0.615407,2.56028,-0.00749303,1.94537
+0.616648,2.56758,-0.00748319,1.9524
+0.617889,2.57326,-0.00756193,1.9567
+0.61913,2.57951,-0.0076237,1.96307
+0.62037,2.58623,-0.00768797,1.96862
+0.621611,2.59143,-0.00782556,1.97328
+0.622852,2.5984,-0.00790025,1.98
+0.624093,2.60382,-0.00805687,1.98411
+0.625333,2.60979,-0.0081899,1.99022
+0.626574,2.61618,-0.00831913,1.99551
+0.627815,2.62115,-0.00850716,1.99998
+0.629056,2.62778,-0.00862365,2.00639
+0.630296,2.63292,-0.00880549,2.01031
+0.631537,2.63862,-0.00894854,2.01614
+0.632778,2.64467,-0.00907359,2.02116
+0.634019,2.64939,-0.00923681,2.02544
+0.635259,2.65569,-0.00931528,2.03152
+0.6365,2.66055,-0.0094396,2.03525
+0.637741,2.66595,-0.00950946,2.04081
+0.638981,2.67168,-0.00954829,2.