PerturbedEquilibrium - BUGFIX! - Reject a perturbed equilibrium solved with a wall

.claude/agents/regression-guardian.md

@@ -36,8 +36,8 @@ Useful invocations:
   cases are added over time, so never assume the list. (At time of writing: `diiid_n1`,
   `diiid_n1_riccati`, `diiid_slayer_n1`, `efit_fixedbdy_separatrix`, `gal_resistive_diiid`,
   `gal_resistive_pe`, `ggj_ray_q500i`, `ggj_reference`, `solovev_n1`, `solovev_multi_n`,
-  `solovev_kinetic_calculated`, `solovev_kinetic_multiion`, `solovev_kinetic_ntv`,
-  `solovev_kinetic_nuzero`.)
+  `solovev_kinetic_nuzero`, `tj_n1`, `tj_kinetic_calculated`, `tj_kinetic_multiion`,
+  `tj_kinetic_ntv`.)
 - `--cases <case[,case]> --refs <refA,refB>` — compare two refs (branches/commits), e.g.
   `--refs develop,my-feature-branch`. `local` names the uncommitted working tree; prefer a branch
   ref (see "Which refs to compare" below).
@@ -49,7 +49,7 @@ Useful invocations:
 Pick the case(s) by the module touched:
 - Equilibrium / ideal stability / perturbed equilibrium → `diiid_n1` (broadest, ~31 quantities).
 - Analytical equilibrium + ideal stability → `solovev_n1`, `solovev_multi_n`.
-- KineticForces (NTV) changes → `solovev_kinetic_calculated`.
+- KineticForces (NTV) changes → `tj_kinetic_calculated`.
 - InnerLayer (GGJ resistive layer) changes → `ggj_reference`.
 
 When in doubt, run `diiid_n1` plus whichever case targets the module you changed.

examples/Solovev_ideal_example/analyze_example.jl

@@ -9,7 +9,6 @@ h5path = joinpath(@__DIR__, "gpec.h5")
 # Summary plots
 p_eq = Analysis.Equilibrium.plot_equilibrium_summary(h5path)
 p_ffs = Analysis.ForceFreeStates.plot_ffs_summary(h5path)
-p_pe  = Analysis.PerturbedEquilibrium.plot_perturbed_equilibrium_summary(h5path)
 
 display(p_eq);  Plots.savefig(p_eq,  joinpath(@__DIR__, "equilibrium_summary.png"))
 display(p_ffs); Plots.savefig(p_ffs, joinpath(@__DIR__, "ffs_summary.png"))

examples/Solovev_ideal_example/gpec.toml

@@ -1,4 +1,4 @@
-# Solovev analytical equilibrium — n=1 ideal stability + perturbed equilibrium.
+# Solovev analytical equilibrium — n=1 ideal stability.
 # Near-marginal n=1 external kink stabilized by a close conformal wall (see [Wall] note).
 # Integration-pipeline smoke-test and et[1] regression fixture; the equilibrium is
 # generated analytically from the embedded [Solovev] section.
@@ -27,14 +27,6 @@ shape = "conformal"                     # Wall shape (conformal, elliptical, dee
 standoff_fraction = 0.2415              # Wall standoff in plasma minor radii
 equal_arc_wall = true                   # Equal arc length distribution of nodes on wall
 
-[ForcingTerms]
-forcing_data_file = "forcing.dat"       # Manual mode table (n, m, amplitude_real, amplitude_imag)
-forcing_data_format = "ascii"           # Format: "ascii", "hdf5", or "coil" (Biot-Savart from 3D wires)
-
-[PerturbedEquilibrium]
-compute_singular_coupling = true        # Compute singular layer coupling metrics
-regularization_width = 0.0              # Regularization width for singular surfaces (0 = disabled)
-
 # Note: SLAYER tearing analysis is not run on the Solovev analytic equilibrium
 # — its Δ' / inner-layer dispersion does not yield meaningful tearing roots.
 # The SLAYER example and regression case use the DIII-D-like geqdsk instead