72 parameter ( nmax = 2 )
78 REAL a( nmax, nmax ), af( nmax, nmax ), b( nmax ),
79 $ w( nmax ), x( nmax )
91 COMMON / infoc / infot, nout, ok, lerr
92 COMMON / srnamc / srnamt
100 WRITE( nout, fmt = * )
106 a( i, j ) = 1. /
REAL( i+j )
107 af( i, j ) = 1. /
REAL( i+j )
121 CALL sgerqf( -1, 0, a, 1, b, w, 1, info )
122 CALL chkxer(
'SGERQF', infot, nout, lerr, ok )
124 CALL sgerqf( 0, -1, a, 1, b, w, 1, info )
125 CALL chkxer(
'SGERQF', infot, nout, lerr, ok )
127 CALL sgerqf( 2, 1, a, 1, b, w, 2, info )
128 CALL chkxer(
'SGERQF', infot, nout, lerr, ok )
130 CALL sgerqf( 2, 1, a, 2, b, w, 1, info )
131 CALL chkxer(
'SGERQF', infot, nout, lerr, ok )
137 CALL sgerq2( -1, 0, a, 1, b, w, info )
138 CALL chkxer(
'SGERQ2', infot, nout, lerr, ok )
140 CALL sgerq2( 0, -1, a, 1, b, w, info )
141 CALL chkxer(
'SGERQ2', infot, nout, lerr, ok )
143 CALL sgerq2( 2, 1, a, 1, b, w, info )
144 CALL chkxer(
'SGERQ2', infot, nout, lerr, ok )
150 CALL sgerqs( -1, 0, 0, a, 1, x, b, 1, w, 1, info )
151 CALL chkxer(
'SGERQS', infot, nout, lerr, ok )
153 CALL sgerqs( 0, -1, 0, a, 1, x, b, 1, w, 1, info )
154 CALL chkxer(
'SGERQS', infot, nout, lerr, ok )
156 CALL sgerqs( 2, 1, 0, a, 2, x, b, 1, w, 1, info )
157 CALL chkxer(
'SGERQS', infot, nout, lerr, ok )
159 CALL sgerqs( 0, 0, -1, a, 1, x, b, 1, w, 1, info )
160 CALL chkxer(
'SGERQS', infot, nout, lerr, ok )
162 CALL sgerqs( 2, 2, 0, a, 1, x, b, 2, w, 1, info )
163 CALL chkxer(
'SGERQS', infot, nout, lerr, ok )
165 CALL sgerqs( 2, 2, 0, a, 2, x, b, 1, w, 1, info )
166 CALL chkxer(
'SGERQS', infot, nout, lerr, ok )
168 CALL sgerqs( 1, 1, 2, a, 1, x, b, 1, w, 1, info )
169 CALL chkxer(
'SGERQS', infot, nout, lerr, ok )
175 CALL sorgrq( -1, 0, 0, a, 1, x, w, 1, info )
176 CALL chkxer(
'SORGRQ', infot, nout, lerr, ok )
178 CALL sorgrq( 0, -1, 0, a, 1, x, w, 1, info )
179 CALL chkxer(
'SORGRQ', infot, nout, lerr, ok )
181 CALL sorgrq( 2, 1, 0, a, 2, x, w, 2, info )
182 CALL chkxer(
'SORGRQ', infot, nout, lerr, ok )
184 CALL sorgrq( 0, 0, -1, a, 1, x, w, 1, info )
185 CALL chkxer(
'SORGRQ', infot, nout, lerr, ok )
187 CALL sorgrq( 1, 2, 2, a, 1, x, w, 1, info )
188 CALL chkxer(
'SORGRQ', infot, nout, lerr, ok )
190 CALL sorgrq( 2, 2, 0, a, 1, x, w, 2, info )
191 CALL chkxer(
'SORGRQ', infot, nout, lerr, ok )
193 CALL sorgrq( 2, 2, 0, a, 2, x, w, 1, info )
194 CALL chkxer(
'SORGRQ', infot, nout, lerr, ok )
200 CALL sorgr2( -1, 0, 0, a, 1, x, w, info )
201 CALL chkxer(
'SORGR2', infot, nout, lerr, ok )
203 CALL sorgr2( 0, -1, 0, a, 1, x, w, info )
204 CALL chkxer(
'SORGR2', infot, nout, lerr, ok )
206 CALL sorgr2( 2, 1, 0, a, 2, x, w, info )
207 CALL chkxer(
'SORGR2', infot, nout, lerr, ok )
209 CALL sorgr2( 0, 0, -1, a, 1, x, w, info )
210 CALL chkxer(
'SORGR2', infot, nout, lerr, ok )
212 CALL sorgr2( 1, 2, 2, a, 2, x, w, info )
213 CALL chkxer(
'SORGR2', infot, nout, lerr, ok )
215 CALL sorgr2( 2, 2, 0, a, 1, x, w, info )
216 CALL chkxer(
'SORGR2', infot, nout, lerr, ok )
222 CALL sormrq(
'/',
'N', 0, 0, 0, a, 1, x, af, 1, w, 1, info )
223 CALL chkxer(
'SORMRQ', infot, nout, lerr, ok )
225 CALL sormrq(
'L',
'/', 0, 0, 0, a, 1, x, af, 1, w, 1, info )
226 CALL chkxer(
'SORMRQ', infot, nout, lerr, ok )
228 CALL sormrq(
'L',
'N', -1, 0, 0, a, 1, x, af, 1, w, 1, info )
229 CALL chkxer(
'SORMRQ', infot, nout, lerr, ok )
231 CALL sormrq(
'L',
'N', 0, -1, 0, a, 1, x, af, 1, w, 1, info )
232 CALL chkxer(
'SORMRQ', infot, nout, lerr, ok )
234 CALL sormrq(
'L',
'N', 0, 0, -1, a, 1, x, af, 1, w, 1, info )
235 CALL chkxer(
'SORMRQ', infot, nout, lerr, ok )
237 CALL sormrq(
'L',
'N', 0, 1, 1, a, 1, x, af, 1, w, 1, info )
238 CALL chkxer(
'SORMRQ', infot, nout, lerr, ok )
240 CALL sormrq(
'R',
'N', 1, 0, 1, a, 1, x, af, 1, w, 1, info )
241 CALL chkxer(
'SORMRQ', infot, nout, lerr, ok )
243 CALL sormrq(
'L',
'N', 2, 1, 2, a, 1, x, af, 2, w, 1, info )
244 CALL chkxer(
'SORMRQ', infot, nout, lerr, ok )
246 CALL sormrq(
'R',
'N', 1, 2, 2, a, 1, x, af, 1, w, 1, info )
247 CALL chkxer(
'SORMRQ', infot, nout, lerr, ok )
249 CALL sormrq(
'L',
'N', 2, 1, 0, a, 1, x, af, 1, w, 1, info )
250 CALL chkxer(
'SORMRQ', infot, nout, lerr, ok )
252 CALL sormrq(
'L',
'N', 1, 2, 0, a, 1, x, af, 1, w, 1, info )
253 CALL chkxer(
'SORMRQ', infot, nout, lerr, ok )
255 CALL sormrq(
'R',
'N', 2, 1, 0, a, 1, x, af, 2, w, 1, info )
256 CALL chkxer(
'SORMRQ', infot, nout, lerr, ok )
262 CALL sormr2(
'/',
'N', 0, 0, 0, a, 1, x, af, 1, w, info )
263 CALL chkxer(
'SORMR2', infot, nout, lerr, ok )
265 CALL sormr2(
'L',
'/', 0, 0, 0, a, 1, x, af, 1, w, info )
266 CALL chkxer(
'SORMR2', infot, nout, lerr, ok )
268 CALL sormr2(
'L',
'N', -1, 0, 0, a, 1, x, af, 1, w, info )
269 CALL chkxer(
'SORMR2', infot, nout, lerr, ok )
271 CALL sormr2(
'L',
'N', 0, -1, 0, a, 1, x, af, 1, w, info )
272 CALL chkxer(
'SORMR2', infot, nout, lerr, ok )
274 CALL sormr2(
'L',
'N', 0, 0, -1, a, 1, x, af, 1, w, info )
275 CALL chkxer(
'SORMR2', infot, nout, lerr, ok )
277 CALL sormr2(
'L',
'N', 0, 1, 1, a, 1, x, af, 1, w, info )
278 CALL chkxer(
'SORMR2', infot, nout, lerr, ok )
280 CALL sormr2(
'R',
'N', 1, 0, 1, a, 1, x, af, 1, w, info )
281 CALL chkxer(
'SORMR2', infot, nout, lerr, ok )
283 CALL sormr2(
'L',
'N', 2, 1, 2, a, 1, x, af, 2, w, info )
284 CALL chkxer(
'SORMR2', infot, nout, lerr, ok )
286 CALL sormr2(
'R',
'N', 1, 2, 2, a, 1, x, af, 1, w, info )
287 CALL chkxer(
'SORMR2', infot, nout, lerr, ok )
289 CALL sormr2(
'L',
'N', 2, 1, 0, a, 1, x, af, 1, w, info )
290 CALL chkxer(
'SORMR2', infot, nout, lerr, ok )
294 CALL alaesm( path, ok, nout )
subroutine sorgr2(M, N, K, A, LDA, TAU, WORK, INFO)
SORGR2 generates all or part of the orthogonal matrix Q from an RQ factorization determined by sgerqf...
subroutine alaesm(PATH, OK, NOUT)
ALAESM
subroutine sgerqs(M, N, NRHS, A, LDA, TAU, B, LDB, WORK, LWORK, INFO)
SGERQS
subroutine chkxer(SRNAMT, INFOT, NOUT, LERR, OK)
subroutine sormrq(SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC, WORK, LWORK, INFO)
SORMRQ
subroutine sgerq2(M, N, A, LDA, TAU, WORK, INFO)
SGERQ2 computes the RQ factorization of a general rectangular matrix using an unblocked algorithm...
subroutine sgerqf(M, N, A, LDA, TAU, WORK, LWORK, INFO)
SGERQF
subroutine sorgrq(M, N, K, A, LDA, TAU, WORK, LWORK, INFO)
SORGRQ
subroutine sormr2(SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC, WORK, INFO)
SORMR2 multiplies a general matrix by the orthogonal matrix from a RQ factorization determined by sge...