RMNLIB |
(ROSSR3) |
BASE |
SUBROUTINE ROSSR3 (P, a, delta, c, d, m)
OUT
IN
A(1) and C(M) are not used. The general tri-diagonal problem is reducible to the above form by suitably scaling: A(i)P(i-1)+B(i)P(i)+C(i)P(i+1)=D(i) should be divided by B(i) (assuming B(i).ne.0) where A(1)=C(M)=0 If any of A, C or D need not be preserved, they may be used as actual arguments for the work vectors.