def lnLikelihood(image, flux, x, y, psf, invvar):
star = MixtureOfGaussians(np.array([
flux,
]), np.array([x, y]), np.array([
0.,
])) # delta-function
cstar = star.convolve(psf)
modelimage = cstar.evaluate_grid(-5, 7, -5, 7)
return -0.5 * np.sum(invvar * (image - modelimage)**2)
def lnLikelihood(image, flux, x, y, psf, invvar):
star = MixtureOfGaussians(np.array([flux, ]), np.array([x, y]), np.array([0., ])) # delta-function
cstar = star.convolve(psf)
modelimage = cstar.evaluate_grid(-5, 7, -5, 7)
return -0.5 * np.sum(invvar * (image - modelimage)**2)
