Is there a way to dynamically add/remove of reactions from ReactionSystem in catalyst

Specific Domains Modelling & Simulations
,
1

Hi I am trying to figure out a way to dynamically add/remove reactions once I have defined a ReactionSystem programmatically.
Suppose I define the system using this function.

function generate_creation(n,MAT)
    # Creates a vector `X` with n species.
    @species X(t)[1:n]
	@species A(t),B(t)
	@parameters k
    rxs = [
		[Reaction(k, [A,B], [X[i]]) for i in 1:n];
                [ADD REACTIONS USING MATRIX MAT]
    ]
    # Assembly and return a complete `ReactionSystem` 
    @named lp_n = ReactionSystem(rxs, t)
    return complete(lp_n)
end

I let the system reach the equilibrium. Then change the matrix mat. This changes the set of reactions. Can i update the set of reactions without needing to call this function once again.

2

Can you maybe elaborate on what you wish to do exactly? Do you want the system to reach equilibrium, and then add some reactions and continue from where you were before? Do you know beforehand already what reactions will be added and when? Depending on the answers, I see different approaches.

3

Yes It is close to the former. The reactions that will be removed will depend on the state at the time the system has reached the equilibrium. I am basically trying to implement an evolutionary algorithm. So the matrix (which creates the set of reactions) changes depending on which molecule exist at the end of the jump process/ODE solution.