# Tau-leaping Algorithm for Stochastic Simulation

**URL:** https://discourse.julialang.org/t/tau-leaping-algorithm-for-stochastic-simulation/128237
**Category:** New to Julia
**Tags:** question, package
**Created:** [April 20, 2025, 3:05am UTC](https://discourse.julialang.org/t/tau-leaping-algorithm-for-stochastic-simulation/128237 "2025-04-20T03:05:29Z")
**Posts on this page:** 5
**Page:** 1

<div class="post-metadata">

### Author: ![jameslawkc](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/jameslawkc/32/216416_2.png) [@jameslawkc](https://discourse.julialang.org/u/jameslawkc)
#### Post date: [April 20, 2025, 3:05am UTC](https://discourse.julialang.org/t/tau-leaping-algorithm-for-stochastic-simulation/128237/1 "2025-04-20T03:05:29Z")

</div>

I know that there is a package called JumpProcesses.jl but there is no option to implement the tau-leaping algorithm. Is there any existing package that can implement this method?

---

<div class="post-metadata">

### Author: ![arnold-c](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/arnold-c/32/222741_2.png) [@arnold-c](https://discourse.julialang.org/u/arnold-c)
#### Post date: [April 20, 2025, 3:30am UTC](https://discourse.julialang.org/t/tau-leaping-algorithm-for-stochastic-simulation/128237/2 "2025-04-20T03:30:33Z")

</div>

> [@jameslawkc](#):
>
> ere is a package called [JumpProcesses.jl](https://juliaregistries.github.io/General/packages/redirect_to_repo/JumpProcesses) but there is no option to implement the tau-leaping algorithm. Is there any existing package that can implement this method?

Welcome to the discourse @jameslawkc. I believe there is an implementation of Tau Leaping in JumpProcesses.jl. See [this section](https://github.com/SciML/JumpProcesses.jl/blob/a404c82e7146b4db3479742e61659109d3f684ac/src/simple_regular_solve.jl#L2) of the package code, and [this page](https://docs.sciml.ai/JumpProcesses/stable/jump_solve/#StochasticDiffEq.jl) of the docs (which also point to StochasticDiffEq.jl for an adaptive implementation of a Tau Leaping algorithm). I’m afraid I can’t help much with any implementation questions you may have as I’ve always just written my own versions, but hopefully [this test file](https://github.com/SciML/StochasticDiffEq.jl/blob/899da89d5069e68374f9edb6805830285a7acb98/test/tau_leaping.jl#L4) from StochasticDiffEq.jl can get you most of the way!

---

<div class="post-metadata">

### Author: ![jameslawkc](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/jameslawkc/32/216416_2.png) [@jameslawkc](https://discourse.julialang.org/u/jameslawkc)
#### Post date: [April 20, 2025, 3:55am UTC](https://discourse.julialang.org/t/tau-leaping-algorithm-for-stochastic-simulation/128237/3 "2025-04-20T03:55:13Z")

</div>

Thank you for the reply. I notice the tau leaping algorithm is called from StochasticDiffEq.jl. Is that the same algorithm in the context for jump processes? It just a bit weird we are calling this from a package that deals with continuous stochastic process.

---

<div class="post-metadata">

### Author: ![isaacsas](https://avatars.discourse-cdn.com/v4/letter/i/f6c823/32.png) [@isaacsas](https://discourse.julialang.org/u/isaacsas)
#### Post date: [April 20, 2025, 11:38am UTC](https://discourse.julialang.org/t/tau-leaping-algorithm-for-stochastic-simulation/128237/4 "2025-04-20T11:38:58Z")

</div>

Yes, the tau-leaping integrators just live in StochasticDiffEq but take inputs from JumpProcesses. This is because they reuse components from other integrators in StochasticDiffEq.

---

<div class="post-metadata">

### Author: ![ChrisRackauckas](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/chrisrackauckas/32/77_2.png) [@ChrisRackauckas](https://discourse.julialang.org/u/ChrisRackauckas)
#### Post date: [April 20, 2025, 3:24pm UTC](https://discourse.julialang.org/t/tau-leaping-algorithm-for-stochastic-simulation/128237/5 "2025-04-20T15:24:54Z")

</div>

Note that there should be a lot more work on it this summer, so things should improve in this space “soon”
