151 SUBROUTINE zlsets( M, P, N, A, AF, LDA, B, BF, LDB, C, CF, D, DF,
152 $ x, work, lwork, rwork, result )
160 INTEGER LDA, LDB, LWORK, M, N, P
166 DOUBLE PRECISION RESULT( 2 ), RWORK( * )
167 COMPLEX*16 A( lda, * ), AF( lda, * ), B( ldb, * ),
168 $ bf( ldb, * ), c( * ), cf( * ), d( * ), df( * ),
169 $ work( lwork ), x( * )
182 CALL
zlacpy(
'Full', m, n, a, lda, af, lda )
183 CALL
zlacpy(
'Full', p, n, b, ldb, bf, ldb )
184 CALL
zcopy( m, c, 1, cf, 1 )
185 CALL
zcopy( p, d, 1, df, 1 )
189 CALL
zgglse( m, n, p, af, lda, bf, ldb, cf, df, x, work, lwork,
196 CALL
zcopy( m, c, 1, cf, 1 )
197 CALL
zcopy( p, d, 1, df, 1 )
198 CALL
zget02(
'No transpose', m, n, 1, a, lda, x, n, cf, m, rwork,
203 CALL
zget02(
'No transpose', p, n, 1, b, ldb, x, n, df, p, rwork,
subroutine zlsets(M, P, N, A, AF, LDA, B, BF, LDB, C, CF, D, DF, X, WORK, LWORK, RWORK, RESULT)
ZLSETS
subroutine zlacpy(UPLO, M, N, A, LDA, B, LDB)
ZLACPY copies all or part of one two-dimensional array to another.
subroutine zcopy(N, ZX, INCX, ZY, INCY)
ZCOPY
subroutine zget02(TRANS, M, N, NRHS, A, LDA, X, LDX, B, LDB, RWORK, RESID)
ZGET02
subroutine zgglse(M, N, P, A, LDA, B, LDB, C, D, X, WORK, LWORK, INFO)
ZGGLSE solves overdetermined or underdetermined systems for OTHER matrices