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)
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 zgels(TRANS, M, N, NRHS, A, LDA, B, LDB, WORK, LWORK, INFO)
ZGELS solves overdetermined or underdetermined systems for GE matrices