56 SUBROUTINE derrls( PATH, NUNIT )
77 DOUBLE PRECISION RCOND
81 DOUBLE PRECISION 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 dgels(
'/', 0, 0, 0, a, 1, b, 1, w, 1, info )
121 CALL chkxer(
'DGELS ', infot, nout, lerr, ok )
123 CALL dgels(
'N', -1, 0, 0, a, 1, b, 1, w, 1, info )
124 CALL chkxer(
'DGELS ', infot, nout, lerr, ok )
126 CALL dgels(
'N', 0, -1, 0, a, 1, b, 1, w, 1, info )
127 CALL chkxer(
'DGELS ', infot, nout, lerr, ok )
129 CALL dgels(
'N', 0, 0, -1, a, 1, b, 1, w, 1, info )
130 CALL chkxer(
'DGELS ', infot, nout, lerr, ok )
132 CALL dgels(
'N', 2, 0, 0, a, 1, b, 2, w, 2, info )
133 CALL chkxer(
'DGELS ', infot, nout, lerr, ok )
135 CALL dgels(
'N', 2, 0, 0, a, 2, b, 1, w, 2, info )
136 CALL chkxer(
'DGELS ', infot, nout, lerr, ok )
138 CALL dgels(
'N', 1, 1, 0, a, 1, b, 1, w, 1, info )
139 CALL chkxer(
'DGELS ', infot, nout, lerr, ok )
145 CALL dgelss( -1, 0, 0, a, 1, b, 1, s, rcond, irnk, w, 1, info )
146 CALL chkxer(
'DGELSS', infot, nout, lerr, ok )
148 CALL dgelss( 0, -1, 0, a, 1, b, 1, s, rcond, irnk, w, 1, info )
149 CALL chkxer(
'DGELSS', infot, nout, lerr, ok )
151 CALL dgelss( 0, 0, -1, a, 1, b, 1, s, rcond, irnk, w, 1, info )
152 CALL chkxer(
'DGELSS', infot, nout, lerr, ok )
154 CALL dgelss( 2, 0, 0, a, 1, b, 2, s, rcond, irnk, w, 2, info )
155 CALL chkxer(
'DGELSS', infot, nout, lerr, ok )
157 CALL dgelss( 2, 0, 0, a, 2, b, 1, s, rcond, irnk, w, 2, info )
158 CALL chkxer(
'DGELSS', infot, nout, lerr, ok )
164 CALL dgelsx( -1, 0, 0, a, 1, b, 1, ip, rcond, irnk, w, info )
165 CALL chkxer(
'DGELSX', infot, nout, lerr, ok )
167 CALL dgelsx( 0, -1, 0, a, 1, b, 1, ip, rcond, irnk, w, info )
168 CALL chkxer(
'DGELSX', infot, nout, lerr, ok )
170 CALL dgelsx( 0, 0, -1, a, 1, b, 1, ip, rcond, irnk, w, info )
171 CALL chkxer(
'DGELSX', infot, nout, lerr, ok )
173 CALL dgelsx( 2, 0, 0, a, 1, b, 2, ip, rcond, irnk, w, info )
174 CALL chkxer(
'DGELSX', infot, nout, lerr, ok )
176 CALL dgelsx( 2, 0, 0, a, 2, b, 1, ip, rcond, irnk, w, info )
177 CALL chkxer(
'DGELSX', infot, nout, lerr, ok )
183 CALL dgelsy( -1, 0, 0, a, 1, b, 1, ip, rcond, irnk, w, 10,
185 CALL chkxer(
'DGELSY', infot, nout, lerr, ok )
187 CALL dgelsy( 0, -1, 0, a, 1, b, 1, ip, rcond, irnk, w, 10,
189 CALL chkxer(
'DGELSY', infot, nout, lerr, ok )
191 CALL dgelsy( 0, 0, -1, a, 1, b, 1, ip, rcond, irnk, w, 10,
193 CALL chkxer(
'DGELSY', infot, nout, lerr, ok )
195 CALL dgelsy( 2, 0, 0, a, 1, b, 2, ip, rcond, irnk, w, 10,
197 CALL chkxer(
'DGELSY', infot, nout, lerr, ok )
199 CALL dgelsy( 2, 0, 0, a, 2, b, 1, ip, rcond, irnk, w, 10,
201 CALL chkxer(
'DGELSY', infot, nout, lerr, ok )
203 CALL dgelsy( 2, 2, 1, a, 2, b, 2, ip, rcond, irnk, w, 1, info )
204 CALL chkxer(
'DGELSY', infot, nout, lerr, ok )
210 CALL dgelsd( -1, 0, 0, a, 1, b, 1, s, rcond, irnk, w, 10, ip,
212 CALL chkxer(
'DGELSD', infot, nout, lerr, ok )
214 CALL dgelsd( 0, -1, 0, a, 1, b, 1, s, rcond, irnk, w, 10, ip,
216 CALL chkxer(
'DGELSD', infot, nout, lerr, ok )
218 CALL dgelsd( 0, 0, -1, a, 1, b, 1, s, rcond, irnk, w, 10, ip,
220 CALL chkxer(
'DGELSD', infot, nout, lerr, ok )
222 CALL dgelsd( 2, 0, 0, a, 1, b, 2, s, rcond, irnk, w, 10, ip,
224 CALL chkxer(
'DGELSD', infot, nout, lerr, ok )
226 CALL dgelsd( 2, 0, 0, a, 2, b, 1, s, rcond, irnk, w, 10, ip,
228 CALL chkxer(
'DGELSD', infot, nout, lerr, ok )
230 CALL dgelsd( 2, 2, 1, a, 2, b, 2, s, rcond, irnk, w, 1, ip,
232 CALL chkxer(
'DGELSD', infot, nout, lerr, ok )
237 CALL alaesm( path, ok, nout )
subroutine dgelsx(M, N, NRHS, A, LDA, B, LDB, JPVT, RCOND, RANK, WORK, INFO)
DGELSX solves overdetermined or underdetermined systems for GE matrices
subroutine dgels(TRANS, M, N, NRHS, A, LDA, B, LDB, WORK, LWORK, INFO)
DGELS solves overdetermined or underdetermined systems for GE matrices
subroutine chkxer(SRNAMT, INFOT, NOUT, LERR, OK)
subroutine dgelsy(M, N, NRHS, A, LDA, B, LDB, JPVT, RCOND, RANK, WORK, LWORK, INFO)
DGELSY solves overdetermined or underdetermined systems for GE matrices
subroutine dgelss(M, N, NRHS, A, LDA, B, LDB, S, RCOND, RANK, WORK, LWORK, INFO)
DGELSS solves overdetermined or underdetermined systems for GE matrices
subroutine derrls(PATH, NUNIT)
DERRLS
subroutine alaesm(PATH, OK, NOUT)
ALAESM
subroutine dgelsd(M, N, NRHS, A, LDA, B, LDB, S, RCOND, RANK, WORK, LWORK, IWORK, INFO)
DGELSD computes the minimum-norm solution to a linear least squares problem for GE matrices ...