l1tv 0.1.0__tar.gz

l1tv-0.1.0/.github/workflows/ci.yml

@@ -0,0 +1,70 @@
+name: CI
+
+on:
+  push:
+    branches: [master]
+  pull_request:
+
+jobs:
+  rust:
+    name: Rust ${{ matrix.os }}
+    runs-on: ${{ matrix.os }}
+    strategy:
+      fail-fast: false
+      matrix:
+        os: [ubuntu-latest, macos-latest, windows-latest]
+    steps:
+      - uses: actions/checkout@v4
+      # Pin Python so pyo3-ffi's build script doesn't pick up the runner's
+      # default Python (now 3.14 on GHA's latest images), which is newer than
+      # pyo3 0.24.2's supported range.
+      - uses: actions/setup-python@v5
+        with:
+          python-version: "3.12"
+      - uses: dtolnay/rust-toolchain@stable
+        with:
+          components: clippy, rustfmt
+      - uses: Swatinem/rust-cache@v2
+      - name: cargo fmt
+        run: cargo fmt --all -- --check
+      # clippy is platform-agnostic in intent; running it on one platform is
+      # enough and avoids macOS/Windows-specific lint flakes.
+      - name: cargo clippy
+        if: matrix.os == 'ubuntu-latest'
+        run: cargo clippy --all-targets -- -D warnings
+      - name: cargo test
+        run: cargo test --release
+
+  python:
+    name: Python ${{ matrix.python-version }} ${{ matrix.os }}
+    runs-on: ${{ matrix.os }}
+    strategy:
+      fail-fast: false
+      matrix:
+        os: [ubuntu-latest, macos-latest, windows-latest]
+        python-version: ["3.9", "3.12"]
+    steps:
+      - uses: actions/checkout@v4
+      - uses: actions/setup-python@v5
+        with:
+          python-version: ${{ matrix.python-version }}
+      - uses: dtolnay/rust-toolchain@stable
+      - uses: Swatinem/rust-cache@v2
+      # PEP 517 install — pulls maturin from build-system.requires, builds the
+      # Rust extension via maturin, and installs the package + test extras
+      # (pytest, scipy, matplotlib). Avoids the venv-required `maturin develop`
+      # path that fails on stock GitHub-hosted runners.
+      - name: install l1tv (PEP 517)
+        run: pip install ".[test]"
+      # Run pytest and demos from a directory that does NOT contain `l1tv/`,
+      # otherwise Python's cwd-on-sys.path rule shadows the installed package
+      # with the source-tree copy (which lacks the compiled `_core` extension).
+      - name: pytest
+        working-directory: ${{ runner.temp }}
+        run: pytest -v ${{ github.workspace }}/tests
+      - name: smoke demos
+        working-directory: ${{ runner.temp }}
+        shell: bash
+        run: |
+          python "${GITHUB_WORKSPACE}/demos_python/demo_real.py"
+          python "${GITHUB_WORKSPACE}/demos_python/demo_circ.py"

l1tv-0.1.0/.github/workflows/release.yml

@@ -0,0 +1,111 @@
+name: Release
+
+on:
+  push:
+    tags:
+      - "v*"
+  workflow_dispatch:
+
+permissions:
+  contents: read
+  id-token: write  # for trusted publishing to PyPI
+
+jobs:
+  linux:
+    runs-on: ubuntu-latest
+    strategy:
+      fail-fast: false
+      matrix:
+        target: [x86_64, aarch64]
+    steps:
+      - uses: actions/checkout@v4
+      - uses: actions/setup-python@v5
+        with: { python-version: "3.12" }
+      - uses: PyO3/maturin-action@v1
+        with:
+          target: ${{ matrix.target }}
+          # Drop --zig and rely on `manylinux: auto`. With --zig, maturin
+          # cross-links on the host and needs an aarch64-targeted Python on
+          # PATH (which there isn't), causing the aarch64 build to fail with
+          # "Couldn't find any python interpreters". `manylinux: auto` runs
+          # the build inside the official quay.io/pypa/manylinux*_<arch>
+          # container, which ships /opt/python/cp3X-cp3X/bin/python for
+          # every supported version; --find-interpreter picks them up so we
+          # get a wheel per Python per arch.
+          args: --release --out dist --find-interpreter
+          manylinux: auto
+      - uses: actions/upload-artifact@v4
+        with:
+          name: wheels-linux-${{ matrix.target }}
+          path: dist
+
+  macos:
+    # macos-latest is Apple Silicon (arm64). Phase 4 cross-cutting decision:
+    # drop the macos-13 (x86) runner family-wide due to 24-hour queue hangs.
+    # Intel-mac users fall back to source builds.
+    runs-on: macos-latest
+    strategy:
+      fail-fast: false
+      matrix:
+        target: [x86_64, aarch64]
+    steps:
+      - uses: actions/checkout@v4
+      - uses: actions/setup-python@v5
+        with: { python-version: "3.12" }
+      - uses: PyO3/maturin-action@v1
+        with:
+          target: ${{ matrix.target }}
+          args: --release --out dist
+      - uses: actions/upload-artifact@v4
+        with:
+          name: wheels-macos-${{ matrix.target }}
+          path: dist
+
+  windows:
+    runs-on: windows-latest
+    strategy:
+      fail-fast: false
+      matrix:
+        target: [x64]
+    steps:
+      - uses: actions/checkout@v4
+      - uses: actions/setup-python@v5
+        with: { python-version: "3.12" }
+      - uses: PyO3/maturin-action@v1
+        with:
+          target: ${{ matrix.target }}
+          args: --release --out dist
+      - uses: actions/upload-artifact@v4
+        with:
+          name: wheels-windows-${{ matrix.target }}
+          path: dist
+
+  sdist:
+    runs-on: ubuntu-latest
+    steps:
+      - uses: actions/checkout@v4
+      - uses: PyO3/maturin-action@v1
+        with:
+          command: sdist
+          args: --out dist
+      - uses: actions/upload-artifact@v4
+        with:
+          name: wheels-sdist
+          path: dist
+
+  publish:
+    name: Publish to PyPI
+    runs-on: ubuntu-latest
+    needs: [linux, macos, windows, sdist]
+    if: startsWith(github.ref, 'refs/tags/v')
+    # Tied to the `pypi` GitHub Actions environment — must match the
+    # "Environment name" in PyPI's trusted-publisher config for `l1tv`.
+    # Required reviewers / branch rules configured at the env level act as
+    # the gate before this job mints an OIDC token and uploads.
+    environment: pypi
+    steps:
+      - uses: actions/download-artifact@v4
+        with:
+          path: dist
+          merge-multiple: true
+      - uses: pypa/gh-action-pypi-publish@release/v1

l1tv-0.1.0/.gitignore

@@ -0,0 +1,50 @@
+
+unitTests/swUnitTestDistTransformL1.m
+
+Docs/create_images.m
+
+unitTests/swDistTransformL1FH.m
+
+unitTests/swUnitTestDistAngles.m
+
+unitTests/swUnitTestDistTransformNE.m
+
+unitTests/swUnitTestMinL1TVCirc.m
+
+# Rust build artifacts
+target/
+Cargo.lock.bak
+
+# Compiled Python extension (platform-specific; built locally or by CI)
+l1tv/_core*.so
+l1tv/_core*.pyd
+
+# Python cache
+__pycache__/
+*.pyc
+*.pyo
+.pytest_cache/
+
+# macOS
+.DS_Store
+
+# Distribution / packaging
+dist/
+build/
+*.egg-info/
+.eggs/
+
+# Editors / IDEs
+.vscode/
+.idea/
+*.swp
+
+# Jupyter
+.ipynb_checkpoints/
+
+# Demo plot artifacts (regenerated on demand)
+demos_python/*.png
+
+# Virtualenvs
+.venv/
+venv/

l1tv-0.1.0/Auxiliary/deltaSNR.m

@@ -0,0 +1,14 @@
+function dsnr = deltaSNR( groundTruth, data, estimate, dataSpace)
+%deltaSNR Computes the signal to noise ratio improvement of a restoration
+
+switch dataSpace
+    case 'circ'
+        dist = @(x,y) distAngle(x,y);
+    case 'real'
+        dist = @(x,y) abs(x - y);
+end
+
+dsnr = 10 * log10(sum( dist(groundTruth, data).^2 ) / sum( dist(groundTruth, estimate).^2 ));
+
+end
+

l1tv-0.1.0/Auxiliary/distAngle.m

@@ -0,0 +1,14 @@
+function d = distAngle( phi, psi )
+%distAngle Distance between two angles phi and psi
+
+if (all(abs(phi) <= pi)) && (all(abs(psi) <= pi))
+    % cheap computation for angles in the interval [-pi, pi]
+    aux = phi - psi;
+    d = min(abs([aux(:) + 2 * pi, aux(:), aux(:) - 2 * pi]), [], 2);
+    d = reshape(d, size(aux));
+else
+    % more expensive calculation for angles outside the interval [-pi, pi]
+    aux = abs(angle(exp(1i * (phi - psi))));
+    d = min(aux, 2*pi-aux);
+end
+

l1tv-0.1.0/Auxiliary/distTransformL1.m

@@ -0,0 +1,15 @@
+function B = distTransformL1( B, R, alpha )
+%distTransformL1 Fast computation of L1 distance transform
+
+K = numel(B);
+% forward pass
+for k=2:K
+    B(k) = min( B(k), B(k-1) + alpha * (R(k) - R(k-1)));
+end
+% backward pass
+for k=K-1:-1:1
+    B(k) = min( B(k), B(k+1) + alpha * (R(k+1) - R(k)));
+end
+
+end
+

l1tv-0.1.0/Auxiliary/randCP.m

@@ -0,0 +1,13 @@
+function r = randCP( r, lambda)
+%randCP Bernoulli-thinned innovations approximating a compound Poisson.
+% Each entry of r is independently zeroed with probability exp(-lambda);
+% otherwise it is left at the supplied value. This matches the marginal
+% P(N = 0) of a Poisson(lambda) point process and is used in the demos as
+% a cheap surrogate for compound-Poisson innovations (good for small
+% lambda; not the actual compound Poisson distribution).
+
+x = rand(size(r)); % uniform distr. vector
+r(x(:) <= exp(-lambda)) = 0;
+
+end
+

l1tv-0.1.0/Auxiliary/randl.m

@@ -0,0 +1,10 @@
+function y = randl( varargin )
+%randl Random normal Laplacian noise, mean = 0, variance sigma^2 = 1
+% pdf: $\frac{1}{\sqrt{2}} e^{- \sqrt{2} |x|}$
+% Same input arguments as function rand
+
+x = rand( varargin{:} ) - 0.5;
+y = -sign(x) .* log(1 - 2 * abs(x)) / sqrt(2);
+
+end
+

l1tv-0.1.0/Auxiliary/wrapAngle.m

@@ -0,0 +1,7 @@
+function x = wrapAngle( y )
+%wrapAngle Wraps the angle y to the interval [-pi, pi]
+
+x = angle(exp(1i * y));
+
+end
+

l1tv-0.1.0/CHANGELOG.md

@@ -0,0 +1,59 @@
+# Changelog
+
+All notable changes to the `l1tv` Python package are documented in this file.
+Format follows [Keep a Changelog](https://keepachangelog.com/en/1.1.0/); this
+project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0.html).
+
+The MATLAB reference implementation is tracked separately in this same repo's
+git history; algorithmic semantics are unchanged across the port.
+
+## [0.1.0] — Unreleased
+
+Initial Python/Rust port of the MATLAB L1TV library by Storath, Weinmann
+(2016).
+
+### Added
+
+- Rust core (`_l1tv_core` cdylib, pyo3 0.24.1, numpy 0.24, ndarray 0.16):
+  - `solve_l1_tv_real(y, alpha, weights, use_dist_transform)` — exact
+    L1-TV solver for real-valued 1-D signals.
+  - `solve_l1_tv_circ(y, alpha, weights, use_dist_transform)` — exact
+    L1-TV solver for circle-valued 1-D signals with shortest-arc distance.
+  - `dist_transform.l1_dt_with_argmin` — forward + backward L1 distance
+    transform with argmin tracking (used by both solvers).
+- Python facade (`l1tv` package):
+  - `min_l1_tv` and `min_l1_tv_circ` public functions.
+  - Boundary handling: any 1-D shape (list, tuple, row, column) is
+    accepted via `np.ascontiguousarray + .ravel()`. Validates finite y,
+    non-negative alpha, shape-matched non-negative weights.
+- Tests:
+  - 25 Rust unit tests covering DT correctness, candidate-set theorem,
+    edge cases, and naive-vs-DT objective parity for both solvers.
+  - 28 Python tests covering the public API, input validation, and shape
+    handling.
+  - 16 MATLAB-vs-Rust objective parity tests (parametrised) using
+    pre-generated `.mat` fixtures; relerr < 1e-10 on every fixture.
+- Demos: `demos_python/demo_real.py` and `demo_circ.py` mirror the
+  original MATLAB demos and run end-to-end with optional matplotlib
+  plotting.
+
+### Memory layout
+
+The Rust port uses a rolling pair of K-cost columns plus a `u32` argmin-
+pointer table for O(N) backtracking, replacing MATLAB's full `(K, N)`
+double-precision cost table. Net memory roughly halves at typical sizes
+and can drop further if argmin pointers shrink to `u16` (when `K < 65536`).
+
+### Differences from the MATLAB reference
+
+- Input shape: silently produced wrong results for row-vector input in
+  `L1TV_Circ.m` (the antipodal vertical concat became a 2×K matrix). The
+  Python facade rules this out by construction.
+- Numerical type: `f64` only. Cast at the boundary if you need `f32`.
+- `setPath.m`-style invasive path setup is not relevant; just
+  `pip install l1tv`.
+
+These differences do not affect algorithmic semantics for column-vector
+MATLAB inputs, which is the only shape the MATLAB demos used.
+
+[0.1.0]: https://github.com/mstorath/L1TV/releases/tag/v0.1.0

l1tv-0.1.0/CITATION.cff

@@ -0,0 +1,41 @@
+cff-version: 1.2.0
+message: "If you use this software, please cite both the software and the associated paper below."
+title: "L1TV: Exact L1-TV regularization of real- or circle-valued signals"
+abstract: "Reference MATLAB implementation of exact, non-iterative algorithms for L1-TV regularization of real-valued and circle-valued (phase / orientation) signals. Solves the discrete L1-TV minimisation problem to global optimum via dynamic programming with a distance-transform speedup."
+type: software
+url: "https://github.com/mstorath/L1TV"
+repository-code: "https://github.com/mstorath/L1TV"
+license: MIT
+authors:
+  - family-names: Storath
+    given-names: Martin
+    email: martin.storath@thws.de
+    orcid: "https://orcid.org/0000-0003-1427-0776"
+    affiliation: "Lab for Mathematical Methods in Computer Vision and Machine Learning, Technische Hochschule Würzburg-Schweinfurt"
+  - family-names: Weinmann
+    given-names: Andreas
+    email: andreas.weinmann@thws.de
+    orcid: "https://orcid.org/0000-0002-4969-7609"
+    affiliation: "Department of Mathematics and Natural Sciences, Hochschule Darmstadt"
+preferred-citation:
+  type: article
+  title: "Exact algorithms for L1-TV regularization of real-valued or circle-valued signals"
+  authors:
+    - family-names: Storath
+      given-names: Martin
+      orcid: "https://orcid.org/0000-0003-1427-0776"
+    - family-names: Weinmann
+      given-names: Andreas
+      orcid: "https://orcid.org/0000-0002-4969-7609"
+    - family-names: Unser
+      given-names: Michael
+      orcid: "https://orcid.org/0000-0003-1248-2513"
+  journal: "SIAM Journal on Scientific Computing"
+  volume: 38
+  issue: 1
+  start: "A614"
+  end: "A630"
+  year: 2016
+  publisher:
+    name: "Society for Industrial and Applied Mathematics"
+  doi: "10.1137/15M101796X"

l1tv-0.1.0/Cargo.lock

@@ -0,0 +1,230 @@
+# This file is automatically @generated by Cargo.
+# It is not intended for manual editing.
+version = 4
+
+[[package]]
+name = "_l1tv_core"
+version = "0.1.0"
+dependencies = [
+ "ndarray",
+ "numpy",
+ "pyo3",
+]
+
+[[package]]
+name = "autocfg"
+version = "1.5.0"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "c08606f8c3cbf4ce6ec8e28fb0014a2c086708fe954eaa885384a6165172e7e8"
+
+[[package]]
+name = "heck"
+version = "0.5.0"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "2304e00983f87ffb38b55b444b5e3b60a884b5d30c0fca7d82fe33449bbe55ea"
+
+[[package]]
+name = "libc"
+version = "0.2.186"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "68ab91017fe16c622486840e4c83c9a37afeff978bd239b5293d61ece587de66"
+
+[[package]]
+name = "matrixmultiply"
+version = "0.3.10"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "a06de3016e9fae57a36fd14dba131fccf49f74b40b7fbdb472f96e361ec71a08"
+dependencies = [
+ "autocfg",
+ "rawpointer",
+]
+
+[[package]]
+name = "ndarray"
+version = "0.16.1"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "882ed72dce9365842bf196bdeedf5055305f11fc8c03dee7bb0194a6cad34841"
+dependencies = [
+ "matrixmultiply",
+ "num-complex",
+ "num-integer",
+ "num-traits",
+ "portable-atomic",
+ "portable-atomic-util",
+ "rawpointer",
+]
+
+[[package]]
+name = "num-complex"
+version = "0.4.6"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "73f88a1307638156682bada9d7604135552957b7818057dcef22705b4d509495"
+dependencies = [
+ "num-traits",
+]
+
+[[package]]
+name = "num-integer"
+version = "0.1.46"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "7969661fd2958a5cb096e56c8e1ad0444ac2bbcd0061bd28660485a44879858f"
+dependencies = [
+ "num-traits",
+]
+
+[[package]]
+name = "num-traits"
+version = "0.2.19"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "071dfc062690e90b734c0b2273ce72ad0ffa95f0c74596bc250dcfd960262841"
+dependencies = [
+ "autocfg",
+]
+
+[[package]]
+name = "numpy"
+version = "0.28.0"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "778da78c64ddc928ebf5ad9df5edf0789410ff3bdbf3619aed51cd789a6af1e2"
+dependencies = [
+ "libc",
+ "ndarray",
+ "num-complex",
+ "num-integer",
+ "num-traits",
+ "pyo3",
+ "pyo3-build-config",
+ "rustc-hash",
+]
+
+[[package]]
+name = "once_cell"
+version = "1.21.4"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "9f7c3e4beb33f85d45ae3e3a1792185706c8e16d043238c593331cc7cd313b50"
+
+[[package]]
+name = "portable-atomic"
+version = "1.13.1"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "c33a9471896f1c69cecef8d20cbe2f7accd12527ce60845ff44c153bb2a21b49"
+
+[[package]]
+name = "portable-atomic-util"
+version = "0.2.7"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "c2a106d1259c23fac8e543272398ae0e3c0b8d33c88ed73d0cc71b0f1d902618"
+dependencies = [
+ "portable-atomic",
+]
+
+[[package]]
+name = "proc-macro2"
+version = "1.0.106"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "8fd00f0bb2e90d81d1044c2b32617f68fcb9fa3bb7640c23e9c748e53fb30934"
+dependencies = [
+ "unicode-ident",
+]
+
+[[package]]
+name = "pyo3"
+version = "0.28.3"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "91fd8e38a3b50ed1167fb981cd6fd60147e091784c427b8f7183a7ee32c31c12"
+dependencies = [
+ "libc",
+ "once_cell",
+ "portable-atomic",
+ "pyo3-build-config",
+ "pyo3-ffi",
+ "pyo3-macros",
+]
+
+[[package]]
+name = "pyo3-build-config"
+version = "0.28.3"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "e368e7ddfdeb98c9bca7f8383be1648fd84ab466bf2bc015e94008db6d35611e"
+dependencies = [
+ "target-lexicon",
+]
+
+[[package]]
+name = "pyo3-ffi"
+version = "0.28.3"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "7f29e10af80b1f7ccaf7f69eace800a03ecd13e883acfacc1e5d0988605f651e"
+dependencies = [
+ "libc",
+ "pyo3-build-config",
+]
+
+[[package]]
+name = "pyo3-macros"
+version = "0.28.3"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "df6e520eff47c45997d2fc7dd8214b25dd1310918bbb2642156ef66a67f29813"
+dependencies = [
+ "proc-macro2",
+ "pyo3-macros-backend",
+ "quote",
+ "syn",
+]
+
+[[package]]
+name = "pyo3-macros-backend"
+version = "0.28.3"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "c4cdc218d835738f81c2338f822078af45b4afdf8b2e33cbb5916f108b813acb"
+dependencies = [
+ "heck",
+ "proc-macro2",
+ "pyo3-build-config",
+ "quote",
+ "syn",
+]
+
+[[package]]
+name = "quote"
+version = "1.0.45"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "41f2619966050689382d2b44f664f4bc593e129785a36d6ee376ddf37259b924"
+dependencies = [
+ "proc-macro2",
+]
+
+[[package]]
+name = "rawpointer"
+version = "0.2.1"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "60a357793950651c4ed0f3f52338f53b2f809f32d83a07f72909fa13e4c6c1e3"
+
+[[package]]
+name = "rustc-hash"
+version = "2.1.2"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "94300abf3f1ae2e2b8ffb7b58043de3d399c73fa6f4b73826402a5c457614dbe"
+
+[[package]]
+name = "syn"
+version = "2.0.117"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "e665b8803e7b1d2a727f4023456bbbbe74da67099c585258af0ad9c5013b9b99"
+dependencies = [
+ "proc-macro2",
+ "quote",
+ "unicode-ident",
+]
+
+[[package]]
+name = "target-lexicon"
+version = "0.13.5"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "adb6935a6f5c20170eeceb1a3835a49e12e19d792f6dd344ccc76a985ca5a6ca"
+
+[[package]]
+name = "unicode-ident"
+version = "1.0.24"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "e6e4313cd5fcd3dad5cafa179702e2b244f760991f45397d14d4ebf38247da75"

l1tv-0.1.0/Cargo.toml

@@ -0,0 +1,41 @@
+[package]
+name = "_l1tv_core"
+version = "0.1.0"
+edition = "2021"
+description = "Rust core for the l1tv Python package — exact L1-TV regularisation of real- or circle-valued 1-D signals (Storath, Weinmann, Unser 2016). Automated port of the original MATLAB by Claude Opus coding agent (Anthropic, 2026)."
+authors = [
+    "Martin Storath <storath@gmail.com> (original MATLAB)",
+    "Andreas Weinmann (original MATLAB)",
+    "Claude Opus coding agent, Anthropic (Rust port, 2026)"
+]
+license = "MIT"
+repository = "https://github.com/mstorath/L1TV"
+readme = "README_PYTHON.md"
+
+[lib]
+name = "_l1tv_core"
+crate-type = ["cdylib"]
+
+[dependencies]
+# pyo3 0.28 (June 2025+) is the first family of versions that supports
+# CPython 3.14. The previous family-wide pin at 0.24.1 (a security-fix pin
+# for GHSA-pph8-gcv7-4qj5) does not support 3.14 and broke the linux release
+# build with `--find-interpreter` after the manylinux containers added 3.14.
+# numpy crate version tracks pyo3 minor; ndarray version tracks numpy's.
+pyo3 = { version = "0.28" }
+numpy = "0.28"
+ndarray = "0.16"
+
+# `extension-module` is gated behind a Cargo feature instead of being on
+# unconditionally. Reason: with the feature ON, pyo3 instructs the linker to
+# leave Python symbols unresolved (resolved at runtime by the Python
+# interpreter that loads the .so). That works for the wheel target but breaks
+# `cargo test`, whose test binary has no Python interpreter to satisfy the
+# undefined symbols at runtime — macOS's strict linker rejects the binary.
+# Maturin enables the feature via [tool.maturin].features in pyproject.toml.
+[features]
+extension-module = ["pyo3/extension-module"]
+
+[profile.release]
+opt-level = 3
+lto = true