78 REAL RW( NMAX ), S( NMAX )
79 COMPLEX A( NMAX, NMAX ), B( NMAX, NMAX ), W( NMAX )
94 COMMON / infoc / infot, nout, ok, lerr
95 COMMON / srnamc / srnamt
101 a( 1, 1 ) = ( 1.0e+0, 0.0e+0 )
102 a( 1, 2 ) = ( 2.0e+0, 0.0e+0 )
103 a( 2, 2 ) = ( 3.0e+0, 0.0e+0 )
104 a( 2, 1 ) = ( 4.0e+0, 0.0e+0 )
106 WRITE( nout, fmt = * )
110 IF(
lsamen( 2, c2,
'LS' ) )
THEN
116 CALL cgels(
'/', 0, 0, 0, a, 1, b, 1, w, 1, info )
117 CALL chkxer(
'CGELS ', infot, nout, lerr, ok )
119 CALL cgels(
'N', -1, 0, 0, a, 1, b, 1, w, 1, info )
120 CALL chkxer(
'CGELS ', infot, nout, lerr, ok )
122 CALL cgels(
'N', 0, -1, 0, a, 1, b, 1, w, 1, info )
123 CALL chkxer(
'CGELS ', infot, nout, lerr, ok )
125 CALL cgels(
'N', 0, 0, -1, a, 1, b, 1, w, 1, info )
126 CALL chkxer(
'CGELS ', infot, nout, lerr, ok )
128 CALL cgels(
'N', 2, 0, 0, a, 1, b, 2, w, 2, info )
129 CALL chkxer(
'CGELS ', infot, nout, lerr, ok )
131 CALL cgels(
'N', 2, 0, 0, a, 2, b, 1, w, 2, info )
132 CALL chkxer(
'CGELS ', infot, nout, lerr, ok )
134 CALL cgels(
'N', 1, 1, 0, a, 1, b, 1, w, 1, info )
135 CALL chkxer(
'CGELS ', infot, nout, lerr, ok )
141 CALL cgelss( -1, 0, 0, a, 1, b, 1, s, rcond, irnk, w, 1, rw,
143 CALL chkxer(
'CGELSS', infot, nout, lerr, ok )
145 CALL cgelss( 0, -1, 0, a, 1, b, 1, s, rcond, irnk, w, 1, rw,
147 CALL chkxer(
'CGELSS', infot, nout, lerr, ok )
149 CALL cgelss( 0, 0, -1, a, 1, b, 1, s, rcond, irnk, w, 1, rw,
151 CALL chkxer(
'CGELSS', infot, nout, lerr, ok )
153 CALL cgelss( 2, 0, 0, a, 1, b, 2, s, rcond, irnk, w, 2, rw,
155 CALL chkxer(
'CGELSS', infot, nout, lerr, ok )
157 CALL cgelss( 2, 0, 0, a, 2, b, 1, s, rcond, irnk, w, 2, rw,
159 CALL chkxer(
'CGELSS', infot, nout, lerr, ok )
165 CALL cgelsy( -1, 0, 0, a, 1, b, 1, ip, rcond, irnk, w, 10, rw,
167 CALL chkxer(
'CGELSY', infot, nout, lerr, ok )
169 CALL cgelsy( 0, -1, 0, a, 1, b, 1, ip, rcond, irnk, w, 10, rw,
171 CALL chkxer(
'CGELSY', infot, nout, lerr, ok )
173 CALL cgelsy( 0, 0, -1, a, 1, b, 1, ip, rcond, irnk, w, 10, rw,
175 CALL chkxer(
'CGELSY', infot, nout, lerr, ok )
177 CALL cgelsy( 2, 0, 0, a, 1, b, 2, ip, rcond, irnk, w, 10, rw,
179 CALL chkxer(
'CGELSY', infot, nout, lerr, ok )
181 CALL cgelsy( 2, 0, 0, a, 2, b, 1, ip, rcond, irnk, w, 10, rw,
183 CALL chkxer(
'CGELSY', infot, nout, lerr, ok )
185 CALL cgelsy( 0, 3, 0, a, 1, b, 3, ip, rcond, irnk, w, 1, rw,
187 CALL chkxer(
'CGELSY', infot, nout, lerr, ok )
193 CALL cgelsd( -1, 0, 0, a, 1, b, 1, s, rcond, irnk, w, 10,
195 CALL chkxer(
'CGELSD', infot, nout, lerr, ok )
197 CALL cgelsd( 0, -1, 0, a, 1, b, 1, s, rcond, irnk, w, 10,
199 CALL chkxer(
'CGELSD', infot, nout, lerr, ok )
201 CALL cgelsd( 0, 0, -1, a, 1, b, 1, s, rcond, irnk, w, 10,
203 CALL chkxer(
'CGELSD', infot, nout, lerr, ok )
205 CALL cgelsd( 2, 0, 0, a, 1, b, 2, s, rcond, irnk, w, 10,
207 CALL chkxer(
'CGELSD', infot, nout, lerr, ok )
209 CALL cgelsd( 2, 0, 0, a, 2, b, 1, s, rcond, irnk, w, 10,
211 CALL chkxer(
'CGELSD', infot, nout, lerr, ok )
213 CALL cgelsd( 2, 2, 1, a, 2, b, 2, s, rcond, irnk, w, 1,
215 CALL chkxer(
'CGELSD', infot, nout, lerr, ok )
220 CALL alaesm( path, ok, nout )
subroutine chkxer(SRNAMT, INFOT, NOUT, LERR, OK)
logical function lsamen(N, CA, CB)
LSAMEN
subroutine alaesm(PATH, OK, NOUT)
ALAESM
subroutine cgels(TRANS, M, N, NRHS, A, LDA, B, LDB, WORK, LWORK, INFO)
CGELS solves overdetermined or underdetermined systems for GE matrices
subroutine cgelss(M, N, NRHS, A, LDA, B, LDB, S, RCOND, RANK, WORK, LWORK, RWORK, INFO)
CGELSS solves overdetermined or underdetermined systems for GE matrices
subroutine cgelsy(M, N, NRHS, A, LDA, B, LDB, JPVT, RCOND, RANK, WORK, LWORK, RWORK, INFO)
CGELSY solves overdetermined or underdetermined systems for GE matrices
subroutine cgelsd(M, N, NRHS, A, LDA, B, LDB, S, RCOND, RANK, WORK, LWORK, RWORK, IWORK, INFO)
CGELSD computes the minimum-norm solution to a linear least squares problem for GE matrices