74 DOUBLE PRECISION RCOND
78 DOUBLE PRECISION 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 dgels(
'/', 0, 0, 0, a, 1, b, 1, w, 1, info )
117 CALL chkxer(
'DGELS ', infot, nout, lerr, ok )
119 CALL dgels(
'N', -1, 0, 0, a, 1, b, 1, w, 1, info )
120 CALL chkxer(
'DGELS ', infot, nout, lerr, ok )
122 CALL dgels(
'N', 0, -1, 0, a, 1, b, 1, w, 1, info )
123 CALL chkxer(
'DGELS ', infot, nout, lerr, ok )
125 CALL dgels(
'N', 0, 0, -1, a, 1, b, 1, w, 1, info )
126 CALL chkxer(
'DGELS ', infot, nout, lerr, ok )
128 CALL dgels(
'N', 2, 0, 0, a, 1, b, 2, w, 2, info )
129 CALL chkxer(
'DGELS ', infot, nout, lerr, ok )
131 CALL dgels(
'N', 2, 0, 0, a, 2, b, 1, w, 2, info )
132 CALL chkxer(
'DGELS ', infot, nout, lerr, ok )
134 CALL dgels(
'N', 1, 1, 0, a, 1, b, 1, w, 1, info )
135 CALL chkxer(
'DGELS ', infot, nout, lerr, ok )
141 CALL dgelss( -1, 0, 0, a, 1, b, 1, s, rcond, irnk, w, 1, info )
142 CALL chkxer(
'DGELSS', infot, nout, lerr, ok )
144 CALL dgelss( 0, -1, 0, a, 1, b, 1, s, rcond, irnk, w, 1, info )
145 CALL chkxer(
'DGELSS', infot, nout, lerr, ok )
147 CALL dgelss( 0, 0, -1, a, 1, b, 1, s, rcond, irnk, w, 1, info )
148 CALL chkxer(
'DGELSS', infot, nout, lerr, ok )
150 CALL dgelss( 2, 0, 0, a, 1, b, 2, s, rcond, irnk, w, 2, info )
151 CALL chkxer(
'DGELSS', infot, nout, lerr, ok )
153 CALL dgelss( 2, 0, 0, a, 2, b, 1, s, rcond, irnk, w, 2, info )
154 CALL chkxer(
'DGELSS', infot, nout, lerr, ok )
160 CALL dgelsy( -1, 0, 0, a, 1, b, 1, ip, rcond, irnk, w, 10,
162 CALL chkxer(
'DGELSY', infot, nout, lerr, ok )
164 CALL dgelsy( 0, -1, 0, a, 1, b, 1, ip, rcond, irnk, w, 10,
166 CALL chkxer(
'DGELSY', infot, nout, lerr, ok )
168 CALL dgelsy( 0, 0, -1, a, 1, b, 1, ip, rcond, irnk, w, 10,
170 CALL chkxer(
'DGELSY', infot, nout, lerr, ok )
172 CALL dgelsy( 2, 0, 0, a, 1, b, 2, ip, rcond, irnk, w, 10,
174 CALL chkxer(
'DGELSY', infot, nout, lerr, ok )
176 CALL dgelsy( 2, 0, 0, a, 2, b, 1, ip, rcond, irnk, w, 10,
178 CALL chkxer(
'DGELSY', infot, nout, lerr, ok )
180 CALL dgelsy( 2, 2, 1, a, 2, b, 2, ip, rcond, irnk, w, 1, info )
181 CALL chkxer(
'DGELSY', infot, nout, lerr, ok )
187 CALL dgelsd( -1, 0, 0, a, 1, b, 1, s, rcond, irnk, w, 10, ip,
189 CALL chkxer(
'DGELSD', infot, nout, lerr, ok )
191 CALL dgelsd( 0, -1, 0, a, 1, b, 1, s, rcond, irnk, w, 10, ip,
193 CALL chkxer(
'DGELSD', infot, nout, lerr, ok )
195 CALL dgelsd( 0, 0, -1, a, 1, b, 1, s, rcond, irnk, w, 10, ip,
197 CALL chkxer(
'DGELSD', infot, nout, lerr, ok )
199 CALL dgelsd( 2, 0, 0, a, 1, b, 2, s, rcond, irnk, w, 10, ip,
201 CALL chkxer(
'DGELSD', infot, nout, lerr, ok )
203 CALL dgelsd( 2, 0, 0, a, 2, b, 1, s, rcond, irnk, w, 10, ip,
205 CALL chkxer(
'DGELSD', infot, nout, lerr, ok )
207 CALL dgelsd( 2, 2, 1, a, 2, b, 2, s, rcond, irnk, w, 1, ip,
209 CALL chkxer(
'DGELSD', infot, nout, lerr, ok )
214 CALL alaesm( path, ok, nout )
subroutine chkxer(SRNAMT, INFOT, NOUT, LERR, OK)
subroutine alaesm(PATH, OK, NOUT)
ALAESM
subroutine derrls(PATH, NUNIT)
DERRLS
subroutine dgels(TRANS, M, N, NRHS, A, LDA, B, LDB, WORK, LWORK, INFO)
DGELS solves overdetermined or underdetermined systems for GE matrices
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 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
subroutine dgelss(M, N, NRHS, A, LDA, B, LDB, S, RCOND, RANK, WORK, LWORK, INFO)
DGELSS solves overdetermined or underdetermined systems for GE matrices