155 SUBROUTINE clsets( M, P, N, A, AF, LDA, B, BF, LDB, C, CF,
156 $ d, df, x, work, lwork, rwork, result )
164 INTEGER lda, ldb, lwork, m, p, n
167 REAL result( 2 ), rwork( * )
168 COMPLEX a( lda, * ), af( lda, * ), b( ldb, * ),
169 $ bf( ldb, * ), c( * ), d( * ), cf( * ),
170 $ df( * ), work( lwork ), x( * )
186 CALL
clacpy(
'Full', m, n, a, lda, af, lda )
187 CALL
clacpy(
'Full', p, n, b, ldb, bf, ldb )
188 CALL
ccopy( m, c, 1, cf, 1 )
189 CALL
ccopy( p, d, 1, df, 1 )
193 CALL
cgglse( m, n, p, af, lda, bf, ldb, cf, df, x,
194 $ work, lwork, info )
200 CALL
ccopy( m, c, 1, cf, 1 )
201 CALL
ccopy( p, d, 1, df, 1 )
202 CALL
cget02(
'No transpose', m, n, 1, a, lda, x, n, cf, m,
203 $ rwork, result( 1 ) )
207 CALL
cget02(
'No transpose', p, n, 1, b, ldb, x, n, df, p,
208 $ rwork, result( 2 ) )