lorentz_lorene1: Lorentz release compiled with LORENElorentz_lorene2: … LORENE2 (the one giving the Carpet error)
lorentz_lorene2_amend: … LORENE2, but with some changes in
unites.h | Reporter: | Ken Hui |
| Status: | open |
| Milestone: | ET_2021_05 |
| Version: | ET_2021_05 |
| Type: | bug |
| Priority: | minor |
| Component: | Cactus |
Comment (by Ken Hui):
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.
lorentz_lorene1: Lorentz release compiled with LORENE
lorentz_lorene2: … LORENE2 (the one giving the Carpet error)
lorentz_lorene2_amend: … LORENE2, but with some changes in unites.h
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 unites.h of LORENE2 to match that of LORENE.
--- lorene2_unites.h.ori 2022-03-11 02:08:27.660497092 +0800 +++ lorene2_unites.h.amend 2022-03-11 02:03:28.053007609 +0800 @@ -78,13 +78,13 @@ * and black holes.\ingroup (unites) */ namespace Unites { - const double g_si = 6.6738E-11 ; ///< Newton gravitational constant [SI] + const double g_si = 6.6726E-11 ; ///< Newton gravitational constant [SI] 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] - const double msol_si = 1.9885E+30 ; ///< Solar mass [kg] - const double mev_si = 1.602176565E-13 ; ///< One MeV [J] + const double msol_si = 1.989E+30 ; ///< Solar mass [kg] + const double mev_si = 1.6021892E-13 ; ///< One MeV [J] const--
Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2599/nsnstohmns-cannot-be-reproduced-using