122 SUBROUTINE sgeqls( M, N, NRHS, A, LDA, TAU, B, LDB, WORK, LWORK,
131 INTEGER info, lda, ldb, lwork, m, n, nrhs
134 REAL a( lda, * ),
b( ldb, * ), tau( * ),
142 parameter( one = 1.0e+0 )
157 ELSE IF( n.LT.0 .OR. n.GT.m )
THEN
159 ELSE IF( nrhs.LT.0 )
THEN
161 ELSE IF( lda.LT.max( 1, m ) )
THEN
163 ELSE IF( ldb.LT.max( 1, m ) )
THEN
165 ELSE IF( lwork.LT.1 .OR. lwork.LT.nrhs .AND. m.GT.0 .AND. n.GT.0 )
170 CALL
xerbla(
'SGEQLS', -info )
176 IF( n.EQ.0 .OR. nrhs.EQ.0 .OR. m.EQ.0 )
181 CALL
sormql(
'Left',
'Transpose', m, nrhs, n, a, lda, tau,
b, ldb,
182 $ work, lwork, info )
186 CALL
strsm(
'Left',
'Lower',
'No transpose',
'Non-unit', n, nrhs,
187 $ one, a( m-n+1, 1 ), lda,
b( m-n+1, 1 ), ldb )
subroutine sormql(SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC, WORK, LWORK, INFO)
SORMQL
subroutine xerbla(SRNAME, INFO)
XERBLA
subroutine sgeqls(M, N, NRHS, A, LDA, TAU, B, LDB, WORK, LWORK, INFO)
SGEQLS
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real b(3) integer i
subroutine strsm(SIDE, UPLO, TRANSA, DIAG, M, N, ALPHA, A, LDA, B, LDB)
STRSM