Feature request: vectorized finite difference evaluation
評価
- 難易度
- 5/5
- 見積もり時間
- 1週間以上
- 初心者へのやさしさ
- 45/100
- issue の種類
- 機能追加
- 明瞭さ
- おおむね明確
- 活発さ
- 静か
- 技術スタック
- julia
- 領域
- performance
調査の方向性
まず、既存の FiniteDiff および AutoFiniteDiff API と、issue にある手動でベクトル化された GPU のユースケースを読みます。dims 引数で行単位の有限差分接ベクトルをどのように表現するか、また出力形状をヤコビアンにどのように対応付けるべきかを判断します。一般化された API でベクトル化評価を 1 回の関数呼び出しで実行できれば完了です。
索引モデルが issue の本文から書いたものです。
説明
I have an implicit function that is extremely expensive to evaluate. However, it is fully GPU parallelized, in that if a matrix is provided, then each row is evaluated in parallel. For this particular function, the output is a matrix of size = (something, number of rows)
Doing finite differences through this is then easily parallelizable by making each row be a finite difference "tangent" (sorry if I'm misusing that word), so that the full Jacobian can be constructed in one single function evaluation as opposed to many evaluations.
I am doing this now manually, but it would be nice if this functionality could be generalized into FiniteDiff. perhaps a dims argument could be provided to AutoFiniteDiff specifying the dimension about which the evaluation is vectorized (in my case across the rows, so dims=1).
- 主要言語
- Julia
- スター
- 305
- フォーク
- 42
- 平均マージ
- 24分
- マージ済み PR(30日)
- 2
環境構築
このプロジェクトには開発コンテナ、Dockerfile、コントリビューションガイドがありません。まず README を読み、一般的な手順ははじめてのコントリビューションガイドを参照してください。
はじめの一歩
- issue を最後まで読み、次にプロジェクトのコントリビューションガイドを読みます。
- 着手することを issue にコメントします — 二人が同じ作業をするのを防げます。
- リポジトリをフォークし、ブランチを切って変更します。
- issue 番号を参照したプルリクエストを送ります。
JuliaDiff/FiniteDiff.jl のほかの issue
-
tagbot-manual
難易度 2/5 1〜3時間 初心者へのやさしさ 55/100
JuliaDiff/FiniteDiff.jl#218 ·
-
難易度 3/5 1〜2日 初心者へのやさしさ 45/100
JuliaDiff/FiniteDiff.jl#199 · コメント 3 件 ·
-
StridedArraysオープン
難易度 4/5 3〜5日 初心者へのやさしさ 25/100
JuliaDiff/FiniteDiff.jl#158 · コメント 1 件 ·
-
難易度 5/5 1週間以上 初心者へのやさしさ 25/100
JuliaDiff/FiniteDiff.jl#157 · コメント 1 件 ·
-
難易度 5/5 1週間以上 初心者へのやさしさ 20/100
JuliaDiff/FiniteDiff.jl#150 · コメント 1 件 ·
JuliaDiff/FiniteDiff.jl の issue をすべて見る
似ている issue
-
難易度 2/5 1〜3時間 初心者へのやさしさ 72/100
oscar-system/GAP.jl#1455 ·
メンテナーはふだん 1 日以内に返信
-
難易度 2/5 1〜3時間 初心者へのやさしさ 65/100
JuliaAI/MLJTransforms.jl#85 ·
-
難易度 2/5 1〜3時間 初心者へのやさしさ 70/100
SciML/OrdinaryDiffEq.jl#4875 ·
メンテナーはふだん 3 日以内に返信
-
documentation
難易度 2/5 1〜3時間 初心者へのやさしさ 62/100
メンテナーはふだん 1 日以内に返信
-
難易度 2/5 1〜3時間 初心者へのやさしさ 66/100
SciML/NonlinearSolve.jl#1372 ·
メンテナーはふだん 1 日以内に返信