inverse to linsolve for AS quadratic

ZIB-IOL · Oct 11, 2024 · e2b2012 · e2b2012
1 parent d203d56
commit e2b2012
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/active_set_quadratic.jl b/src/active_set_quadratic.jl
@@ -36,7 +36,7 @@ function detect_quadratic_function(grad!, x0; test=true)
         X[:, i] .-= x0
         G[:, i] .= storage .- g0
     end
-    A = G * inv(X)
+    A = G \ X
     b = g0 - A * x0
     if test
         x_test = randn(T, n)