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,