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

Aperta
#5,486 0 commenti 0 reazioni 0 assegnatari Vedi su GitHub

I maintainer di solito rispondono entro 1 giorno

Nessuno ha ancora preso questa issue.

Valutazione

Difficoltà
4/5
Tempo stimato
3-5 giorni
Idoneità per principianti
48/100
Tipo di issue
Bug
Chiarezza
Abbastanza chiara
Stato di attività
Attiva
Stack tecnologico
python
Ambito
hpc

Direzione di ricerca

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.

Scritto dal modello di indicizzazione a partire dal testo della issue.

Descrizione

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

Lingua principale
Python
Stelle
682
Fork
204
Merge medio
4g 14h
PR unite (30g)
46

Preparare l'ambiente

  • Nessun Dockerfile né file Docker Compose
  • Ha un modello di pull request
  • Nessuna guida per i contributori

Come iniziare

  1. Leggi tutta la issue e poi la guida ai contributi del progetto.
  2. Commenta sulla issue per dire che te ne occupi tu — evita che due persone facciano lo stesso lavoro.
  3. Fai un fork del repository e lavora su un branch.
  4. Apri una pull request che faccia riferimento al numero della issue.

Altre issue di firedrakeproject/firedrake

Tutte le issue di firedrakeproject/firedrake

Issue simili

Altre issue su Python

Ricevi le nuove issue nella tua casella

Un breve riepilogo di issue GitHub adatte ai principianti.