We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
There was an error while loading. Please reload this page.
2 parents 668866a + b210c02 commit 15a7d07Copy full SHA for 15a7d07
test/jumpsystem.jl
@@ -197,3 +197,25 @@ OBS2 = OBS
197
@test isequal(OBS2, @nonamespace js5.OBS)
198
@unpack OBS = js5
199
@test isequal(OBS2, OBS)
200
+
201
+# test to make sure dep graphs are correct
202
+let
203
+ # A + 2X --> 3X
204
+ # 3X --> A + 2X
205
+ # B --> X
206
+ # X --> B
207
+ @variables t A(t) X(t) B(t)
208
+ jumps = [MassActionJump(1.0, [A => 1, X => 2], [A => -1, X => 1]),
209
+ MassActionJump(1.0, [X => 3], [A => 1, X => -1]),
210
+ MassActionJump(1.0, [B => 1], [B => -1, X => 1]),
211
+ MassActionJump(1.0, [X => 1], [B => 1, X => -1])]
212
+ @named js = JumpSystem(jumps, t, [A, X, B], [])
213
+ jdeps = asgraph(js)
214
+ vdeps = variable_dependencies(js)
215
+ vtoj = jdeps.badjlist
216
+ @test vtoj == [[1], [1, 2, 4], [3]]
217
+ jtov = vdeps.badjlist
218
+ @test jtov == [[1, 2], [1, 2], [2, 3], [2, 3]]
219
+ jtoj = eqeq_dependencies(jdeps, vdeps).fadjlist
220
+ @test jtoj == [[1, 2, 4], [1, 2, 4], [1, 2, 3, 4], [1, 2, 3, 4]]
221
+end
0 commit comments