print 'Initial position:', src.pos x,y = src.pos.x, src.pos.y tr.freezeParam('images') tr.printThawedParams() for step in range(50): dlnp,X,alpha = tr.optimize(priors=False, shared_params=False, alphas=alphas) print 'dlnp', dlnp print 'pos', src.pos.x, src.pos.y print 'Delta position:', src.pos.x - x, src.pos.y - y #if dlnp < 0.1: if dlnp == 0.: break print 'Final position:', src.pos pos = src.getPosition() stars.xfit[i] = stars.x0[i] + pos.x stars.yfit[i] = stars.y0[i] + pos.y mod = tr.getModelImage(0) mods2.append(mod) subtim = tr.images[0] resids2.append((subtim.data - mod) * (subtim.inverr > 0)) chis2.append((subtim.data - mod) * subtim.inverr) plt.clf() for i,mod in enumerate(mods2): plt.subplot(R, C, i+1) dimshow(mod, **ima) plt.suptitle('Fit Models: ' + tim.name) ps.savefig()