Automatic standalone structural SubDomainer

.github/workflows/leapfrog-parity.yml

@@ -2,6 +2,10 @@ name: Leapfrog parity
 
 on:
   push:
+    branches:
+      - feature/leapfrog-parity
+      - feature/automatic-subdomainer
+  pull_request:
     branches:
       - feature/leapfrog-parity
   workflow_dispatch:
@@ -19,6 +23,20 @@ jobs:
         with:
           python-version: "3.12"
 
+      - name: Validate standalone notebook
+        run: |
+          python - <<'PY'
+          import json
+          from pathlib import Path
+
+          path = Path("examples/standalone_automatic_lva.ipynb")
+          notebook = json.loads(path.read_text(encoding="utf-8"))
+          for index, cell in enumerate(notebook.get("cells", [])):
+              if cell.get("cell_type") == "code":
+                  compile("".join(cell.get("source", [])), f"{path}:cell-{index}", "exec")
+          print(f"Validated {path}")
+          PY
+
       - name: Install native build tools
         run: |
           sudo apt-get update
@@ -36,5 +54,6 @@ jobs:
         run: |
           python -m pip install pytest scipy
           python -m pytest -q \
+            python/tests/test_automatic_domain_builder.py \
             python/tests/test_labeled_domain_builder.py \
             python/tests/test_leapfrog_values.py

.github/workflows/polatory-benchmark.yml

@@ -0,0 +1,91 @@
+name: Polatory benchmark lab
+
+on:
+  workflow_dispatch:
+    inputs:
+      surface_resolution:
+        description: "Polatory surface resolution in metres"
+        required: false
+        default: "10"
+  push:
+    branches:
+      - feature/automatic-subdomainer
+
+jobs:
+  build-polatory:
+    runs-on: windows-latest
+    timeout-minutes: 300
+    env:
+      POLATORY_BENCHMARK_RESOLUTION: ${{ github.event.inputs.surface_resolution || '10' }}
+      PYTHONPATH: ${{ github.workspace }}\benchmarks\headless_shims
+
+    steps:
+      - name: Check out repository
+        uses: actions/checkout@v4
+        with:
+          submodules: recursive
+
+      - name: Set up Python 3.12
+        uses: actions/setup-python@v5
+        with:
+          python-version: "3.12"
+
+      - name: Install Python build tools
+        run: |
+          python -m pip install --upgrade pip
+          python -m pip install setuptools wheel cmake ninja build
+
+      - name: Bootstrap vcpkg
+        shell: powershell
+        run: |
+          .\vcpkg\bootstrap-vcpkg.bat -disableMetrics
+
+      - name: Build and install Polatory
+        shell: powershell
+        env:
+          VCPKG_ROOT: ${{ github.workspace }}\vcpkg
+        run: |
+          python -m pip install . --no-build-isolation --force-reinstall
+
+      - name: Install benchmark dependencies
+        run: |
+          python -m pip install scipy scikit-image pandas pyvista PyQt6 matplotlib
+
+      - name: Test Polatory import
+        run: |
+          python -c "import polatory; print('Polatory imported successfully')"
+
+      - name: Download Leapfrog gold benchmark release
+        shell: powershell
+        run: |
+          $archive = "benchmark-data\leapfrog-benchmark-data-v1.zip"
+          $folder = "benchmark-data\leapfrog-benchmark-data-v1"
+          New-Item -ItemType Directory -Force -Path "benchmark-data" | Out-Null
+          Invoke-WebRequest `
+            -Uri "https://github.com/Mining-Geologist/RBF/releases/download/leapfrog-benchmark-data-v1/leapfrog-benchmark-data-v1.zip" `
+            -OutFile $archive
+          New-Item -ItemType Directory -Force -Path $folder | Out-Null
+          Expand-Archive -Path $archive -DestinationPath $folder -Force
+          Get-ChildItem $folder
+
+      - name: Run synthetic structural-LVA suite
+        env:
+          QT_QPA_PLATFORM: offscreen
+        run: |
+          python benchmarks/leapfrog_blending/run_headless_suite.py
+
+      - name: Run real Leapfrog gold suite
+        env:
+          QT_QPA_PLATFORM: offscreen
+          LEAPFROG_GOLD_DIR: ${{ github.workspace }}\benchmark-data\leapfrog-benchmark-data-v1
+        run: |
+          python benchmarks/leapfrog_gold/run_real_gold_suite_ci.py
+
+      - name: Upload benchmark meshes, overlays, and metrics
+        if: always()
+        uses: actions/upload-artifact@v4
+        with:
+          name: polatory-leapfrog-benchmark-results
+          path: benchmark-results/
+          if-no-files-found: warn
+          retention-days: 30

benchmarks/headless_shims/polatory_lva_pyqt_app_v2.py

@@ -0,0 +1,108 @@
+"""Minimal non-interactive compatibility surface for CI benchmark imports.
+
+The production launcher chain historically imports the original PyQt v2 application
+at module import time, even when a benchmark only needs the automatic domain builder
+and chunk-safe mesher.  The full v2 GUI is intentionally not duplicated in the
+repository.  GitHub Actions prepends this directory to ``PYTHONPATH`` so the launcher
+patch modules can be imported headlessly without creating widgets.
+
+No benchmark calls these UI methods.  They exist only so v3-v8 can capture and patch
+the same attributes they patch in the interactive application.
+"""
+
+from __future__ import annotations
+
+from typing import Any
+
+import numpy as np
+import pandas as pd
+import pyvista as pv
+from PyQt6 import QtCore, QtWidgets
+
+
+APP_TITLE = "Polatory Automatic Structural LVA"
+ROLE_INSIDE = "Inside"
+ROLE_OUTSIDE = "Outside"
+ROLE_IGNORE = "Ignore"
+ROLE_OPTIONS = (ROLE_INSIDE, ROLE_OUTSIDE, ROLE_IGNORE)
+
+
+def category_keys(series: pd.Series) -> pd.Series:
+    return series.astype("string").fillna("<NA>").astype(str)
+
+
+def default_role_for_category(value: str) -> str:
+    text = str(value).strip().lower()
+    try:
+        number = float(text)
+    except ValueError:
+        number = np.nan
+    if np.isfinite(number):
+        if number > 0.0:
+            return ROLE_OUTSIDE
+        if number < 0.0:
+            return ROLE_INSIDE
+        return ROLE_IGNORE
+    if any(token in text for token in ("inside", "ore", "deposit", "target")):
+        return ROLE_INSIDE
+    if any(token in text for token in ("outside", "waste", "background")):
+        return ROLE_OUTSIDE
+    return ROLE_IGNORE
+
+
+def structured_points(
+    minimum: np.ndarray,
+    maximum: np.ndarray,
+    dimensions: tuple[int, int, int],
+) -> np.ndarray:
+    """Return the same regular Fortran-ordered point grid used by the GUI worker."""
+    minimum = np.asarray(minimum, dtype=float)
+    maximum = np.asarray(maximum, dtype=float)
+    dimensions = tuple(int(value) for value in dimensions)
+    axes = [
+        np.linspace(minimum[index], maximum[index], dimensions[index])
+        for index in range(3)
+    ]
+    x, y, z = np.meshgrid(*axes, indexing="ij")
+    return np.column_stack(
+        [x.ravel(order="F"), y.ravel(order="F"), z.ravel(order="F")]
+    )
+
+
+class ModelWorker:
+    def run(self) -> None:
+        raise RuntimeError("The CI compatibility worker must be patched before use.")
+
+
+class MainWindow:
+    """Attribute shell required only while the launcher patch modules are imported."""
+
+    def __init__(self, *args: Any, **kwargs: Any) -> None:
+        del args, kwargs
+
+    def model_parameters(self) -> dict[str, Any]:
+        return {}
+
+    def model_finished(self, result: dict[str, Any]) -> None:
+        del result
+
+    def load_csv(self) -> None:
+        return None
+
+    def run_model(self) -> None:
+        return None
+
+    def populate_category_roles(self, column: str) -> None:
+        del column
+
+    def _mapped_arrays(self):
+        raise RuntimeError("The headless compatibility shell has no mapped UI data.")
+
+    def apply_data_mapping(self) -> None:
+        return None
+
+
+def main() -> int:
+    raise RuntimeError(
+        "This module is a CI import shim, not an interactive application launcher."
+    )

benchmarks/leapfrog_blending/01_two_domain_blend_binary.csv

@@ -0,0 +1,25 @@
+X,Y,Z,Category,Indicator,DomainHint,ExpectedBoundaryZ
+-160,-60,60,Inside,-1,Left,80
+-160,-60,100,Outside,1,Left,80
+-160,0,60,Inside,-1,Left,80
+-160,0,100,Outside,1,Left,80
+-160,60,60,Inside,-1,Left,80
+-160,60,100,Outside,1,Left,80
+-80,-60,20,Inside,-1,Left,40
+-80,-60,60,Outside,1,Left,40
+-80,0,20,Inside,-1,Left,40
+-80,0,60,Outside,1,Left,40
+-80,60,20,Inside,-1,Left,40
+-80,60,60,Outside,1,Left,40
+80,-60,120,Inside,-1,Right,140
+80,-60,160,Outside,1,Right,140
+80,0,120,Inside,-1,Right,140
+80,0,160,Outside,1,Right,140
+80,60,120,Inside,-1,Right,140
+80,60,160,Outside,1,Right,140
+160,-60,160,Inside,-1,Right,180
+160,-60,200,Outside,1,Right,180
+160,0,160,Inside,-1,Right,180
+160,0,200,Outside,1,Right,180
+160,60,160,Inside,-1,Right,180
+160,60,200,Outside,1,Right,180

benchmarks/leapfrog_blending/01_two_domain_blend_points.csv

@@ -0,0 +1,37 @@
+X,Y,Z,Role,Indicator,DomainHint,ContactSurfaceZ
+-160,-60,60,Inside,-1,Left,80
+-160,-60,80,Contact,0,Left,80
+-160,-60,100,Outside,1,Left,80
+-160,0,60,Inside,-1,Left,80
+-160,0,80,Contact,0,Left,80
+-160,0,100,Outside,1,Left,80
+-160,60,60,Inside,-1,Left,80
+-160,60,80,Contact,0,Left,80
+-160,60,100,Outside,1,Left,80
+-80,-60,20,Inside,-1,Left,40
+-80,-60,40,Contact,0,Left,40
+-80,-60,60,Outside,1,Left,40
+-80,0,20,Inside,-1,Left,40
+-80,0,40,Contact,0,Left,40
+-80,0,60,Outside,1,Left,40
+-80,60,20,Inside,-1,Left,40
+-80,60,40,Contact,0,Left,40
+-80,60,60,Outside,1,Left,40
+80,-60,120,Inside,-1,Right,140
+80,-60,140,Contact,0,Right,140
+80,-60,160,Outside,1,Right,140
+80,0,120,Inside,-1,Right,140
+80,0,140,Contact,0,Right,140
+80,0,160,Outside,1,Right,140
+80,60,120,Inside,-1,Right,140
+80,60,140,Contact,0,Right,140
+80,60,160,Outside,1,Right,140
+160,-60,160,Inside,-1,Right,180
+160,-60,180,Contact,0,Right,180
+160,-60,200,Outside,1,Right,180
+160,0,160,Inside,-1,Right,180
+160,0,180,Contact,0,Right,180
+160,0,200,Outside,1,Right,180
+160,60,160,Inside,-1,Right,180
+160,60,180,Contact,0,Right,180
+160,60,200,Outside,1,Right,180

benchmarks/leapfrog_blending/02_background_decay_binary.csv

@@ -0,0 +1,31 @@
+X,Y,Z,Category,Indicator,ExpectedBoundaryZ
+-160,-60,60,Inside,-1,80
+-160,-60,100,Outside,1,80
+-160,0,60,Inside,-1,80
+-160,0,100,Outside,1,80
+-160,60,60,Inside,-1,80
+-160,60,100,Outside,1,80
+-80,-60,20,Inside,-1,40
+-80,-60,60,Outside,1,40
+-80,0,20,Inside,-1,40
+-80,0,60,Outside,1,40
+-80,60,20,Inside,-1,40
+-80,60,60,Outside,1,40
+0,-60,-20,Inside,-1,0
+0,-60,20,Outside,1,0
+0,0,-20,Inside,-1,0
+0,0,20,Outside,1,0
+0,60,-20,Inside,-1,0
+0,60,20,Outside,1,0
+80,-60,20,Inside,-1,40
+80,-60,60,Outside,1,40
+80,0,20,Inside,-1,40
+80,0,60,Outside,1,40
+80,60,20,Inside,-1,40
+80,60,60,Outside,1,40
+160,-60,60,Inside,-1,80
+160,-60,100,Outside,1,80
+160,0,60,Inside,-1,80
+160,0,100,Outside,1,80
+160,60,60,Inside,-1,80
+160,60,100,Outside,1,80

benchmarks/leapfrog_blending/02_background_decay_points.csv

@@ -0,0 +1,46 @@
+X,Y,Z,Role,Indicator,ContactSurfaceZ
+-160,-60,60,Inside,-1,80
+-160,-60,80,Contact,0,80
+-160,-60,100,Outside,1,80
+-160,0,60,Inside,-1,80
+-160,0,80,Contact,0,80
+-160,0,100,Outside,1,80
+-160,60,60,Inside,-1,80
+-160,60,80,Contact,0,80
+-160,60,100,Outside,1,80
+-80,-60,20,Inside,-1,40
+-80,-60,40,Contact,0,40
+-80,-60,60,Outside,1,40
+-80,0,20,Inside,-1,40
+-80,0,40,Contact,0,40
+-80,0,60,Outside,1,40
+-80,60,20,Inside,-1,40
+-80,60,40,Contact,0,40
+-80,60,60,Outside,1,40
+0,-60,-20,Inside,-1,0
+0,-60,0,Contact,0,0
+0,-60,20,Outside,1,0
+0,0,-20,Inside,-1,0
+0,0,0,Contact,0,0
+0,0,20,Outside,1,0
+0,60,-20,Inside,-1,0
+0,60,0,Contact,0,0
+0,60,20,Outside,1,0
+80,-60,20,Inside,-1,40
+80,-60,40,Contact,0,40
+80,-60,60,Outside,1,40
+80,0,20,Inside,-1,40
+80,0,40,Contact,0,40
+80,0,60,Outside,1,40
+80,60,20,Inside,-1,40
+80,60,40,Contact,0,40
+80,60,60,Outside,1,40
+160,-60,60,Inside,-1,80
+160,-60,80,Contact,0,80
+160,-60,100,Outside,1,80
+160,0,60,Inside,-1,80
+160,0,80,Contact,0,80
+160,0,100,Outside,1,80
+160,60,60,Inside,-1,80
+160,60,80,Contact,0,80
+160,60,100,Outside,1,80