Tearing.Dispersion - Update the root polish's Jacobian by Broyden instead of two differences per step

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

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

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

評価

難易度
4/5
見積もり時間
3〜5日
初心者へのやさしさ
52/100
issue の種類
リファクタリング
明瞭さ
明確に書かれている
活発さ
活発
技術スタック
julia
領域
performance

調査の方向性

Read _polish_root in Tearing/Dispersion/GrowthRateExtraction.jl and trace how it evaluates the residual, accepts trial steps, and handles up to three backtracks. Update the Jacobian after accepted steps using the step and residual change, and rebuild it from differences after a failed backtrack. Run the harness and compare polished roots with the existing tolerance.

索引モデルが issue の本文から書いたものです。

説明

The root polish (_polish_root in Tearing/Dispersion/GrowthRateExtraction.jl) builds the real 2×2 Jacobian of (Re f, Im f) from two forward differences on every iteration. With the trial step and up to three backtracks that is 3–6 residual evaluations per iteration, and on the SLAYER path each evaluation is a stiff Riccati solve.

The cost is small today (diiid_slayer_n1 127 → 128 s), but it scales with the number of intersections polished and with anything that re-polishes, such as more surfaces or a finer scan.

Proposal: build the Jacobian once from two differences, then update it by Broyden's rank-one rule with the step and residual change already computed, so each iteration costs one evaluation. Rebuild it from differences only after a failed backtrack. Convergence and stopping criteria stay as they are, so roots should agree with the current polish to its tolerance; check with the harness.

主要言語
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 を短くまとめたダイジェスト。