****** |
i n p u t ******
w, nroww,nrow,nbandl,nbandu.....describe the lu-factorization of a
|
banded matrix | |
a of order nrow as constructed in banfac .
|
for details, see | |
banfac .
b.....right side of the system to be solved .
|
|
|
****** |
o u t p u t ******
|
b.....contains the solution | |
x , of order nrow .
|
|
|
****** |
m e t h o d ******
|
(with |
a = l*u, as stored in w,) the unit lower triangular system
|
|
|
l(u*x) = b |
|
is solved for y = u*x, and y stored in b . then the
|
upper triangular system | |
u*x = y is solved for x . the calcul-
ations are so arranged that the innermost loops stay within columns.
|
|
|
|