<div dir="ltr">All evaluation scripts of that paper are available in the supplemental material at <a href="http://dx.doi.org/10.5281/zenodo.31003">http://dx.doi.org/10.5281/zenodo.31003</a>.<div>Section III.C explains how the data is extrapolated for energy measurements or other observables. For the local density and correlations we used a second order polynomial, some of the extrapolations are shown in the panels in Fig. 6.<div><br></div><div><br></div><div>Michele<br><div><br></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Jul 13, 2017 at 11:37 AM, Robertson Esperanza <span dir="ltr"><<a href="mailto:robbie.esperanza@gmail.com" target="_blank">robbie.esperanza@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><span style="font-size:12.8px">I have seen that you have extrapolated to obtain quantities as M approaches infinity. I would like to know the extrapolation method used (especially in Figures 3 and 4).</span><div><br style="font-size:12.8px"><span style="font-size:12.8px">Regards,</span><div style="font-size:12.8px">Robertson Esperanza</div></div></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Sun, Mar 26, 2017 at 8:33 PM, 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">We have some example in our paper at <a href="https://doi.org/10.1103/PhysRevB.92.195139" target="_blank">https://doi.org/10.1103/Phy<wbr>sRevB.92.195139</a>.<div><br></div><div>Basically you should always extrapolate in the bond dimension.</div><div><br></div><div><br></div><div>Best,</div><div>Michele</div><div><div class="m_-8900819910704678448h5"><div><br></div><div><br><div><blockquote type="cite"><div>On 26 Mar 2017, at 11:36, Robertson Esperanza <<a href="mailto:robbie.esperanza@gmail.com" target="_blank">robbie.esperanza@gmail.com</a>> wrote:</div><br class="m_-8900819910704678448m_-6932479936982201518Apple-interchange-newline"><div><div dir="ltr">Now it makes sense, thank you for the clarification. You said earlier that I should <span style="font-size:12.8px">address the convergence in the bond dimension separately, how is it done?</span></div><div class="gmail_extra"><br><div class="gmail_quote">On Sun, Mar 26, 2017 at 3:50 PM, 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>I suppose you are just printing the values of the first sweep, then the math is as follow:</div><div><br></div><div>N=96*2-1 (number of bonds that gets optimized)</div><div>since one sweep is left-to-right and right-to-left: 2*N=382</div><div><br></div><div><br></div><br><div><blockquote type="cite"><div><div class="m_-8900819910704678448m_-6932479936982201518h5"><div>On 26 Mar 2017, at 01:45, Robertson Esperanza <<a href="mailto:robbie.esperanza@gmail.com" target="_blank">robbie.esperanza@gmail.com</a>> wrote:</div><br class="m_-8900819910704678448m_-6932479936982201518m_2338771395665549078Apple-interchange-newline"></div></div><div><div><div class="m_-8900819910704678448m_-6932479936982201518h5"><div dir="ltr">To clarify, here is a text file showing an iteration ('iteration': 0) showing 'observable': 'TruncatedWeight':<br><br>from here, it shows 382 values for that observable. I'm wondering why there are 382 truncation error values.</div><div class="gmail_extra"><br><div class="gmail_quote">On Sun, Mar 26, 2017 at 12:50 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">First note that this refers only to convergence of a single simulation. One should then address the convergence in the bond dimension separately.<br>
<br>
I would say that there is no exact rule for convergence, but I will give you some example of what we usually check.<br>
- Convergence of energy.<br>
- Truncation error. This is the sum of all discarded eigenvalues in the reduced density matrix.<br>
- Any other observable of interest, for example the local density is good indicator, because sometime you fall in a state which is not symmetric and you are pretty sure that something wrong happened.<br>
For this values the iteration values we store one number per optimization, i.e. one sweep are 2*N with N the system size.<br>
<br>
It is a bit arbitrary how to aggregate the truncation errors. It should just be clear outlined in the publications, so that readers know what you are talking about. Usually you either sum them of take the maximum value.<br>
<br>
<br>
Regards,<br>
Michele<br>
<div><div class="m_-8900819910704678448m_-6932479936982201518m_2338771395665549078h5"><br>
<br>
> On 24 Mar 2017, at 22:24, Robertson Esperanza <<a href="mailto:robbie.esperanza@gmail.com" target="_blank">robbie.esperanza@gmail.com</a>> wrote:<br>
><br>
> Thank you for all the help so far, now I'm trying to check for convergence of my results.<br>
><br>
> First, is it enough to just check how the truncation error changes with sweeps (at a certain MAXSTATE value) or I need to both check the iteration history of the ground state energy and the truncation error?<br>
><br>
> Second, checking the supplementary materials provided (Computer Physics Communications 185 (2014) 3430–3440) for the Hubbard ladder, why are there 381 values of energy and truncation weights for every iteration? Why is the total truncation error for that iteration is the sum of all the 381 truncation weight values (based on iteration.py)? How is truncation error defined (e.g. is it 1 - sum of all reduced density eigenvalues corresponding to all the picked states?) ?<br>
><br>
>  As far as I understand the DMRG algorithm every iteration/sweep has a corresponding bond dimension in which as you progress through iterations (until you reach the "SWEEPS"th iteration), bond dimension increases until you reach "MAXSTATES" at the "SWEEPS"th iteration.<br>
><br>
> Regards,<br>
> Robertson Esperanza<br>
><br>
><br>
</div></div>> ----<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//li<wbr>stinfo/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//pi<wbr>permail/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@lis<wbr>ts.phys.ethz.ch</a>.<br>
<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//li<wbr>stinfo/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//pi<wbr>permail/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@lis<wbr>ts.phys.ethz.ch</a>.</blockquote></div><br></div>
</div></div><span id="m_-8900819910704678448m_-6932479936982201518m_2338771395665549078cid:E7D0604D-FF9D-46F4-AC82-4E15332DBE62@ADS.NKSA.CH"><sample_result.rtf></span><span><br><br>----<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//li<wbr>stinfo/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//pi<wbr>permail/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@lis<wbr>ts.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//li<wbr>stinfo/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//pi<wbr>permail/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@lis<wbr>ts.phys.ethz.ch</a>.<br></blockquote></div><br></div>
<br><br>----<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//li<wbr>stinfo/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//pi<wbr>permail/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@lis<wbr>ts.phys.ethz.ch</a>.</div></blockquote></div><br></div></div></div></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//li<wbr>stinfo/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//pi<wbr>permail/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@lis<wbr>ts.phys.ethz.ch</a>.<br></blockquote></div><br></div>
</div></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>