From e5291b3425888f4bf86313d94e5dd2d070dc1a33 Mon Sep 17 00:00:00 2001 From: Iago Mendes Date: Fri, 15 Nov 2024 11:29:18 -0800 Subject: [PATCH] Damp initial iterations of Broyden's method used in ID parameter control --- support/Pipelines/Bbh/ControlId.py | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/support/Pipelines/Bbh/ControlId.py b/support/Pipelines/Bbh/ControlId.py index 1fee4da015f66..f6d46cdf25675 100644 --- 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()