Puts the empty interval test first

grame-cncm · Apr 4, 2024 · 1c593b3 · 1c593b3
1 parent 26b639d
commit 1c593b3
Show file tree

Hide file tree

Showing 2 changed files with 6 additions and 5 deletions.
diff --git a/interval/intervalCos.cpp b/interval/intervalCos.cpp
@@ -34,12 +34,13 @@ static double cosPi(double x)
 
 interval interval_algebra::Cos(const interval& x)
 {
-    int precision = exactPrecisionUnary(cos, 0, pow(2, x.lsb()));
-    if (precision == INT_MIN or taylor_lsb) precision = 2*x.lsb() - 1; // if x.lsb() is so small that the automatic computation doesn't work
-
     if (x.isEmpty()) {
         return empty();
     }
+
+    int precision = exactPrecisionUnary(cos, 0, pow(2, x.lsb()));
+    if (precision == INT_MIN or taylor_lsb) precision = 2*x.lsb() - 1; // if x.lsb() is so small that the automatic computation doesn't work
+
     if (x.size() >= 2*M_PI) {
         return {-1, 1, precision};
     }

diff --git a/interval/intervalCosh.cpp b/interval/intervalCosh.cpp
@@ -28,12 +28,12 @@ namespace itv {
 
 interval interval_algebra::Cosh(const interval& x)
 {
-    double v = 0;  // absolute lowest slope is at zero
-
     if (x.isEmpty()) {
         return empty();
     }
 
+    double v = 0;  // absolute lowest slope is at zero
+
     if (x.hasZero()) {
         int precision = exactPrecisionUnary(cosh, v, pow(2, x.lsb()));
         if (precision == INT_MIN or taylor_lsb)