78 REAL A( NMAX, NMAX ), B( NMAX, NMAX ), S( NMAX ),
94 COMMON / infoc / infot, nout, ok, lerr
95 COMMON / srnamc / srnamt
100 WRITE( nout, fmt = * )
108 IF(
lsamen( 2, c2,
'LS' ) )
THEN
116 CALL sgels(
'/', 0, 0, 0, a, 1, b, 1, w, 1, info )
117 CALL chkxer(
'SGELS ', infot, nout, lerr, ok )
119 CALL sgels(
'N', -1, 0, 0, a, 1, b, 1, w, 1, info )
120 CALL chkxer(
'SGELS ', infot, nout, lerr, ok )
122 CALL sgels(
'N', 0, -1, 0, a, 1, b, 1, w, 1, info )
123 CALL chkxer(
'SGELS ', infot, nout, lerr, ok )
125 CALL sgels(
'N', 0, 0, -1, a, 1, b, 1, w, 1, info )
126 CALL chkxer(
'SGELS ', infot, nout, lerr, ok )
128 CALL sgels(
'N', 2, 0, 0, a, 1, b, 2, w, 2, info )
129 CALL chkxer(
'SGELS ', infot, nout, lerr, ok )
131 CALL sgels(
'N', 2, 0, 0, a, 2, b, 1, w, 2, info )
132 CALL chkxer(
'SGELS ', infot, nout, lerr, ok )
134 CALL sgels(
'N', 1, 1, 0, a, 1, b, 1, w, 1, info )
135 CALL chkxer(
'SGELS ', infot, nout, lerr, ok )
141 CALL sgelss( -1, 0, 0, a, 1, b, 1, s, rcond, irnk, w, 1, info )
142 CALL chkxer(
'SGELSS', infot, nout, lerr, ok )
144 CALL sgelss( 0, -1, 0, a, 1, b, 1, s, rcond, irnk, w, 1, info )
145 CALL chkxer(
'SGELSS', infot, nout, lerr, ok )
147 CALL sgelss( 0, 0, -1, a, 1, b, 1, s, rcond, irnk, w, 1, info )
148 CALL chkxer(
'SGELSS', infot, nout, lerr, ok )
150 CALL sgelss( 2, 0, 0, a, 1, b, 2, s, rcond, irnk, w, 2, info )
151 CALL chkxer(
'SGELSS', infot, nout, lerr, ok )
153 CALL sgelss( 2, 0, 0, a, 2, b, 1, s, rcond, irnk, w, 2, info )
154 CALL chkxer(
'SGELSS', infot, nout, lerr, ok )
160 CALL sgelsy( -1, 0, 0, a, 1, b, 1, ip, rcond, irnk, w, 10,
162 CALL chkxer(
'SGELSY', infot, nout, lerr, ok )
164 CALL sgelsy( 0, -1, 0, a, 1, b, 1, ip, rcond, irnk, w, 10,
166 CALL chkxer(
'SGELSY', infot, nout, lerr, ok )
168 CALL sgelsy( 0, 0, -1, a, 1, b, 1, ip, rcond, irnk, w, 10,
170 CALL chkxer(
'SGELSY', infot, nout, lerr, ok )
172 CALL sgelsy( 2, 0, 0, a, 1, b, 2, ip, rcond, irnk, w, 10,
174 CALL chkxer(
'SGELSY', infot, nout, lerr, ok )
176 CALL sgelsy( 2, 0, 0, a, 2, b, 1, ip, rcond, irnk, w, 10,
178 CALL chkxer(
'SGELSY', infot, nout, lerr, ok )
180 CALL sgelsy( 2, 2, 1, a, 2, b, 2, ip, rcond, irnk, w, 1, info )
181 CALL chkxer(
'SGELSY', infot, nout, lerr, ok )
187 CALL sgelsd( -1, 0, 0, a, 1, b, 1, s, rcond, irnk, w, 10,
189 CALL chkxer(
'SGELSD', infot, nout, lerr, ok )
191 CALL sgelsd( 0, -1, 0, a, 1, b, 1, s, rcond, irnk, w, 10,
193 CALL chkxer(
'SGELSD', infot, nout, lerr, ok )
195 CALL sgelsd( 0, 0, -1, a, 1, b, 1, s, rcond, irnk, w, 10,
197 CALL chkxer(
'SGELSD', infot, nout, lerr, ok )
199 CALL sgelsd( 2, 0, 0, a, 1, b, 2, s, rcond, irnk, w, 10,
201 CALL chkxer(
'SGELSD', infot, nout, lerr, ok )
203 CALL sgelsd( 2, 0, 0, a, 2, b, 1, s, rcond, irnk, w, 10,
205 CALL chkxer(
'SGELSD', infot, nout, lerr, ok )
207 CALL sgelsd( 2, 2, 1, a, 2, b, 2, s, rcond, irnk, w, 1, ip,
209 CALL chkxer(
'SGELSD', infot, nout, lerr, ok )
214 CALL alaesm( path, ok, nout )
subroutine chkxer(SRNAMT, INFOT, NOUT, LERR, OK)
logical function lsamen(N, CA, CB)
LSAMEN
subroutine alaesm(PATH, OK, NOUT)
ALAESM
subroutine sgels(TRANS, M, N, NRHS, A, LDA, B, LDB, WORK, LWORK, INFO)
SGELS solves overdetermined or underdetermined systems 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 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 sgelsy(M, N, NRHS, A, LDA, B, LDB, JPVT, RCOND, RANK, WORK, LWORK, INFO)
SGELSY solves overdetermined or underdetermined systems for GE matrices