155 SUBROUTINE slsets( 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 a( lda, * ), af( lda, * ),
b( ldb, * ),
168 $ bf( ldb, * ), result( 2 ), rwork( * ),
169 $ c( * ), d( * ), cf( * ), df( * ),
170 $ work( lwork ),
x( * )
186 CALL
slacpy(
'Full', m, n, a, lda, af, lda )
187 CALL
slacpy(
'Full', p, n,
b, ldb, bf, ldb )
188 CALL
scopy( m, c, 1, cf, 1 )
189 CALL
scopy( p, d, 1, df, 1 )
193 CALL
sgglse( m, n, p, af, lda, bf, ldb, cf, df,
x,
194 $ work, lwork, info )
200 CALL
scopy( m, c, 1, cf, 1 )
201 CALL
scopy( p, d, 1, df, 1 )
202 CALL
sget02(
'No transpose', m, n, 1, a, lda,
x, n, cf, m,
203 $ rwork, result( 1 ) )
207 CALL
sget02(
'No transpose', p, n, 1,
b, ldb,
x, n, df, p,
208 $ rwork, result( 2 ) )
subroutine sget02(TRANS, M, N, NRHS, A, LDA, X, LDX, B, LDB, RWORK, RESID)
SGET02
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real b(3) integer i
subroutine scopy(N, SX, INCX, SY, INCY)
SCOPY
subroutine sgglse(M, N, P, A, LDA, B, LDB, C, D, X, WORK, LWORK, INFO)
SGGLSE solves overdetermined or underdetermined systems for OTHER matrices
subroutine slacpy(UPLO, M, N, A, LDA, B, LDB)
SLACPY copies all or part of one two-dimensional array to another.
subroutine slsets(M, P, N, A, AF, LDA, B, BF, LDB, C, CF, D, DF, X, WORK, LWORK, RWORK, RESULT)
SLSETS