init

Author: TorkelE

test/reactionsystem_core/custom_crn_functions.jl

@@ -49,7 +49,7 @@ let
         ps = rnd_ps(custom_function_network_1, rng; factor)
         t = rand(rng)
         @test f_eval(custom_function_network_1, u0, ps, t) ≈ f_eval(custom_function_network_2, u0, ps, t)
-        @test jac_eval(custom_function_network_1, u0, ps, t) ≈ jac_eval(custom_function_network_2, u0, ps, t)
+        @test_broken jac_eval(custom_function_network_1, u0, ps, t) ≈ jac_eval(custom_function_network_2, u0, ps, t) # Weird error dur to some SciML update. Reported in https://github.com/SciML/ModelingToolkit.jl/issues/3447.
         @test g_eval(custom_function_network_1, u0, ps, t) ≈ g_eval(custom_function_network_2, u0, ps, t)
     end
 end
@@ -101,20 +101,20 @@ end
 # Tests `ReactionSystem`s.
 let
     @species x(t) y(t)
-    @parameters k v n 
+    @parameters k v n
     rs1 = @reaction_network rs begin
         mm(x, v, k), 0 --> x
         mmr(x, v, k), 0 --> x
         hill(x, v, k, n), 0 --> x
         hillr(x, v, k, n), 0 --> x
-        hillar(x, y, v, k, n), 0 --> x    
+        hillar(x, y, v, k, n), 0 --> x
     end
     rs2 = @reaction_network rs begin
         v * x / (x + k), 0 --> x
         v * k / (x + k), 0 --> x
         v * (x^n) / (x^n + k^n), 0 --> x
         v * (k^n) / (x^n + k^n), 0 --> x
-        v * (x^n) / (x^n + y^n + k^n), 0 --> x    
+        v * (x^n) / (x^n + y^n + k^n), 0 --> x
     end
 
     @test Catalyst.expand_registered_functions(rs1) == rs2
@@ -123,14 +123,14 @@ end
 # Tests `Reaction`s.
 let
     @species x(t) y(t)
-    @parameters k v n 
-    
+    @parameters k v n
+
     r1 = @reaction mm(x, v, k), 0 --> x
     r2 = @reaction mmr(x, v, k), 0 --> x
     r3 = @reaction hill(x, v, k, n), 0 --> x
     r4 = @reaction hillr(x, v, k, n), 0 --> x
     r5 = @reaction hillar(x, y, v, k, n), 0 --> x + y
-    
+
     @test isequal(Catalyst.expand_registered_functions(r1).rate, v * x / (x + k))
     @test isequal(Catalyst.expand_registered_functions(r2).rate, v * k / (x + k))
     @test isequal(Catalyst.expand_registered_functions(r3).rate, v * (x^n) / (x^n + k^n))
@@ -143,14 +143,14 @@ end
 let
     @parameters T
     @variables X(T) Y(T)
-    @parameters K V N 
-    
+    @parameters K V N
+
     eq1 = 0 ~ mm(X, V, K)
     eq2 = 0 ~ mmr(X, V, K)
     eq3 = 0 ~ hill(X, V, K, N)
     eq4 = 0 ~ hillr(X, V, K, N)
     eq5 = 0 ~ hillar(X, Y, V, K, N)
-    
+
     @test isequal(Catalyst.expand_registered_functions(eq1), 0 ~ V * X / (X + K))
     @test isequal(Catalyst.expand_registered_functions(eq2), 0 ~ V * K / (X + K))
     @test isequal(Catalyst.expand_registered_functions(eq3), 0 ~ V * (X^N) / (X^N + K^N))
@@ -166,7 +166,7 @@ let
     @parameters v K
     eqs = [
         Reaction(mm(X,v,K), [], [X]),
-        mm(V,v,K) ~ V + 1    
+        mm(V,v,K) ~ V + 1
     ]
     @named rs = ReactionSystem(eqs, t)
 
@@ -211,7 +211,7 @@ let
         (v * (X^n) / (X^n + K^n) > 1000.0) => [X ~ v * (K^n) / (X^n + K^n) + 2]
     ]
     continuous_events = ModelingToolkit.SymbolicContinuousCallback.(continuous_events)
-    discrete_events = ModelingToolkit.SymbolicDiscreteCallback.(discrete_events)    
+    discrete_events = ModelingToolkit.SymbolicDiscreteCallback.(discrete_events)
     @test isequal(only(Catalyst.get_rxs(rs_expanded)).rate, v0 + v * (X^n) / (X^n + Y^n + K^n))
     @test isequal(get_continuous_events(rs_expanded), continuous_events)
     @test isequal(get_discrete_events(rs_expanded), discrete_events)