pynamicalsys 1.3.0__py3-none-any.whl → 1.3.1__py3-none-any.whl

pynamicalsys/__version__.py

@@ -28,7 +28,7 @@ version_tuple: VERSION_TUPLE
 commit_id: COMMIT_ID
 __commit_id__: COMMIT_ID
 
-__version__ = version = '1.3.0'
-__version_tuple__ = version_tuple = (1, 3, 0)
+__version__ = version = '1.3.1'
+__version_tuple__ = version_tuple = (1, 3, 1)
 
 __commit_id__ = commit_id = None

pynamicalsys/common/recurrence_quantification_analysis.py

@@ -91,7 +91,7 @@ class RTEConfig:
             raise ValueError("metric must be 'supremum', 'euclidean' or 'manhattan'")
 
 
-@njit(cache=True)
+@njit
 def _recurrence_matrix(
     arr: NDArray[np.float64], threshold: float, metric_id: int
 ) -> NDArray[np.uint8]:

pynamicalsys/common/utils.py

@@ -21,7 +21,7 @@ from numpy.typing import NDArray
 from numba import njit
 
 
-@njit(cache=True)
+@njit
 def qr(M: NDArray[np.float64]) -> Tuple[NDArray[np.float64], NDArray[np.float64]]:
     """
     Perform numerically stable QR decomposition using modified Gram-Schmidt with reorthogonalization.
@@ -92,7 +92,7 @@ def qr(M: NDArray[np.float64]) -> Tuple[NDArray[np.float64], NDArray[np.float64]
     return Q, R
 
 
-@njit(cache=True)
+@njit
 def householder_qr(
     M: NDArray[np.float64],
 ) -> Tuple[NDArray[np.float64], NDArray[np.float64]]:
@@ -184,7 +184,7 @@ def householder_qr(
     return Q, R
 
 
-@njit(cache=True)
+@njit
 def finite_difference_jacobian(
     u: NDArray[np.float64],
     parameters: NDArray[np.float64],

pynamicalsys/continuous_time/chaotic_indicators.py

@@ -25,7 +25,7 @@ from pynamicalsys.continuous_time.trajectory_analysis import step, evolve_system
 from pynamicalsys.continuous_time.numerical_integrators import rk4_step_wrapped
 
 
-# @njit(cache=True)
+# @njit
 def lyapunov_exponents(
     u: NDArray[np.float64],
     parameters: NDArray[np.float64],
@@ -136,7 +136,7 @@ def lyapunov_exponents(
         return [result]
 
 
-@njit(cache=True)
+@njit
 def SALI(
     u: NDArray[np.float64],
     parameters: NDArray[np.float64],
@@ -238,7 +238,7 @@ def SALI(
         return [[time, sali]]
 
 
-# @njit(cache=True)
+# @njit
 def LDI(
     u: NDArray[np.float64],
     parameters: NDArray[np.float64],

pynamicalsys/continuous_time/numerical_integrators.py

@@ -23,7 +23,7 @@ from numba import njit, prange
 from pynamicalsys.continuous_time.models import variational_equations
 
 
-@njit(cache=True)
+@njit
 def rk4_step(
     t: float,
     u: NDArray[np.float64],
@@ -43,7 +43,7 @@ def rk4_step(
     return u_next
 
 
-@njit(cache=True)
+@njit
 def variational_rk4_step(
     t: float,
     u: NDArray[np.float64],
@@ -115,7 +115,7 @@ RK45_B4 = np.array(
 )
 
 
-@njit(cache=True)
+@njit
 def rk45_step(t, u, parameters, equations_of_motion, time_step, atol=1e-6, rtol=1e-3):
     """Single adaptive step of RK45 (Dormand-Prince).
 
@@ -163,7 +163,7 @@ def rk45_step(t, u, parameters, equations_of_motion, time_step, atol=1e-6, rtol=
     return u5, t + time_step, time_step_new, accept
 
 
-@njit(cache=True)
+@njit
 def variational_rk45_step(
     t,
     u,
@@ -227,7 +227,7 @@ def variational_rk45_step(
     return u5, t + time_step, time_step_new, accept
 
 
-@njit(cache=True)
+@njit
 def rk4_step_wrapped(
     t: float,
     u: NDArray[np.float64],
@@ -265,7 +265,7 @@ def rk4_step_wrapped(
     return u_next, t_next, h_next, accept
 
 
-@njit(cache=True)
+@njit
 def rk45_step_wrapped(
     t: float,
     u: NDArray[np.float64],
@@ -302,7 +302,7 @@ def rk45_step_wrapped(
         )
 
 
-@njit(cache=True)
+@njit
 def estimate_initial_step(
     t0: float,
     u0: np.ndarray,

pynamicalsys/continuous_time/trajectory_analysis.py

@@ -24,7 +24,7 @@ from numpy.typing import NDArray
 from pynamicalsys.continuous_time.numerical_integrators import rk4_step_wrapped
 
 
-@njit(cache=True)
+@njit
 def step(
     time: np.float64,
     u: NDArray[np.float64],
@@ -68,7 +68,7 @@ def step(
     return u_new, time_new, time_step_new
 
 
-@njit(cache=True)
+@njit
 def evolve_system(
     u: NDArray[np.float64],
     parameters: NDArray[np.float64],
@@ -102,7 +102,7 @@ def evolve_system(
     return u
 
 
-@njit(cache=True)
+@njit
 def generate_trajectory(
     u: NDArray[np.float64],
     parameters: NDArray[np.float64],

pynamicalsys/discrete_time/dynamical_indicators.py

@@ -32,7 +32,7 @@ from pynamicalsys.discrete_time.trajectory_analysis import (
 from pynamicalsys.common.utils import qr, householder_qr, fit_poly, wedge_norm
 
 
-@njit(cache=True)
+@njit
 def lyapunov_1D(
     u: NDArray[np.float64],
     parameters: NDArray[np.float64],
@@ -120,7 +120,7 @@ def lyapunov_1D(
     return history if return_history else np.array([exponent / sample_size])
 
 
-@njit(cache=True)
+@njit
 def lyapunov_er(
     u: NDArray[np.float64],
     parameters: NDArray[np.float64],
@@ -246,7 +246,7 @@ def lyapunov_er(
         return aux_exponents, u_contig
 
 
-@njit(cache=True)
+@njit
 def maximum_lyapunov_er(
     u: NDArray[np.float64],
     parameters: NDArray[np.float64],
@@ -367,7 +367,7 @@ def maximum_lyapunov_er(
         return np.array([exponent / sample_size]), u_contig
 
 
-@njit(cache=True)
+@njit
 def lyapunov_qr(
     u: NDArray[np.float64],
     parameters: NDArray[np.float64],
@@ -701,7 +701,7 @@ def dig(
     return -np.log10(abs(WB0 - WB1))
 
 
-@njit(cache=True)
+@njit
 def SALI(
     u: NDArray[np.float64],
     parameters: NDArray[np.float64],
@@ -761,7 +761,7 @@ def SALI(
     - Uses QR decomposition to initialize orthonormal deviation vectors.
     - Computes both Parallel (PAI) and Antiparallel (AAI) Alignment Indices.
     - Early termination occurs if SALI < `tol` (indicating chaotic behavior).
-    - Optimized with `@njit(cache=True)` for performance.
+    - Optimized with `@njit` for performance.
     """
 
     np.random.seed(seed)  # For reproducibility
@@ -1047,7 +1047,7 @@ def GALI_k(
         return np.array([gali])
 
 
-@njit(cache=True)
+@njit
 def hurst_exponent(
     u: NDArray[np.float64],
     parameters: NDArray[np.float64],
@@ -1224,13 +1224,13 @@ def finite_time_hurst_exponent(
         return H_values
 
 
-@njit(cache=True)
+@njit
 def lyapunov_vectors():
     # ! To be implemented...
     pass
 
 
-@njit(cache=True)
+@njit
 def lagrangian_descriptors(
     u: NDArray[np.float64],
     parameters: NDArray[np.float64],

pynamicalsys/discrete_time/models.py

@@ -408,7 +408,7 @@ def logistic_map_jacobian(
 # ! ------------------ !
 
 
-@njit(cache=True)
+@njit
 def rulkov_map(
     u: NDArray[np.float64], parameters: Union[NDArray[np.float64], Sequence[float]]
 ) -> NDArray[np.float64]:
@@ -419,7 +419,7 @@ def rulkov_map(
     return np.array([x_new, y_new])
 
 
-@njit(cache=True)
+@njit
 def rulkov_map_jacobian(
     u: NDArray[np.float64],
     parameters: Union[NDArray[np.float64], Sequence[float]],

pynamicalsys/discrete_time/trajectory_analysis.py

@@ -22,7 +22,7 @@ from numba import njit, prange
 from numpy.typing import NDArray
 
 
-@njit(cache=True)
+@njit
 def iterate_mapping(
     u: NDArray[np.float64],
     parameters: NDArray[np.float64],
@@ -78,7 +78,7 @@ def iterate_mapping(
     return u
 
 
-@njit(cache=True)
+@njit
 def generate_trajectory(
     u: NDArray[np.float64],
     parameters: NDArray[np.float64],
@@ -331,7 +331,7 @@ def bifurcation_diagram(
         return param_values, results
 
 
-@njit(cache=True)
+@njit
 def period_counter(
     u: NDArray[np.float64],
     parameters: NDArray[np.float64],
@@ -415,7 +415,7 @@ def period_counter(
     return -1
 
 
-@njit(cache=True)
+@njit
 def rotation_number(
     u: Union[NDArray[np.float64], Sequence[float], float],
     parameters: Union[NDArray[np.float64], Sequence[float], float],
@@ -438,7 +438,7 @@ def rotation_number(
     return rn
 
 
-@njit(cache=True)
+@njit
 def escape_basin_and_time_entering(
     u: NDArray[np.float64],
     parameters: NDArray[np.float64],
@@ -515,7 +515,7 @@ def escape_basin_and_time_entering(
     return -1, max_time
 
 
-@njit(cache=True)
+@njit
 def escape_time_exiting(
     u: NDArray[np.float64],
     parameters: NDArray[np.float64],
@@ -584,7 +584,7 @@ def escape_time_exiting(
     return -1, max_time  # No escape
 
 
-@njit(cache=True)
+@njit
 def survival_probability(
     escape_times: NDArray[np.int32],
     max_time: np.int32,
@@ -659,7 +659,7 @@ def survival_probability(
     return t_values, survival_probs
 
 
-@njit(cache=True)
+@njit
 def is_periodic(
     u: NDArray[np.float64],
     parameters: NDArray[np.float64],
@@ -811,7 +811,7 @@ def scan_phase_space(
     return result
 
 
-@njit(cache=True)
+@njit
 def scan_symmetry_line(
     points: NDArray[np.float64],
     parameters: NDArray[np.float64],
@@ -1095,7 +1095,7 @@ def find_periodic_orbit(
     return periodic_orbit
 
 
-@njit(cache=True)
+@njit
 def eigenvalues_and_eigenvectors(
     u: NDArray[np.float64],
     parameters: NDArray[np.float64],

pynamicalsys/discrete_time/transport.py

@@ -390,7 +390,7 @@ def mean_squared_displacement(
     return output
 
 
-@njit(cache=True)
+@njit
 def recurrence_times(
     u: NDArray[np.float64],
     parameters: NDArray[np.float64],

{pynamicalsys-1.3.0.dist-info → pynamicalsys-1.3.1.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: pynamicalsys
-Version: 1.3.0
+Version: 1.3.1
 Summary: A Python toolkit for the analysis of dynamical systems
 Author-email: Matheus Rolim Sales <rolim.sales.m@gmail.com>
 License:                     GNU GENERAL PUBLIC LICENSE

{pynamicalsys-1.3.0.dist-info → pynamicalsys-1.3.1.dist-info}/RECORD

@@ -1,14 +1,14 @@
 pynamicalsys/__init__.py,sha256=IBHDcCD7wX8pULFXLk0yvF3uuNodvDKGirJektm1nR4,1215
-pynamicalsys/__version__.py,sha256=0Oc4EBzGTJOvXX0Vym4evglW1NQPpe8RLn8TdxsKzfs,704
+pynamicalsys/__version__.py,sha256=-4u7pjQ4caDQqa-1Qgms81j5hpkXjmjUYRCVEaLmb88,704
 pynamicalsys/common/__init__.py,sha256=W4OESm7TA-EUTGUOtIgBgvSd1n80Mgl0KA2dCkw8yIw,721
 pynamicalsys/common/basin_analysis.py,sha256=teGEOMPPOdzKNZsN0_OvkzBlcV-UOANRjcdFyI2-FLM,5473
-pynamicalsys/common/recurrence_quantification_analysis.py,sha256=neBolgoWrGJYOQta71Jf2ZMwNK1RoaZb_wNHqb1mM7I,14356
-pynamicalsys/common/utils.py,sha256=-sPpnoC3G8fJmdGZHzyQvbYU3Ym3BUMBlklaTgcz_Ac,10304
+pynamicalsys/common/recurrence_quantification_analysis.py,sha256=WQDFwM1aCJJuaATTVNsHzCXode53lzG_PqPoHmBcprs,14344
+pynamicalsys/common/utils.py,sha256=KexTurkehh_O7TR38iiBEnUVAAvqGRpwbxVKsqLDA4A,10268
 pynamicalsys/continuous_time/__init__.py,sha256=W4OESm7TA-EUTGUOtIgBgvSd1n80Mgl0KA2dCkw8yIw,721
-pynamicalsys/continuous_time/chaotic_indicators.py,sha256=UNZL1yJvUAvFC-RD-Xr3xZIEFCXgdpgw8KzmTh7U9OU,12542
+pynamicalsys/continuous_time/chaotic_indicators.py,sha256=kwoI259MoGntEBO7-LBnb-SG6JuTq0Ixgz8lxXNl-tk,12506
 pynamicalsys/continuous_time/models.py,sha256=YJtmsVZqSn3GDNqV0IKSQBHWZSCDY_EJ5Sgh4ulDc9I,5339
-pynamicalsys/continuous_time/numerical_integrators.py,sha256=CtP8GWNWGE5_2FYnrpbKFZMKklJdZBDWq3FkeXoIG2U,9684
-pynamicalsys/continuous_time/trajectory_analysis.py,sha256=hpoEw0EZLRV3UoWEGM7P7rH6LLY8lRNZXJ2ymQnNNDw,5256
+pynamicalsys/continuous_time/numerical_integrators.py,sha256=loIzJenp04ZWiukCK3oCs8vn9BeEcxckDO58B7QRKTo,9600
+pynamicalsys/continuous_time/trajectory_analysis.py,sha256=GwAMFuvH0kLqCaROXSFoXmj1PYauwwB3h5M_gR17FLY,5220
 pynamicalsys/continuous_time/validators.py,sha256=dvtlcWbn2w5_vNUB68f7ywFPyFrPKe3cUDqGCk52mhw,3834
 pynamicalsys/core/__init__.py,sha256=W4OESm7TA-EUTGUOtIgBgvSd1n80Mgl0KA2dCkw8yIw,721
 pynamicalsys/core/basin_metrics.py,sha256=sqSj-UCzBQdRaTvOBtzo8SCgJN90T3GiWy4nlLsbg5c,7604
@@ -17,12 +17,12 @@ pynamicalsys/core/discrete_dynamical_systems.py,sha256=JcDL71M3kaQJLnkWOVTG28ylo
 pynamicalsys/core/plot_styler.py,sha256=sVX_rb1HSt340U3zyPlXojylhma9cKCOgYirMIge18M,5647
 pynamicalsys/core/time_series_metrics.py,sha256=7fwXO1I6lIjmv_MzYxBn-9k2gjjnhm_y_Tp3diE1z68,5102
 pynamicalsys/discrete_time/__init__.py,sha256=W4OESm7TA-EUTGUOtIgBgvSd1n80Mgl0KA2dCkw8yIw,721
-pynamicalsys/discrete_time/dynamical_indicators.py,sha256=XDZ6e-nFCDO7UgkkBnuqsr70c3LLJ8bvew82pXrV1no,51759
-pynamicalsys/discrete_time/models.py,sha256=u9MlaCzBzR_OomWRYbEC2sS2qpfDg7qzZjX8YDijVe0,10756
-pynamicalsys/discrete_time/trajectory_analysis.py,sha256=W3nqLJgqJYMV9DbBbSSOpzNDQ7XlroP3lC0sgSoFztY,48669
-pynamicalsys/discrete_time/transport.py,sha256=mp6jqI8zfqwVrV5ZXqyOaFKzAjuCPXl6pQtzmX-gDok,14429
+pynamicalsys/discrete_time/dynamical_indicators.py,sha256=2Dz-ubootvvey2Hx9g3bGDSVUty3kVSXWvosAlwysHU,51651
+pynamicalsys/discrete_time/models.py,sha256=8mw7jrjMLDxKF7qtRyOAIf6zRu7koX43epdzSbqwi1o,10732
+pynamicalsys/discrete_time/trajectory_analysis.py,sha256=z9obMDO62pvUQI1xVikAVo5m6BiyMu7Rmvt3uhKiuYg,48549
+pynamicalsys/discrete_time/transport.py,sha256=CyRiF5ZNLQVUgHbFANH9wIGJuMTOyPRg9Q67_3C8zFM,14417
 pynamicalsys/discrete_time/validators.py,sha256=sk5PuQ6vxNj6ReqIfBTb2iQqlgY0SUkndS56xcp3w_w,10013
-pynamicalsys-1.3.0.dist-info/METADATA,sha256=_A4OpfJ39pd-eGH-HY08lyuq-KcDbDZB6ktdNb-uWwk,45897
-pynamicalsys-1.3.0.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
-pynamicalsys-1.3.0.dist-info/top_level.txt,sha256=1AqoHXye_hGRkvqu4KY0SdEuNTgqf3xyTlV8xhiJJFU,13
-pynamicalsys-1.3.0.dist-info/RECORD,,
+pynamicalsys-1.3.1.dist-info/METADATA,sha256=yTQx9DYbgYQygq4dfqDmUmdIETMXrha58TswTb8_A-w,45897
+pynamicalsys-1.3.1.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
+pynamicalsys-1.3.1.dist-info/top_level.txt,sha256=1AqoHXye_hGRkvqu4KY0SdEuNTgqf3xyTlV8xhiJJFU,13
+pynamicalsys-1.3.1.dist-info/RECORD,,