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 ) )
subroutine clsets(M, P, N, A, AF, LDA, B, BF, LDB, C, CF, D, DF, X, WORK, LWORK, RWORK, RESULT)
CLSETS
subroutine cgglse(M, N, P, A, LDA, B, LDB, C, D, X, WORK, LWORK, INFO)
CGGLSE solves overdetermined or underdetermined systems for OTHER matrices
subroutine clacpy(UPLO, M, N, A, LDA, B, LDB)
CLACPY copies all or part of one two-dimensional array to another.
subroutine ccopy(N, CX, INCX, CY, INCY)
CCOPY
subroutine cget02(TRANS, M, N, NRHS, A, LDA, X, LDX, B, LDB, RWORK, RESID)
CGET02