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 dlacpy(UPLO, M, N, A, LDA, B, LDB)
DLACPY copies all or part of one two-dimensional array to another.
subroutine dcopy(N, DX, INCX, DY, INCY)
DCOPY
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 dget02(TRANS, M, N, NRHS, A, LDA, X, LDX, B, LDB, RWORK, RESID)
DGET02
subroutine dlsets(M, P, N, A, AF, LDA, B, BF, LDB, C, CF, D, DF, X, WORK, LWORK, RWORK, RESULT)
DLSETS