fixed bug.

diff --git a/SRC/ilaclr.f b/SRC/ilaclr.f
index 3d303c1..0028191 100644 (file)
--- a/SRC/ilaclr.f
+++ b/SRC/ilaclr.f
@@ -111,12 +111,9 @@
          ILACLR = 0
          DO J = 1, N
             I=M
-            DO WHILE ((A(I,J).NE.ZERO).AND.(I.GE.1))
-              I=I-1
-              IF (I.EQ.0) THEN
-                 EXIT
-              END IF
-            ENDDO         
+            DO WHILE((A(MAX(I,1),J).EQ.ZERO).AND.(I.GE.1))
+               I=I-1
+            ENDDO
             ILACLR = MAX( ILACLR, I )
          END DO
       END IF

diff --git a/SRC/iladlr.f b/SRC/iladlr.f
index ef1aef5..218c711 100644 (file)
--- a/SRC/iladlr.f
+++ b/SRC/iladlr.f
@@ -111,11 +111,8 @@
          ILADLR = 0
          DO J = 1, N
             I=M
-            DO WHILE ((A(I,J).NE.ZERO).AND.(I.GE.1))
-              I=I-1
-              IF (I.EQ.0) THEN
-                 EXIT
-              END IF
+            DO WHILE((A(MAX(I,1),J).EQ.ZERO).AND.(I.GE.1))
+               I=I-1
             ENDDO
             ILADLR = MAX( ILADLR, I )
          END DO

diff --git a/SRC/ilaslr.f b/SRC/ilaslr.f
index 285e51e..2233f01 100644 (file)
--- a/SRC/ilaslr.f
+++ b/SRC/ilaslr.f
@@ -111,11 +111,8 @@
          ILASLR = 0
          DO J = 1, N
             I=M
-            DO WHILE ((A(I,J).NE.ZERO).AND.(I.GE.1))
-              I=I-1
-              IF (I.EQ.0) THEN
-                 EXIT
-              END IF
+            DO WHILE((A(MAX(I,1),J).EQ.ZERO).AND.(I.GE.1))
+               I=I-1
             ENDDO
             ILASLR = MAX( ILASLR, I )
          END DO

diff --git a/SRC/ilazlr.f b/SRC/ilazlr.f
index 7fc1db2..dc5ceef 100644 (file)
--- a/SRC/ilazlr.f
+++ b/SRC/ilazlr.f
@@ -111,12 +111,9 @@
          ILAZLR = 0
          DO J = 1, N
             I=M
-            DO WHILE ((A(I,J).NE.ZERO).AND.(I.GE.1))
-              I=I-1
-              IF (I.EQ.0) THEN
-                 EXIT
-              END IF
-            ENDDO         
+            DO WHILE((A(MAX(I,1),J).EQ.ZERO).AND.(I.GE.1))
+               I=I-1
+            ENDDO
             ILAZLR = MAX( ILAZLR, I )
          END DO
       END IF