Skip to content

Commit 091c95b

Browse files
committed
started
1 parent 8bf2adc commit 091c95b

File tree

1 file changed

+30
-0
lines changed

1 file changed

+30
-0
lines changed

src/systems/jumps/jumpsystem.jl

Lines changed: 30 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -1,5 +1,35 @@
11
JumpType = Union{VariableRateJump, ConstantRateJump, MassActionJump}
22

3+
"""
4+
$(TYPEDEF)
5+
6+
A system of jumps processes.
7+
8+
# Fields
9+
$(FIELDS)
10+
11+
# Example
12+
13+
```
14+
using ModelingToolKit, DiffEqJump
15+
16+
@parameters β γ t
17+
@variables S I R
18+
rate₁ = β*S*I
19+
affect₁ = [S ~ S - 1, I ~ I + 1]
20+
rate₂ = γ*I+t
21+
affect₂ = [I ~ I - 1, R ~ R + 1]
22+
j₁ = ConstantRateJump(rate₁,affect₁)
23+
j₂ = VariableRateJump(rate₂,affect₂)
24+
j₃ = MassActionJump(2*β+γ, [R => 1], [S => 1, R => -1])
25+
js = JumpSystem([j₁,j₂,j₃], t, [S,I,R], [β,γ])
26+
u₀map = [S => 999, I => 1, R => 0]
27+
parammap = [β => .1/1000, γ => .01]
28+
dprob = DiscreteProblem(js2, u₀map, tspan, parammap)
29+
jprob = JumpProblem(js2, dprob, Direct(), save_
30+
31+
```
32+
"""
333
struct JumpSystem{U <: ArrayPartition} <: AbstractSystem
434
eqs::U
535
iv::Variable

0 commit comments

Comments
 (0)