Catalyst.jl Files

Catalyst v17.0.0

Diff since v16.5.0

Breaking changes

Catalyst 17.0 updates to JumpProcesses 10, and requires ModelingToolkitBase 1.78 or later, the first release that supports it.

• The random number generator moves from JumpProblem to solve/init. Passing rng to JumpProblem, or to a HybridProblem with jumps, now throws an ArgumentError. Write solve(jprob, SSAStepper(); rng = StableRNG(1234)) instead. solve and init also accept a seed.
• Jumps draw their random numbers from the generator passed to solve. Creating a new generator with the same seed for every solve therefore repeats one trajectory. To draw independent, reproducible samples, create one generator and pass it to every solve.
• A JumpProblem supports one solve at a time. Solves reuse the problem's jump state instead of copying it, and problems created with remake share it. For parallel simulations, use an EnsembleProblem or a deepcopy of the problem per thread or task.
• Mass-action rate constants are recomputed from the current parameters at every solve/init and every reset_aggregated_jumps!. Callbacks written directly with DiffEqCallbacks or SciMLBase that change parameters must still call reset_aggregated_jumps!(integrator), whose update_jump_params keyword is removed. jprob.massaction_jump.scaled_rates no longer holds the rate constants.

See HISTORY.md (github.com) for details and examples, and the JumpProcesses [HISTORY.md](https://github.com/SciML/JumpProcesses.jl/blob/master/HISTORY.md#100 for the complete list of JumpProcesses 10 changes.

Merged pull requests:

• Support JumpProcesses 10 in tests, benchmarks, and docs (#1572) (@isaacsas)