<html>#2599: nsnstohmns cannot be reproduced using LORENE2
<table style='border-spacing: 1ex 0pt; '>
<tr><td style='text-align:right'> Reporter:</td><td>Ken Hui</td></tr>
<tr><td style='text-align:right'> Status:</td><td>open</td></tr>
<tr><td style='text-align:right'>Milestone:</td><td>ET_2021_05</td></tr>
<tr><td style='text-align:right'> Version:</td><td>ET_2021_05</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>Cactus</td></tr>
</table>
<p>Comment (by Ken Hui):</p>
<p>I followed the particular ticket you mentioned to plot the Psi4 and maximum density, but I don’t think they give any information related to the failure of the simulation.</p>
<p><img alt="" src="https://bitbucket.org/repo/8zBLME5/images/2417778790-Psi4_compare.png">
<img alt="" src="https://bitbucket.org/repo/8zBLME5/images/2427778410-dens_compare.png">
lorentz_lorene1: Lorentz release compiled with LORENE<br>
lorentz_lorene2: … LORENE2 (the one giving the Carpet error)<br>
lorentz_lorene2_amend: … LORENE2, but with some changes in <code>unites.h</code> </p>
<p>The values of the density differ even at the first step between lorentz_lorene1 and lorentz_lorene2, which is caused by differences in some physical constants, so I amended some numbers in <code>unites.h</code> of LORENE2 to match that of LORENE. </p>
<div class="codehilite"><pre><span></span><span class="gd">--- lorene2_unites.h.ori 2022-03-11 02:08:27.660497092 +0800</span>
<span class="gi">+++ lorene2_unites.h.amend 2022-03-11 02:03:28.053007609 +0800</span>
<span class="gu">@@ -78,13 +78,13 @@</span>
* and black holes.\ingroup (unites)
*/
namespace Unites {
<span class="gd">- const double g_si = 6.6738E-11 ; ///< Newton gravitational constant [SI]</span>
<span class="gi">+ const double g_si = 6.6726E-11 ; ///< Newton gravitational constant [SI]</span>
const double c_si = 2.99792458E+8 ; ///< Velocity of light [m/s]
const double kB_si = 1.3806488E-23 ; ///< Boltzmann constant [J/K]
const double rhonuc_si = 1.66E+17 ; ///< Nuclear density [kg/m3] (arbitrary)
const double km_si = 1.E+3 ; ///< One kilometer [m]
<span class="gd">- const double msol_si = 1.9885E+30 ; ///< Solar mass [kg]</span>
<span class="gd">- const double mev_si = 1.602176565E-13 ; ///< One MeV [J]</span>
<span class="gi">+ const double msol_si = 1.989E+30 ; ///< Solar mass [kg]</span>
<span class="gi">+ const double mev_si = 1.6021892E-13 ; ///< One MeV [J]</span>
const
<p>--<br/>
Ticket URL: <a href='https://bitbucket.org/einsteintoolkit/tickets/issues/2599/nsnstohmns-cannot-be-reproduced-using'>https://bitbucket.org/einsteintoolkit/tickets/issues/2599/nsnstohmns-cannot-be-reproduced-using</a></p>
</html>