JumpProcesses.jl
JumpProcesses.jl copied to clipboard
Build and simulate jump equations like Gillespie simulations and jump diffusions with constant and state-dependent rates and mix with differential equations and scientific machine learning (SciML)
`CartesianGridRej` performs significantly better than `CartesianGridIter`.
This is in response to a discourse question: https://discourse.julialang.org/t/jumpproblem-hopping-constants-callback/86045?u=robertgregg Essentially, it would be nice to be able to modify `hopping_constants` in a callback to simulate a moving boundary separates two...
This summer I am working on implementing spatial SSAs. There are currently two issues open about it: https://github.com/SciML/DiffEqJump.jl/issues/107 and https://github.com/SciML/DiffEqJump.jl/issues/130. The first PR https://github.com/SciML/DiffEqJump.jl/pull/183 implements the NSM, adds two tests...
Extends the VariableRateJump interface to allow for user-defined bounds on time-dependent jump rates and look-ahead horizons for which the bound is valid. Implements a fictitious jump rejection algorithm based on...
I am creating a new issue to help track the work on migrating some of the `DiffEqJump.jl` documentation to `Documenter.jl.` In talking with @isaacsas , he indicated that he wanted...
This PR updates the specific docstrings below. This is a first PR, to also confirm that the template and updates look good. Once everything looks good, I can continue to...
```julia prob = ODEProblem((du, u, p, t) -> du .= 0, ArrayPartition([0.]), (0, 10.)) jump = VariableRateJump((u, p, t) -> 1, int -> (int.u += 1)) jprob = JumpProblem(prob, Direct(),...
Here is the code: ``` using DifferentialEquations function f(du,u,p,t) du[1] = u[1] end function affect!(integrator) #print(typeof(integrator.u)) integrator.u = [0.0, 1.0] end function main() x₀ = [1.0, 2.0] prob = ODEProblem(f,...
In https://github.com/SciML/DiffEqJump.jl/pull/225 I couldn't consistently use the user-provided rng to JumpProblem as some places only have integrators as input (for example in callbacks). It would be nice to have a...
The DiffEq ecosystem already has some support for DDEs, so I thought extending to DSSAs would be a natural extension for DiffEqJump. Is there any interest in this?