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