add non-autonomous errors to bifurcation computations and homotopy continuation computation.

Author: TorkelE

ext/CatalystBifurcationKitExtension/bifurcation_kit_extension.jl

@@ -3,6 +3,7 @@
 # Creates a BifurcationProblem, using a ReactionSystem as an input.
 function BK.BifurcationProblem(rs::ReactionSystem, u0_bif, ps, bif_par, args...;
                                plot_var=nothing, record_from_solution=BK.record_sol_default, jac=true, u0=[], kwargs...)
+    is_autonomous(rs) || error("Attempting to compute steady state for non-autonomous system (e.g. does some rate depend on $(rs.iv)?), this is not possible.")
 
     # Converts symbols to symbolics.
     (bif_par isa Symbol) && (bif_par = ModelingToolkit.get_var_to_name(rs)[bif_par])

ext/CatalystHomotopyContinuationExtension/homotopy_continuation_extension.jl

@@ -36,6 +36,7 @@ Notes:
 ```
   """
 function Catalyst.hc_steady_states(rs::ReactionSystem, ps; filter_negative=true, neg_thres=-1e-20, u0=[], kwargs...)
+    is_autonomous(rs) || error("Attempting to compute steady state for non-autonomous system (e.g. does some rate depend on $(rs.iv)?), this is not possible.")
     ss_poly = steady_state_polynomial(rs, ps, u0)
     sols = HC.real_solutions(HC.solve(ss_poly; kwargs...))
     reorder_sols!(sols, ss_poly, rs)

test/extensions/homotopy_continuation.jl

@@ -138,4 +138,13 @@ let
 
     # Computes bifurcation diagram.
     @test_throws Exception hc_steady_states(incomplete_network, p_start)
+end
+
+# Tests that non-autonomous system throws an error
+let 
+    rs = @reaction_network begin
+        (k,t), 0 <--> X
+    end
+    ps = [k => 1.0]
+    @test_throws Exception hc_steady_states(rs, ps)
 end