fix: Avoid non-converging flow due to large perturbation

danieledler · danieledler · commit 1a0747a0e867 · 2024-12-11T10:38:49.000+01:00
diff --git a/src/utils/FlowCalculator.cpp b/src/utils/FlowCalculator.cpp
@@ -301,8 +301,8 @@ IterationResult powerIterate(double alpha, Iteration&& iter)
     err = iter(iterations, alpha, beta);
 
     // Perturb the system if equilibrium
-    if (std::abs(err - oldErr) < 1e-15) {
-      alpha += 1.0e-10;
+    if (std::abs(err - oldErr) < 1e-17) {
+      alpha += 1.0e-12;
       beta = 1.0 - alpha;
     }
 

Original file line number	Diff line number	Diff line change
`@@ -301,8 +301,8 @@ IterationResult powerIterate(double alpha, Iteration&& iter)`
`301`	`301`	`err = iter(iterations, alpha, beta);`
`302`	`302`
`303`	`303`	`// Perturb the system if equilibrium`
`304`		`- if (std::abs(err - oldErr) < 1e-15) {`
`305`		`- alpha += 1.0e-10;`
	`304`	`+ if (std::abs(err - oldErr) < 1e-17) {`
	`305`	`+ alpha += 1.0e-12;`
`306`	`306`	`beta = 1.0 - alpha;`
`307`	`307`	`}`
`308`	`308`