look at E-M fit of PSF -- try initializing from center of patch