mirror of
https://github.com/fjosw/pyerrors.git
synced 2025-05-14 19:43:41 +02:00
Documentation updated
This commit is contained in:
fjosw
parent
983a6095db
commit
9442fda122
3 changed files with 67 additions and 67 deletions
|
|
@ -336,7 +336,7 @@
|
|||
</span><span id="L-119"><a href="#L-119"><span class="linenos"> 119</span></a> <span class="k">def</span> <span class="nf">gamma_method</span><span class="p">(</span><span class="bp">self</span><span class="p">,</span> <span class="o">**</span><span class="n">kwargs</span><span class="p">):</span>
|
||||
</span><span id="L-120"><a href="#L-120"><span class="linenos"> 120</span></a> <span class="sd">"""Apply the gamma method to the content of the Corr."""</span>
|
||||
</span><span id="L-121"><a href="#L-121"><span class="linenos"> 121</span></a> <span class="k">for</span> <span class="n">item</span> <span class="ow">in</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">:</span>
|
||||
</span><span id="L-122"><a href="#L-122"><span class="linenos"> 122</span></a> <span class="k">if</span> <span class="ow">not</span><span class="p">(</span><span class="n">item</span> <span class="ow">is</span> <span class="kc">None</span><span class="p">):</span>
|
||||
</span><span id="L-122"><a href="#L-122"><span class="linenos"> 122</span></a> <span class="k">if</span> <span class="ow">not</span> <span class="p">(</span><span class="n">item</span> <span class="ow">is</span> <span class="kc">None</span><span class="p">):</span>
|
||||
</span><span id="L-123"><a href="#L-123"><span class="linenos"> 123</span></a> <span class="k">if</span> <span class="bp">self</span><span class="o">.</span><span class="n">N</span> <span class="o">==</span> <span class="mi">1</span><span class="p">:</span>
|
||||
</span><span id="L-124"><a href="#L-124"><span class="linenos"> 124</span></a> <span class="n">item</span><span class="p">[</span><span class="mi">0</span><span class="p">]</span><span class="o">.</span><span class="n">gamma_method</span><span class="p">(</span><span class="o">**</span><span class="n">kwargs</span><span class="p">)</span>
|
||||
</span><span id="L-125"><a href="#L-125"><span class="linenos"> 125</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
|
|
@ -356,7 +356,7 @@
|
|||
</span><span id="L-139"><a href="#L-139"><span class="linenos"> 139</span></a>
|
||||
</span><span id="L-140"><a href="#L-140"><span class="linenos"> 140</span></a> <span class="k">if</span> <span class="n">vector_l</span> <span class="ow">is</span> <span class="kc">None</span><span class="p">:</span>
|
||||
</span><span id="L-141"><a href="#L-141"><span class="linenos"> 141</span></a> <span class="n">vector_l</span><span class="p">,</span> <span class="n">vector_r</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">asarray</span><span class="p">([</span><span class="mf">1.</span><span class="p">]</span> <span class="o">+</span> <span class="p">(</span><span class="bp">self</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="p">[</span><span class="mf">0.</span><span class="p">]),</span> <span class="n">np</span><span class="o">.</span><span class="n">asarray</span><span class="p">([</span><span class="mf">1.</span><span class="p">]</span> <span class="o">+</span> <span class="p">(</span><span class="bp">self</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="p">[</span><span class="mf">0.</span><span class="p">])</span>
|
||||
</span><span id="L-142"><a href="#L-142"><span class="linenos"> 142</span></a> <span class="k">elif</span><span class="p">(</span><span class="n">vector_r</span> <span class="ow">is</span> <span class="kc">None</span><span class="p">):</span>
|
||||
</span><span id="L-142"><a href="#L-142"><span class="linenos"> 142</span></a> <span class="k">elif</span> <span class="p">(</span><span class="n">vector_r</span> <span class="ow">is</span> <span class="kc">None</span><span class="p">):</span>
|
||||
</span><span id="L-143"><a href="#L-143"><span class="linenos"> 143</span></a> <span class="n">vector_r</span> <span class="o">=</span> <span class="n">vector_l</span>
|
||||
</span><span id="L-144"><a href="#L-144"><span class="linenos"> 144</span></a> <span class="k">if</span> <span class="nb">isinstance</span><span class="p">(</span><span class="n">vector_l</span><span class="p">,</span> <span class="nb">list</span><span class="p">)</span> <span class="ow">and</span> <span class="ow">not</span> <span class="nb">isinstance</span><span class="p">(</span><span class="n">vector_r</span><span class="p">,</span> <span class="nb">list</span><span class="p">):</span>
|
||||
</span><span id="L-145"><a href="#L-145"><span class="linenos"> 145</span></a> <span class="k">if</span> <span class="nb">len</span><span class="p">(</span><span class="n">vector_l</span><span class="p">)</span> <span class="o">!=</span> <span class="bp">self</span><span class="o">.</span><span class="n">T</span><span class="p">:</span>
|
||||
|
|
@ -395,7 +395,7 @@
|
|||
</span><span id="L-178"><a href="#L-178"><span class="linenos"> 178</span></a><span class="sd"> """</span>
|
||||
</span><span id="L-179"><a href="#L-179"><span class="linenos"> 179</span></a> <span class="k">if</span> <span class="bp">self</span><span class="o">.</span><span class="n">N</span> <span class="o">==</span> <span class="mi">1</span><span class="p">:</span>
|
||||
</span><span id="L-180"><a href="#L-180"><span class="linenos"> 180</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"Trying to pick item from projected Corr"</span><span class="p">)</span>
|
||||
</span><span id="L-181"><a href="#L-181"><span class="linenos"> 181</span></a> <span class="n">newcontent</span> <span class="o">=</span> <span class="p">[</span><span class="kc">None</span> <span class="k">if</span><span class="p">(</span><span class="n">item</span> <span class="ow">is</span> <span class="kc">None</span><span class="p">)</span> <span class="k">else</span> <span class="n">item</span><span class="p">[</span><span class="n">i</span><span class="p">,</span> <span class="n">j</span><span class="p">]</span> <span class="k">for</span> <span class="n">item</span> <span class="ow">in</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">]</span>
|
||||
</span><span id="L-181"><a href="#L-181"><span class="linenos"> 181</span></a> <span class="n">newcontent</span> <span class="o">=</span> <span class="p">[</span><span class="kc">None</span> <span class="k">if</span> <span class="p">(</span><span class="n">item</span> <span class="ow">is</span> <span class="kc">None</span><span class="p">)</span> <span class="k">else</span> <span class="n">item</span><span class="p">[</span><span class="n">i</span><span class="p">,</span> <span class="n">j</span><span class="p">]</span> <span class="k">for</span> <span class="n">item</span> <span class="ow">in</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">]</span>
|
||||
</span><span id="L-182"><a href="#L-182"><span class="linenos"> 182</span></a> <span class="k">return</span> <span class="n">Corr</span><span class="p">(</span><span class="n">newcontent</span><span class="p">)</span>
|
||||
</span><span id="L-183"><a href="#L-183"><span class="linenos"> 183</span></a>
|
||||
</span><span id="L-184"><a href="#L-184"><span class="linenos"> 184</span></a> <span class="k">def</span> <span class="nf">plottable</span><span class="p">(</span><span class="bp">self</span><span class="p">):</span>
|
||||
|
|
@ -428,7 +428,7 @@
|
|||
</span><span id="L-211"><a href="#L-211"><span class="linenos"> 211</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="kc">None</span><span class="p">)</span>
|
||||
</span><span id="L-212"><a href="#L-212"><span class="linenos"> 212</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="L-213"><a href="#L-213"><span class="linenos"> 213</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="mf">0.5</span> <span class="o">*</span> <span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span><span class="p">]</span> <span class="o">+</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="bp">self</span><span class="o">.</span><span class="n">T</span> <span class="o">-</span> <span class="n">t</span><span class="p">]))</span>
|
||||
</span><span id="L-214"><a href="#L-214"><span class="linenos"> 214</span></a> <span class="k">if</span><span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="L-214"><a href="#L-214"><span class="linenos"> 214</span></a> <span class="k">if</span> <span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="L-215"><a href="#L-215"><span class="linenos"> 215</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"Corr could not be symmetrized: No redundant values"</span><span class="p">)</span>
|
||||
</span><span id="L-216"><a href="#L-216"><span class="linenos"> 216</span></a> <span class="k">return</span> <span class="n">Corr</span><span class="p">(</span><span class="n">newcontent</span><span class="p">,</span> <span class="n">prange</span><span class="o">=</span><span class="bp">self</span><span class="o">.</span><span class="n">prange</span><span class="p">)</span>
|
||||
</span><span id="L-217"><a href="#L-217"><span class="linenos"> 217</span></a>
|
||||
|
|
@ -450,7 +450,7 @@
|
|||
</span><span id="L-233"><a href="#L-233"><span class="linenos"> 233</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="kc">None</span><span class="p">)</span>
|
||||
</span><span id="L-234"><a href="#L-234"><span class="linenos"> 234</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="L-235"><a href="#L-235"><span class="linenos"> 235</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="mf">0.5</span> <span class="o">*</span> <span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span><span class="p">]</span> <span class="o">-</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="bp">self</span><span class="o">.</span><span class="n">T</span> <span class="o">-</span> <span class="n">t</span><span class="p">]))</span>
|
||||
</span><span id="L-236"><a href="#L-236"><span class="linenos"> 236</span></a> <span class="k">if</span><span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="L-236"><a href="#L-236"><span class="linenos"> 236</span></a> <span class="k">if</span> <span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="L-237"><a href="#L-237"><span class="linenos"> 237</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"Corr could not be symmetrized: No redundant values"</span><span class="p">)</span>
|
||||
</span><span id="L-238"><a href="#L-238"><span class="linenos"> 238</span></a> <span class="k">return</span> <span class="n">Corr</span><span class="p">(</span><span class="n">newcontent</span><span class="p">,</span> <span class="n">prange</span><span class="o">=</span><span class="bp">self</span><span class="o">.</span><span class="n">prange</span><span class="p">)</span>
|
||||
</span><span id="L-239"><a href="#L-239"><span class="linenos"> 239</span></a>
|
||||
|
|
@ -751,7 +751,7 @@
|
|||
</span><span id="L-534"><a href="#L-534"><span class="linenos"> 534</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="kc">None</span><span class="p">)</span>
|
||||
</span><span id="L-535"><a href="#L-535"><span class="linenos"> 535</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="L-536"><a href="#L-536"><span class="linenos"> 536</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="mf">0.5</span> <span class="o">*</span> <span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">+</span> <span class="mi">1</span><span class="p">]</span> <span class="o">-</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">-</span> <span class="mi">1</span><span class="p">]))</span>
|
||||
</span><span id="L-537"><a href="#L-537"><span class="linenos"> 537</span></a> <span class="k">if</span><span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="L-537"><a href="#L-537"><span class="linenos"> 537</span></a> <span class="k">if</span> <span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="L-538"><a href="#L-538"><span class="linenos"> 538</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s1">'Derivative is undefined at all timeslices'</span><span class="p">)</span>
|
||||
</span><span id="L-539"><a href="#L-539"><span class="linenos"> 539</span></a> <span class="k">return</span> <span class="n">Corr</span><span class="p">(</span><span class="n">newcontent</span><span class="p">,</span> <span class="n">padding</span><span class="o">=</span><span class="p">[</span><span class="mi">1</span><span class="p">,</span> <span class="mi">1</span><span class="p">])</span>
|
||||
</span><span id="L-540"><a href="#L-540"><span class="linenos"> 540</span></a> <span class="k">elif</span> <span class="n">variant</span> <span class="o">==</span> <span class="s2">"forward"</span><span class="p">:</span>
|
||||
|
|
@ -761,7 +761,7 @@
|
|||
</span><span id="L-544"><a href="#L-544"><span class="linenos"> 544</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="kc">None</span><span class="p">)</span>
|
||||
</span><span id="L-545"><a href="#L-545"><span class="linenos"> 545</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="L-546"><a href="#L-546"><span class="linenos"> 546</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">+</span> <span class="mi">1</span><span class="p">]</span> <span class="o">-</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span><span class="p">])</span>
|
||||
</span><span id="L-547"><a href="#L-547"><span class="linenos"> 547</span></a> <span class="k">if</span><span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="L-547"><a href="#L-547"><span class="linenos"> 547</span></a> <span class="k">if</span> <span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="L-548"><a href="#L-548"><span class="linenos"> 548</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"Derivative is undefined at all timeslices"</span><span class="p">)</span>
|
||||
</span><span id="L-549"><a href="#L-549"><span class="linenos"> 549</span></a> <span class="k">return</span> <span class="n">Corr</span><span class="p">(</span><span class="n">newcontent</span><span class="p">,</span> <span class="n">padding</span><span class="o">=</span><span class="p">[</span><span class="mi">0</span><span class="p">,</span> <span class="mi">1</span><span class="p">])</span>
|
||||
</span><span id="L-550"><a href="#L-550"><span class="linenos"> 550</span></a> <span class="k">elif</span> <span class="n">variant</span> <span class="o">==</span> <span class="s2">"backward"</span><span class="p">:</span>
|
||||
|
|
@ -771,7 +771,7 @@
|
|||
</span><span id="L-554"><a href="#L-554"><span class="linenos"> 554</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="kc">None</span><span class="p">)</span>
|
||||
</span><span id="L-555"><a href="#L-555"><span class="linenos"> 555</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="L-556"><a href="#L-556"><span class="linenos"> 556</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span><span class="p">]</span> <span class="o">-</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">-</span> <span class="mi">1</span><span class="p">])</span>
|
||||
</span><span id="L-557"><a href="#L-557"><span class="linenos"> 557</span></a> <span class="k">if</span><span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="L-557"><a href="#L-557"><span class="linenos"> 557</span></a> <span class="k">if</span> <span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="L-558"><a href="#L-558"><span class="linenos"> 558</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"Derivative is undefined at all timeslices"</span><span class="p">)</span>
|
||||
</span><span id="L-559"><a href="#L-559"><span class="linenos"> 559</span></a> <span class="k">return</span> <span class="n">Corr</span><span class="p">(</span><span class="n">newcontent</span><span class="p">,</span> <span class="n">padding</span><span class="o">=</span><span class="p">[</span><span class="mi">1</span><span class="p">,</span> <span class="mi">0</span><span class="p">])</span>
|
||||
</span><span id="L-560"><a href="#L-560"><span class="linenos"> 560</span></a> <span class="k">elif</span> <span class="n">variant</span> <span class="o">==</span> <span class="s2">"improved"</span><span class="p">:</span>
|
||||
|
|
@ -781,7 +781,7 @@
|
|||
</span><span id="L-564"><a href="#L-564"><span class="linenos"> 564</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="kc">None</span><span class="p">)</span>
|
||||
</span><span id="L-565"><a href="#L-565"><span class="linenos"> 565</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="L-566"><a href="#L-566"><span class="linenos"> 566</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">((</span><span class="mi">1</span> <span class="o">/</span> <span class="mi">12</span><span class="p">)</span> <span class="o">*</span> <span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">-</span> <span class="mi">2</span><span class="p">]</span> <span class="o">-</span> <span class="mi">8</span> <span class="o">*</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">-</span> <span class="mi">1</span><span class="p">]</span> <span class="o">+</span> <span class="mi">8</span> <span class="o">*</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">+</span> <span class="mi">1</span><span class="p">]</span> <span class="o">-</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">+</span> <span class="mi">2</span><span class="p">]))</span>
|
||||
</span><span id="L-567"><a href="#L-567"><span class="linenos"> 567</span></a> <span class="k">if</span><span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="L-567"><a href="#L-567"><span class="linenos"> 567</span></a> <span class="k">if</span> <span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="L-568"><a href="#L-568"><span class="linenos"> 568</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s1">'Derivative is undefined at all timeslices'</span><span class="p">)</span>
|
||||
</span><span id="L-569"><a href="#L-569"><span class="linenos"> 569</span></a> <span class="k">return</span> <span class="n">Corr</span><span class="p">(</span><span class="n">newcontent</span><span class="p">,</span> <span class="n">padding</span><span class="o">=</span><span class="p">[</span><span class="mi">2</span><span class="p">,</span> <span class="mi">2</span><span class="p">])</span>
|
||||
</span><span id="L-570"><a href="#L-570"><span class="linenos"> 570</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
|
|
@ -805,7 +805,7 @@
|
|||
</span><span id="L-588"><a href="#L-588"><span class="linenos"> 588</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="kc">None</span><span class="p">)</span>
|
||||
</span><span id="L-589"><a href="#L-589"><span class="linenos"> 589</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="L-590"><a href="#L-590"><span class="linenos"> 590</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">((</span><span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</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="o">*</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span><span class="p">]</span> <span class="o">+</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">-</span> <span class="mi">1</span><span class="p">]))</span>
|
||||
</span><span id="L-591"><a href="#L-591"><span class="linenos"> 591</span></a> <span class="k">if</span><span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="L-591"><a href="#L-591"><span class="linenos"> 591</span></a> <span class="k">if</span> <span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="L-592"><a href="#L-592"><span class="linenos"> 592</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"Derivative is undefined at all timeslices"</span><span class="p">)</span>
|
||||
</span><span id="L-593"><a href="#L-593"><span class="linenos"> 593</span></a> <span class="k">return</span> <span class="n">Corr</span><span class="p">(</span><span class="n">newcontent</span><span class="p">,</span> <span class="n">padding</span><span class="o">=</span><span class="p">[</span><span class="mi">1</span><span class="p">,</span> <span class="mi">1</span><span class="p">])</span>
|
||||
</span><span id="L-594"><a href="#L-594"><span class="linenos"> 594</span></a> <span class="k">elif</span> <span class="n">variant</span> <span class="o">==</span> <span class="s2">"improved"</span><span class="p">:</span>
|
||||
|
|
@ -815,7 +815,7 @@
|
|||
</span><span id="L-598"><a href="#L-598"><span class="linenos"> 598</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="kc">None</span><span class="p">)</span>
|
||||
</span><span id="L-599"><a href="#L-599"><span class="linenos"> 599</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="L-600"><a href="#L-600"><span class="linenos"> 600</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">((</span><span class="mi">1</span> <span class="o">/</span> <span class="mi">12</span><span class="p">)</span> <span class="o">*</span> <span class="p">(</span><span class="o">-</span><span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">+</span> <span class="mi">2</span><span class="p">]</span> <span class="o">+</span> <span class="mi">16</span> <span class="o">*</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">+</span> <span class="mi">1</span><span class="p">]</span> <span class="o">-</span> <span class="mi">30</span> <span class="o">*</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span><span class="p">]</span> <span class="o">+</span> <span class="mi">16</span> <span class="o">*</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">-</span> <span class="mi">1</span><span class="p">]</span> <span class="o">-</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">-</span> <span class="mi">2</span><span class="p">]))</span>
|
||||
</span><span id="L-601"><a href="#L-601"><span class="linenos"> 601</span></a> <span class="k">if</span><span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="L-601"><a href="#L-601"><span class="linenos"> 601</span></a> <span class="k">if</span> <span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="L-602"><a href="#L-602"><span class="linenos"> 602</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"Derivative is undefined at all timeslices"</span><span class="p">)</span>
|
||||
</span><span id="L-603"><a href="#L-603"><span class="linenos"> 603</span></a> <span class="k">return</span> <span class="n">Corr</span><span class="p">(</span><span class="n">newcontent</span><span class="p">,</span> <span class="n">padding</span><span class="o">=</span><span class="p">[</span><span class="mi">2</span><span class="p">,</span> <span class="mi">2</span><span class="p">])</span>
|
||||
</span><span id="L-604"><a href="#L-604"><span class="linenos"> 604</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
|
|
@ -846,7 +846,7 @@
|
|||
</span><span id="L-629"><a href="#L-629"><span class="linenos"> 629</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="kc">None</span><span class="p">)</span>
|
||||
</span><span id="L-630"><a href="#L-630"><span class="linenos"> 630</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="L-631"><a href="#L-631"><span class="linenos"> 631</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span><span class="p">]</span> <span class="o">/</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">+</span> <span class="mi">1</span><span class="p">])</span>
|
||||
</span><span id="L-632"><a href="#L-632"><span class="linenos"> 632</span></a> <span class="k">if</span><span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="L-632"><a href="#L-632"><span class="linenos"> 632</span></a> <span class="k">if</span> <span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="L-633"><a href="#L-633"><span class="linenos"> 633</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s1">'m_eff is undefined at all timeslices'</span><span class="p">)</span>
|
||||
</span><span id="L-634"><a href="#L-634"><span class="linenos"> 634</span></a>
|
||||
</span><span id="L-635"><a href="#L-635"><span class="linenos"> 635</span></a> <span class="k">return</span> <span class="n">np</span><span class="o">.</span><span class="n">log</span><span class="p">(</span><span class="n">Corr</span><span class="p">(</span><span class="n">newcontent</span><span class="p">,</span> <span class="n">padding</span><span class="o">=</span><span class="p">[</span><span class="mi">0</span><span class="p">,</span> <span class="mi">1</span><span class="p">]))</span>
|
||||
|
|
@ -871,7 +871,7 @@
|
|||
</span><span id="L-654"><a href="#L-654"><span class="linenos"> 654</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="kc">None</span><span class="p">)</span>
|
||||
</span><span id="L-655"><a href="#L-655"><span class="linenos"> 655</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="L-656"><a href="#L-656"><span class="linenos"> 656</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">np</span><span class="o">.</span><span class="n">abs</span><span class="p">(</span><span class="n">find_root</span><span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span><span class="p">][</span><span class="mi">0</span><span class="p">]</span> <span class="o">/</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">+</span> <span class="mi">1</span><span class="p">][</span><span class="mi">0</span><span class="p">],</span> <span class="n">root_function</span><span class="p">,</span> <span class="n">guess</span><span class="o">=</span><span class="n">guess</span><span class="p">)))</span>
|
||||
</span><span id="L-657"><a href="#L-657"><span class="linenos"> 657</span></a> <span class="k">if</span><span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="L-657"><a href="#L-657"><span class="linenos"> 657</span></a> <span class="k">if</span> <span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="L-658"><a href="#L-658"><span class="linenos"> 658</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s1">'m_eff is undefined at all timeslices'</span><span class="p">)</span>
|
||||
</span><span id="L-659"><a href="#L-659"><span class="linenos"> 659</span></a>
|
||||
</span><span id="L-660"><a href="#L-660"><span class="linenos"> 660</span></a> <span class="k">return</span> <span class="n">Corr</span><span class="p">(</span><span class="n">newcontent</span><span class="p">,</span> <span class="n">padding</span><span class="o">=</span><span class="p">[</span><span class="mi">0</span><span class="p">,</span> <span class="mi">1</span><span class="p">])</span>
|
||||
|
|
@ -883,7 +883,7 @@
|
|||
</span><span id="L-666"><a href="#L-666"><span class="linenos"> 666</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="kc">None</span><span class="p">)</span>
|
||||
</span><span id="L-667"><a href="#L-667"><span class="linenos"> 667</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="L-668"><a href="#L-668"><span class="linenos"> 668</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">((</span><span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">+</span> <span class="mi">1</span><span class="p">]</span> <span class="o">+</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">-</span> <span class="mi">1</span><span class="p">])</span> <span class="o">/</span> <span class="p">(</span><span class="mi">2</span> <span class="o">*</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span><span class="p">]))</span>
|
||||
</span><span id="L-669"><a href="#L-669"><span class="linenos"> 669</span></a> <span class="k">if</span><span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="L-669"><a href="#L-669"><span class="linenos"> 669</span></a> <span class="k">if</span> <span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="L-670"><a href="#L-670"><span class="linenos"> 670</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"m_eff is undefined at all timeslices"</span><span class="p">)</span>
|
||||
</span><span id="L-671"><a href="#L-671"><span class="linenos"> 671</span></a> <span class="k">return</span> <span class="n">np</span><span class="o">.</span><span class="n">arccosh</span><span class="p">(</span><span class="n">Corr</span><span class="p">(</span><span class="n">newcontent</span><span class="p">,</span> <span class="n">padding</span><span class="o">=</span><span class="p">[</span><span class="mi">1</span><span class="p">,</span> <span class="mi">1</span><span class="p">]))</span>
|
||||
</span><span id="L-672"><a href="#L-672"><span class="linenos"> 672</span></a>
|
||||
|
|
@ -946,7 +946,7 @@
|
|||
</span><span id="L-729"><a href="#L-729"><span class="linenos"> 729</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"no plateau range provided"</span><span class="p">)</span>
|
||||
</span><span id="L-730"><a href="#L-730"><span class="linenos"> 730</span></a> <span class="k">if</span> <span class="bp">self</span><span class="o">.</span><span class="n">N</span> <span class="o">!=</span> <span class="mi">1</span><span class="p">:</span>
|
||||
</span><span id="L-731"><a href="#L-731"><span class="linenos"> 731</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"Correlator must be projected before getting a plateau."</span><span class="p">)</span>
|
||||
</span><span id="L-732"><a href="#L-732"><span class="linenos"> 732</span></a> <span class="k">if</span><span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span><span class="p">]</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">t</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="n">plateau_range</span><span class="p">[</span><span class="mi">0</span><span class="p">],</span> <span class="n">plateau_range</span><span class="p">[</span><span class="mi">1</span><span class="p">]</span> <span class="o">+</span> <span class="mi">1</span><span class="p">)])):</span>
|
||||
</span><span id="L-732"><a href="#L-732"><span class="linenos"> 732</span></a> <span class="k">if</span> <span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span><span class="p">]</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">t</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="n">plateau_range</span><span class="p">[</span><span class="mi">0</span><span class="p">],</span> <span class="n">plateau_range</span><span class="p">[</span><span class="mi">1</span><span class="p">]</span> <span class="o">+</span> <span class="mi">1</span><span class="p">)])):</span>
|
||||
</span><span id="L-733"><a href="#L-733"><span class="linenos"> 733</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"plateau is undefined at all timeslices in plateaurange."</span><span class="p">)</span>
|
||||
</span><span id="L-734"><a href="#L-734"><span class="linenos"> 734</span></a> <span class="k">if</span> <span class="n">auto_gamma</span><span class="p">:</span>
|
||||
</span><span id="L-735"><a href="#L-735"><span class="linenos"> 735</span></a> <span class="bp">self</span><span class="o">.</span><span class="n">gamma_method</span><span class="p">()</span>
|
||||
|
|
@ -1211,7 +1211,7 @@
|
|||
</span><span id="L-994"><a href="#L-994"><span class="linenos"> 994</span></a>
|
||||
</span><span id="L-995"><a href="#L-995"><span class="linenos"> 995</span></a> <span class="k">def</span> <span class="fm">__mul__</span><span class="p">(</span><span class="bp">self</span><span class="p">,</span> <span class="n">y</span><span class="p">):</span>
|
||||
</span><span id="L-996"><a href="#L-996"><span class="linenos"> 996</span></a> <span class="k">if</span> <span class="nb">isinstance</span><span class="p">(</span><span class="n">y</span><span class="p">,</span> <span class="n">Corr</span><span class="p">):</span>
|
||||
</span><span id="L-997"><a href="#L-997"><span class="linenos"> 997</span></a> <span class="k">if</span> <span class="ow">not</span><span class="p">((</span><span class="bp">self</span><span class="o">.</span><span class="n">N</span> <span class="o">==</span> <span class="mi">1</span> <span class="ow">or</span> <span class="n">y</span><span class="o">.</span><span class="n">N</span> <span class="o">==</span> <span class="mi">1</span> <span class="ow">or</span> <span class="bp">self</span><span class="o">.</span><span class="n">N</span> <span class="o">==</span> <span class="n">y</span><span class="o">.</span><span class="n">N</span><span class="p">)</span> <span class="ow">and</span> <span class="bp">self</span><span class="o">.</span><span class="n">T</span> <span class="o">==</span> <span class="n">y</span><span class="o">.</span><span class="n">T</span><span class="p">):</span>
|
||||
</span><span id="L-997"><a href="#L-997"><span class="linenos"> 997</span></a> <span class="k">if</span> <span class="ow">not</span> <span class="p">((</span><span class="bp">self</span><span class="o">.</span><span class="n">N</span> <span class="o">==</span> <span class="mi">1</span> <span class="ow">or</span> <span class="n">y</span><span class="o">.</span><span class="n">N</span> <span class="o">==</span> <span class="mi">1</span> <span class="ow">or</span> <span class="bp">self</span><span class="o">.</span><span class="n">N</span> <span class="o">==</span> <span class="n">y</span><span class="o">.</span><span class="n">N</span><span class="p">)</span> <span class="ow">and</span> <span class="bp">self</span><span class="o">.</span><span class="n">T</span> <span class="o">==</span> <span class="n">y</span><span class="o">.</span><span class="n">T</span><span class="p">):</span>
|
||||
</span><span id="L-998"><a href="#L-998"><span class="linenos"> 998</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"Multiplication of Corr object requires N=N or N=1 and T=T"</span><span class="p">)</span>
|
||||
</span><span id="L-999"><a href="#L-999"><span class="linenos"> 999</span></a> <span class="n">newcontent</span> <span class="o">=</span> <span class="p">[]</span>
|
||||
</span><span id="L-1000"><a href="#L-1000"><span class="linenos">1000</span></a> <span class="k">for</span> <span class="n">t</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">T</span><span class="p">):</span>
|
||||
|
|
@ -1239,7 +1239,7 @@
|
|||
</span><span id="L-1022"><a href="#L-1022"><span class="linenos">1022</span></a>
|
||||
</span><span id="L-1023"><a href="#L-1023"><span class="linenos">1023</span></a> <span class="k">def</span> <span class="fm">__truediv__</span><span class="p">(</span><span class="bp">self</span><span class="p">,</span> <span class="n">y</span><span class="p">):</span>
|
||||
</span><span id="L-1024"><a href="#L-1024"><span class="linenos">1024</span></a> <span class="k">if</span> <span class="nb">isinstance</span><span class="p">(</span><span class="n">y</span><span class="p">,</span> <span class="n">Corr</span><span class="p">):</span>
|
||||
</span><span id="L-1025"><a href="#L-1025"><span class="linenos">1025</span></a> <span class="k">if</span> <span class="ow">not</span><span class="p">((</span><span class="bp">self</span><span class="o">.</span><span class="n">N</span> <span class="o">==</span> <span class="mi">1</span> <span class="ow">or</span> <span class="n">y</span><span class="o">.</span><span class="n">N</span> <span class="o">==</span> <span class="mi">1</span> <span class="ow">or</span> <span class="bp">self</span><span class="o">.</span><span class="n">N</span> <span class="o">==</span> <span class="n">y</span><span class="o">.</span><span class="n">N</span><span class="p">)</span> <span class="ow">and</span> <span class="bp">self</span><span class="o">.</span><span class="n">T</span> <span class="o">==</span> <span class="n">y</span><span class="o">.</span><span class="n">T</span><span class="p">):</span>
|
||||
</span><span id="L-1025"><a href="#L-1025"><span class="linenos">1025</span></a> <span class="k">if</span> <span class="ow">not</span> <span class="p">((</span><span class="bp">self</span><span class="o">.</span><span class="n">N</span> <span class="o">==</span> <span class="mi">1</span> <span class="ow">or</span> <span class="n">y</span><span class="o">.</span><span class="n">N</span> <span class="o">==</span> <span class="mi">1</span> <span class="ow">or</span> <span class="bp">self</span><span class="o">.</span><span class="n">N</span> <span class="o">==</span> <span class="n">y</span><span class="o">.</span><span class="n">N</span><span class="p">)</span> <span class="ow">and</span> <span class="bp">self</span><span class="o">.</span><span class="n">T</span> <span class="o">==</span> <span class="n">y</span><span class="o">.</span><span class="n">T</span><span class="p">):</span>
|
||||
</span><span id="L-1026"><a href="#L-1026"><span class="linenos">1026</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"Multiplication of Corr object requires N=N or N=1 and T=T"</span><span class="p">)</span>
|
||||
</span><span id="L-1027"><a href="#L-1027"><span class="linenos">1027</span></a> <span class="n">newcontent</span> <span class="o">=</span> <span class="p">[]</span>
|
||||
</span><span id="L-1028"><a href="#L-1028"><span class="linenos">1028</span></a> <span class="k">for</span> <span class="n">t</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">T</span><span class="p">):</span>
|
||||
|
|
@ -1622,7 +1622,7 @@
|
|||
</span><span id="Corr-120"><a href="#Corr-120"><span class="linenos"> 120</span></a> <span class="k">def</span> <span class="nf">gamma_method</span><span class="p">(</span><span class="bp">self</span><span class="p">,</span> <span class="o">**</span><span class="n">kwargs</span><span class="p">):</span>
|
||||
</span><span id="Corr-121"><a href="#Corr-121"><span class="linenos"> 121</span></a> <span class="sd">"""Apply the gamma method to the content of the Corr."""</span>
|
||||
</span><span id="Corr-122"><a href="#Corr-122"><span class="linenos"> 122</span></a> <span class="k">for</span> <span class="n">item</span> <span class="ow">in</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">:</span>
|
||||
</span><span id="Corr-123"><a href="#Corr-123"><span class="linenos"> 123</span></a> <span class="k">if</span> <span class="ow">not</span><span class="p">(</span><span class="n">item</span> <span class="ow">is</span> <span class="kc">None</span><span class="p">):</span>
|
||||
</span><span id="Corr-123"><a href="#Corr-123"><span class="linenos"> 123</span></a> <span class="k">if</span> <span class="ow">not</span> <span class="p">(</span><span class="n">item</span> <span class="ow">is</span> <span class="kc">None</span><span class="p">):</span>
|
||||
</span><span id="Corr-124"><a href="#Corr-124"><span class="linenos"> 124</span></a> <span class="k">if</span> <span class="bp">self</span><span class="o">.</span><span class="n">N</span> <span class="o">==</span> <span class="mi">1</span><span class="p">:</span>
|
||||
</span><span id="Corr-125"><a href="#Corr-125"><span class="linenos"> 125</span></a> <span class="n">item</span><span class="p">[</span><span class="mi">0</span><span class="p">]</span><span class="o">.</span><span class="n">gamma_method</span><span class="p">(</span><span class="o">**</span><span class="n">kwargs</span><span class="p">)</span>
|
||||
</span><span id="Corr-126"><a href="#Corr-126"><span class="linenos"> 126</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
|
|
@ -1642,7 +1642,7 @@
|
|||
</span><span id="Corr-140"><a href="#Corr-140"><span class="linenos"> 140</span></a>
|
||||
</span><span id="Corr-141"><a href="#Corr-141"><span class="linenos"> 141</span></a> <span class="k">if</span> <span class="n">vector_l</span> <span class="ow">is</span> <span class="kc">None</span><span class="p">:</span>
|
||||
</span><span id="Corr-142"><a href="#Corr-142"><span class="linenos"> 142</span></a> <span class="n">vector_l</span><span class="p">,</span> <span class="n">vector_r</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">asarray</span><span class="p">([</span><span class="mf">1.</span><span class="p">]</span> <span class="o">+</span> <span class="p">(</span><span class="bp">self</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="p">[</span><span class="mf">0.</span><span class="p">]),</span> <span class="n">np</span><span class="o">.</span><span class="n">asarray</span><span class="p">([</span><span class="mf">1.</span><span class="p">]</span> <span class="o">+</span> <span class="p">(</span><span class="bp">self</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="p">[</span><span class="mf">0.</span><span class="p">])</span>
|
||||
</span><span id="Corr-143"><a href="#Corr-143"><span class="linenos"> 143</span></a> <span class="k">elif</span><span class="p">(</span><span class="n">vector_r</span> <span class="ow">is</span> <span class="kc">None</span><span class="p">):</span>
|
||||
</span><span id="Corr-143"><a href="#Corr-143"><span class="linenos"> 143</span></a> <span class="k">elif</span> <span class="p">(</span><span class="n">vector_r</span> <span class="ow">is</span> <span class="kc">None</span><span class="p">):</span>
|
||||
</span><span id="Corr-144"><a href="#Corr-144"><span class="linenos"> 144</span></a> <span class="n">vector_r</span> <span class="o">=</span> <span class="n">vector_l</span>
|
||||
</span><span id="Corr-145"><a href="#Corr-145"><span class="linenos"> 145</span></a> <span class="k">if</span> <span class="nb">isinstance</span><span class="p">(</span><span class="n">vector_l</span><span class="p">,</span> <span class="nb">list</span><span class="p">)</span> <span class="ow">and</span> <span class="ow">not</span> <span class="nb">isinstance</span><span class="p">(</span><span class="n">vector_r</span><span class="p">,</span> <span class="nb">list</span><span class="p">):</span>
|
||||
</span><span id="Corr-146"><a href="#Corr-146"><span class="linenos"> 146</span></a> <span class="k">if</span> <span class="nb">len</span><span class="p">(</span><span class="n">vector_l</span><span class="p">)</span> <span class="o">!=</span> <span class="bp">self</span><span class="o">.</span><span class="n">T</span><span class="p">:</span>
|
||||
|
|
@ -1681,7 +1681,7 @@
|
|||
</span><span id="Corr-179"><a href="#Corr-179"><span class="linenos"> 179</span></a><span class="sd"> """</span>
|
||||
</span><span id="Corr-180"><a href="#Corr-180"><span class="linenos"> 180</span></a> <span class="k">if</span> <span class="bp">self</span><span class="o">.</span><span class="n">N</span> <span class="o">==</span> <span class="mi">1</span><span class="p">:</span>
|
||||
</span><span id="Corr-181"><a href="#Corr-181"><span class="linenos"> 181</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"Trying to pick item from projected Corr"</span><span class="p">)</span>
|
||||
</span><span id="Corr-182"><a href="#Corr-182"><span class="linenos"> 182</span></a> <span class="n">newcontent</span> <span class="o">=</span> <span class="p">[</span><span class="kc">None</span> <span class="k">if</span><span class="p">(</span><span class="n">item</span> <span class="ow">is</span> <span class="kc">None</span><span class="p">)</span> <span class="k">else</span> <span class="n">item</span><span class="p">[</span><span class="n">i</span><span class="p">,</span> <span class="n">j</span><span class="p">]</span> <span class="k">for</span> <span class="n">item</span> <span class="ow">in</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">]</span>
|
||||
</span><span id="Corr-182"><a href="#Corr-182"><span class="linenos"> 182</span></a> <span class="n">newcontent</span> <span class="o">=</span> <span class="p">[</span><span class="kc">None</span> <span class="k">if</span> <span class="p">(</span><span class="n">item</span> <span class="ow">is</span> <span class="kc">None</span><span class="p">)</span> <span class="k">else</span> <span class="n">item</span><span class="p">[</span><span class="n">i</span><span class="p">,</span> <span class="n">j</span><span class="p">]</span> <span class="k">for</span> <span class="n">item</span> <span class="ow">in</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">]</span>
|
||||
</span><span id="Corr-183"><a href="#Corr-183"><span class="linenos"> 183</span></a> <span class="k">return</span> <span class="n">Corr</span><span class="p">(</span><span class="n">newcontent</span><span class="p">)</span>
|
||||
</span><span id="Corr-184"><a href="#Corr-184"><span class="linenos"> 184</span></a>
|
||||
</span><span id="Corr-185"><a href="#Corr-185"><span class="linenos"> 185</span></a> <span class="k">def</span> <span class="nf">plottable</span><span class="p">(</span><span class="bp">self</span><span class="p">):</span>
|
||||
|
|
@ -1714,7 +1714,7 @@
|
|||
</span><span id="Corr-212"><a href="#Corr-212"><span class="linenos"> 212</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="kc">None</span><span class="p">)</span>
|
||||
</span><span id="Corr-213"><a href="#Corr-213"><span class="linenos"> 213</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="Corr-214"><a href="#Corr-214"><span class="linenos"> 214</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="mf">0.5</span> <span class="o">*</span> <span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span><span class="p">]</span> <span class="o">+</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="bp">self</span><span class="o">.</span><span class="n">T</span> <span class="o">-</span> <span class="n">t</span><span class="p">]))</span>
|
||||
</span><span id="Corr-215"><a href="#Corr-215"><span class="linenos"> 215</span></a> <span class="k">if</span><span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="Corr-215"><a href="#Corr-215"><span class="linenos"> 215</span></a> <span class="k">if</span> <span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="Corr-216"><a href="#Corr-216"><span class="linenos"> 216</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"Corr could not be symmetrized: No redundant values"</span><span class="p">)</span>
|
||||
</span><span id="Corr-217"><a href="#Corr-217"><span class="linenos"> 217</span></a> <span class="k">return</span> <span class="n">Corr</span><span class="p">(</span><span class="n">newcontent</span><span class="p">,</span> <span class="n">prange</span><span class="o">=</span><span class="bp">self</span><span class="o">.</span><span class="n">prange</span><span class="p">)</span>
|
||||
</span><span id="Corr-218"><a href="#Corr-218"><span class="linenos"> 218</span></a>
|
||||
|
|
@ -1736,7 +1736,7 @@
|
|||
</span><span id="Corr-234"><a href="#Corr-234"><span class="linenos"> 234</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="kc">None</span><span class="p">)</span>
|
||||
</span><span id="Corr-235"><a href="#Corr-235"><span class="linenos"> 235</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="Corr-236"><a href="#Corr-236"><span class="linenos"> 236</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="mf">0.5</span> <span class="o">*</span> <span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span><span class="p">]</span> <span class="o">-</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="bp">self</span><span class="o">.</span><span class="n">T</span> <span class="o">-</span> <span class="n">t</span><span class="p">]))</span>
|
||||
</span><span id="Corr-237"><a href="#Corr-237"><span class="linenos"> 237</span></a> <span class="k">if</span><span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="Corr-237"><a href="#Corr-237"><span class="linenos"> 237</span></a> <span class="k">if</span> <span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="Corr-238"><a href="#Corr-238"><span class="linenos"> 238</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"Corr could not be symmetrized: No redundant values"</span><span class="p">)</span>
|
||||
</span><span id="Corr-239"><a href="#Corr-239"><span class="linenos"> 239</span></a> <span class="k">return</span> <span class="n">Corr</span><span class="p">(</span><span class="n">newcontent</span><span class="p">,</span> <span class="n">prange</span><span class="o">=</span><span class="bp">self</span><span class="o">.</span><span class="n">prange</span><span class="p">)</span>
|
||||
</span><span id="Corr-240"><a href="#Corr-240"><span class="linenos"> 240</span></a>
|
||||
|
|
@ -2037,7 +2037,7 @@
|
|||
</span><span id="Corr-535"><a href="#Corr-535"><span class="linenos"> 535</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="kc">None</span><span class="p">)</span>
|
||||
</span><span id="Corr-536"><a href="#Corr-536"><span class="linenos"> 536</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="Corr-537"><a href="#Corr-537"><span class="linenos"> 537</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="mf">0.5</span> <span class="o">*</span> <span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">+</span> <span class="mi">1</span><span class="p">]</span> <span class="o">-</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">-</span> <span class="mi">1</span><span class="p">]))</span>
|
||||
</span><span id="Corr-538"><a href="#Corr-538"><span class="linenos"> 538</span></a> <span class="k">if</span><span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="Corr-538"><a href="#Corr-538"><span class="linenos"> 538</span></a> <span class="k">if</span> <span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="Corr-539"><a href="#Corr-539"><span class="linenos"> 539</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s1">'Derivative is undefined at all timeslices'</span><span class="p">)</span>
|
||||
</span><span id="Corr-540"><a href="#Corr-540"><span class="linenos"> 540</span></a> <span class="k">return</span> <span class="n">Corr</span><span class="p">(</span><span class="n">newcontent</span><span class="p">,</span> <span class="n">padding</span><span class="o">=</span><span class="p">[</span><span class="mi">1</span><span class="p">,</span> <span class="mi">1</span><span class="p">])</span>
|
||||
</span><span id="Corr-541"><a href="#Corr-541"><span class="linenos"> 541</span></a> <span class="k">elif</span> <span class="n">variant</span> <span class="o">==</span> <span class="s2">"forward"</span><span class="p">:</span>
|
||||
|
|
@ -2047,7 +2047,7 @@
|
|||
</span><span id="Corr-545"><a href="#Corr-545"><span class="linenos"> 545</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="kc">None</span><span class="p">)</span>
|
||||
</span><span id="Corr-546"><a href="#Corr-546"><span class="linenos"> 546</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="Corr-547"><a href="#Corr-547"><span class="linenos"> 547</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">+</span> <span class="mi">1</span><span class="p">]</span> <span class="o">-</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span><span class="p">])</span>
|
||||
</span><span id="Corr-548"><a href="#Corr-548"><span class="linenos"> 548</span></a> <span class="k">if</span><span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="Corr-548"><a href="#Corr-548"><span class="linenos"> 548</span></a> <span class="k">if</span> <span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="Corr-549"><a href="#Corr-549"><span class="linenos"> 549</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"Derivative is undefined at all timeslices"</span><span class="p">)</span>
|
||||
</span><span id="Corr-550"><a href="#Corr-550"><span class="linenos"> 550</span></a> <span class="k">return</span> <span class="n">Corr</span><span class="p">(</span><span class="n">newcontent</span><span class="p">,</span> <span class="n">padding</span><span class="o">=</span><span class="p">[</span><span class="mi">0</span><span class="p">,</span> <span class="mi">1</span><span class="p">])</span>
|
||||
</span><span id="Corr-551"><a href="#Corr-551"><span class="linenos"> 551</span></a> <span class="k">elif</span> <span class="n">variant</span> <span class="o">==</span> <span class="s2">"backward"</span><span class="p">:</span>
|
||||
|
|
@ -2057,7 +2057,7 @@
|
|||
</span><span id="Corr-555"><a href="#Corr-555"><span class="linenos"> 555</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="kc">None</span><span class="p">)</span>
|
||||
</span><span id="Corr-556"><a href="#Corr-556"><span class="linenos"> 556</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="Corr-557"><a href="#Corr-557"><span class="linenos"> 557</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span><span class="p">]</span> <span class="o">-</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">-</span> <span class="mi">1</span><span class="p">])</span>
|
||||
</span><span id="Corr-558"><a href="#Corr-558"><span class="linenos"> 558</span></a> <span class="k">if</span><span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="Corr-558"><a href="#Corr-558"><span class="linenos"> 558</span></a> <span class="k">if</span> <span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="Corr-559"><a href="#Corr-559"><span class="linenos"> 559</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"Derivative is undefined at all timeslices"</span><span class="p">)</span>
|
||||
</span><span id="Corr-560"><a href="#Corr-560"><span class="linenos"> 560</span></a> <span class="k">return</span> <span class="n">Corr</span><span class="p">(</span><span class="n">newcontent</span><span class="p">,</span> <span class="n">padding</span><span class="o">=</span><span class="p">[</span><span class="mi">1</span><span class="p">,</span> <span class="mi">0</span><span class="p">])</span>
|
||||
</span><span id="Corr-561"><a href="#Corr-561"><span class="linenos"> 561</span></a> <span class="k">elif</span> <span class="n">variant</span> <span class="o">==</span> <span class="s2">"improved"</span><span class="p">:</span>
|
||||
|
|
@ -2067,7 +2067,7 @@
|
|||
</span><span id="Corr-565"><a href="#Corr-565"><span class="linenos"> 565</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="kc">None</span><span class="p">)</span>
|
||||
</span><span id="Corr-566"><a href="#Corr-566"><span class="linenos"> 566</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="Corr-567"><a href="#Corr-567"><span class="linenos"> 567</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">((</span><span class="mi">1</span> <span class="o">/</span> <span class="mi">12</span><span class="p">)</span> <span class="o">*</span> <span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">-</span> <span class="mi">2</span><span class="p">]</span> <span class="o">-</span> <span class="mi">8</span> <span class="o">*</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">-</span> <span class="mi">1</span><span class="p">]</span> <span class="o">+</span> <span class="mi">8</span> <span class="o">*</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">+</span> <span class="mi">1</span><span class="p">]</span> <span class="o">-</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">+</span> <span class="mi">2</span><span class="p">]))</span>
|
||||
</span><span id="Corr-568"><a href="#Corr-568"><span class="linenos"> 568</span></a> <span class="k">if</span><span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="Corr-568"><a href="#Corr-568"><span class="linenos"> 568</span></a> <span class="k">if</span> <span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="Corr-569"><a href="#Corr-569"><span class="linenos"> 569</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s1">'Derivative is undefined at all timeslices'</span><span class="p">)</span>
|
||||
</span><span id="Corr-570"><a href="#Corr-570"><span class="linenos"> 570</span></a> <span class="k">return</span> <span class="n">Corr</span><span class="p">(</span><span class="n">newcontent</span><span class="p">,</span> <span class="n">padding</span><span class="o">=</span><span class="p">[</span><span class="mi">2</span><span class="p">,</span> <span class="mi">2</span><span class="p">])</span>
|
||||
</span><span id="Corr-571"><a href="#Corr-571"><span class="linenos"> 571</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
|
|
@ -2091,7 +2091,7 @@
|
|||
</span><span id="Corr-589"><a href="#Corr-589"><span class="linenos"> 589</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="kc">None</span><span class="p">)</span>
|
||||
</span><span id="Corr-590"><a href="#Corr-590"><span class="linenos"> 590</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="Corr-591"><a href="#Corr-591"><span class="linenos"> 591</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">((</span><span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</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="o">*</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span><span class="p">]</span> <span class="o">+</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">-</span> <span class="mi">1</span><span class="p">]))</span>
|
||||
</span><span id="Corr-592"><a href="#Corr-592"><span class="linenos"> 592</span></a> <span class="k">if</span><span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="Corr-592"><a href="#Corr-592"><span class="linenos"> 592</span></a> <span class="k">if</span> <span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="Corr-593"><a href="#Corr-593"><span class="linenos"> 593</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"Derivative is undefined at all timeslices"</span><span class="p">)</span>
|
||||
</span><span id="Corr-594"><a href="#Corr-594"><span class="linenos"> 594</span></a> <span class="k">return</span> <span class="n">Corr</span><span class="p">(</span><span class="n">newcontent</span><span class="p">,</span> <span class="n">padding</span><span class="o">=</span><span class="p">[</span><span class="mi">1</span><span class="p">,</span> <span class="mi">1</span><span class="p">])</span>
|
||||
</span><span id="Corr-595"><a href="#Corr-595"><span class="linenos"> 595</span></a> <span class="k">elif</span> <span class="n">variant</span> <span class="o">==</span> <span class="s2">"improved"</span><span class="p">:</span>
|
||||
|
|
@ -2101,7 +2101,7 @@
|
|||
</span><span id="Corr-599"><a href="#Corr-599"><span class="linenos"> 599</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="kc">None</span><span class="p">)</span>
|
||||
</span><span id="Corr-600"><a href="#Corr-600"><span class="linenos"> 600</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="Corr-601"><a href="#Corr-601"><span class="linenos"> 601</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">((</span><span class="mi">1</span> <span class="o">/</span> <span class="mi">12</span><span class="p">)</span> <span class="o">*</span> <span class="p">(</span><span class="o">-</span><span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">+</span> <span class="mi">2</span><span class="p">]</span> <span class="o">+</span> <span class="mi">16</span> <span class="o">*</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">+</span> <span class="mi">1</span><span class="p">]</span> <span class="o">-</span> <span class="mi">30</span> <span class="o">*</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span><span class="p">]</span> <span class="o">+</span> <span class="mi">16</span> <span class="o">*</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">-</span> <span class="mi">1</span><span class="p">]</span> <span class="o">-</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">-</span> <span class="mi">2</span><span class="p">]))</span>
|
||||
</span><span id="Corr-602"><a href="#Corr-602"><span class="linenos"> 602</span></a> <span class="k">if</span><span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="Corr-602"><a href="#Corr-602"><span class="linenos"> 602</span></a> <span class="k">if</span> <span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="Corr-603"><a href="#Corr-603"><span class="linenos"> 603</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"Derivative is undefined at all timeslices"</span><span class="p">)</span>
|
||||
</span><span id="Corr-604"><a href="#Corr-604"><span class="linenos"> 604</span></a> <span class="k">return</span> <span class="n">Corr</span><span class="p">(</span><span class="n">newcontent</span><span class="p">,</span> <span class="n">padding</span><span class="o">=</span><span class="p">[</span><span class="mi">2</span><span class="p">,</span> <span class="mi">2</span><span class="p">])</span>
|
||||
</span><span id="Corr-605"><a href="#Corr-605"><span class="linenos"> 605</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
|
|
@ -2132,7 +2132,7 @@
|
|||
</span><span id="Corr-630"><a href="#Corr-630"><span class="linenos"> 630</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="kc">None</span><span class="p">)</span>
|
||||
</span><span id="Corr-631"><a href="#Corr-631"><span class="linenos"> 631</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="Corr-632"><a href="#Corr-632"><span class="linenos"> 632</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span><span class="p">]</span> <span class="o">/</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">+</span> <span class="mi">1</span><span class="p">])</span>
|
||||
</span><span id="Corr-633"><a href="#Corr-633"><span class="linenos"> 633</span></a> <span class="k">if</span><span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="Corr-633"><a href="#Corr-633"><span class="linenos"> 633</span></a> <span class="k">if</span> <span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="Corr-634"><a href="#Corr-634"><span class="linenos"> 634</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s1">'m_eff is undefined at all timeslices'</span><span class="p">)</span>
|
||||
</span><span id="Corr-635"><a href="#Corr-635"><span class="linenos"> 635</span></a>
|
||||
</span><span id="Corr-636"><a href="#Corr-636"><span class="linenos"> 636</span></a> <span class="k">return</span> <span class="n">np</span><span class="o">.</span><span class="n">log</span><span class="p">(</span><span class="n">Corr</span><span class="p">(</span><span class="n">newcontent</span><span class="p">,</span> <span class="n">padding</span><span class="o">=</span><span class="p">[</span><span class="mi">0</span><span class="p">,</span> <span class="mi">1</span><span class="p">]))</span>
|
||||
|
|
@ -2157,7 +2157,7 @@
|
|||
</span><span id="Corr-655"><a href="#Corr-655"><span class="linenos"> 655</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="kc">None</span><span class="p">)</span>
|
||||
</span><span id="Corr-656"><a href="#Corr-656"><span class="linenos"> 656</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="Corr-657"><a href="#Corr-657"><span class="linenos"> 657</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">np</span><span class="o">.</span><span class="n">abs</span><span class="p">(</span><span class="n">find_root</span><span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span><span class="p">][</span><span class="mi">0</span><span class="p">]</span> <span class="o">/</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">+</span> <span class="mi">1</span><span class="p">][</span><span class="mi">0</span><span class="p">],</span> <span class="n">root_function</span><span class="p">,</span> <span class="n">guess</span><span class="o">=</span><span class="n">guess</span><span class="p">)))</span>
|
||||
</span><span id="Corr-658"><a href="#Corr-658"><span class="linenos"> 658</span></a> <span class="k">if</span><span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="Corr-658"><a href="#Corr-658"><span class="linenos"> 658</span></a> <span class="k">if</span> <span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="Corr-659"><a href="#Corr-659"><span class="linenos"> 659</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s1">'m_eff is undefined at all timeslices'</span><span class="p">)</span>
|
||||
</span><span id="Corr-660"><a href="#Corr-660"><span class="linenos"> 660</span></a>
|
||||
</span><span id="Corr-661"><a href="#Corr-661"><span class="linenos"> 661</span></a> <span class="k">return</span> <span class="n">Corr</span><span class="p">(</span><span class="n">newcontent</span><span class="p">,</span> <span class="n">padding</span><span class="o">=</span><span class="p">[</span><span class="mi">0</span><span class="p">,</span> <span class="mi">1</span><span class="p">])</span>
|
||||
|
|
@ -2169,7 +2169,7 @@
|
|||
</span><span id="Corr-667"><a href="#Corr-667"><span class="linenos"> 667</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="kc">None</span><span class="p">)</span>
|
||||
</span><span id="Corr-668"><a href="#Corr-668"><span class="linenos"> 668</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="Corr-669"><a href="#Corr-669"><span class="linenos"> 669</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">((</span><span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">+</span> <span class="mi">1</span><span class="p">]</span> <span class="o">+</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">-</span> <span class="mi">1</span><span class="p">])</span> <span class="o">/</span> <span class="p">(</span><span class="mi">2</span> <span class="o">*</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span><span class="p">]))</span>
|
||||
</span><span id="Corr-670"><a href="#Corr-670"><span class="linenos"> 670</span></a> <span class="k">if</span><span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="Corr-670"><a href="#Corr-670"><span class="linenos"> 670</span></a> <span class="k">if</span> <span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="Corr-671"><a href="#Corr-671"><span class="linenos"> 671</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"m_eff is undefined at all timeslices"</span><span class="p">)</span>
|
||||
</span><span id="Corr-672"><a href="#Corr-672"><span class="linenos"> 672</span></a> <span class="k">return</span> <span class="n">np</span><span class="o">.</span><span class="n">arccosh</span><span class="p">(</span><span class="n">Corr</span><span class="p">(</span><span class="n">newcontent</span><span class="p">,</span> <span class="n">padding</span><span class="o">=</span><span class="p">[</span><span class="mi">1</span><span class="p">,</span> <span class="mi">1</span><span class="p">]))</span>
|
||||
</span><span id="Corr-673"><a href="#Corr-673"><span class="linenos"> 673</span></a>
|
||||
|
|
@ -2232,7 +2232,7 @@
|
|||
</span><span id="Corr-730"><a href="#Corr-730"><span class="linenos"> 730</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"no plateau range provided"</span><span class="p">)</span>
|
||||
</span><span id="Corr-731"><a href="#Corr-731"><span class="linenos"> 731</span></a> <span class="k">if</span> <span class="bp">self</span><span class="o">.</span><span class="n">N</span> <span class="o">!=</span> <span class="mi">1</span><span class="p">:</span>
|
||||
</span><span id="Corr-732"><a href="#Corr-732"><span class="linenos"> 732</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"Correlator must be projected before getting a plateau."</span><span class="p">)</span>
|
||||
</span><span id="Corr-733"><a href="#Corr-733"><span class="linenos"> 733</span></a> <span class="k">if</span><span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span><span class="p">]</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">t</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="n">plateau_range</span><span class="p">[</span><span class="mi">0</span><span class="p">],</span> <span class="n">plateau_range</span><span class="p">[</span><span class="mi">1</span><span class="p">]</span> <span class="o">+</span> <span class="mi">1</span><span class="p">)])):</span>
|
||||
</span><span id="Corr-733"><a href="#Corr-733"><span class="linenos"> 733</span></a> <span class="k">if</span> <span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span><span class="p">]</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">t</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="n">plateau_range</span><span class="p">[</span><span class="mi">0</span><span class="p">],</span> <span class="n">plateau_range</span><span class="p">[</span><span class="mi">1</span><span class="p">]</span> <span class="o">+</span> <span class="mi">1</span><span class="p">)])):</span>
|
||||
</span><span id="Corr-734"><a href="#Corr-734"><span class="linenos"> 734</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"plateau is undefined at all timeslices in plateaurange."</span><span class="p">)</span>
|
||||
</span><span id="Corr-735"><a href="#Corr-735"><span class="linenos"> 735</span></a> <span class="k">if</span> <span class="n">auto_gamma</span><span class="p">:</span>
|
||||
</span><span id="Corr-736"><a href="#Corr-736"><span class="linenos"> 736</span></a> <span class="bp">self</span><span class="o">.</span><span class="n">gamma_method</span><span class="p">()</span>
|
||||
|
|
@ -2497,7 +2497,7 @@
|
|||
</span><span id="Corr-995"><a href="#Corr-995"><span class="linenos"> 995</span></a>
|
||||
</span><span id="Corr-996"><a href="#Corr-996"><span class="linenos"> 996</span></a> <span class="k">def</span> <span class="fm">__mul__</span><span class="p">(</span><span class="bp">self</span><span class="p">,</span> <span class="n">y</span><span class="p">):</span>
|
||||
</span><span id="Corr-997"><a href="#Corr-997"><span class="linenos"> 997</span></a> <span class="k">if</span> <span class="nb">isinstance</span><span class="p">(</span><span class="n">y</span><span class="p">,</span> <span class="n">Corr</span><span class="p">):</span>
|
||||
</span><span id="Corr-998"><a href="#Corr-998"><span class="linenos"> 998</span></a> <span class="k">if</span> <span class="ow">not</span><span class="p">((</span><span class="bp">self</span><span class="o">.</span><span class="n">N</span> <span class="o">==</span> <span class="mi">1</span> <span class="ow">or</span> <span class="n">y</span><span class="o">.</span><span class="n">N</span> <span class="o">==</span> <span class="mi">1</span> <span class="ow">or</span> <span class="bp">self</span><span class="o">.</span><span class="n">N</span> <span class="o">==</span> <span class="n">y</span><span class="o">.</span><span class="n">N</span><span class="p">)</span> <span class="ow">and</span> <span class="bp">self</span><span class="o">.</span><span class="n">T</span> <span class="o">==</span> <span class="n">y</span><span class="o">.</span><span class="n">T</span><span class="p">):</span>
|
||||
</span><span id="Corr-998"><a href="#Corr-998"><span class="linenos"> 998</span></a> <span class="k">if</span> <span class="ow">not</span> <span class="p">((</span><span class="bp">self</span><span class="o">.</span><span class="n">N</span> <span class="o">==</span> <span class="mi">1</span> <span class="ow">or</span> <span class="n">y</span><span class="o">.</span><span class="n">N</span> <span class="o">==</span> <span class="mi">1</span> <span class="ow">or</span> <span class="bp">self</span><span class="o">.</span><span class="n">N</span> <span class="o">==</span> <span class="n">y</span><span class="o">.</span><span class="n">N</span><span class="p">)</span> <span class="ow">and</span> <span class="bp">self</span><span class="o">.</span><span class="n">T</span> <span class="o">==</span> <span class="n">y</span><span class="o">.</span><span class="n">T</span><span class="p">):</span>
|
||||
</span><span id="Corr-999"><a href="#Corr-999"><span class="linenos"> 999</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"Multiplication of Corr object requires N=N or N=1 and T=T"</span><span class="p">)</span>
|
||||
</span><span id="Corr-1000"><a href="#Corr-1000"><span class="linenos">1000</span></a> <span class="n">newcontent</span> <span class="o">=</span> <span class="p">[]</span>
|
||||
</span><span id="Corr-1001"><a href="#Corr-1001"><span class="linenos">1001</span></a> <span class="k">for</span> <span class="n">t</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">T</span><span class="p">):</span>
|
||||
|
|
@ -2525,7 +2525,7 @@
|
|||
</span><span id="Corr-1023"><a href="#Corr-1023"><span class="linenos">1023</span></a>
|
||||
</span><span id="Corr-1024"><a href="#Corr-1024"><span class="linenos">1024</span></a> <span class="k">def</span> <span class="fm">__truediv__</span><span class="p">(</span><span class="bp">self</span><span class="p">,</span> <span class="n">y</span><span class="p">):</span>
|
||||
</span><span id="Corr-1025"><a href="#Corr-1025"><span class="linenos">1025</span></a> <span class="k">if</span> <span class="nb">isinstance</span><span class="p">(</span><span class="n">y</span><span class="p">,</span> <span class="n">Corr</span><span class="p">):</span>
|
||||
</span><span id="Corr-1026"><a href="#Corr-1026"><span class="linenos">1026</span></a> <span class="k">if</span> <span class="ow">not</span><span class="p">((</span><span class="bp">self</span><span class="o">.</span><span class="n">N</span> <span class="o">==</span> <span class="mi">1</span> <span class="ow">or</span> <span class="n">y</span><span class="o">.</span><span class="n">N</span> <span class="o">==</span> <span class="mi">1</span> <span class="ow">or</span> <span class="bp">self</span><span class="o">.</span><span class="n">N</span> <span class="o">==</span> <span class="n">y</span><span class="o">.</span><span class="n">N</span><span class="p">)</span> <span class="ow">and</span> <span class="bp">self</span><span class="o">.</span><span class="n">T</span> <span class="o">==</span> <span class="n">y</span><span class="o">.</span><span class="n">T</span><span class="p">):</span>
|
||||
</span><span id="Corr-1026"><a href="#Corr-1026"><span class="linenos">1026</span></a> <span class="k">if</span> <span class="ow">not</span> <span class="p">((</span><span class="bp">self</span><span class="o">.</span><span class="n">N</span> <span class="o">==</span> <span class="mi">1</span> <span class="ow">or</span> <span class="n">y</span><span class="o">.</span><span class="n">N</span> <span class="o">==</span> <span class="mi">1</span> <span class="ow">or</span> <span class="bp">self</span><span class="o">.</span><span class="n">N</span> <span class="o">==</span> <span class="n">y</span><span class="o">.</span><span class="n">N</span><span class="p">)</span> <span class="ow">and</span> <span class="bp">self</span><span class="o">.</span><span class="n">T</span> <span class="o">==</span> <span class="n">y</span><span class="o">.</span><span class="n">T</span><span class="p">):</span>
|
||||
</span><span id="Corr-1027"><a href="#Corr-1027"><span class="linenos">1027</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"Multiplication of Corr object requires N=N or N=1 and T=T"</span><span class="p">)</span>
|
||||
</span><span id="Corr-1028"><a href="#Corr-1028"><span class="linenos">1028</span></a> <span class="n">newcontent</span> <span class="o">=</span> <span class="p">[]</span>
|
||||
</span><span id="Corr-1029"><a href="#Corr-1029"><span class="linenos">1029</span></a> <span class="k">for</span> <span class="n">t</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">T</span><span class="p">):</span>
|
||||
|
|
@ -2888,7 +2888,7 @@ region indentified for this correlator.</li>
|
|||
<div class="pdoc-code codehilite"><pre><span></span><span id="Corr.gamma_method-120"><a href="#Corr.gamma_method-120"><span class="linenos">120</span></a> <span class="k">def</span> <span class="nf">gamma_method</span><span class="p">(</span><span class="bp">self</span><span class="p">,</span> <span class="o">**</span><span class="n">kwargs</span><span class="p">):</span>
|
||||
</span><span id="Corr.gamma_method-121"><a href="#Corr.gamma_method-121"><span class="linenos">121</span></a> <span class="sd">"""Apply the gamma method to the content of the Corr."""</span>
|
||||
</span><span id="Corr.gamma_method-122"><a href="#Corr.gamma_method-122"><span class="linenos">122</span></a> <span class="k">for</span> <span class="n">item</span> <span class="ow">in</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">:</span>
|
||||
</span><span id="Corr.gamma_method-123"><a href="#Corr.gamma_method-123"><span class="linenos">123</span></a> <span class="k">if</span> <span class="ow">not</span><span class="p">(</span><span class="n">item</span> <span class="ow">is</span> <span class="kc">None</span><span class="p">):</span>
|
||||
</span><span id="Corr.gamma_method-123"><a href="#Corr.gamma_method-123"><span class="linenos">123</span></a> <span class="k">if</span> <span class="ow">not</span> <span class="p">(</span><span class="n">item</span> <span class="ow">is</span> <span class="kc">None</span><span class="p">):</span>
|
||||
</span><span id="Corr.gamma_method-124"><a href="#Corr.gamma_method-124"><span class="linenos">124</span></a> <span class="k">if</span> <span class="bp">self</span><span class="o">.</span><span class="n">N</span> <span class="o">==</span> <span class="mi">1</span><span class="p">:</span>
|
||||
</span><span id="Corr.gamma_method-125"><a href="#Corr.gamma_method-125"><span class="linenos">125</span></a> <span class="n">item</span><span class="p">[</span><span class="mi">0</span><span class="p">]</span><span class="o">.</span><span class="n">gamma_method</span><span class="p">(</span><span class="o">**</span><span class="n">kwargs</span><span class="p">)</span>
|
||||
</span><span id="Corr.gamma_method-126"><a href="#Corr.gamma_method-126"><span class="linenos">126</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
|
|
@ -2926,7 +2926,7 @@ region indentified for this correlator.</li>
|
|||
</span><span id="Corr.projected-140"><a href="#Corr.projected-140"><span class="linenos">140</span></a>
|
||||
</span><span id="Corr.projected-141"><a href="#Corr.projected-141"><span class="linenos">141</span></a> <span class="k">if</span> <span class="n">vector_l</span> <span class="ow">is</span> <span class="kc">None</span><span class="p">:</span>
|
||||
</span><span id="Corr.projected-142"><a href="#Corr.projected-142"><span class="linenos">142</span></a> <span class="n">vector_l</span><span class="p">,</span> <span class="n">vector_r</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">asarray</span><span class="p">([</span><span class="mf">1.</span><span class="p">]</span> <span class="o">+</span> <span class="p">(</span><span class="bp">self</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="p">[</span><span class="mf">0.</span><span class="p">]),</span> <span class="n">np</span><span class="o">.</span><span class="n">asarray</span><span class="p">([</span><span class="mf">1.</span><span class="p">]</span> <span class="o">+</span> <span class="p">(</span><span class="bp">self</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="p">[</span><span class="mf">0.</span><span class="p">])</span>
|
||||
</span><span id="Corr.projected-143"><a href="#Corr.projected-143"><span class="linenos">143</span></a> <span class="k">elif</span><span class="p">(</span><span class="n">vector_r</span> <span class="ow">is</span> <span class="kc">None</span><span class="p">):</span>
|
||||
</span><span id="Corr.projected-143"><a href="#Corr.projected-143"><span class="linenos">143</span></a> <span class="k">elif</span> <span class="p">(</span><span class="n">vector_r</span> <span class="ow">is</span> <span class="kc">None</span><span class="p">):</span>
|
||||
</span><span id="Corr.projected-144"><a href="#Corr.projected-144"><span class="linenos">144</span></a> <span class="n">vector_r</span> <span class="o">=</span> <span class="n">vector_l</span>
|
||||
</span><span id="Corr.projected-145"><a href="#Corr.projected-145"><span class="linenos">145</span></a> <span class="k">if</span> <span class="nb">isinstance</span><span class="p">(</span><span class="n">vector_l</span><span class="p">,</span> <span class="nb">list</span><span class="p">)</span> <span class="ow">and</span> <span class="ow">not</span> <span class="nb">isinstance</span><span class="p">(</span><span class="n">vector_r</span><span class="p">,</span> <span class="nb">list</span><span class="p">):</span>
|
||||
</span><span id="Corr.projected-146"><a href="#Corr.projected-146"><span class="linenos">146</span></a> <span class="k">if</span> <span class="nb">len</span><span class="p">(</span><span class="n">vector_l</span><span class="p">)</span> <span class="o">!=</span> <span class="bp">self</span><span class="o">.</span><span class="n">T</span><span class="p">:</span>
|
||||
|
|
@ -2987,7 +2987,7 @@ By default it will return the lowest source, which usually means unsmeared-unsme
|
|||
</span><span id="Corr.item-179"><a href="#Corr.item-179"><span class="linenos">179</span></a><span class="sd"> """</span>
|
||||
</span><span id="Corr.item-180"><a href="#Corr.item-180"><span class="linenos">180</span></a> <span class="k">if</span> <span class="bp">self</span><span class="o">.</span><span class="n">N</span> <span class="o">==</span> <span class="mi">1</span><span class="p">:</span>
|
||||
</span><span id="Corr.item-181"><a href="#Corr.item-181"><span class="linenos">181</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"Trying to pick item from projected Corr"</span><span class="p">)</span>
|
||||
</span><span id="Corr.item-182"><a href="#Corr.item-182"><span class="linenos">182</span></a> <span class="n">newcontent</span> <span class="o">=</span> <span class="p">[</span><span class="kc">None</span> <span class="k">if</span><span class="p">(</span><span class="n">item</span> <span class="ow">is</span> <span class="kc">None</span><span class="p">)</span> <span class="k">else</span> <span class="n">item</span><span class="p">[</span><span class="n">i</span><span class="p">,</span> <span class="n">j</span><span class="p">]</span> <span class="k">for</span> <span class="n">item</span> <span class="ow">in</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">]</span>
|
||||
</span><span id="Corr.item-182"><a href="#Corr.item-182"><span class="linenos">182</span></a> <span class="n">newcontent</span> <span class="o">=</span> <span class="p">[</span><span class="kc">None</span> <span class="k">if</span> <span class="p">(</span><span class="n">item</span> <span class="ow">is</span> <span class="kc">None</span><span class="p">)</span> <span class="k">else</span> <span class="n">item</span><span class="p">[</span><span class="n">i</span><span class="p">,</span> <span class="n">j</span><span class="p">]</span> <span class="k">for</span> <span class="n">item</span> <span class="ow">in</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">]</span>
|
||||
</span><span id="Corr.item-183"><a href="#Corr.item-183"><span class="linenos">183</span></a> <span class="k">return</span> <span class="n">Corr</span><span class="p">(</span><span class="n">newcontent</span><span class="p">)</span>
|
||||
</span></pre></div>
|
||||
|
||||
|
|
@ -3068,7 +3068,7 @@ timeslice and the error on each timeslice.</p>
|
|||
</span><span id="Corr.symmetric-212"><a href="#Corr.symmetric-212"><span class="linenos">212</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="kc">None</span><span class="p">)</span>
|
||||
</span><span id="Corr.symmetric-213"><a href="#Corr.symmetric-213"><span class="linenos">213</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="Corr.symmetric-214"><a href="#Corr.symmetric-214"><span class="linenos">214</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="mf">0.5</span> <span class="o">*</span> <span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span><span class="p">]</span> <span class="o">+</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="bp">self</span><span class="o">.</span><span class="n">T</span> <span class="o">-</span> <span class="n">t</span><span class="p">]))</span>
|
||||
</span><span id="Corr.symmetric-215"><a href="#Corr.symmetric-215"><span class="linenos">215</span></a> <span class="k">if</span><span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="Corr.symmetric-215"><a href="#Corr.symmetric-215"><span class="linenos">215</span></a> <span class="k">if</span> <span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="Corr.symmetric-216"><a href="#Corr.symmetric-216"><span class="linenos">216</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"Corr could not be symmetrized: No redundant values"</span><span class="p">)</span>
|
||||
</span><span id="Corr.symmetric-217"><a href="#Corr.symmetric-217"><span class="linenos">217</span></a> <span class="k">return</span> <span class="n">Corr</span><span class="p">(</span><span class="n">newcontent</span><span class="p">,</span> <span class="n">prange</span><span class="o">=</span><span class="bp">self</span><span class="o">.</span><span class="n">prange</span><span class="p">)</span>
|
||||
</span></pre></div>
|
||||
|
|
@ -3108,7 +3108,7 @@ timeslice and the error on each timeslice.</p>
|
|||
</span><span id="Corr.anti_symmetric-234"><a href="#Corr.anti_symmetric-234"><span class="linenos">234</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="kc">None</span><span class="p">)</span>
|
||||
</span><span id="Corr.anti_symmetric-235"><a href="#Corr.anti_symmetric-235"><span class="linenos">235</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="Corr.anti_symmetric-236"><a href="#Corr.anti_symmetric-236"><span class="linenos">236</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="mf">0.5</span> <span class="o">*</span> <span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span><span class="p">]</span> <span class="o">-</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="bp">self</span><span class="o">.</span><span class="n">T</span> <span class="o">-</span> <span class="n">t</span><span class="p">]))</span>
|
||||
</span><span id="Corr.anti_symmetric-237"><a href="#Corr.anti_symmetric-237"><span class="linenos">237</span></a> <span class="k">if</span><span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="Corr.anti_symmetric-237"><a href="#Corr.anti_symmetric-237"><span class="linenos">237</span></a> <span class="k">if</span> <span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="Corr.anti_symmetric-238"><a href="#Corr.anti_symmetric-238"><span class="linenos">238</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"Corr could not be symmetrized: No redundant values"</span><span class="p">)</span>
|
||||
</span><span id="Corr.anti_symmetric-239"><a href="#Corr.anti_symmetric-239"><span class="linenos">239</span></a> <span class="k">return</span> <span class="n">Corr</span><span class="p">(</span><span class="n">newcontent</span><span class="p">,</span> <span class="n">prange</span><span class="o">=</span><span class="bp">self</span><span class="o">.</span><span class="n">prange</span><span class="p">)</span>
|
||||
</span></pre></div>
|
||||
|
|
@ -3724,7 +3724,7 @@ Parity quantum number of the correlator, can be +1 or -1</li>
|
|||
</span><span id="Corr.deriv-535"><a href="#Corr.deriv-535"><span class="linenos">535</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="kc">None</span><span class="p">)</span>
|
||||
</span><span id="Corr.deriv-536"><a href="#Corr.deriv-536"><span class="linenos">536</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="Corr.deriv-537"><a href="#Corr.deriv-537"><span class="linenos">537</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="mf">0.5</span> <span class="o">*</span> <span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">+</span> <span class="mi">1</span><span class="p">]</span> <span class="o">-</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">-</span> <span class="mi">1</span><span class="p">]))</span>
|
||||
</span><span id="Corr.deriv-538"><a href="#Corr.deriv-538"><span class="linenos">538</span></a> <span class="k">if</span><span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="Corr.deriv-538"><a href="#Corr.deriv-538"><span class="linenos">538</span></a> <span class="k">if</span> <span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="Corr.deriv-539"><a href="#Corr.deriv-539"><span class="linenos">539</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s1">'Derivative is undefined at all timeslices'</span><span class="p">)</span>
|
||||
</span><span id="Corr.deriv-540"><a href="#Corr.deriv-540"><span class="linenos">540</span></a> <span class="k">return</span> <span class="n">Corr</span><span class="p">(</span><span class="n">newcontent</span><span class="p">,</span> <span class="n">padding</span><span class="o">=</span><span class="p">[</span><span class="mi">1</span><span class="p">,</span> <span class="mi">1</span><span class="p">])</span>
|
||||
</span><span id="Corr.deriv-541"><a href="#Corr.deriv-541"><span class="linenos">541</span></a> <span class="k">elif</span> <span class="n">variant</span> <span class="o">==</span> <span class="s2">"forward"</span><span class="p">:</span>
|
||||
|
|
@ -3734,7 +3734,7 @@ Parity quantum number of the correlator, can be +1 or -1</li>
|
|||
</span><span id="Corr.deriv-545"><a href="#Corr.deriv-545"><span class="linenos">545</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="kc">None</span><span class="p">)</span>
|
||||
</span><span id="Corr.deriv-546"><a href="#Corr.deriv-546"><span class="linenos">546</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="Corr.deriv-547"><a href="#Corr.deriv-547"><span class="linenos">547</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">+</span> <span class="mi">1</span><span class="p">]</span> <span class="o">-</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span><span class="p">])</span>
|
||||
</span><span id="Corr.deriv-548"><a href="#Corr.deriv-548"><span class="linenos">548</span></a> <span class="k">if</span><span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="Corr.deriv-548"><a href="#Corr.deriv-548"><span class="linenos">548</span></a> <span class="k">if</span> <span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="Corr.deriv-549"><a href="#Corr.deriv-549"><span class="linenos">549</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"Derivative is undefined at all timeslices"</span><span class="p">)</span>
|
||||
</span><span id="Corr.deriv-550"><a href="#Corr.deriv-550"><span class="linenos">550</span></a> <span class="k">return</span> <span class="n">Corr</span><span class="p">(</span><span class="n">newcontent</span><span class="p">,</span> <span class="n">padding</span><span class="o">=</span><span class="p">[</span><span class="mi">0</span><span class="p">,</span> <span class="mi">1</span><span class="p">])</span>
|
||||
</span><span id="Corr.deriv-551"><a href="#Corr.deriv-551"><span class="linenos">551</span></a> <span class="k">elif</span> <span class="n">variant</span> <span class="o">==</span> <span class="s2">"backward"</span><span class="p">:</span>
|
||||
|
|
@ -3744,7 +3744,7 @@ Parity quantum number of the correlator, can be +1 or -1</li>
|
|||
</span><span id="Corr.deriv-555"><a href="#Corr.deriv-555"><span class="linenos">555</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="kc">None</span><span class="p">)</span>
|
||||
</span><span id="Corr.deriv-556"><a href="#Corr.deriv-556"><span class="linenos">556</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="Corr.deriv-557"><a href="#Corr.deriv-557"><span class="linenos">557</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span><span class="p">]</span> <span class="o">-</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">-</span> <span class="mi">1</span><span class="p">])</span>
|
||||
</span><span id="Corr.deriv-558"><a href="#Corr.deriv-558"><span class="linenos">558</span></a> <span class="k">if</span><span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="Corr.deriv-558"><a href="#Corr.deriv-558"><span class="linenos">558</span></a> <span class="k">if</span> <span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="Corr.deriv-559"><a href="#Corr.deriv-559"><span class="linenos">559</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"Derivative is undefined at all timeslices"</span><span class="p">)</span>
|
||||
</span><span id="Corr.deriv-560"><a href="#Corr.deriv-560"><span class="linenos">560</span></a> <span class="k">return</span> <span class="n">Corr</span><span class="p">(</span><span class="n">newcontent</span><span class="p">,</span> <span class="n">padding</span><span class="o">=</span><span class="p">[</span><span class="mi">1</span><span class="p">,</span> <span class="mi">0</span><span class="p">])</span>
|
||||
</span><span id="Corr.deriv-561"><a href="#Corr.deriv-561"><span class="linenos">561</span></a> <span class="k">elif</span> <span class="n">variant</span> <span class="o">==</span> <span class="s2">"improved"</span><span class="p">:</span>
|
||||
|
|
@ -3754,7 +3754,7 @@ Parity quantum number of the correlator, can be +1 or -1</li>
|
|||
</span><span id="Corr.deriv-565"><a href="#Corr.deriv-565"><span class="linenos">565</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="kc">None</span><span class="p">)</span>
|
||||
</span><span id="Corr.deriv-566"><a href="#Corr.deriv-566"><span class="linenos">566</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="Corr.deriv-567"><a href="#Corr.deriv-567"><span class="linenos">567</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">((</span><span class="mi">1</span> <span class="o">/</span> <span class="mi">12</span><span class="p">)</span> <span class="o">*</span> <span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">-</span> <span class="mi">2</span><span class="p">]</span> <span class="o">-</span> <span class="mi">8</span> <span class="o">*</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">-</span> <span class="mi">1</span><span class="p">]</span> <span class="o">+</span> <span class="mi">8</span> <span class="o">*</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">+</span> <span class="mi">1</span><span class="p">]</span> <span class="o">-</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">+</span> <span class="mi">2</span><span class="p">]))</span>
|
||||
</span><span id="Corr.deriv-568"><a href="#Corr.deriv-568"><span class="linenos">568</span></a> <span class="k">if</span><span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="Corr.deriv-568"><a href="#Corr.deriv-568"><span class="linenos">568</span></a> <span class="k">if</span> <span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="Corr.deriv-569"><a href="#Corr.deriv-569"><span class="linenos">569</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s1">'Derivative is undefined at all timeslices'</span><span class="p">)</span>
|
||||
</span><span id="Corr.deriv-570"><a href="#Corr.deriv-570"><span class="linenos">570</span></a> <span class="k">return</span> <span class="n">Corr</span><span class="p">(</span><span class="n">newcontent</span><span class="p">,</span> <span class="n">padding</span><span class="o">=</span><span class="p">[</span><span class="mi">2</span><span class="p">,</span> <span class="mi">2</span><span class="p">])</span>
|
||||
</span><span id="Corr.deriv-571"><a href="#Corr.deriv-571"><span class="linenos">571</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
|
|
@ -3804,7 +3804,7 @@ Available choice: symmetric, forward, backward, improved, default: symmetric</li
|
|||
</span><span id="Corr.second_deriv-589"><a href="#Corr.second_deriv-589"><span class="linenos">589</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="kc">None</span><span class="p">)</span>
|
||||
</span><span id="Corr.second_deriv-590"><a href="#Corr.second_deriv-590"><span class="linenos">590</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="Corr.second_deriv-591"><a href="#Corr.second_deriv-591"><span class="linenos">591</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">((</span><span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</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="o">*</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span><span class="p">]</span> <span class="o">+</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">-</span> <span class="mi">1</span><span class="p">]))</span>
|
||||
</span><span id="Corr.second_deriv-592"><a href="#Corr.second_deriv-592"><span class="linenos">592</span></a> <span class="k">if</span><span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="Corr.second_deriv-592"><a href="#Corr.second_deriv-592"><span class="linenos">592</span></a> <span class="k">if</span> <span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="Corr.second_deriv-593"><a href="#Corr.second_deriv-593"><span class="linenos">593</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"Derivative is undefined at all timeslices"</span><span class="p">)</span>
|
||||
</span><span id="Corr.second_deriv-594"><a href="#Corr.second_deriv-594"><span class="linenos">594</span></a> <span class="k">return</span> <span class="n">Corr</span><span class="p">(</span><span class="n">newcontent</span><span class="p">,</span> <span class="n">padding</span><span class="o">=</span><span class="p">[</span><span class="mi">1</span><span class="p">,</span> <span class="mi">1</span><span class="p">])</span>
|
||||
</span><span id="Corr.second_deriv-595"><a href="#Corr.second_deriv-595"><span class="linenos">595</span></a> <span class="k">elif</span> <span class="n">variant</span> <span class="o">==</span> <span class="s2">"improved"</span><span class="p">:</span>
|
||||
|
|
@ -3814,7 +3814,7 @@ Available choice: symmetric, forward, backward, improved, default: symmetric</li
|
|||
</span><span id="Corr.second_deriv-599"><a href="#Corr.second_deriv-599"><span class="linenos">599</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="kc">None</span><span class="p">)</span>
|
||||
</span><span id="Corr.second_deriv-600"><a href="#Corr.second_deriv-600"><span class="linenos">600</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="Corr.second_deriv-601"><a href="#Corr.second_deriv-601"><span class="linenos">601</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">((</span><span class="mi">1</span> <span class="o">/</span> <span class="mi">12</span><span class="p">)</span> <span class="o">*</span> <span class="p">(</span><span class="o">-</span><span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">+</span> <span class="mi">2</span><span class="p">]</span> <span class="o">+</span> <span class="mi">16</span> <span class="o">*</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">+</span> <span class="mi">1</span><span class="p">]</span> <span class="o">-</span> <span class="mi">30</span> <span class="o">*</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span><span class="p">]</span> <span class="o">+</span> <span class="mi">16</span> <span class="o">*</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">-</span> <span class="mi">1</span><span class="p">]</span> <span class="o">-</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">-</span> <span class="mi">2</span><span class="p">]))</span>
|
||||
</span><span id="Corr.second_deriv-602"><a href="#Corr.second_deriv-602"><span class="linenos">602</span></a> <span class="k">if</span><span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="Corr.second_deriv-602"><a href="#Corr.second_deriv-602"><span class="linenos">602</span></a> <span class="k">if</span> <span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="Corr.second_deriv-603"><a href="#Corr.second_deriv-603"><span class="linenos">603</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"Derivative is undefined at all timeslices"</span><span class="p">)</span>
|
||||
</span><span id="Corr.second_deriv-604"><a href="#Corr.second_deriv-604"><span class="linenos">604</span></a> <span class="k">return</span> <span class="n">Corr</span><span class="p">(</span><span class="n">newcontent</span><span class="p">,</span> <span class="n">padding</span><span class="o">=</span><span class="p">[</span><span class="mi">2</span><span class="p">,</span> <span class="mi">2</span><span class="p">])</span>
|
||||
</span><span id="Corr.second_deriv-605"><a href="#Corr.second_deriv-605"><span class="linenos">605</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
|
|
@ -3871,7 +3871,7 @@ Available choice: symmetric, improved, default: symmetric</li>
|
|||
</span><span id="Corr.m_eff-630"><a href="#Corr.m_eff-630"><span class="linenos">630</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="kc">None</span><span class="p">)</span>
|
||||
</span><span id="Corr.m_eff-631"><a href="#Corr.m_eff-631"><span class="linenos">631</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="Corr.m_eff-632"><a href="#Corr.m_eff-632"><span class="linenos">632</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span><span class="p">]</span> <span class="o">/</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">+</span> <span class="mi">1</span><span class="p">])</span>
|
||||
</span><span id="Corr.m_eff-633"><a href="#Corr.m_eff-633"><span class="linenos">633</span></a> <span class="k">if</span><span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="Corr.m_eff-633"><a href="#Corr.m_eff-633"><span class="linenos">633</span></a> <span class="k">if</span> <span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="Corr.m_eff-634"><a href="#Corr.m_eff-634"><span class="linenos">634</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s1">'m_eff is undefined at all timeslices'</span><span class="p">)</span>
|
||||
</span><span id="Corr.m_eff-635"><a href="#Corr.m_eff-635"><span class="linenos">635</span></a>
|
||||
</span><span id="Corr.m_eff-636"><a href="#Corr.m_eff-636"><span class="linenos">636</span></a> <span class="k">return</span> <span class="n">np</span><span class="o">.</span><span class="n">log</span><span class="p">(</span><span class="n">Corr</span><span class="p">(</span><span class="n">newcontent</span><span class="p">,</span> <span class="n">padding</span><span class="o">=</span><span class="p">[</span><span class="mi">0</span><span class="p">,</span> <span class="mi">1</span><span class="p">]))</span>
|
||||
|
|
@ -3896,7 +3896,7 @@ Available choice: symmetric, improved, default: symmetric</li>
|
|||
</span><span id="Corr.m_eff-655"><a href="#Corr.m_eff-655"><span class="linenos">655</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="kc">None</span><span class="p">)</span>
|
||||
</span><span id="Corr.m_eff-656"><a href="#Corr.m_eff-656"><span class="linenos">656</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="Corr.m_eff-657"><a href="#Corr.m_eff-657"><span class="linenos">657</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">np</span><span class="o">.</span><span class="n">abs</span><span class="p">(</span><span class="n">find_root</span><span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span><span class="p">][</span><span class="mi">0</span><span class="p">]</span> <span class="o">/</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">+</span> <span class="mi">1</span><span class="p">][</span><span class="mi">0</span><span class="p">],</span> <span class="n">root_function</span><span class="p">,</span> <span class="n">guess</span><span class="o">=</span><span class="n">guess</span><span class="p">)))</span>
|
||||
</span><span id="Corr.m_eff-658"><a href="#Corr.m_eff-658"><span class="linenos">658</span></a> <span class="k">if</span><span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="Corr.m_eff-658"><a href="#Corr.m_eff-658"><span class="linenos">658</span></a> <span class="k">if</span> <span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="Corr.m_eff-659"><a href="#Corr.m_eff-659"><span class="linenos">659</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s1">'m_eff is undefined at all timeslices'</span><span class="p">)</span>
|
||||
</span><span id="Corr.m_eff-660"><a href="#Corr.m_eff-660"><span class="linenos">660</span></a>
|
||||
</span><span id="Corr.m_eff-661"><a href="#Corr.m_eff-661"><span class="linenos">661</span></a> <span class="k">return</span> <span class="n">Corr</span><span class="p">(</span><span class="n">newcontent</span><span class="p">,</span> <span class="n">padding</span><span class="o">=</span><span class="p">[</span><span class="mi">0</span><span class="p">,</span> <span class="mi">1</span><span class="p">])</span>
|
||||
|
|
@ -3908,7 +3908,7 @@ Available choice: symmetric, improved, default: symmetric</li>
|
|||
</span><span id="Corr.m_eff-667"><a href="#Corr.m_eff-667"><span class="linenos">667</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="kc">None</span><span class="p">)</span>
|
||||
</span><span id="Corr.m_eff-668"><a href="#Corr.m_eff-668"><span class="linenos">668</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="Corr.m_eff-669"><a href="#Corr.m_eff-669"><span class="linenos">669</span></a> <span class="n">newcontent</span><span class="o">.</span><span class="n">append</span><span class="p">((</span><span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">+</span> <span class="mi">1</span><span class="p">]</span> <span class="o">+</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span> <span class="o">-</span> <span class="mi">1</span><span class="p">])</span> <span class="o">/</span> <span class="p">(</span><span class="mi">2</span> <span class="o">*</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span><span class="p">]))</span>
|
||||
</span><span id="Corr.m_eff-670"><a href="#Corr.m_eff-670"><span class="linenos">670</span></a> <span class="k">if</span><span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="Corr.m_eff-670"><a href="#Corr.m_eff-670"><span class="linenos">670</span></a> <span class="k">if</span> <span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="n">x</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">newcontent</span><span class="p">])):</span>
|
||||
</span><span id="Corr.m_eff-671"><a href="#Corr.m_eff-671"><span class="linenos">671</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"m_eff is undefined at all timeslices"</span><span class="p">)</span>
|
||||
</span><span id="Corr.m_eff-672"><a href="#Corr.m_eff-672"><span class="linenos">672</span></a> <span class="k">return</span> <span class="n">np</span><span class="o">.</span><span class="n">arccosh</span><span class="p">(</span><span class="n">Corr</span><span class="p">(</span><span class="n">newcontent</span><span class="p">,</span> <span class="n">padding</span><span class="o">=</span><span class="p">[</span><span class="mi">1</span><span class="p">,</span> <span class="mi">1</span><span class="p">]))</span>
|
||||
</span><span id="Corr.m_eff-673"><a href="#Corr.m_eff-673"><span class="linenos">673</span></a>
|
||||
|
|
@ -4034,7 +4034,7 @@ Decides whether output is printed to the standard output.</li>
|
|||
</span><span id="Corr.plateau-730"><a href="#Corr.plateau-730"><span class="linenos">730</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"no plateau range provided"</span><span class="p">)</span>
|
||||
</span><span id="Corr.plateau-731"><a href="#Corr.plateau-731"><span class="linenos">731</span></a> <span class="k">if</span> <span class="bp">self</span><span class="o">.</span><span class="n">N</span> <span class="o">!=</span> <span class="mi">1</span><span class="p">:</span>
|
||||
</span><span id="Corr.plateau-732"><a href="#Corr.plateau-732"><span class="linenos">732</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"Correlator must be projected before getting a plateau."</span><span class="p">)</span>
|
||||
</span><span id="Corr.plateau-733"><a href="#Corr.plateau-733"><span class="linenos">733</span></a> <span class="k">if</span><span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span><span class="p">]</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">t</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="n">plateau_range</span><span class="p">[</span><span class="mi">0</span><span class="p">],</span> <span class="n">plateau_range</span><span class="p">[</span><span class="mi">1</span><span class="p">]</span> <span class="o">+</span> <span class="mi">1</span><span class="p">)])):</span>
|
||||
</span><span id="Corr.plateau-733"><a href="#Corr.plateau-733"><span class="linenos">733</span></a> <span class="k">if</span> <span class="p">(</span><span class="nb">all</span><span class="p">([</span><span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t</span><span class="p">]</span> <span class="ow">is</span> <span class="kc">None</span> <span class="k">for</span> <span class="n">t</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="n">plateau_range</span><span class="p">[</span><span class="mi">0</span><span class="p">],</span> <span class="n">plateau_range</span><span class="p">[</span><span class="mi">1</span><span class="p">]</span> <span class="o">+</span> <span class="mi">1</span><span class="p">)])):</span>
|
||||
</span><span id="Corr.plateau-734"><a href="#Corr.plateau-734"><span class="linenos">734</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"plateau is undefined at all timeslices in plateaurange."</span><span class="p">)</span>
|
||||
</span><span id="Corr.plateau-735"><a href="#Corr.plateau-735"><span class="linenos">735</span></a> <span class="k">if</span> <span class="n">auto_gamma</span><span class="p">:</span>
|
||||
</span><span id="Corr.plateau-736"><a href="#Corr.plateau-736"><span class="linenos">736</span></a> <span class="bp">self</span><span class="o">.</span><span class="n">gamma_method</span><span class="p">()</span>
|
||||
|
|
|
|||
|
|
@ -245,7 +245,7 @@
|
|||
</span><span id="L-157"><a href="#L-157"><span class="linenos">157</span></a> <span class="n">config_no</span> <span class="o">=</span> <span class="n">struct</span><span class="o">.</span><span class="n">unpack</span><span class="p">(</span><span class="s1">'i'</span><span class="p">,</span> <span class="n">t</span><span class="p">)[</span><span class="mi">0</span><span class="p">]</span>
|
||||
</span><span id="L-158"><a href="#L-158"><span class="linenos">158</span></a> <span class="n">configlist</span><span class="p">[</span><span class="o">-</span><span class="mi">1</span><span class="p">]</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">config_no</span><span class="p">)</span>
|
||||
</span><span id="L-159"><a href="#L-159"><span class="linenos">159</span></a> <span class="k">for</span> <span class="n">i</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="n">nrw</span><span class="p">):</span>
|
||||
</span><span id="L-160"><a href="#L-160"><span class="linenos">160</span></a> <span class="k">if</span><span class="p">(</span><span class="n">version</span> <span class="o">==</span> <span class="s1">'2.0'</span><span class="p">):</span>
|
||||
</span><span id="L-160"><a href="#L-160"><span class="linenos">160</span></a> <span class="k">if</span> <span class="p">(</span><span class="n">version</span> <span class="o">==</span> <span class="s1">'2.0'</span><span class="p">):</span>
|
||||
</span><span id="L-161"><a href="#L-161"><span class="linenos">161</span></a> <span class="n">tmpd</span> <span class="o">=</span> <span class="n">_read_array_openQCD2</span><span class="p">(</span><span class="n">fp</span><span class="p">)</span>
|
||||
</span><span id="L-162"><a href="#L-162"><span class="linenos">162</span></a> <span class="n">tmpd</span> <span class="o">=</span> <span class="n">_read_array_openQCD2</span><span class="p">(</span><span class="n">fp</span><span class="p">)</span>
|
||||
</span><span id="L-163"><a href="#L-163"><span class="linenos">163</span></a> <span class="n">tmp_rw</span> <span class="o">=</span> <span class="n">tmpd</span><span class="p">[</span><span class="s1">'arr'</span><span class="p">]</span>
|
||||
|
|
@ -452,7 +452,7 @@
|
|||
</span><span id="L-364"><a href="#L-364"><span class="linenos">364</span></a> <span class="n">configlist</span><span class="o">.</span><span class="n">append</span><span class="p">([])</span>
|
||||
</span><span id="L-365"><a href="#L-365"><span class="linenos">365</span></a> <span class="k">while</span> <span class="kc">True</span><span class="p">:</span>
|
||||
</span><span id="L-366"><a href="#L-366"><span class="linenos">366</span></a> <span class="n">t</span> <span class="o">=</span> <span class="n">fp</span><span class="o">.</span><span class="n">read</span><span class="p">(</span><span class="mi">4</span><span class="p">)</span>
|
||||
</span><span id="L-367"><a href="#L-367"><span class="linenos">367</span></a> <span class="k">if</span><span class="p">(</span><span class="nb">len</span><span class="p">(</span><span class="n">t</span><span class="p">)</span> <span class="o"><</span> <span class="mi">4</span><span class="p">):</span>
|
||||
</span><span id="L-367"><a href="#L-367"><span class="linenos">367</span></a> <span class="k">if</span> <span class="p">(</span><span class="nb">len</span><span class="p">(</span><span class="n">t</span><span class="p">)</span> <span class="o"><</span> <span class="mi">4</span><span class="p">):</span>
|
||||
</span><span id="L-368"><a href="#L-368"><span class="linenos">368</span></a> <span class="k">break</span>
|
||||
</span><span id="L-369"><a href="#L-369"><span class="linenos">369</span></a> <span class="n">nc</span> <span class="o">=</span> <span class="n">struct</span><span class="o">.</span><span class="n">unpack</span><span class="p">(</span><span class="s1">'i'</span><span class="p">,</span> <span class="n">t</span><span class="p">)[</span><span class="mi">0</span><span class="p">]</span>
|
||||
</span><span id="L-370"><a href="#L-370"><span class="linenos">370</span></a> <span class="n">configlist</span><span class="p">[</span><span class="o">-</span><span class="mi">1</span><span class="p">]</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">nc</span><span class="p">)</span>
|
||||
|
|
@ -842,7 +842,7 @@
|
|||
</span><span id="L-754"><a href="#L-754"><span class="linenos">754</span></a>
|
||||
</span><span id="L-755"><a href="#L-755"><span class="linenos">755</span></a> <span class="n">t</span> <span class="o">=</span> <span class="n">fp</span><span class="o">.</span><span class="n">read</span><span class="p">(</span><span class="mi">12</span><span class="p">)</span>
|
||||
</span><span id="L-756"><a href="#L-756"><span class="linenos">756</span></a> <span class="n">Ls</span> <span class="o">=</span> <span class="n">struct</span><span class="o">.</span><span class="n">unpack</span><span class="p">(</span><span class="s1">'<iii'</span><span class="p">,</span> <span class="n">t</span><span class="p">)</span>
|
||||
</span><span id="L-757"><a href="#L-757"><span class="linenos">757</span></a> <span class="k">if</span><span class="p">(</span><span class="n">Ls</span><span class="p">[</span><span class="mi">0</span><span class="p">]</span> <span class="o">==</span> <span class="n">Ls</span><span class="p">[</span><span class="mi">1</span><span class="p">]</span> <span class="ow">and</span> <span class="n">Ls</span><span class="p">[</span><span class="mi">1</span><span class="p">]</span> <span class="o">==</span> <span class="n">Ls</span><span class="p">[</span><span class="mi">2</span><span class="p">]):</span>
|
||||
</span><span id="L-757"><a href="#L-757"><span class="linenos">757</span></a> <span class="k">if</span> <span class="p">(</span><span class="n">Ls</span><span class="p">[</span><span class="mi">0</span><span class="p">]</span> <span class="o">==</span> <span class="n">Ls</span><span class="p">[</span><span class="mi">1</span><span class="p">]</span> <span class="ow">and</span> <span class="n">Ls</span><span class="p">[</span><span class="mi">1</span><span class="p">]</span> <span class="o">==</span> <span class="n">Ls</span><span class="p">[</span><span class="mi">2</span><span class="p">]):</span>
|
||||
</span><span id="L-758"><a href="#L-758"><span class="linenos">758</span></a> <span class="n">L</span> <span class="o">=</span> <span class="n">Ls</span><span class="p">[</span><span class="mi">0</span><span class="p">]</span>
|
||||
</span><span id="L-759"><a href="#L-759"><span class="linenos">759</span></a> <span class="k">if</span> <span class="ow">not</span> <span class="p">(</span><span class="n">supposed_L</span> <span class="o">==</span> <span class="n">L</span><span class="p">)</span> <span class="ow">and</span> <span class="n">supposed_L</span><span class="p">:</span>
|
||||
</span><span id="L-760"><a href="#L-760"><span class="linenos">760</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"It seems the length given in the header and by you contradict each other"</span><span class="p">)</span>
|
||||
|
|
@ -857,7 +857,7 @@
|
|||
</span><span id="L-769"><a href="#L-769"><span class="linenos">769</span></a> <span class="k">if</span> <span class="n">c</span> <span class="o">></span> <span class="n">cmax</span><span class="p">:</span>
|
||||
</span><span id="L-770"><a href="#L-770"><span class="linenos">770</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s1">'Flow has been determined between c=0 and c=</span><span class="si">%lf</span><span class="s1"> with tolerance </span><span class="si">%lf</span><span class="s1">'</span> <span class="o">%</span> <span class="p">(</span><span class="n">cmax</span><span class="p">,</span> <span class="n">tol</span><span class="p">))</span>
|
||||
</span><span id="L-771"><a href="#L-771"><span class="linenos">771</span></a>
|
||||
</span><span id="L-772"><a href="#L-772"><span class="linenos">772</span></a> <span class="k">if</span><span class="p">(</span><span class="n">zthfl</span> <span class="o">==</span> <span class="mi">2</span><span class="p">):</span>
|
||||
</span><span id="L-772"><a href="#L-772"><span class="linenos">772</span></a> <span class="k">if</span> <span class="p">(</span><span class="n">zthfl</span> <span class="o">==</span> <span class="mi">2</span><span class="p">):</span>
|
||||
</span><span id="L-773"><a href="#L-773"><span class="linenos">773</span></a> <span class="n">nfl</span> <span class="o">=</span> <span class="mi">2</span> <span class="c1"># number of flows</span>
|
||||
</span><span id="L-774"><a href="#L-774"><span class="linenos">774</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="L-775"><a href="#L-775"><span class="linenos">775</span></a> <span class="n">nfl</span> <span class="o">=</span> <span class="mi">1</span>
|
||||
|
|
@ -865,7 +865,7 @@
|
|||
</span><span id="L-777"><a href="#L-777"><span class="linenos">777</span></a>
|
||||
</span><span id="L-778"><a href="#L-778"><span class="linenos">778</span></a> <span class="k">while</span> <span class="kc">True</span><span class="p">:</span>
|
||||
</span><span id="L-779"><a href="#L-779"><span class="linenos">779</span></a> <span class="n">t</span> <span class="o">=</span> <span class="n">fp</span><span class="o">.</span><span class="n">read</span><span class="p">(</span><span class="mi">4</span><span class="p">)</span>
|
||||
</span><span id="L-780"><a href="#L-780"><span class="linenos">780</span></a> <span class="k">if</span><span class="p">(</span><span class="nb">len</span><span class="p">(</span><span class="n">t</span><span class="p">)</span> <span class="o"><</span> <span class="mi">4</span><span class="p">):</span>
|
||||
</span><span id="L-780"><a href="#L-780"><span class="linenos">780</span></a> <span class="k">if</span> <span class="p">(</span><span class="nb">len</span><span class="p">(</span><span class="n">t</span><span class="p">)</span> <span class="o"><</span> <span class="mi">4</span><span class="p">):</span>
|
||||
</span><span id="L-781"><a href="#L-781"><span class="linenos">781</span></a> <span class="k">break</span>
|
||||
</span><span id="L-782"><a href="#L-782"><span class="linenos">782</span></a> <span class="n">traj_list</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">struct</span><span class="o">.</span><span class="n">unpack</span><span class="p">(</span><span class="s1">'i'</span><span class="p">,</span> <span class="n">t</span><span class="p">)[</span><span class="mi">0</span><span class="p">])</span> <span class="c1"># trajectory number when measurement was done</span>
|
||||
</span><span id="L-783"><a href="#L-783"><span class="linenos">783</span></a>
|
||||
|
|
@ -890,7 +890,7 @@
|
|||
</span><span id="L-802"><a href="#L-802"><span class="linenos">802</span></a>
|
||||
</span><span id="L-803"><a href="#L-803"><span class="linenos">803</span></a> <span class="k">while</span> <span class="kc">True</span><span class="p">:</span>
|
||||
</span><span id="L-804"><a href="#L-804"><span class="linenos">804</span></a> <span class="n">t</span> <span class="o">=</span> <span class="n">fp</span><span class="o">.</span><span class="n">read</span><span class="p">(</span><span class="mi">4</span><span class="p">)</span>
|
||||
</span><span id="L-805"><a href="#L-805"><span class="linenos">805</span></a> <span class="k">if</span><span class="p">(</span><span class="nb">len</span><span class="p">(</span><span class="n">t</span><span class="p">)</span> <span class="o"><</span> <span class="mi">4</span><span class="p">):</span>
|
||||
</span><span id="L-805"><a href="#L-805"><span class="linenos">805</span></a> <span class="k">if</span> <span class="p">(</span><span class="nb">len</span><span class="p">(</span><span class="n">t</span><span class="p">)</span> <span class="o"><</span> <span class="mi">4</span><span class="p">):</span>
|
||||
</span><span id="L-806"><a href="#L-806"><span class="linenos">806</span></a> <span class="k">break</span>
|
||||
</span><span id="L-807"><a href="#L-807"><span class="linenos">807</span></a> <span class="n">traj_list</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">struct</span><span class="o">.</span><span class="n">unpack</span><span class="p">(</span><span class="s1">'i'</span><span class="p">,</span> <span class="n">t</span><span class="p">)[</span><span class="mi">0</span><span class="p">])</span>
|
||||
</span><span id="L-808"><a href="#L-808"><span class="linenos">808</span></a> <span class="c1"># Wsl</span>
|
||||
|
|
@ -1224,7 +1224,7 @@
|
|||
</span><span id="read_rwms-158"><a href="#read_rwms-158"><span class="linenos">158</span></a> <span class="n">config_no</span> <span class="o">=</span> <span class="n">struct</span><span class="o">.</span><span class="n">unpack</span><span class="p">(</span><span class="s1">'i'</span><span class="p">,</span> <span class="n">t</span><span class="p">)[</span><span class="mi">0</span><span class="p">]</span>
|
||||
</span><span id="read_rwms-159"><a href="#read_rwms-159"><span class="linenos">159</span></a> <span class="n">configlist</span><span class="p">[</span><span class="o">-</span><span class="mi">1</span><span class="p">]</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">config_no</span><span class="p">)</span>
|
||||
</span><span id="read_rwms-160"><a href="#read_rwms-160"><span class="linenos">160</span></a> <span class="k">for</span> <span class="n">i</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="n">nrw</span><span class="p">):</span>
|
||||
</span><span id="read_rwms-161"><a href="#read_rwms-161"><span class="linenos">161</span></a> <span class="k">if</span><span class="p">(</span><span class="n">version</span> <span class="o">==</span> <span class="s1">'2.0'</span><span class="p">):</span>
|
||||
</span><span id="read_rwms-161"><a href="#read_rwms-161"><span class="linenos">161</span></a> <span class="k">if</span> <span class="p">(</span><span class="n">version</span> <span class="o">==</span> <span class="s1">'2.0'</span><span class="p">):</span>
|
||||
</span><span id="read_rwms-162"><a href="#read_rwms-162"><span class="linenos">162</span></a> <span class="n">tmpd</span> <span class="o">=</span> <span class="n">_read_array_openQCD2</span><span class="p">(</span><span class="n">fp</span><span class="p">)</span>
|
||||
</span><span id="read_rwms-163"><a href="#read_rwms-163"><span class="linenos">163</span></a> <span class="n">tmpd</span> <span class="o">=</span> <span class="n">_read_array_openQCD2</span><span class="p">(</span><span class="n">fp</span><span class="p">)</span>
|
||||
</span><span id="read_rwms-164"><a href="#read_rwms-164"><span class="linenos">164</span></a> <span class="n">tmp_rw</span> <span class="o">=</span> <span class="n">tmpd</span><span class="p">[</span><span class="s1">'arr'</span><span class="p">]</span>
|
||||
|
|
@ -1478,7 +1478,7 @@ Print additional information that is useful for debugging.</li>
|
|||
</span><span id="extract_t0-365"><a href="#extract_t0-365"><span class="linenos">365</span></a> <span class="n">configlist</span><span class="o">.</span><span class="n">append</span><span class="p">([])</span>
|
||||
</span><span id="extract_t0-366"><a href="#extract_t0-366"><span class="linenos">366</span></a> <span class="k">while</span> <span class="kc">True</span><span class="p">:</span>
|
||||
</span><span id="extract_t0-367"><a href="#extract_t0-367"><span class="linenos">367</span></a> <span class="n">t</span> <span class="o">=</span> <span class="n">fp</span><span class="o">.</span><span class="n">read</span><span class="p">(</span><span class="mi">4</span><span class="p">)</span>
|
||||
</span><span id="extract_t0-368"><a href="#extract_t0-368"><span class="linenos">368</span></a> <span class="k">if</span><span class="p">(</span><span class="nb">len</span><span class="p">(</span><span class="n">t</span><span class="p">)</span> <span class="o"><</span> <span class="mi">4</span><span class="p">):</span>
|
||||
</span><span id="extract_t0-368"><a href="#extract_t0-368"><span class="linenos">368</span></a> <span class="k">if</span> <span class="p">(</span><span class="nb">len</span><span class="p">(</span><span class="n">t</span><span class="p">)</span> <span class="o"><</span> <span class="mi">4</span><span class="p">):</span>
|
||||
</span><span id="extract_t0-369"><a href="#extract_t0-369"><span class="linenos">369</span></a> <span class="k">break</span>
|
||||
</span><span id="extract_t0-370"><a href="#extract_t0-370"><span class="linenos">370</span></a> <span class="n">nc</span> <span class="o">=</span> <span class="n">struct</span><span class="o">.</span><span class="n">unpack</span><span class="p">(</span><span class="s1">'i'</span><span class="p">,</span> <span class="n">t</span><span class="p">)[</span><span class="mi">0</span><span class="p">]</span>
|
||||
</span><span id="extract_t0-371"><a href="#extract_t0-371"><span class="linenos">371</span></a> <span class="n">configlist</span><span class="p">[</span><span class="o">-</span><span class="mi">1</span><span class="p">]</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">nc</span><span class="p">)</span>
|
||||
|
|
|
|||
|
|
@ -147,11 +147,11 @@
|
|||
</span><span id="L-74"><a href="#L-74"><span class="linenos"> 74</span></a>
|
||||
</span><span id="L-75"><a href="#L-75"><span class="linenos"> 75</span></a> <span class="k">if</span> <span class="n">version</span> <span class="ow">not</span> <span class="ow">in</span> <span class="n">known_versions</span><span class="p">:</span>
|
||||
</span><span id="L-76"><a href="#L-76"><span class="linenos"> 76</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"This version is not known!"</span><span class="p">)</span>
|
||||
</span><span id="L-77"><a href="#L-77"><span class="linenos"> 77</span></a> <span class="k">if</span><span class="p">(</span><span class="n">version</span><span class="p">[</span><span class="o">-</span><span class="mi">1</span><span class="p">]</span> <span class="o">==</span> <span class="s2">"c"</span><span class="p">):</span>
|
||||
</span><span id="L-77"><a href="#L-77"><span class="linenos"> 77</span></a> <span class="k">if</span> <span class="p">(</span><span class="n">version</span><span class="p">[</span><span class="o">-</span><span class="mi">1</span><span class="p">]</span> <span class="o">==</span> <span class="s2">"c"</span><span class="p">):</span>
|
||||
</span><span id="L-78"><a href="#L-78"><span class="linenos"> 78</span></a> <span class="n">appended</span> <span class="o">=</span> <span class="kc">False</span>
|
||||
</span><span id="L-79"><a href="#L-79"><span class="linenos"> 79</span></a> <span class="n">compact</span> <span class="o">=</span> <span class="kc">True</span>
|
||||
</span><span id="L-80"><a href="#L-80"><span class="linenos"> 80</span></a> <span class="n">version</span> <span class="o">=</span> <span class="n">version</span><span class="p">[:</span><span class="o">-</span><span class="mi">1</span><span class="p">]</span>
|
||||
</span><span id="L-81"><a href="#L-81"><span class="linenos"> 81</span></a> <span class="k">elif</span><span class="p">(</span><span class="n">version</span><span class="p">[</span><span class="o">-</span><span class="mi">1</span><span class="p">]</span> <span class="o">==</span> <span class="s2">"a"</span><span class="p">):</span>
|
||||
</span><span id="L-81"><a href="#L-81"><span class="linenos"> 81</span></a> <span class="k">elif</span> <span class="p">(</span><span class="n">version</span><span class="p">[</span><span class="o">-</span><span class="mi">1</span><span class="p">]</span> <span class="o">==</span> <span class="s2">"a"</span><span class="p">):</span>
|
||||
</span><span id="L-82"><a href="#L-82"><span class="linenos"> 82</span></a> <span class="n">appended</span> <span class="o">=</span> <span class="kc">True</span>
|
||||
</span><span id="L-83"><a href="#L-83"><span class="linenos"> 83</span></a> <span class="n">compact</span> <span class="o">=</span> <span class="kc">False</span>
|
||||
</span><span id="L-84"><a href="#L-84"><span class="linenos"> 84</span></a> <span class="n">version</span> <span class="o">=</span> <span class="n">version</span><span class="p">[:</span><span class="o">-</span><span class="mi">1</span><span class="p">]</span>
|
||||
|
|
@ -325,7 +325,7 @@
|
|||
</span><span id="L-252"><a href="#L-252"><span class="linenos">252</span></a> <span class="n">lines</span> <span class="o">=</span> <span class="n">fp</span><span class="o">.</span><span class="n">readlines</span><span class="p">()</span>
|
||||
</span><span id="L-253"><a href="#L-253"><span class="linenos">253</span></a> <span class="c1"># check, if the correlator is in fact</span>
|
||||
</span><span id="L-254"><a href="#L-254"><span class="linenos">254</span></a> <span class="c1"># printed completely</span>
|
||||
</span><span id="L-255"><a href="#L-255"><span class="linenos">255</span></a> <span class="k">if</span><span class="p">(</span><span class="n">start_read</span> <span class="o">+</span> <span class="n">T</span> <span class="o">></span> <span class="nb">len</span><span class="p">(</span><span class="n">lines</span><span class="p">)):</span>
|
||||
</span><span id="L-255"><a href="#L-255"><span class="linenos">255</span></a> <span class="k">if</span> <span class="p">(</span><span class="n">start_read</span> <span class="o">+</span> <span class="n">T</span> <span class="o">></span> <span class="nb">len</span><span class="p">(</span><span class="n">lines</span><span class="p">)):</span>
|
||||
</span><span id="L-256"><a href="#L-256"><span class="linenos">256</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"EOF before end of correlator data! Maybe "</span> <span class="o">+</span> <span class="n">path</span> <span class="o">+</span> <span class="s1">'/'</span> <span class="o">+</span> <span class="n">item</span> <span class="o">+</span> <span class="s1">'/'</span> <span class="o">+</span> <span class="n">sub_ls</span><span class="p">[</span><span class="n">cfg</span><span class="p">]</span> <span class="o">+</span> <span class="s2">" is corrupted?"</span><span class="p">)</span>
|
||||
</span><span id="L-257"><a href="#L-257"><span class="linenos">257</span></a> <span class="c1"># and start to read the correlator.</span>
|
||||
</span><span id="L-258"><a href="#L-258"><span class="linenos">258</span></a> <span class="c1"># the range here is chosen like this,</span>
|
||||
|
|
@ -335,7 +335,7 @@
|
|||
</span><span id="L-262"><a href="#L-262"><span class="linenos">262</span></a> <span class="k">if</span> <span class="n">k</span> <span class="o">==</span> <span class="n">start_read</span> <span class="o">-</span> <span class="mi">5</span> <span class="o">-</span> <span class="n">b2b</span><span class="p">:</span>
|
||||
</span><span id="L-263"><a href="#L-263"><span class="linenos">263</span></a> <span class="k">if</span> <span class="n">lines</span><span class="p">[</span><span class="n">k</span><span class="p">]</span><span class="o">.</span><span class="n">strip</span><span class="p">()</span> <span class="o">!=</span> <span class="s1">'name '</span> <span class="o">+</span> <span class="n">name</span><span class="p">:</span>
|
||||
</span><span id="L-264"><a href="#L-264"><span class="linenos">264</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s1">'Wrong format'</span><span class="p">,</span> <span class="n">sub_ls</span><span class="p">[</span><span class="n">cfg</span><span class="p">])</span>
|
||||
</span><span id="L-265"><a href="#L-265"><span class="linenos">265</span></a> <span class="k">if</span><span class="p">(</span><span class="n">k</span> <span class="o">>=</span> <span class="n">start_read</span> <span class="ow">and</span> <span class="n">k</span> <span class="o"><</span> <span class="n">start_read</span> <span class="o">+</span> <span class="n">T</span><span class="p">):</span>
|
||||
</span><span id="L-265"><a href="#L-265"><span class="linenos">265</span></a> <span class="k">if</span> <span class="p">(</span><span class="n">k</span> <span class="o">>=</span> <span class="n">start_read</span> <span class="ow">and</span> <span class="n">k</span> <span class="o"><</span> <span class="n">start_read</span> <span class="o">+</span> <span class="n">T</span><span class="p">):</span>
|
||||
</span><span id="L-266"><a href="#L-266"><span class="linenos">266</span></a> <span class="n">floats</span> <span class="o">=</span> <span class="nb">list</span><span class="p">(</span><span class="nb">map</span><span class="p">(</span><span class="nb">float</span><span class="p">,</span> <span class="n">lines</span><span class="p">[</span><span class="n">k</span><span class="p">]</span><span class="o">.</span><span class="n">split</span><span class="p">()))</span>
|
||||
</span><span id="L-267"><a href="#L-267"><span class="linenos">267</span></a> <span class="n">deltas</span><span class="p">[</span><span class="n">k</span> <span class="o">-</span> <span class="n">start_read</span><span class="p">][</span><span class="n">i</span><span class="p">][</span><span class="n">cfg</span><span class="p">]</span> <span class="o">=</span> <span class="n">floats</span><span class="p">[</span><span class="o">-</span><span class="mi">2</span><span class="p">:][</span><span class="n">im</span><span class="p">]</span>
|
||||
</span><span id="L-268"><a href="#L-268"><span class="linenos">268</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
|
|
@ -346,7 +346,7 @@
|
|||
</span><span id="L-273"><a href="#L-273"><span class="linenos">273</span></a> <span class="c1"># we can iterate over the whole file.</span>
|
||||
</span><span id="L-274"><a href="#L-274"><span class="linenos">274</span></a> <span class="c1"># here one can also implement the chekc from above.</span>
|
||||
</span><span id="L-275"><a href="#L-275"><span class="linenos">275</span></a> <span class="k">for</span> <span class="n">k</span><span class="p">,</span> <span class="n">line</span> <span class="ow">in</span> <span class="nb">enumerate</span><span class="p">(</span><span class="n">fp</span><span class="p">):</span>
|
||||
</span><span id="L-276"><a href="#L-276"><span class="linenos">276</span></a> <span class="k">if</span><span class="p">(</span><span class="n">k</span> <span class="o">>=</span> <span class="n">start_read</span> <span class="ow">and</span> <span class="n">k</span> <span class="o"><</span> <span class="n">start_read</span> <span class="o">+</span> <span class="n">T</span><span class="p">):</span>
|
||||
</span><span id="L-276"><a href="#L-276"><span class="linenos">276</span></a> <span class="k">if</span> <span class="p">(</span><span class="n">k</span> <span class="o">>=</span> <span class="n">start_read</span> <span class="ow">and</span> <span class="n">k</span> <span class="o"><</span> <span class="n">start_read</span> <span class="o">+</span> <span class="n">T</span><span class="p">):</span>
|
||||
</span><span id="L-277"><a href="#L-277"><span class="linenos">277</span></a> <span class="n">floats</span> <span class="o">=</span> <span class="nb">list</span><span class="p">(</span><span class="nb">map</span><span class="p">(</span><span class="nb">float</span><span class="p">,</span> <span class="n">line</span><span class="o">.</span><span class="n">split</span><span class="p">()))</span>
|
||||
</span><span id="L-278"><a href="#L-278"><span class="linenos">278</span></a> <span class="k">if</span> <span class="n">version</span> <span class="o">==</span> <span class="s2">"0.0"</span><span class="p">:</span>
|
||||
</span><span id="L-279"><a href="#L-279"><span class="linenos">279</span></a> <span class="n">deltas</span><span class="p">[</span><span class="n">k</span> <span class="o">-</span> <span class="n">start</span><span class="p">][</span><span class="n">i</span><span class="p">][</span><span class="n">cnfg</span><span class="p">]</span> <span class="o">=</span> <span class="n">floats</span><span class="p">[</span><span class="n">im</span> <span class="o">-</span> <span class="n">single</span><span class="p">]</span>
|
||||
|
|
@ -511,11 +511,11 @@
|
|||
</span><span id="read_sfcf-75"><a href="#read_sfcf-75"><span class="linenos"> 75</span></a>
|
||||
</span><span id="read_sfcf-76"><a href="#read_sfcf-76"><span class="linenos"> 76</span></a> <span class="k">if</span> <span class="n">version</span> <span class="ow">not</span> <span class="ow">in</span> <span class="n">known_versions</span><span class="p">:</span>
|
||||
</span><span id="read_sfcf-77"><a href="#read_sfcf-77"><span class="linenos"> 77</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"This version is not known!"</span><span class="p">)</span>
|
||||
</span><span id="read_sfcf-78"><a href="#read_sfcf-78"><span class="linenos"> 78</span></a> <span class="k">if</span><span class="p">(</span><span class="n">version</span><span class="p">[</span><span class="o">-</span><span class="mi">1</span><span class="p">]</span> <span class="o">==</span> <span class="s2">"c"</span><span class="p">):</span>
|
||||
</span><span id="read_sfcf-78"><a href="#read_sfcf-78"><span class="linenos"> 78</span></a> <span class="k">if</span> <span class="p">(</span><span class="n">version</span><span class="p">[</span><span class="o">-</span><span class="mi">1</span><span class="p">]</span> <span class="o">==</span> <span class="s2">"c"</span><span class="p">):</span>
|
||||
</span><span id="read_sfcf-79"><a href="#read_sfcf-79"><span class="linenos"> 79</span></a> <span class="n">appended</span> <span class="o">=</span> <span class="kc">False</span>
|
||||
</span><span id="read_sfcf-80"><a href="#read_sfcf-80"><span class="linenos"> 80</span></a> <span class="n">compact</span> <span class="o">=</span> <span class="kc">True</span>
|
||||
</span><span id="read_sfcf-81"><a href="#read_sfcf-81"><span class="linenos"> 81</span></a> <span class="n">version</span> <span class="o">=</span> <span class="n">version</span><span class="p">[:</span><span class="o">-</span><span class="mi">1</span><span class="p">]</span>
|
||||
</span><span id="read_sfcf-82"><a href="#read_sfcf-82"><span class="linenos"> 82</span></a> <span class="k">elif</span><span class="p">(</span><span class="n">version</span><span class="p">[</span><span class="o">-</span><span class="mi">1</span><span class="p">]</span> <span class="o">==</span> <span class="s2">"a"</span><span class="p">):</span>
|
||||
</span><span id="read_sfcf-82"><a href="#read_sfcf-82"><span class="linenos"> 82</span></a> <span class="k">elif</span> <span class="p">(</span><span class="n">version</span><span class="p">[</span><span class="o">-</span><span class="mi">1</span><span class="p">]</span> <span class="o">==</span> <span class="s2">"a"</span><span class="p">):</span>
|
||||
</span><span id="read_sfcf-83"><a href="#read_sfcf-83"><span class="linenos"> 83</span></a> <span class="n">appended</span> <span class="o">=</span> <span class="kc">True</span>
|
||||
</span><span id="read_sfcf-84"><a href="#read_sfcf-84"><span class="linenos"> 84</span></a> <span class="n">compact</span> <span class="o">=</span> <span class="kc">False</span>
|
||||
</span><span id="read_sfcf-85"><a href="#read_sfcf-85"><span class="linenos"> 85</span></a> <span class="n">version</span> <span class="o">=</span> <span class="n">version</span><span class="p">[:</span><span class="o">-</span><span class="mi">1</span><span class="p">]</span>
|
||||
|
|
@ -689,7 +689,7 @@
|
|||
</span><span id="read_sfcf-253"><a href="#read_sfcf-253"><span class="linenos">253</span></a> <span class="n">lines</span> <span class="o">=</span> <span class="n">fp</span><span class="o">.</span><span class="n">readlines</span><span class="p">()</span>
|
||||
</span><span id="read_sfcf-254"><a href="#read_sfcf-254"><span class="linenos">254</span></a> <span class="c1"># check, if the correlator is in fact</span>
|
||||
</span><span id="read_sfcf-255"><a href="#read_sfcf-255"><span class="linenos">255</span></a> <span class="c1"># printed completely</span>
|
||||
</span><span id="read_sfcf-256"><a href="#read_sfcf-256"><span class="linenos">256</span></a> <span class="k">if</span><span class="p">(</span><span class="n">start_read</span> <span class="o">+</span> <span class="n">T</span> <span class="o">></span> <span class="nb">len</span><span class="p">(</span><span class="n">lines</span><span class="p">)):</span>
|
||||
</span><span id="read_sfcf-256"><a href="#read_sfcf-256"><span class="linenos">256</span></a> <span class="k">if</span> <span class="p">(</span><span class="n">start_read</span> <span class="o">+</span> <span class="n">T</span> <span class="o">></span> <span class="nb">len</span><span class="p">(</span><span class="n">lines</span><span class="p">)):</span>
|
||||
</span><span id="read_sfcf-257"><a href="#read_sfcf-257"><span class="linenos">257</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"EOF before end of correlator data! Maybe "</span> <span class="o">+</span> <span class="n">path</span> <span class="o">+</span> <span class="s1">'/'</span> <span class="o">+</span> <span class="n">item</span> <span class="o">+</span> <span class="s1">'/'</span> <span class="o">+</span> <span class="n">sub_ls</span><span class="p">[</span><span class="n">cfg</span><span class="p">]</span> <span class="o">+</span> <span class="s2">" is corrupted?"</span><span class="p">)</span>
|
||||
</span><span id="read_sfcf-258"><a href="#read_sfcf-258"><span class="linenos">258</span></a> <span class="c1"># and start to read the correlator.</span>
|
||||
</span><span id="read_sfcf-259"><a href="#read_sfcf-259"><span class="linenos">259</span></a> <span class="c1"># the range here is chosen like this,</span>
|
||||
|
|
@ -699,7 +699,7 @@
|
|||
</span><span id="read_sfcf-263"><a href="#read_sfcf-263"><span class="linenos">263</span></a> <span class="k">if</span> <span class="n">k</span> <span class="o">==</span> <span class="n">start_read</span> <span class="o">-</span> <span class="mi">5</span> <span class="o">-</span> <span class="n">b2b</span><span class="p">:</span>
|
||||
</span><span id="read_sfcf-264"><a href="#read_sfcf-264"><span class="linenos">264</span></a> <span class="k">if</span> <span class="n">lines</span><span class="p">[</span><span class="n">k</span><span class="p">]</span><span class="o">.</span><span class="n">strip</span><span class="p">()</span> <span class="o">!=</span> <span class="s1">'name '</span> <span class="o">+</span> <span class="n">name</span><span class="p">:</span>
|
||||
</span><span id="read_sfcf-265"><a href="#read_sfcf-265"><span class="linenos">265</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s1">'Wrong format'</span><span class="p">,</span> <span class="n">sub_ls</span><span class="p">[</span><span class="n">cfg</span><span class="p">])</span>
|
||||
</span><span id="read_sfcf-266"><a href="#read_sfcf-266"><span class="linenos">266</span></a> <span class="k">if</span><span class="p">(</span><span class="n">k</span> <span class="o">>=</span> <span class="n">start_read</span> <span class="ow">and</span> <span class="n">k</span> <span class="o"><</span> <span class="n">start_read</span> <span class="o">+</span> <span class="n">T</span><span class="p">):</span>
|
||||
</span><span id="read_sfcf-266"><a href="#read_sfcf-266"><span class="linenos">266</span></a> <span class="k">if</span> <span class="p">(</span><span class="n">k</span> <span class="o">>=</span> <span class="n">start_read</span> <span class="ow">and</span> <span class="n">k</span> <span class="o"><</span> <span class="n">start_read</span> <span class="o">+</span> <span class="n">T</span><span class="p">):</span>
|
||||
</span><span id="read_sfcf-267"><a href="#read_sfcf-267"><span class="linenos">267</span></a> <span class="n">floats</span> <span class="o">=</span> <span class="nb">list</span><span class="p">(</span><span class="nb">map</span><span class="p">(</span><span class="nb">float</span><span class="p">,</span> <span class="n">lines</span><span class="p">[</span><span class="n">k</span><span class="p">]</span><span class="o">.</span><span class="n">split</span><span class="p">()))</span>
|
||||
</span><span id="read_sfcf-268"><a href="#read_sfcf-268"><span class="linenos">268</span></a> <span class="n">deltas</span><span class="p">[</span><span class="n">k</span> <span class="o">-</span> <span class="n">start_read</span><span class="p">][</span><span class="n">i</span><span class="p">][</span><span class="n">cfg</span><span class="p">]</span> <span class="o">=</span> <span class="n">floats</span><span class="p">[</span><span class="o">-</span><span class="mi">2</span><span class="p">:][</span><span class="n">im</span><span class="p">]</span>
|
||||
</span><span id="read_sfcf-269"><a href="#read_sfcf-269"><span class="linenos">269</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
|
|
@ -710,7 +710,7 @@
|
|||
</span><span id="read_sfcf-274"><a href="#read_sfcf-274"><span class="linenos">274</span></a> <span class="c1"># we can iterate over the whole file.</span>
|
||||
</span><span id="read_sfcf-275"><a href="#read_sfcf-275"><span class="linenos">275</span></a> <span class="c1"># here one can also implement the chekc from above.</span>
|
||||
</span><span id="read_sfcf-276"><a href="#read_sfcf-276"><span class="linenos">276</span></a> <span class="k">for</span> <span class="n">k</span><span class="p">,</span> <span class="n">line</span> <span class="ow">in</span> <span class="nb">enumerate</span><span class="p">(</span><span class="n">fp</span><span class="p">):</span>
|
||||
</span><span id="read_sfcf-277"><a href="#read_sfcf-277"><span class="linenos">277</span></a> <span class="k">if</span><span class="p">(</span><span class="n">k</span> <span class="o">>=</span> <span class="n">start_read</span> <span class="ow">and</span> <span class="n">k</span> <span class="o"><</span> <span class="n">start_read</span> <span class="o">+</span> <span class="n">T</span><span class="p">):</span>
|
||||
</span><span id="read_sfcf-277"><a href="#read_sfcf-277"><span class="linenos">277</span></a> <span class="k">if</span> <span class="p">(</span><span class="n">k</span> <span class="o">>=</span> <span class="n">start_read</span> <span class="ow">and</span> <span class="n">k</span> <span class="o"><</span> <span class="n">start_read</span> <span class="o">+</span> <span class="n">T</span><span class="p">):</span>
|
||||
</span><span id="read_sfcf-278"><a href="#read_sfcf-278"><span class="linenos">278</span></a> <span class="n">floats</span> <span class="o">=</span> <span class="nb">list</span><span class="p">(</span><span class="nb">map</span><span class="p">(</span><span class="nb">float</span><span class="p">,</span> <span class="n">line</span><span class="o">.</span><span class="n">split</span><span class="p">()))</span>
|
||||
</span><span id="read_sfcf-279"><a href="#read_sfcf-279"><span class="linenos">279</span></a> <span class="k">if</span> <span class="n">version</span> <span class="o">==</span> <span class="s2">"0.0"</span><span class="p">:</span>
|
||||
</span><span id="read_sfcf-280"><a href="#read_sfcf-280"><span class="linenos">280</span></a> <span class="n">deltas</span><span class="p">[</span><span class="n">k</span> <span class="o">-</span> <span class="n">start</span><span class="p">][</span><span class="n">i</span><span class="p">][</span><span class="n">cnfg</span><span class="p">]</span> <span class="o">=</span> <span class="n">floats</span><span class="p">[</span><span class="n">im</span> <span class="o">-</span> <span class="n">single</span><span class="p">]</span>
|
||||
|
|
|
|||
Loading…
Add table
Add a link
Reference in a new issue