Tearing.Dispersion - Scale the polish's finite-difference step to |Q|

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

维护者通常 3 天内回复

还没有人认领这个 Issue。

评估

难度
3/5
预计耗时
1-2 天
新手友好度
54/100
Issue 类型
缺陷
描述清晰度
基本清楚
活跃度
活跃
技术栈
julia
领域
performance

调研方向

Find _polish_root and inspect how its finite-difference step h is used to compute the Jacobian near small trust radii and large |Q|. Run the harness before and after any proposed step scaling, then compare and report changes in accepted roots.

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

描述

_polish_root uses h = max(1e-3 R, eps(Float64)), where R is the trust radius. The step does not scale with |Q|.

  • When R is small (near-coincident intersections), fx - fk loses digits to cancellation.
  • If h ever falls below about eps·|Q|, Q + h == Q, the Jacobian determinant is zero, and the polish silently does nothing.

Proposal: h = clamp(1e-3 R, sqrt(eps)·max(|Q|, 1), R/2), or an equivalent relative step. This changes the polish numerics slightly, so it needs a harness run, and any change in accepted roots should be reported.

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

环境准备

  • 没有 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 摘要。