73 INTEGER I, INFO, J, MB, NB
76 DOUBLE PRECISION A( NMAX, NMAX ), T( NMAX, NMAX ), W( NMAX ),
77 $ C( NMAX, NMAX ), TAU(NMAX*2)
89 COMMON / infoc / infot, nout, ok, lerr
90 COMMON / srnamc / srnamt
98 WRITE( nout, fmt = * )
104 a( i, j ) = 1.d0 / dble( i+j )
105 c( i, j ) = 1.d0 / dble( i+j )
106 t( i, j ) = 1.d0 / dble( i+j )
118 CALL dgeqr( -1, 0, a, 1, tau, 1, w, 1, info )
119 CALL chkxer(
'DGEQR', infot, nout, lerr, ok )
121 CALL dgeqr( 0, -1, a, 1, tau, 1, w, 1, info )
122 CALL chkxer(
'DGEQR', infot, nout, lerr, ok )
124 CALL dgeqr( 1, 1, a, 0, tau, 1, w, 1, info )
125 CALL chkxer(
'DGEQR', infot, nout, lerr, ok )
127 CALL dgeqr( 3, 2, a, 3, tau, 1, w, 1, info )
128 CALL chkxer(
'DGEQR', infot, nout, lerr, ok )
130 CALL dgeqr( 3, 2, a, 3, tau, 7, w, 0, info )
131 CALL chkxer(
'DGEQR', infot, nout, lerr, ok )
139 CALL dlatsqr( -1, 0, mb, nb, a, 1, tau, 1, w, 1, info )
140 CALL chkxer(
'DLATSQR', infot, nout, lerr, ok )
142 CALL dlatsqr( 1, 2, mb, nb, a, 1, tau, 1, w, 1, info )
143 CALL chkxer(
'DLATSQR', infot, nout, lerr, ok )
144 CALL dlatsqr( 0, -1, mb, nb, a, 1, tau, 1, w, 1, info )
145 CALL chkxer(
'DLATSQR', infot, nout, lerr, ok )
147 CALL dlatsqr( 2, 1, -1, nb, a, 2, tau, 1, w, 1, info )
148 CALL chkxer(
'DLATSQR', infot, nout, lerr, ok )
150 CALL dlatsqr( 2, 1, mb, 2, a, 2, tau, 1, w, 1, info )
151 CALL chkxer(
'DLATSQR', infot, nout, lerr, ok )
153 CALL dlatsqr( 2, 1, mb, nb, a, 1, tau, 1, w, 1, info )
154 CALL chkxer(
'DLATSQR', infot, nout, lerr, ok )
156 CALL dlatsqr( 2, 1, mb, nb, a, 2, tau, 0, w, 1, info )
157 CALL chkxer(
'DLATSQR', infot, nout, lerr, ok )
159 CALL dlatsqr( 2, 1, mb, nb, a, 2, tau, 2, w, 0, info )
160 CALL chkxer(
'DLATSQR', infot, nout, lerr, ok )
171 CALL dgemqr(
'/',
'N', 0, 0, 0, a, 1, tau, 1, c, 1, w, 1,info)
172 CALL chkxer(
'DGEMQR', infot, nout, lerr, ok )
174 CALL dgemqr(
'L',
'/', 0, 0, 0, a, 1, tau, 1, c, 1, w, 1,info)
175 CALL chkxer(
'DGEMQR', infot, nout, lerr, ok )
177 CALL dgemqr(
'L',
'N', -1, 0, 0, a, 1, tau, 1, c, 1, w,1,info)
178 CALL chkxer(
'DGEMQR', infot, nout, lerr, ok )
180 CALL dgemqr(
'L',
'N', 0, -1, 0, a, 1, tau, 1, c, 1, w,1,info)
181 CALL chkxer(
'DGEMQR', infot, nout, lerr, ok )
183 CALL dgemqr(
'L',
'N', 0, 0, -1, a, 1, tau, 1, c, 1, w,1,info)
184 CALL chkxer(
'DGEMQR', infot, nout, lerr, ok )
186 CALL dgemqr(
'R',
'N', 0, 0, -1, a, 1, tau, 1, c, 1, w,1,info)
187 CALL chkxer(
'DGEMQR', infot, nout, lerr, ok )
189 CALL dgemqr(
'L',
'N', 2, 1, 0, a, 0, tau, 1, c, 1, w, 1,info)
190 CALL chkxer(
'DGEMQR', infot, nout, lerr, ok )
192 CALL dgemqr(
'R',
'N', 2, 2, 1, a, 2, tau, 0, c, 1, w, 1,info)
193 CALL chkxer(
'DGEMQR', infot, nout, lerr, ok )
195 CALL dgemqr(
'L',
'N', 2, 2, 1, a, 2, tau, 0, c, 1, w, 1,info)
196 CALL chkxer(
'DGEMQR', infot, nout, lerr, ok )
198 CALL dgemqr(
'L',
'N', 2, 1, 1, a, 2, tau, 6, c, 0, w, 1,info)
199 CALL chkxer(
'DGEMQR', infot, nout, lerr, ok )
201 CALL dgemqr(
'L',
'N', 2, 2, 1, a, 2, tau, 6, c, 2, w, 0,info)
202 CALL chkxer(
'DGEMQR', infot, nout, lerr, ok )
208 CALL dgelq( -1, 0, a, 1, tau, 1, w, 1, info )
209 CALL chkxer(
'DGELQ', infot, nout, lerr, ok )
211 CALL dgelq( 0, -1, a, 1, tau, 1, w, 1, info )
212 CALL chkxer(
'DGELQ', infot, nout, lerr, ok )
214 CALL dgelq( 1, 1, a, 0, tau, 1, w, 1, info )
215 CALL chkxer(
'DGELQ', infot, nout, lerr, ok )
217 CALL dgelq( 2, 3, a, 3, tau, 1, w, 1, info )
218 CALL chkxer(
'DGELQ', infot, nout, lerr, ok )
220 CALL dgelq( 2, 3, a, 3, tau, 7, w, 0, info )
221 CALL chkxer(
'DGELQ', infot, nout, lerr, ok )
229 CALL dlaswlq( -1, 0, mb, nb, a, 1, tau, 1, w, 1, info )
230 CALL chkxer(
'DLASWLQ', infot, nout, lerr, ok )
232 CALL dlaswlq( 2, 1, mb, nb, a, 1, tau, 1, w, 1, info )
233 CALL chkxer(
'DLASWLQ', infot, nout, lerr, ok )
234 CALL dlaswlq( 0, -1, mb, nb, a, 1, tau, 1, w, 1, info )
235 CALL chkxer(
'DLASWLQ', infot, nout, lerr, ok )
237 CALL dlaswlq( 1, 2, -1, nb, a, 1, tau, 1, w, 1, info )
238 CALL chkxer(
'DLASWLQ', infot, nout, lerr, ok )
239 CALL dlaswlq( 1, 1, 2, nb, a, 1, tau, 1, w, 1, info )
240 CALL chkxer(
'DLASWLQ', infot, nout, lerr, ok )
242 CALL dlaswlq( 1, 2, mb, -1, a, 1, tau, 1, w, 1, info )
243 CALL chkxer(
'DLASWLQ', infot, nout, lerr, ok )
245 CALL dlaswlq( 1, 2, mb, nb, a, 0, tau, 1, w, 1, info )
246 CALL chkxer(
'DLASWLQ', infot, nout, lerr, ok )
248 CALL dlaswlq( 1, 2, mb, nb, a, 1, tau, 0, w, 1, info )
249 CALL chkxer(
'DLASWLQ', infot, nout, lerr, ok )
251 CALL dlaswlq( 1, 2, mb, nb, a, 1, tau, 1, w, 0, info )
252 CALL chkxer(
'DLASWLQ', infot, nout, lerr, ok )
261 CALL dgemlq(
'/',
'N', 0, 0, 0, a, 1, tau, 1, c, 1, w, 1,info)
262 CALL chkxer(
'DGEMLQ', infot, nout, lerr, ok )
264 CALL dgemlq(
'L',
'/', 0, 0, 0, a, 1, tau, 1, c, 1, w, 1,info)
265 CALL chkxer(
'DGEMLQ', infot, nout, lerr, ok )
267 CALL dgemlq(
'L',
'N', -1, 0, 0, a, 1, tau, 1, c, 1, w,1,info)
268 CALL chkxer(
'DGEMLQ', infot, nout, lerr, ok )
270 CALL dgemlq(
'L',
'N', 0, -1, 0, a, 1, tau, 1, c, 1, w,1,info)
271 CALL chkxer(
'DGEMLQ', infot, nout, lerr, ok )
273 CALL dgemlq(
'L',
'N', 0, 0, -1, a, 1, tau, 1, c, 1, w,1,info)
274 CALL chkxer(
'DGEMLQ', infot, nout, lerr, ok )
276 CALL dgemlq(
'R',
'N', 0, 0, -1, a, 1, tau, 1, c, 1, w,1,info)
277 CALL chkxer(
'DGEMLQ', infot, nout, lerr, ok )
279 CALL dgemlq(
'L',
'N', 1, 2, 0, a, 0, tau, 1, c, 1, w, 1,info)
280 CALL chkxer(
'DGEMLQ', infot, nout, lerr, ok )
282 CALL dgemlq(
'R',
'N', 2, 2, 1, a, 1, tau, 0, c, 1, w, 1,info)
283 CALL chkxer(
'DGEMLQ', infot, nout, lerr, ok )
285 CALL dgemlq(
'L',
'N', 2, 2, 1, a, 1, tau, 0, c, 1, w, 1,info)
286 CALL chkxer(
'DGEMLQ', infot, nout, lerr, ok )
288 CALL dgemlq(
'L',
'N', 1, 2, 1, a, 1, tau, 6, c, 0, w, 1,info)
289 CALL chkxer(
'DGEMLQ', infot, nout, lerr, ok )
291 CALL dgemlq(
'L',
'N', 2, 2, 1, a, 2, tau, 6, c, 2, w, 0,info)
292 CALL chkxer(
'DGEMLQ', infot, nout, lerr, ok )
296 CALL alaesm( path, ok, nout )