up

Author: TorkelE

Project.toml

@@ -64,6 +64,7 @@ HomotopyContinuation = "f213a82b-91d6-5c5d-acf7-10f1c761b327"
 LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
 NonlinearSolve = "8913a72c-1f9b-4ce2-8d82-65094dcecaec"
 OrdinaryDiffEq = "1dea7af3-3e70-54e6-95c3-0bf5283fa5ed"
+Plots = "91a5bcdd-55d7-5caf-9e0b-520d859cae80"
 Random = "9a3f8284-a2c9-5f02-9a11-845980a1fd5c"
 SafeTestsets = "1bc83da4-3b8d-516f-aca4-4fe02f6d838f"
 SciMLBase = "0bca4576-84f4-4d90-8ffe-ffa030f20462"

src/reaction_network.jl

@@ -699,6 +699,7 @@ function read_observed_options(options)
         for obs_eq in observed_eqs.args
             (obs_eq.args[1] != :~) && error("Malformed observable formula :\"$(obs_eq)\". Formula should contain two expressions separated by a '~'.")
             (obs_eq.args[2] isa Symbol) || error("Malformed observable formula :\"$(obs_eq)\". Left hand side should be a single symbol.")
+            in(obs_eq.args[2], forbidden_symbols_error) && error("A forbidden symbol ($(obs_eq.args[2])) was used as an observable name.")
             push!(observed_vars.args, obs_eq.args[2])
         end
     else  

test/dsl/dsl_options.jl

@@ -447,4 +447,12 @@ let
         end
         (p,d), 0 <--> X1 + X2
     end
+
+    # A forbidden symbol used as observable name.
+    @test_throws Exception @eval @reaction_network begin
+        @observables begin
+            t ~ X1 + X2
+        end
+        (p,d), 0 <--> X1 + X2
+    end
 end