Most iterative techniques for solving the symmetric positive-definite systemAx=b involve approximating the matrixA by another symmetric positive-definite matrixM and then solving a system of the formMz=d at each iteration. On a vector machine such as the CDC-STAR-100, the solution of this new system can be very time consuming. If, however, an approximationM−1 can be given toA−1, the solutionz=M−1d can be computed rapidly by matrix multiplication, a fast operation on the STAR. Approximations using the Neumann expansion of the inverse ofA give reasonable forms forM−1 and are presented. Computational results using the conjugate gradient method for the “5-point” matrixA are given.ZusammenfassungDie meisten iterativen Methoden zur Lösung des symmetrischen positiv-definitiven SystemsAx=b enthalten die Näherung der MatrixA durch eine andere symmetrische positiv-definitive MatrixM und anschließend daran die Lösung eines Systems der ArtMz=d bei jeder Wiederholung. Auf einer Vektor-Maschine wie der CDC-STAR-100 kann die Lösung dieses neuen Systems sehr zeitraubend sein. Wenn jedoch eine NäherungM−1 zuA−1 gegeben werden kann, so kann die Lösungz=M−1d sehr schnell durch Matrixmultiplikation errechnet werden. Diese Kalkulation kann auf dem STAR schnell ausgeführt werden. Näherungen, bei denen die Neumann-Entwicklung der Inversen vonA verwendet wird, ergeben angemessene Ausdrücke fürM−1. Diese Ausdrücke sind angeführt. Die mit Hilfe der Konjugierten-Gradienten-Methode errechneten Resultate für die „5-Punk”-MatrixA sind angegeben.
[1]
G. Rodrigue,et al.
Incomplete Cholesky conjugate gradient method for the STAR (5-point operator)
,
1977
.
[2]
Louis A. Hageman,et al.
Iterative Solution of Large Linear Systems.
,
1971
.
[3]
Gene H. Golub,et al.
A generalized conjugate gradient method for the numerical solution of elliptic partial differential equations
,
2007,
Milestones in Matrix Computation.
[4]
A. K. Cline.
Several Observations on the Use of Conjugate Gradient Methods
,
1978
.
[5]
J. Gillis,et al.
Matrix Iterative Analysis
,
1961
.
[6]
J. Meijerink,et al.
An iterative solution method for linear systems of which the coefficient matrix is a symmetric -matrix
,
1977
.