pylibstats 0.1.5__tar.gz

pylibstats-0.1.5/.github/workflows/ci.yml

@@ -0,0 +1,33 @@
+name: CI
+
+on:
+  push:
+    branches: [main]
+  pull_request:
+    branches: [main]
+
+concurrency:
+  group: ci-${{ github.ref }}
+  cancel-in-progress: true
+
+jobs:
+  build-and-test:
+    name: ${{ matrix.os }} / Python ${{ matrix.python }}
+    runs-on: ${{ matrix.os }}
+    strategy:
+      fail-fast: false
+      matrix:
+        os: [ubuntu-latest, macos-latest, windows-latest]
+        python: ["3.11", "3.12", "3.13"]
+    steps:
+      - uses: actions/checkout@v4
+
+      - uses: actions/setup-python@v5
+        with:
+          python-version: ${{ matrix.python }}
+
+      - name: Build and install
+        run: pip install ".[test]" -v
+
+      - name: Run tests
+        run: pytest tests -v --tb=short

pylibstats-0.1.5/.github/workflows/wheels.yml

@@ -0,0 +1,104 @@
+name: Build wheels
+
+on:
+  push:
+    tags: ["v*"]
+  workflow_dispatch:
+
+concurrency:
+  group: wheels-${{ github.ref }}
+  cancel-in-progress: true
+
+jobs:
+  build-wheels:
+    name: Wheels / ${{ matrix.name }}
+    runs-on: ${{ matrix.os }}
+    strategy:
+      fail-fast: false
+      matrix:
+        include:
+          - name: linux-x86_64
+            os: ubuntu-latest
+            cibw_archs_linux: x86_64
+          - name: linux-aarch64
+            os: ubuntu-24.04-arm
+            cibw_archs_linux: aarch64
+          - name: macos
+            os: macos-latest
+            cibw_archs_linux: ""
+          - name: windows
+            os: windows-latest
+            cibw_archs_linux: ""
+
+    steps:
+      - uses: actions/checkout@v4
+
+      - uses: actions/setup-python@v5
+        with:
+          python-version: "3.12"
+
+      - name: Install cibuildwheel
+        run: pip install cibuildwheel
+
+      - name: Build wheels
+        run: cibuildwheel --output-dir wheelhouse
+        env:
+          # macOS: build universal2 (x86_64 + arm64)
+          CIBW_ARCHS_MACOS: "x86_64 arm64"
+          # nanobind uses aligned deallocation, which requires macOS 10.13+
+          CIBW_ENVIRONMENT_MACOS: "MACOSX_DEPLOYMENT_TARGET=10.13"
+          # Linux: set per-matrix (x86_64 on x64 runner, aarch64 on native ARM runner)
+          CIBW_ARCHS_LINUX: ${{ matrix.cibw_archs_linux }}
+          # Windows: x86_64 only
+          CIBW_ARCHS_WINDOWS: "AMD64"
+          # Skip PyPy, musllinux, and older Pythons
+          CIBW_SKIP: "*-musllinux_* cp39-* cp310-*"
+          # Build dependencies (installed in the build environment)
+          CIBW_BUILD_FRONTEND: "pip"
+          # Test the built wheel
+          CIBW_TEST_REQUIRES: "pytest scipy"
+          CIBW_TEST_COMMAND: "pytest {project}/tests -q"
+          # Ensure CMake is available in manylinux images
+          CIBW_BEFORE_ALL_LINUX: "yum install -y cmake || apt-get install -y cmake || true"
+
+      - uses: actions/upload-artifact@v4
+        with:
+          name: wheels-${{ matrix.name }}
+          path: wheelhouse/*.whl
+
+  build-sdist:
+    name: Source distribution
+    runs-on: ubuntu-latest
+    steps:
+      - uses: actions/checkout@v4
+
+      - uses: actions/setup-python@v5
+        with:
+          python-version: "3.12"
+
+      - name: Build sdist
+        run: pip install build && python -m build --sdist
+
+      - uses: actions/upload-artifact@v4
+        with:
+          name: sdist
+          path: dist/*.tar.gz
+
+  publish:
+    name: Publish to PyPI
+    needs: [build-wheels, build-sdist]
+    runs-on: ubuntu-latest
+    if: startsWith(github.ref, 'refs/tags/v')
+    environment:
+      name: pypi
+    permissions:
+      id-token: write
+    steps:
+      - uses: actions/download-artifact@v4
+        with:
+          path: dist
+          merge-multiple: true
+
+      - uses: pypa/gh-action-pypi-publish@release/v1
+        with:
+          packages-dir: dist

pylibstats-0.1.5/.gitignore

@@ -0,0 +1,29 @@
+# Build artifacts
+build/
+_skbuild/
+dist/
+*.egg-info/
+*.so
+*.pyd
+*.dll
+
+# Python
+__pycache__/
+*.pyc
+.venv/
+venv/
+
+# CMake
+CMakeCache.txt
+CMakeFiles/
+cmake_install.cmake
+Makefile
+
+# IDE
+.vscode/
+.idea/
+*.swp
+
+# OS
+.DS_Store
+Thumbs.db

pylibstats-0.1.5/CMakeLists.txt

@@ -0,0 +1,79 @@
+cmake_minimum_required(VERSION 3.20)
+project(pylibstats LANGUAGES CXX)
+
+set(CMAKE_CXX_STANDARD 20)
+set(CMAKE_CXX_STANDARD_REQUIRED ON)
+set(CMAKE_CXX_EXTENSIONS OFF)
+
+# ---------------------------------------------------------------------------
+# Python + nanobind
+# ---------------------------------------------------------------------------
+find_package(Python REQUIRED COMPONENTS Interpreter Development.Module)
+
+# scikit-build-core adds site-packages to CMAKE_PREFIX_PATH automatically.
+# For direct CMake invocation, discover nanobind via Python as a fallback.
+if(NOT nanobind_ROOT AND NOT nanobind_DIR)
+    execute_process(
+        COMMAND "${Python_EXECUTABLE}" -m nanobind --cmake_dir
+        OUTPUT_STRIP_TRAILING_WHITESPACE OUTPUT_VARIABLE nanobind_ROOT
+        RESULT_VARIABLE _nanobind_rc)
+    if(NOT _nanobind_rc EQUAL 0)
+        set(nanobind_ROOT "")
+    endif()
+endif()
+
+find_package(nanobind CONFIG REQUIRED)
+
+# ---------------------------------------------------------------------------
+# libstats — prefer installed package, fall back to FetchContent
+# ---------------------------------------------------------------------------
+find_package(libstats 1.1.5 QUIET)
+
+if(NOT libstats_FOUND)
+    message(STATUS "libstats not found via find_package — using FetchContent")
+    include(FetchContent)
+    FetchContent_Declare(
+        libstats
+        GIT_REPOSITORY https://github.com/OldCrow/libstats.git
+        GIT_TAG        v1.1.5
+    )
+    # Disable libstats tests/tools/examples when building as a subdirectory.
+    set(LIBSTATS_BUILD_TESTS OFF CACHE BOOL "" FORCE)
+    set(LIBSTATS_BUILD_TOOLS OFF CACHE BOOL "" FORCE)
+    set(LIBSTATS_BUILD_EXAMPLES OFF CACHE BOOL "" FORCE)
+    set(BUILD_DOCS OFF CACHE BOOL "" FORCE)
+    FetchContent_MakeAvailable(libstats)
+    set(LIBSTATS_TARGET libstats_static)
+else()
+    message(STATUS "Using installed libstats ${LIBSTATS_VERSION}")
+    # Installed exports use libstats::libstats_static (not libstats::static)
+    if(TARGET libstats::libstats_static)
+        set(LIBSTATS_TARGET libstats::libstats_static)
+    elseif(TARGET libstats::static)
+        set(LIBSTATS_TARGET libstats::static)
+    elseif(DEFINED LIBSTATS_LIBRARIES)
+        set(LIBSTATS_TARGET ${LIBSTATS_LIBRARIES})
+    else()
+        message(FATAL_ERROR "libstats found but no usable static target")
+    endif()
+endif()
+
+# ---------------------------------------------------------------------------
+# Extension module
+# ---------------------------------------------------------------------------
+nanobind_add_module(
+    _core
+    STABLE_ABI
+    NB_STATIC
+    src/pylibstats/_core.cpp
+)
+
+target_link_libraries(_core PRIVATE ${LIBSTATS_TARGET})
+
+# Windows: prevent min/max macros from <windows.h> colliding with std::min/max
+if(WIN32)
+    target_compile_definitions(_core PRIVATE NOMINMAX _USE_MATH_DEFINES)
+endif()
+
+# Let nanobind / scikit-build-core handle the install destination
+install(TARGETS _core LIBRARY DESTINATION pylibstats)

pylibstats-0.1.5/PKG-INFO

@@ -0,0 +1,111 @@
+Metadata-Version: 2.2
+Name: pylibstats
+Version: 0.1.5
+Summary: Python bindings for libstats — SIMD-accelerated statistical distributions via NumPy
+License: MIT
+Requires-Python: >=3.11
+Requires-Dist: numpy>=1.20
+Provides-Extra: test
+Requires-Dist: pytest>=7.0; extra == "test"
+Requires-Dist: scipy; extra == "test"
+Description-Content-Type: text/markdown
+
+# pylibstats
+
+Python bindings for [libstats](https://github.com/OldCrow/libstats) — a C++20 statistical distributions library with SIMD batch operations.
+
+## Features
+
+- **9 distributions**: Gaussian, Exponential, Uniform, Poisson, Discrete Uniform, Gamma, Beta, Chi-Squared, Student's t
+- **NumPy integration**: pass arrays directly to `pdf()`, `cdf()`, `log_pdf()` — the SIMD/parallel batch path runs automatically
+- **GIL-releasing**: batch operations release the Python GIL for concurrent workloads
+- **SciPy-compatible naming**: `pdf`, `cdf`, `ppf`, `fit`, `sample`
+- **Input validation**: all constructor, setter, and `fit()` parameters are validated in Python with clear `ValueError` messages
+
+## Quick start
+
+```python
+import numpy as np
+import pylibstats
+
+dist = pylibstats.Gaussian(mu=0.0, sigma=1.0)
+
+# Scalar
+dist.pdf(1.0)
+dist.cdf(0.0)        # 0.5
+dist.ppf(0.975)      # ~1.96
+
+# Batch (SIMD-accelerated)
+x = np.linspace(-4, 4, 100_000)
+densities = dist.pdf(x)
+
+# Sampling
+samples = dist.sample(n=10_000, seed=42)
+
+# Fitting
+dist.fit(samples)
+```
+
+## Building from source
+
+Requires Python ≥3.11, CMake ≥3.20, and a C++20 compiler.
+
+```bash
+pip install .
+```
+
+This fetches libstats v1.1.5 via CMake FetchContent if not already installed.
+
+### Building against a local libstats
+
+To link against a locally built libstats (e.g. a development branch), install
+libstats to a prefix and point `pip` at it:
+
+```bash
+# In the libstats repo
+cmake --install build --prefix /path/to/libstats/install
+
+# In this repo — use libstats_DIR, not CMAKE_PREFIX_PATH
+# (overriding CMAKE_PREFIX_PATH breaks nanobind discovery)
+pip install --no-build-isolation -ve . \
+    -Ccmake.define.libstats_DIR=/path/to/libstats/install/lib/cmake/libstats
+```
+
+`--no-build-isolation` requires build deps in the active environment:
+
+```bash
+pip install "scikit-build-core>=0.10" "nanobind>=2.0"
+```
+
+## Running tests
+
+```bash
+pip install ".[test]"
+pytest
+```
+
+## Examples
+
+See the `examples/` directory:
+
+- `basic_usage.py` — scalar/batch operations, sampling, and fitting
+- `benchmark.py` — wall-clock comparison against SciPy (PDF and CDF)
+- `scipy_comparison.py` — numerical accuracy verification across all 9 distributions
+
+## Known limitations
+
+- **Beta CDF performance**: the regularised incomplete beta function in libstats is slower than SciPy's implementation (~0.5× speedup). All other distribution/operation combinations are faster.
+
+## Contributing
+
+### macOS ABI note
+
+On macOS, libstats may be compiled with Homebrew LLVM while Python extensions use Apple clang. These ship different `libc++` versions whose exception-handling ABIs are incompatible — C++ exceptions thrown from libstats segfault during stack unwinding instead of propagating normally.
+
+pylibstats works around this by validating all parameters in pure Python (in `__init__.py`) *before* calling into the C++ layer, so the error path never crosses the ABI boundary. If you add new parameters or distribution classes, follow the same pattern: validate in Python, then delegate to `_core`.
+
+See `libstats/include/core/error_handling.h` for the upstream discussion.
+
+## License
+
+MIT

pylibstats-0.1.5/README.md

@@ -0,0 +1,99 @@
+# pylibstats
+
+Python bindings for [libstats](https://github.com/OldCrow/libstats) — a C++20 statistical distributions library with SIMD batch operations.
+
+## Features
+
+- **9 distributions**: Gaussian, Exponential, Uniform, Poisson, Discrete Uniform, Gamma, Beta, Chi-Squared, Student's t
+- **NumPy integration**: pass arrays directly to `pdf()`, `cdf()`, `log_pdf()` — the SIMD/parallel batch path runs automatically
+- **GIL-releasing**: batch operations release the Python GIL for concurrent workloads
+- **SciPy-compatible naming**: `pdf`, `cdf`, `ppf`, `fit`, `sample`
+- **Input validation**: all constructor, setter, and `fit()` parameters are validated in Python with clear `ValueError` messages
+
+## Quick start
+
+```python
+import numpy as np
+import pylibstats
+
+dist = pylibstats.Gaussian(mu=0.0, sigma=1.0)
+
+# Scalar
+dist.pdf(1.0)
+dist.cdf(0.0)        # 0.5
+dist.ppf(0.975)      # ~1.96
+
+# Batch (SIMD-accelerated)
+x = np.linspace(-4, 4, 100_000)
+densities = dist.pdf(x)
+
+# Sampling
+samples = dist.sample(n=10_000, seed=42)
+
+# Fitting
+dist.fit(samples)
+```
+
+## Building from source
+
+Requires Python ≥3.11, CMake ≥3.20, and a C++20 compiler.
+
+```bash
+pip install .
+```
+
+This fetches libstats v1.1.5 via CMake FetchContent if not already installed.
+
+### Building against a local libstats
+
+To link against a locally built libstats (e.g. a development branch), install
+libstats to a prefix and point `pip` at it:
+
+```bash
+# In the libstats repo
+cmake --install build --prefix /path/to/libstats/install
+
+# In this repo — use libstats_DIR, not CMAKE_PREFIX_PATH
+# (overriding CMAKE_PREFIX_PATH breaks nanobind discovery)
+pip install --no-build-isolation -ve . \
+    -Ccmake.define.libstats_DIR=/path/to/libstats/install/lib/cmake/libstats
+```
+
+`--no-build-isolation` requires build deps in the active environment:
+
+```bash
+pip install "scikit-build-core>=0.10" "nanobind>=2.0"
+```
+
+## Running tests
+
+```bash
+pip install ".[test]"
+pytest
+```
+
+## Examples
+
+See the `examples/` directory:
+
+- `basic_usage.py` — scalar/batch operations, sampling, and fitting
+- `benchmark.py` — wall-clock comparison against SciPy (PDF and CDF)
+- `scipy_comparison.py` — numerical accuracy verification across all 9 distributions
+
+## Known limitations
+
+- **Beta CDF performance**: the regularised incomplete beta function in libstats is slower than SciPy's implementation (~0.5× speedup). All other distribution/operation combinations are faster.
+
+## Contributing
+
+### macOS ABI note
+
+On macOS, libstats may be compiled with Homebrew LLVM while Python extensions use Apple clang. These ship different `libc++` versions whose exception-handling ABIs are incompatible — C++ exceptions thrown from libstats segfault during stack unwinding instead of propagating normally.
+
+pylibstats works around this by validating all parameters in pure Python (in `__init__.py`) *before* calling into the C++ layer, so the error path never crosses the ABI boundary. If you add new parameters or distribution classes, follow the same pattern: validate in Python, then delegate to `_core`.
+
+See `libstats/include/core/error_handling.h` for the upstream discussion.
+
+## License
+
+MIT

pylibstats-0.1.5/WARP.md

@@ -0,0 +1,50 @@
+# WARP.md
+
+Guidance for working in `pylibstats`.
+
+## Project purpose
+
+`pylibstats` provides Python bindings for `libstats` via `nanobind` and `scikit-build-core`.
+
+Core goals:
+
+- expose the `libstats` distribution API in Python
+- keep NumPy-based batch paths fast and simple
+- keep Python wrappers, native bindings, and stubs synchronized
+
+## Key files
+
+- `CMakeLists.txt` — native extension build and `libstats` dependency wiring
+- `pyproject.toml` — package metadata and build backend config
+- `src/pylibstats/_core.cpp` — nanobind bindings
+- `src/pylibstats/_common.h` — NumPy conversion helpers
+- `src/pylibstats/__init__.py` — Python wrappers and validation
+- `src/pylibstats/__init__.pyi`, `src/pylibstats/_core.pyi` — typing stubs
+- `tests/` — pytest suite
+
+## Dependency notes
+
+- Build first tries `find_package(libstats)`.
+- If not found, CMake fetches `libstats` from GitHub tag `v1.1.5`.
+- For local development against a custom `libstats` install, pass `libstats_DIR` (do not override `CMAKE_PREFIX_PATH`, which can break nanobind discovery).
+
+## Common commands
+
+```powershell
+python -m pip install -e C:\Users\gdwol\Development\pylibstats
+python -m pip install ".[test]"
+python -m pytest C:\Users\gdwol\Development\pylibstats\tests -q
+```
+
+Local `libstats` build override example:
+
+```powershell
+python -m pip install --no-build-isolation -ve C:\Users\gdwol\Development\pylibstats `
+  -Ccmake.define.libstats_DIR=C:\path\to\libstats\install\lib\cmake\libstats
+```
+
+## Editing rules
+
+1. Keep `_core.cpp`, `__init__.py`, and `.pyi` stubs consistent.
+2. Add or update tests for any behavior/API change.
+3. Keep docs concise and accurate.

pylibstats-0.1.5/examples/basic_usage.py

@@ -0,0 +1,78 @@
+"""Basic pylibstats usage examples."""
+
+import numpy as np
+
+import pylibstats
+
+# ── Construction ─────────────────────────────────────────────────────────────
+
+gaussian = pylibstats.Gaussian(mu=0.0, sigma=1.0)
+gamma = pylibstats.Gamma(alpha=2.0, beta=0.5)
+beta = pylibstats.Beta(alpha=2.0, beta=5.0)
+poisson = pylibstats.Poisson(lam=4.0)
+
+# Normal is an alias for Gaussian
+normal = pylibstats.Normal(mu=100.0, sigma=15.0)
+
+print("── Distributions ──")
+print(f"  {gaussian!r}")
+print(f"  {gamma!r}")
+print(f"  {beta!r}")
+print(f"  {poisson!r}")
+print(f"  {normal!r}")
+
+# ── Scalar operations ────────────────────────────────────────────────────────
+
+print("\n── Scalar PDF / CDF / PPF ──")
+print(f"  Gaussian PDF(0)   = {gaussian.pdf(0.0):.6f}")
+print(f"  Gaussian CDF(0)   = {gaussian.cdf(0.0):.6f}")
+print(f"  Gaussian PPF(0.975) = {gaussian.ppf(0.975):.4f}")
+print(f"  Gamma PDF(2)      = {gamma.pdf(2.0):.6f}")
+print(f"  Beta CDF(0.3)     = {beta.cdf(0.3):.6f}")
+print(f"  Poisson PMF(4)    = {poisson.pdf(4.0):.6f}")
+
+# ── Batch operations (SIMD-accelerated) ──────────────────────────────────────
+
+x = np.linspace(-4, 4, 10)
+pdf_values = gaussian.pdf(x)
+cdf_values = gaussian.cdf(x)
+
+print("\n── Batch operations (10 points) ──")
+print(f"  x   = {x}")
+print(f"  pdf = {pdf_values}")
+print(f"  cdf = {cdf_values}")
+
+# ── Moment properties ────────────────────────────────────────────────────────
+
+print("\n── Moment properties ──")
+for name, dist in [("Gaussian(0,1)", gaussian), ("Gamma(2,0.5)", gamma),
+                    ("Beta(2,5)", beta), ("Poisson(4)", poisson)]:
+    print(f"  {name:16s}  mean={dist.mean:.4f}  var={dist.variance:.4f}  "
+          f"skew={dist.skewness:.4f}  kurt={dist.kurtosis:.4f}")
+
+# ── Sampling ─────────────────────────────────────────────────────────────────
+
+samples = gaussian.sample(n=10_000, seed=42)
+print(f"\n── Sampling ──")
+print(f"  10,000 Gaussian samples: mean={samples.mean():.4f}, std={samples.std():.4f}")
+
+beta_samples = beta.sample(n=10_000, seed=42)
+print(f"  10,000 Beta(2,5) samples: mean={beta_samples.mean():.4f}, "
+      f"min={beta_samples.min():.4f}, max={beta_samples.max():.4f}")
+
+# ── Fitting ──────────────────────────────────────────────────────────────────
+
+data = np.random.default_rng(42).normal(loc=5.0, scale=2.0, size=5000)
+fitted = pylibstats.Gaussian()
+fitted.fit(data)
+print(f"\n── Fitting ──")
+print(f"  Fitted Gaussian to 5,000 samples from N(5, 2):")
+print(f"  mu={fitted.mu:.4f}  sigma={fitted.sigma:.4f}")
+
+# ── Parameter modification ───────────────────────────────────────────────────
+
+print(f"\n── Parameter modification ──")
+dist = pylibstats.Exponential(lam=1.0)
+print(f"  Before: lam={dist.lam}, mean={dist.mean:.4f}")
+dist.lam = 3.0
+print(f"  After:  lam={dist.lam}, mean={dist.mean:.4f}")

pylibstats-0.1.5/examples/benchmark.py

@@ -0,0 +1,110 @@
+"""Benchmark pylibstats batch operations against SciPy.
+
+Measures wall-clock time for PDF and CDF evaluation on arrays of 100k and 1M
+elements across all continuous distributions. Prints a formatted table with
+per-distribution speedup ratios.
+
+Usage:
+    python examples/benchmark.py
+"""
+
+import time
+
+import numpy as np
+from scipy import stats as sp
+
+import pylibstats
+
+
+def bench(fn, *, warmup: int = 2, repeats: int = 5) -> float:
+    """Return median wall-clock seconds for fn()."""
+    for _ in range(warmup):
+        fn()
+    times = []
+    for _ in range(repeats):
+        t0 = time.perf_counter()
+        fn()
+        times.append(time.perf_counter() - t0)
+    return float(np.median(times))
+
+
+# ── Distribution configs ─────────────────────────────────────────────────────
+# Each entry: (name, pylibstats instance, scipy frozen instance, x_range)
+
+DISTRIBUTIONS = [
+    ("Gaussian(0,1)",
+     pylibstats.Gaussian(0.0, 1.0),
+     sp.norm(0.0, 1.0),
+     (-4.0, 4.0)),
+
+    ("Exponential(1)",
+     pylibstats.Exponential(1.0),
+     sp.expon(scale=1.0),
+     (0.01, 10.0)),
+
+    ("Uniform(0,1)",
+     pylibstats.Uniform(0.0, 1.0),
+     sp.uniform(0.0, 1.0),
+     (0.0, 1.0)),
+
+    ("Gamma(2,1)",
+     pylibstats.Gamma(2.0, 1.0),
+     sp.gamma(a=2.0, scale=1.0),
+     (0.01, 10.0)),
+
+    ("Beta(2,5)",
+     pylibstats.Beta(2.0, 5.0),
+     sp.beta(a=2.0, b=5.0),
+     (0.01, 0.99)),
+
+    ("ChiSquared(5)",
+     pylibstats.ChiSquared(5.0),
+     sp.chi2(df=5),
+     (0.1, 20.0)),
+
+    ("StudentT(10)",
+     pylibstats.StudentT(10.0),
+     sp.t(df=10),
+     (-5.0, 5.0)),
+]
+
+SIZES = [100_000, 1_000_000]
+
+
+def run_benchmarks() -> None:
+    header = (f"{'Distribution':>20s}  {'N':>10s}  {'Op':>4s}  "
+              f"{'pylibstats (ms)':>15s}  {'SciPy (ms)':>12s}  {'Speedup':>8s}")
+    sep = "─" * len(header)
+
+    print()
+    print(header)
+    print(sep)
+
+    for name, pl_dist, sc_dist, (lo, hi) in DISTRIBUTIONS:
+        for n in SIZES:
+            x = np.linspace(lo, hi, n)
+
+            # PDF — bind loop vars via default args to avoid late-binding
+            t_pl = bench(lambda d=pl_dist, arr=x: d.pdf(arr)) * 1000
+            t_sc = bench(lambda d=sc_dist, arr=x: d.pdf(arr)) * 1000
+            ratio = t_sc / t_pl if t_pl > 0 else float("inf")
+            print(f"{name:>20s}  {n:>10,d}  {'PDF':>4s}  "
+                  f"{t_pl:>15.2f}  {t_sc:>12.2f}  {ratio:>7.1f}x")
+
+            # CDF
+            t_pl = bench(lambda d=pl_dist, arr=x: d.cdf(arr)) * 1000
+            t_sc = bench(lambda d=sc_dist, arr=x: d.cdf(arr)) * 1000
+            ratio = t_sc / t_pl if t_pl > 0 else float("inf")
+            print(f"{'':>20s}  {'':>10s}  {'CDF':>4s}  "
+                  f"{t_pl:>15.2f}  {t_sc:>12.2f}  {ratio:>7.1f}x")
+
+    print(sep)
+    print()
+
+
+if __name__ == "__main__":
+    import scipy
+    print("pylibstats vs SciPy batch performance benchmark")
+    print(f"NumPy {np.__version__}, SciPy {scipy.__version__}")
+    print(f"Repeats: 5 (median), Warmup: 2")
+    run_benchmarks()