79 COMPLEX*16 A( nmax, nmax ), T( nmax, nmax ), W( nmax ),
92 COMMON / infoc / infot, nout, ok, lerr
93 COMMON / srnamc / srnamt
96 INTRINSIC dble, dcmplx
101 WRITE( nout, fmt = * )
107 a( i, j ) = 1.d0 / dcmplx( dble( i+j ), 0.d0 )
108 c( i, j ) = 1.d0 / dcmplx( dble( i+j ), 0.d0 )
109 t( i, j ) = 1.d0 / dcmplx( dble( i+j ), 0.d0 )
121 CALL
zgeqrt( -1, 0, 1, a, 1, t, 1, w, info )
122 CALL
chkxer(
'ZGEQRT', infot, nout, lerr, ok )
124 CALL
zgeqrt( 0, -1, 1, a, 1, t, 1, w, info )
125 CALL
chkxer(
'ZGEQRT', infot, nout, lerr, ok )
127 CALL
zgeqrt( 0, 0, 0, a, 1, t, 1, w, info )
128 CALL
chkxer(
'ZGEQRT', infot, nout, lerr, ok )
130 CALL
zgeqrt( 2, 1, 1, a, 1, t, 1, w, info )
131 CALL
chkxer(
'ZGEQRT', infot, nout, lerr, ok )
133 CALL
zgeqrt( 2, 2, 2, a, 2, t, 1, w, info )
134 CALL
chkxer(
'ZGEQRT', infot, nout, lerr, ok )
140 CALL
zgeqrt2( -1, 0, a, 1, t, 1, info )
141 CALL
chkxer(
'ZGEQRT2', infot, nout, lerr, ok )
143 CALL
zgeqrt2( 0, -1, a, 1, t, 1, info )
144 CALL
chkxer(
'ZGEQRT2', infot, nout, lerr, ok )
146 CALL
zgeqrt2( 2, 1, a, 1, t, 1, info )
147 CALL
chkxer(
'ZGEQRT2', infot, nout, lerr, ok )
149 CALL
zgeqrt2( 2, 2, a, 2, t, 1, info )
150 CALL
chkxer(
'ZGEQRT2', infot, nout, lerr, ok )
156 CALL
zgeqrt3( -1, 0, a, 1, t, 1, info )
157 CALL
chkxer(
'ZGEQRT3', infot, nout, lerr, ok )
159 CALL
zgeqrt3( 0, -1, a, 1, t, 1, info )
160 CALL
chkxer(
'ZGEQRT3', infot, nout, lerr, ok )
162 CALL
zgeqrt3( 2, 1, a, 1, t, 1, info )
163 CALL
chkxer(
'ZGEQRT3', infot, nout, lerr, ok )
165 CALL
zgeqrt3( 2, 2, a, 2, t, 1, info )
166 CALL
chkxer(
'ZGEQRT3', infot, nout, lerr, ok )
172 CALL
zgemqrt(
'/',
'N', 0, 0, 0, 1, a, 1, t, 1, c, 1, w, info )
173 CALL
chkxer(
'ZGEMQRT', infot, nout, lerr, ok )
175 CALL
zgemqrt(
'L',
'/', 0, 0, 0, 1, a, 1, t, 1, c, 1, w, info )
176 CALL
chkxer(
'ZGEMQRT', infot, nout, lerr, ok )
178 CALL
zgemqrt(
'L',
'N', -1, 0, 0, 1, a, 1, t, 1, c, 1, w, info )
179 CALL
chkxer(
'ZGEMQRT', infot, nout, lerr, ok )
181 CALL
zgemqrt(
'L',
'N', 0, -1, 0, 1, a, 1, t, 1, c, 1, w, info )
182 CALL
chkxer(
'ZGEMQRT', infot, nout, lerr, ok )
184 CALL
zgemqrt(
'L',
'N', 0, 0, -1, 1, a, 1, t, 1, c, 1, w, info )
185 CALL
chkxer(
'ZGEMQRT', infot, nout, lerr, ok )
187 CALL
zgemqrt(
'R',
'N', 0, 0, -1, 1, a, 1, t, 1, c, 1, w, info )
188 CALL
chkxer(
'ZGEMQRT', infot, nout, lerr, ok )
190 CALL
zgemqrt(
'L',
'N', 0, 0, 0, 0, a, 1, t, 1, c, 1, w, info )
191 CALL
chkxer(
'ZGEMQRT', infot, nout, lerr, ok )
193 CALL
zgemqrt(
'R',
'N', 1, 2, 1, 1, a, 1, t, 1, c, 1, w, info )
194 CALL
chkxer(
'ZGEMQRT', infot, nout, lerr, ok )
196 CALL
zgemqrt(
'L',
'N', 2, 1, 1, 1, a, 1, t, 1, c, 1, w, info )
197 CALL
chkxer(
'ZGEMQRT', infot, nout, lerr, ok )
199 CALL
zgemqrt(
'R',
'N', 1, 1, 1, 1, a, 1, t, 0, c, 1, w, info )
200 CALL
chkxer(
'ZGEMQRT', infot, nout, lerr, ok )
202 CALL
zgemqrt(
'L',
'N', 1, 1, 1, 1, a, 1, t, 1, c, 0, w, info )
203 CALL
chkxer(
'ZGEMQRT', infot, nout, lerr, ok )
207 CALL
alaesm( path, ok, nout )
subroutine zgeqrt2(M, N, A, LDA, T, LDT, INFO)
ZGEQRT2 computes a QR factorization of a general real or complex matrix using the compact WY represen...
subroutine zgemqrt(SIDE, TRANS, M, N, K, NB, V, LDV, T, LDT, C, LDC, WORK, INFO)
ZGEMQRT
subroutine chkxer(SRNAMT, INFOT, NOUT, LERR, OK)
recursive subroutine zgeqrt3(M, N, A, LDA, T, LDT, INFO)
ZGEQRT3 recursively computes a QR factorization of a general real or complex matrix using the compact...
subroutine zerrqrt(PATH, NUNIT)
ZERRQRT
subroutine alaesm(PATH, OK, NOUT)
ALAESM
subroutine zgeqrt(M, N, NB, A, LDA, T, LDT, WORK, INFO)
ZGEQRT