diff options
Diffstat (limited to 'stokes.hpp')
-rw-r--r-- | stokes.hpp | 8 |
1 files changed, 4 insertions, 4 deletions
@@ -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++; } } |