jorgensd / jorgensd/dolfinx-tutorial
Rethink nullspace tutorials
Nessuno ha ancora preso questa issue.
- Lingua principale
- Jupyter Notebook
- Stelle
- 167
- Fork
- 89
- Merge medio
- 3h 25m
- PR unite (30g)
- 2
Descrizione
Direct solver
PETSc does magic after the fact with nullspace information, which is "correct", but mumps might fail if
icntl_24:1 and inctl_25: 0 is not set for more complex spaces.
Iterative solver
Hypre/Boomeramg ignores the near nullspace if
"pc_hypre_boomeramg_nodal_coarsen": (1-6)
"pc_hypre_boomeramg_vec_interp_variant": (1-3)
is not set.
Thus currently has no effect.
Should consider a two stage approach, stating that constant nullspace is "trivial", and make another one for elasticity,
following guidelines form either:
https://colab.research.google.com/github/firedrakeproject/notebooks/blob/main/03-elasticity.ipynb
or
https://docs.mfem.org/4.8/hypre_8cpp_source.html#l05296
for option tests.
Since the problems are symmetric, Im not sure why I chose gmres and not cg ...
Discussed with @MariusCausemann
Guida per i contributori
Nessuna guida per i contributori indicizzata per questo repository
Come iniziare
- Leggi tutta la issue e poi la guida ai contributi del progetto.
- Commenta sulla issue per dire che te ne occupi tu — evita che due persone facciano lo stesso lavoro.
- Fai un fork del repository e lavora su un branch.
- Apri una pull request che faccia riferimento al numero della issue.
Direzione di ricerca
Inizia esaminando il materiale sul nullspace nel tutorial e confrontalo con il notebook di elasticità di Firedrake collegato e con il codice sorgente di MFEM hypre. Verifica quindi i casi relativi alle opzioni dei solver diretti e iterativi descritti nell’issue, poi definisci le modifiche necessarie al tutorial per fornire indicazioni separate sul nullspace costante e sull’elasticità, includendo la discussione sulla scelta del solver.
Scritto dal modello di indicizzazione a partire dal testo della issue.
Valutazione
- Stack tecnologico
- jupyter-notebook
- Ambito
- documentation
- Tipo di issue
- Documentazione
- Difficoltà
- 5/5
- Tempo stimato
- Più di una settimana
- Stato di attività
- Ferma
- Chiarezza
- Da chiarire
- Idoneità per principianti
- 30/100