rustima 0.2.0__tar.gz → 0.3.0__tar.gz

rustima-0.3.0/CHANGELOG.md

@@ -0,0 +1,112 @@
+# Changelog
+
+All notable changes to rustima are documented here.
+Versions follow `0.MINOR.PATCH`, where breaking changes bump MINOR.
+
+## 0.3.0 — 2026-07-25
+
+A diagnosis sweep of the engine (~1,800 adversarial inputs, 9 numerical
+regimes against statsmodels, and a structural audit) found and fixed four
+defects that returned **wrong numbers without raising**. Full report:
+`planner/ver9/DIAGNOSIS_V9.md`.
+
+> **Results from 0.2.x are not comparable to 0.3.0.** The log-likelihood,
+> AIC and BIC of every model **without differencing** (`d=0, D=0`) were
+> biased, and `inference="hessian"` standard errors were inflated. If you
+> have published or stored numbers from 0.2.x, re-run them.
+
+### Breaking — corrected numerical results
+
+- **Initial state covariance for non-differenced models.** `P₀` for the
+  stationary ARMA block used the DARE fixed point (the steady-state
+  *filtering* covariance) instead of the unconditional *Lyapunov*
+  covariance. For AR(1) this used `σ²` where the correct prior is
+  `σ²/(1-φ²)`. The error grows like `1/(1-φ²)` and does not vanish with
+  sample size: `Δloglike` ≈ −0.1 at φ=0.5, −62 at φ=0.995, and `ΔAIC`
+  +156 for a seasonal SAR=0.99 model. Fixed-parameter log-likelihoods now
+  match statsmodels to 1e-10. Models with `d≥1` were always correct.
+- **Standard errors from `inference="hessian"`.** The Hessian was mapped
+  from unconstrained to constrained space with the Jacobian instead of its
+  inverse, applying the transform twice and inflating AR/MA standard
+  errors by up to 6.6×. Now matches statsmodels `cov_params_approx` to
+  five decimals. `inference="opg"` was unaffected.
+- **False convergence with trend terms.** A deterministic trend made the
+  CSS/Burg start-parameter estimator see near-integration, saturating the
+  AR start on the constraint boundary where L-BFGS-B stalls and reports
+  `converged=True` with a gradient of order 100 — up to 192 nats below the
+  optimum. Start parameters are now detrended, and out-of-box estimates
+  fall back to zeros instead of clamping.
+- **Near-cancellation detection.** AR and MA roots were compared against
+  the sign-flipped MA polynomial, so genuine cancellation (`θ = −φ`) was
+  never detected — the warning and the multi-start restart filter both ran
+  on the wrong condition.
+- **Post-estimation ran a differently-initialized filter.** `forecast`,
+  `rolling_forecast`, `residuals`, `batch_forecast` and `diagnostics`
+  ignored the model's enforcement flags and silently used an approximate
+  diffuse initialization, so `res.llf` and `res.resid` / `res.diagnostics()`
+  described different models — Ljung-Box tested residuals the likelihood
+  never saw.
+
+### Breaking — stricter inputs and API
+
+- Masked arrays are rejected by `SARIMAXModel` and `auto_arima` (previously
+  the mask was silently dropped and masked entries used as data).
+- `sarimax_inference` validates the parameter-vector length instead of
+  returning empty arrays.
+- `SARIMAXResult.param_names` raises if the Rust layout and the Python
+  naming disagree, instead of padding with `param_N` and printing values
+  under the wrong names.
+- `sarimax_rolling_forecast` refuses configurations whose origin snapshots
+  would exceed 2 GiB (previously unbounded — 4.4 GB measured at n=1000,
+  s=365, and an allocation failure aborts the process).
+- `sarimax_residuals` returns an additional `prediction_variances` key.
+  Code asserting an exact key set must be updated.
+- Rust: the never-supported `SarimaxConfig.measurement_error` field is
+  removed.
+
+### Added
+
+- `PredictionResult.conf_int(alpha=...)` and `.se_mean` —
+  `get_prediction()` now produces intervals across the whole range, not
+  only the out-of-sample tail.
+- `sarimax_residuals` returns `prediction_variances` (one-step-ahead
+  prediction variance).
+- `sarimax_forecast`, `sarimax_rolling_forecast`, `sarimax_residuals`,
+  `sarimax_batch_forecast` and `sarimax_diagnostics` accept
+  `enforce_stationarity` / `enforce_invertibility` (appended to the
+  signature; positional calls are unaffected).
+
+### Fixed
+
+- **Deadlock in forked processes.** Rayon's global pool is not fork-safe:
+  any batch or grid call in a `fork()`ed child hung forever with no error
+  (reachable through `multiprocessing`'s fork context after a single
+  earlier batch call). Parallel entry points now fail fast with a clear
+  message, and the optimizer's internal parallel sections fall back to
+  sequential execution.
+- Batch workers isolate panics per series, so the documented per-series
+  error isolation holds even if a worker panics.
+- The analytical gradient no longer disappears near unit roots: the score
+  pass now shares the likelihood pass's `F_t ≤ 0` fallback and
+  steady-state adoption guard, keeping the gradient finite wherever the
+  log-likelihood is.
+- The analytical score accounts for `∂P₀/∂θ` when the initialization
+  depends on the parameters.
+
+### Internal
+
+- Single-sourced the trend basis, the flat-parameter count, and
+  single-series input validation (the last was duplicated into
+  `sarimax_grid_search`, the `auto_arima` hot path). Model flags moved from
+  adjacent positional booleans into a named struct.
+- Removed dead code: `auto_arima`'s unused evaluation path, the
+  `measurement_error` plumbing, and the redundant Python numerical-Hessian
+  implementation.
+- Python test suite 407 → 429, including regression tests for every defect
+  above and eight previously uncovered option-matrix combinations. CI is
+  green for the first time since 2026-05-23.
+
+## 0.2.0
+
+Initial tagged release with the SARIMAX engine, `auto_arima`, batch
+processing and the statsmodels-compatible Python API.

{rustima-0.2.0 → rustima-0.3.0}/Cargo.lock

@@ -880,7 +880,7 @@ checksum = "357703d41365b4b27c590e3ed91eabb1b663f07c4c084095e60cbed4362dff0d"
 
 [[package]]
 name = "rustima"
-version = "0.2.0"
+version = "0.3.0"
 dependencies = [
  "anyhow",
  "approx",

{rustima-0.2.0 → rustima-0.3.0}/Cargo.toml

@@ -1,6 +1,6 @@
 [package]
 name = "rustima"
-version = "0.2.0"
+version = "0.3.0"
 edition = "2021"
 rust-version = "1.83"
 description = "High-performance SARIMAX engine with PyO3 bindings"

{rustima-0.2.0 → rustima-0.3.0}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: rustima
-Version: 0.2.0
+Version: 0.3.0
 Classifier: Development Status :: 4 - Beta
 Classifier: Intended Audience :: Science/Research
 Classifier: Intended Audience :: Developers

{rustima-0.2.0 → rustima-0.3.0}/benches/bench_fit.rs

@@ -26,7 +26,6 @@ fn bench_fit_ar1(c: &mut Criterion) {
         enforce_invertibility: true,
         concentrate_scale: true,
         simple_differencing: false,
-        measurement_error: false,
     };
 
     c.bench_function("fit_ar1_n500", |b| {
@@ -44,7 +43,6 @@ fn bench_fit_arima111(c: &mut Criterion) {
         enforce_invertibility: true,
         concentrate_scale: true,
         simple_differencing: false,
-        measurement_error: false,
     };
 
     c.bench_function("fit_arima111_n500", |b| {
@@ -66,7 +64,6 @@ fn bench_fit_sarima_111_111_12(c: &mut Criterion) {
         enforce_invertibility: false,
         concentrate_scale: true,
         simple_differencing: false,
-        measurement_error: false,
     };
 
     c.bench_function("fit_sarima111_111_12_n500", |b| {

{rustima-0.2.0 → rustima-0.3.0}/benches/bench_kalman.rs

@@ -30,7 +30,6 @@ fn bench_kalman_ar1(c: &mut Criterion) {
         enforce_invertibility: false,
         concentrate_scale: true,
         simple_differencing: false,
-        measurement_error: false,
     };
     let params = SarimaxParams {
         trend_coeffs: vec![],
@@ -60,7 +59,6 @@ fn bench_kalman_arima111(c: &mut Criterion) {
         enforce_invertibility: false,
         concentrate_scale: true,
         simple_differencing: false,
-        measurement_error: false,
     };
     let params = SarimaxParams {
         trend_coeffs: vec![],

{rustima-0.2.0 → rustima-0.3.0}/examples/high_order_bench.rs

@@ -28,7 +28,6 @@ fn run_test(label: &str, order: SarimaxOrder, n: usize) {
         enforce_invertibility: true,
         concentrate_scale: true,
         simple_differencing: false,
-        measurement_error: false,
     };
 
     let start = Instant::now();

{rustima-0.2.0 → rustima-0.3.0}/lbfgsb_c/lbfgsb.c

@@ -60,7 +60,7 @@ static integer c__1 = 1;
 
 
     /* Local variables */
-    _Thread_local static integer ld, lr, lt, lz, lwa, lwn, lss, lxp, lws, lwt, lsy, lwy,
+    static _Thread_local integer ld, lr, lt, lz, lwa, lwn, lss, lxp, lws, lwt, lsy, lwy,
 	    lsnd;
 
 /* -jlm-jn */
@@ -343,45 +343,45 @@ static double c_b7 = 0.;
     fileType o__1=NULL;
 
     /* Local variables */
-    _Thread_local static integer i__, k;
-    _Thread_local static double gd, dr, rr, dtd;
-    _Thread_local static integer col;
-    _Thread_local static double tol;
-    _Thread_local static logical wrk;
-    _Thread_local static double stp, cpu1, cpu2;
-    _Thread_local static integer head;
-    _Thread_local static double fold;
-    _Thread_local static integer nact;
-    _Thread_local static double ddum;
-    _Thread_local static integer info, nseg;
-    _Thread_local static double time;
-    _Thread_local static integer nfgv, ifun, iter;
-    _Thread_local static integer wordTemp;
-    _Thread_local static integer *word;
+    static _Thread_local integer i__, k;
+    static _Thread_local double gd, dr, rr, dtd;
+    static _Thread_local integer col;
+    static _Thread_local double tol;
+    static _Thread_local logical wrk;
+    static _Thread_local double stp, cpu1, cpu2;
+    static _Thread_local integer head;
+    static _Thread_local double fold;
+    static _Thread_local integer nact;
+    static _Thread_local double ddum;
+    static _Thread_local integer info, nseg;
+    static _Thread_local double time;
+    static _Thread_local integer nfgv, ifun, iter;
+    static _Thread_local integer wordTemp;
+    static _Thread_local integer *word;
     word = &wordTemp;
-    _Thread_local static double time1, time2;
-    _Thread_local static integer iback;
-    _Thread_local static double gdold;
-    _Thread_local static integer nfree;
-    _Thread_local static logical boxed;
-    _Thread_local static integer itail;
-    _Thread_local static double theta;
-    _Thread_local static double dnorm;
-    _Thread_local static integer nskip, iword;
-    _Thread_local static double xstep, stpmx;
-    _Thread_local static integer ileave;
-    _Thread_local static double cachyt;
-    _Thread_local static integer itfile;
-    _Thread_local static double epsmch;
-    _Thread_local static logical updatd;
-    _Thread_local static double sbtime;
-    _Thread_local static logical prjctd;
-    _Thread_local static integer iupdat;
-    _Thread_local static double sbgnrm;
-    _Thread_local static logical cnstnd;
-    _Thread_local static integer nenter;
-    _Thread_local static double lnscht;
-    _Thread_local static integer nintol;
+    static _Thread_local double time1, time2;
+    static _Thread_local integer iback;
+    static _Thread_local double gdold;
+    static _Thread_local integer nfree;
+    static _Thread_local logical boxed;
+    static _Thread_local integer itail;
+    static _Thread_local double theta;
+    static _Thread_local double dnorm;
+    static _Thread_local integer nskip, iword;
+    static _Thread_local double xstep, stpmx;
+    static _Thread_local integer ileave;
+    static _Thread_local double cachyt;
+    static _Thread_local integer itfile;
+    static _Thread_local double epsmch;
+    static _Thread_local logical updatd;
+    static _Thread_local double sbtime;
+    static _Thread_local logical prjctd;
+    static _Thread_local integer iupdat;
+    static _Thread_local double sbgnrm;
+    static _Thread_local logical cnstnd;
+    static _Thread_local integer nenter;
+    static _Thread_local double lnscht;
+    static _Thread_local integer nintol;
 
 /* -jlm-jn */
 /*     ************ */

{rustima-0.2.0 → rustima-0.3.0}/lbfgsb_c/lbfgsb.h

@@ -5,10 +5,20 @@
 #ifndef lbfgsb_h
 #define lbfgsb_h
 
-/* MSVC's C mode has no C11 _Thread_local keyword — map to its equivalent
- * so Windows wheel builds compile. gcc/clang paths are unaffected. */
+/* Portable thread-local shim. The C11 _Thread_local keyword is missing on:
+ *  - MSVC C mode                          -> __declspec(thread)
+ *  - gcc < 4.9 (manylinux2014 aarch64
+ *    cross toolchain is gcc 4.8; it even
+ *    accepts -std=c11 without the keyword) -> GNU __thread
+ *  - any pre-C11 mode                      -> GNU __thread
+ * Modern clang/gcc with C11 keep the native keyword. */
 #if defined(_MSC_VER) && !defined(__clang__)
 #define _Thread_local __declspec(thread)
+#elif defined(__GNUC__) && !defined(__clang__) && \
+    (__GNUC__ < 4 || (__GNUC__ == 4 && __GNUC_MINOR__ < 9))
+#define _Thread_local __thread
+#elif !defined(__STDC_VERSION__) || (__STDC_VERSION__ < 201112L)
+#define _Thread_local __thread
 #endif
 
 

{rustima-0.2.0 → rustima-0.3.0}/lbfgsb_c/linesearch.c

@@ -43,18 +43,18 @@ static integer c__1 = 1;
 
 
     /* Table of constant values */
-    _Thread_local static double c_b14 = FTOL;
-    _Thread_local static double c_b15 = GTOL;
-    _Thread_local static double c_b16 = XTOL;
-    _Thread_local static double c_b17 = STEPMIN;
+    static _Thread_local double c_b14 = FTOL;
+    static _Thread_local double c_b15 = GTOL;
+    static _Thread_local double c_b16 = XTOL;
+    static _Thread_local double c_b17 = STEPMIN;
     /* System generated locals */
     integer i__1;
     double d__1;
 
 
     /* Local variables */
-    _Thread_local static integer i__;
-    _Thread_local static double a1, a2;
+    static _Thread_local integer i__;
+    static _Thread_local double a1, a2;
 
     /* Parameter adjustments */
     --z__;
@@ -170,10 +170,10 @@ int dcsrch(double *f, double *g, double *stp,
 
 
     /* Local variables */
-    _Thread_local static double fm, gm, fx, fy, gx, gy, fxm, fym, gxm, gym, stx, sty;
-    _Thread_local static integer stage;
-    _Thread_local static double finit, ginit, width, ftest, gtest, stmin, stmax, width1;
-    _Thread_local static logical brackt;
+    static _Thread_local double fm, gm, fx, fy, gx, gy, fxm, fym, gxm, gym, stx, sty;
+    static _Thread_local integer stage;
+    static _Thread_local double finit, ginit, width, ftest, gtest, stmin, stmax, width1;
+    static _Thread_local logical brackt;
 
     /*
      ********** 
@@ -496,7 +496,7 @@ L1000:
     double sqrt(double);
 
     /* Local variables */
-    _Thread_local static double p, q, r__, s, sgnd, stpc, stpf, stpq, gamma, theta;
+    static _Thread_local double p, q, r__, s, sgnd, stpc, stpf, stpq, gamma, theta;
 
     /*
      ********** 

{rustima-0.2.0 → rustima-0.3.0}/lbfgsb_c/linpack.c

@@ -18,9 +18,9 @@ int dpofa(double *a, integer *lda, integer *n, integer *
     double sqrt(double);
 
     /* Local variables */
-    _Thread_local static integer j, k;
-    _Thread_local static double s, t;
-    _Thread_local static integer jm1;
+    static _Thread_local integer j, k;
+    static _Thread_local double s, t;
+    static _Thread_local integer jm1;
 
 /*
     dpofa factors a double precision symmetric positive definite 
@@ -115,8 +115,8 @@ int dtrsl(double *t, integer *ldt, integer *n,
     integer t_dim1, t_offset, i__1, i__2;
 
     /* Local variables */
-    _Thread_local static integer j, jj, case__;
-    _Thread_local static double temp;
+    static _Thread_local integer j, jj, case__;
+    static _Thread_local double temp;
     /*
     extern double ddot(integer *, double *, integer *, double *, 
 	    integer *);

{rustima-0.2.0 → rustima-0.3.0}/lbfgsb_c/miniCBLAS.c

@@ -52,7 +52,7 @@
     integer i__1;
 
     /* Local variables */
-    _Thread_local static integer i, m, ix, iy, mp1;
+    static _Thread_local integer i, m, ix, iy, mp1;
 
 
 /*     constant times a vector plus a vector.   
@@ -139,7 +139,7 @@ L40:
     integer i__1;
 
     /* Local variables */
-    _Thread_local static integer i, m, ix, iy, mp1;
+    static _Thread_local integer i, m, ix, iy, mp1;
 
 
     /*     copies a vector, x, to a vector, y.   
@@ -226,9 +226,9 @@ double ddotRef(integer *n, double *dx, integer *incx, double *dy,
     double ret_val;
 
     /* Local variables */
-    _Thread_local static integer i, m;
-    _Thread_local static double dtemp;
-    _Thread_local static integer ix, iy, mp1;
+    static _Thread_local integer i, m;
+    static _Thread_local double dtemp;
+    static _Thread_local integer ix, iy, mp1;
 
 
     /*     forms the dot product of two vectors.   
@@ -314,7 +314,7 @@ L60:
     integer i__1, i__2;
 
     /* Local variables */
-    _Thread_local static integer i, m, nincx, mp1;
+    static _Thread_local integer i, m, nincx, mp1;
 
 
     /*     scales a vector by a constant.   

{rustima-0.2.0 → rustima-0.3.0}/lbfgsb_c/print.c

@@ -30,7 +30,7 @@ int prn1lb(integer *n, integer *m, double *l,
     integer i__1;
 
     /* Local variables */
-    _Thread_local static integer i__;
+    static _Thread_local integer i__;
 
     /* Parameter adjustments */
     --x;
@@ -103,7 +103,7 @@ int prn1lb(integer *n, integer *m, double *l,
     integer i__1;
 
     /* Local variables */
-    _Thread_local static integer i__, imod;
+    static _Thread_local integer i__, imod;
 
     /* Parameter adjustments */
     --g;
@@ -186,7 +186,7 @@ int prn3lb(integer *n, double *x, double *f, integer *
     integer i__1;
 
     /* Local variables */
-    _Thread_local static integer i__;
+    static _Thread_local integer i__;
     /* Parameter adjustments */
     --x;
 
@@ -304,7 +304,7 @@ int errclb(integer *n, integer *m, double *factr,
     integer i__1;
 
     /* Local variables */
-    _Thread_local static integer i__;
+    static _Thread_local integer i__;
 
     /* Parameter adjustments */
     --nbd;

{rustima-0.2.0 → rustima-0.3.0}/lbfgsb_c/subalgorithms.c

@@ -12,7 +12,7 @@ static integer c__11 = 11;
     integer i__1;
 
     /* Local variables */
-    _Thread_local static integer i__, nbdd;
+    static _Thread_local integer i__, nbdd;
 
 
 /*     ************ */
@@ -127,8 +127,8 @@ static integer c__11 = 11;
     double sqrt(double);
 
     /* Local variables */
-    _Thread_local static integer i__, k, i2;
-    _Thread_local static double sum;
+    static _Thread_local integer i__, k, i2;
+    static _Thread_local double sum;
 
 /*     ************ */
 
@@ -275,28 +275,28 @@ static integer c__11 = 11;
 
 
     /* Local variables */
-    _Thread_local static integer i__, j;
-    _Thread_local static double f1, f2, dt, tj, tl, tu, tj0;
-    _Thread_local static integer ibp;
-    _Thread_local static double dtm;
+    static _Thread_local integer i__, j;
+    static _Thread_local double f1, f2, dt, tj, tl, tu, tj0;
+    static _Thread_local integer ibp;
+    static _Thread_local double dtm;
     extern /* Subroutine */ int bmv(integer *, double *, double *,
 	    integer *, double *, double *, integer *);
-    _Thread_local static double wmc, wmp, wmw;
-    _Thread_local static integer col2;
-    _Thread_local static double dibp;
-    _Thread_local static integer iter;
-    _Thread_local static double zibp, tsum, dibp2;
-    _Thread_local static logical bnded;
-    _Thread_local static double neggi;
-    _Thread_local static integer nfree;
-    _Thread_local static double bkmin;
-    _Thread_local static integer nleft;
-    _Thread_local static double f2_org__;
-    _Thread_local static integer nbreak, ibkmin;
+    static _Thread_local double wmc, wmp, wmw;
+    static _Thread_local integer col2;
+    static _Thread_local double dibp;
+    static _Thread_local integer iter;
+    static _Thread_local double zibp, tsum, dibp2;
+    static _Thread_local logical bnded;
+    static _Thread_local double neggi;
+    static _Thread_local integer nfree;
+    static _Thread_local double bkmin;
+    static _Thread_local integer nleft;
+    static _Thread_local double f2_org__;
+    static _Thread_local integer nbreak, ibkmin;
     extern /* Subroutine */ int hpsolb(integer *, double *, integer *,
 	    integer *);
-    _Thread_local static integer pointr;
-    _Thread_local static logical xlower, xupper;
+    static _Thread_local integer pointr;
+    static _Thread_local logical xlower, xupper;
 
 /*     ************ */
 
@@ -828,11 +828,11 @@ L999:
 	    wt_dim1, wt_offset, i__1, i__2;
 
     /* Local variables */
-    _Thread_local static integer i__, j, k;
-    _Thread_local static double a1, a2;
+    static _Thread_local integer i__, j, k;
+    static _Thread_local double a1, a2;
     extern /* Subroutine */ int bmv(integer *, double *, double *,
 	    integer *, double *, double *, integer *);
-    _Thread_local static integer pointr;
+    static _Thread_local integer pointr;
 
 /*     ************ */
 
@@ -928,10 +928,10 @@ L999:
 	    wy_dim1, wy_offset, sy_dim1, sy_offset, i__1, i__2, i__3;
 
     /* Local variables */
-    _Thread_local static integer i__, k, k1, m2, is, js, iy, jy, is1, js1, col2, dend, pend;
-    _Thread_local static integer upcl;
-    _Thread_local static double temp1, temp2, temp3, temp4;
-    _Thread_local static integer ipntr, jpntr, dbegin, pbegin;
+    static _Thread_local integer i__, k, k1, m2, is, js, iy, jy, is1, js1, col2, dend, pend;
+    static _Thread_local integer upcl;
+    static _Thread_local double temp1, temp2, temp3, temp4;
+    static _Thread_local integer ipntr, jpntr, dbegin, pbegin;
 
 /*     ************ */
 
@@ -1310,8 +1310,8 @@ L999:
 	    i__2, i__3;
 
     /* Local variables */
-    _Thread_local static integer i__, j, k, k1;
-    _Thread_local static double ddum;
+    static _Thread_local integer i__, j, k, k1;
+    static _Thread_local double ddum;
 
 /*     ************ */
 
@@ -1399,7 +1399,7 @@ L999:
     integer i__1;
 
     /* Local variables */
-    _Thread_local static integer i__, k, iact;
+    static _Thread_local integer i__, k, iact;
 
 
 /*     ************ */
@@ -1517,9 +1517,9 @@ L999:
     integer i__1;
 
     /* Local variables */
-    _Thread_local static integer i__, j, k;
-    _Thread_local static double out, ddum;
-    _Thread_local static integer indxin, indxou;
+    static _Thread_local integer i__, j, k;
+    static _Thread_local double out, ddum;
+    static _Thread_local integer indxin, indxou;
 
 /*     ************ */
 
@@ -1639,8 +1639,8 @@ L30:
 	    ss_dim1, ss_offset, i__1, i__2;
 
     /* Local variables */
-    _Thread_local static integer j;
-    _Thread_local static integer pointr;
+    static _Thread_local integer j;
+    static _Thread_local integer pointr;
 
 /*     ************ */
 
@@ -1746,8 +1746,8 @@ L30:
     double d__1, d__2;
 
     /* Local variables */
-    _Thread_local static integer i__;
-    _Thread_local static double gi;
+    static _Thread_local integer i__;
+    static _Thread_local double gi;
 
 /*     ************ */
 
@@ -1827,13 +1827,13 @@ L30:
     double d__1, d__2;
 
     /* Local variables */
-    _Thread_local static integer i__, j, k, m2;
-    _Thread_local static double dk;
-    _Thread_local static integer js, jy;
-    _Thread_local static double xk;
-    _Thread_local static integer ibd, col2;
-    _Thread_local static double dd_p__, temp1, temp2, alpha;
-    _Thread_local static integer pointr;
+    static _Thread_local integer i__, j, k, m2;
+    static _Thread_local double dk;
+    static _Thread_local integer js, jy;
+    static _Thread_local double xk;
+    static _Thread_local integer ibd, col2;
+    static _Thread_local double dd_p__, temp1, temp2, alpha;
+    static _Thread_local integer pointr;
 
 /*     ********************************************************************** */
 

{rustima-0.2.0 → rustima-0.3.0}/pyproject.toml

@@ -4,7 +4,7 @@ build-backend = "maturin"
 
 [project]
 name = "rustima"
-version = "0.2.0"
+version = "0.3.0"
 requires-python = ">=3.10"
 description = "High-performance SARIMAX (Seasonal ARIMA with eXogenous regressors) engine in Rust with statsmodels-compatible Python API"
 readme = "README.md"

{rustima-0.2.0 → rustima-0.3.0}/python/rustima/__init__.py

@@ -14,15 +14,15 @@ import os as _os
 if "RAYON_NUM_THREADS" not in _os.environ:
     import sys as _sys
     if _sys.platform == "darwin":
+        import subprocess as _sp
         try:
-            import subprocess as _sp
             _n = _sp.check_output(
                 ["sysctl", "-n", "hw.perflevel0.physicalcpu"],
                 stderr=_sp.DEVNULL, timeout=2,
             ).decode().strip()
             if _n.isdigit() and int(_n) > 0:
                 _os.environ["RAYON_NUM_THREADS"] = _n
-        except Exception:
+        except (OSError, _sp.SubprocessError, UnicodeDecodeError):
             pass
 
 # Native extension module — maturin places the .so here automatically.