|
[Sponsors] |
March 13, 2021, 17:23 |
ILU(0) preconditioner in SU2
|
#1 |
New Member
Bruan
Join Date: Mar 2021
Posts: 1
Rep Power: 0 |
Does anyone know how parallelism is achieved in highly serial ILU and still gives relatively good performance?
|
|
March 14, 2021, 09:56 |
|
#2 |
Senior Member
Pedro Gomes
Join Date: Dec 2017
Posts: 466
Rep Power: 14 |
In SU2 (and many codes) ILU (and LU-SGS) operate only on the local matrix of each MPI rank.
This is called "Additive Schwarz" (makes it sound less hacky if you write a paper). That is, it is not the ILU of the entire global problem. But linear systems are still solved for the global problem that is why we get away with this approach. However, for the same tolerance the systems take a few more iterations to solve once you have many cores. The higher the CFL the worse this gets, and for elasticity problems it becomes very noticeable. There are ways to compute and apply the full ILU in parallel by transforming it in an iterative problem (a kind of Jacobi iteration to compute it, and then to form the solves with L and U parts). These approaches are common for massive parallel architectures (GPU's) but they would not perform well on CPU due to much less memory bandwidth (unless you use a number of CPU comparable to the number of threads in a GPU). |
|
Tags |
ilu, parallel, preconditioner, su2 |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
SU2 7.0.7 Built on CentOS 7, parallel computation pyscript mpi exit error? | EternalSeekerX | SU2 | 3 | October 9, 2020 19:28 |
SU2 Divergence Problem | Gui_AP | SU2 | 5 | February 21, 2020 12:10 |
Introducing SU2 International Developers Society (IDS) | fpalacios | SU2 News & Announcements | 1 | June 17, 2019 23:38 |
Official release of SU2 V3.0 and SU2 Educational V1.0 | fpalacios | SU2 News & Announcements | 2 | January 22, 2014 06:28 |
Welcome to the Stanford University Unstructured (SU2) forum! | economon | SU2 | 0 | January 7, 2013 03:48 |