<div dir="ltr">Regarding the BNS issues, is it something due to when you recover from checkpoint?<div><br></div><div>In Whisky we had an issue related to that due to Con2Prim running after recovery from checkpoint.</div><div><br></div><div>E.g., you run a simulation from iteration 0 to iteration 24 (made up numbers just to explain myself).</div><div>You try the same simulation, but this time you run from 0 to 12, checkpoint, recover from checkpoint and run from 12 to 24.</div><div>The results between  these two simulations at iteration 24 were slightly different</div><div><br></div><div>The issue was in the Newton Raphson being called by Con2Prim and using a different initial guess when recovering from checkpoint. During a run it uses the previous values of the primitive variables as an initial guess, but these are not available when restarting from checkpoint. The solution, as far as I remember, was not to run Con2Prim during recovery from checkpoint (since the primitive variables are anyway recovered from the checkpoint files).</div><div><br></div><div>Cheers,</div><div>Bruno</div><div><br></div></div><br><div class="gmail_quote gmail_quote_container"><div dir="ltr" class="gmail_attr">Il giorno gio 16 gen 2025 alle ore 19:59 Roland Haas <<a href="mailto:rhaas@illinois.edu">rhaas@illinois.edu</a>> ha scritto:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hello all,<br>
<br>
> # Reproducibility issues with BNS<br>
>   * Wolfgang Tichy has obtained different results with different MPI ranks on<br>
Sorry, my bad. Correction: Michal Pirog (who is working with Wolfgang)<br>
noticed these.<br>
<br>
Yours,<br>
Roland<br>
<br>
-- <br>
My email is as private as my paper mail. I therefore support encrypting<br>
and signing email messages. Get my PGP key from <a href="http://pgp.mit.edu" rel="noreferrer" target="_blank">http://pgp.mit.edu</a> .<br>
_______________________________________________<br>
Users mailing list<br>
<a href="mailto:Users@einsteintoolkit.org" target="_blank">Users@einsteintoolkit.org</a><br>
<a href="http://lists.einsteintoolkit.org/mailman/listinfo/users" rel="noreferrer" target="_blank">http://lists.einsteintoolkit.org/mailman/listinfo/users</a><br>
</blockquote></div><div><br clear="all"></div><div><br></div><span class="gmail_signature_prefix">-- </span><br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><p><font color="#000000">Pr<font face="arial, sans-serif">of. Bruno Giacomazzo<br>Deputy Director of the Department of Physics<br></font></font><span style="font-family:arial,sans-serif;color:rgb(0,0,0)">University of Milano-Bicocca<br></span><span style="font-family:arial,sans-serif;color:rgb(0,0,0)">Piazza della Scienza 3<br></span><span style="color:rgb(0,0,0)"><font face="arial, sans-serif">20126 Milano<br></font></span><span style="color:rgb(0,0,0);font-size:12.8px">Italy</span></p><p><font color="#000000"><span style="font-size:12.8px">email: </span><span style="font-size:12.8px"><a href="mailto:bruno.giacomazzo@unimib.it" target="_blank">bruno.giacomazzo@unimib.it</a><br></span><span style="font-size:12.8px">phone: (+39) 02 6448 2321</span></font><br><font color="#000000"><span style="font-size:12.8px">web: </span></font><span style="color:rgb(0,0,0);font-size:12.8px"><a href="http://www.brunogiacomazzo.org/" style="font-size:12.8px" target="_blank">http://www.brunogiacomazzo.org</a></span></p><p><font color="#000000">---------------------------------------------------------------------<br><span style="font-size:12.8px">There are only 10 types of people in the world:<br></span><span style="font-size:12.8px">Those who understand binary, and those who don't<br></span><span style="font-size:12.8px">----------------------------------------------------------------------<br><br></span></font></p></div></div></div></div></div>