77 DOUBLE PRECISION rcond
81 DOUBLE PRECISION rw( nmax ), s( nmax )
82 COMPLEX*16 a( nmax, nmax ),
b( nmax, nmax ), w( nmax )
98 COMMON / infoc / infot, nout, ok, lerr
99 COMMON / srnamc / srnamt
105 a( 1, 1 ) = ( 1.0d+0, 0.0d+0 )
106 a( 1, 2 ) = ( 2.0d+0, 0.0d+0 )
107 a( 2, 2 ) = ( 3.0d+0, 0.0d+0 )
108 a( 2, 1 ) = ( 4.0d+0, 0.0d+0 )
110 WRITE( nout, fmt = * )
114 IF(
lsamen( 2, c2,
'LS' ) )
THEN
120 CALL
zgels(
'/', 0, 0, 0, a, 1,
b, 1, w, 1, info )
121 CALL
chkxer(
'ZGELS ', infot, nout, lerr, ok )
123 CALL
zgels(
'N', -1, 0, 0, a, 1,
b, 1, w, 1, info )
124 CALL
chkxer(
'ZGELS ', infot, nout, lerr, ok )
126 CALL
zgels(
'N', 0, -1, 0, a, 1,
b, 1, w, 1, info )
127 CALL
chkxer(
'ZGELS ', infot, nout, lerr, ok )
129 CALL
zgels(
'N', 0, 0, -1, a, 1,
b, 1, w, 1, info )
130 CALL
chkxer(
'ZGELS ', infot, nout, lerr, ok )
132 CALL
zgels(
'N', 2, 0, 0, a, 1,
b, 2, w, 2, info )
133 CALL
chkxer(
'ZGELS ', infot, nout, lerr, ok )
135 CALL
zgels(
'N', 2, 0, 0, a, 2,
b, 1, w, 2, info )
136 CALL
chkxer(
'ZGELS ', infot, nout, lerr, ok )
138 CALL
zgels(
'N', 1, 1, 0, a, 1,
b, 1, w, 1, info )
139 CALL
chkxer(
'ZGELS ', infot, nout, lerr, ok )
145 CALL
zgelss( -1, 0, 0, a, 1,
b, 1, s, rcond, irnk, w, 1, rw,
147 CALL
chkxer(
'ZGELSS', infot, nout, lerr, ok )
149 CALL
zgelss( 0, -1, 0, a, 1,
b, 1, s, rcond, irnk, w, 1, rw,
151 CALL
chkxer(
'ZGELSS', infot, nout, lerr, ok )
153 CALL
zgelss( 0, 0, -1, a, 1,
b, 1, s, rcond, irnk, w, 1, rw,
155 CALL
chkxer(
'ZGELSS', infot, nout, lerr, ok )
157 CALL
zgelss( 2, 0, 0, a, 1,
b, 2, s, rcond, irnk, w, 2, rw,
159 CALL
chkxer(
'ZGELSS', infot, nout, lerr, ok )
161 CALL
zgelss( 2, 0, 0, a, 2,
b, 1, s, rcond, irnk, w, 2, rw,
163 CALL
chkxer(
'ZGELSS', infot, nout, lerr, ok )
169 CALL
zgelsx( -1, 0, 0, a, 1,
b, 1, ip, rcond, irnk, w, rw,
171 CALL
chkxer(
'ZGELSX', infot, nout, lerr, ok )
173 CALL
zgelsx( 0, -1, 0, a, 1,
b, 1, ip, rcond, irnk, w, rw,
175 CALL
chkxer(
'ZGELSX', infot, nout, lerr, ok )
177 CALL
zgelsx( 0, 0, -1, a, 1,
b, 1, ip, rcond, irnk, w, rw,
179 CALL
chkxer(
'ZGELSX', infot, nout, lerr, ok )
181 CALL
zgelsx( 2, 0, 0, a, 1,
b, 2, ip, rcond, irnk, w, rw,
183 CALL
chkxer(
'ZGELSX', infot, nout, lerr, ok )
185 CALL
zgelsx( 2, 0, 0, a, 2,
b, 1, ip, rcond, irnk, w, rw,
187 CALL
chkxer(
'ZGELSX', infot, nout, lerr, ok )
193 CALL
zgelsy( -1, 0, 0, a, 1,
b, 1, ip, rcond, irnk, w, 10, rw,
195 CALL
chkxer(
'ZGELSY', infot, nout, lerr, ok )
197 CALL
zgelsy( 0, -1, 0, a, 1,
b, 1, ip, rcond, irnk, w, 10, rw,
199 CALL
chkxer(
'ZGELSY', infot, nout, lerr, ok )
201 CALL
zgelsy( 0, 0, -1, a, 1,
b, 1, ip, rcond, irnk, w, 10, rw,
203 CALL
chkxer(
'ZGELSY', infot, nout, lerr, ok )
205 CALL
zgelsy( 2, 0, 0, a, 1,
b, 2, ip, rcond, irnk, w, 10, rw,
207 CALL
chkxer(
'ZGELSY', infot, nout, lerr, ok )
209 CALL
zgelsy( 2, 0, 0, a, 2,
b, 1, ip, rcond, irnk, w, 10, rw,
211 CALL
chkxer(
'ZGELSY', infot, nout, lerr, ok )
213 CALL
zgelsy( 0, 3, 0, a, 1,
b, 3, ip, rcond, irnk, w, 1, rw,
215 CALL
chkxer(
'ZGELSY', infot, nout, lerr, ok )
221 CALL
zgelsd( -1, 0, 0, a, 1,
b, 1, s, rcond, irnk, w, 10, rw,
223 CALL
chkxer(
'ZGELSD', infot, nout, lerr, ok )
225 CALL
zgelsd( 0, -1, 0, a, 1,
b, 1, s, rcond, irnk, w, 10, rw,
227 CALL
chkxer(
'ZGELSD', infot, nout, lerr, ok )
229 CALL
zgelsd( 0, 0, -1, a, 1,
b, 1, s, rcond, irnk, w, 10, rw,
231 CALL
chkxer(
'ZGELSD', infot, nout, lerr, ok )
233 CALL
zgelsd( 2, 0, 0, a, 1,
b, 2, s, rcond, irnk, w, 10, rw,
235 CALL
chkxer(
'ZGELSD', infot, nout, lerr, ok )
237 CALL
zgelsd( 2, 0, 0, a, 2,
b, 1, s, rcond, irnk, w, 10, rw,
239 CALL
chkxer(
'ZGELSD', infot, nout, lerr, ok )
241 CALL
zgelsd( 2, 2, 1, a, 2,
b, 2, s, rcond, irnk, w, 1, rw, ip,
243 CALL
chkxer(
'ZGELSD', infot, nout, lerr, ok )
248 CALL
alaesm( path, ok, nout )
subroutine zgelsx(M, N, NRHS, A, LDA, B, LDB, JPVT, RCOND, RANK, WORK, RWORK, INFO)
ZGELSX solves overdetermined or underdetermined systems for GE matrices
subroutine zgelsy(M, N, NRHS, A, LDA, B, LDB, JPVT, RCOND, RANK, WORK, LWORK, RWORK, INFO)
ZGELSY solves overdetermined or underdetermined systems for GE matrices
subroutine chkxer(SRNAMT, INFOT, NOUT, LERR, OK)
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real b(3) integer i
subroutine zgelss(M, N, NRHS, A, LDA, B, LDB, S, RCOND, RANK, WORK, LWORK, RWORK, INFO)
ZGELSS solves overdetermined or underdetermined systems for GE matrices
logical function lsamen(N, CA, CB)
LSAMEN
subroutine zgelsd(M, N, NRHS, A, LDA, B, LDB, S, RCOND, RANK, WORK, LWORK, RWORK, IWORK, INFO)
ZGELSD computes the minimum-norm solution to a linear least squares problem for GE matrices ...
subroutine alaesm(PATH, OK, NOUT)
ALAESM
subroutine zerrls(PATH, NUNIT)
ZERRLS
subroutine zgels(TRANS, M, N, NRHS, A, LDA, B, LDB, WORK, LWORK, INFO)
ZGELS solves overdetermined or underdetermined systems for GE matrices