Skip to content

Commit 22b7b8e

Browse files
committed
fix substituter calls
1 parent ead3edb commit 22b7b8e

File tree

2 files changed

+2
-2
lines changed

2 files changed

+2
-2
lines changed

src/systems/jumps/jumpsystem.jl

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -233,7 +233,7 @@ function DiffEqJump.JumpProblem(js::JumpSystem, prob, aggregator; kwargs...)
233233
# handling parameter substition and empty param vecs
234234
p = (prob.p == DiffEqBase.NullParameters()) ? Operation[] : prob.p
235235
parammap = map((x,y)->Pair(x(),y), parameters(js), p)
236-
subber = substituter(first.(parammap), last.(parammap))
236+
subber = substituter(parammap)
237237

238238
majs = MassActionJump[assemble_maj(js, j, statetoid, subber, invttype) for j in eqs.x[1]]
239239
crjs = ConstantRateJump[assemble_crj(js, j, statetoid) for j in eqs.x[2]]

test/reactionsystem.jl

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -117,7 +117,7 @@ jumps[20] = VariableRateJump((u,p,t) -> p[20]*t*u[1]*binomial(u[2],2)*u[3], inte
117117
statetoid = Dict(convert(Variable,state) => i for (i,state) in enumerate(states(js)))
118118
parammap = map((x,y)->Pair(x(),y),parameters(js),pars)
119119
for i = 1:14
120-
maj = MT.assemble_maj(js, js.eqs[i], statetoid, ModelingToolkit.substituter(first.(parammap), last.(parammap)),eltype(pars))
120+
maj = MT.assemble_maj(js, js.eqs[i], statetoid, ModelingToolkit.substituter(parammap),eltype(pars))
121121
@test abs(jumps[i].scaled_rates - maj.scaled_rates) < 100*eps()
122122
@test jumps[i].reactant_stoch == maj.reactant_stoch
123123
@test jumps[i].net_stoch == maj.net_stoch

0 commit comments

Comments
 (0)