Algorithmen auf Graphen und dünn besetzte Matrizen - Bergische ...

1.1. DÜNN BESETZTE VEKTOREN1.1.6 Algorithmus (b)SAXP Y (n, ind x , ind y , val x , val y , nnz(x), nnz(y), α){Berechnet bei gepackten Vektoren x, y ein SAXPY y := αx + y. DasResultat steht in dem vollen Vektor w zur Verfügung und wird erst amEnde in y geschrieben.}{innerhalb eines GAXPYs nur beim ersten SAXPY notwendig}scatter (w, n, ind y , val y , nnz(y)) {w enthält y}j := nnz(y)for i = 1 to nnz(x) do {Lauf über x}if w indx(i) ≠ 0 thenw indx(i) := α · val x (i) + w indx(i)elsew indx(i) = α · val x (i)j + +ind y (j) = ind x (i)end ifend fornnz(y) := j{struct(y) ist bekannt, die Werte nicht}{innerhalb eines GAXPYs nur beim letzten SAXPY notwendig}for i = 1 to nnz(y) do{Lauf über y, sammelt die Werte ein}val y (i) := w indy(i)w indy(i) := 0end forAufwand: O(nnz(x) + nnz(y)), bei einfachem SAXPY etwas aufwendigerals Algorithmus 1.1.6.9