On entry: the upper triangle of the
symmetric band matrix
, with the diagonal of the matrix stored in the
th row of the array, and the
superdiagonals within the band stored in the first
rows of the array. Each column of the matrix is stored in the corresponding column of the array. For example, if
and
, the storage scheme is
Elements in the top left corner of the array are not used. The matrix elements within the band can be assigned to the correct elements of the array using the following
code:
Do 20 j = 1, n
Do 10 i = max(1,j-m1+1), j
a(i-j+m1,j) = matrix(i,j)
End Do
End Do