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

オープン
#605 コメント 0 件 リアクション 0 件 担当者 0 名 GitHub で見る

メンテナーはふだん 3 日以内に返信

まだ誰も着手していません。

評価

難易度
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
平均マージ
4日 3時間
マージ済み PR(30日)
70

環境構築

  • Dockerfile・Docker Compose ファイルなし
  • プルリクエストのテンプレートあり
  • コントリビューションガイドなし

はじめの一歩

  1. issue を最後まで読み、次にプロジェクトのコントリビューションガイドを読みます。
  2. 着手することを issue にコメントします — 二人が同じ作業をするのを防げます。
  3. リポジトリをフォークし、ブランチを切って変更します。
  4. issue 番号を参照したプルリクエストを送ります。

OpenFUSIONToolkit/GPEC のほかの issue

OpenFUSIONToolkit/GPEC の issue をすべて見る

似ている issue

Julia の issue をもっと見る

新しい issue をメールで受け取る

初心者向けの GitHub issue を短くまとめたダイジェスト。