tsdynamics 2.2.1__tar.gz → 2.2.2__tar.gz

{tsdynamics-2.2.1 → tsdynamics-2.2.2}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: tsdynamics
-Version: 2.2.1
+Version: 2.2.2
 Summary: Compiled ODE/DDE integration, discrete maps, and chaos analysis for dynamical systems.
 Project-URL: Source, https://github.com/El3ssar/TSDynamics
 Project-URL: Documentation, https://el3ssar.github.io/TSDynamics/

{tsdynamics-2.2.1 → tsdynamics-2.2.2}/src/tsdynamics/__init__.py

@@ -48,7 +48,7 @@ from .systems import continuous as _continuous
 from .systems import discrete as _discrete
 
 # Single source of truth for the package version; rewritten by python-semantic-release.
-__version__ = "2.2.1"
+__version__ = "2.2.2"
 
 # Re-export every system class at the top level so ``from tsdynamics import Lorenz``
 # works without users having to remember which submodule a system lives in.

{tsdynamics-2.2.1 → tsdynamics-2.2.2}/src/tsdynamics/systems/continuous/chaotic_attractors.py

@@ -547,6 +547,9 @@ class MultiChua(ContinuousSystem):
 class Duffing(ContinuousSystem):
     params = {"alpha": 1.0, "beta": -1.0, "delta": 0.1, "gamma": 0.35, "omega": 1.4}
     dim = 3
+    # The explicit default (dopri5) fails to integrate this system; an implicit
+    # solver handles it robustly, so make that the default.
+    _default_method = "LSODA"
 
     @staticmethod
     def _equations(Y, t, *, alpha, beta, delta, gamma, omega):
@@ -889,6 +892,7 @@ class SprottK(ContinuousSystem):
 class SprottL(ContinuousSystem):
     params = {"a": 0.9, "b": 3.9}
     dim = 3
+    _default_method = "LSODA"  # explicit default solver fails; use an implicit one
 
     @staticmethod
     def _equations(Y, t, *, a, b):
@@ -975,6 +979,7 @@ class SprottO(ContinuousSystem):
 class SprottP(ContinuousSystem):
     params = {"a": 2.7}
     dim = 3
+    _default_method = "LSODA"  # explicit default solver fails; use an implicit one
 
     @staticmethod
     def _equations(Y, t, *, a):
@@ -1080,6 +1085,7 @@ class SprottMore(ContinuousSystem):
 class SprottJerk(ContinuousSystem):
     params = {"mu": 2.017}
     dim = 3
+    _default_method = "LSODA"  # explicit default solver fails; use an implicit one
 
     @staticmethod
     def _equations(Y, t, *, mu):

{tsdynamics-2.2.1 → tsdynamics-2.2.2}/src/tsdynamics/systems/continuous/chem_bio_systems.py

@@ -37,6 +37,9 @@ class Oregonator(ContinuousSystem):
         "epsilon": 1e-2,
     }
     dim = 3  # Three variables: X, Y, Z (reduced forms of the chemical species)
+    # Classic stiff system (Field–Noyes); an explicit solver cannot integrate
+    # it, so default to an implicit one.
+    _default_method = "LSODA"
 
     @staticmethod
     def _equations(Y, t, *, q, f, mu, epsilon):

{tsdynamics-2.2.1 → tsdynamics-2.2.2}/tests/test_ode_systems.py

@@ -134,10 +134,9 @@ def test_ode_full_integration_sweep(ode_entry) -> None:
     # get a fixed (seeded) draw so the sweep is reproducible, not flaky.
     np.random.seed(zlib.crc32(ode_entry.name.encode()) & 0xFFFFFFFF)
     ic = sys.resolve_ic(None)
-    if ode_entry.name == "Oregonator":  # stiff — needs an implicit solver
-        traj = sys.integrate(ic=ic, final_time=1.0, dt=0.1, method="LSODA", rtol=1e-8, atol=1e-10)
-    else:
-        traj = sys.integrate(ic=ic, final_time=2.0, dt=0.1, rtol=1e-5, atol=1e-7)
+    # Each system carries its own _default_method (stiff systems default to an
+    # implicit solver), so the plain default path must integrate them all.
+    traj = sys.integrate(ic=ic, final_time=2.0, dt=0.1, rtol=1e-5, atol=1e-7)
     assert traj.y.shape[1] == sys.dim
     assert np.all(np.isfinite(traj.y))