Merge pull request #733 from SciML/1.6_to_1.9

Test on 1.9 instead of 1.6

.github/workflows/CI.yml

@@ -15,7 +15,7 @@ jobs:
           - Core
         version:
           - '1'
-          - '1.6'
+          - '1.9'
     steps:
       - uses: actions/checkout@v4
       - uses: julia-actions/setup-julia@v1

Project.toml

@@ -43,7 +43,7 @@ Requires = "1.0"
 SymbolicUtils = "1.0.3"
 Symbolics = "5.0.3"
 Unitful = "1.12.4"
-julia = "1.6"
+julia = "1.9"
 
 [extras]
 DomainSets = "5b8099bc-c8ec-5219-889f-1d9e522a28bf"

test/miscellaneous_tests/nonlinear_solve.jl

@@ -16,7 +16,7 @@ let
     steady_state_network_1 = @reaction_network begin
         (k1, k2), ∅ ↔ X1
         (k3, k4), ∅ ↔ 3X2
-        (k5, k6), ∅ ↔ X3 + X4
+        (k5, k6*X1), ∅ ↔ X3
     end
 
     # Creates NonlinearProblem.
@@ -31,10 +31,10 @@ let
     # Tests solutions are correct.
     @test isapprox(sol1[1], p[1] / p[2]; atol=1e-10)
     @test isapprox(sol1[2]^3 / factorial(3), p[3] / p[4]; atol=1e-10)
-    @test isapprox(sol1[3] * sol1[4], p[5] / p[6]; atol=1e-10)
+    @test isapprox(sol1[3], (p[5] * p[2]) / (p[6] * p[1]); atol=1e-10)
     @test isapprox(sol2[1], p[1] / p[2]; atol=1e-10)
     @test isapprox(sol2[2]^3 / factorial(3), p[3] / p[4]; atol=1e-10)
-    @test isapprox(sol2[3] * sol2[4], p[5] / p[6]; atol=1e-10)
+    @test isapprox(sol2[3], (p[5] * p[2]) / (p[6] * p[1]); atol=1e-10)
 end
 
 # Creates a system with multiple steady states.