def run2():
pygsl.set_debug_level(10)
x = 2.0
step = odeiv.step_rk8pd(2, func, None, x)
control = odeiv.control_y_new(step, 1e-6, 1e-6)
evolve = odeiv.evolve(step, control)
pass
def run2():
pygsl.set_debug_level(10)
x = 2.0
step = odeiv.step_rk8pd(2, func, None, x)
control = odeiv.control_y_new(step, 1e-6, 1e-6)
evolve = odeiv.evolve(step, control)
pass
def run():
mu = 2.0
dimension = 2
# The different possible steppers for the function
# Comment your favourite one out to test it.
#stepper = odeiv.step_rk2
#stepper = odeiv.step_rk4
#stepper = odeiv.step_rkf45
#stepper = odeiv.step_rkck
stepper = odeiv.step_rk8pd
#stepper = odeiv.step_rk2imp
#stepper = odeiv.step_rk4imp
#stepper = odeiv.step_gear1
#stepper = odeiv.step_gear2
#stepper = odeiv.step_bsimp
step = stepper(dimension, func, jac, mu)
# All above steppers exept the odeiv.step_bsimp (Buerlisch - Stoer Improved
# method can calculate without an jacobian ...
# step = stepper(dimension, func)
#print step
#print step.type()
control = odeiv.control_y_new(step, 1e-6, 1e-6)
evolve = odeiv.evolve(step, control)
print "# Using stepper %s with order %d" %(step.name(), step.order())
print "# Using Control ", control.name()
print "# %5s %9s %9s %9s " % ("iter", "t", "y[0]", "y[1]")
h = 1.
tstart = 0.0
t1 = 5
ystart = (1.0, 0.0)
t = tstart
y = ystart
stamp = time.time()
nsteps = 1000
for i in xrange(nsteps):
if t >= t1:
break
t, h, y = evolve.apply(t, t1, h, y, mu)
#print t
##y = y[-1]
print " %5d % 10.6f % 10.6f % 10.6f " %(i, t, y[0], y[1])
else:
raise ValueError, "Maximum number of steps exceeded!"
print "Needed %f seconds" %( time.time() - stamp,)
print " % 10.6f % 10.6f % 10.6f " %(t, y[0], y[1])
stamp = time.time()
def run():
mu = 2.0
dimension = 2
# The different possible steppers for the function
# Comment your favourite one out to test it.
#stepper = odeiv.step_rk2
#stepper = odeiv.step_rk4
#stepper = odeiv.step_rkf45
#stepper = odeiv.step_rkck
stepper = odeiv.step_rk8pd
#stepper = odeiv.step_rk2imp
#stepper = odeiv.step_rk4imp
#stepper = odeiv.step_gear1
#stepper = odeiv.step_gear2
#stepper = odeiv.step_bsimp
step = stepper(dimension, func, jac, mu)
# All above steppers exept the odeiv.step_bsimp (Buerlisch - Stoer Improved
# method can calculate without an jacobian ...
# step = stepper(dimension, func)
#print step
#print step.type()
control = odeiv.control_y_new(step, 1e-6, 1e-6)
evolve = odeiv.evolve(step, control)
print "# Using stepper %s with order %d" % (step.name(), step.order())
print "# Using Control ", control.name()
print "# %5s %9s %9s %9s " % ("iter", "t", "y[0]", "y[1]")
h = 1.
tstart = 0.0
t1 = 5
ystart = (1.0, 0.0)
t = tstart
y = ystart
stamp = time.time()
nsteps = 1000
for i in xrange(nsteps):
if t >= t1:
break
t, h, y = evolve.apply(t, t1, h, y, mu)
#print t
##y = y[-1]
print " %5d % 10.6f % 10.6f % 10.6f " % (i, t, y[0], y[1])
else:
raise ValueError, "Maximum number of steps exceeded!"
print "Needed %f seconds" % (time.time() - stamp, )
print " % 10.6f % 10.6f % 10.6f " % (t, y[0], y[1])
stamp = time.time()
