Does `PCDPC` actually implement DCP rather than PCD??

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

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

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

評価

難易度
4/5
見積もり時間
3〜5日
初心者へのやさしさ
48/100
issue の種類
バグ
明瞭さ
おおむね明確
活発さ
活発
技術スタック
python
領域
hpc

調査の方向性

Read firedrake/preconditioners/pcd.py, especially the docstring near lines 13–22 and the implementation near lines 124–128. Compare both operator orderings with equation 4.13 in the cited Elman et al. reference, then determine whether the code, documentation, or both are inconsistent. Done means the PCD definition and its implementation agree with the reference.

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

説明

bug straightforward

The PCD preconditioner implements the pressure-convection-diffusion preconditioner from Elman et al 2006

$$
S \approx K F^{-1} M
$$

Where $K$ is a diffusion operator, $F$ is an advection-diffusion operator on the pressure space, and $M$ is the pressure mass matrix.

The preconditioner should apply the inverse (equation 4.13 in the reference above):

$$
S^{-1} \approx M^{-1} F K^{-1}
$$

However, both the docstring and the code write the inverse as:

$$
S^{-1} \approx K^{-1} F M^{-1}
$$

https://github.com/firedrakeproject/firedrake/blob/2bfed39d4dfa617dd677ceafe2a880d7057c0d69/firedrake/preconditioners/pcd.py#L13-L22

https://github.com/firedrakeproject/firedrake/blob/2bfed39d4dfa617dd677ceafe2a880d7057c0d69/firedrake/preconditioners/pcd.py#L124-L128

主要言語
Python
スター
684
フォーク
204
平均マージ
4日 14時間
マージ済み PR(30日)
46

環境構築

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

はじめの一歩

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

firedrakeproject/firedrake のほかの issue

firedrakeproject/firedrake の issue をすべて見る

似ている issue

Python の issue をもっと見る

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

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