Reduce template code bloat

Author: kimwalisch

src/LogarithmicIntegral.cpp

@@ -315,10 +315,8 @@ int64_t Li(int64_t x)
 #endif
   if (x > 1e8)
     return (int64_t) ::Li((long double) x);
-  if (x > 100)
-    return (int64_t) ::Li((double) x);
   else
-    return (int64_t) ::Li((float) x);
+    return (int64_t) ::Li((double) x);
 }
 
 int64_t Li_inverse(int64_t x)
@@ -329,10 +327,8 @@ int64_t Li_inverse(int64_t x)
 #endif
   if (x > 1e8)
     return Li_inverse_overflow_check<long double>(x);
-  if (x > 100)
-    return Li_inverse_overflow_check<double>(x);
   else
-    return Li_inverse_overflow_check<float>(x);
+    return Li_inverse_overflow_check<double>(x);
 }
 
 #ifdef HAVE_INT128_T
@@ -345,10 +341,8 @@ int128_t Li(int128_t x)
 #endif
   if (x > 1e8)
     return (int128_t) ::Li((long double) x);
-  if (x > 100)
-    return (int128_t) ::Li((double) x);
   else
-    return (int128_t) ::Li((float) x);
+    return (int128_t) ::Li((double) x);
 }
 
 int128_t Li_inverse(int128_t x)
@@ -359,10 +353,8 @@ int128_t Li_inverse(int128_t x)
 #endif
   if (x > 1e8)
     return Li_inverse_overflow_check<long double>(x);
-  if (x > 100)
-    return Li_inverse_overflow_check<double>(x);
   else
-    return Li_inverse_overflow_check<float>(x);
+    return Li_inverse_overflow_check<double>(x);
 }
 
 #endif

src/RiemannR.cpp

@@ -549,10 +549,8 @@ int64_t RiemannR(int64_t x)
 #endif
   if (x > 1e8)
     return (int64_t) ::RiemannR((long double) x);
-  if (x > 100)
-    return (int64_t) ::RiemannR((double) x);
   else
-    return (int64_t) ::RiemannR((float) x);
+    return (int64_t) ::RiemannR((double) x);
 }
 
 int64_t RiemannR_inverse(int64_t x)
@@ -563,10 +561,8 @@ int64_t RiemannR_inverse(int64_t x)
 #endif
   if (x > 1e8)
     return RiemannR_inverse_overflow_check<long double>(x);
-  if (x > 100)
-    return RiemannR_inverse_overflow_check<double>(x);
   else
-    return RiemannR_inverse_overflow_check<float>(x);
+    return RiemannR_inverse_overflow_check<double>(x);
 }
 
 #ifdef HAVE_INT128_T
@@ -579,10 +575,8 @@ int128_t RiemannR(int128_t x)
 #endif
   if (x > 1e8)
     return (int128_t) ::RiemannR((long double) x);
-  if (x > 100)
-    return (int128_t) ::RiemannR((double) x);
   else
-    return (int128_t) ::RiemannR((float) x);
+    return (int128_t) ::RiemannR((double) x);
 }
 
 int128_t RiemannR_inverse(int128_t x)
@@ -593,10 +587,8 @@ int128_t RiemannR_inverse(int128_t x)
 #endif
   if (x > 1e8)
     return RiemannR_inverse_overflow_check<long double>(x);
-  if (x > 100)
-    return RiemannR_inverse_overflow_check<double>(x);
   else
-    return RiemannR_inverse_overflow_check<float>(x);
+    return RiemannR_inverse_overflow_check<double>(x);
 }
 
 #endif