73 parameter ( nmax = 2 )
76 INTEGER I, INFO, J, NB
79 COMPLEX*16 A( nmax, nmax ), T( nmax, nmax ), W( nmax ),
80 $ c( nmax, nmax ), tau(nmax)
92 COMMON / infoc / infot, nout, ok, lerr
93 COMMON / srnamc / srnamt
101 WRITE( nout, fmt = * )
107 a( i, j ) = 1.d0 / dble( i+j )
108 c( i, j ) = 1.d0 / dble( i+j )
109 t( i, j ) = 1.d0 / dble( i+j )
121 CALL zgeqr( -1, 0, a, 1, tau, 1, w, 1, info )
122 CALL chkxer(
'ZGEQR', infot, nout, lerr, ok )
124 CALL zgeqr( 0, -1, a, 1, tau, 1, w, 1, info )
125 CALL chkxer(
'ZGEQR', infot, nout, lerr, ok )
127 CALL zgeqr( 1, 1, a, 0, tau, 1, w, 1, info )
128 CALL chkxer(
'ZGEQR', infot, nout, lerr, ok )
130 CALL zgeqr( 3, 2, a, 3, tau, 1, w, 1, info )
131 CALL chkxer(
'ZGEQR', infot, nout, lerr, ok )
133 CALL zgeqr( 3, 2, a, 3, tau, 8, w, 0, info )
134 CALL chkxer(
'ZGEQR', infot, nout, lerr, ok )
143 CALL zgemqr(
'/',
'N', 0, 0, 0, a, 1, tau, 1, c, 1, w, 1,info)
144 CALL chkxer(
'ZGEMQR', infot, nout, lerr, ok )
146 CALL zgemqr(
'L',
'/', 0, 0, 0, a, 1, tau, 1, c, 1, w, 1,info)
147 CALL chkxer(
'ZGEMQR', infot, nout, lerr, ok )
149 CALL zgemqr(
'L',
'N', -1, 0, 0, a, 1, tau, 1, c, 1, w,1,info)
150 CALL chkxer(
'ZGEMQR', infot, nout, lerr, ok )
152 CALL zgemqr(
'L',
'N', 0, -1, 0, a, 1, tau, 1, c, 1, w,1,info)
153 CALL chkxer(
'ZGEMQR', infot, nout, lerr, ok )
155 CALL zgemqr(
'L',
'N', 0, 0, -1, a, 1, tau, 1, c, 1, w,1,info)
156 CALL chkxer(
'ZGEMQR', infot, nout, lerr, ok )
158 CALL zgemqr(
'R',
'N', 0, 0, -1, a, 1, tau, 1, c, 1, w,1,info)
159 CALL chkxer(
'ZGEMQR', infot, nout, lerr, ok )
161 CALL zgemqr(
'L',
'N', 2, 1, 0, a, 0, tau, 1, c, 1, w, 1,info)
162 CALL chkxer(
'ZGEMQR', infot, nout, lerr, ok )
164 CALL zgemqr(
'R',
'N', 2, 2, 1, a, 2, tau, 0, c, 1, w, 1,info)
165 CALL chkxer(
'ZGEMQR', infot, nout, lerr, ok )
167 CALL zgemqr(
'L',
'N', 2, 2, 1, a, 2, tau, 0, c, 1, w, 1,info)
168 CALL chkxer(
'ZGEMQR', infot, nout, lerr, ok )
170 CALL zgemqr(
'L',
'N', 2, 1, 1, a, 2, tau, 6, c, 0, w, 1,info)
171 CALL chkxer(
'ZGEMQR', infot, nout, lerr, ok )
173 CALL zgemqr(
'L',
'N', 2, 2, 1, a, 2, tau, 6, c, 2, w, 0,info)
174 CALL chkxer(
'ZGEMQR', infot, nout, lerr, ok )
180 CALL zgelq( -1, 0, a, 1, tau, 1, w, 1, info )
181 CALL chkxer(
'ZGELQ', infot, nout, lerr, ok )
183 CALL zgelq( 0, -1, a, 1, tau, 1, w, 1, info )
184 CALL chkxer(
'ZGELQ', infot, nout, lerr, ok )
186 CALL zgelq( 1, 1, a, 0, tau, 1, w, 1, info )
187 CALL chkxer(
'ZGELQ', infot, nout, lerr, ok )
189 CALL zgelq( 2, 3, a, 3, tau, 1, w, 1, info )
190 CALL chkxer(
'ZGELQ', infot, nout, lerr, ok )
192 CALL zgelq( 2, 3, a, 3, tau, 8, w, 0, info )
193 CALL chkxer(
'ZGELQ', infot, nout, lerr, ok )
202 CALL zgemlq(
'/',
'N', 0, 0, 0, a, 1, tau, 1, c, 1, w, 1,info)
203 CALL chkxer(
'ZGEMLQ', infot, nout, lerr, ok )
205 CALL zgemlq(
'L',
'/', 0, 0, 0, a, 1, tau, 1, c, 1, w, 1,info)
206 CALL chkxer(
'ZGEMLQ', infot, nout, lerr, ok )
208 CALL zgemlq(
'L',
'N', -1, 0, 0, a, 1, tau, 1, c, 1, w,1,info)
209 CALL chkxer(
'ZGEMLQ', infot, nout, lerr, ok )
211 CALL zgemlq(
'L',
'N', 0, -1, 0, a, 1, tau, 1, c, 1, w,1,info)
212 CALL chkxer(
'ZGEMLQ', infot, nout, lerr, ok )
214 CALL zgemlq(
'L',
'N', 0, 0, -1, a, 1, tau, 1, c, 1, w,1,info)
215 CALL chkxer(
'ZGEMLQ', infot, nout, lerr, ok )
217 CALL zgemlq(
'R',
'N', 0, 0, -1, a, 1, tau, 1, c, 1, w,1,info)
218 CALL chkxer(
'ZGEMLQ', infot, nout, lerr, ok )
220 CALL zgemlq(
'L',
'N', 1, 2, 0, a, 0, tau, 1, c, 1, w, 1,info)
221 CALL chkxer(
'ZGEMLQ', infot, nout, lerr, ok )
223 CALL zgemlq(
'R',
'N', 2, 2, 1, a, 1, tau, 0, c, 1, w, 1,info)
224 CALL chkxer(
'ZGEMLQ', infot, nout, lerr, ok )
226 CALL zgemlq(
'L',
'N', 2, 2, 1, a, 1, tau, 0, c, 1, w, 1,info)
227 CALL chkxer(
'ZGEMLQ', infot, nout, lerr, ok )
229 CALL zgemlq(
'L',
'N', 1, 2, 1, a, 1, tau, 6, c, 0, w, 1,info)
230 CALL chkxer(
'ZGEMLQ', infot, nout, lerr, ok )
232 CALL zgemlq(
'L',
'N', 2, 2, 1, a, 2, tau, 6, c, 2, w, 0,info)
233 CALL chkxer(
'ZGEMLQ', infot, nout, lerr, ok )
237 CALL alaesm( path, ok, nout )
subroutine zgeqr(M, N, A, LDA, T, TSIZE, WORK, LWORK, INFO)
subroutine alaesm(PATH, OK, NOUT)
ALAESM
subroutine zerrtsqr(PATH, NUNIT)
ZERRTSQR
subroutine chkxer(SRNAMT, INFOT, NOUT, LERR, OK)
subroutine zgemqr(SIDE, TRANS, M, N, K, A, LDA, T, TSIZE, C, LDC, WORK, LWORK, INFO)
subroutine zgelq(M, N, A, LDA, T, TSIZE, WORK, LWORK, INFO)
subroutine zgemlq(SIDE, TRANS, M, N, K, A, LDA, T, TSIZE, C, LDC, WORK, LWORK, INFO)