... for real.

Author: cschreib

src/egg-gennoise.cpp

@@ -112,18 +112,13 @@ int phypp_main(int argc, char* argv[]) {
 
     // Compute PSF FWHM
     double fwhm; {
-        vec1u idok = where(psf > 0.2*max(psf));
         vec2d r2 = generate_img(psf.dims, [&](double ix, double iy) {
             return sqr(ix - hsize) + sqr(iy - hsize);
         });
 
-        auto res = linfit(log(psf[idok]/max(psf)), 1, r2[idok]);
-        if (res.params[0] >= 0) {
-            error("could not estimate FWHM of the PSF");
-            return 1;
-        }
-
-        fwhm = 2*sqrt(-log(2)/res.params[0]);
+        vec1u idok = where(psf > 0.2*max(psf) && r2 > 0);
+        double res = mean(log(psf[idok]/max(psf))/r2[idok]);
+        fwhm = 2*sqrt(-log(2)/res);
 
         if (verbose) note("FWHM=", fwhm, " pixels");
     }