end = 2*nx//3
kf = IteratedMSKalmanFilter(distmesh, frame, flowframe, cuda)
#Now perturb the positions a bit...
kf.state.X = kf.state.X*1.2
kf.state.refresh()
rend = kf.state.renderer
cuda = kf.state.renderer.cudagl
kf.state.render()
state = kf.state
nx = nx
deltaX = -2
nF = video.shape[2]
nI = 10
count = 0
for i in range(500):
count += 1
print 'Frame %d' % count
#frame = video[:,:,i]
#mask = (frame > 0).astype(np.uint8)
#flowframe = flow[:,:,:,i]
kf.predict()
kf.state.refresh()
kf.state.render()
#Wait for user input:
raw_input("Press enter to continue, ctrl-d to quit")
