5.5.3 Die parallele IUL-Faktorisierung

Next: 5.5.4 Vergleich von ILU- Up: 5.5 Unvollständige Zerlegungen Previous: 5.5.2 Die parallele ILU-Faktorisierung

5.5.3 Die parallele IUL-Faktorisierung

Zur Überwindung der doppelten Kommunikation in Algorithmus 5.15 wird eine unvollständige ${\ensuremath{\color{red} \mathfrak{U} } }{\ensuremath{\color{red} \mathfrak{L} } }$ -Faktorisierung benutzt.

Der Durchlaufsinn kehrt sich gegenüber der ILU-Faktorisierung um. Dadurch werden die redundant gespeicherten Matrizen ${\ensuremath{\color{red} \mathfrak{K} } }_V$ , ${\ensuremath{\color{red} \mathfrak{K} } }_E$ , ${\ensuremath{\color{red} \mathfrak{K} } }_E$ , ${\ensuremath{\color{red} \mathfrak{K} } }_{EV}$ und ${\ensuremath{\color{red} \mathfrak{K} } }_{VE}$ während der Faktorisierung lokal verändert, so daß ihre Akkumulation erst kurz vor der entsprechenden Faktorisierung ausgeführt werden darf !
$\begin{algorithmus}% latex2html id marker 20292 [H]\caption{Parallelisierte IUL-... ...ce & mesh $\rightarrow$\space diagonal matrix \end{tabular} \end{algorithmus}$
Der Auflösungsschritt ist nur noch eine Anwendung von (4.12a).
$\begin{algorithmus}% latex2html id marker 20424 [H]\caption{Paralleler Aufl\uml ... ...ath{\color{red} \mathfrak{w} } }}_V)$\space \\ \end{tabular} \end{algorithmus}$

Gundolf Haase
1998-12-22