diff options
author | Julien Langou <julien.langou@ucdenver.edu> | 2016-12-04 11:29:15 +0100 |
---|---|---|
committer | Julien Langou <julien.langou@ucdenver.edu> | 2016-12-04 11:29:15 +0100 |
commit | a9fcb8e2b98ff5aedbf17aa7be51ad15c617b14e (patch) | |
tree | c4513b54cd19163a82649b268d14ef25fb63146d | |
parent | f0875b93c4322544212af8d800110a7f7681e3ab (diff) |
remove MB and NB from the code, they are not used and this routine should be agnostic to them
-rw-r--r-- | SRC/dgetsls.f | 12 |
1 files changed, 4 insertions, 8 deletions
diff --git a/SRC/dgetsls.f b/SRC/dgetsls.f index c61938da..bf3aa1bf 100644 --- a/SRC/dgetsls.f +++ b/SRC/dgetsls.f @@ -159,7 +159,7 @@ * * .. Scalar Arguments .. CHARACTER TRANS - INTEGER INFO, LDA, LDB, LWORK, M, N, NRHS, MB + INTEGER INFO, LDA, LDB, LWORK, M, N, NRHS * .. * .. Array Arguments .. DOUBLE PRECISION A( LDA, * ), B( LDB, * ), WORK( * ) @@ -176,7 +176,7 @@ LOGICAL LQUERY, TRAN INTEGER I, IASCL, IBSCL, J, MINMN, MAXMN, BROW, LW, $ SCLLEN, MNK, WSIZEO, WSIZEM, LW1, LW2, - $ INFO2, NB + $ INFO2 DOUBLE PRECISION ANRM, BIGNUM, BNRM, SMLNUM * .. * .. External Functions .. @@ -225,21 +225,17 @@ IF ( M.GE.N ) THEN CALL DGEQR( M, N, A, LDA, WORK(1), -1, WORK(6), -1, $ INFO2) - MB = INT(WORK(4)) - NB = INT(WORK(5)) LW = INT(WORK(6)) CALL DGEMQR( 'L', TRANS, M, NRHS, N, A, LDA, WORK(1), - $ INT(WORK(2)), B, LDB, WORK(6), -1 , INFO2 ) + $ -1, B, LDB, WORK(6), -1 , INFO2 ) WSIZEO = INT(WORK(2))+MAX(LW,INT(WORK(6))) WSIZEM = INT(WORK(3))+MAX(LW,INT(WORK(6))) ELSE CALL DGELQ( M, N, A, LDA, WORK(1), -1, WORK(6), -1, $ INFO2) - MB = INT(WORK(4)) - NB = INT(WORK(5)) LW = INT(WORK(6)) CALL DGEMLQ( 'L', TRANS, N, NRHS, M, A, LDA, WORK(1), - $ INT(WORK(2)), B, LDB, WORK(6), -1 , INFO2 ) + $ -1, B, LDB, WORK(6), -1 , INFO2 ) WSIZEO = INT(WORK(2))+MAX(LW,INT(WORK(6))) WSIZEM = INT(WORK(3))+MAX(LW,INT(WORK(6))) END IF |