Performance advice for Jump simulation tutorial#750
Conversation
|
Also, @isaacsas, If you know a good reference book/article on Jump simulation it would be useful to know (so I could put a reference to it here, for further reading). |
|
I think for bio readers this is a nice one that is focused on algorithms: Luca Marchetti, Corrado Priami, Vo Hong Thanh, |
|
I tried writing the tutorial assuming that someone wasn't familiar with rates, propensities, and rate laws, and hence basically just used rate instead of propensities. I don't min making it more rigorous, but then I want to add a proper introduction of what these things are. I think if we just start using propensity people not familiar with it will get confused really quickly. (long term I want a better introduction to this somewhere in the docs as well, just so that it can be cross-referenced in situations like this) |
|
Is it worth mentioning small case static array stuff? |
|
We should probably benchmark that case and see if there is a speed up. I don’t think we’ve ever done so. |
|
I have, you can see it here: https://github.com/sdwfrost/Gillespie.jl?tab=readme-ov-file#benchmarks It's a bit old but it works well for small problems. In fact, for the smallest of small problems we should probably make a different codegen approach, but that's a whole different story. |
baafeb4 to
5ffc323
Compare
Co-authored-by: Sam Isaacson <isaacsas@users.noreply.github.com>
Co-authored-by: Sam Isaacson <isaacsas@users.noreply.github.com>
Co-authored-by: Sam Isaacson <isaacsas@users.noreply.github.com>
Co-authored-by: Sam Isaacson <isaacsas@users.noreply.github.com>
Co-authored-by: Sam Isaacson <isaacsas@users.noreply.github.com>
Co-authored-by: Sam Isaacson <isaacsas@users.noreply.github.com>
Co-authored-by: Sam Isaacson <isaacsas@users.noreply.github.com>
Co-authored-by: Sam Isaacson <isaacsas@users.noreply.github.com>
Co-authored-by: Sam Isaacson <isaacsas@users.noreply.github.com>
Co-authored-by: Sam Isaacson <isaacsas@users.noreply.github.com>
Co-authored-by: Sam Isaacson <isaacsas@users.noreply.github.com>
5ffc323 to
bdf964a
Compare
Adds a single tutorial, providing some advice on Jump simulation performance:
save_positionsandsaveat.k, E + S --> E + Pis better thank*E, S --> P.RSSAis probably better thanDirect(andRSSACRshould be preferred for big models).