f01ref
f01ref
© Numerical Algorithms Group, 2002.
Purpose
F01REF Operations with unitary matrices, form columns of Q after
factorization by F01RCF
Synopsis
[a,ifail] = f01ref(a,theta<,ncolq,n,wheret,ifail>)
Description
The unitary matrix Q is assumed to be given by
H
Q=(Q Q ...Q ) ,
n n-1 1
Q being given in the form
k
(I 0 )
Q =(0 T ),
k ( k)
where
H
T =I-(gamma) u u ,
k k k k
((zeta) )
( k)
u =(z ),
k ( k )
(gamma) is a scalar for which Re (gamma) =1.0, (zeta) is a real
k k k
scalar and z is an (m-k) element vector.
k
z must be supplied in the kth column of A in elements
k
a ,...,a and (theta) , given by
k+1,k m,k k
(theta) =((zeta) ,Im (gamma) ),
k k k
must be supplied either in a or in THETA(k) depending upon the
k,k
parameter WHERET.
Parameters
f01ref
Required Input Arguments:
a (:,:) complex
theta (:) complex
Optional Input Arguments: <Default>
ncolq integer size(a,2)
n integer size(a,2)
wheret (1) string 's'
ifail integer -1
Output Arguments:
a (:,:) complex
ifail integer