--- rpl/lapack/lapack/iladlr.f 2010/12/21 13:53:41 1.8 +++ rpl/lapack/lapack/iladlr.f 2011/07/22 07:38:13 1.9 @@ -1,12 +1,10 @@ INTEGER FUNCTION ILADLR( M, N, A, LDA ) IMPLICIT NONE * -* -- LAPACK auxiliary routine (version 3.2.2) -- -* -* -- June 2010 -- -* +* -- LAPACK auxiliary routine (version 3.3.1) -- * -- LAPACK is a software package provided by Univ. of Tennessee, -- * -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- +* -- April 2011 -- * * .. Scalar Arguments .. INTEGER M, N, LDA @@ -55,9 +53,10 @@ * Scan up each column tracking the last zero row seen. ILADLR = 0 DO J = 1, N - DO I = M, 1, -1 - IF( A(I, J).NE.ZERO ) EXIT - END DO + I=M + DO WHILE ((A(I,J).NE.ZERO).AND.(I.GE.1)) + I=I-1 + ENDDO ILADLR = MAX( ILADLR, I ) END DO END IF