test: test initial values dependent on constants

SciML · Oct 14, 2024 · 34a8c8e · 34a8c8e
1 parent 486ae04
commit 34a8c8e
Showing 1 changed file with 13 additions and 0 deletions.
diff --git a/test/constants.jl b/test/constants.jl
@@ -37,3 +37,16 @@ eqs = [D(x) ~ β]
 simp = structural_simplify(sys)
 
 @test isempty(MT.collect_constants(nothing))
+
+@testset "Issue#3044" begin
+    @constants h = 1
+    @parameters τ = 0.5 * h
+    @variables x(t) = h
+    eqs = [D(x) ~ (h - x) / τ]
+
+    @mtkbuild fol_model = ODESystem(eqs, t)
+
+    prob = ODEProblem(fol_model, [], (0.0, 10.0))
+    @test prob[x] ≈ 1
+    @test prob.ps[τ] ≈ 0.5
+end