-
-
Notifications
You must be signed in to change notification settings - Fork 77
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Detailed balance #965
Detailed balance #965
Conversation
This needs to be updated to master. |
@vyudu when you have a chance can you update to master? (No rush though as I’m traveling this week as we discussed.) |
bump |
Sorry about the delay - will be active again in a week or so |
This one should be ready for review @isaacsas |
src/network_analysis.jl
Outdated
if deficiency(rs) > 0 | ||
rxn_idxs = [edgeindex(D, Graphs.src(e), Graphs.dst(e)) for e in spanning_forest] | ||
S_F = netstoichmat(rs)[:, rxn_idxs] | ||
sols = nullspace(S_F) |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Should this use ModelingToolkit's integer nullspace?
Code and tests for testing whether a reaction networks with a given set of rate parameters is detailed balanced.