jorgensd / jorgensd/dolfinx-tutorial

Rethink nullspace tutorials

Aperta
#286 1 commento 0 reazioni 0 assegnatari Vedi su GitHub

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

  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.

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

Ricevi le nuove issue nella tua casella

Un breve riepilogo di issue GitHub adatte ai principianti.