Fixed other calls exceding column 72 with -DLAPACK_64 flag and corrected comments in dorgx family

Author: Johnathan Rhyne

SRC/dlumm.f

@@ -276,15 +276,15 @@ RECURSIVE SUBROUTINE DLUMM(SIDEL, DIAGL, DIAGU, N, ALPHA,
 *
 *        A_{12} = \alpha*L_{11}*U_{12}
 *
-            CALL DTRMM('Left', 'Lower', 'No Transpose', DIAGL, K, N-K,
-     $               ALPHA, A, LDA, A(1,K+1), LDA)
+            CALL DTRMM('Left', 'Lower', 'No Transpose', DIAGL, K,
+     $               N-K, ALPHA, A, LDA, A(1,K+1), LDA)
 *
 *        Compute A_{21}
 *
 *        A_{21} = \alpha*L_{21}*U_{11}
 *
-            CALL DTRMM('Right', 'Upper', 'No Transpose', DIAGU, N-K, K,
-     $               ALPHA, A, LDA, A(K+1,1), LDA)
+            CALL DTRMM('Right', 'Upper', 'No Transpose', DIAGU, N-K,
+     $               K, ALPHA, A, LDA, A(K+1,1), LDA)
 *
 *        Compute A_{11}
 *
@@ -329,15 +329,15 @@ RECURSIVE SUBROUTINE DLUMM(SIDEL, DIAGL, DIAGU, N, ALPHA,
 *
 *        A_{12} = \alpha*U_{12}*L_{22}
 *
-            CALL DTRMM('Right', 'Lower', 'No Transpose', DIAGL, K, N-K,
-     $               ALPHA, A(K+1,K+1), LDA, A(1,K+1), LDA)
+            CALL DTRMM('Right', 'Lower', 'No Transpose', DIAGL, K,
+     $               N-K, ALPHA, A(K+1,K+1), LDA, A(1,K+1), LDA)
 *
 *        Compute A_{21}
 *
 *        A_{21} = \alpha*U_{22}*L_{21}
 *
-            CALL DTRMM('Left', 'Upper', 'No Transpose', DIAGU, N-K, K,
-     $               ALPHA, A(K+1, K+1), LDA, A(K+1,1), LDA)
+            CALL DTRMM('Left', 'Upper', 'No Transpose', DIAGU, N-K,
+     $               K, ALPHA, A(K+1, K+1), LDA, A(K+1,1), LDA)
 *
 *        Compute A_{22}
 *

SRC/dorglq.f

@@ -144,7 +144,7 @@ SUBROUTINE DORGLQ( M, N, K, A, LDA, TAU, WORK, LWORK, INFO )
      $                   NB, NBMIN, NX
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           DLARFB0C2, DLARFT, DORGL2, XERBLA
+      EXTERNAL           DLARFB0C2, DLARFT, DORGL2, DORGLK, XERBLA
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          MAX, MIN
@@ -159,7 +159,8 @@ SUBROUTINE DORGLQ( M, N, K, A, LDA, TAU, WORK, LWORK, INFO )
 *
       INFO = 0
       NB = ILAENV( 1, 'DORGLQ', ' ', M, N, K, -1 )
-      LWKOPT = MAX( 1, M ) * NB
+      ! Only need a workspace for calls to dorgl2
+      LWKOPT = MAX( 1, M )
       WORK( 1 ) = LWKOPT
       LQUERY = ( LWORK.EQ.-1 )
       IF( M.LT.0 ) THEN
@@ -268,7 +269,7 @@ SUBROUTINE DORGLQ( M, N, K, A, LDA, TAU, WORK, LWORK, INFO )
      $            LDA, A(I,I), LDA, A(I+IB,I), LDA)
 *
 *           Apply H to columns i:n of current block
-
+*
             CALL DORGLK( IB, N-I+1, A( I, I ), LDA)
          END DO
 *
@@ -283,19 +284,18 @@ SUBROUTINE DORGLQ( M, N, K, A, LDA, TAU, WORK, LWORK, INFO )
 *           Form the triangular factor of the block reflector
 *           H = H(i) H(i+1) . . . H(i+ib-1)
 *
-            CALL DLARFT( 'Forward', 'Rowwise', N-I+1, IB, A( I, I ),
-     $                  LDA, TAU( I ), WORK, LDWORK )
+            CALL DLARFT( 'Forward', 'Transpose', N-I+1, IB, A(I,I),
+     $                  LDA, TAU( I ), A( I, I ), LDA )
+*
+*           Apply H to A(i+ib:m,i:n) from the right
 *
-*           Apply H**T to A(i+ib:m,i:n) from the right
+            CALL DLARFB0C2(.FALSE., 'Right', 'No Transpose',
+     $            'Forward', 'Rowwise', M-I-IB+1, N-I+1, IB, A(I,I),
+     $            LDA, A(I,I), LDA, A(I+IB,I), LDA)
 *
-            CALL DLARFB0C2(.FALSE., 'Right', 'Transpose', 'Forward',
-     $            'Rowwise', M-I-IB+1, N-I+1, IB, A(I,I), LDA, WORK, 
-     $            LDWORK, A(I+IB,I), LDA)
+*           Apply H to columns i:n of current block
 *
-*           Apply H**T to columns i:n of current block
-
-            CALL DORGL2( IB, N-I+1, IB, A( I, I ), LDA, TAU( I ),
-     $                   WORK, IINFO )
+            CALL DORGLK( IB, N-I+1, A( I, I ), LDA)
          END IF
       END IF
 *

SRC/dorgql.f

@@ -142,7 +142,7 @@ SUBROUTINE DORGQL( M, N, K, A, LDA, TAU, WORK, LWORK, INFO )
       INTEGER            I, IB, IINFO, IWS, KK, LWKOPT, NB, NBMIN
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           DLARFB0C2, DLARFT, DORG2L, XERBLA
+      EXTERNAL           DLARFB0C2, DLARFT, DORG2L, DORGKL, XERBLA
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          MAX, MIN

SRC/dorgrq.f

@@ -145,7 +145,7 @@ SUBROUTINE DORGRQ( M, N, K, A, LDA, TAU, WORK, LWORK, INFO )
      $                   LWKOPT, NB, NBMIN, NX
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           DLARFB0C2, DLARFT, DORGR2, XERBLA
+      EXTERNAL           DLARFB0C2, DLARFT, DORGR2, DORGRK, XERBLA
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          MAX, MIN
@@ -175,7 +175,8 @@ SUBROUTINE DORGRQ( M, N, K, A, LDA, TAU, WORK, LWORK, INFO )
             LWKOPT = 1
          ELSE
             NB = ILAENV( 1, 'DORGRQ', ' ', M, N, K, -1 )
-            LWKOPT = M*NB
+            ! Only need a workspace for calls to dorgr2
+            LWKOPT = M
          END IF
          WORK( 1 ) = LWKOPT
 *
@@ -277,13 +278,13 @@ SUBROUTINE DORGRQ( M, N, K, A, LDA, TAU, WORK, LWORK, INFO )
             CALL DLARFT( 'Transpose', 'Rowwise', N-K+I+IB-1, IB,
      $                A( II, 1 ), LDA, TAU( I ), A( II, N-K+I ), LDA )
 *
-*           Apply H**T to A(1:m-k+i-1,1:n-k+i+ib-1) from the right
+*           Apply H to A(1:m-k+i-1,1:n-k+i+ib-1) from the right
 *
             CALL DLARFB0C2(.FALSE., 'Right', 'No Transpose',
      $            'Backward', 'Rowwise', II-1, N-K+I+IB-1, IB, A(II,1),
      $             LDA, A( II, N-K+I ), LDA, A, LDA)
 *
-*           Apply H**T to columns 1:n-k+i+ib-1 of current block
+*           Apply H to columns 1:n-k+i+ib-1 of current block
 *
             CALL DORGRK( IB, N-K+I+IB-1, A( II, 1 ), LDA )
          END DO