Solution of Linear Algebraic Equations part 10

void dsprstx(double sa[], unsigned long ija[], double x[], double b[], unsigned long n); These are double versions of sprsax and sprstx. if (itrnsp) dsprstx(sa,ija,x,r,n); else dsprsax(sa,ija,x,r,n); } | Sparse Linear Systems 89 void dsprstx double sa unsigned long ija double x double b unsigned long n These are double versions of sprsax and sprstx. if itrnsp dsprstx sa ija x r n else dsprsax sa ija x r n extern unsigned long ija extern double sa The matrix is stored somewhere. void asolve unsigned long n double b double x int itrnsp unsigned long i for i 1 i n i x i sa i b i sa i b i The matrix A is the diagonal part of A stored in the first n elements of sa. Since the transpose matrix has the same diagonal the flag itrnsp is not used. CITED REFERENCES AND FURTHER READING Tewarson . 1973 Sparse Matrices New York Academic Press . 1 Jacobs . ed. 1977 The State of the Art in Numerical Analysis London Academic Press Chapter by . Reid . 2 George A. and Liu . 1981 Computer Solution ofLarge Sparse Positive Definite Systems Englewood Cliffs NJ Prentice-Hall . 3 NAG Fortran Library Numerical Algorithms Group 256 Banbury Road Oxford OX27DE . . 4 IMSL Math Library UsersManual IMSL Inc. 2500 CityWest Boulevard Houston TX 77042 . 5 Eisenstat . Gursky . Schultz . and Sherman . 1977 Yale Sparse Matrix Package Technical Reports 112 and 114 Yale University Department of Computer Science . 6 Knuth . 1968 FundamentalAlgorithms vol. 1 of The Art of ComputerProgramming Reading MA Addison-Wesley . 7 Kincaid . Respess . Young . and Grimes . 1982 ACM Transactions on Mathematical Software vol. 8 pp. 302-322. 8 PCGPAK User s Guide New Haven Scientific Computing Associates Inc. . 9 Bentley J. 1986 Programming Pearls Reading MA Addison-Wesley 9. 10 Golub . and Van Loan . 1989 Matrix Computations 2nd ed. Baltimore Johns Hopkins University Press Chapters 4 and 10 particularly . 11 Stoer J. and Bulirsch R. 1980 Introduction to NumericalAnalysis New York Springer-Verlag Chapter 8. 12 Baker L. 1991 More C Tools for Scientists and Engineers New York McGraw-Hill . 13 Fletcher R. 1976 in Numerical Analysis Dundee 1975 .