© 1976 by British Computer Society
| ||||||||||||||||||||||||||||||||||||||||||||||||||||
A direct method for the solution of large sparse systems of linear equations
Department of Chemical Engineering, Ben Gurion University of the Negev, Beer-Sheva, Israel
A new method for the solution of large sparse systems of linear equations that offers considerable savings of computation time and storage is presented. The proposed method is easily programmable for linear systems with a band type matrix of coefficients and can be used for large sparse systems, where the non-zero elements in the matrix of coefficients are randomly spaced.
The system is restructured by tearing and the whole system is solved as a smaller linear system of the tear variables, using a modification of the Newton method, for which the derivatives are calculated directly from the equations.
Received March 1974. Revised April 1975.
* Present Address: Department of Chemical Engineering, Ben Gurion University of the Negev, Beer-Sheva 84-120, Israel.
Department of Chemical Engineering, TerchnionIsrael Institute of Technology, Haifa, Israel