def test_scalar_variable(): ''' Test the correct handling of scalar variables ''' tau = 10*ms G = NeuronGroup(10, '''E_L : volt (shared) s2 : 1 (shared) dv/dt = (E_L - v) / tau : volt''') # Setting should work in these ways G.E_L = -70*mV assert_allclose(G.E_L[:], -70*mV) G.E_L[:] = -60*mV assert_allclose(G.E_L[:], -60*mV) G.E_L = 'E_L + s2*mV - 10*mV' assert_allclose(G.E_L[:], -70*mV) G.E_L[:] = '-75*mV' assert_allclose(G.E_L[:], -75*mV) net = Network(G) net.run(defaultclock.dt)