Refactor

kimwalisch · Feb 10, 2024 · 96505d5 · 96505d5
1 parent 83e0e92
commit 96505d5
Showing 1 changed file with 6 additions and 10 deletions.
diff --git a/src/RiemannR.cpp b/src/RiemannR.cpp
@@ -212,17 +212,15 @@ long double RiemannR_prime(long double x)
 
   for (int k = 1; k < 128 && std::abs(old_sum - sum) >= epsilon; k++)
   {
-    long double k_inv = 1.0L / (long double) k;
-    term *= logx * k_inv;
+    term *= logx / k;
     old_sum = sum;
     sum += term * zetaInv[k];
   }
 
   // For k >= 128, approximate zeta(k + 1) by 1
   for (int k = 128; std::abs(old_sum - sum) >= epsilon; k++)
   {
-    long double k_inv = 1.0L / (long double) k;
-    term *= logx * k_inv;
+    term *= logx / k;
     old_sum = sum;
     sum += term;
   }
@@ -259,19 +257,17 @@ long double RiemannR(long double x)
 
   for (int k = 1; k < 128 && std::abs(old_sum - sum) >= epsilon; k++)
   {
-    long double k_inv = 1.0L / (long double) k;
-    term *= logx * k_inv;
+    term *= logx / k;
     old_sum = sum;
-    sum += term * k_inv * zetaInv[k];
+    sum += term / k * zetaInv[k];
   }
 
   // For k >= 128, approximate zeta(k + 1) by 1
   for (int k = 128; std::abs(old_sum - sum) >= epsilon; k++)
   {
-    long double k_inv = 1.0L / (long double) k;
-    term *= logx * k_inv;
+    term *= logx / k;
     old_sum = sum;
-    sum += term * k_inv;
+    sum += term / k;
   }
 
   return sum;