asteroid_spinprops 1.3.4__tar.gz → 1.3.6__tar.gz

{asteroid_spinprops-1.3.4 → asteroid_spinprops-1.3.6}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: asteroid_spinprops
-Version: 1.3.4
+Version: 1.3.6
 Summary: Collection of tools used for fitting sHG1G2 and SOCCA photometric models to sparse asteroid photometry
 License: MIT
 Author: Odysseas

{asteroid_spinprops-1.3.4 → asteroid_spinprops-1.3.6}/asteroid_spinprops/ssolib/modelfit.py

@@ -237,7 +237,6 @@ def get_fit_params(
                         rms_landscape, factor=4, sigma=1.0
                     )
                     ny, nx = interp_vals.shape
-                    [rarange[0], rarange[-1], decrange[0], decrange[-1]]
                     ra_vals = np.linspace(rarange.min(), rarange.max(), nx)
                     dec_vals = np.linspace(decrange.min(), decrange.max(), ny)
                     ys, xs = utils.detect_local_minima(interp_vals)
@@ -247,16 +246,25 @@ def get_fit_params(
                     ra_init = ra_minima
                     dec_init = dec_minima
 
+                    # Add near-pole initialization points
                     ra_init = np.append(ra_init, 220)
                     ra_init = np.append(ra_init, 140)
 
-                    dec_init = np.append(ra_init, 70)
-                    dec_init = np.append(ra_init, -70)
+                    dec_init = np.append(dec_init, 70)
+                    dec_init = np.append(dec_init, -70)
+
+                    # Remove pairs at the parameter space border
+                    RA_MARGIN = 1.0  # degrees
+
+                    ra_mask = (ra_init > RA_MARGIN) & (ra_init < 360 - RA_MARGIN)
+
+                    ra_init = ra_init[ra_mask]
+                    dec_init = dec_init[ra_mask]
 
                 H_key = next(
                     (
                         f"H_{i}" for i in range(1, 7) if f"H_{i}" in shg1g2_params
-                    ),  # Harcoded number of filters? Grr...
+                    ),  # FIXME: Harcoded N of bands, won't throw error if N>6, but to be recosnidered
                     None,
                 )
 
@@ -324,10 +332,6 @@ def get_fit_params(
                         SOCCA_opt["prms"] = p_rms
                         SOCCA_opt["k_terms"] = k_val
                 except Exception:
-                    if time_me:
-                        t2 = time.time()
-                        etime = t2 - t1
-                        SOCCA_opt["invtime"] = etime
                     SOCCA_opt = {"Failed at SOCCA inversion": 1}
                 if time_me:
                     t2 = time.time()

{asteroid_spinprops-1.3.4 → asteroid_spinprops-1.3.6}/pyproject.toml

@@ -1,6 +1,6 @@
 [project]
 name = "asteroid_spinprops"
-version = "1.3.4"
+version = "1.3.6"
 description = "Collection of tools used for fitting sHG1G2 and SOCCA photometric models to sparse asteroid photometry"
 authors = [
     {name = "Odysseas",email = "odysseas.xenos@proton.me"}