Some tweaking of the Stokes finder, how the lines are recorded, and fewer executables.

1 files changed, 4 insertions, 4 deletions

diff --git a/stokes.hpp b/stokes.hpp
index 8c0c1fb..60a3f55 100644
--- a/stokes.hpp
+++ b/stokes.hpp
@@ -237,14 +237,14 @@ public:
 
       newCost = cNew.totalCost(nt);
 
-      δ *= 1.5;
+      δ *= 2;
     }
 
     std::cerr << newCost << " " << stepSize << " " << δ << std::endl;
 
     gs = cNew.gs;
 
-    return {δ / 1.5, stepSize};
+    return {δ / 2, stepSize};
   }
 
   void relaxNewton(unsigned nt, Real δ₀, Real ε, unsigned maxSteps) {
@@ -254,9 +254,9 @@ public:
     while (stepSize > ε && steps < maxSteps) {
       std::tie(δ, stepSize) = relaxStepNewton(nt, δ);
       if (δ > 100) {
-        nt++;
+        break;
       }
-      δ /= 2;
+      δ /= 3;
       steps++;
     }
   }