From e823f8e53c36d30e9a0e806ab358d31494ba2b32 Mon Sep 17 00:00:00 2001
From: Philippe WAUTELET <philippe.wautelet@aero.obs-mip.fr>
Date: Mon, 9 Dec 2019 17:38:09 +0100
Subject: [PATCH] Philippe 09/12/2019: BITREP without floor
---
src/LIB/BITREP/br_transcendentals.cpp | 6 ++++++
1 file changed, 6 insertions(+)
diff --git a/src/LIB/BITREP/br_transcendentals.cpp b/src/LIB/BITREP/br_transcendentals.cpp
index 8b1fdfc6c..2bbbb4dae 100644
--- a/src/LIB/BITREP/br_transcendentals.cpp
+++ b/src/LIB/BITREP/br_transcendentals.cpp
@@ -823,7 +823,13 @@ double __internal_exp_kernel(double x, int scale)
double t, z;
int i;
+#if 0
t = std::floor (x*1.4426950408889634e+0 + 4.99999999999999945e-1);
+#else
+ t = x*1.4426950408889634e+0 + 4.99999999999999945e-1;
+ if (t<0) {t=t-1.;}
+ t=(int)t;
+#endif
i = (int)t;
z = __BITREPFMA(t, -6.9314718055994529e-1, x);
z = __BITREPFMA(t, -2.3190468138462996e-17, z);
--
GitLab