handle non-majs too

Author: isaacsas

src/systems/jumps/jumpsystem.jl

@@ -256,10 +256,13 @@ end
 
 
 ### Functions to determine which states a jump depends on
-get_variables!(dep, jump::Union{ConstantRateJump,VariableRateJump}, variables) = get_variables!(dep, jump.rate, variables)
+function get_variables!(dep, jump::Union{ConstantRateJump,VariableRateJump}, variables) 
+    (jump.rate isa Operation) && get_variables!(dep, jump.rate, variables)
+end
 
 function get_variables!(dep, jump::MassActionJump, variables)
-    !(jump.scaled_rates isa Number) && get_variables!(dep, jump.scaled_rates, variables)
+    sr = jump.scaled_rates
+    (sr isa Operation) && get_variables!(dep, sr, variables)
     for varasop in jump.reactant_stoch
         (varasop[1].op in variables) && push!(dep, varasop[1])
     end