summaryrefslogtreecommitdiff
path: root/stokes.hpp
diff options
context:
space:
mode:
Diffstat (limited to 'stokes.hpp')
-rw-r--r--stokes.hpp8
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++;
}
}