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()
