rustima 0.1.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.1.0 → rustima-0.3.0}/Cargo.lock

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

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

@@ -1,6 +1,6 @@
 [package]
 name = "rustima"
-version = "0.1.0"
+version = "0.3.0"
 edition = "2021"
 rust-version = "1.83"
 description = "High-performance SARIMAX engine with PyO3 bindings"
@@ -27,7 +27,9 @@ statrs = "0.18"
 rayon = "1.10"
 
 # Python bindings
-pyo3 = { version = "0.28", features = ["auto-initialize"] }
+# auto-initialize 제거(2026-07-11): 크레이트 내 Python::with_gil 사용처 0 — 임베딩 불필요.
+# 켜두면 manylinux(정적 python) 휠 빌드가 실패해 PyPI 릴리스를 막는다.
+pyo3 = { version = "0.28" }
 numpy = "0.28"
 
 # Error handling

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

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: rustima
-Version: 0.1.0
+Version: 0.3.0
 Classifier: Development Status :: 4 - Beta
 Classifier: Intended Audience :: Science/Research
 Classifier: Intended Audience :: Developers
@@ -13,6 +13,7 @@ Classifier: Programming Language :: Python :: 3.10
 Classifier: Programming Language :: Python :: 3.11
 Classifier: Programming Language :: Python :: 3.12
 Classifier: Programming Language :: Python :: 3.13
+Classifier: Programming Language :: Python :: 3.14
 Classifier: Programming Language :: Rust
 Classifier: Topic :: Scientific/Engineering
 Classifier: Topic :: Scientific/Engineering :: Mathematics

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

@@ -1,6 +1,6 @@
 use criterion::{criterion_group, criterion_main, Criterion};
-use sarimax_rs::optimizer;
-use sarimax_rs::types::{SarimaxConfig, SarimaxOrder, Trend};
+use rustima::optimizer;
+use rustima::types::{SarimaxConfig, SarimaxOrder, Trend};
 
 /// Deterministic LCG data generator (fixed seed → reproducible data).
 fn generate_ar1_data(n: usize, phi: f64, seed: u64) -> Vec<f64> {
@@ -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.1.0 → rustima-0.3.0}/benches/bench_kalman.rs

@@ -1,9 +1,9 @@
 use criterion::{criterion_group, criterion_main, Criterion};
-use sarimax_rs::initialization::KalmanInit;
-use sarimax_rs::kalman::kalman_loglike;
-use sarimax_rs::params::SarimaxParams;
-use sarimax_rs::state_space::StateSpace;
-use sarimax_rs::types::{SarimaxConfig, SarimaxOrder, Trend};
+use rustima::initialization::KalmanInit;
+use rustima::kalman::kalman_loglike;
+use rustima::params::SarimaxParams;
+use rustima::state_space::StateSpace;
+use rustima::types::{SarimaxConfig, SarimaxOrder, Trend};
 
 /// Deterministic LCG data generator (fixed seed → reproducible data).
 fn generate_ar1_data(n: usize, phi: f64, seed: u64) -> Vec<f64> {
@@ -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.1.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.3.0/lbfgsb_c/License.txt

@@ -0,0 +1,60 @@
+L-BFGS-B License
+================
+
+This directory contains a C translation of L-BFGS-B (version 3.0):
+
+  - Original FORTRAN L-BFGS-B: Ciyou Zhu, Richard Byrd, Jorge Nocedal,
+    Jose Luis Morales (Northwestern University).
+  - C translation (L-BFGS-B-C): Stephen Becker, 2015
+    (stephen.becker@colorado.edu), https://github.com/stephenbeckr/L-BFGS-B-C
+
+L-BFGS-B is released under the "New BSD License" (aka "Modified BSD License"
+or "3-clause license").
+
+New BSD License
+---------------
+
+Redistribution and use in source and binary forms, with or without
+modification, are permitted provided that the following conditions are met:
+
+1. Redistributions of source code must retain the above copyright notice,
+   this list of conditions and the following disclaimer.
+
+2. Redistributions in binary form must reproduce the above copyright notice,
+   this list of conditions and the following disclaimer in the documentation
+   and/or other materials provided with the distribution.
+
+3. Neither the name of the copyright holder nor the names of its contributors
+   may be used to endorse or promote products derived from this software
+   without specific prior written permission.
+
+THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
+AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
+IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
+ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE
+LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
+CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
+SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
+INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
+CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
+ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
+POSSIBILITY OF SUCH DAMAGE.
+
+Citation request (from the original authors)
+--------------------------------------------
+
+The authors of L-BFGS-B expect that all publications describing work using
+this software, or all commercial products using it, quote at least one of
+the references below:
+
+  R. H. Byrd, P. Lu and J. Nocedal. A Limited Memory Algorithm for Bound
+  Constrained Optimization (1995), SIAM Journal on Scientific and
+  Statistical Computing, 16, 5, pp. 1190-1208.
+
+  C. Zhu, R. H. Byrd and J. Nocedal. Algorithm 778: L-BFGS-B, FORTRAN
+  routines for large scale bound constrained optimization (1997), ACM
+  Transactions on Mathematical Software, 23, 4, pp. 550-560.
+
+  J.L. Morales and J. Nocedal. Remark on "Algorithm 778: L-BFGS-B, FORTRAN
+  routines for large scale bound constrained optimization" (2011), ACM
+  Transactions on Mathematical Software, 38, 1.

{rustima-0.1.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.1.0 → rustima-0.3.0}/lbfgsb_c/lbfgsb.h

@@ -5,6 +5,22 @@
 #ifndef lbfgsb_h
 #define lbfgsb_h
 
+/* 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
+
 
 
  /* Force 64-bit integers on all platforms (including Windows LLP64

{rustima-0.1.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.1.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.1.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.1.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;