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 cgglse(M, N, P, A, LDA, B, LDB, C, D, X, WORK, LWORK, INFO)
CGGLSE solves overdetermined or underdetermined systems for OTHER matrices
subroutine cget02(TRANS, M, N, NRHS, A, LDA, X, LDX, B, LDB, RWORK, RESID)
CGET02
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 clsets(M, P, N, A, AF, LDA, B, BF, LDB, C, CF, D, DF, X, WORK, LWORK, RWORK, RESULT)
CLSETS