EQUIL - NEW FEATURE - Pin located kinetic-resonance surfaces into the auto psi grid - #422
Draft
logan-nc wants to merge 2 commits into
Draft
EQUIL - NEW FEATURE - Pin located kinetic-resonance surfaces into the auto psi grid#422logan-nc wants to merge 2 commits into
logan-nc wants to merge 2 commits into
Conversation
… auto psi grid The two-pass auto grid's criterion is ideal-driven and knows nothing about kinetic resonance locations. When a run builds calculated kinetic matrices, locate the Omega_l = 0 surfaces (same locator as the NTV quadrature paneling) and insert them as plain knots via merge_mandatory_nodes -- knot-at-node, no cleared zone, inserted before rational bracketing so the Delta-prime clean-interval treatment wins locally. Nodes inside the near-axis validity region are suppressed anyway and not pinned. DIII-D: +2 net knots, et[1] unchanged to 2e-6, EL steps drop 14%. Co-Authored-By: Claude Fable 5 <noreply@anthropic.com> Claude-Session: https://claude.ai/code/session_01LzbLFQKyuRE5DYZmLokKmk
…onstructions Co-Authored-By: Claude Fable 5 <noreply@anthropic.com> Claude-Session: https://claude.ai/code/session_01LzbLFQKyuRE5DYZmLokKmk
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Add this suggestion to a batch that can be applied as a single commit.This suggestion is invalid because no changes were made to the code.Suggestions cannot be applied while the pull request is closed.Suggestions cannot be applied while viewing a subset of changes.Only one suggestion per line can be applied in a batch.Add this suggestion to a batch that can be applied as a single commit.Applying suggestions on deleted lines is not supported.You must change the existing code in this line in order to create a valid suggestion.Outdated suggestions cannot be applied.This suggestion has been applied or marked resolved.Suggestions cannot be applied from pending reviews.Suggestions cannot be applied on multi-line comments.Suggestions cannot be applied while the pull request is queued to merge.Suggestion cannot be applied right now. Please check back later.
Stacked on #414. The kinetic-grid follow-up chosen after the certified-grid postmortem: instead of adaptive compression, make the auto grid resonance-aware statically — the cheapest layer of the safety-net design, and per the tests below, the only one currently justified.
What
When a run builds calculated kinetic matrices, locate the Ω_ℓ = 0 kinetic-resonance surfaces (same
kinetic_resonance_psi_nodeslocator the NTV quadrature panels at — one source of truth) and pin them into the two-pass auto ψ grid as plain knots viamerge_mandatory_nodes— knot-at-node, no cleared zone, inserted before rational bracketing so the Δ′ clean-interval treatment wins wherever a resonance sits inside a rational's bracket. Nodes inside the near-axis validity region (ψ < ψ_c, #414) are suppressed anyway and are not pinned. ~30 lines, no new knobs.Evidence
toroidal_rotation_factor = 0.2parks Ω_ℓ=0 nodes at ψ = 0.578, 0.812): naive auto grid vs 1024-knot gold agrees to 1.24e-3 — the profile-curvature term already tracks the reshaped ω_E — and pinning reproduces it to 1e-6. The feared naive-grid physics risk is bounded at the convergence-scatter level for this machine; pinning is cheap insurance plus the step-count win.Design provenance
This replaces the parked adaptive approach (#410): detection of possible resonances is essentially free from the frequency profiles, so the grid can be made resonance-aware without any kernel-driven refinement loop. The remaining safety-net layers (a free post-hoc residual detector; one bounded probe-and-insert round) are designed and recorded in
handoff/issue376/RESULTS.md§31 but not implemented — every stress test so far says the seeded auto grid is sufficient.🤖 Generated with Claude Code
https://claude.ai/code/session_01LzbLFQKyuRE5DYZmLokKmk