From 65a136a7a49f200ca66703ee931e97d172c88278 Mon Sep 17 00:00:00 2001 From: Mike Jarvis Date: Fri, 16 Oct 2026 13:58:29 -0400 Subject: [PATCH 1/2] Replace broken link with link to Moffat's paper (#1368) --- galsim/moffat.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/galsim/moffat.py b/galsim/moffat.py index 189d1b5605..d8b32bd936 100644 --- a/galsim/moffat.py +++ b/galsim/moffat.py @@ -43,7 +43,7 @@ class Moffat(GSObject): For more information, refer to - http://home.fnal.gov/~neilsen/notebook/astroPSF/astroPSF.html + https://ui.adsabs.harvard.edu/abs/1969A%26A.....3..455M/abstract A Moffat can be initialized using one (and only one) of three possible size parameters: ``scale_radius``, ``fwhm``, or ``half_light_radius``. Exactly one of these three is required. From c7b5b354203dc1f032ae36de0e7ab96c895f0c34 Mon Sep 17 00:00:00 2001 From: Mike Jarvis Date: Thu, 15 Oct 2026 20:59:48 -0400 Subject: [PATCH 2/2] Use zeros rather than empty for Sigma, since not all array values are filled in. (#1365) --- galsim/real.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/galsim/real.py b/galsim/real.py index 8f12ad23a4..78cebd4e35 100644 --- a/galsim/real.py +++ b/galsim/real.py @@ -1265,7 +1265,7 @@ def _initialize(self, imgs, bands, xis, PSFs, # Note: put NSED axis last, since significantly faster to compute them this way, # even though we eventually convert to images which are strided in this format. coef = np.zeros((nk, nk, NSED), dtype=np.complex128) - Sigma = np.empty((nk, nk, NSED, NSED), dtype=np.complex128) + Sigma = np.zeros((nk, nk, NSED, NSED), dtype=np.complex128) # Solve the weighted linear least squares problem for each Fourier mode. This is # effectively a constrained chromatic deconvolution. Take advantage of symmetries.