diff options
| author | jason <jason@8a072113-8704-0410-8d35-dd094bca7971> | 2008-10-28 01:38:50 +0000 |
|---|---|---|
| committer | jason <jason@8a072113-8704-0410-8d35-dd094bca7971> | 2008-10-28 01:38:50 +0000 |
| commit | baba851215b44ac3b60b9248eb02bcce7eb76247 (patch) | |
| tree | 8c0f5c006875532a30d4409f5e94b0f310ff00a7 /SRC/iladlc.f | |
| download | lapack-baba851215b44ac3b60b9248eb02bcce7eb76247.tar.gz lapack-baba851215b44ac3b60b9248eb02bcce7eb76247.tar.bz2 lapack-baba851215b44ac3b60b9248eb02bcce7eb76247.zip | |
Move LAPACK trunk into position.
Diffstat (limited to 'SRC/iladlc.f')
| -rw-r--r-- | SRC/iladlc.f | 58 |
1 files changed, 58 insertions, 0 deletions
diff --git a/SRC/iladlc.f b/SRC/iladlc.f new file mode 100644 index 00000000..fb983d6b --- /dev/null +++ b/SRC/iladlc.f @@ -0,0 +1,58 @@ + INTEGER FUNCTION ILADLC(M, N, A, LDA) + IMPLICIT NONE +! +! -- LAPACK auxiliary routine (version 3.1) -- +! Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd.. +! December 2007 +! +! .. Scalar Arguments .. + INTEGER M, N, LDA +! .. +! .. Array Arguments .. + DOUBLE PRECISION A( LDA, * ) +! .. +! +! Purpose +! ======= +! +! ILADLC scans A for its last non-zero column. +! +! Arguments +! ========= +! +! M (input) INTEGER +! The number of rows of the matrix A. +! +! N (input) INTEGER +! The number of columns of the matrix A. +! +! A (input) DOUBLE PRECISION array, dimension (LDA,N) +! The m by n matrix A. +! +! LDA (input) INTEGER +! The leading dimension of the array A. LDA >= max(1,M). +! +! ===================================================================== +! +! .. Parameters .. + DOUBLE PRECISION ZERO + PARAMETER ( ZERO = 0.0D+0 ) +! .. +! .. Local Scalars .. + INTEGER I, J +! .. +! .. Executable Statements .. +! +! Quick test for the common case where one corner is non-zero. + IF( N.EQ.0 .OR. A(1, N).NE.ZERO .OR. A(M, N).NE.ZERO ) THEN + ILADLC = N + ELSE +! Now scan each column from the end, returning with the first non-zero. + DO ILADLC = N, 1, -1 + DO I = 1, M + IF( A(I, ILADLC).NE.ZERO ) RETURN + END DO + END DO + END IF + RETURN + END FUNCTION |