processforge 0.3.4__tar.gz → 0.3.8__tar.gz

{processforge-0.3.4/src/processforge.egg-info → processforge-0.3.8}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: processforge
-Version: 0.3.4
+Version: 0.3.8
 Summary: A Python-based process simulation framework for chemical engineering applications.
 Author-email: Processforge Team <team@processforge.dev>
 License-Expression: BSD-3-Clause

{processforge-0.3.4 → processforge-0.3.8}/flowsheets/openmc/msre_eigenvalue.json

@@ -17,7 +17,7 @@
       "id": 1,
       "density": 2.26,
       "density_units": "g/cm3",
-      "temperature": 0.0,
+      "temperature": 300.0,
       "elements": [
         { "element": "C", "percent": 1.0, "percent_type": "ao" }
       ]
@@ -26,7 +26,7 @@
       "id": 2,
       "density": 9.0,
       "density_units": "g/cm3",
-      "temperature": 0.0,
+      "temperature": 400.0,
       "elements": [
         { "element": "Ni", "percent": 0.72, "percent_type": "wo" },
         { "element": "Mo", "percent": 0.16, "percent_type": "wo" },
@@ -38,7 +38,7 @@
       "id": 3,
       "density": 2.2,
       "density_units": "g/cm3",
-      "temperature": 0.0,
+      "temperature": 300.0,
       "depletable": true,
       "nuclides": [
         { "name": "Li7",  "percent": 0.335000, "percent_type": "ao" },
@@ -55,7 +55,7 @@
       "id": 4,
       "density": 8.9,
       "density_units": "g/cm3",
-      "temperature": 0.0,
+      "temperature": 300.0,
       "elements": [
         { "element": "Ni", "percent": 0.72, "percent_type": "ao" },
         { "element": "Cr", "percent": 0.20, "percent_type": "ao" },
@@ -66,7 +66,7 @@
       "id": 5,
       "density": 1.0e-4,
       "density_units": "g/cm3",
-      "temperature": 0.0,
+      "temperature": 300.0,
       "nuclides": [
         { "name": "He4", "percent": 1.0, "percent_type": "ao" }
       ]

{processforge-0.3.4 → processforge-0.3.8}/pyproject.toml

@@ -4,7 +4,7 @@ build-backend = "setuptools.build_meta"
 
 [project]
 name = "processforge"
-version = "0.3.4"
+version = "0.3.8"
 description = "A Python-based process simulation framework for chemical engineering applications."
 readme = "README.md"
 license = "BSD-3-Clause"

{processforge-0.3.4 → processforge-0.3.8}/src/processforge/cli/__init__.py

@@ -31,7 +31,7 @@ def register_commands(app: typer.Typer) -> None:
     )(apply)
     app.command(
         "plan",
-        help="Validate a flowsheet, run DOF analysis, structural diff, and generate a Mermaid diagram",
+        help="Validate a flowsheet, run DOF analysis, structural diff, container health check, and generate a Mermaid diagram",
     )(plan)
     app.command("diagram", help="Generate a flowsheet diagram")(diagram)
     app.command(

processforge-0.3.8/src/processforge/cli/apply.py

@@ -0,0 +1,210 @@
+"""``pf apply`` — state-based warm start with homotopy fallback."""
+
+from __future__ import annotations
+
+import os
+import time
+from typing import Literal
+
+import typer
+from loguru import logger
+
+from ..eo import EOFlowsheet
+from ..provenance import build_run_info
+from ..result import save_results_zarr
+from ..state import StateManager
+from .common import (
+    build_divergence_report,
+    build_run_metadata,
+    check_providers,
+    display_backend,
+    flowsheet_basename,
+    load_state_manager,
+    log_residual_breakdown,
+    output_root,
+    require_existing_file,
+    save_snapshot,
+    validate_runtime_flowsheet,
+    write_divergence_report,
+)
+
+from .display import print_structural_diff
+
+
+def apply(
+    flowsheet: str = typer.Argument(help="Path to the flowsheet JSON file"),
+    backend: Literal["scipy", "pyomo", "casadi"] | None = typer.Option(
+        None,
+        "--backend",
+        help="Override the flowsheet's simulation.backend",
+    ),
+    tolerance: float = typer.Option(
+        1e-6,
+        "--tolerance",
+        help="Newton solver convergence tolerance (default: 1e-6)",
+    ),
+    max_iter: int = typer.Option(
+        50,
+        "--max-iter",
+        help="Max Newton iterations (default: 50)",
+    ),
+    skip_homotopy: bool = typer.Option(
+        False,
+        "--skip-homotopy",
+        help="Disable homotopy fallback; cold-start only",
+    ),
+) -> None:
+    """Apply flowsheet: drift detection, warm-start, homotopy fallback, convergence guardrails."""
+    require_existing_file(flowsheet)
+    config = validate_runtime_flowsheet(flowsheet)
+
+    # Check provider availability (assumes any containers are already running)
+    check_providers(config, flowsheet)
+
+    base_name = flowsheet_basename(flowsheet)
+
+    sim_cfg = config.get("simulation", {})
+    mode = sim_cfg.get("mode", "steady")
+    if mode != "steady":
+        logger.error("pf apply is only supported for steady-state EO flowsheets.")
+        raise SystemExit(1)
+
+    outputs_dir = output_root()
+    os.makedirs(outputs_dir, exist_ok=True)
+    sm, state = load_state_manager(outputs_dir, base_name)
+
+    # Structural diff: detect topology changes
+    topology_changed = False
+    if state is not None:
+        diff = sm.detect_structural_diff(config, state)
+        print_structural_diff(diff)
+        topology_changed = bool(diff.get("added") or diff.get("removed"))
+        if topology_changed:
+            logger.warning(
+                "Topology changed (units added/removed). "
+                "Homotopy requires identical topology — falling back to cold start."
+            )
+
+    # Parameter drift (only meaningful when topology is unchanged)
+    drifted: list[str] = []
+    current_metadata = build_run_metadata(config, tolerance, max_iter, backend or "scipy")
+    if state is not None and not topology_changed:
+        mismatches = sm.validate_metadata(current_metadata, state)
+        if mismatches:
+            logger.warning(f"Metadata mismatch: {mismatches}")
+        drifted = sm.detect_drift(config, state)
+        if not drifted:
+            logger.info("No drift detected. System is already at the desired state.")
+            return
+        stream_drifts = [d for d in drifted if d.startswith("streams.")]
+        unit_drifts = [d for d in drifted if d.startswith("units.")]
+        logger.warning("Drift detected:")
+        if stream_drifts:
+            logger.warning(f"  Stream drift : {stream_drifts}")
+        if unit_drifts:
+            logger.warning(f"  Unit drift   : {unit_drifts}")
+
+    # Build flowsheet; attach saved state for warm-start unless topology changed
+    fs = EOFlowsheet(config, backend=backend)
+    fs.saved_state = state if not topology_changed else None
+    fs.solver_tol = tolerance
+    fs.solver_max_iter = max_iter
+
+    logger.info("=== Running Apply (Steady-State EO) ===")
+    t0 = time.perf_counter()
+    results = fs.run()
+    elapsed = time.perf_counter() - t0
+
+    logger.info(
+        f"Direct solve completed in {elapsed:.2f}s (converged={fs.converged})."
+    )
+
+    if fs.converged:
+        snapshot_id = save_snapshot(
+            sm, config, fs.x_converged, fs.var_names,
+            metadata=current_metadata,
+            parent_snapshot_id=state.snapshot_id if state is not None and not topology_changed else None,
+            label="converged state",
+        )
+        run_info = build_run_info(config, x0=fs.x0, var_names=fs.var_names)
+        zarr_path = os.path.join(outputs_dir, f"{base_name}_results.zarr")
+        save_results_zarr(results, zarr_path, run_info=run_info)
+        logger.info("=== Apply Summary ===")
+        logger.info("  Status       : CONVERGED")
+        logger.info(f"  Final ||F||  : {fs.solver_stats.get('final_norm', '?'):.3e}")
+        logger.info(f"  Iterations   : {fs.solver_stats.get('iterations', '?')}")
+        logger.info(f"  Backend      : {display_backend(config, fs.backend)}")
+        logger.info(f"  Snapshot ID  : {snapshot_id}")
+        logger.info(f"  Results zarr : {zarr_path}")
+        logger.info(f"  Elapsed (s)  : {elapsed:.2f}")
+        return
+
+    # Direct solve failed — try homotopy (only when topology is same and state exists)
+    if state is not None and not topology_changed and drifted and not skip_homotopy:
+        logger.warning("Direct solve failed. Attempting homotopy continuation...")
+        from ..eo.solver import EOSolver, solve_with_homotopy
+        from ..eo.flowsheet import EOFlowsheet as _EO
+        from ..providers.manager import teardown_providers
+
+        solver = EOSolver(backend=fs.backend, tol=tolerance, max_iter=max_iter)
+        tmp_fs = _EO(config, backend=backend)
+        manager = tmp_fs._build()
+        try:
+            x_hom, hom_converged, hom_stats = solve_with_homotopy(
+                tmp_fs, manager, solver, state, drifted
+            )
+        finally:
+            teardown_providers(tmp_fs._provider_map)
+
+        if hom_converged:
+            logger.info(
+                f"Homotopy converged: ||F||={hom_stats.get('final_norm', '?'):.3e}, "
+                f"iterations={hom_stats.get('iterations', '?')}"
+            )
+            save_snapshot(
+                sm, config, x_hom, fs.var_names,
+                metadata=current_metadata,
+                parent_snapshot_id=state.snapshot_id if state is not None and not topology_changed else None,
+                label="homotopy solution",
+            )
+            run_info = build_run_info(config, x0=fs.x0, var_names=fs.var_names)
+            save_results_zarr(
+                results,
+                os.path.join(outputs_dir, f"{base_name}_results.zarr"),
+                run_info=run_info,
+            )
+            logger.info("Homotopy apply succeeded. New snapshot saved.")
+            return
+
+        # Both failed — auto-revert and write divergence report
+        logger.error("Homotopy also failed to converge.")
+        prev_id = state.snapshot_id if state is not None else "unknown"
+        sm.rollback(1)
+        logger.warning(f"Reverted .pfstate to snapshot before {prev_id}.")
+        breakdown = log_residual_breakdown(fs)
+        divergence = build_divergence_report(
+            drifted_params=drifted,
+            solver_stats=hom_stats,
+            x_last=x_hom,
+            var_names=fs.var_names,
+            breakdown=breakdown,
+        )
+    else:
+        # Cold start also failed; no rollback (nothing to revert to)
+        if not drifted and state is None:
+            logger.error("Cold-start solve failed to converge (no prior snapshot).")
+        elif skip_homotopy:
+            logger.error("Cold-start solve failed to converge (--skip-homotopy).")
+        else:
+            logger.error("Cold-start solve failed to converge.")
+        breakdown = log_residual_breakdown(fs)
+        divergence = build_divergence_report(
+            drifted_params=[],
+            solver_stats=fs.solver_stats,
+            x_last=getattr(fs, "x_converged", []),
+            var_names=getattr(fs, "var_names", []),
+            breakdown=breakdown,
+        )
+
+    write_divergence_report(outputs_dir, base_name, divergence)
+    raise SystemExit(1)

{processforge-0.3.4 → processforge-0.3.8}/src/processforge/cli/common.py

@@ -98,6 +98,23 @@ def build_run_metadata(config: dict, solver_tol: float, solver_max_iter: int, ba
     }
 
 
+def display_backend(config: dict, eo_backend: str) -> str:
+    """Resolve the backend label shown in run/apply summaries.
+
+    When a flowsheet delegates its computation to one or more providers
+    (e.g. ``openmc``, ``cantera``), the provider type is the engine that
+    actually ran the simulation — report that rather than the EO solver
+    backend (which is a no-op for provider-driven flowsheets).  Falls back
+    to the EO solver backend when no providers are declared.
+    """
+    provider_types = sorted(
+        {p.get("type") for p in config.get("providers", {}).values() if p.get("type")}
+    )
+    if provider_types:
+        return "+".join(provider_types)
+    return eo_backend
+
+
 def extract_providers(flowsheet_path: str) -> dict:
     """Read flowsheet JSON and return the raw providers dict."""
     try:
@@ -128,6 +145,21 @@ def load_state_manager(outputs_dir: str, base_name: str) -> tuple[StateManager,
 # Provider checking
 # ---------------------------------------------------------------------------
 
+def is_local_provider_url(url: str | None) -> bool:
+    """Return True if *url* points at a locally-managed provider.
+
+    A provider is considered local — and thus a candidate for ``pf init`` to
+    generate docker-compose and pull its image — when it has no explicit URL
+    or its URL points at localhost / 127.0.0.1 / 0.0.0.0. Any other URL is
+    treated as an externally-managed (cloud / remote) provider that
+    ``pf init`` must not try to spin up or pull an image for.
+    """
+    if not url:
+        return True
+    lowered = url.lower()
+    return any(token in lowered for token in ("localhost", "127.0.0.1", "0.0.0.0"))
+
+
 def _resolve_provider_url(cfg: dict, ptype: str) -> str:
     """Build the provider URL, falling back to the default port."""
     from ..providers.registry import get_provider_default_port
@@ -139,6 +171,19 @@ def _resolve_provider_url(cfg: dict, ptype: str) -> str:
     return url
 
 
+def _ping_provider_health(url: str, timeout: int = 5) -> tuple[bool, "dict | str"]:
+    """GET ``{url}/health`` and return ``(ok, payload_or_error)``.
+
+    Returns the parsed JSON health payload on success, or the error string on
+    failure (unreachable service, timeout, or non-JSON response).
+    """
+    try:
+        with urllib.request.urlopen(f"{url.rstrip('/')}/health", timeout=timeout) as resp:
+            return True, json.loads(resp.read().decode())
+    except (urllib.error.URLError, OSError, TimeoutError, ValueError) as exc:
+        return False, str(exc)
+
+
 def check_providers(
     config: dict,
     flowsheet_path: str,
@@ -203,50 +248,6 @@ def check_providers(
         raise SystemExit(1)
 
 
-def has_containerized_providers(config: dict) -> bool:
-    """Return True if any declared provider runs in a Docker container."""
-    from ..providers.registry import is_containerized
-
-    return any(
-        is_containerized(cfg.get("type", ""))
-        for cfg in config.get("providers", {}).values()
-    )
-
-
-def wait_for_provider_health(config: dict, timeout: int = 30) -> None:
-    """Poll ``/health`` for each localhost containerized provider until ready.
-
-    Only providers whose resolved URL host is ``localhost``/``127.0.0.1`` are
-    waited on — remote providers are assumed to be managed elsewhere.
-    """
-    import time
-    import urllib.parse
-
-    from ..providers.registry import is_containerized
-
-    providers = config.get("providers", {})
-    deadline = time.time() + timeout
-    for name, cfg in providers.items():
-        if not is_containerized(cfg.get("type", "")):
-            continue
-        url = _resolve_provider_url(cfg, cfg.get("type", ""))
-        host = urllib.parse.urlparse(url).hostname
-        if host not in ("localhost", "127.0.0.1"):
-            continue
-        health = f"{url}/health"
-        while time.time() < deadline:
-            try:
-                urllib.request.urlopen(health, timeout=2)
-                logger.info(f"Provider '{name}' ready at {url}")
-                break
-            except (urllib.error.URLError, OSError):
-                time.sleep(1)
-        else:
-            logger.warning(
-                f"Provider '{name}' not ready at {url} after {timeout}s."
-            )
-
-
 # ---------------------------------------------------------------------------
 # Divergence reports
 # ---------------------------------------------------------------------------

{processforge-0.3.4 → processforge-0.3.8}/src/processforge/cli/display.py

@@ -53,6 +53,61 @@ def print_unit_mismatches(mismatches: list[Any]) -> None:
             logger.warning(f"Unit annotation — stream '{m.stream_name}'.{m.property_name}: {m.message}")
 
 
+def print_provider_health(config: dict) -> list[str]:
+    """Log per-provider reachability for flowsheets using Docker containers.
+
+    For each provider declared in *config*:
+
+    * containerized providers are probed via ``GET /health`` and reported with
+      their URL, ``status``, and ``provider_type`` from the health payload;
+    * pip-installable providers are reported as importable or not.
+
+    Returns a list of failure messages — one per unreachable containerized
+    provider (empty means every provider is healthy/available).
+    """
+    import importlib
+
+    from ..providers.registry import is_containerized, _PROVIDER_CATALOG
+
+    providers = config.get("providers", {})
+    logger.info("=== Provider / Container Health ===")
+    failures: list[str] = []
+    if not providers:
+        logger.info("  No providers declared.")
+        return failures
+
+    for name, cfg in providers.items():
+        ptype = cfg.get("type", "")
+        if is_containerized(ptype):
+            from .common import _resolve_provider_url, _ping_provider_health
+
+            url = _resolve_provider_url(cfg, ptype)
+            ok, info = _ping_provider_health(url, timeout=5)
+            if ok:
+                payload = info if isinstance(info, dict) else {}
+                status = payload.get("status", "?")
+                provider_type = payload.get("provider_type", "?")
+                logger.info(
+                    f"  [OK] {name} [{ptype}] {url} — status={status} provider_type={provider_type}"
+                )
+            else:
+                msg = f"  [ERR] {name} [{ptype}] {url} — unreachable: {info}"
+                logger.error(msg)
+                failures.append(f"Provider '{name}' unreachable at {url}: {info}")
+        else:
+            catalog = _PROVIDER_CATALOG.get(ptype, {})
+            module = catalog.get("module", "")
+            try:
+                importlib.util.find_spec(module)
+                logger.info(f"  [OK] {name} [{ptype}] (pip — importable)")
+            except (ModuleNotFoundError, ValueError):
+                dep = catalog.get("optional_dep")
+                hint = f"pip install 'processforge[{dep}]'" if dep else "built-in"
+                logger.warning(f"  [WARN] {name} [{ptype}] — not installed ({hint})")
+
+    return failures
+
+
 def print_structural_diff(diff: dict) -> None:
     """Print a +/~/- structural diff of units."""
     logger.info("=== Structural Diff vs. Saved State ===")

{processforge-0.3.4 → processforge-0.3.8}/src/processforge/cli/init.py

@@ -11,7 +11,56 @@ import subprocess
 import typer
 from loguru import logger
 
-from .common import extract_providers
+from .common import extract_providers, is_local_provider_url
+from ..lock import flowsheet_env_dir, read_lock
+
+
+def _migrate_legacy_env(pf_dir: str) -> None:
+    """Move a legacy root-level lock.json / docker-compose.yml into a per-flowsheet dir.
+
+    Older processforge stored a single environment at ``.processforge/lock.json``
+    and ``.processforge/docker-compose.yml``. This relocates those into the hashed
+    env dir derived from the flowsheet recorded in the legacy lock, so an existing
+    repo keeps its provider environment after upgrading.
+    """
+    legacy_lock = os.path.join(pf_dir, "lock.json")
+    legacy_compose = os.path.join(pf_dir, "docker-compose.yml")
+    if not (os.path.exists(legacy_lock) or os.path.exists(legacy_compose)):
+        return
+
+    recorded = None
+    if os.path.exists(legacy_lock):
+        try:
+            recorded = read_lock(pf_dir).get("flowsheet")
+        except Exception:
+            recorded = None
+
+    if recorded and os.path.exists(recorded):
+        env_dir = flowsheet_env_dir(pf_dir, recorded)
+    else:
+        base = os.path.splitext(os.path.basename(recorded or "legacy"))[0] or "legacy"
+        env_dir = os.path.join(pf_dir, base)
+        logger.warning(
+            "Migrating legacy .processforge/lock.json + docker-compose.yml into "
+            f"'{os.path.relpath(env_dir, pf_dir)}/' (recorded flowsheet not found on "
+            "disk — re-run `pf init <flowsheet.json>` to restore the correct env dir)."
+        )
+
+    if os.path.exists(env_dir):
+        logger.warning(
+            f"Migration target {os.path.relpath(env_dir, pf_dir)}/ already exists — "
+            "skipping legacy migration."
+        )
+        return
+
+    os.makedirs(env_dir, exist_ok=True)
+    for src in (legacy_lock, legacy_compose):
+        if os.path.exists(src):
+            shutil.move(src, os.path.join(env_dir, os.path.basename(src)))
+    logger.info(
+        "Migrated legacy .processforge/ environment into per-flowsheet dir "
+        f"'{os.path.relpath(env_dir, pf_dir)}/'."
+    )
 
 
 def init(
@@ -42,6 +91,11 @@ def init(
     os.makedirs(pf_dir, exist_ok=True)
     os.makedirs(outputs_dir, exist_ok=True)
 
+    # Migrate a legacy single-environment layout (.processforge/lock.json and
+    # .processforge/docker-compose.yml at the root) into a per-flowsheet env
+    # dir so existing repos aren't silently broken by the new structure.
+    _migrate_legacy_env(pf_dir)
+
     # Remove stale .pfstate snapshot directories from outputs/
     stale_count = 0
     for entry in os.listdir(outputs_dir):
@@ -84,21 +138,32 @@ def init(
 
     # Categorize providers
     pip_providers: dict[str, dict] = {}
-    docker_providers: dict[str, dict] = {}
+    local_docker_providers: dict[str, dict] = {}
+    remote_docker_providers: dict[str, dict] = {}
     for name, cfg in providers.items():
         ptype = cfg.get("type", "")
         if is_containerized(ptype):
             url = cfg.get("url")
-            if not url:
+            if is_local_provider_url(url):
                 port = get_provider_default_port(ptype) or 9000
-                url = f"http://localhost:{port}"
-            docker_providers[name] = {
-                "type": ptype,
-                "url": url,
-                "docker_image": cfg.get("docker_image") or get_provider_docker_image(ptype),
-                "port": get_provider_default_port(ptype),
-            }
-            logger.info(f"  {name}: type={ptype}, url={url} (Docker)")
+                if not url:
+                    url = f"http://localhost:{port}"
+                local_docker_providers[name] = {
+                    "type": ptype,
+                    "url": url,
+                    "docker_image": cfg.get("docker_image") or get_provider_docker_image(ptype),
+                    "port": port,
+                }
+                logger.info(f"  {name}: type={ptype}, url={url} (Docker, local)")
+            else:
+                remote_docker_providers[name] = {
+                    "type": ptype,
+                    "url": url,
+                    "docker_image": cfg.get("docker_image") or get_provider_docker_image(ptype),
+                }
+                logger.info(
+                    f"  {name}: type={ptype}, url={url} (Docker, remote — skipping compose)"
+                )
         else:
             pip_providers[name] = {"type": ptype}
             logger.info(f"  {name}: type={ptype} (pip)")
@@ -116,15 +181,19 @@ def init(
             hint = f"pip install 'processforge[{dep}]'" if dep else "built-in"
             logger.warning(f"  {name} — not installed. Install with: {hint}")
 
-    # Generate compose for Docker providers
-    if docker_providers:
-        compose_path = os.path.join(pf_dir, "docker-compose.yml")
+    # Generate compose + pull images only for locally-managed Docker providers.
+    # Providers with an explicit remote URL are assumed to be running elsewhere
+    # (e.g. a cloud deployment of the ghcr.io image) and are not touched here.
+    if local_docker_providers:
+        compose_path = os.path.join(
+            flowsheet_env_dir(pf_dir, flowsheet_path), "docker-compose.yml"
+        )
         if os.path.exists(compose_path):
             logger.warning(
                 f"Environment already initialized — reinitializing from {flowsheet_path}"
             )
 
-        generate_compose(pf_dir, docker_providers, outputs_dir)
+        generate_compose(pf_dir, local_docker_providers, outputs_dir, flowsheet=flowsheet_path)
         logger.info(f"Generated {compose_path}")
 
         # Attempt docker compose pull — stream progress live
@@ -155,6 +224,10 @@ def init(
             )
         except subprocess.TimeoutExpired:
             logger.warning("docker compose pull timed out after 600s.")
+    elif remote_docker_providers:
+        logger.info(
+            "All containerized providers use remote URLs — skipping Docker setup."
+        )
     else:
         logger.info("No containerized providers — skipping Docker setup.")
 
@@ -177,5 +250,5 @@ def init(
     from .. import __version__ as pf_version
 
     write_lock(pf_dir, flowsheet_path, lock_providers, pf_version)
-    logger.info(f"Wrote {os.path.join(pf_dir, 'lock.json')}")
+    logger.info(f"Wrote {os.path.join(flowsheet_env_dir(pf_dir, flowsheet_path), 'lock.json')}")
     logger.info(".processforge/ initialised successfully.")

{processforge-0.3.4 → processforge-0.3.8}/src/processforge/cli/plan.py

@@ -18,7 +18,12 @@ from .common import (
     require_existing_file,
     validate_snapshot_config,
 )
-from .display import print_dof_report, print_structural_diff, print_unit_mismatches
+from .display import (
+    print_dof_report,
+    print_provider_health,
+    print_structural_diff,
+    print_unit_mismatches,
+)
 
 
 def plan(
@@ -116,7 +121,10 @@ def plan(
         logger.info("  Warm-start available : No snapshot found")
         logger.info("  Homotopy eligible    : No (cold start)")
 
-    # Step 7: Mermaid diagram
+    # Step 7: Provider / container health check
+    health_failures = print_provider_health(config)
+
+    # Step 8: Mermaid diagram
     if not no_diagram:
         diagram_output_dir = output_dir or "diagrams"
         os.makedirs(diagram_output_dir, exist_ok=True)
@@ -129,5 +137,5 @@ def plan(
 
     # Exit non-zero on hard errors
     hard_errors = [m for m in mismatches if not m.compatible]
-    if hard_errors or report.system_dof < 0:
+    if hard_errors or report.system_dof < 0 or health_failures:
         raise SystemExit(1)