From e56493d71c855e532d0bb47191c5a703521226a5 Mon Sep 17 00:00:00 2001
From: Aayush Sabharwal <aayush.sabharwal@gmail.com>
Date: Sat, 7 Jun 2025 13:42:47 +0530
Subject: [PATCH] fix: properly retain defaults in `generate_initializesystem`

---
 src/systems/nonlinear/initializesystem.jl | 4 ++++
 1 file changed, 4 insertions(+)

diff --git a/src/systems/nonlinear/initializesystem.jl b/src/systems/nonlinear/initializesystem.jl
index b1ca187e95..f6fdd6eeda 100644
--- a/src/systems/nonlinear/initializesystem.jl
+++ b/src/systems/nonlinear/initializesystem.jl
@@ -167,6 +167,8 @@ function generate_initializesystem_timevarying(sys::AbstractSystem;
     for k in keys(defs)
         defs[k] = substitute(defs[k], paramsubs)
     end
+    initials = Dict(k => v for (k, v) in pmap if isinitial(k))
+    merge!(defs, initials)
     isys = System(Vector{Equation}(eqs_ics),
         vars,
         pars;
@@ -280,6 +282,8 @@ function generate_initializesystem_timeindependent(sys::AbstractSystem;
     for k in keys(defs)
         defs[k] = substitute(defs[k], paramsubs)
     end
+    initials = Dict(k => v for (k, v) in pmap if isinitial(k))
+    merge!(defs, initials)
     isys = System(Vector{Equation}(eqs_ics),
         vars,
         pars;