Purpose
To compute the suitable maps for Hessenberg index H and signature array S. Auxiliary routine for the periodic QZ algorithms.Specification
SUBROUTINE MB03BA( K, H, S, SMULT, AMAP, QMAP )
C .. Scalar Arguments ..
INTEGER K, H, SMULT
C .. Array Arguments ..
INTEGER AMAP(*), QMAP(*), S(*)
Arguments
Input/Output Parameters
K (input) INTEGER
The number of factors. K >= 1.
H (input) INTEGER
Index which corresponds to A_1.
S (input) INTEGER array, dimension (K)
The signature array. Each entry of S must be 1 or -1.
SMULT (output) INTEGER
Signature multiplier. Entries of S are virtually
multiplied by SMULT.
AMAP (output) INTEGER array, dimension (K)
The map for accessing the factors, that is,
if AMAP(I) = J, then the factor A_I is stored at the J-th
position in A.
QMAP (output) INTEGER array, dimension (K)
The map for accessing the orthognal transformation
matrices, that is, if QMAP(I) = J, then the matrix Q_I is
stored at the J-th position in Q.
Further Comments
NoneExample
Program Text
NoneProgram Data
NoneProgram Results
None