151 SUBROUTINE dlsets( 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 A( lda, * ), AF( lda, * ), B( ldb, * ),
167 $ bf( ldb, * ), c( * ), cf( * ), d( * ), df( * ),
168 $ result( 2 ), rwork( * ), work( lwork ), x( * )
181 CALL
dlacpy(
'Full', m, n, a, lda, af, lda )
182 CALL
dlacpy(
'Full', p, n, b, ldb, bf, ldb )
183 CALL
dcopy( m, c, 1, cf, 1 )
184 CALL
dcopy( p, d, 1, df, 1 )
188 CALL
dgglse( m, n, p, af, lda, bf, ldb, cf, df, x, work, lwork,
195 CALL
dcopy( m, c, 1, cf, 1 )
196 CALL
dcopy( p, d, 1, df, 1 )
197 CALL
dget02(
'No transpose', m, n, 1, a, lda, x, n, cf, m, rwork,
202 CALL
dget02(
'No transpose', p, n, 1, b, ldb, x, n, df, p, rwork,
subroutine dcopy(N, DX, INCX, DY, INCY)
DCOPY
subroutine dlacpy(UPLO, M, N, A, LDA, B, LDB)
DLACPY copies all or part of one two-dimensional array to another.
subroutine dgglse(M, N, P, A, LDA, B, LDB, C, D, X, WORK, LWORK, INFO)
DGGLSE solves overdetermined or underdetermined systems for OTHER matrices
subroutine dlsets(M, P, N, A, AF, LDA, B, BF, LDB, C, CF, D, DF, X, WORK, LWORK, RWORK, RESULT)
DLSETS
subroutine dget02(TRANS, M, N, NRHS, A, LDA, X, LDX, B, LDB, RWORK, RESID)
DGET02