Die Verfügbarkeit der Ähnlichkeitstransformationsmatrizen
erlaubt die Anwendung einer Variante des klassischen
CG-Algorithmus [39][40].
In diesem Verfahren
treten sehr große oder sehr kleine Zahlen nur bei der
Berechnung innerer Produkte auf. Die Berechnung dieser
inneren Produkte kann durch eine geeignete Skalierung
auch in doppelter Genauigkeit erfolgen.
Die Iteration ist auf eine Dreitermrekursion aufgebaut.
Linksseitige Vorkonditionierung wird verwendet.
Dieser Algorithmus minimiert bei jeder Iteration
die -Norm
des Lösungsfehlers
( die exakte Lösung):