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

[Bug fix] Off-by-one in FPI counting

parent 7d16074e
No related branches found
No related tags found
No related merge requests found
...@@ -41,7 +41,7 @@ int FixedPointIterator<Factory, StateUpdater, VelocityUpdater>::run( ...@@ -41,7 +41,7 @@ int FixedPointIterator<Factory, StateUpdater, VelocityUpdater>::run(
Vector previousVelocityIterate = velocityIterate; Vector previousVelocityIterate = velocityIterate;
size_t fixedPointIteration; size_t fixedPointIteration;
for (fixedPointIteration = 1; fixedPointIteration <= fixedPointMaxIterations_; for (fixedPointIteration = 0; fixedPointIteration < fixedPointMaxIterations_;
++fixedPointIteration) { ++fixedPointIteration) {
Vector v_m; Vector v_m;
velocityUpdater->extractOldVelocity(v_m); velocityUpdater->extractOldVelocity(v_m);
...@@ -63,8 +63,10 @@ int FixedPointIterator<Factory, StateUpdater, VelocityUpdater>::run( ...@@ -63,8 +63,10 @@ int FixedPointIterator<Factory, StateUpdater, VelocityUpdater>::run(
velocityProblemSolver.solve(); velocityProblemSolver.solve();
if (velocityMatrixNorm.diff(previousVelocityIterate, velocityIterate) < if (velocityMatrixNorm.diff(previousVelocityIterate, velocityIterate) <
fixedPointTolerance_) fixedPointTolerance_) {
fixedPointIteration++;
break; break;
}
previousVelocityIterate = velocityIterate; previousVelocityIterate = velocityIterate;
} }
......
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