Equilibrium - Make the magnetic-axis search's singular-Jacobian threshold relative

未关闭
#608 0 条评论 0 个 reaction 已指派 0 人 在 GitHub 查看

维护者通常 3 天内回复

还没有人认领这个 Issue。

评估

难度
5/5
预计耗时
一周以上
新手友好度
38/100
Issue 类型
功能
描述清晰度
基本清楚
活跃度
活跃
技术栈
julia

调研方向

Start with Equilibrium/Axisymmetric/Direct.jl to inspect the Newton search and its singular_det check. Then check FastInterpolations' implementation or documentation for whether 2-D cubic_interp is C² or C¹. Done means the singular-Jacobian threshold scales with the equilibrium, and the fallback accounts for interpolation behavior that can make Newton cycle.

由索引模型根据 Issue 内容生成。

描述

The Newton search for the magnetic axis (Equilibrium/Axisymmetric/Direct.jl) treats the Hessian of ψ as singular below singular_det = 1e-20, an absolute threshold in the units of ψ_RR·ψ_ZZ. Its meaning therefore depends on the device size and the field strength.

Proposal: compare the determinant with a scale built from the equilibrium, for example (ψ_axis − ψ_boundary)²/L⁴ with L the minor radius, so the same relative threshold applies to a compact device and to an ITER-sized one. Also confirm whether FastInterpolations' 2-D cubic_interp is C² or only C¹. If only C¹, the Hessian jumps at grid nodes and an axis near a node can make Newton cycle, which the fallback should then expect.

主要语言
Julia
星标
11
派生
7
平均合并
3 天 11 小时
30 天内合并 PR
68

环境准备

  • 没有 Dockerfile 或 Docker Compose 文件
  • 有 Pull Request 模板
  • 没有贡献指南

从这里开始

  1. 先读完整个 Issue,再读项目的贡献指南。
  2. 在 Issue 下留言说明你要接手 —— 这能避免两个人做同样的事。
  3. Fork 仓库,在一个分支上完成修改。
  4. 提交 Pull Request,并在描述里引用这个 Issue 编号。

OpenFUSIONToolkit/GPEC 的其他 Issue

查看 OpenFUSIONToolkit/GPEC 的全部 Issue

相似的 Issue

更多 Julia Issue

把新 issue 发到你的邮箱

精选适合新手参与的 GitHub issue 摘要。