up

Author: TorkelE

src/reactionsystem.jl

@@ -207,6 +207,13 @@ end
 
 ### ReactionSystem Structure ###
 
+# Constant storing all reaction system fields (in order). Used to check whether the `ReactionSystem`
+# structure have been updated (in the `reactionsystem_uptodate` function).
+const reactionsystem_fields = [:eqs, :rxs, :iv, :sivs, :unknowns, :species, :ps, :var_to_name, 
+                               :observed, :name, :systems, :defaults, :connection_type, 
+                               :networkproperties, :combinatoric_ratelaws, :continuous_events, 
+                               :discrete_events, :metadata, :complete]
+
 """
 $(TYPEDEF)
 
@@ -1025,6 +1032,14 @@ end
 
 ### General `ReactionSystem`-specific Functions ###
 
+# Checks if the `ReactionSystem` structure have been updated without also updating the 
+# `reactionsystem_fields` constant. If this is the case, returns `false`. This is used in 
+# certain functionalities which would break if the `ReactionSystem` structure is updated without
+# also updating tehse functionalities.
+function reactionsystem_uptodate()
+    fieldnames(ReactionSystem) == reactionsystem_fields
+end
+
 # used in the `__unpacksys` function.
 function __unpacksys(rn)
     ex = :(begin end)

test/reactionsystem_core/reactionsystem.jl

@@ -740,4 +740,12 @@ let
     @test nameof(ModelingToolkit.get_iv(empty_network)) == :t
     @test length(ModelingToolkit.get_unknowns(empty_network)) == 0
     @test length(ModelingToolkit.get_ps(empty_network)) == 0
+end
+
+# Checks that the `reactionsystem_uptodate` function work. If it does not, the ReactionSystem
+# strcuture's fields have been updated, without updating the `reactionsystem_fields` costant. If so,
+# there are several places in the code where the `reactionsystem_uptodate` function is called, here
+# the code might need adaptation to take the updated reaction system into account.
+let
+    @test Catalyst.reactionsystem_uptodate()
 end