81 REAL a( nmax, nmax ), b( nmax, nmax ), s( nmax ),
98 COMMON / infoc / infot, nout, ok, lerr
99 COMMON / srnamc / srnamt
104 WRITE( nout, fmt = * )
112 IF(
lsamen( 2, c2,
'LS' ) )
THEN
120 CALL sgels(
'/', 0, 0, 0, a, 1, b, 1, w, 1, info )
121 CALL chkxer(
'SGELS ', infot, nout, lerr, ok )
123 CALL sgels(
'N', -1, 0, 0, a, 1, b, 1, w, 1, info )
124 CALL chkxer(
'SGELS ', infot, nout, lerr, ok )
126 CALL sgels(
'N', 0, -1, 0, a, 1, b, 1, w, 1, info )
127 CALL chkxer(
'SGELS ', infot, nout, lerr, ok )
129 CALL sgels(
'N', 0, 0, -1, a, 1, b, 1, w, 1, info )
130 CALL chkxer(
'SGELS ', infot, nout, lerr, ok )
132 CALL sgels(
'N', 2, 0, 0, a, 1, b, 2, w, 2, info )
133 CALL chkxer(
'SGELS ', infot, nout, lerr, ok )
135 CALL sgels(
'N', 2, 0, 0, a, 2, b, 1, w, 2, info )
136 CALL chkxer(
'SGELS ', infot, nout, lerr, ok )
138 CALL sgels(
'N', 1, 1, 0, a, 1, b, 1, w, 1, info )
139 CALL chkxer(
'SGELS ', infot, nout, lerr, ok )
145 CALL sgelss( -1, 0, 0, a, 1, b, 1, s, rcond, irnk, w, 1, info )
146 CALL chkxer(
'SGELSS', infot, nout, lerr, ok )
148 CALL sgelss( 0, -1, 0, a, 1, b, 1, s, rcond, irnk, w, 1, info )
149 CALL chkxer(
'SGELSS', infot, nout, lerr, ok )
151 CALL sgelss( 0, 0, -1, a, 1, b, 1, s, rcond, irnk, w, 1, info )
152 CALL chkxer(
'SGELSS', infot, nout, lerr, ok )
154 CALL sgelss( 2, 0, 0, a, 1, b, 2, s, rcond, irnk, w, 2, info )
155 CALL chkxer(
'SGELSS', infot, nout, lerr, ok )
157 CALL sgelss( 2, 0, 0, a, 2, b, 1, s, rcond, irnk, w, 2, info )
158 CALL chkxer(
'SGELSS', infot, nout, lerr, ok )
164 CALL sgelsx( -1, 0, 0, a, 1, b, 1, ip, rcond, irnk, w, info )
165 CALL chkxer(
'SGELSX', infot, nout, lerr, ok )
167 CALL sgelsx( 0, -1, 0, a, 1, b, 1, ip, rcond, irnk, w, info )
168 CALL chkxer(
'SGELSX', infot, nout, lerr, ok )
170 CALL sgelsx( 0, 0, -1, a, 1, b, 1, ip, rcond, irnk, w, info )
171 CALL chkxer(
'SGELSX', infot, nout, lerr, ok )
173 CALL sgelsx( 2, 0, 0, a, 1, b, 2, ip, rcond, irnk, w, info )
174 CALL chkxer(
'SGELSX', infot, nout, lerr, ok )
176 CALL sgelsx( 2, 0, 0, a, 2, b, 1, ip, rcond, irnk, w, info )
177 CALL chkxer(
'SGELSX', infot, nout, lerr, ok )
183 CALL sgelsy( -1, 0, 0, a, 1, b, 1, ip, rcond, irnk, w, 10,
185 CALL chkxer(
'SGELSY', infot, nout, lerr, ok )
187 CALL sgelsy( 0, -1, 0, a, 1, b, 1, ip, rcond, irnk, w, 10,
189 CALL chkxer(
'SGELSY', infot, nout, lerr, ok )
191 CALL sgelsy( 0, 0, -1, a, 1, b, 1, ip, rcond, irnk, w, 10,
193 CALL chkxer(
'SGELSY', infot, nout, lerr, ok )
195 CALL sgelsy( 2, 0, 0, a, 1, b, 2, ip, rcond, irnk, w, 10,
197 CALL chkxer(
'SGELSY', infot, nout, lerr, ok )
199 CALL sgelsy( 2, 0, 0, a, 2, b, 1, ip, rcond, irnk, w, 10,
201 CALL chkxer(
'SGELSY', infot, nout, lerr, ok )
203 CALL sgelsy( 2, 2, 1, a, 2, b, 2, ip, rcond, irnk, w, 1, info )
204 CALL chkxer(
'SGELSY', infot, nout, lerr, ok )
210 CALL sgelsd( -1, 0, 0, a, 1, b, 1, s, rcond, irnk, w, 10,
212 CALL chkxer(
'SGELSD', infot, nout, lerr, ok )
214 CALL sgelsd( 0, -1, 0, a, 1, b, 1, s, rcond, irnk, w, 10,
216 CALL chkxer(
'SGELSD', infot, nout, lerr, ok )
218 CALL sgelsd( 0, 0, -1, a, 1, b, 1, s, rcond, irnk, w, 10,
220 CALL chkxer(
'SGELSD', infot, nout, lerr, ok )
222 CALL sgelsd( 2, 0, 0, a, 1, b, 2, s, rcond, irnk, w, 10,
224 CALL chkxer(
'SGELSD', infot, nout, lerr, ok )
226 CALL sgelsd( 2, 0, 0, a, 2, b, 1, s, rcond, irnk, w, 10,
228 CALL chkxer(
'SGELSD', infot, nout, lerr, ok )
230 CALL sgelsd( 2, 2, 1, a, 2, b, 2, s, rcond, irnk, w, 1, ip,
232 CALL chkxer(
'SGELSD', infot, nout, lerr, ok )
237 CALL alaesm( path, ok, nout )
subroutine sgelsd(M, N, NRHS, A, LDA, B, LDB, S, RCOND, RANK, WORK, LWORK, IWORK, INFO)
SGELSD computes the minimum-norm solution to a linear least squares problem for GE matrices ...
subroutine sgelss(M, N, NRHS, A, LDA, B, LDB, S, RCOND, RANK, WORK, LWORK, INFO)
SGELSS solves overdetermined or underdetermined systems for GE matrices
subroutine sgels(TRANS, M, N, NRHS, A, LDA, B, LDB, WORK, LWORK, INFO)
SGELS solves overdetermined or underdetermined systems for GE matrices
subroutine chkxer(SRNAMT, INFOT, NOUT, LERR, OK)
subroutine sgelsy(M, N, NRHS, A, LDA, B, LDB, JPVT, RCOND, RANK, WORK, LWORK, INFO)
SGELSY solves overdetermined or underdetermined systems for GE matrices
logical function lsamen(N, CA, CB)
LSAMEN
subroutine sgelsx(M, N, NRHS, A, LDA, B, LDB, JPVT, RCOND, RANK, WORK, INFO)
SGELSX solves overdetermined or underdetermined systems for GE matrices
subroutine alaesm(PATH, OK, NOUT)
ALAESM