numerical conditioning issue

scottshambaugh · Feb 19, 2024 · 673e767 · 673e767
1 parent 6f4ec03
commit 673e767
Show file tree

Hide file tree

Showing 2 changed files with 6 additions and 1 deletion.
diff --git a/CHANGELOG.md b/CHANGELOG.md
@@ -27,6 +27,11 @@
 ### Removed    
 
 
+## [0.11.1] - 2024-02-19
+### Changed    
+* Fix underflow issue in dvars
+
+
 ## [0.11.0] - 2023-11-05
 ### Added    
 * Python 3.12 support

diff --git a/src/monaco/dvars_sensitivity.py b/src/monaco/dvars_sensitivity.py
@@ -269,7 +269,7 @@ def calc_L(phi : np.ndarray,
     M = np.ones((m, 1))
     R = calc_R(phi, X)
     Rinv = np.linalg.inv(R)
-    Rdet = max(np.linalg.det(R), 1e-12)  # Protect for poor conditioning
+    Rdet = max(scipy.linalg.det(R), 1e-12)  # Protect for poor conditioning
 
     mu = np.linalg.inv(M.T @ Rinv @ M) @ (M.T @ Rinv @ Y)