fix: fix InitializationProblems that have vectors in u0map

SebastianM-C · SebastianM-C · commit e3def3881ba8 · 2024-11-23T00:17:33.000+02:00
diff --git a/src/systems/diffeqs/abstractodesystem.jl b/src/systems/diffeqs/abstractodesystem.jl
@@ -1357,7 +1357,7 @@ function InitializationProblem{iip, specialize}(sys::AbstractODESystem,
     end
     if u0T != Union{}
         u0map = Dict(k => if symbolic_type(v) == NotSymbolic() && !is_array_of_symbolics(v)
-                         v isa AbstractArray ? u0T.(v) : u0T(v)
+                         v isa AbstractArray ? eltype(u0T).(v) : u0T(v)
                      else
                          v
                      end

Original file line number	Diff line number	Diff line change
`@@ -1357,7 +1357,7 @@ function InitializationProblem{iip, specialize}(sys::AbstractODESystem,`
`1357`	`1357`	`end`
`1358`	`1358`	`if u0T != Union{}`
`1359`	`1359`	`u0map = Dict(k => if symbolic_type(v) == NotSymbolic() && !is_array_of_symbolics(v)`
`1360`		`- v isa AbstractArray ? u0T.(v) : u0T(v)`
	`1360`	`+ v isa AbstractArray ? eltype(u0T).(v) : u0T(v)`
`1361`	`1361`	`else`
`1362`	`1362`	`v`
`1363`	`1363`	`end`