Merge pull request #378 from vladimir-ch/lapack/lantr-avoid-outofbound-write

LAPACK: avoid out-of-bound write in ?LANTR

Author: langou

SRC/clantr.f

@@ -287,7 +287,7 @@ REAL             FUNCTION CLANTR( NORM, UPLO, DIAG, M, N, A, LDA,
             END IF
          ELSE
             IF( LSAME( DIAG, 'U' ) ) THEN
-               DO 210 I = 1, N
+               DO 210 I = 1, MIN( M, N )
                   WORK( I ) = ONE
   210          CONTINUE
                DO 220 I = N + 1, M

SRC/dlantr.f

@@ -285,7 +285,7 @@ DOUBLE PRECISION FUNCTION DLANTR( NORM, UPLO, DIAG, M, N, A, LDA,
             END IF
          ELSE
             IF( LSAME( DIAG, 'U' ) ) THEN
-               DO 210 I = 1, N
+               DO 210 I = 1, MIN( M, N )
                   WORK( I ) = ONE
   210          CONTINUE
                DO 220 I = N + 1, M

SRC/slantr.f

@@ -285,7 +285,7 @@ REAL             FUNCTION SLANTR( NORM, UPLO, DIAG, M, N, A, LDA,
             END IF
          ELSE
             IF( LSAME( DIAG, 'U' ) ) THEN
-               DO 210 I = 1, N
+               DO 210 I = 1, MIN( M, N )
                   WORK( I ) = ONE
   210          CONTINUE
                DO 220 I = N + 1, M

SRC/zlantr.f

@@ -287,7 +287,7 @@ DOUBLE PRECISION FUNCTION ZLANTR( NORM, UPLO, DIAG, M, N, A, LDA,
             END IF
          ELSE
             IF( LSAME( DIAG, 'U' ) ) THEN
-               DO 210 I = 1, N
+               DO 210 I = 1, MIN( M, N )
                   WORK( I ) = ONE
   210          CONTINUE
                DO 220 I = N + 1, M