Sam Isaacson

Ignore the algebraic equation, that was from my misunderstanding of Chris’s reply. I thought the algebraic equation would have precedence and lead to the elimination of one ODE and conversion...

@YingboMa Once we get @kaandocal's PR merged (i.e. adding tests for accuracy/correctness, cleaning up and such), it would be great to work on making the conservation law detection and reduction...

Another reference discussing methods for calculating conservation laws: https://www.sciencedirect.com/science/article/pii/S0301462203002503?casa_token=lwJ3ukNRX6cAAAAA:oA1T3kKDFf9ebVYhPkO-j3v7XneNP2YhnfEYBmrOhAinjkUpmOEDGVqmp03--EnyjMJndAMg#BIB8

A library with mostly BSD code for doing structural analysis (including detecting conservation laws): Paper: https://doi.org/10.1016/j.biosystems.2018.05.008 Library: https://github.com/sys-bio/Libstructural

I think the network API stuff needs more testing before we start advertising it. I'd be hesistant right now to say more about it in a paper than that we...

I guess where this would get tricky is something like `2A + B --> C` with `A` constant and a rate `k`. Should this be rewritten to `B --> C`...

Also need to consider boundary condition species (which do not get ODEs from reactions, but may have constraint equations or events).

Closed by https://github.com/SciML/Catalyst.jl/pull/514

Reopening as I realized this won't work beyond ODESystems and NonlinearSystems until if/when events and/or constraints can be supported with SDESystems and/or JumpSystems. Needed components in MT: - [ ]...

This would also allow users to put together stoichiometry matrix + rate law vector representations for the ODE rhs... I think such a network representation would be useful for many...