add validation

Author: vyudu

src/network_analysis.jl

@@ -215,6 +215,13 @@ end
     **Warning**: Unlike other Catalyst functions, the `fluxmat` function will return a `Matrix{Num}` in the symbolic case. This is to allow easier computation of the matrix decomposition of the ODEs, and to ensure that multiplying the sparse form of the matrix will work.
 """
 function fluxmat(rn::ReactionSystem, pmap::Dict = Dict(); sparse=false)
+    deps = Set()
+    for (i, rx) in enumerate(reactions(rn))
+        empty!(deps)
+        get_variables!(deps, rx.rate, species(rn))
+        (!isempty(deps)) && (error("Reaction $rx's rate constant depends on species $(join(deps, ", ")). `adjacencymat` cannot support rate constants of this form."))
+    end
+
     rates = if isempty(pmap)
         reactionrates(rn)
     else
@@ -290,6 +297,10 @@ function massactionvector(rn::ReactionSystem, scmap::Dict = Dict(); combinatoric
     rxs = reactions(rn)
     sm = speciesmap(rn)
 
+    for rx in rxs
+        !ismassaction(rx, rn) && error("The supplied ReactionSystem has non-mass action reaction $rx. The `massactionvector` can only be constructed for mass action networks.")
+    end
+
     specs = if isempty(scmap) 
         species(rn)
     else
@@ -1001,7 +1012,6 @@ function iscomplexbalanced(rs::ReactionSystem, parametermap::Dict; sparse = fals
         error("The supplied ReactionSystem has reactions that are not ismassaction. Testing for being complex balanced is currently only supported for pure mass action networks.")
     end
 
-    L = laplacianmat(rs, parametermap; sparse)
     D = incidencemat(rs; sparse)
     S = netstoichmat(rs; sparse)
 
@@ -1053,6 +1063,13 @@ end
         - In `adjacencymat`, the rows and columns both represent complexes, and an entry (c1, c2) is non-zero if there is a reaction c1 --> c2.
 """
 function adjacencymat(rn::ReactionSystem, pmap::Dict = Dict(); sparse = false)
+    deps = Set()
+    for (i, rx) in enumerate(reactions(rn))
+        empty!(deps)
+        get_variables!(deps, rx.rate, species(rn))
+        (!isempty(deps)) && (error("Reaction $rx's rate constant depends on species $(join(deps, ", ")). `adjacencymat` cannot support rate constants of this form."))
+    end
+
     rates = if isempty(pmap)
         reactionrates(rn)
     else