updating documentation, using xLARF1y where applicable, and removing some extraneous variables

Author: jprhyne

SRC/dlarf1f.f

@@ -36,7 +36,7 @@
 *>
 *> \verbatim
 *>
-*> DLARF applies a real elementary reflector H to a real m by n matrix
+*> DLARF1F applies a real elementary reflector H to a real m by n matrix
 *> C, from either the left or the right. H is represented in the form
 *>
 *>       H = I - tau * v * v**T
@@ -193,7 +193,7 @@ SUBROUTINE DLARF1F( SIDE, M, N, V, INCV, TAU, C, LDC, WORK )
 *     .. Executable Statements ..
 *
       APPLYLEFT = LSAME( SIDE, 'L' )
-      LASTV = 0
+      LASTV = 1
       LASTC = 0
       IF( TAU.NE.ZERO ) THEN
 !     Set up variables for scanning V.  LASTV begins pointing to the end
@@ -222,17 +222,17 @@ SUBROUTINE DLARF1F( SIDE, M, N, V, INCV, TAU, C, LDC, WORK )
 !     Scan for the last non-zero row in C(:,1:lastv).
             LASTC = ILADLR(M, LASTV, C, LDC)
          END IF
-      ELSE
-!        TAU is 0, so H = I. Meaning HC = C = CH.
-         RETURN
       END IF
       IF( APPLYLEFT ) THEN
 *
 *        Form  H * C
 *
          ! Check if lastv = 1. This means v = 1, So we just need to compute
          ! C := HC = (1-\tau)C.
-         IF( LASTV.LE.1 ) THEN
+         IF( LASTV.EQ.1 ) THEN
+*
+*           C(1,1:lastc) := ( 1 - tau ) * C(1,1:lastc)
+*
             CALL DSCAL(LASTC, ONE - TAU, C, LDC)
          ELSE
 *
@@ -260,6 +260,9 @@ SUBROUTINE DLARF1F( SIDE, M, N, V, INCV, TAU, C, LDC, WORK )
          ! Check if n = 1. This means v = 1, so we just need to compute
          ! C := CH = C(1-\tau).
          IF( LASTV.EQ.1 ) THEN
+*
+*           C(1:lastc,1) := ( 1 - tau ) * C(1:lastc,1)
+*
             CALL DSCAL(LASTC, ONE - TAU, C, 1)
          ELSE
 *

SRC/dlarf1l.f

@@ -8,7 +8,7 @@
 *            http://www.netlib.org/lapack/explore-html/
 *
 *> \htmlonly
-*> Download DLARF + dependencies
+*> Download DLARF1L + dependencies
 *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/dlarf1l.f">
 *> [TGZ]</a>
 *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/dlarf1l.f">
@@ -147,7 +147,7 @@ SUBROUTINE DLARF1L( SIDE, M, N, V, INCV, TAU, C, LDC, WORK )
 *     ..
 *     .. Local Scalars ..
       LOGICAL            APPLYLEFT
-      INTEGER            I, J, FIRSTV, LASTV, LASTC
+      INTEGER            I, FIRSTV, LASTV, LASTC
 *     ..
 *     .. External Subroutines ..
       EXTERNAL           DGEMV, DGER
@@ -183,8 +183,6 @@ SUBROUTINE DLARF1L( SIDE, M, N, V, INCV, TAU, C, LDC, WORK )
 !     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.