Vincent Du
Vincent Du
Thank you both!
Also, need to add a test checking conditional complex balance, i.e. a reaction network that is complex balance for some set of rate parameters but not for a different set.
Methodology for the persistence check: 1. Detect the set of minimal siphons in the reaction network, following [Shiu et al., 2009](https://people.tamu.edu/~annejls/SiphonsInCRNs.pdf) 2. Check necessary and sufficient conditions for persistence outlined...
[Absolute concentration robustness](https://arxiv.org/pdf/2401.00078) [Higher deficiency algorithm and others](https://etd.ohiolink.edu/acprod/odb_etd/ws/send_file/send?accession=osu1307122057&disposition=inline) [Monotone Systems will not have multistationarity](https://arxiv.org/abs/1411.2896)
``` rn = @reaction_network begin @parameters k ksq = k^2 (k, ksq), A B end rn = complete(rn) u0_guess = [:A => 1., :B => 1.] p_start = [:k =>...
Sorry, forgot about this, will look this afternoon
Will look at this tonight
It seems like it is hanging due to some kind of type inference. Cutting the compiler shows that the compiler is hanging in `specialize_method`, which does some type inference. When...
What's the expected behavior here? should `isequal(rn1.X1, rn1.X2)` return true?
Hmm what would be the way to check that? I don't think there's a way that they could generate the same code when expanded if the observable IVs are inferred...