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 sormlq(SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC, WORK, LWORK, INFO)
SORMLQ
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 chkxer(SRNAMT, INFOT, NOUT, LERR, OK)
subroutine sgelqs(M, N, NRHS, A, LDA, TAU, B, LDB, WORK, LWORK, INFO)
SGELQS
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real b(3) integer i
subroutine sorgl2(M, N, K, A, LDA, TAU, WORK, INFO)
SORGL2
subroutine sgelqf(M, N, A, LDA, TAU, WORK, LWORK, INFO)
SGELQF
subroutine sgelq2(M, N, A, LDA, TAU, WORK, INFO)
SGELQ2 computes the LQ factorization of a general rectangular matrix using an unblocked algorithm...
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j
subroutine alaesm(PATH, OK, NOUT)
ALAESM
subroutine sorglq(M, N, K, A, LDA, TAU, WORK, LWORK, INFO)
SORGLQ
subroutine serrlq(PATH, NUNIT)
SERRLQ