From a14ed8e1bbf105b8b66e64faa6fb78141f528b3e Mon Sep 17 00:00:00 2001
From: Matthieu Schaller <schaller@strw.leidenuniv.nl>
Date: Mon, 4 Feb 2019 22:36:15 +0100
Subject: [PATCH] Correct initial guess for the Newton scheme in the EAGLE
 cooling code.

---
 src/cooling/EAGLE/cooling.c | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/cooling/EAGLE/cooling.c b/src/cooling/EAGLE/cooling.c
index 1b93417dde..7193a12f37 100644
--- a/src/cooling/EAGLE/cooling.c
+++ b/src/cooling/EAGLE/cooling.c
@@ -57,7 +57,7 @@ static const float newton_tolerance = 1.0e-4;
 static const float bisection_tolerance = 1.0e-6;
 static const float rounding_tolerance = 1.0e-4;
 static const double bracket_factor = 1.0488088481701515;    /* sqrt(1.1) */
-static const double newton_log_u_guess_cgs = 1.414213562e6; /* log10(2e12) */
+static const double newton_log_u_guess_cgs = 12.30103; /* log10(2e12) */
 
 /**
  * @brief Find the index of the current redshift along the redshift dimension
-- 
GitLab