|
823 | 823 | @equations 2X ~ $c - X |
824 | 824 | end) |
825 | 825 |
|
826 | | - u0 = [] |
| 826 | + u0 = [rn.X => 0.0] |
827 | 827 | ps = [] |
828 | 828 | oprob = ODEProblem(rn, u0, (0.0, 100.0), ps; structural_simplify=true) |
829 | 829 | sol = solve(oprob, Tsit5(); abstol=1e-9, reltol=1e-9) |
|
851 | 851 |
|
852 | 852 | rn = complete(compose(base_rn, [internal_rn])) |
853 | 853 |
|
854 | | - u0 = [X => 3.0, internal_rn.X => 4.0] |
| 854 | + u0 = [V1 => 1.0, X => 3.0, internal_rn.V2 => 2.0, internal_rn.X => 4.0] |
855 | 855 | ps = [p => 1.0, d => 0.2, internal_rn.p => 2.0, internal_rn.d => 0.5] |
856 | 856 | oprob = ODEProblem(rn, u0, (0.0, 1000.0), ps; structural_simplify=true) |
857 | 857 | sol = solve(oprob, Rosenbrock23(); abstol=1e-9, reltol=1e-9) |
|
895 | 895 | @test equations(rn)[3] isa Equation |
896 | 896 |
|
897 | 897 | # Checks that simulations has the correct output |
898 | | - u0 = Dict([S => 1 + rand(rng), Y => 1 + rand(rng)]) |
| 898 | + u0 = Dict([S => 1 + rand(rng), X => 1 + rand(rng), Y => 1 + rand(rng)]) |
899 | 899 | ps = Dict([p => 1 + rand(rng), d => 1 + rand(rng), k => 1 + rand(rng)]) |
900 | 900 | oprob = ODEProblem(rn, u0, (0.0, 10000.0), ps; structural_simplify=true) |
901 | 901 | sol = solve(oprob, Tsit5(); abstol=1e-9, reltol=1e-9) |
@@ -1025,4 +1025,3 @@ let |
1025 | 1025 | @parameters v n |
1026 | 1026 | @test isequal(Catalyst.expand_registered_functions(equations(rn4)[1]), D(A) ~ v*(A^n)) |
1027 | 1027 | end |
1028 | | - |
0 commit comments