cqrt11()

real function cqrt11	(	integer	m,
		integer	k,
		complex, dimension( lda, * )	a,
		integer	lda,
		complex, dimension( * )	tau,
		complex, dimension( lwork )	work,
		integer	lwork )

CQRT11

Purpose:

!>
!> CQRT11 computes the test ratio
!>
!>       || Q'*Q - I || / (eps * m)
!>
!> where the orthogonal matrix Q is represented as a product of
!> elementary transformations.  Each transformation has the form
!>
!>    H(k) = I - tau(k) v(k) v(k)'
!>
!> where tau(k) is stored in TAU(k) and v(k) is an m-vector of the form
!> [ 0 ... 0 1 x(k) ]', where x(k) is a vector of length m-k stored
!> in A(k+1:m,k).
!> 

Parameters

[in]	M	!> M is INTEGER !> The number of rows of the matrix A. !>
[in]	K	!> K is INTEGER !> The number of columns of A whose subdiagonal entries !> contain information about orthogonal transformations. !>
[in]	A	!> A is COMPLEX array, dimension (LDA,K) !> The (possibly partial) output of a QR reduction routine. !>
[in]	LDA	!> LDA is INTEGER !> The leading dimension of the array A. !>
[in]	TAU	!> TAU is COMPLEX array, dimension (K) !> The scaling factors tau for the elementary transformations as !> computed by the QR factorization routine. !>
[out]	WORK	!> WORK is COMPLEX array, dimension (LWORK) !>
[in]	LWORK	!> LWORK is INTEGER !> The length of the array WORK. LWORK >= M*M + M. !>

Author

Univ. of Tennessee

Univ. of California Berkeley

Univ. of Colorado Denver

NAG Ltd.

Definition at line 97 of file cqrt11.f.

*
*  -- LAPACK test routine --
*  -- LAPACK is a software package provided by Univ. of Tennessee,    --
*  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
*
*     .. Scalar Arguments ..
      INTEGER            K, LDA, LWORK, M
*     ..
*     .. Array Arguments ..
      COMPLEX            A( LDA, * ), TAU( * ), WORK( LWORK )
*     ..
*
*  =====================================================================
*
*     .. Parameters ..
      REAL               ZERO, ONE
      parameter( zero = 0.0e0, one = 1.0e0 )
*     ..
*     .. Local Scalars ..
      INTEGER            INFO, J
*     ..
*     .. External Functions ..
      REAL               CLANGE, SLAMCH
      EXTERNAL           clange, slamch
*     ..
*     .. External Subroutines ..
      EXTERNAL           claset, cunm2r, xerbla
*     ..
*     .. Intrinsic Functions ..
      INTRINSIC          cmplx, real
*     ..
*     .. Local Arrays ..
      REAL               RDUMMY( 1 )
*     ..
*     .. Executable Statements ..
*
      cqrt11 = zero
*
*     Test for sufficient workspace
*
      IF( lwork.LT.m*m+m ) THEN
         CALL xerbla( 'CQRT11', 7 )
         RETURN
      END IF
*
*     Quick return if possible
*
      IF( m.LE.0 )
     $   RETURN
*
      CALL claset( 'Full', m, m, cmplx( zero ), cmplx( one ), work, m )
*
*     Form Q
*
      CALL cunm2r( 'Left', 'No transpose', m, m, k, a, lda, tau, work,
     $             m, work( m*m+1 ), info )
*
*     Form Q'*Q
*
      CALL cunm2r( 'Left', 'Conjugate transpose', m, m, k, a, lda, tau,
     $             work, m, work( m*m+1 ), info )
*
      DO j = 1, m
         work( ( j-1 )*m+j ) = work( ( j-1 )*m+j ) - one
      END DO
*
      cqrt11 = clange( 'One-norm', m, m, work, m, rdummy ) /
     $         ( real( m )*slamch( 'Epsilon' ) )
*
      RETURN
*
*     End of CQRT11
*

Here is the call graph for this function:

Here is the caller graph for this function: