1010
1111class TestIntegratorRunnerForODEs (TestCase ):
1212 def test_ode (self ):
13+
1314 sigma = 10
1415 beta = 8 / 3
1516 rho = 28
@@ -21,16 +22,16 @@ def lorenz(x, y, z, t):
2122 dz = x * y - beta * z
2223 return dx , dy , dz
2324
24- runner = bp .integrators . IntegratorRunner (lorenz , monitors = ['x' , 'y' , 'z' ], inits = [1. , 1. , 1. ])
25+ runner = bp .IntegratorRunner (lorenz , monitors = ['x' , 'y' , 'z' ], inits = [1. , 1. , 1. ])
2526 runner .run (100. )
2627 fig = plt .figure ()
2728 fig .add_subplot (111 , projection = '3d' )
2829 plt .plot (runner .mon .x [:, 0 ], runner .mon .y [:, 0 ], runner .mon .z [:, 0 ], )
2930 plt .show ()
3031
31- runner = bp .integrators . IntegratorRunner (lorenz ,
32- monitors = ['x' , 'y' , 'z' ],
33- inits = [1. , (1. , 0. ), (1. , 0. )])
32+ runner = bp .IntegratorRunner (lorenz ,
33+ monitors = ['x' , 'y' , 'z' ],
34+ inits = [1. , (1. , 0. ), (1. , 0. )])
3435 runner .run (100. )
3536 for i in range (2 ):
3637 fig = plt .figure ()
@@ -47,7 +48,7 @@ def test_ode2(self):
4748 dw = lambda w , t , V : (V + a - b * w ) / tau
4849 fhn = bp .odeint (bp .JointEq ([dV , dw ]), method = 'rk4' , dt = 0.1 )
4950
50- runner = bp .integrators . IntegratorRunner (fhn , monitors = ['V' , 'w' ], inits = [1. , 1. ])
51+ runner = bp .IntegratorRunner (fhn , monitors = ['V' , 'w' ], inits = [1. , 1. ])
5152 runner .run (100. , args = dict (Iext = 1.5 ))
5253 bp .visualize .line_plot (runner .mon .ts , runner .mon ['V' ], legend = 'V' )
5354 bp .visualize .line_plot (runner .mon .ts , runner .mon ['w' ], legend = 'w' , show = True )
@@ -61,9 +62,9 @@ def test_ode3(self):
6162 fhn = bp .odeint (bp .JointEq ([dV , dw ]), method = 'rk4' , dt = 0.1 )
6263
6364 Iext , duration = bp .inputs .section_input ([0. , 1. , 0.5 ], [200 , 500 , 200 ], return_length = True )
64- runner = bp .integrators . IntegratorRunner (fhn ,
65- monitors = ['V' , 'w' ],
66- inits = [1. , 1. ])
65+ runner = bp .IntegratorRunner (fhn ,
66+ monitors = ['V' , 'w' ],
67+ inits = [1. , 1. ])
6768 runner .run (duration , dyn_args = dict (Iext = Iext ))
6869 bp .visualize .line_plot (runner .mon .ts , runner .mon ['V' ], legend = 'V' )
6970 bp .visualize .line_plot (runner .mon .ts , runner .mon ['w' ], legend = 'w' , show = True )
@@ -76,9 +77,9 @@ def test_ode_continuous_run(self):
7677 dw = lambda w , t , V : (V + a - b * w ) / tau
7778 fhn = bp .odeint (bp .JointEq ([dV , dw ]), method = 'rk4' , dt = 0.1 )
7879
79- runner = bp .integrators . IntegratorRunner (fhn ,
80- monitors = ['V' , 'w' ],
81- inits = [1. , 1. ])
80+ runner = bp .IntegratorRunner (fhn ,
81+ monitors = ['V' , 'w' ],
82+ inits = [1. , 1. ])
8283 Iext , duration = bp .inputs .section_input ([0. , 1. , 0.5 ], [200 , 200 , 200 ], return_length = True )
8384 runner .run (duration , dyn_args = dict (Iext = Iext ))
8485 bp .visualize .line_plot (runner .mon .ts , runner .mon ['V' ], legend = 'V' )
@@ -100,9 +101,9 @@ def test_ode_dyn_args(self):
100101 Iext , duration = bp .inputs .section_input ([0. , 1. , 0.5 ],
101102 [200 , 500 , 199 ],
102103 return_length = True )
103- runner = bp .integrators . IntegratorRunner (fhn ,
104- monitors = ['V' , 'w' ],
105- inits = [1. , 1. ])
104+ runner = bp .IntegratorRunner (fhn ,
105+ monitors = ['V' , 'w' ],
106+ inits = [1. , 1. ])
106107 with self .assertRaises (ValueError ):
107108 runner .run (duration + 1 , dyn_args = dict (Iext = Iext ))
108109
0 commit comments