updating dlarf1f and dlarf1l to fix a bug found within dorg2l

Author: jprhyne

SRC/dlarf1f.f

@@ -211,7 +211,7 @@ SUBROUTINE DLARF1F( SIDE, M, N, V, INCV, TAU, C, LDC, WORK )
 !     Look for the last non-zero row in V.
 !        Since we are assuming that V(1) = 1, and it is not stored, so we
 !        shouldn't access it.
-         DO WHILE( LASTV.GT.1 .AND. V( I ).EQ.ZERO )
+         DO WHILE( LASTV.GE.2 .AND. V( I ).EQ.ZERO )
             LASTV = LASTV - 1
             I = I - INCV
          END DO
@@ -232,7 +232,7 @@ SUBROUTINE DLARF1F( SIDE, M, N, V, INCV, TAU, C, LDC, WORK )
 *
          ! Check if lastv = 1. This means v = 1, So we just need to compute
          ! C := HC = (1-\tau)C.
-         IF( LASTV.EQ.1 ) THEN
+         IF( LASTV.LE.1 ) THEN
             CALL DSCAL(LASTC, ONE - TAU, C, LDC)
          ELSE
 *

SRC/dlarf1l.f

@@ -170,23 +170,16 @@ SUBROUTINE DLARF1L( SIDE, M, N, V, INCV, TAU, C, LDC, WORK )
          ELSE
             LASTV = N
          END IF
-         IF( INCV.GT.0 ) THEN
-            I = 1 + (LASTV-1) * INCV
-         ELSE
-            I = 1
-         END IF
 !     Look for the last non-zero row in V.
-         DO WHILE( LASTV.GT.0 .AND. V( I ).EQ.ZERO )
-            LASTV = LASTV - 1
-            I = I - INCV
-         END DO
          IF( APPLYLEFT ) THEN
 !     Scan for the last non-zero column in C(1:lastv,:).
             LASTC = ILADLC(LASTV, N, C, LDC)
          ELSE
 !     Scan for the last non-zero row in C(:,1:lastv).
             LASTC = ILADLR(M, LASTV, C, LDC)
          END IF
+      ELSE
+         RETURN
       END IF
 !     Note that lastc.eq.0 renders the BLAS operations null; no special
 !     case is needed at this level.

SRC/dorg2l.f

@@ -176,6 +176,7 @@ SUBROUTINE DORG2L( M, N, K, A, LDA, TAU, WORK, INFO )
 *
 *        Apply H(i) to A(1:m-k+i,1:n-k+i) from the left
 *
+         !A(M-N+II, II) = ONE
          CALL DLARF1L( 'Left', M-N+II, II-1, A( 1, II ), 1, TAU( I ),
      $               A,
      $               LDA, WORK )