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
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real b(3) integer i
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
LOGICAL function lsamen(N, CA, CB)
LSAMEN
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 ...