From 09ba6851a712af711f1b3410a227bfd38c288027 Mon Sep 17 00:00:00 2001
From: Mladen Ivkovic <mladen.ivkovic@hotmail.com>
Date: Fri, 9 Feb 2024 14:33:39 +0000
Subject: [PATCH] use doubles to compute F^2, as on some occations it may
exceed the single precision limits
---
src/rt/GEAR/rt_unphysical.h | 8 ++++----
1 file changed, 4 insertions(+), 4 deletions(-)
diff --git a/src/rt/GEAR/rt_unphysical.h b/src/rt/GEAR/rt_unphysical.h
index 39bb0086d4..af689dd990 100644
--- a/src/rt/GEAR/rt_unphysical.h
+++ b/src/rt/GEAR/rt_unphysical.h
@@ -72,11 +72,11 @@ __attribute__((always_inline)) INLINE static void rt_check_unphysical_state(
}
/* Check for too high fluxes */
- const float flux2 = flux[0] * flux[0] + flux[1] * flux[1] + flux[2] * flux[2];
- const float flux_norm = sqrtf(flux2);
- const float flux_max = rt_params.reduced_speed_of_light * *energy_density;
+ const double flux2 = flux[0] * flux[0] + flux[1] * flux[1] + flux[2] * flux[2];
+ const double flux_norm = sqrt(flux2);
+ const double flux_max = rt_params.reduced_speed_of_light * *energy_density;
if (flux_norm > flux_max) {
- const float correct = flux_max / flux_norm;
+ const double correct = flux_max / flux_norm;
flux[0] *= correct;
flux[1] *= correct;
flux[2] *= correct;
--
GitLab