<div dir="ltr">Thank you.<div><br></div><div>For a mere 10 timesteps on a 3 by 3 system,</div><div>the time for 'spinless fermions' model is ~1s</div><div>while the time for 'fermion Hubbard' model is ~475s</div><div>Is this within your expectation? The source codes are attached</div><div><br></div><div>I tried changing 'mpo_compression' to 'singlesite' but the result becomes incorrect.</div><div><br></div><div>Sorry for keep bothering you. I just want to get the most out of this excellent library.</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sun, Jul 9, 2017 at 4:26 AM, Michele Dolfi <span dir="ltr"><<a href="mailto:dolfim@phys.ethz.ch" target="_blank">dolfim@phys.ethz.ch</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word"><div>Well, there are multiple factors that makes it more computationally expensive, but not really thousand times…</div><div><br></div><div>The relevant factor are:</div><div>- hamiltonian terms are applied on two-sites, i.e. each will cost M^2*d^2. when you change from spinless d=2 to spinfull d=4, you are at least 4 times slower</div><div>- because of the fermionic sign, each term in the trotter decomposition requires multiple operations</div><div>- the compression to keep the MPS of a controllable dimension is done using an iterative scheme. depending on how strong you are compressing, it might take more iterations</div><div><br></div><div>The algorithm used is described in section 2.3 of our paper <a href="https://doi.org/10.1016/j.cpc.2014.08.019" target="_blank">https://doi.org/10.1016/<wbr>j.cpc.2014.08.019</a>, or more a bit more extensively in section 2.5 of my thesis <a href="http://dx.doi.org/10.3929/ethz-a-010735455" target="_blank">http://dx.doi.org/10.3929/<wbr>ethz-a-010735455</a>.</div><div><br></div><div><br></div><div>Best,</div><div>Michele</div><div><br></div><div><br></div><br><div><blockquote type="cite"><div><div class="h5"><div>On 8 Jul 2017, at 21:39, H <<a href="mailto:w.boson@gmail.com" target="_blank">w.boson@gmail.com</a>> wrote:</div><br class="m_7157877933521856715Apple-interchange-newline"></div></div><div><div><div class="h5"><div dir="ltr">Thank you. I have indeed overlooked this important parameter.<div><br></div><div>It now works smoothly with spinless fermions or bosons. However for spinful fermions the time it takes is more than a thousand times (or more) longer. Even for tiny 2D systems the time it takes is huge. Is it a known limitation and is there a physical reason?</div><div><br></div><div><br></div><div class="gmail_extra"><br><div class="gmail_quote">On Sat, Jul 8, 2017 at 2:48 AM, Michele Dolfi <span dir="ltr"><<a href="mailto:dolfim@phys.ethz.ch" target="_blank">dolfim@phys.ethz.ch</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word"><div>Note that by default mps_evolve only works for nearest neighbors lattices.</div><div>In higher dimensions this is not the case, since the lattice is unrolled into a one dimensional chain with longer range interactions.</div><div><br></div><div>We have a second flavor of the time evolution algorithm for these cases, you can use it via the parameter:</div><div>te_type=mpo</div><div>Quite some improvements to this variant was introduced in ALPS 2.3, so I suggest you to use that version.</div><div><br></div><div><br></div><div>Best,</div><div>Michele</div><div><br></div></div></blockquote></div></div></div>
<br><br></div></div><span class="">----<br>Comp-phys-alps-users Mailing List for the ALPS Project<br><a href="http://alps.comp-phys.org/" target="_blank">http://alps.comp-phys.org/</a><br><br>List info: <a href="https://lists.phys.ethz.ch//listinfo/comp-phys-alps-users" target="_blank">https://lists.phys.ethz.ch//<wbr>listinfo/comp-phys-alps-users</a><br>Archive: <a href="https://lists.phys.ethz.ch//pipermail/comp-phys-alps-users" target="_blank">https://lists.phys.ethz.ch//<wbr>pipermail/comp-phys-alps-users</a><br><br>Unsubscribe by writing a mail to <a href="mailto:comp-phys-alps-users-leave@lists.phys.ethz.ch" target="_blank">comp-phys-alps-users-leave@<wbr>lists.phys.ethz.ch</a>.</span></div></blockquote></div><br></div><br><br>
<br>
----<br>
Comp-phys-alps-users Mailing List for the ALPS Project<br>
<a href="http://alps.comp-phys.org/" rel="noreferrer" target="_blank">http://alps.comp-phys.org/</a><br>
<br>
List info: <a href="https://lists.phys.ethz.ch//listinfo/comp-phys-alps-users" rel="noreferrer" target="_blank">https://lists.phys.ethz.ch//<wbr>listinfo/comp-phys-alps-users</a><br>
Archive: <a href="https://lists.phys.ethz.ch//pipermail/comp-phys-alps-users" rel="noreferrer" target="_blank">https://lists.phys.ethz.ch//<wbr>pipermail/comp-phys-alps-users</a><br>
<br>
Unsubscribe by writing a mail to <a href="mailto:comp-phys-alps-users-leave@lists.phys.ethz.ch">comp-phys-alps-users-leave@<wbr>lists.phys.ethz.ch</a>.<br></blockquote></div><br></div>