next up previous contents
Next: 5.4.3 Red-Black-Gauß-Seidel iteration Up: 5.4 Gauß-Seidel iteration Previous: 5.4.1 Sequential algorithm   Contents


5.4.2 Data flow of the Gauß-Seidel forward iteration

As presented in Alg. 5.7 and Fig. 5.2 the components of  $ \underline{u}^k$ are coupled, e.g., $ u_i^k$ cannot be calculated before $ u_{i-1}^k$ was determined. In general, all $ u_{s}^k$ ($ s<i$) for which $ K_{i,s}\neq 0$ (graph/pattern of matrix) holds have to be calculated before $ u_i^k$ can be updated.

Figure 5.2: Data flow of Gauß-Seidel forward iteration
\begin{figure}\unitlength0.035\textwidth
\mbox{}\hfill
\begin{picture}(14,6)
%
...
...13,1){\makebox(0,0){$u_{j}^{k}$}}
%
\par\end{picture}\hfill\mbox{}
\end{figure}


Due to that coupling of the components of $ \underline{u}^k$ the direct parallelization will become fine grain and the vectorization makes no sense.

Gundolf Haase 2000-03-20