74 DOUBLE PRECISION RCOND
78 DOUBLE PRECISION RW( NMAX ), S( NMAX )
79 COMPLEX*16 A( NMAX, NMAX ), B( NMAX, NMAX ), W( NMAX )
94 COMMON / infoc / infot, nout, ok, lerr
95 COMMON / srnamc / srnamt
101 a( 1, 1 ) = ( 1.0d+0, 0.0d+0 )
102 a( 1, 2 ) = ( 2.0d+0, 0.0d+0 )
103 a( 2, 2 ) = ( 3.0d+0, 0.0d+0 )
104 a( 2, 1 ) = ( 4.0d+0, 0.0d+0 )
106 WRITE( nout, fmt = * )
110 IF( lsamen( 2, c2,
'LS' ) )
THEN
116 CALL zgels(
'/', 0, 0, 0, a, 1, b, 1, w, 1, info )
117 CALL chkxer(
'ZGELS ', infot, nout, lerr, ok )
119 CALL zgels(
'N', -1, 0, 0, a, 1, b, 1, w, 1, info )
120 CALL chkxer(
'ZGELS ', infot, nout, lerr, ok )
122 CALL zgels(
'N', 0, -1, 0, a, 1, b, 1, w, 1, info )
123 CALL chkxer(
'ZGELS ', infot, nout, lerr, ok )
125 CALL zgels(
'N', 0, 0, -1, a, 1, b, 1, w, 1, info )
126 CALL chkxer(
'ZGELS ', infot, nout, lerr, ok )
128 CALL zgels(
'N', 2, 0, 0, a, 1, b, 2, w, 2, info )
129 CALL chkxer(
'ZGELS ', infot, nout, lerr, ok )
131 CALL zgels(
'N', 2, 0, 0, a, 2, b, 1, w, 2, info )
132 CALL chkxer(
'ZGELS ', infot, nout, lerr, ok )
134 CALL zgels(
'N', 1, 1, 0, a, 1, b, 1, w, 1, info )
135 CALL chkxer(
'ZGELS ', infot, nout, lerr, ok )
141 CALL zgelss( -1, 0, 0, a, 1, b, 1, s, rcond, irnk, w, 1, rw,
143 CALL chkxer(
'ZGELSS', infot, nout, lerr, ok )
145 CALL zgelss( 0, -1, 0, a, 1, b, 1, s, rcond, irnk, w, 1, rw,
147 CALL chkxer(
'ZGELSS', infot, nout, lerr, ok )
149 CALL zgelss( 0, 0, -1, a, 1, b, 1, s, rcond, irnk, w, 1, rw,
151 CALL chkxer(
'ZGELSS', infot, nout, lerr, ok )
153 CALL zgelss( 2, 0, 0, a, 1, b, 2, s, rcond, irnk, w, 2, rw,
155 CALL chkxer(
'ZGELSS', infot, nout, lerr, ok )
157 CALL zgelss( 2, 0, 0, a, 2, b, 1, s, rcond, irnk, w, 2, rw,
159 CALL chkxer(
'ZGELSS', infot, nout, lerr, ok )
165 CALL zgelsy( -1, 0, 0, a, 1, b, 1, ip, rcond, irnk, w, 10, rw,
167 CALL chkxer(
'ZGELSY', infot, nout, lerr, ok )
169 CALL zgelsy( 0, -1, 0, a, 1, b, 1, ip, rcond, irnk, w, 10, rw,
171 CALL chkxer(
'ZGELSY', infot, nout, lerr, ok )
173 CALL zgelsy( 0, 0, -1, a, 1, b, 1, ip, rcond, irnk, w, 10, rw,
175 CALL chkxer(
'ZGELSY', infot, nout, lerr, ok )
177 CALL zgelsy( 2, 0, 0, a, 1, b, 2, ip, rcond, irnk, w, 10, rw,
179 CALL chkxer(
'ZGELSY', infot, nout, lerr, ok )
181 CALL zgelsy( 2, 0, 0, a, 2, b, 1, ip, rcond, irnk, w, 10, rw,
183 CALL chkxer(
'ZGELSY', infot, nout, lerr, ok )
185 CALL zgelsy( 0, 3, 0, a, 1, b, 3, ip, rcond, irnk, w, 1, rw,
187 CALL chkxer(
'ZGELSY', infot, nout, lerr, ok )
193 CALL zgelsd( -1, 0, 0, a, 1, b, 1, s, rcond, irnk, w, 10, rw,
195 CALL chkxer(
'ZGELSD', infot, nout, lerr, ok )
197 CALL zgelsd( 0, -1, 0, a, 1, b, 1, s, rcond, irnk, w, 10, rw,
199 CALL chkxer(
'ZGELSD', infot, nout, lerr, ok )
201 CALL zgelsd( 0, 0, -1, a, 1, b, 1, s, rcond, irnk, w, 10, rw,
203 CALL chkxer(
'ZGELSD', infot, nout, lerr, ok )
205 CALL zgelsd( 2, 0, 0, a, 1, b, 2, s, rcond, irnk, w, 10, rw,
207 CALL chkxer(
'ZGELSD', infot, nout, lerr, ok )
209 CALL zgelsd( 2, 0, 0, a, 2, b, 1, s, rcond, irnk, w, 10, rw,
211 CALL chkxer(
'ZGELSD', infot, nout, lerr, ok )
213 CALL zgelsd( 2, 2, 1, a, 2, b, 2, s, rcond, irnk, w, 1, rw, ip,
215 CALL chkxer(
'ZGELSD', infot, nout, lerr, ok )
220 CALL alaesm( path, ok, nout )
subroutine chkxer(SRNAMT, INFOT, NOUT, LERR, OK)
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
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
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 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