cont-r.py

from sector import sector
from finger import finger
from f2s import f2s
from ucont import ucont

# s = sector.load('objects/step.object')

f = finger.load('fingers/ic/e-0.1.finger')

# determine the ``right'' b value
branch, f = ucont(f, 'dummy', 'b', 'vy',     1,    1.0)

s = f2s(f, 1.0e6, {'dummy':0.0})
# s = f2s(f, 10.0**5.0, {'dummy': 0.0})

# determine the ``right'' epsilon value

branch, s = ucont(s, 'dummy', 'e', 'omega',     1,    1.0)

def z0(x, p):
  return p - 10.0

branch, s = ucont(s, 'r', 'e', 'omega', 1000, -100.0, [z0])