75 DOUBLE PRECISION A( NMAX, NMAX ), AF( NMAX, NMAX ), B( NMAX ),
76 $ W( NMAX ), X( NMAX )
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 af( i, j ) = 1.d0 / dble( i+j )
119 CALL dgeqrf( -1, 0, a, 1, b, w, 1, info )
120 CALL chkxer(
'DGEQRF', infot, nout, lerr, ok )
122 CALL dgeqrf( 0, -1, a, 1, b, w, 1, info )
123 CALL chkxer(
'DGEQRF', infot, nout, lerr, ok )
125 CALL dgeqrf( 2, 1, a, 1, b, w, 1, info )
126 CALL chkxer(
'DGEQRF', infot, nout, lerr, ok )
128 CALL dgeqrf( 1, 2, a, 1, b, w, 1, info )
129 CALL chkxer(
'DGEQRF', infot, nout, lerr, ok )
135 CALL dgeqrfp( -1, 0, a, 1, b, w, 1, info )
136 CALL chkxer(
'DGEQRFP', infot, nout, lerr, ok )
138 CALL dgeqrfp( 0, -1, a, 1, b, w, 1, info )
139 CALL chkxer(
'DGEQRFP', infot, nout, lerr, ok )
141 CALL dgeqrfp( 2, 1, a, 1, b, w, 1, info )
142 CALL chkxer(
'DGEQRFP', infot, nout, lerr, ok )
144 CALL dgeqrfp( 1, 2, a, 1, b, w, 1, info )
145 CALL chkxer(
'DGEQRFP', infot, nout, lerr, ok )
151 CALL dgeqr2( -1, 0, a, 1, b, w, info )
152 CALL chkxer(
'DGEQR2', infot, nout, lerr, ok )
154 CALL dgeqr2( 0, -1, a, 1, b, w, info )
155 CALL chkxer(
'DGEQR2', infot, nout, lerr, ok )
157 CALL dgeqr2( 2, 1, a, 1, b, w, info )
158 CALL chkxer(
'DGEQR2', infot, nout, lerr, ok )
164 CALL dgeqr2p( -1, 0, a, 1, b, w, info )
165 CALL chkxer(
'DGEQR2P', infot, nout, lerr, ok )
167 CALL dgeqr2p( 0, -1, a, 1, b, w, info )
168 CALL chkxer(
'DGEQR2P', infot, nout, lerr, ok )
170 CALL dgeqr2p( 2, 1, a, 1, b, w, info )
171 CALL chkxer(
'DGEQR2P', infot, nout, lerr, ok )
177 CALL dorgqr( -1, 0, 0, a, 1, x, w, 1, info )
178 CALL chkxer(
'DORGQR', infot, nout, lerr, ok )
180 CALL dorgqr( 0, -1, 0, a, 1, x, w, 1, info )
181 CALL chkxer(
'DORGQR', infot, nout, lerr, ok )
183 CALL dorgqr( 1, 2, 0, a, 1, x, w, 2, info )
184 CALL chkxer(
'DORGQR', infot, nout, lerr, ok )
186 CALL dorgqr( 0, 0, -1, a, 1, x, w, 1, info )
187 CALL chkxer(
'DORGQR', infot, nout, lerr, ok )
189 CALL dorgqr( 1, 1, 2, a, 1, x, w, 1, info )
190 CALL chkxer(
'DORGQR', infot, nout, lerr, ok )
192 CALL dorgqr( 2, 2, 0, a, 1, x, w, 2, info )
193 CALL chkxer(
'DORGQR', infot, nout, lerr, ok )
195 CALL dorgqr( 2, 2, 0, a, 2, x, w, 1, info )
196 CALL chkxer(
'DORGQR', infot, nout, lerr, ok )
202 CALL dorg2r( -1, 0, 0, a, 1, x, w, info )
203 CALL chkxer(
'DORG2R', infot, nout, lerr, ok )
205 CALL dorg2r( 0, -1, 0, a, 1, x, w, info )
206 CALL chkxer(
'DORG2R', infot, nout, lerr, ok )
208 CALL dorg2r( 1, 2, 0, a, 1, x, w, info )
209 CALL chkxer(
'DORG2R', infot, nout, lerr, ok )
211 CALL dorg2r( 0, 0, -1, a, 1, x, w, info )
212 CALL chkxer(
'DORG2R', infot, nout, lerr, ok )
214 CALL dorg2r( 2, 1, 2, a, 2, x, w, info )
215 CALL chkxer(
'DORG2R', infot, nout, lerr, ok )
217 CALL dorg2r( 2, 1, 0, a, 1, x, w, info )
218 CALL chkxer(
'DORG2R', infot, nout, lerr, ok )
224 CALL dormqr(
'/',
'N', 0, 0, 0, a, 1, x, af, 1, w, 1, info )
225 CALL chkxer(
'DORMQR', infot, nout, lerr, ok )
227 CALL dormqr(
'L',
'/', 0, 0, 0, a, 1, x, af, 1, w, 1, info )
228 CALL chkxer(
'DORMQR', infot, nout, lerr, ok )
230 CALL dormqr(
'L',
'N', -1, 0, 0, a, 1, x, af, 1, w, 1, info )
231 CALL chkxer(
'DORMQR', infot, nout, lerr, ok )
233 CALL dormqr(
'L',
'N', 0, -1, 0, a, 1, x, af, 1, w, 1, info )
234 CALL chkxer(
'DORMQR', infot, nout, lerr, ok )
236 CALL dormqr(
'L',
'N', 0, 0, -1, a, 1, x, af, 1, w, 1, info )
237 CALL chkxer(
'DORMQR', infot, nout, lerr, ok )
239 CALL dormqr(
'L',
'N', 0, 1, 1, a, 1, x, af, 1, w, 1, info )
240 CALL chkxer(
'DORMQR', infot, nout, lerr, ok )
242 CALL dormqr(
'R',
'N', 1, 0, 1, a, 1, x, af, 1, w, 1, info )
243 CALL chkxer(
'DORMQR', infot, nout, lerr, ok )
245 CALL dormqr(
'L',
'N', 2, 1, 0, a, 1, x, af, 2, w, 1, info )
246 CALL chkxer(
'DORMQR', infot, nout, lerr, ok )
248 CALL dormqr(
'R',
'N', 1, 2, 0, a, 1, x, af, 1, w, 1, info )
249 CALL chkxer(
'DORMQR', infot, nout, lerr, ok )
251 CALL dormqr(
'L',
'N', 2, 1, 0, a, 2, x, af, 1, w, 1, info )
252 CALL chkxer(
'DORMQR', infot, nout, lerr, ok )
254 CALL dormqr(
'L',
'N', 1, 2, 0, a, 1, x, af, 1, w, 1, info )
255 CALL chkxer(
'DORMQR', infot, nout, lerr, ok )
257 CALL dormqr(
'R',
'N', 2, 1, 0, a, 1, x, af, 2, w, 1, info )
258 CALL chkxer(
'DORMQR', infot, nout, lerr, ok )
264 CALL dorm2r(
'/',
'N', 0, 0, 0, a, 1, x, af, 1, w, info )
265 CALL chkxer(
'DORM2R', infot, nout, lerr, ok )
267 CALL dorm2r(
'L',
'/', 0, 0, 0, a, 1, x, af, 1, w, info )
268 CALL chkxer(
'DORM2R', infot, nout, lerr, ok )
270 CALL dorm2r(
'L',
'N', -1, 0, 0, a, 1, x, af, 1, w, info )
271 CALL chkxer(
'DORM2R', infot, nout, lerr, ok )
273 CALL dorm2r(
'L',
'N', 0, -1, 0, a, 1, x, af, 1, w, info )
274 CALL chkxer(
'DORM2R', infot, nout, lerr, ok )
276 CALL dorm2r(
'L',
'N', 0, 0, -1, a, 1, x, af, 1, w, info )
277 CALL chkxer(
'DORM2R', infot, nout, lerr, ok )
279 CALL dorm2r(
'L',
'N', 0, 1, 1, a, 1, x, af, 1, w, info )
280 CALL chkxer(
'DORM2R', infot, nout, lerr, ok )
282 CALL dorm2r(
'R',
'N', 1, 0, 1, a, 1, x, af, 1, w, info )
283 CALL chkxer(
'DORM2R', infot, nout, lerr, ok )
285 CALL dorm2r(
'L',
'N', 2, 1, 0, a, 1, x, af, 2, w, info )
286 CALL chkxer(
'DORM2R', infot, nout, lerr, ok )
288 CALL dorm2r(
'R',
'N', 1, 2, 0, a, 1, x, af, 1, w, info )
289 CALL chkxer(
'DORM2R', infot, nout, lerr, ok )
291 CALL dorm2r(
'L',
'N', 2, 1, 0, a, 2, x, af, 1, w, info )
292 CALL chkxer(
'DORM2R', infot, nout, lerr, ok )
296 CALL alaesm( path, ok, nout )
subroutine alaesm(path, ok, nout)
ALAESM
subroutine chkxer(srnamt, infot, nout, lerr, ok)
subroutine derrqr(path, nunit)
DERRQR
subroutine dgeqr2(m, n, a, lda, tau, work, info)
DGEQR2 computes the QR factorization of a general rectangular matrix using an unblocked algorithm.
subroutine dgeqr2p(m, n, a, lda, tau, work, info)
DGEQR2P computes the QR factorization of a general rectangular matrix with non-negative diagonal elem...
subroutine dgeqrf(m, n, a, lda, tau, work, lwork, info)
DGEQRF
subroutine dgeqrfp(m, n, a, lda, tau, work, lwork, info)
DGEQRFP
subroutine dorg2r(m, n, k, a, lda, tau, work, info)
DORG2R generates all or part of the orthogonal matrix Q from a QR factorization determined by sgeqrf ...
subroutine dorgqr(m, n, k, a, lda, tau, work, lwork, info)
DORGQR
subroutine dorm2r(side, trans, m, n, k, a, lda, tau, c, ldc, work, info)
DORM2R multiplies a general matrix by the orthogonal matrix from a QR factorization determined by sge...
subroutine dormqr(side, trans, m, n, k, a, lda, tau, c, ldc, work, lwork, info)
DORMQR