fix: flux = 0 was not handled correctly

galsci · Dec 23, 2024 · 21dc598 · 21dc598
1 parent 8d160f4
commit 21dc598
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/src/pysm3/models/catalog.py b/src/pysm3/models/catalog.py
@@ -65,8 +65,8 @@ def flux2amp(flux, fwhm):
     sigma = fwhm2sigma(fwhm)
     amp = flux / (2 * np.pi * sigma**2)
     # sim_objects fails if amp is zero
-    c = 1e-9 # minimum amplitude
-    amp[np.logical_and(amp < c, amp > 0)] = c
+    c = 1e-9 # clip
+    amp[np.logical_and(amp < c, amp >= 0)] = c
     amp[np.logical_and(amp > -c, amp < 0)] = -c
     return amp