56 SUBROUTINE serrlq( PATH, NUNIT )
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 sgelqf( -1, 0, a, 1, b, w, 1, info )
122 CALL chkxer(
'SGELQF', infot, nout, lerr, ok )
124 CALL sgelqf( 0, -1, a, 1, b, w, 1, info )
125 CALL chkxer(
'SGELQF', infot, nout, lerr, ok )
127 CALL sgelqf( 2, 1, a, 1, b, w, 2, info )
128 CALL chkxer(
'SGELQF', infot, nout, lerr, ok )
130 CALL sgelqf( 2, 1, a, 2, b, w, 1, info )
131 CALL chkxer(
'SGELQF', infot, nout, lerr, ok )
137 CALL sgelq2( -1, 0, a, 1, b, w, info )
138 CALL chkxer(
'SGELQ2', infot, nout, lerr, ok )
140 CALL sgelq2( 0, -1, a, 1, b, w, info )
141 CALL chkxer(
'SGELQ2', infot, nout, lerr, ok )
143 CALL sgelq2( 2, 1, a, 1, b, w, info )
144 CALL chkxer(
'SGELQ2', infot, nout, lerr, ok )
150 CALL sgelqs( -1, 0, 0, a, 1, x, b, 1, w, 1, info )
151 CALL chkxer(
'SGELQS', infot, nout, lerr, ok )
153 CALL sgelqs( 0, -1, 0, a, 1, x, b, 1, w, 1, info )
154 CALL chkxer(
'SGELQS', infot, nout, lerr, ok )
156 CALL sgelqs( 2, 1, 0, a, 2, x, b, 1, w, 1, info )
157 CALL chkxer(
'SGELQS', infot, nout, lerr, ok )
159 CALL sgelqs( 0, 0, -1, a, 1, x, b, 1, w, 1, info )
160 CALL chkxer(
'SGELQS', infot, nout, lerr, ok )
162 CALL sgelqs( 2, 2, 0, a, 1, x, b, 2, w, 1, info )
163 CALL chkxer(
'SGELQS', infot, nout, lerr, ok )
165 CALL sgelqs( 1, 2, 0, a, 1, x, b, 1, w, 1, info )
166 CALL chkxer(
'SGELQS', infot, nout, lerr, ok )
168 CALL sgelqs( 1, 1, 2, a, 1, x, b, 1, w, 1, info )
169 CALL chkxer(
'SGELQS', infot, nout, lerr, ok )
175 CALL sorglq( -1, 0, 0, a, 1, x, w, 1, info )
176 CALL chkxer(
'SORGLQ', infot, nout, lerr, ok )
178 CALL sorglq( 0, -1, 0, a, 1, x, w, 1, info )
179 CALL chkxer(
'SORGLQ', infot, nout, lerr, ok )
181 CALL sorglq( 2, 1, 0, a, 2, x, w, 2, info )
182 CALL chkxer(
'SORGLQ', infot, nout, lerr, ok )
184 CALL sorglq( 0, 0, -1, a, 1, x, w, 1, info )
185 CALL chkxer(
'SORGLQ', infot, nout, lerr, ok )
187 CALL sorglq( 1, 1, 2, a, 1, x, w, 1, info )
188 CALL chkxer(
'SORGLQ', infot, nout, lerr, ok )
190 CALL sorglq( 2, 2, 0, a, 1, x, w, 2, info )
191 CALL chkxer(
'SORGLQ', infot, nout, lerr, ok )
193 CALL sorglq( 2, 2, 0, a, 2, x, w, 1, info )
194 CALL chkxer(
'SORGLQ', infot, nout, lerr, ok )
200 CALL sorgl2( -1, 0, 0, a, 1, x, w, info )
201 CALL chkxer(
'SORGL2', infot, nout, lerr, ok )
203 CALL sorgl2( 0, -1, 0, a, 1, x, w, info )
204 CALL chkxer(
'SORGL2', infot, nout, lerr, ok )
206 CALL sorgl2( 2, 1, 0, a, 2, x, w, info )
207 CALL chkxer(
'SORGL2', infot, nout, lerr, ok )
209 CALL sorgl2( 0, 0, -1, a, 1, x, w, info )
210 CALL chkxer(
'SORGL2', infot, nout, lerr, ok )
212 CALL sorgl2( 1, 1, 2, a, 1, x, w, info )
213 CALL chkxer(
'SORGL2', infot, nout, lerr, ok )
215 CALL sorgl2( 2, 2, 0, a, 1, x, w, info )
216 CALL chkxer(
'SORGL2', infot, nout, lerr, ok )
222 CALL sormlq(
'/',
'N', 0, 0, 0, a, 1, x, af, 1, w, 1, info )
223 CALL chkxer(
'SORMLQ', infot, nout, lerr, ok )
225 CALL sormlq(
'L',
'/', 0, 0, 0, a, 1, x, af, 1, w, 1, info )
226 CALL chkxer(
'SORMLQ', infot, nout, lerr, ok )
228 CALL sormlq(
'L',
'N', -1, 0, 0, a, 1, x, af, 1, w, 1, info )
229 CALL chkxer(
'SORMLQ', infot, nout, lerr, ok )
231 CALL sormlq(
'L',
'N', 0, -1, 0, a, 1, x, af, 1, w, 1, info )
232 CALL chkxer(
'SORMLQ', infot, nout, lerr, ok )
234 CALL sormlq(
'L',
'N', 0, 0, -1, a, 1, x, af, 1, w, 1, info )
235 CALL chkxer(
'SORMLQ', infot, nout, lerr, ok )
237 CALL sormlq(
'L',
'N', 0, 1, 1, a, 1, x, af, 1, w, 1, info )
238 CALL chkxer(
'SORMLQ', infot, nout, lerr, ok )
240 CALL sormlq(
'R',
'N', 1, 0, 1, a, 1, x, af, 1, w, 1, info )
241 CALL chkxer(
'SORMLQ', infot, nout, lerr, ok )
243 CALL sormlq(
'L',
'N', 2, 0, 2, a, 1, x, af, 2, w, 1, info )
244 CALL chkxer(
'SORMLQ', infot, nout, lerr, ok )
246 CALL sormlq(
'R',
'N', 0, 2, 2, a, 1, x, af, 1, w, 1, info )
247 CALL chkxer(
'SORMLQ', infot, nout, lerr, ok )
249 CALL sormlq(
'L',
'N', 2, 1, 0, a, 2, x, af, 1, w, 1, info )
250 CALL chkxer(
'SORMLQ', infot, nout, lerr, ok )
252 CALL sormlq(
'L',
'N', 1, 2, 0, a, 1, x, af, 1, w, 1, info )
253 CALL chkxer(
'SORMLQ', infot, nout, lerr, ok )
255 CALL sormlq(
'R',
'N', 2, 1, 0, a, 1, x, af, 2, w, 1, info )
256 CALL chkxer(
'SORMLQ', infot, nout, lerr, ok )
262 CALL sorml2(
'/',
'N', 0, 0, 0, a, 1, x, af, 1, w, info )
263 CALL chkxer(
'SORML2', infot, nout, lerr, ok )
265 CALL sorml2(
'L',
'/', 0, 0, 0, a, 1, x, af, 1, w, info )
266 CALL chkxer(
'SORML2', infot, nout, lerr, ok )
268 CALL sorml2(
'L',
'N', -1, 0, 0, a, 1, x, af, 1, w, info )
269 CALL chkxer(
'SORML2', infot, nout, lerr, ok )
271 CALL sorml2(
'L',
'N', 0, -1, 0, a, 1, x, af, 1, w, info )
272 CALL chkxer(
'SORML2', infot, nout, lerr, ok )
274 CALL sorml2(
'L',
'N', 0, 0, -1, a, 1, x, af, 1, w, info )
275 CALL chkxer(
'SORML2', infot, nout, lerr, ok )
277 CALL sorml2(
'L',
'N', 0, 1, 1, a, 1, x, af, 1, w, info )
278 CALL chkxer(
'SORML2', infot, nout, lerr, ok )
280 CALL sorml2(
'R',
'N', 1, 0, 1, a, 1, x, af, 1, w, info )
281 CALL chkxer(
'SORML2', infot, nout, lerr, ok )
283 CALL sorml2(
'L',
'N', 2, 1, 2, a, 1, x, af, 2, w, info )
284 CALL chkxer(
'SORML2', infot, nout, lerr, ok )
286 CALL sorml2(
'R',
'N', 1, 2, 2, a, 1, x, af, 1, w, info )
287 CALL chkxer(
'SORML2', infot, nout, lerr, ok )
289 CALL sorml2(
'L',
'N', 2, 1, 0, a, 2, x, af, 1, w, info )
290 CALL chkxer(
'SORML2', infot, nout, lerr, ok )
294 CALL alaesm( path, ok, nout )
subroutine sgelqs(M, N, NRHS, A, LDA, TAU, B, LDB, WORK, LWORK, INFO)
SGELQS
subroutine alaesm(PATH, OK, NOUT)
ALAESM
subroutine sorglq(M, N, K, A, LDA, TAU, WORK, LWORK, INFO)
SORGLQ
subroutine sorml2(SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC, WORK, INFO)
SORML2 multiplies a general matrix by the orthogonal matrix from a LQ factorization determined by sge...
subroutine sorgl2(M, N, K, A, LDA, TAU, WORK, INFO)
SORGL2
subroutine chkxer(SRNAMT, INFOT, NOUT, LERR, OK)
subroutine sormlq(SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC, WORK, LWORK, INFO)
SORMLQ
subroutine sgelq2(M, N, A, LDA, TAU, WORK, INFO)
SGELQ2 computes the LQ factorization of a general rectangular matrix using an unblocked algorithm...
subroutine serrlq(PATH, NUNIT)
SERRLQ
subroutine sgelqf(M, N, A, LDA, TAU, WORK, LWORK, INFO)
SGELQF