Scalarize initial conditions of vectors

Author: hersle

src/systems/nonlinear/initializesystem.jl

@@ -45,7 +45,7 @@ function generate_initializesystem(sys::ODESystem;
                     push!(defs, y => x)
                 elseif y isa Symbolics.Arr
                     # TODO: don't scalarize arrays
-                    push!(defs, (collect(y) .=> x)...)
+                    merge!(defs, Dict(scalarize(y .=> x)))
                 elseif y isa Symbolics.BasicSymbolic
                     # y is a derivative expression expanded; add it to the initialization equations
                     push!(eqs_ics, y ~ x)