mirror of
https://github.com/fjosw/pyerrors.git
synced 2025-03-16 07:10:24 +01:00
Documentation updated
This commit is contained in:
fjosw
parent
53cfbecb59
commit
ac6b711bec
1 changed files with 30 additions and 30 deletions
|
|
@ -471,16 +471,16 @@
|
|||
</span><span id="L-258"><a href="#L-258"><span class="linenos"> 258</span></a><span class="sd"> "Eigenvector" - Use the method described in arXiv:2004.10472 [hep-lat] to find the set of v(t) belonging to the state.</span>
|
||||
</span><span id="L-259"><a href="#L-259"><span class="linenos"> 259</span></a><span class="sd"> The reference state is identified by its eigenvalue at t=ts</span>
|
||||
</span><span id="L-260"><a href="#L-260"><span class="linenos"> 260</span></a><span class="sd"> """</span>
|
||||
</span><span id="L-261"><a href="#L-261"><span class="linenos"> 261</span></a> <span class="k">if</span> <span class="n">sorted_list</span> <span class="ow">is</span> <span class="kc">None</span><span class="p">:</span>
|
||||
</span><span id="L-262"><a href="#L-262"><span class="linenos"> 262</span></a> <span class="k">if</span> <span class="p">(</span><span class="n">ts</span> <span class="ow">is</span> <span class="kc">None</span><span class="p">):</span>
|
||||
</span><span id="L-263"><a href="#L-263"><span class="linenos"> 263</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"ts is required if sorted_list=None"</span><span class="p">)</span>
|
||||
</span><span id="L-264"><a href="#L-264"><span class="linenos"> 264</span></a> <span class="k">if</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">t0</span><span class="p">]</span> <span class="ow">is</span> <span class="kc">None</span><span class="p">)</span> <span class="ow">or</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">ts</span><span class="p">]</span> <span class="ow">is</span> <span class="kc">None</span><span class="p">):</span>
|
||||
</span><span id="L-265"><a href="#L-265"><span class="linenos"> 265</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"Corr not defined at t0/ts"</span><span class="p">)</span>
|
||||
</span><span id="L-266"><a href="#L-266"><span class="linenos"> 266</span></a> <span class="n">G0</span><span class="p">,</span> <span class="n">Gt</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">empty</span><span class="p">([</span><span class="bp">self</span><span class="o">.</span><span class="n">N</span><span class="p">,</span> <span class="bp">self</span><span class="o">.</span><span class="n">N</span><span class="p">],</span> <span class="n">dtype</span><span class="o">=</span><span class="s2">"double"</span><span class="p">),</span> <span class="n">np</span><span class="o">.</span><span class="n">empty</span><span class="p">([</span><span class="bp">self</span><span class="o">.</span><span class="n">N</span><span class="p">,</span> <span class="bp">self</span><span class="o">.</span><span class="n">N</span><span class="p">],</span> <span class="n">dtype</span><span class="o">=</span><span class="s2">"double"</span><span class="p">)</span>
|
||||
</span><span id="L-267"><a href="#L-267"><span class="linenos"> 267</span></a> <span class="n">symmetric_corr</span> <span class="o">=</span> <span class="bp">self</span><span class="o">.</span><span class="n">matrix_symmetric</span><span class="p">()</span>
|
||||
</span><span id="L-261"><a href="#L-261"><span class="linenos"> 261</span></a> <span class="n">symmetric_corr</span> <span class="o">=</span> <span class="bp">self</span><span class="o">.</span><span class="n">matrix_symmetric</span><span class="p">()</span>
|
||||
</span><span id="L-262"><a href="#L-262"><span class="linenos"> 262</span></a> <span class="k">if</span> <span class="n">sorted_list</span> <span class="ow">is</span> <span class="kc">None</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="p">(</span><span class="n">ts</span> <span class="ow">is</span> <span class="kc">None</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="s2">"ts is required if sorted_list=None"</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="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t0</span><span class="p">]</span> <span class="ow">is</span> <span class="kc">None</span><span class="p">)</span> <span class="ow">or</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">ts</span><span class="p">]</span> <span class="ow">is</span> <span class="kc">None</span><span class="p">):</span>
|
||||
</span><span id="L-266"><a href="#L-266"><span class="linenos"> 266</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"Corr not defined at t0/ts"</span><span class="p">)</span>
|
||||
</span><span id="L-267"><a href="#L-267"><span class="linenos"> 267</span></a> <span class="n">G0</span><span class="p">,</span> <span class="n">Gt</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">empty</span><span class="p">([</span><span class="bp">self</span><span class="o">.</span><span class="n">N</span><span class="p">,</span> <span class="bp">self</span><span class="o">.</span><span class="n">N</span><span class="p">],</span> <span class="n">dtype</span><span class="o">=</span><span class="s2">"double"</span><span class="p">),</span> <span class="n">np</span><span class="o">.</span><span class="n">empty</span><span class="p">([</span><span class="bp">self</span><span class="o">.</span><span class="n">N</span><span class="p">,</span> <span class="bp">self</span><span class="o">.</span><span class="n">N</span><span class="p">],</span> <span class="n">dtype</span><span class="o">=</span><span class="s2">"double"</span><span class="p">)</span>
|
||||
</span><span id="L-268"><a href="#L-268"><span class="linenos"> 268</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="bp">self</span><span class="o">.</span><span class="n">N</span><span class="p">):</span>
|
||||
</span><span id="L-269"><a href="#L-269"><span class="linenos"> 269</span></a> <span class="k">for</span> <span class="n">j</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">N</span><span class="p">):</span>
|
||||
</span><span id="L-270"><a href="#L-270"><span class="linenos"> 270</span></a> <span class="n">G0</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="o">=</span> <span class="n">symmetric_corr</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t0</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="o">.</span><span class="n">value</span>
|
||||
</span><span id="L-270"><a href="#L-270"><span class="linenos"> 270</span></a> <span class="n">G0</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="o">=</span> <span class="n">symmetric_corr</span><span class="p">[</span><span class="n">t0</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="o">.</span><span class="n">value</span>
|
||||
</span><span id="L-271"><a href="#L-271"><span class="linenos"> 271</span></a> <span class="n">Gt</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="o">=</span> <span class="n">symmetric_corr</span><span class="p">[</span><span class="n">ts</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="o">.</span><span class="n">value</span>
|
||||
</span><span id="L-272"><a href="#L-272"><span class="linenos"> 272</span></a>
|
||||
</span><span id="L-273"><a href="#L-273"><span class="linenos"> 273</span></a> <span class="n">sp_vecs</span> <span class="o">=</span> <span class="n">_GEVP_solver</span><span class="p">(</span><span class="n">Gt</span><span class="p">,</span> <span class="n">G0</span><span class="p">)</span>
|
||||
|
|
@ -493,8 +493,8 @@
|
|||
</span><span id="L-280"><a href="#L-280"><span class="linenos"> 280</span></a> <span class="n">G0</span><span class="p">,</span> <span class="n">Gt</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">empty</span><span class="p">([</span><span class="bp">self</span><span class="o">.</span><span class="n">N</span><span class="p">,</span> <span class="bp">self</span><span class="o">.</span><span class="n">N</span><span class="p">],</span> <span class="n">dtype</span><span class="o">=</span><span class="s2">"double"</span><span class="p">),</span> <span class="n">np</span><span class="o">.</span><span class="n">empty</span><span class="p">([</span><span class="bp">self</span><span class="o">.</span><span class="n">N</span><span class="p">,</span> <span class="bp">self</span><span class="o">.</span><span class="n">N</span><span class="p">],</span> <span class="n">dtype</span><span class="o">=</span><span class="s2">"double"</span><span class="p">)</span>
|
||||
</span><span id="L-281"><a href="#L-281"><span class="linenos"> 281</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="bp">self</span><span class="o">.</span><span class="n">N</span><span class="p">):</span>
|
||||
</span><span id="L-282"><a href="#L-282"><span class="linenos"> 282</span></a> <span class="k">for</span> <span class="n">j</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">N</span><span class="p">):</span>
|
||||
</span><span id="L-283"><a href="#L-283"><span class="linenos"> 283</span></a> <span class="n">G0</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="o">=</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t0</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="o">.</span><span class="n">value</span>
|
||||
</span><span id="L-284"><a href="#L-284"><span class="linenos"> 284</span></a> <span class="n">Gt</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="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="n">i</span><span class="p">,</span> <span class="n">j</span><span class="p">]</span><span class="o">.</span><span class="n">value</span>
|
||||
</span><span id="L-283"><a href="#L-283"><span class="linenos"> 283</span></a> <span class="n">G0</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="o">=</span> <span class="n">symmetric_corr</span><span class="p">[</span><span class="n">t0</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="o">.</span><span class="n">value</span>
|
||||
</span><span id="L-284"><a href="#L-284"><span class="linenos"> 284</span></a> <span class="n">Gt</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="o">=</span> <span class="n">symmetric_corr</span><span class="p">[</span><span class="n">t</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="o">.</span><span class="n">value</span>
|
||||
</span><span id="L-285"><a href="#L-285"><span class="linenos"> 285</span></a>
|
||||
</span><span id="L-286"><a href="#L-286"><span class="linenos"> 286</span></a> <span class="n">sp_vecs</span> <span class="o">=</span> <span class="n">_GEVP_solver</span><span class="p">(</span><span class="n">Gt</span><span class="p">,</span> <span class="n">G0</span><span class="p">)</span>
|
||||
</span><span id="L-287"><a href="#L-287"><span class="linenos"> 287</span></a> <span class="k">if</span> <span class="n">sorted_list</span> <span class="o">==</span> <span class="s2">"Eigenvalue"</span><span class="p">:</span>
|
||||
|
|
@ -1693,16 +1693,16 @@
|
|||
</span><span id="Corr-259"><a href="#Corr-259"><span class="linenos"> 259</span></a><span class="sd"> "Eigenvector" - Use the method described in arXiv:2004.10472 [hep-lat] to find the set of v(t) belonging to the state.</span>
|
||||
</span><span id="Corr-260"><a href="#Corr-260"><span class="linenos"> 260</span></a><span class="sd"> The reference state is identified by its eigenvalue at t=ts</span>
|
||||
</span><span id="Corr-261"><a href="#Corr-261"><span class="linenos"> 261</span></a><span class="sd"> """</span>
|
||||
</span><span id="Corr-262"><a href="#Corr-262"><span class="linenos"> 262</span></a> <span class="k">if</span> <span class="n">sorted_list</span> <span class="ow">is</span> <span class="kc">None</span><span class="p">:</span>
|
||||
</span><span id="Corr-263"><a href="#Corr-263"><span class="linenos"> 263</span></a> <span class="k">if</span> <span class="p">(</span><span class="n">ts</span> <span class="ow">is</span> <span class="kc">None</span><span class="p">):</span>
|
||||
</span><span id="Corr-264"><a href="#Corr-264"><span class="linenos"> 264</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"ts is required if sorted_list=None"</span><span class="p">)</span>
|
||||
</span><span id="Corr-265"><a href="#Corr-265"><span class="linenos"> 265</span></a> <span class="k">if</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">t0</span><span class="p">]</span> <span class="ow">is</span> <span class="kc">None</span><span class="p">)</span> <span class="ow">or</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">ts</span><span class="p">]</span> <span class="ow">is</span> <span class="kc">None</span><span class="p">):</span>
|
||||
</span><span id="Corr-266"><a href="#Corr-266"><span class="linenos"> 266</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"Corr not defined at t0/ts"</span><span class="p">)</span>
|
||||
</span><span id="Corr-267"><a href="#Corr-267"><span class="linenos"> 267</span></a> <span class="n">G0</span><span class="p">,</span> <span class="n">Gt</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">empty</span><span class="p">([</span><span class="bp">self</span><span class="o">.</span><span class="n">N</span><span class="p">,</span> <span class="bp">self</span><span class="o">.</span><span class="n">N</span><span class="p">],</span> <span class="n">dtype</span><span class="o">=</span><span class="s2">"double"</span><span class="p">),</span> <span class="n">np</span><span class="o">.</span><span class="n">empty</span><span class="p">([</span><span class="bp">self</span><span class="o">.</span><span class="n">N</span><span class="p">,</span> <span class="bp">self</span><span class="o">.</span><span class="n">N</span><span class="p">],</span> <span class="n">dtype</span><span class="o">=</span><span class="s2">"double"</span><span class="p">)</span>
|
||||
</span><span id="Corr-268"><a href="#Corr-268"><span class="linenos"> 268</span></a> <span class="n">symmetric_corr</span> <span class="o">=</span> <span class="bp">self</span><span class="o">.</span><span class="n">matrix_symmetric</span><span class="p">()</span>
|
||||
</span><span id="Corr-262"><a href="#Corr-262"><span class="linenos"> 262</span></a> <span class="n">symmetric_corr</span> <span class="o">=</span> <span class="bp">self</span><span class="o">.</span><span class="n">matrix_symmetric</span><span class="p">()</span>
|
||||
</span><span id="Corr-263"><a href="#Corr-263"><span class="linenos"> 263</span></a> <span class="k">if</span> <span class="n">sorted_list</span> <span class="ow">is</span> <span class="kc">None</span><span class="p">:</span>
|
||||
</span><span id="Corr-264"><a href="#Corr-264"><span class="linenos"> 264</span></a> <span class="k">if</span> <span class="p">(</span><span class="n">ts</span> <span class="ow">is</span> <span class="kc">None</span><span class="p">):</span>
|
||||
</span><span id="Corr-265"><a href="#Corr-265"><span class="linenos"> 265</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"ts is required if sorted_list=None"</span><span class="p">)</span>
|
||||
</span><span id="Corr-266"><a href="#Corr-266"><span class="linenos"> 266</span></a> <span class="k">if</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">t0</span><span class="p">]</span> <span class="ow">is</span> <span class="kc">None</span><span class="p">)</span> <span class="ow">or</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">ts</span><span class="p">]</span> <span class="ow">is</span> <span class="kc">None</span><span class="p">):</span>
|
||||
</span><span id="Corr-267"><a href="#Corr-267"><span class="linenos"> 267</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"Corr not defined at t0/ts"</span><span class="p">)</span>
|
||||
</span><span id="Corr-268"><a href="#Corr-268"><span class="linenos"> 268</span></a> <span class="n">G0</span><span class="p">,</span> <span class="n">Gt</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">empty</span><span class="p">([</span><span class="bp">self</span><span class="o">.</span><span class="n">N</span><span class="p">,</span> <span class="bp">self</span><span class="o">.</span><span class="n">N</span><span class="p">],</span> <span class="n">dtype</span><span class="o">=</span><span class="s2">"double"</span><span class="p">),</span> <span class="n">np</span><span class="o">.</span><span class="n">empty</span><span class="p">([</span><span class="bp">self</span><span class="o">.</span><span class="n">N</span><span class="p">,</span> <span class="bp">self</span><span class="o">.</span><span class="n">N</span><span class="p">],</span> <span class="n">dtype</span><span class="o">=</span><span class="s2">"double"</span><span class="p">)</span>
|
||||
</span><span id="Corr-269"><a href="#Corr-269"><span class="linenos"> 269</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="bp">self</span><span class="o">.</span><span class="n">N</span><span class="p">):</span>
|
||||
</span><span id="Corr-270"><a href="#Corr-270"><span class="linenos"> 270</span></a> <span class="k">for</span> <span class="n">j</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">N</span><span class="p">):</span>
|
||||
</span><span id="Corr-271"><a href="#Corr-271"><span class="linenos"> 271</span></a> <span class="n">G0</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="o">=</span> <span class="n">symmetric_corr</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t0</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="o">.</span><span class="n">value</span>
|
||||
</span><span id="Corr-271"><a href="#Corr-271"><span class="linenos"> 271</span></a> <span class="n">G0</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="o">=</span> <span class="n">symmetric_corr</span><span class="p">[</span><span class="n">t0</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="o">.</span><span class="n">value</span>
|
||||
</span><span id="Corr-272"><a href="#Corr-272"><span class="linenos"> 272</span></a> <span class="n">Gt</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="o">=</span> <span class="n">symmetric_corr</span><span class="p">[</span><span class="n">ts</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="o">.</span><span class="n">value</span>
|
||||
</span><span id="Corr-273"><a href="#Corr-273"><span class="linenos"> 273</span></a>
|
||||
</span><span id="Corr-274"><a href="#Corr-274"><span class="linenos"> 274</span></a> <span class="n">sp_vecs</span> <span class="o">=</span> <span class="n">_GEVP_solver</span><span class="p">(</span><span class="n">Gt</span><span class="p">,</span> <span class="n">G0</span><span class="p">)</span>
|
||||
|
|
@ -1715,8 +1715,8 @@
|
|||
</span><span id="Corr-281"><a href="#Corr-281"><span class="linenos"> 281</span></a> <span class="n">G0</span><span class="p">,</span> <span class="n">Gt</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">empty</span><span class="p">([</span><span class="bp">self</span><span class="o">.</span><span class="n">N</span><span class="p">,</span> <span class="bp">self</span><span class="o">.</span><span class="n">N</span><span class="p">],</span> <span class="n">dtype</span><span class="o">=</span><span class="s2">"double"</span><span class="p">),</span> <span class="n">np</span><span class="o">.</span><span class="n">empty</span><span class="p">([</span><span class="bp">self</span><span class="o">.</span><span class="n">N</span><span class="p">,</span> <span class="bp">self</span><span class="o">.</span><span class="n">N</span><span class="p">],</span> <span class="n">dtype</span><span class="o">=</span><span class="s2">"double"</span><span class="p">)</span>
|
||||
</span><span id="Corr-282"><a href="#Corr-282"><span class="linenos"> 282</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="bp">self</span><span class="o">.</span><span class="n">N</span><span class="p">):</span>
|
||||
</span><span id="Corr-283"><a href="#Corr-283"><span class="linenos"> 283</span></a> <span class="k">for</span> <span class="n">j</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">N</span><span class="p">):</span>
|
||||
</span><span id="Corr-284"><a href="#Corr-284"><span class="linenos"> 284</span></a> <span class="n">G0</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="o">=</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t0</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="o">.</span><span class="n">value</span>
|
||||
</span><span id="Corr-285"><a href="#Corr-285"><span class="linenos"> 285</span></a> <span class="n">Gt</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="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="n">i</span><span class="p">,</span> <span class="n">j</span><span class="p">]</span><span class="o">.</span><span class="n">value</span>
|
||||
</span><span id="Corr-284"><a href="#Corr-284"><span class="linenos"> 284</span></a> <span class="n">G0</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="o">=</span> <span class="n">symmetric_corr</span><span class="p">[</span><span class="n">t0</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="o">.</span><span class="n">value</span>
|
||||
</span><span id="Corr-285"><a href="#Corr-285"><span class="linenos"> 285</span></a> <span class="n">Gt</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="o">=</span> <span class="n">symmetric_corr</span><span class="p">[</span><span class="n">t</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="o">.</span><span class="n">value</span>
|
||||
</span><span id="Corr-286"><a href="#Corr-286"><span class="linenos"> 286</span></a>
|
||||
</span><span id="Corr-287"><a href="#Corr-287"><span class="linenos"> 287</span></a> <span class="n">sp_vecs</span> <span class="o">=</span> <span class="n">_GEVP_solver</span><span class="p">(</span><span class="n">Gt</span><span class="p">,</span> <span class="n">G0</span><span class="p">)</span>
|
||||
</span><span id="Corr-288"><a href="#Corr-288"><span class="linenos"> 288</span></a> <span class="k">if</span> <span class="n">sorted_list</span> <span class="o">==</span> <span class="s2">"Eigenvalue"</span><span class="p">:</span>
|
||||
|
|
@ -3042,16 +3042,16 @@ timeslice and the error on each timeslice.</p>
|
|||
</span><span id="Corr.GEVP-259"><a href="#Corr.GEVP-259"><span class="linenos">259</span></a><span class="sd"> "Eigenvector" - Use the method described in arXiv:2004.10472 [hep-lat] to find the set of v(t) belonging to the state.</span>
|
||||
</span><span id="Corr.GEVP-260"><a href="#Corr.GEVP-260"><span class="linenos">260</span></a><span class="sd"> The reference state is identified by its eigenvalue at t=ts</span>
|
||||
</span><span id="Corr.GEVP-261"><a href="#Corr.GEVP-261"><span class="linenos">261</span></a><span class="sd"> """</span>
|
||||
</span><span id="Corr.GEVP-262"><a href="#Corr.GEVP-262"><span class="linenos">262</span></a> <span class="k">if</span> <span class="n">sorted_list</span> <span class="ow">is</span> <span class="kc">None</span><span class="p">:</span>
|
||||
</span><span id="Corr.GEVP-263"><a href="#Corr.GEVP-263"><span class="linenos">263</span></a> <span class="k">if</span> <span class="p">(</span><span class="n">ts</span> <span class="ow">is</span> <span class="kc">None</span><span class="p">):</span>
|
||||
</span><span id="Corr.GEVP-264"><a href="#Corr.GEVP-264"><span class="linenos">264</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"ts is required if sorted_list=None"</span><span class="p">)</span>
|
||||
</span><span id="Corr.GEVP-265"><a href="#Corr.GEVP-265"><span class="linenos">265</span></a> <span class="k">if</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">t0</span><span class="p">]</span> <span class="ow">is</span> <span class="kc">None</span><span class="p">)</span> <span class="ow">or</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">ts</span><span class="p">]</span> <span class="ow">is</span> <span class="kc">None</span><span class="p">):</span>
|
||||
</span><span id="Corr.GEVP-266"><a href="#Corr.GEVP-266"><span class="linenos">266</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"Corr not defined at t0/ts"</span><span class="p">)</span>
|
||||
</span><span id="Corr.GEVP-267"><a href="#Corr.GEVP-267"><span class="linenos">267</span></a> <span class="n">G0</span><span class="p">,</span> <span class="n">Gt</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">empty</span><span class="p">([</span><span class="bp">self</span><span class="o">.</span><span class="n">N</span><span class="p">,</span> <span class="bp">self</span><span class="o">.</span><span class="n">N</span><span class="p">],</span> <span class="n">dtype</span><span class="o">=</span><span class="s2">"double"</span><span class="p">),</span> <span class="n">np</span><span class="o">.</span><span class="n">empty</span><span class="p">([</span><span class="bp">self</span><span class="o">.</span><span class="n">N</span><span class="p">,</span> <span class="bp">self</span><span class="o">.</span><span class="n">N</span><span class="p">],</span> <span class="n">dtype</span><span class="o">=</span><span class="s2">"double"</span><span class="p">)</span>
|
||||
</span><span id="Corr.GEVP-268"><a href="#Corr.GEVP-268"><span class="linenos">268</span></a> <span class="n">symmetric_corr</span> <span class="o">=</span> <span class="bp">self</span><span class="o">.</span><span class="n">matrix_symmetric</span><span class="p">()</span>
|
||||
</span><span id="Corr.GEVP-262"><a href="#Corr.GEVP-262"><span class="linenos">262</span></a> <span class="n">symmetric_corr</span> <span class="o">=</span> <span class="bp">self</span><span class="o">.</span><span class="n">matrix_symmetric</span><span class="p">()</span>
|
||||
</span><span id="Corr.GEVP-263"><a href="#Corr.GEVP-263"><span class="linenos">263</span></a> <span class="k">if</span> <span class="n">sorted_list</span> <span class="ow">is</span> <span class="kc">None</span><span class="p">:</span>
|
||||
</span><span id="Corr.GEVP-264"><a href="#Corr.GEVP-264"><span class="linenos">264</span></a> <span class="k">if</span> <span class="p">(</span><span class="n">ts</span> <span class="ow">is</span> <span class="kc">None</span><span class="p">):</span>
|
||||
</span><span id="Corr.GEVP-265"><a href="#Corr.GEVP-265"><span class="linenos">265</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"ts is required if sorted_list=None"</span><span class="p">)</span>
|
||||
</span><span id="Corr.GEVP-266"><a href="#Corr.GEVP-266"><span class="linenos">266</span></a> <span class="k">if</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">t0</span><span class="p">]</span> <span class="ow">is</span> <span class="kc">None</span><span class="p">)</span> <span class="ow">or</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">ts</span><span class="p">]</span> <span class="ow">is</span> <span class="kc">None</span><span class="p">):</span>
|
||||
</span><span id="Corr.GEVP-267"><a href="#Corr.GEVP-267"><span class="linenos">267</span></a> <span class="k">raise</span> <span class="ne">Exception</span><span class="p">(</span><span class="s2">"Corr not defined at t0/ts"</span><span class="p">)</span>
|
||||
</span><span id="Corr.GEVP-268"><a href="#Corr.GEVP-268"><span class="linenos">268</span></a> <span class="n">G0</span><span class="p">,</span> <span class="n">Gt</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">empty</span><span class="p">([</span><span class="bp">self</span><span class="o">.</span><span class="n">N</span><span class="p">,</span> <span class="bp">self</span><span class="o">.</span><span class="n">N</span><span class="p">],</span> <span class="n">dtype</span><span class="o">=</span><span class="s2">"double"</span><span class="p">),</span> <span class="n">np</span><span class="o">.</span><span class="n">empty</span><span class="p">([</span><span class="bp">self</span><span class="o">.</span><span class="n">N</span><span class="p">,</span> <span class="bp">self</span><span class="o">.</span><span class="n">N</span><span class="p">],</span> <span class="n">dtype</span><span class="o">=</span><span class="s2">"double"</span><span class="p">)</span>
|
||||
</span><span id="Corr.GEVP-269"><a href="#Corr.GEVP-269"><span class="linenos">269</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="bp">self</span><span class="o">.</span><span class="n">N</span><span class="p">):</span>
|
||||
</span><span id="Corr.GEVP-270"><a href="#Corr.GEVP-270"><span class="linenos">270</span></a> <span class="k">for</span> <span class="n">j</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">N</span><span class="p">):</span>
|
||||
</span><span id="Corr.GEVP-271"><a href="#Corr.GEVP-271"><span class="linenos">271</span></a> <span class="n">G0</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="o">=</span> <span class="n">symmetric_corr</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t0</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="o">.</span><span class="n">value</span>
|
||||
</span><span id="Corr.GEVP-271"><a href="#Corr.GEVP-271"><span class="linenos">271</span></a> <span class="n">G0</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="o">=</span> <span class="n">symmetric_corr</span><span class="p">[</span><span class="n">t0</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="o">.</span><span class="n">value</span>
|
||||
</span><span id="Corr.GEVP-272"><a href="#Corr.GEVP-272"><span class="linenos">272</span></a> <span class="n">Gt</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="o">=</span> <span class="n">symmetric_corr</span><span class="p">[</span><span class="n">ts</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="o">.</span><span class="n">value</span>
|
||||
</span><span id="Corr.GEVP-273"><a href="#Corr.GEVP-273"><span class="linenos">273</span></a>
|
||||
</span><span id="Corr.GEVP-274"><a href="#Corr.GEVP-274"><span class="linenos">274</span></a> <span class="n">sp_vecs</span> <span class="o">=</span> <span class="n">_GEVP_solver</span><span class="p">(</span><span class="n">Gt</span><span class="p">,</span> <span class="n">G0</span><span class="p">)</span>
|
||||
|
|
@ -3064,8 +3064,8 @@ timeslice and the error on each timeslice.</p>
|
|||
</span><span id="Corr.GEVP-281"><a href="#Corr.GEVP-281"><span class="linenos">281</span></a> <span class="n">G0</span><span class="p">,</span> <span class="n">Gt</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">empty</span><span class="p">([</span><span class="bp">self</span><span class="o">.</span><span class="n">N</span><span class="p">,</span> <span class="bp">self</span><span class="o">.</span><span class="n">N</span><span class="p">],</span> <span class="n">dtype</span><span class="o">=</span><span class="s2">"double"</span><span class="p">),</span> <span class="n">np</span><span class="o">.</span><span class="n">empty</span><span class="p">([</span><span class="bp">self</span><span class="o">.</span><span class="n">N</span><span class="p">,</span> <span class="bp">self</span><span class="o">.</span><span class="n">N</span><span class="p">],</span> <span class="n">dtype</span><span class="o">=</span><span class="s2">"double"</span><span class="p">)</span>
|
||||
</span><span id="Corr.GEVP-282"><a href="#Corr.GEVP-282"><span class="linenos">282</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="bp">self</span><span class="o">.</span><span class="n">N</span><span class="p">):</span>
|
||||
</span><span id="Corr.GEVP-283"><a href="#Corr.GEVP-283"><span class="linenos">283</span></a> <span class="k">for</span> <span class="n">j</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">N</span><span class="p">):</span>
|
||||
</span><span id="Corr.GEVP-284"><a href="#Corr.GEVP-284"><span class="linenos">284</span></a> <span class="n">G0</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="o">=</span> <span class="bp">self</span><span class="o">.</span><span class="n">content</span><span class="p">[</span><span class="n">t0</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="o">.</span><span class="n">value</span>
|
||||
</span><span id="Corr.GEVP-285"><a href="#Corr.GEVP-285"><span class="linenos">285</span></a> <span class="n">Gt</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="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="n">i</span><span class="p">,</span> <span class="n">j</span><span class="p">]</span><span class="o">.</span><span class="n">value</span>
|
||||
</span><span id="Corr.GEVP-284"><a href="#Corr.GEVP-284"><span class="linenos">284</span></a> <span class="n">G0</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="o">=</span> <span class="n">symmetric_corr</span><span class="p">[</span><span class="n">t0</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="o">.</span><span class="n">value</span>
|
||||
</span><span id="Corr.GEVP-285"><a href="#Corr.GEVP-285"><span class="linenos">285</span></a> <span class="n">Gt</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="o">=</span> <span class="n">symmetric_corr</span><span class="p">[</span><span class="n">t</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="o">.</span><span class="n">value</span>
|
||||
</span><span id="Corr.GEVP-286"><a href="#Corr.GEVP-286"><span class="linenos">286</span></a>
|
||||
</span><span id="Corr.GEVP-287"><a href="#Corr.GEVP-287"><span class="linenos">287</span></a> <span class="n">sp_vecs</span> <span class="o">=</span> <span class="n">_GEVP_solver</span><span class="p">(</span><span class="n">Gt</span><span class="p">,</span> <span class="n">G0</span><span class="p">)</span>
|
||||
</span><span id="Corr.GEVP-288"><a href="#Corr.GEVP-288"><span class="linenos">288</span></a> <span class="k">if</span> <span class="n">sorted_list</span> <span class="o">==</span> <span class="s2">"Eigenvalue"</span><span class="p">:</span>
|
||||
|
|
|
|||
Loading…
Add table
Reference in a new issue