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