perturbation-sim: eigenvalue-perturbation + edge-propagation outage simulator

Cargo.toml

@@ -39,6 +39,12 @@ exclude = [
     "crates/learning",
     "crates/jc",
     "crates/sigker",
+    # Outage perturbation-shape simulator — standalone zero-dep research crate.
+    # The applied companion to jc: supplies the genuine eigenvalue-perturbation
+    # result (Weyl/Davis-Kahan on the weighted Laplacian) that jc::weyl
+    # (equidistribution) is NOT, composed with a DC-power-flow/LODF cascade.
+    # Verify via `cargo test --manifest-path crates/perturbation-sim/Cargo.toml`.
+    "crates/perturbation-sim",
     # Place/Residue golden-spiral codec — standalone, with a MANDATORY git dep on
     # the AdaWorldAPI ndarray fork (the SIMD foundation). Verified via
     # `cargo test --manifest-path crates/helix/Cargo.toml`.

crates/perturbation-sim/Cargo.toml

@@ -0,0 +1,62 @@
+[package]
+name = "perturbation-sim"
+version = "0.1.0"
+edition = "2021"
+license = "Apache-2.0"
+description = "Spectral + edge-propagation outage simulator: models the perturbation SHAPE of a cascading grid failure. Weyl/Davis-Kahan eigenvalue-perturbation of the weighted graph Laplacian (a line trip = low-rank perturbation E) composed with a DC-power-flow / LODF cascade. Standalone, zero-dep, deterministic — the applied companion to lance-graph crates/jc (jc::weyl = equidistribution; this crate is the missing eigenvalue-perturbation pillar) and jc::ewa_sandwich (edge propagation)."
+
+# Standalone research crate (excluded from the lance-graph workspace, like jc /
+# sigker / helix). Verify with:
+#   cargo test  --manifest-path crates/perturbation-sim/Cargo.toml
+#   cargo run   --manifest-path crates/perturbation-sim/Cargo.toml --example simulate
+[features]
+# Route the Morton/Walsh pyramid transform through ndarray's SIMD
+# (AVX-512 under x86-64-v4) Walsh–Hadamard. AMX is int8/bf16 tile-GEMM and is
+# NOT used (this crate's WHT is f32, its field tier f64). Default OFF → stays
+# zero-dep with the scalar `fwht` fallback. Build accelerated with:
+#   RUSTFLAGS='-C target-cpu=x86-64-v4' cargo … --features ndarray-simd
+# (or target-cpu=native locally). All SIMD comes from `ndarray::simd` per the
+# workspace rule — never raw intrinsics here.
+ndarray-simd = ["dep:ndarray"]
+
+[dependencies]
+# The AdaWorldAPI fork ("The Foundation"), optional + behind `ndarray-simd`.
+# Sourced by GIT (not a local path): an optional *path* dep is read at manifest
+# resolution, so a clean checkout without the sibling `../../../ndarray` fails
+# even with the feature OFF — breaking the zero-dep default (codex P2 on #504,
+# the helix-#460 lesson). A git source resolves remotely and is only fetched/
+# built when `ndarray-simd` activates it. Local-dev alternative: swap to
+#   ndarray = { path = "../../../ndarray", optional = true, default-features = false, features = ["std"] }
+ndarray = { git = "https://github.com/AdaWorldAPI/ndarray.git", branch = "master", optional = true, default-features = false, features = ["std"] }
+
+[lib]
+name = "perturbation_sim"
+path = "src/lib.rs"
+
+[[example]]
+name = "simulate"
+path = "examples/simulate.rs"
+
+[[example]]
+name = "iberian"
+path = "examples/iberian.rs"
+
+[[example]]
+name = "validate"
+path = "examples/validate.rs"
+
+[[example]]
+name = "weakest_links"
+path = "examples/weakest_links.rs"
+
+[[example]]
+name = "hhtl_levels"
+path = "examples/hhtl_levels.rs"
+
+[[example]]
+name = "reinforce"
+path = "examples/reinforce.rs"
+
+[[example]]
+name = "hhtl_resident"
+path = "examples/hhtl_resident.rs"