Damp initial iterations of Broyden's method used in ID parameter control

sxs-collaboration · Nov 15, 2024 · e5291b3 · e5291b3
1 parent a0a5f64
commit e5291b3
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/support/Pipelines/Bbh/ControlId.py b/support/Pipelines/Bbh/ControlId.py
@@ -298,7 +298,8 @@ def Residual(u):
         F = Residual(u)
 
         # Update the Jacobian using Broyden's method
-        J += np.outer(F, Delta_u) / np.dot(Delta_u, Delta_u)
+        damping_factor = 1.0 - np.exp(-0.5 * iteration)
+        J += np.outer(F, Delta_u) / np.dot(Delta_u, Delta_u) * damping_factor
 
     data_file.close()