75 REAL 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. / real( i+j )
105 af( i, j ) = 1. / real( i+j )
119 CALL sgeqrf( -1, 0, a, 1, b, w, 1, info )
120 CALL chkxer(
'SGEQRF', infot, nout, lerr, ok )
122 CALL sgeqrf( 0, -1, a, 1, b, w, 1, info )
123 CALL chkxer(
'SGEQRF', infot, nout, lerr, ok )
125 CALL sgeqrf( 2, 1, a, 1, b, w, 1, info )
126 CALL chkxer(
'SGEQRF', infot, nout, lerr, ok )
128 CALL sgeqrf( 1, 2, a, 1, b, w, 1, info )
129 CALL chkxer(
'SGEQRF', infot, nout, lerr, ok )
135 CALL sgeqrfp( -1, 0, a, 1, b, w, 1, info )
136 CALL chkxer(
'SGEQRFP', infot, nout, lerr, ok )
138 CALL sgeqrfp( 0, -1, a, 1, b, w, 1, info )
139 CALL chkxer(
'SGEQRFP', infot, nout, lerr, ok )
141 CALL sgeqrfp( 2, 1, a, 1, b, w, 1, info )
142 CALL chkxer(
'SGEQRFP', infot, nout, lerr, ok )
144 CALL sgeqrfp( 1, 2, a, 1, b, w, 1, info )
145 CALL chkxer(
'SGEQRFP', infot, nout, lerr, ok )
151 CALL sgeqr2( -1, 0, a, 1, b, w, info )
152 CALL chkxer(
'SGEQR2', infot, nout, lerr, ok )
154 CALL sgeqr2( 0, -1, a, 1, b, w, info )
155 CALL chkxer(
'SGEQR2', infot, nout, lerr, ok )
157 CALL sgeqr2( 2, 1, a, 1, b, w, info )
158 CALL chkxer(
'SGEQR2', infot, nout, lerr, ok )
164 CALL sgeqr2p( -1, 0, a, 1, b, w, info )
165 CALL chkxer(
'SGEQR2P', infot, nout, lerr, ok )
167 CALL sgeqr2p( 0, -1, a, 1, b, w, info )
168 CALL chkxer(
'SGEQR2P', infot, nout, lerr, ok )
170 CALL sgeqr2p( 2, 1, a, 1, b, w, info )
171 CALL chkxer(
'SGEQR2P', infot, nout, lerr, ok )
177 CALL sorgqr( -1, 0, 0, a, 1, x, w, 1, info )
178 CALL chkxer(
'SORGQR', infot, nout, lerr, ok )
180 CALL sorgqr( 0, -1, 0, a, 1, x, w, 1, info )
181 CALL chkxer(
'SORGQR', infot, nout, lerr, ok )
183 CALL sorgqr( 1, 2, 0, a, 1, x, w, 2, info )
184 CALL chkxer(
'SORGQR', infot, nout, lerr, ok )
186 CALL sorgqr( 0, 0, -1, a, 1, x, w, 1, info )
187 CALL chkxer(
'SORGQR', infot, nout, lerr, ok )
189 CALL sorgqr( 1, 1, 2, a, 1, x, w, 1, info )
190 CALL chkxer(
'SORGQR', infot, nout, lerr, ok )
192 CALL sorgqr( 2, 2, 0, a, 1, x, w, 2, info )
193 CALL chkxer(
'SORGQR', infot, nout, lerr, ok )
195 CALL sorgqr( 2, 2, 0, a, 2, x, w, 1, info )
196 CALL chkxer(
'SORGQR', infot, nout, lerr, ok )
202 CALL sorg2r( -1, 0, 0, a, 1, x, w, info )
203 CALL chkxer(
'SORG2R', infot, nout, lerr, ok )
205 CALL sorg2r( 0, -1, 0, a, 1, x, w, info )
206 CALL chkxer(
'SORG2R', infot, nout, lerr, ok )
208 CALL sorg2r( 1, 2, 0, a, 1, x, w, info )
209 CALL chkxer(
'SORG2R', infot, nout, lerr, ok )
211 CALL sorg2r( 0, 0, -1, a, 1, x, w, info )
212 CALL chkxer(
'SORG2R', infot, nout, lerr, ok )
214 CALL sorg2r( 2, 1, 2, a, 2, x, w, info )
215 CALL chkxer(
'SORG2R', infot, nout, lerr, ok )
217 CALL sorg2r( 2, 1, 0, a, 1, x, w, info )
218 CALL chkxer(
'SORG2R', infot, nout, lerr, ok )
224 CALL sormqr(
'/',
'N', 0, 0, 0, a, 1, x, af, 1, w, 1, info )
225 CALL chkxer(
'SORMQR', infot, nout, lerr, ok )
227 CALL sormqr(
'L',
'/', 0, 0, 0, a, 1, x, af, 1, w, 1, info )
228 CALL chkxer(
'SORMQR', infot, nout, lerr, ok )
230 CALL sormqr(
'L',
'N', -1, 0, 0, a, 1, x, af, 1, w, 1, info )
231 CALL chkxer(
'SORMQR', infot, nout, lerr, ok )
233 CALL sormqr(
'L',
'N', 0, -1, 0, a, 1, x, af, 1, w, 1, info )
234 CALL chkxer(
'SORMQR', infot, nout, lerr, ok )
236 CALL sormqr(
'L',
'N', 0, 0, -1, a, 1, x, af, 1, w, 1, info )
237 CALL chkxer(
'SORMQR', infot, nout, lerr, ok )
239 CALL sormqr(
'L',
'N', 0, 1, 1, a, 1, x, af, 1, w, 1, info )
240 CALL chkxer(
'SORMQR', infot, nout, lerr, ok )
242 CALL sormqr(
'R',
'N', 1, 0, 1, a, 1, x, af, 1, w, 1, info )
243 CALL chkxer(
'SORMQR', infot, nout, lerr, ok )
245 CALL sormqr(
'L',
'N', 2, 1, 0, a, 1, x, af, 2, w, 1, info )
246 CALL chkxer(
'SORMQR', infot, nout, lerr, ok )
248 CALL sormqr(
'R',
'N', 1, 2, 0, a, 1, x, af, 1, w, 1, info )
249 CALL chkxer(
'SORMQR', infot, nout, lerr, ok )
251 CALL sormqr(
'L',
'N', 2, 1, 0, a, 2, x, af, 1, w, 1, info )
252 CALL chkxer(
'SORMQR', infot, nout, lerr, ok )
254 CALL sormqr(
'L',
'N', 1, 2, 0, a, 1, x, af, 1, w, 1, info )
255 CALL chkxer(
'SORMQR', infot, nout, lerr, ok )
257 CALL sormqr(
'R',
'N', 2, 1, 0, a, 1, x, af, 2, w, 1, info )
258 CALL chkxer(
'SORMQR', infot, nout, lerr, ok )
264 CALL sorm2r(
'/',
'N', 0, 0, 0, a, 1, x, af, 1, w, info )
265 CALL chkxer(
'SORM2R', infot, nout, lerr, ok )
267 CALL sorm2r(
'L',
'/', 0, 0, 0, a, 1, x, af, 1, w, info )
268 CALL chkxer(
'SORM2R', infot, nout, lerr, ok )
270 CALL sorm2r(
'L',
'N', -1, 0, 0, a, 1, x, af, 1, w, info )
271 CALL chkxer(
'SORM2R', infot, nout, lerr, ok )
273 CALL sorm2r(
'L',
'N', 0, -1, 0, a, 1, x, af, 1, w, info )
274 CALL chkxer(
'SORM2R', infot, nout, lerr, ok )
276 CALL sorm2r(
'L',
'N', 0, 0, -1, a, 1, x, af, 1, w, info )
277 CALL chkxer(
'SORM2R', infot, nout, lerr, ok )
279 CALL sorm2r(
'L',
'N', 0, 1, 1, a, 1, x, af, 1, w, info )
280 CALL chkxer(
'SORM2R', infot, nout, lerr, ok )
282 CALL sorm2r(
'R',
'N', 1, 0, 1, a, 1, x, af, 1, w, info )
283 CALL chkxer(
'SORM2R', infot, nout, lerr, ok )
285 CALL sorm2r(
'L',
'N', 2, 1, 0, a, 1, x, af, 2, w, info )
286 CALL chkxer(
'SORM2R', infot, nout, lerr, ok )
288 CALL sorm2r(
'R',
'N', 1, 2, 0, a, 1, x, af, 1, w, info )
289 CALL chkxer(
'SORM2R', infot, nout, lerr, ok )
291 CALL sorm2r(
'L',
'N', 2, 1, 0, a, 2, x, af, 1, w, info )
292 CALL chkxer(
'SORM2R', infot, nout, lerr, ok )
296 CALL alaesm( path, ok, nout )
subroutine alaesm(path, ok, nout)
ALAESM
subroutine chkxer(srnamt, infot, nout, lerr, ok)
subroutine sgeqr2(m, n, a, lda, tau, work, info)
SGEQR2 computes the QR factorization of a general rectangular matrix using an unblocked algorithm.
subroutine sgeqr2p(m, n, a, lda, tau, work, info)
SGEQR2P computes the QR factorization of a general rectangular matrix with non-negative diagonal elem...
subroutine sgeqrf(m, n, a, lda, tau, work, lwork, info)
SGEQRF
subroutine sgeqrfp(m, n, a, lda, tau, work, lwork, info)
SGEQRFP
subroutine sorg2r(m, n, k, a, lda, tau, work, info)
SORG2R generates all or part of the orthogonal matrix Q from a QR factorization determined by sgeqrf ...
subroutine sorgqr(m, n, k, a, lda, tau, work, lwork, info)
SORGQR
subroutine sorm2r(side, trans, m, n, k, a, lda, tau, c, ldc, work, info)
SORM2R multiplies a general matrix by the orthogonal matrix from a QR factorization determined by sge...
subroutine sormqr(side, trans, m, n, k, a, lda, tau, c, ldc, work, lwork, info)
SORMQR
subroutine serrqr(path, nunit)
SERRQR