--- rpl/lapack/lapack/dgesvj.f 2012/08/22 09:48:14 1.10 +++ rpl/lapack/lapack/dgesvj.f 2012/12/14 12:30:20 1.11 @@ -258,7 +258,7 @@ *> \author Univ. of Colorado Denver *> \author NAG Ltd. * -*> \date November 2011 +*> \date September 2012 * *> \ingroup doubleGEcomputational * @@ -335,10 +335,10 @@ SUBROUTINE DGESVJ( JOBA, JOBU, JOBV, M, N, A, LDA, SVA, MV, V, $ LDV, WORK, LWORK, INFO ) * -* -- LAPACK computational routine (version 3.4.0) -- +* -- LAPACK computational routine (version 3.4.2) -- * -- LAPACK is a software package provided by Univ. of Tennessee, -- * -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- -* November 2011 +* September 2012 * * .. Scalar Arguments .. INTEGER INFO, LDA, LDV, LWORK, M, MV, N @@ -1573,11 +1573,11 @@ END IF * * Undo scaling, if necessary (and possible). - IF( ( ( SKL.GT.ONE ) .AND. ( SVA( 1 ).LT.( BIG / - $ SKL) ) ) .OR. ( ( SKL.LT.ONE ) .AND. ( SVA( N2 ).GT. + IF( ( ( SKL.GT.ONE ) .AND. ( SVA( 1 ).LT.( BIG / SKL) ) ) + $ .OR. ( ( SKL.LT.ONE ) .AND. ( SVA( MAX( N2, 1 ) ) .GT. $ ( SFMIN / SKL) ) ) ) THEN DO 2400 p = 1, N - SVA( p ) = SKL*SVA( p ) + SVA( P ) = SKL*SVA( P ) 2400 CONTINUE SKL= ONE END IF