56 SUBROUTINE zerrlqt( PATH, NUNIT )
73 parameter ( nmax = 2 )
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 zgelqt( -1, 0, 1, a, 1, t, 1, w, info )
122 CALL chkxer(
'ZGELQT', infot, nout, lerr, ok )
124 CALL zgelqt( 0, -1, 1, a, 1, t, 1, w, info )
125 CALL chkxer(
'ZGELQT', infot, nout, lerr, ok )
127 CALL zgelqt( 0, 0, 0, a, 1, t, 1, w, info )
128 CALL chkxer(
'ZGELQT', infot, nout, lerr, ok )
130 CALL zgelqt( 2, 1, 1, a, 1, t, 1, w, info )
131 CALL chkxer(
'ZGELQT', infot, nout, lerr, ok )
133 CALL zgelqt( 2, 2, 2, a, 2, t, 1, w, info )
134 CALL chkxer(
'ZGELQT', infot, nout, lerr, ok )
140 CALL zgelqt3( -1, 0, a, 1, t, 1, info )
141 CALL chkxer(
'ZGELQT3', infot, nout, lerr, ok )
143 CALL zgelqt3( 0, -1, a, 1, t, 1, info )
144 CALL chkxer(
'ZGELQT3', infot, nout, lerr, ok )
146 CALL zgelqt3( 2, 2, a, 1, t, 1, info )
147 CALL chkxer(
'ZGELQT3', infot, nout, lerr, ok )
149 CALL zgelqt3( 2, 2, a, 2, t, 1, info )
150 CALL chkxer(
'ZGELQT3', infot, nout, lerr, ok )
156 CALL zgemlqt(
'/',
'N', 0, 0, 0, 1, a, 1, t, 1, c, 1, w, info )
157 CALL chkxer(
'ZGEMLQT', infot, nout, lerr, ok )
159 CALL zgemlqt(
'L',
'/', 0, 0, 0, 1, a, 1, t, 1, c, 1, w, info )
160 CALL chkxer(
'ZGEMLQT', infot, nout, lerr, ok )
162 CALL zgemlqt(
'L',
'N', -1, 0, 0, 1, a, 1, t, 1, c, 1, w, info )
163 CALL chkxer(
'ZGEMLQT', infot, nout, lerr, ok )
165 CALL zgemlqt(
'L',
'N', 0, -1, 0, 1, a, 1, t, 1, c, 1, w, info )
166 CALL chkxer(
'ZGEMLQT', infot, nout, lerr, ok )
168 CALL zgemlqt(
'L',
'N', 0, 0, -1, 1, a, 1, t, 1, c, 1, w, info )
169 CALL chkxer(
'ZGEMLQT', infot, nout, lerr, ok )
171 CALL zgemlqt(
'R',
'N', 0, 0, -1, 1, a, 1, t, 1, c, 1, w, info )
172 CALL chkxer(
'ZGEMLQT', infot, nout, lerr, ok )
174 CALL zgemlqt(
'L',
'N', 0, 0, 0, 0, a, 1, t, 1, c, 1, w, info )
175 CALL chkxer(
'ZGEMLQT', infot, nout, lerr, ok )
177 CALL zgemlqt(
'R',
'N', 2, 2, 2, 1, a, 1, t, 1, c, 1, w, info )
178 CALL chkxer(
'ZGEMLQT', infot, nout, lerr, ok )
180 CALL zgemlqt(
'L',
'N', 2, 2, 2, 1, a, 1, t, 1, c, 1, w, info )
181 CALL chkxer(
'ZGEMLQT', infot, nout, lerr, ok )
183 CALL zgemlqt(
'R',
'N', 1, 1, 1, 1, a, 1, t, 0, c, 1, w, info )
184 CALL chkxer(
'ZGEMLQT', infot, nout, lerr, ok )
186 CALL zgemlqt(
'L',
'N', 1, 1, 1, 1, a, 1, t, 1, c, 0, w, info )
187 CALL chkxer(
'ZGEMLQT', infot, nout, lerr, ok )
191 CALL alaesm( path, ok, nout )
subroutine zgelqt(M, N, MB, A, LDA, T, LDT, WORK, INFO)
ZGELQT
subroutine zgemlqt(SIDE, TRANS, M, N, K, MB, V, LDV, T, LDT, C, LDC, WORK, INFO)
ZGEMLQT
subroutine zerrlqt(PATH, NUNIT)
ZERLQT
subroutine alaesm(PATH, OK, NOUT)
ALAESM
recursive subroutine zgelqt3(M, N, A, LDA, T, LDT, INFO)
ZGELQT3 recursively computes a LQ factorization of a general real or complex matrix using the compact...
subroutine chkxer(SRNAMT, INFOT, NOUT, LERR, OK)