Remove Symbolics.scalarize call for reduced equations

Author: bowenszhu

src/deim.jl

@@ -161,7 +161,7 @@ function deim(sys::ODESystem, snapshot::AbstractMatrix, pod_dim::Integer;
     reduced_rhss, linear_projection_eqs = deim(dvs, A, g, F, ŷ, V, U; kwargs...)
 
     reduced_deqs = D.(ŷ) ~ reduced_rhss
-    @set! sys.eqs = [Symbolics.scalarize(reduced_deqs); eqs]
+    @set! sys.eqs = [reduced_deqs; eqs]
 
     old_observed = ModelingToolkit.get_observed(sys)
     fullstates = [map(eq -> eq.lhs, old_observed); dvs; ModelingToolkit.get_unknowns(sys)]