<html>#2119: The binary neutron star gallery example gives different results with the release candiate.
<table style='border-spacing: 1ex 0pt; '>
<tr><td style='text-align:right'> Reporter:</td><td>Peter Diener</td></tr>
<tr><td style='text-align:right'>   Status:</td><td>open</td></tr>
<tr><td style='text-align:right'>Milestone:</td><td>ET_2019_10</td></tr>
<tr><td style='text-align:right'>  Version:</td><td>development version</td></tr>
<tr><td style='text-align:right'>     Type:</td><td>bug</td></tr>
<tr><td style='text-align:right'> Priority:</td><td>minor</td></tr>
<tr><td style='text-align:right'>Component:</td><td>Other</td></tr>
</table>

<p>Comment (by Roland Haas):</p>
<p>I reran the same code, same ID, on the same workstation using gcc’s <code>-Ofast</code> rather than <code>-O3</code> which should give an indication on how much any part of the simulation is affected by relaxed IEEE conformance (in the C code) or term reordering (in the Fortran code where rules are not as strict as in C anyway). </p>
<p>The results are rather spectacular in the post merger waveform and density evolution:</p>
<p><img alt="" src="https://bitbucket.org/repo/8zBLME5/images/2404527760-psi4_Zwicky.png" />
where the waves shut off after t~2400 when <code>-Ofast</code> is used. Similarly for density (similar to what Peter had already found) the result is quite noticable:</p>
<p><img alt="" src="https://bitbucket.org/repo/8zBLME5/images/1420980466-rho_Optimize.png" />
in the post-merger phase and already somewhat visible during merger when density peaks.</p>
<p>This simulation had a rather coarse resolution for the stars (0.28 Msun where one would normally want have at least 0.25 Msun). </p>
<p>This would mean that certainly the residual differences in post-merger evolution that Peter saw could be caused by just running the code on different clusters with different compilers and different optimization settings compared to the cluster originally used when producing the gallery data. </p>
<p>I will re-run once more, this time using <code>-O3</code> and the <em>new</em> Meudon-BNS thorn ie <a data-is-external-link="true" href="https://bitbucket.org/einsteintoolkit/tickets/commits/1b9d6a3eff34224a357a61e3b74d5468f8efa82b" rel="nofollow">1b9d6a3eff34224a357a61e3b74d5468f8efa82b</a>.</p>
<p>--<br/>
Ticket URL: <a href='https://bitbucket.org/einsteintoolkit/tickets/issues/2119/the-binary-neutron-star-gallery-example'>https://bitbucket.org/einsteintoolkit/tickets/issues/2119/the-binary-neutron-star-gallery-example</a></p>
</html>