def METHOD_NAME(self, data, g, velx, tvals, beta = -.9, nIterates = 1):
    for i in range(nIterates):
        lh0 = ndimage.gaussian_filter(self.dlHd2(data, g, velx, tvals), 50)
        lh1 = ndimage.gaussian_filter(self.dlHd2(data, g, velx + .1, tvals), 50)
        upd = beta*.1*lh0/(lh1-lh0 + .001)
        #pl.plot(lh1 - lh0)
        #plot(lh1)
        #pl.plot(upd)
        velx[:] = ndimage.gaussian_filter(velx + upd, 50)