[ET Trac] #2599: nsnstohmns cannot be reproduced using LORENE2

[Ken Hui]

#2599: nsnstohmns cannot be reproduced using LORENE2

 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.

![](https://bitbucket.org/repo/8zBLME5/images/2417778790-Psi4_compare.png)
![](https://bitbucket.org/repo/8zBLME5/images/2427778410-dens_compare.png)
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 double r_unit = 1.e4 ;  ///< Lorene's unit of length = 10 km
     const double v_unit = c_si ; ///< Lorene's unit of velocity = c 
```

‌

--
Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2599/nsnstohmns-cannot-be-reproduced-using
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.einsteintoolkit.org/pipermail/trac/attachments/20220310/4b29ce3a/attachment.html