LAPACK 3.12.0 LAPACK: Linear Algebra PACKage
Searching...
No Matches
crzt02.f
Go to the documentation of this file.
1*> \brief \b CRZT02
2*
3* =========== DOCUMENTATION ===========
4*
5* Online html documentation available at
6* http://www.netlib.org/lapack/explore-html/
7*
8* Definition:
9* ===========
10*
11* REAL FUNCTION CRZT02( M, N, AF, LDA, TAU, WORK,
12* LWORK )
13*
14* .. Scalar Arguments ..
15* INTEGER LDA, LWORK, M, N
16* ..
17* .. Array Arguments ..
18* COMPLEX AF( LDA, * ), TAU( * ), WORK( LWORK )
19* ..
20*
21*
22*> \par Purpose:
23* =============
24*>
25*> \verbatim
26*>
27*> CRZT02 returns
28*> || I - Q'*Q || / ( M * eps)
29*> where the matrix Q is defined by the Householder transformations
30*> generated by CTZRZF.
31*> \endverbatim
32*
33* Arguments:
34* ==========
35*
36*> \param[in] M
37*> \verbatim
38*> M is INTEGER
39*> The number of rows of the matrix AF.
40*> \endverbatim
41*>
42*> \param[in] N
43*> \verbatim
44*> N is INTEGER
45*> The number of columns of the matrix AF.
46*> \endverbatim
47*>
48*> \param[in] AF
49*> \verbatim
50*> AF is COMPLEX array, dimension (LDA,N)
51*> The output of CTZRZF.
52*> \endverbatim
53*>
54*> \param[in] LDA
55*> \verbatim
56*> LDA is INTEGER
57*> The leading dimension of the array AF.
58*> \endverbatim
59*>
60*> \param[in] TAU
61*> \verbatim
62*> TAU is COMPLEX array, dimension (M)
63*> Details of the Householder transformations as returned by
64*> CTZRZF.
65*> \endverbatim
66*>
67*> \param[out] WORK
68*> \verbatim
69*> WORK is COMPLEX array, dimension (LWORK)
70*> \endverbatim
71*>
72*> \param[in] LWORK
73*> \verbatim
74*> LWORK is INTEGER
75*> Length of WORK array. LWORK >= N*N+N.
76*> \endverbatim
77*
78* Authors:
79* ========
80*
81*> \author Univ. of Tennessee
82*> \author Univ. of California Berkeley
83*> \author Univ. of Colorado Denver
84*> \author NAG Ltd.
85*
86*> \ingroup complex_lin
87*
88* =====================================================================
89 REAL function crzt02( m, n, af, lda, tau, work,
90 \$ lwork )
91*
92* -- LAPACK test routine --
93* -- LAPACK is a software package provided by Univ. of Tennessee, --
94* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
95*
96* .. Scalar Arguments ..
97 INTEGER lda, lwork, m, n
98* ..
99* .. Array Arguments ..
100 COMPLEX af( lda, * ), tau( * ), work( lwork )
101* ..
102*
103* =====================================================================
104*
105* .. Parameters ..
106 REAL zero, one
107 parameter( zero = 0.0e0, one = 1.0e0 )
108* ..
109* .. Local Scalars ..
110 INTEGER i, info
111* ..
112* .. Local Arrays ..
113 REAL rwork( 1 )
114* ..
115* .. External Functions ..
116 REAL clange, slamch
117 EXTERNAL clange, slamch
118* ..
119* .. External Subroutines ..
120 EXTERNAL claset, cunmrz, xerbla
121* ..
122* .. Intrinsic Functions ..
123 INTRINSIC cmplx, max, real
124* ..
125* .. Executable Statements ..
126*
127 crzt02 = zero
128*
129 IF( lwork.LT.n*n+n ) THEN
130 CALL xerbla( 'CRZT02', 7 )
131 RETURN
132 END IF
133*
134* Quick return if possible
135*
136 IF( m.LE.0 .OR. n.LE.0 )
137 \$ RETURN
138*
139* Q := I
140*
141 CALL claset( 'Full', n, n, cmplx( zero ), cmplx( one ), work, n )
142*
143* Q := P(1) * ... * P(m) * Q
144*
145 CALL cunmrz( 'Left', 'No transpose', n, n, m, n-m, af, lda, tau,
146 \$ work, n, work( n*n+1 ), lwork-n*n, info )
147*
148* Q := P(m)' * ... * P(1)' * Q
149*
150 CALL cunmrz( 'Left', 'Conjugate transpose', n, n, m, n-m, af, lda,
151 \$ tau, work, n, work( n*n+1 ), lwork-n*n, info )
152*
153* Q := Q - I
154*
155 DO 10 i = 1, n
156 work( ( i-1 )*n+i ) = work( ( i-1 )*n+i ) - one
157 10 CONTINUE
158*
159 crzt02 = clange( 'One-norm', n, n, work, n, rwork ) /
160 \$ ( slamch( 'Epsilon' )*real( max( m, n ) ) )
161 RETURN
162*
163* End of CRZT02
164*
165 END
subroutine xerbla(srname, info)
Definition cblat2.f:3285
real function crzt02(m, n, af, lda, tau, work, lwork)
CRZT02
Definition crzt02.f:91
real function slamch(cmach)
SLAMCH
Definition slamch.f:68
real function clange(norm, m, n, a, lda, work)
CLANGE returns the value of the 1-norm, Frobenius norm, infinity-norm, or the largest absolute value ...
Definition clange.f:115
subroutine claset(uplo, m, n, alpha, beta, a, lda)
CLASET initializes the off-diagonal elements and the diagonal elements of a matrix to given values.
Definition claset.f:106
subroutine cunmrz(side, trans, m, n, k, l, a, lda, tau, c, ldc, work, lwork, info)
CUNMRZ
Definition cunmrz.f:187