mirror of
https://github.com/fjosw/pyerrors.git
synced 2025-05-15 03:53:41 +02:00
Documentation updated
This commit is contained in:
fjosw
parent
1294f8ffe0
commit
d54f3402fa
2 changed files with 44 additions and 36 deletions
|
|
@ -68,6 +68,9 @@
|
|||
<li>
|
||||
<a class="variable" href="#Obs.N_sigma_global">N_sigma_global</a>
|
||||
</li>
|
||||
<li>
|
||||
<a class="variable" href="#Obs.N_sigma_dict">N_sigma_dict</a>
|
||||
</li>
|
||||
<li>
|
||||
<a class="variable" href="#Obs.filter_eps">filter_eps</a>
|
||||
</li>
|
||||
|
|
@ -339,11 +342,14 @@
|
|||
<span class="sd"> exists this overwrites the standard value for that ensemble.</span>
|
||||
<span class="sd"> tau_exp_global : float</span>
|
||||
<span class="sd"> Standard value for tau_exp (default 0.0)</span>
|
||||
<span class="sd"> tau_exp_dict :dict</span>
|
||||
<span class="sd"> tau_exp_dict : dict</span>
|
||||
<span class="sd"> Dictionary for tau_exp values. If an entry for a given ensemble exists</span>
|
||||
<span class="sd"> this overwrites the standard value for that ensemble.</span>
|
||||
<span class="sd"> N_sigma_global : float</span>
|
||||
<span class="sd"> Standard value for N_sigma (default 1.0)</span>
|
||||
<span class="sd"> N_sigma_dict : dict</span>
|
||||
<span class="sd"> Dictionary for N_sigma values. If an entry for a given ensemble exists</span>
|
||||
<span class="sd"> this overwrites the standard value for that ensemble.</span>
|
||||
<span class="sd"> """</span>
|
||||
<span class="vm">__slots__</span> <span class="o">=</span> <span class="p">[</span><span class="s1">'names'</span><span class="p">,</span> <span class="s1">'shape'</span><span class="p">,</span> <span class="s1">'r_values'</span><span class="p">,</span> <span class="s1">'deltas'</span><span class="p">,</span> <span class="s1">'N'</span><span class="p">,</span> <span class="s1">'_value'</span><span class="p">,</span> <span class="s1">'_dvalue'</span><span class="p">,</span>
|
||||
<span class="s1">'ddvalue'</span><span class="p">,</span> <span class="s1">'reweighted'</span><span class="p">,</span> <span class="s1">'S'</span><span class="p">,</span> <span class="s1">'tau_exp'</span><span class="p">,</span> <span class="s1">'N_sigma'</span><span class="p">,</span>
|
||||
|
|
@ -356,6 +362,7 @@
|
|||
<span class="n">tau_exp_global</span> <span class="o">=</span> <span class="mf">0.0</span>
|
||||
<span class="n">tau_exp_dict</span> <span class="o">=</span> <span class="p">{}</span>
|
||||
<span class="n">N_sigma_global</span> <span class="o">=</span> <span class="mf">1.0</span>
|
||||
<span class="n">N_sigma_dict</span> <span class="o">=</span> <span class="p">{}</span>
|
||||
<span class="n">filter_eps</span> <span class="o">=</span> <span class="mf">1e-10</span>
|
||||
|
||||
<span class="k">def</span> <span class="fm">__init__</span><span class="p">(</span><span class="bp">self</span><span class="p">,</span> <span class="n">samples</span><span class="p">,</span> <span class="n">names</span><span class="p">,</span> <span class="n">idl</span><span class="o">=</span><span class="kc">None</span><span class="p">,</span> <span class="n">means</span><span class="o">=</span><span class="kc">None</span><span class="p">,</span> <span class="o">**</span><span class="n">kwargs</span><span class="p">):</span>
|
||||
|
|
@ -494,6 +501,7 @@
|
|||
|
||||
<span class="bp">self</span><span class="o">.</span><span class="n">S</span> <span class="o">=</span> <span class="p">{}</span>
|
||||
<span class="bp">self</span><span class="o">.</span><span class="n">tau_exp</span> <span class="o">=</span> <span class="p">{}</span>
|
||||
<span class="bp">self</span><span class="o">.</span><span class="n">N_sigma</span> <span class="o">=</span> <span class="p">{}</span>
|
||||
|
||||
<span class="k">if</span> <span class="n">kwargs</span><span class="o">.</span><span class="n">get</span><span class="p">(</span><span class="s1">'fft'</span><span class="p">)</span> <span class="ow">is</span> <span class="kc">False</span><span class="p">:</span>
|
||||
<span class="n">fft</span> <span class="o">=</span> <span class="kc">False</span>
|
||||
|
|
@ -504,8 +512,8 @@
|
|||
<span class="k">if</span> <span class="n">kwarg_name</span> <span class="ow">in</span> <span class="n">kwargs</span><span class="p">:</span>
|
||||
<span class="n">tmp</span> <span class="o">=</span> <span class="n">kwargs</span><span class="o">.</span><span class="n">get</span><span class="p">(</span><span class="n">kwarg_name</span><span class="p">)</span>
|
||||
<span class="k">if</span> <span class="nb">isinstance</span><span class="p">(</span><span class="n">tmp</span><span class="p">,</span> <span class="p">(</span><span class="nb">int</span><span class="p">,</span> <span class="nb">float</span><span class="p">)):</span>
|
||||
<span class="k">if</span> <span class="n">tmp</span> <span class="o"><=</span> <span class="mi">0</span><span class="p">:</span>
|
||||
<span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="n">kwarg_name</span> <span class="o">+</span> <span class="s1">' has to be larger than 0.'</span><span class="p">)</span>
|
||||
<span class="k">if</span> <span class="n">tmp</span> <span class="o"><</span> <span class="mi">0</span><span class="p">:</span>
|
||||
<span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="n">kwarg_name</span> <span class="o">+</span> <span class="s1">' has to be larger or equal to 0.'</span><span class="p">)</span>
|
||||
<span class="k">for</span> <span class="n">e</span><span class="p">,</span> <span class="n">e_name</span> <span class="ow">in</span> <span class="nb">enumerate</span><span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">e_names</span><span class="p">):</span>
|
||||
<span class="nb">getattr</span><span class="p">(</span><span class="bp">self</span><span class="p">,</span> <span class="n">kwarg_name</span><span class="p">)[</span><span class="n">e_name</span><span class="p">]</span> <span class="o">=</span> <span class="n">tmp</span>
|
||||
<span class="k">else</span><span class="p">:</span>
|
||||
|
|
@ -519,13 +527,7 @@
|
|||
|
||||
<span class="n">_parse_kwarg</span><span class="p">(</span><span class="s1">'S'</span><span class="p">)</span>
|
||||
<span class="n">_parse_kwarg</span><span class="p">(</span><span class="s1">'tau_exp'</span><span class="p">)</span>
|
||||
|
||||
<span class="k">if</span> <span class="s1">'N_sigma'</span> <span class="ow">in</span> <span class="n">kwargs</span><span class="p">:</span>
|
||||
<span class="bp">self</span><span class="o">.</span><span class="n">N_sigma</span> <span class="o">=</span> <span class="n">kwargs</span><span class="o">.</span><span class="n">get</span><span class="p">(</span><span class="s1">'N_sigma'</span><span class="p">)</span>
|
||||
<span class="k">if</span> <span class="ow">not</span> <span class="nb">isinstance</span><span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">N_sigma</span><span class="p">,</span> <span class="p">(</span><span class="nb">int</span><span class="p">,</span> <span class="nb">float</span><span class="p">)):</span>
|
||||
<span class="k">raise</span> <span class="ne">TypeError</span><span class="p">(</span><span class="s1">'N_sigma is not a number.'</span><span class="p">)</span>
|
||||
<span class="k">else</span><span class="p">:</span>
|
||||
<span class="bp">self</span><span class="o">.</span><span class="n">N_sigma</span> <span class="o">=</span> <span class="n">Obs</span><span class="o">.</span><span class="n">N_sigma_global</span>
|
||||
<span class="n">_parse_kwarg</span><span class="p">(</span><span class="s1">'N_sigma'</span><span class="p">)</span>
|
||||
|
||||
<span class="k">for</span> <span class="n">e</span><span class="p">,</span> <span class="n">e_name</span> <span class="ow">in</span> <span class="nb">enumerate</span><span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">e_names</span><span class="p">):</span>
|
||||
|
||||
|
|
@ -578,7 +580,7 @@
|
|||
<span class="c1"># Critical slowing down analysis</span>
|
||||
<span class="k">for</span> <span class="n">n</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="mi">1</span><span class="p">,</span> <span class="n">w_max</span> <span class="o">//</span> <span class="mi">2</span><span class="p">):</span>
|
||||
<span class="n">_compute_drho</span><span class="p">(</span><span class="n">n</span> <span class="o">+</span> <span class="mi">1</span><span class="p">)</span>
|
||||
<span class="k">if</span> <span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">e_rho</span><span class="p">[</span><span class="n">e_name</span><span class="p">][</span><span class="n">n</span><span class="p">]</span> <span class="o">-</span> <span class="bp">self</span><span class="o">.</span><span class="n">N_sigma</span> <span class="o">*</span> <span class="bp">self</span><span class="o">.</span><span class="n">e_drho</span><span class="p">[</span><span class="n">e_name</span><span class="p">][</span><span class="n">n</span><span class="p">])</span> <span class="o"><</span> <span class="mi">0</span> <span class="ow">or</span> <span class="n">n</span> <span class="o">>=</span> <span class="n">w_max</span> <span class="o">//</span> <span class="mi">2</span> <span class="o">-</span> <span class="mi">2</span><span class="p">:</span>
|
||||
<span class="k">if</span> <span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">e_rho</span><span class="p">[</span><span class="n">e_name</span><span class="p">][</span><span class="n">n</span><span class="p">]</span> <span class="o">-</span> <span class="bp">self</span><span class="o">.</span><span class="n">N_sigma</span><span class="p">[</span><span class="n">e_name</span><span class="p">]</span> <span class="o">*</span> <span class="bp">self</span><span class="o">.</span><span class="n">e_drho</span><span class="p">[</span><span class="n">e_name</span><span class="p">][</span><span class="n">n</span><span class="p">])</span> <span class="o"><</span> <span class="mi">0</span> <span class="ow">or</span> <span class="n">n</span> <span class="o">>=</span> <span class="n">w_max</span> <span class="o">//</span> <span class="mi">2</span> <span class="o">-</span> <span class="mi">2</span><span class="p">:</span>
|
||||
<span class="c1"># Bias correction hep-lat/0306017 eq. (49) included</span>
|
||||
<span class="bp">self</span><span class="o">.</span><span class="n">e_tauint</span><span class="p">[</span><span class="n">e_name</span><span class="p">]</span> <span class="o">=</span> <span class="bp">self</span><span class="o">.</span><span class="n">e_n_tauint</span><span class="p">[</span><span class="n">e_name</span><span class="p">][</span><span class="n">n</span><span class="p">]</span> <span class="o">*</span> <span class="p">(</span><span class="mi">1</span> <span class="o">+</span> <span class="p">(</span><span class="mi">2</span> <span class="o">*</span> <span class="n">n</span> <span class="o">+</span> <span class="mi">1</span><span class="p">)</span> <span class="o">/</span> <span class="n">e_N</span><span class="p">)</span> <span class="o">/</span> <span class="p">(</span><span class="mi">1</span> <span class="o">+</span> <span class="mi">1</span> <span class="o">/</span> <span class="n">e_N</span><span class="p">)</span> <span class="o">+</span> <span class="n">texp</span> <span class="o">*</span> <span class="n">np</span><span class="o">.</span><span class="n">abs</span><span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">e_rho</span><span class="p">[</span><span class="n">e_name</span><span class="p">][</span><span class="n">n</span> <span class="o">+</span> <span class="mi">1</span><span class="p">])</span> <span class="c1"># The absolute makes sure, that the tail contribution is always positive</span>
|
||||
<span class="bp">self</span><span class="o">.</span><span class="n">e_dtauint</span><span class="p">[</span><span class="n">e_name</span><span class="p">]</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">sqrt</span><span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">e_n_dtauint</span><span class="p">[</span><span class="n">e_name</span><span class="p">][</span><span class="n">n</span><span class="p">]</span> <span class="o">**</span> <span class="mi">2</span> <span class="o">+</span> <span class="n">texp</span> <span class="o">**</span> <span class="mi">2</span> <span class="o">*</span> <span class="bp">self</span><span class="o">.</span><span class="n">e_drho</span><span class="p">[</span><span class="n">e_name</span><span class="p">][</span><span class="n">n</span> <span class="o">+</span> <span class="mi">1</span><span class="p">]</span> <span class="o">**</span> <span class="mi">2</span><span class="p">)</span>
|
||||
|
|
@ -669,7 +671,7 @@
|
|||
<span class="k">if</span> <span class="nb">len</span><span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">e_names</span><span class="p">)</span> <span class="o">></span> <span class="mi">1</span><span class="p">:</span>
|
||||
<span class="nb">print</span><span class="p">(</span><span class="s1">''</span><span class="p">,</span> <span class="n">e_name</span><span class="p">,</span> <span class="s1">'</span><span class="se">\t</span><span class="s1"> </span><span class="si">%3.8e</span><span class="s1"> +/- </span><span class="si">%3.8e</span><span class="s1">'</span> <span class="o">%</span> <span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">e_dvalue</span><span class="p">[</span><span class="n">e_name</span><span class="p">],</span> <span class="bp">self</span><span class="o">.</span><span class="n">e_ddvalue</span><span class="p">[</span><span class="n">e_name</span><span class="p">]))</span>
|
||||
<span class="k">if</span> <span class="bp">self</span><span class="o">.</span><span class="n">tau_exp</span><span class="p">[</span><span class="n">e_name</span><span class="p">]</span> <span class="o">></span> <span class="mi">0</span><span class="p">:</span>
|
||||
<span class="nb">print</span><span class="p">(</span><span class="s1">' t_int</span><span class="se">\t</span><span class="s1"> </span><span class="si">%3.8e</span><span class="s1"> +/- </span><span class="si">%3.8e</span><span class="s1"> tau_exp = </span><span class="si">%3.2f</span><span class="s1">, N_sigma = </span><span class="si">%1.0i</span><span class="s1">'</span> <span class="o">%</span> <span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">e_tauint</span><span class="p">[</span><span class="n">e_name</span><span class="p">],</span> <span class="bp">self</span><span class="o">.</span><span class="n">e_dtauint</span><span class="p">[</span><span class="n">e_name</span><span class="p">],</span> <span class="bp">self</span><span class="o">.</span><span class="n">tau_exp</span><span class="p">[</span><span class="n">e_name</span><span class="p">],</span> <span class="bp">self</span><span class="o">.</span><span class="n">N_sigma</span><span class="p">))</span>
|
||||
<span class="nb">print</span><span class="p">(</span><span class="s1">' t_int</span><span class="se">\t</span><span class="s1"> </span><span class="si">%3.8e</span><span class="s1"> +/- </span><span class="si">%3.8e</span><span class="s1"> tau_exp = </span><span class="si">%3.2f</span><span class="s1">, N_sigma = </span><span class="si">%1.0i</span><span class="s1">'</span> <span class="o">%</span> <span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">e_tauint</span><span class="p">[</span><span class="n">e_name</span><span class="p">],</span> <span class="bp">self</span><span class="o">.</span><span class="n">e_dtauint</span><span class="p">[</span><span class="n">e_name</span><span class="p">],</span> <span class="bp">self</span><span class="o">.</span><span class="n">tau_exp</span><span class="p">[</span><span class="n">e_name</span><span class="p">],</span> <span class="bp">self</span><span class="o">.</span><span class="n">N_sigma</span><span class="p">[</span><span class="n">e_name</span><span class="p">]))</span>
|
||||
<span class="k">else</span><span class="p">:</span>
|
||||
<span class="nb">print</span><span class="p">(</span><span class="s1">' t_int</span><span class="se">\t</span><span class="s1"> </span><span class="si">%3.8e</span><span class="s1"> +/- </span><span class="si">%3.8e</span><span class="s1"> S = </span><span class="si">%3.2f</span><span class="s1">'</span> <span class="o">%</span> <span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">e_tauint</span><span class="p">[</span><span class="n">e_name</span><span class="p">],</span> <span class="bp">self</span><span class="o">.</span><span class="n">e_dtauint</span><span class="p">[</span><span class="n">e_name</span><span class="p">],</span> <span class="bp">self</span><span class="o">.</span><span class="n">S</span><span class="p">[</span><span class="n">e_name</span><span class="p">]))</span>
|
||||
<span class="k">if</span> <span class="bp">self</span><span class="o">.</span><span class="n">tag</span> <span class="ow">is</span> <span class="ow">not</span> <span class="kc">None</span><span class="p">:</span>
|
||||
|
|
@ -1868,11 +1870,14 @@
|
|||
<span class="sd"> exists this overwrites the standard value for that ensemble.</span>
|
||||
<span class="sd"> tau_exp_global : float</span>
|
||||
<span class="sd"> Standard value for tau_exp (default 0.0)</span>
|
||||
<span class="sd"> tau_exp_dict :dict</span>
|
||||
<span class="sd"> tau_exp_dict : dict</span>
|
||||
<span class="sd"> Dictionary for tau_exp values. If an entry for a given ensemble exists</span>
|
||||
<span class="sd"> this overwrites the standard value for that ensemble.</span>
|
||||
<span class="sd"> N_sigma_global : float</span>
|
||||
<span class="sd"> Standard value for N_sigma (default 1.0)</span>
|
||||
<span class="sd"> N_sigma_dict : dict</span>
|
||||
<span class="sd"> Dictionary for N_sigma values. If an entry for a given ensemble exists</span>
|
||||
<span class="sd"> this overwrites the standard value for that ensemble.</span>
|
||||
<span class="sd"> """</span>
|
||||
<span class="vm">__slots__</span> <span class="o">=</span> <span class="p">[</span><span class="s1">'names'</span><span class="p">,</span> <span class="s1">'shape'</span><span class="p">,</span> <span class="s1">'r_values'</span><span class="p">,</span> <span class="s1">'deltas'</span><span class="p">,</span> <span class="s1">'N'</span><span class="p">,</span> <span class="s1">'_value'</span><span class="p">,</span> <span class="s1">'_dvalue'</span><span class="p">,</span>
|
||||
<span class="s1">'ddvalue'</span><span class="p">,</span> <span class="s1">'reweighted'</span><span class="p">,</span> <span class="s1">'S'</span><span class="p">,</span> <span class="s1">'tau_exp'</span><span class="p">,</span> <span class="s1">'N_sigma'</span><span class="p">,</span>
|
||||
|
|
@ -1885,6 +1890,7 @@
|
|||
<span class="n">tau_exp_global</span> <span class="o">=</span> <span class="mf">0.0</span>
|
||||
<span class="n">tau_exp_dict</span> <span class="o">=</span> <span class="p">{}</span>
|
||||
<span class="n">N_sigma_global</span> <span class="o">=</span> <span class="mf">1.0</span>
|
||||
<span class="n">N_sigma_dict</span> <span class="o">=</span> <span class="p">{}</span>
|
||||
<span class="n">filter_eps</span> <span class="o">=</span> <span class="mf">1e-10</span>
|
||||
|
||||
<span class="k">def</span> <span class="fm">__init__</span><span class="p">(</span><span class="bp">self</span><span class="p">,</span> <span class="n">samples</span><span class="p">,</span> <span class="n">names</span><span class="p">,</span> <span class="n">idl</span><span class="o">=</span><span class="kc">None</span><span class="p">,</span> <span class="n">means</span><span class="o">=</span><span class="kc">None</span><span class="p">,</span> <span class="o">**</span><span class="n">kwargs</span><span class="p">):</span>
|
||||
|
|
@ -2023,6 +2029,7 @@
|
|||
|
||||
<span class="bp">self</span><span class="o">.</span><span class="n">S</span> <span class="o">=</span> <span class="p">{}</span>
|
||||
<span class="bp">self</span><span class="o">.</span><span class="n">tau_exp</span> <span class="o">=</span> <span class="p">{}</span>
|
||||
<span class="bp">self</span><span class="o">.</span><span class="n">N_sigma</span> <span class="o">=</span> <span class="p">{}</span>
|
||||
|
||||
<span class="k">if</span> <span class="n">kwargs</span><span class="o">.</span><span class="n">get</span><span class="p">(</span><span class="s1">'fft'</span><span class="p">)</span> <span class="ow">is</span> <span class="kc">False</span><span class="p">:</span>
|
||||
<span class="n">fft</span> <span class="o">=</span> <span class="kc">False</span>
|
||||
|
|
@ -2033,8 +2040,8 @@
|
|||
<span class="k">if</span> <span class="n">kwarg_name</span> <span class="ow">in</span> <span class="n">kwargs</span><span class="p">:</span>
|
||||
<span class="n">tmp</span> <span class="o">=</span> <span class="n">kwargs</span><span class="o">.</span><span class="n">get</span><span class="p">(</span><span class="n">kwarg_name</span><span class="p">)</span>
|
||||
<span class="k">if</span> <span class="nb">isinstance</span><span class="p">(</span><span class="n">tmp</span><span class="p">,</span> <span class="p">(</span><span class="nb">int</span><span class="p">,</span> <span class="nb">float</span><span class="p">)):</span>
|
||||
<span class="k">if</span> <span class="n">tmp</span> <span class="o"><=</span> <span class="mi">0</span><span class="p">:</span>
|
||||
<span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="n">kwarg_name</span> <span class="o">+</span> <span class="s1">' has to be larger than 0.'</span><span class="p">)</span>
|
||||
<span class="k">if</span> <span class="n">tmp</span> <span class="o"><</span> <span class="mi">0</span><span class="p">:</span>
|
||||
<span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="n">kwarg_name</span> <span class="o">+</span> <span class="s1">' has to be larger or equal to 0.'</span><span class="p">)</span>
|
||||
<span class="k">for</span> <span class="n">e</span><span class="p">,</span> <span class="n">e_name</span> <span class="ow">in</span> <span class="nb">enumerate</span><span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">e_names</span><span class="p">):</span>
|
||||
<span class="nb">getattr</span><span class="p">(</span><span class="bp">self</span><span class="p">,</span> <span class="n">kwarg_name</span><span class="p">)[</span><span class="n">e_name</span><span class="p">]</span> <span class="o">=</span> <span class="n">tmp</span>
|
||||
<span class="k">else</span><span class="p">:</span>
|
||||
|
|
@ -2048,13 +2055,7 @@
|
|||
|
||||
<span class="n">_parse_kwarg</span><span class="p">(</span><span class="s1">'S'</span><span class="p">)</span>
|
||||
<span class="n">_parse_kwarg</span><span class="p">(</span><span class="s1">'tau_exp'</span><span class="p">)</span>
|
||||
|
||||
<span class="k">if</span> <span class="s1">'N_sigma'</span> <span class="ow">in</span> <span class="n">kwargs</span><span class="p">:</span>
|
||||
<span class="bp">self</span><span class="o">.</span><span class="n">N_sigma</span> <span class="o">=</span> <span class="n">kwargs</span><span class="o">.</span><span class="n">get</span><span class="p">(</span><span class="s1">'N_sigma'</span><span class="p">)</span>
|
||||
<span class="k">if</span> <span class="ow">not</span> <span class="nb">isinstance</span><span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">N_sigma</span><span class="p">,</span> <span class="p">(</span><span class="nb">int</span><span class="p">,</span> <span class="nb">float</span><span class="p">)):</span>
|
||||
<span class="k">raise</span> <span class="ne">TypeError</span><span class="p">(</span><span class="s1">'N_sigma is not a number.'</span><span class="p">)</span>
|
||||
<span class="k">else</span><span class="p">:</span>
|
||||
<span class="bp">self</span><span class="o">.</span><span class="n">N_sigma</span> <span class="o">=</span> <span class="n">Obs</span><span class="o">.</span><span class="n">N_sigma_global</span>
|
||||
<span class="n">_parse_kwarg</span><span class="p">(</span><span class="s1">'N_sigma'</span><span class="p">)</span>
|
||||
|
||||
<span class="k">for</span> <span class="n">e</span><span class="p">,</span> <span class="n">e_name</span> <span class="ow">in</span> <span class="nb">enumerate</span><span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">e_names</span><span class="p">):</span>
|
||||
|
||||
|
|
@ -2107,7 +2108,7 @@
|
|||
<span class="c1"># Critical slowing down analysis</span>
|
||||
<span class="k">for</span> <span class="n">n</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="mi">1</span><span class="p">,</span> <span class="n">w_max</span> <span class="o">//</span> <span class="mi">2</span><span class="p">):</span>
|
||||
<span class="n">_compute_drho</span><span class="p">(</span><span class="n">n</span> <span class="o">+</span> <span class="mi">1</span><span class="p">)</span>
|
||||
<span class="k">if</span> <span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">e_rho</span><span class="p">[</span><span class="n">e_name</span><span class="p">][</span><span class="n">n</span><span class="p">]</span> <span class="o">-</span> <span class="bp">self</span><span class="o">.</span><span class="n">N_sigma</span> <span class="o">*</span> <span class="bp">self</span><span class="o">.</span><span class="n">e_drho</span><span class="p">[</span><span class="n">e_name</span><span class="p">][</span><span class="n">n</span><span class="p">])</span> <span class="o"><</span> <span class="mi">0</span> <span class="ow">or</span> <span class="n">n</span> <span class="o">>=</span> <span class="n">w_max</span> <span class="o">//</span> <span class="mi">2</span> <span class="o">-</span> <span class="mi">2</span><span class="p">:</span>
|
||||
<span class="k">if</span> <span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">e_rho</span><span class="p">[</span><span class="n">e_name</span><span class="p">][</span><span class="n">n</span><span class="p">]</span> <span class="o">-</span> <span class="bp">self</span><span class="o">.</span><span class="n">N_sigma</span><span class="p">[</span><span class="n">e_name</span><span class="p">]</span> <span class="o">*</span> <span class="bp">self</span><span class="o">.</span><span class="n">e_drho</span><span class="p">[</span><span class="n">e_name</span><span class="p">][</span><span class="n">n</span><span class="p">])</span> <span class="o"><</span> <span class="mi">0</span> <span class="ow">or</span> <span class="n">n</span> <span class="o">>=</span> <span class="n">w_max</span> <span class="o">//</span> <span class="mi">2</span> <span class="o">-</span> <span class="mi">2</span><span class="p">:</span>
|
||||
<span class="c1"># Bias correction hep-lat/0306017 eq. (49) included</span>
|
||||
<span class="bp">self</span><span class="o">.</span><span class="n">e_tauint</span><span class="p">[</span><span class="n">e_name</span><span class="p">]</span> <span class="o">=</span> <span class="bp">self</span><span class="o">.</span><span class="n">e_n_tauint</span><span class="p">[</span><span class="n">e_name</span><span class="p">][</span><span class="n">n</span><span class="p">]</span> <span class="o">*</span> <span class="p">(</span><span class="mi">1</span> <span class="o">+</span> <span class="p">(</span><span class="mi">2</span> <span class="o">*</span> <span class="n">n</span> <span class="o">+</span> <span class="mi">1</span><span class="p">)</span> <span class="o">/</span> <span class="n">e_N</span><span class="p">)</span> <span class="o">/</span> <span class="p">(</span><span class="mi">1</span> <span class="o">+</span> <span class="mi">1</span> <span class="o">/</span> <span class="n">e_N</span><span class="p">)</span> <span class="o">+</span> <span class="n">texp</span> <span class="o">*</span> <span class="n">np</span><span class="o">.</span><span class="n">abs</span><span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">e_rho</span><span class="p">[</span><span class="n">e_name</span><span class="p">][</span><span class="n">n</span> <span class="o">+</span> <span class="mi">1</span><span class="p">])</span> <span class="c1"># The absolute makes sure, that the tail contribution is always positive</span>
|
||||
<span class="bp">self</span><span class="o">.</span><span class="n">e_dtauint</span><span class="p">[</span><span class="n">e_name</span><span class="p">]</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">sqrt</span><span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">e_n_dtauint</span><span class="p">[</span><span class="n">e_name</span><span class="p">][</span><span class="n">n</span><span class="p">]</span> <span class="o">**</span> <span class="mi">2</span> <span class="o">+</span> <span class="n">texp</span> <span class="o">**</span> <span class="mi">2</span> <span class="o">*</span> <span class="bp">self</span><span class="o">.</span><span class="n">e_drho</span><span class="p">[</span><span class="n">e_name</span><span class="p">][</span><span class="n">n</span> <span class="o">+</span> <span class="mi">1</span><span class="p">]</span> <span class="o">**</span> <span class="mi">2</span><span class="p">)</span>
|
||||
|
|
@ -2198,7 +2199,7 @@
|
|||
<span class="k">if</span> <span class="nb">len</span><span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">e_names</span><span class="p">)</span> <span class="o">></span> <span class="mi">1</span><span class="p">:</span>
|
||||
<span class="nb">print</span><span class="p">(</span><span class="s1">''</span><span class="p">,</span> <span class="n">e_name</span><span class="p">,</span> <span class="s1">'</span><span class="se">\t</span><span class="s1"> </span><span class="si">%3.8e</span><span class="s1"> +/- </span><span class="si">%3.8e</span><span class="s1">'</span> <span class="o">%</span> <span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">e_dvalue</span><span class="p">[</span><span class="n">e_name</span><span class="p">],</span> <span class="bp">self</span><span class="o">.</span><span class="n">e_ddvalue</span><span class="p">[</span><span class="n">e_name</span><span class="p">]))</span>
|
||||
<span class="k">if</span> <span class="bp">self</span><span class="o">.</span><span class="n">tau_exp</span><span class="p">[</span><span class="n">e_name</span><span class="p">]</span> <span class="o">></span> <span class="mi">0</span><span class="p">:</span>
|
||||
<span class="nb">print</span><span class="p">(</span><span class="s1">' t_int</span><span class="se">\t</span><span class="s1"> </span><span class="si">%3.8e</span><span class="s1"> +/- </span><span class="si">%3.8e</span><span class="s1"> tau_exp = </span><span class="si">%3.2f</span><span class="s1">, N_sigma = </span><span class="si">%1.0i</span><span class="s1">'</span> <span class="o">%</span> <span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">e_tauint</span><span class="p">[</span><span class="n">e_name</span><span class="p">],</span> <span class="bp">self</span><span class="o">.</span><span class="n">e_dtauint</span><span class="p">[</span><span class="n">e_name</span><span class="p">],</span> <span class="bp">self</span><span class="o">.</span><span class="n">tau_exp</span><span class="p">[</span><span class="n">e_name</span><span class="p">],</span> <span class="bp">self</span><span class="o">.</span><span class="n">N_sigma</span><span class="p">))</span>
|
||||
<span class="nb">print</span><span class="p">(</span><span class="s1">' t_int</span><span class="se">\t</span><span class="s1"> </span><span class="si">%3.8e</span><span class="s1"> +/- </span><span class="si">%3.8e</span><span class="s1"> tau_exp = </span><span class="si">%3.2f</span><span class="s1">, N_sigma = </span><span class="si">%1.0i</span><span class="s1">'</span> <span class="o">%</span> <span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">e_tauint</span><span class="p">[</span><span class="n">e_name</span><span class="p">],</span> <span class="bp">self</span><span class="o">.</span><span class="n">e_dtauint</span><span class="p">[</span><span class="n">e_name</span><span class="p">],</span> <span class="bp">self</span><span class="o">.</span><span class="n">tau_exp</span><span class="p">[</span><span class="n">e_name</span><span class="p">],</span> <span class="bp">self</span><span class="o">.</span><span class="n">N_sigma</span><span class="p">[</span><span class="n">e_name</span><span class="p">]))</span>
|
||||
<span class="k">else</span><span class="p">:</span>
|
||||
<span class="nb">print</span><span class="p">(</span><span class="s1">' t_int</span><span class="se">\t</span><span class="s1"> </span><span class="si">%3.8e</span><span class="s1"> +/- </span><span class="si">%3.8e</span><span class="s1"> S = </span><span class="si">%3.2f</span><span class="s1">'</span> <span class="o">%</span> <span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">e_tauint</span><span class="p">[</span><span class="n">e_name</span><span class="p">],</span> <span class="bp">self</span><span class="o">.</span><span class="n">e_dtauint</span><span class="p">[</span><span class="n">e_name</span><span class="p">],</span> <span class="bp">self</span><span class="o">.</span><span class="n">S</span><span class="p">[</span><span class="n">e_name</span><span class="p">]))</span>
|
||||
<span class="k">if</span> <span class="bp">self</span><span class="o">.</span><span class="n">tag</span> <span class="ow">is</span> <span class="ow">not</span> <span class="kc">None</span><span class="p">:</span>
|
||||
|
|
@ -2572,6 +2573,9 @@ Dictionary for tau_exp values. If an entry for a given ensemble exists
|
|||
this overwrites the standard value for that ensemble.</li>
|
||||
<li><strong>N_sigma_global</strong> (float):
|
||||
Standard value for N_sigma (default 1.0)</li>
|
||||
<li><strong>N_sigma_dict</strong> (dict):
|
||||
Dictionary for N_sigma values. If an entry for a given ensemble exists
|
||||
this overwrites the standard value for that ensemble.</li>
|
||||
</ul>
|
||||
</div>
|
||||
|
||||
|
|
@ -2730,6 +2734,15 @@ already subtracted from the samples</li>
|
|||
|
||||
|
||||
|
||||
</div>
|
||||
<div id="Obs.N_sigma_dict" class="classattr">
|
||||
<div class="attr variable"><a class="headerlink" href="#Obs.N_sigma_dict">#  </a>
|
||||
|
||||
<span class="name">N_sigma_dict</span><span class="default_value"> = {}</span>
|
||||
</div>
|
||||
|
||||
|
||||
|
||||
</div>
|
||||
<div id="Obs.filter_eps" class="classattr">
|
||||
<div class="attr variable"><a class="headerlink" href="#Obs.filter_eps">#  </a>
|
||||
|
|
@ -2912,6 +2925,7 @@ already subtracted from the samples</li>
|
|||
|
||||
<span class="bp">self</span><span class="o">.</span><span class="n">S</span> <span class="o">=</span> <span class="p">{}</span>
|
||||
<span class="bp">self</span><span class="o">.</span><span class="n">tau_exp</span> <span class="o">=</span> <span class="p">{}</span>
|
||||
<span class="bp">self</span><span class="o">.</span><span class="n">N_sigma</span> <span class="o">=</span> <span class="p">{}</span>
|
||||
|
||||
<span class="k">if</span> <span class="n">kwargs</span><span class="o">.</span><span class="n">get</span><span class="p">(</span><span class="s1">'fft'</span><span class="p">)</span> <span class="ow">is</span> <span class="kc">False</span><span class="p">:</span>
|
||||
<span class="n">fft</span> <span class="o">=</span> <span class="kc">False</span>
|
||||
|
|
@ -2922,8 +2936,8 @@ already subtracted from the samples</li>
|
|||
<span class="k">if</span> <span class="n">kwarg_name</span> <span class="ow">in</span> <span class="n">kwargs</span><span class="p">:</span>
|
||||
<span class="n">tmp</span> <span class="o">=</span> <span class="n">kwargs</span><span class="o">.</span><span class="n">get</span><span class="p">(</span><span class="n">kwarg_name</span><span class="p">)</span>
|
||||
<span class="k">if</span> <span class="nb">isinstance</span><span class="p">(</span><span class="n">tmp</span><span class="p">,</span> <span class="p">(</span><span class="nb">int</span><span class="p">,</span> <span class="nb">float</span><span class="p">)):</span>
|
||||
<span class="k">if</span> <span class="n">tmp</span> <span class="o"><=</span> <span class="mi">0</span><span class="p">:</span>
|
||||
<span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="n">kwarg_name</span> <span class="o">+</span> <span class="s1">' has to be larger than 0.'</span><span class="p">)</span>
|
||||
<span class="k">if</span> <span class="n">tmp</span> <span class="o"><</span> <span class="mi">0</span><span class="p">:</span>
|
||||
<span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="n">kwarg_name</span> <span class="o">+</span> <span class="s1">' has to be larger or equal to 0.'</span><span class="p">)</span>
|
||||
<span class="k">for</span> <span class="n">e</span><span class="p">,</span> <span class="n">e_name</span> <span class="ow">in</span> <span class="nb">enumerate</span><span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">e_names</span><span class="p">):</span>
|
||||
<span class="nb">getattr</span><span class="p">(</span><span class="bp">self</span><span class="p">,</span> <span class="n">kwarg_name</span><span class="p">)[</span><span class="n">e_name</span><span class="p">]</span> <span class="o">=</span> <span class="n">tmp</span>
|
||||
<span class="k">else</span><span class="p">:</span>
|
||||
|
|
@ -2937,13 +2951,7 @@ already subtracted from the samples</li>
|
|||
|
||||
<span class="n">_parse_kwarg</span><span class="p">(</span><span class="s1">'S'</span><span class="p">)</span>
|
||||
<span class="n">_parse_kwarg</span><span class="p">(</span><span class="s1">'tau_exp'</span><span class="p">)</span>
|
||||
|
||||
<span class="k">if</span> <span class="s1">'N_sigma'</span> <span class="ow">in</span> <span class="n">kwargs</span><span class="p">:</span>
|
||||
<span class="bp">self</span><span class="o">.</span><span class="n">N_sigma</span> <span class="o">=</span> <span class="n">kwargs</span><span class="o">.</span><span class="n">get</span><span class="p">(</span><span class="s1">'N_sigma'</span><span class="p">)</span>
|
||||
<span class="k">if</span> <span class="ow">not</span> <span class="nb">isinstance</span><span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">N_sigma</span><span class="p">,</span> <span class="p">(</span><span class="nb">int</span><span class="p">,</span> <span class="nb">float</span><span class="p">)):</span>
|
||||
<span class="k">raise</span> <span class="ne">TypeError</span><span class="p">(</span><span class="s1">'N_sigma is not a number.'</span><span class="p">)</span>
|
||||
<span class="k">else</span><span class="p">:</span>
|
||||
<span class="bp">self</span><span class="o">.</span><span class="n">N_sigma</span> <span class="o">=</span> <span class="n">Obs</span><span class="o">.</span><span class="n">N_sigma_global</span>
|
||||
<span class="n">_parse_kwarg</span><span class="p">(</span><span class="s1">'N_sigma'</span><span class="p">)</span>
|
||||
|
||||
<span class="k">for</span> <span class="n">e</span><span class="p">,</span> <span class="n">e_name</span> <span class="ow">in</span> <span class="nb">enumerate</span><span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">e_names</span><span class="p">):</span>
|
||||
|
||||
|
|
@ -2996,7 +3004,7 @@ already subtracted from the samples</li>
|
|||
<span class="c1"># Critical slowing down analysis</span>
|
||||
<span class="k">for</span> <span class="n">n</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="mi">1</span><span class="p">,</span> <span class="n">w_max</span> <span class="o">//</span> <span class="mi">2</span><span class="p">):</span>
|
||||
<span class="n">_compute_drho</span><span class="p">(</span><span class="n">n</span> <span class="o">+</span> <span class="mi">1</span><span class="p">)</span>
|
||||
<span class="k">if</span> <span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">e_rho</span><span class="p">[</span><span class="n">e_name</span><span class="p">][</span><span class="n">n</span><span class="p">]</span> <span class="o">-</span> <span class="bp">self</span><span class="o">.</span><span class="n">N_sigma</span> <span class="o">*</span> <span class="bp">self</span><span class="o">.</span><span class="n">e_drho</span><span class="p">[</span><span class="n">e_name</span><span class="p">][</span><span class="n">n</span><span class="p">])</span> <span class="o"><</span> <span class="mi">0</span> <span class="ow">or</span> <span class="n">n</span> <span class="o">>=</span> <span class="n">w_max</span> <span class="o">//</span> <span class="mi">2</span> <span class="o">-</span> <span class="mi">2</span><span class="p">:</span>
|
||||
<span class="k">if</span> <span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">e_rho</span><span class="p">[</span><span class="n">e_name</span><span class="p">][</span><span class="n">n</span><span class="p">]</span> <span class="o">-</span> <span class="bp">self</span><span class="o">.</span><span class="n">N_sigma</span><span class="p">[</span><span class="n">e_name</span><span class="p">]</span> <span class="o">*</span> <span class="bp">self</span><span class="o">.</span><span class="n">e_drho</span><span class="p">[</span><span class="n">e_name</span><span class="p">][</span><span class="n">n</span><span class="p">])</span> <span class="o"><</span> <span class="mi">0</span> <span class="ow">or</span> <span class="n">n</span> <span class="o">>=</span> <span class="n">w_max</span> <span class="o">//</span> <span class="mi">2</span> <span class="o">-</span> <span class="mi">2</span><span class="p">:</span>
|
||||
<span class="c1"># Bias correction hep-lat/0306017 eq. (49) included</span>
|
||||
<span class="bp">self</span><span class="o">.</span><span class="n">e_tauint</span><span class="p">[</span><span class="n">e_name</span><span class="p">]</span> <span class="o">=</span> <span class="bp">self</span><span class="o">.</span><span class="n">e_n_tauint</span><span class="p">[</span><span class="n">e_name</span><span class="p">][</span><span class="n">n</span><span class="p">]</span> <span class="o">*</span> <span class="p">(</span><span class="mi">1</span> <span class="o">+</span> <span class="p">(</span><span class="mi">2</span> <span class="o">*</span> <span class="n">n</span> <span class="o">+</span> <span class="mi">1</span><span class="p">)</span> <span class="o">/</span> <span class="n">e_N</span><span class="p">)</span> <span class="o">/</span> <span class="p">(</span><span class="mi">1</span> <span class="o">+</span> <span class="mi">1</span> <span class="o">/</span> <span class="n">e_N</span><span class="p">)</span> <span class="o">+</span> <span class="n">texp</span> <span class="o">*</span> <span class="n">np</span><span class="o">.</span><span class="n">abs</span><span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">e_rho</span><span class="p">[</span><span class="n">e_name</span><span class="p">][</span><span class="n">n</span> <span class="o">+</span> <span class="mi">1</span><span class="p">])</span> <span class="c1"># The absolute makes sure, that the tail contribution is always positive</span>
|
||||
<span class="bp">self</span><span class="o">.</span><span class="n">e_dtauint</span><span class="p">[</span><span class="n">e_name</span><span class="p">]</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">sqrt</span><span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">e_n_dtauint</span><span class="p">[</span><span class="n">e_name</span><span class="p">][</span><span class="n">n</span><span class="p">]</span> <span class="o">**</span> <span class="mi">2</span> <span class="o">+</span> <span class="n">texp</span> <span class="o">**</span> <span class="mi">2</span> <span class="o">*</span> <span class="bp">self</span><span class="o">.</span><span class="n">e_drho</span><span class="p">[</span><span class="n">e_name</span><span class="p">][</span><span class="n">n</span> <span class="o">+</span> <span class="mi">1</span><span class="p">]</span> <span class="o">**</span> <span class="mi">2</span><span class="p">)</span>
|
||||
|
|
@ -3106,7 +3114,7 @@ of the autocorrelation function (default True)</li>
|
|||
<span class="k">if</span> <span class="nb">len</span><span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">e_names</span><span class="p">)</span> <span class="o">></span> <span class="mi">1</span><span class="p">:</span>
|
||||
<span class="nb">print</span><span class="p">(</span><span class="s1">''</span><span class="p">,</span> <span class="n">e_name</span><span class="p">,</span> <span class="s1">'</span><span class="se">\t</span><span class="s1"> </span><span class="si">%3.8e</span><span class="s1"> +/- </span><span class="si">%3.8e</span><span class="s1">'</span> <span class="o">%</span> <span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">e_dvalue</span><span class="p">[</span><span class="n">e_name</span><span class="p">],</span> <span class="bp">self</span><span class="o">.</span><span class="n">e_ddvalue</span><span class="p">[</span><span class="n">e_name</span><span class="p">]))</span>
|
||||
<span class="k">if</span> <span class="bp">self</span><span class="o">.</span><span class="n">tau_exp</span><span class="p">[</span><span class="n">e_name</span><span class="p">]</span> <span class="o">></span> <span class="mi">0</span><span class="p">:</span>
|
||||
<span class="nb">print</span><span class="p">(</span><span class="s1">' t_int</span><span class="se">\t</span><span class="s1"> </span><span class="si">%3.8e</span><span class="s1"> +/- </span><span class="si">%3.8e</span><span class="s1"> tau_exp = </span><span class="si">%3.2f</span><span class="s1">, N_sigma = </span><span class="si">%1.0i</span><span class="s1">'</span> <span class="o">%</span> <span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">e_tauint</span><span class="p">[</span><span class="n">e_name</span><span class="p">],</span> <span class="bp">self</span><span class="o">.</span><span class="n">e_dtauint</span><span class="p">[</span><span class="n">e_name</span><span class="p">],</span> <span class="bp">self</span><span class="o">.</span><span class="n">tau_exp</span><span class="p">[</span><span class="n">e_name</span><span class="p">],</span> <span class="bp">self</span><span class="o">.</span><span class="n">N_sigma</span><span class="p">))</span>
|
||||
<span class="nb">print</span><span class="p">(</span><span class="s1">' t_int</span><span class="se">\t</span><span class="s1"> </span><span class="si">%3.8e</span><span class="s1"> +/- </span><span class="si">%3.8e</span><span class="s1"> tau_exp = </span><span class="si">%3.2f</span><span class="s1">, N_sigma = </span><span class="si">%1.0i</span><span class="s1">'</span> <span class="o">%</span> <span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">e_tauint</span><span class="p">[</span><span class="n">e_name</span><span class="p">],</span> <span class="bp">self</span><span class="o">.</span><span class="n">e_dtauint</span><span class="p">[</span><span class="n">e_name</span><span class="p">],</span> <span class="bp">self</span><span class="o">.</span><span class="n">tau_exp</span><span class="p">[</span><span class="n">e_name</span><span class="p">],</span> <span class="bp">self</span><span class="o">.</span><span class="n">N_sigma</span><span class="p">[</span><span class="n">e_name</span><span class="p">]))</span>
|
||||
<span class="k">else</span><span class="p">:</span>
|
||||
<span class="nb">print</span><span class="p">(</span><span class="s1">' t_int</span><span class="se">\t</span><span class="s1"> </span><span class="si">%3.8e</span><span class="s1"> +/- </span><span class="si">%3.8e</span><span class="s1"> S = </span><span class="si">%3.2f</span><span class="s1">'</span> <span class="o">%</span> <span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">e_tauint</span><span class="p">[</span><span class="n">e_name</span><span class="p">],</span> <span class="bp">self</span><span class="o">.</span><span class="n">e_dtauint</span><span class="p">[</span><span class="n">e_name</span><span class="p">],</span> <span class="bp">self</span><span class="o">.</span><span class="n">S</span><span class="p">[</span><span class="n">e_name</span><span class="p">]))</span>
|
||||
<span class="k">if</span> <span class="bp">self</span><span class="o">.</span><span class="n">tag</span> <span class="ow">is</span> <span class="ow">not</span> <span class="kc">None</span><span class="p">:</span>
|
||||
|
|
|
|||
File diff suppressed because one or more lines are too long
Loading…
Add table
Add a link
Reference in a new issue