Skip to content
Snippets Groups Projects
Commit aa0e7c9f authored by Elias Pipping's avatar Elias Pipping
Browse files

[Output] Special handling for lambda = 0

parent 7eb52129
No related branches found
No related tags found
No related merge requests found
......@@ -72,7 +72,8 @@ FixedPointIterator<Factory, StateUpdater, VelocityUpdater, ErrorNorm>::run(
ScalarVector newAlpha;
stateUpdater->extractAlpha(newAlpha);
if (errorNorm_.diff(alpha, newAlpha) < fixedPointTolerance_) {
if (lambda_ < 1e-12 or
errorNorm_.diff(alpha, newAlpha) < fixedPointTolerance_) {
fixedPointIteration++;
break;
}
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment