@@ -146,38 +146,6 @@ function assemble_maj(js, maj::MassActionJump{U,Vector{Pair{V,W}},Vector{Pair{V2
146
146
return maj
147
147
end
148
148
149
- function assemble_maj_expr (js, maj:: MassActionJump{U,Vector{Pair{V,W}},Vector{Pair{V2,W2}}} ,
150
- statetoid, subber, invttype) where {U,V,W,V2,W2}
151
- sr = maj. scaled_rates
152
- if sr isa Operation
153
- pval = subber (sr). value
154
- elseif sr isa Variable
155
- pval = subber (sr ()). value
156
- else
157
- pval = maj. scaled_rates
158
- end
159
-
160
- rs = Vector {Pair{Int,W}} ()
161
- for (spec,stoich) in maj. reactant_stoch
162
- if iszero (spec)
163
- push! (rs, 0 => stoich)
164
- else
165
- push! (rs, statetoid[convert (Variable,spec)] => stoich)
166
- end
167
- end
168
- sort! (rs)
169
-
170
- ns = Vector {Pair{Int,W2}} ()
171
- for (spec,stoich) in maj. net_stoch
172
- iszero (spec) && error (" Net stoichiometry can not have a species labelled 0." )
173
- push! (ns, statetoid[convert (Variable,spec)] => stoich)
174
- end
175
- sort! (ns)
176
-
177
- maj = MassActionJump (convert (invttype, pval), rs, ns, scale_rates = false )
178
- return maj
179
- end
180
-
181
149
"""
182
150
```julia
183
151
function DiffEqBase.DiscreteProblem(sys::JumpSystem, u0map, tspan,
0 commit comments