ISSN:
0029-5981
Keywords:
Engineering
;
Engineering General
Source:
Wiley InterScience Backfile Collection 1832-2000
Topics:
Mathematics
,
Technology
Notes:
FORTRAN subprograms for the solution of systems of linear algebraic equations are evaluated and compared on the basis of execution speed and accuracy. A symmetric, positive definite, banded test matrix is used in each case. The procedures considered are direct solution, iteration and matrix inversion. Both in-core schemes and those requiring the use of auxiliary data storage devices are included. Some of the techniques used require the full coefficient matrix, whereas others account for symmetry, banding, or sparseness of the system.Matrix inversion is found to be an inefficient technique, even if multiple right-hand side constant vectors are to be solved. In such cases, either Gauss elimination with multiple constant vectors treated simultaneously, or decomposition with the retention of the upper and lower triangular matrices, is recommended. Double precision arithmetic is suggested as a means of reducing round-off errors, and should always be employed when permitted by the computer's core capacity.On the basis of the results obtained with the single test system of equations, specific subprograms are recommended for each category of problem.
Additional Material:
4 Tab.
Type of Medium:
Electronic Resource
URL:
http://dx.doi.org/10.1002/nme.1620070406
Permalink
