Merge pull request #3293 from AayushSabharwal/as/hc-poly-transform

fix: minor bug fix in `PolynomialTransformation`
SciML · Jan 9, 2025 · 6143897 · 6143897
2 parents 961e25a + da923ca
commit 6143897
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/src/systems/nonlinear/homotopy_continuation.jl b/src/systems/nonlinear/homotopy_continuation.jl
@@ -348,7 +348,8 @@ function PolynomialTransformation(sys::NonlinearSystem)
     # Is there a better way to check for uniqueness? `simplify` is relatively slow
     # (maybe use the threaded version?) and `expand` can blow up expression size.
     # Could metatheory help?
-    all_non_poly_terms = mapreduce(d -> d.non_polynomial_terms, vcat, polydata)
+    all_non_poly_terms = mapreduce(
+        d -> d.non_polynomial_terms, vcat, polydata; init = BasicSymbolic[])
     unique!(all_non_poly_terms)
 
     # each variable can only be replaced by one non-polynomial expression involving