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
sgeqlf( -1, 0, a, 1,
b, w, 1, info )
122 CALL
chkxer(
'SGEQLF', infot, nout, lerr, ok )
124 CALL
sgeqlf( 0, -1, a, 1,
b, w, 1, info )
125 CALL
chkxer(
'SGEQLF', infot, nout, lerr, ok )
127 CALL
sgeqlf( 2, 1, a, 1,
b, w, 1, info )
128 CALL
chkxer(
'SGEQLF', infot, nout, lerr, ok )
130 CALL
sgeqlf( 1, 2, a, 1,
b, w, 1, info )
131 CALL
chkxer(
'SGEQLF', infot, nout, lerr, ok )
137 CALL
sgeql2( -1, 0, a, 1,
b, w, info )
138 CALL
chkxer(
'SGEQL2', infot, nout, lerr, ok )
140 CALL
sgeql2( 0, -1, a, 1,
b, w, info )
141 CALL
chkxer(
'SGEQL2', infot, nout, lerr, ok )
143 CALL
sgeql2( 2, 1, a, 1,
b, w, info )
144 CALL
chkxer(
'SGEQL2', infot, nout, lerr, ok )
150 CALL
sgeqls( -1, 0, 0, a, 1,
x,
b, 1, w, 1, info )
151 CALL
chkxer(
'SGEQLS', infot, nout, lerr, ok )
153 CALL
sgeqls( 0, -1, 0, a, 1,
x,
b, 1, w, 1, info )
154 CALL
chkxer(
'SGEQLS', infot, nout, lerr, ok )
156 CALL
sgeqls( 1, 2, 0, a, 1,
x,
b, 1, w, 1, info )
157 CALL
chkxer(
'SGEQLS', infot, nout, lerr, ok )
159 CALL
sgeqls( 0, 0, -1, a, 1,
x,
b, 1, w, 1, info )
160 CALL
chkxer(
'SGEQLS', infot, nout, lerr, ok )
162 CALL
sgeqls( 2, 1, 0, a, 1,
x,
b, 2, w, 1, info )
163 CALL
chkxer(
'SGEQLS', infot, nout, lerr, ok )
165 CALL
sgeqls( 2, 1, 0, a, 2,
x,
b, 1, w, 1, info )
166 CALL
chkxer(
'SGEQLS', infot, nout, lerr, ok )
168 CALL
sgeqls( 1, 1, 2, a, 1,
x,
b, 1, w, 1, info )
169 CALL
chkxer(
'SGEQLS', infot, nout, lerr, ok )
175 CALL
sorgql( -1, 0, 0, a, 1,
x, w, 1, info )
176 CALL
chkxer(
'SORGQL', infot, nout, lerr, ok )
178 CALL
sorgql( 0, -1, 0, a, 1,
x, w, 1, info )
179 CALL
chkxer(
'SORGQL', infot, nout, lerr, ok )
181 CALL
sorgql( 1, 2, 0, a, 1,
x, w, 2, info )
182 CALL
chkxer(
'SORGQL', infot, nout, lerr, ok )
184 CALL
sorgql( 0, 0, -1, a, 1,
x, w, 1, info )
185 CALL
chkxer(
'SORGQL', infot, nout, lerr, ok )
187 CALL
sorgql( 1, 1, 2, a, 1,
x, w, 1, info )
188 CALL
chkxer(
'SORGQL', infot, nout, lerr, ok )
190 CALL
sorgql( 2, 1, 0, a, 1,
x, w, 1, info )
191 CALL
chkxer(
'SORGQL', infot, nout, lerr, ok )
193 CALL
sorgql( 2, 2, 0, a, 2,
x, w, 1, info )
194 CALL
chkxer(
'SORGQL', infot, nout, lerr, ok )
200 CALL
sorg2l( -1, 0, 0, a, 1,
x, w, info )
201 CALL
chkxer(
'SORG2L', infot, nout, lerr, ok )
203 CALL
sorg2l( 0, -1, 0, a, 1,
x, w, info )
204 CALL
chkxer(
'SORG2L', infot, nout, lerr, ok )
206 CALL
sorg2l( 1, 2, 0, a, 1,
x, w, info )
207 CALL
chkxer(
'SORG2L', infot, nout, lerr, ok )
209 CALL
sorg2l( 0, 0, -1, a, 1,
x, w, info )
210 CALL
chkxer(
'SORG2L', infot, nout, lerr, ok )
212 CALL
sorg2l( 2, 1, 2, a, 2,
x, w, info )
213 CALL
chkxer(
'SORG2L', infot, nout, lerr, ok )
215 CALL
sorg2l( 2, 1, 0, a, 1,
x, w, info )
216 CALL
chkxer(
'SORG2L', infot, nout, lerr, ok )
222 CALL
sormql(
'/',
'N', 0, 0, 0, a, 1,
x, af, 1, w, 1, info )
223 CALL
chkxer(
'SORMQL', infot, nout, lerr, ok )
225 CALL
sormql(
'L',
'/', 0, 0, 0, a, 1,
x, af, 1, w, 1, info )
226 CALL
chkxer(
'SORMQL', infot, nout, lerr, ok )
228 CALL
sormql(
'L',
'N', -1, 0, 0, a, 1,
x, af, 1, w, 1, info )
229 CALL
chkxer(
'SORMQL', infot, nout, lerr, ok )
231 CALL
sormql(
'L',
'N', 0, -1, 0, a, 1,
x, af, 1, w, 1, info )
232 CALL
chkxer(
'SORMQL', infot, nout, lerr, ok )
234 CALL
sormql(
'L',
'N', 0, 0, -1, a, 1,
x, af, 1, w, 1, info )
235 CALL
chkxer(
'SORMQL', infot, nout, lerr, ok )
237 CALL
sormql(
'L',
'N', 0, 1, 1, a, 1,
x, af, 1, w, 1, info )
238 CALL
chkxer(
'SORMQL', infot, nout, lerr, ok )
240 CALL
sormql(
'R',
'N', 1, 0, 1, a, 1,
x, af, 1, w, 1, info )
241 CALL
chkxer(
'SORMQL', infot, nout, lerr, ok )
243 CALL
sormql(
'L',
'N', 2, 1, 0, a, 1,
x, af, 2, w, 1, info )
244 CALL
chkxer(
'SORMQL', infot, nout, lerr, ok )
246 CALL
sormql(
'R',
'N', 1, 2, 0, a, 1,
x, af, 1, w, 1, info )
247 CALL
chkxer(
'SORMQL', infot, nout, lerr, ok )
249 CALL
sormql(
'L',
'N', 2, 1, 0, a, 2,
x, af, 1, w, 1, info )
250 CALL
chkxer(
'SORMQL', infot, nout, lerr, ok )
252 CALL
sormql(
'L',
'N', 1, 2, 0, a, 1,
x, af, 1, w, 1, info )
253 CALL
chkxer(
'SORMQL', infot, nout, lerr, ok )
255 CALL
sormql(
'R',
'N', 2, 1, 0, a, 1,
x, af, 2, w, 1, info )
256 CALL
chkxer(
'SORMQL', infot, nout, lerr, ok )
262 CALL
sorm2l(
'/',
'N', 0, 0, 0, a, 1,
x, af, 1, w, info )
263 CALL
chkxer(
'SORM2L', infot, nout, lerr, ok )
265 CALL
sorm2l(
'L',
'/', 0, 0, 0, a, 1,
x, af, 1, w, info )
266 CALL
chkxer(
'SORM2L', infot, nout, lerr, ok )
268 CALL
sorm2l(
'L',
'N', -1, 0, 0, a, 1,
x, af, 1, w, info )
269 CALL
chkxer(
'SORM2L', infot, nout, lerr, ok )
271 CALL
sorm2l(
'L',
'N', 0, -1, 0, a, 1,
x, af, 1, w, info )
272 CALL
chkxer(
'SORM2L', infot, nout, lerr, ok )
274 CALL
sorm2l(
'L',
'N', 0, 0, -1, a, 1,
x, af, 1, w, info )
275 CALL
chkxer(
'SORM2L', infot, nout, lerr, ok )
277 CALL
sorm2l(
'L',
'N', 0, 1, 1, a, 1,
x, af, 1, w, info )
278 CALL
chkxer(
'SORM2L', infot, nout, lerr, ok )
280 CALL
sorm2l(
'R',
'N', 1, 0, 1, a, 1,
x, af, 1, w, info )
281 CALL
chkxer(
'SORM2L', infot, nout, lerr, ok )
283 CALL
sorm2l(
'L',
'N', 2, 1, 0, a, 1,
x, af, 2, w, info )
284 CALL
chkxer(
'SORM2L', infot, nout, lerr, ok )
286 CALL
sorm2l(
'R',
'N', 1, 2, 0, a, 1,
x, af, 1, w, info )
287 CALL
chkxer(
'SORM2L', infot, nout, lerr, ok )
289 CALL
sorm2l(
'L',
'N', 2, 1, 0, a, 2,
x, af, 1, w, info )
290 CALL
chkxer(
'SORM2L', infot, nout, lerr, ok )
294 CALL
alaesm( path, ok, nout )
subroutine sormql(SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC, WORK, LWORK, INFO)
SORMQL
subroutine sorm2l(SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC, WORK, INFO)
SORM2L multiplies a general matrix by the orthogonal matrix from a QL factorization determined by sge...
subroutine chkxer(SRNAMT, INFOT, NOUT, LERR, OK)
subroutine sgeqls(M, N, NRHS, A, LDA, TAU, B, LDB, WORK, LWORK, INFO)
SGEQLS
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real b(3) integer i
subroutine sgeqlf(M, N, A, LDA, TAU, WORK, LWORK, INFO)
SGEQLF
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j
subroutine serrql(PATH, NUNIT)
SERRQL
subroutine sgeql2(M, N, A, LDA, TAU, WORK, INFO)
SGEQL2 computes the QL factorization of a general rectangular matrix using an unblocked algorithm...
subroutine sorgql(M, N, K, A, LDA, TAU, WORK, LWORK, INFO)
SORGQL
subroutine alaesm(PATH, OK, NOUT)
ALAESM
subroutine sorg2l(M, N, K, A, LDA, TAU, WORK, INFO)
SORG2L generates all or part of the orthogonal matrix Q from a QL factorization determined by sgeqlf ...