asteroid_spinprops 1.3.9__tar.gz → 1.4.1__tar.gz

{asteroid_spinprops-1.3.9 → asteroid_spinprops-1.4.1}/PKG-INFO

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

asteroid_spinprops-1.4.1/asteroid_spinprops/ssolib/.ruff_cache/.gitignore

@@ -0,0 +1,2 @@
+# Automatically created by ruff.
+*

asteroid_spinprops-1.4.1/asteroid_spinprops/ssolib/.ruff_cache/CACHEDIR.TAG

@@ -0,0 +1 @@
+Signature: 8a477f597d28d172789f06886806bc55

{asteroid_spinprops-1.3.9 → asteroid_spinprops-1.4.1}/asteroid_spinprops/ssolib/modelfit.py

@@ -296,10 +296,11 @@ def get_fit_params(
                     G1, G2 = 0.15, 0.15
 
                 a_b, a_c = shg1g2_params["a_b"], shg1g2_params["a_c"]
-                if not (1 <= a_b <= 5 and 1 <= a_c <= 5):
+                if (not (1 <= a_b <= 5 and 1 <= a_c <= 5)) or np.isclose(
+                    a_b, a_c, rtol=1e-6, atol=1e-9
+                ):
                     a_b = 1.05
                     a_c = 1.5
-
                 for ra, dec in zip(ra_init, dec_init):
                     for period_sc in period_scan:
                         p_in = [
@@ -350,6 +351,9 @@ def get_fit_params(
                             else:
                                 SOCCA_opt["Period_class"] = 0  # Alias
                             SOCCA_opt["Nbs"] = Nbs
+                            SOCCA_opt["f_obs"] = f_obs
+                            SOCCA_opt["signal_peaks"] = signal_peaks
+                            SOCCA_opt["window_peaks"] = window_peaks
                         except Exception:
                             SOCCA_opt["Period_class"] = -1  # Classification error
                     if period_in is None:
@@ -447,6 +451,8 @@ def get_fit_params(
                 if terminator:
                     if remap:
                         p0in = np.concatenate((p0[3:], p0[:3]))
+                        if remap_kwargs["use_angles"] is True:
+                            p0in = np.insert(p0in, 0, 1.0)  # Initiate rho at 1
                         p0_latent = parameter_remapping(
                             p0in, physical_to_latent=True, **remap_kwargs
                         )

{asteroid_spinprops-1.3.9 → asteroid_spinprops-1.4.1}/asteroid_spinprops/ssolib/utils.py

@@ -318,14 +318,14 @@ def oblateness(a_b, a_c):
     Parameters
     ----------
     a_b : float or array_like
-        Intermediate-to-long axis ratio (b/a).
+        Intermediate-to-long axis ratio (a/b).
     a_c : float or array_like
-        Short-to-long axis ratio (c/a).
+        Short-to-long axis ratio (a/c).
 
     Returns
     -------
     float or array_like
-        Oblateness measure defined as 0.5 * (a_b / a_c + 1 / a_b).
+        Oblateness
     """
     return 1 / 2 * a_b / a_c + 1 / 2 * 1 / a_b
 

{asteroid_spinprops-1.3.9 → asteroid_spinprops-1.4.1}/pyproject.toml

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