Merge pull request #2167 from SciML/myb/1.9_regression

Work around Julia 1.9.0 regression

Author: YingboMa

src/structural_transformation/bareiss.jl

@@ -143,6 +143,9 @@ end
     V = M[(k + 1):end, (k + 1):end]
     V .= exactdiv.(V .* pivot .- M[(k + 1):end, k] * M[k, (k + 1):end]', prev_pivot)
     zero!(M, (k + 1):size(M, 1), k)
+    if M isa AbstractSparseMatrixCSC
+        dropzeros!(M)
+    end
 end
 
 function bareiss_update_virtual_colswap!(zero!, M::AbstractMatrix, k, swapto, pivot,

test/state_selection.jl

@@ -250,7 +250,8 @@ let
             D(dw) => ddw]
 
     # equations ------------------------------------------------------------------
-    flow(x, dp) = K1 * abs(dp) * abs(x) + K2 * sqrt(abs(dp)) * abs(x) + K3 * abs(dp) * x^2
+    # sqrt -> log as a hack
+    flow(x, dp) = K1 * abs(dp) * abs(x) + K2 * log(abs(dp)) * abs(x) + K3 * abs(dp) * x^2
     xm = xf / x_f_fullscale
     Δp1 = p_s - p1
     Δp2 = p2