Merge pull request #885 from mjacobse/fix_procedure_args

langou · web-flow · commit 129637ea6fb3 · 2023-07-18T08:33:40.000-06:00
Fix actual arguments in some procedure calls
diff --git a/SRC/cgelqt3.f b/SRC/cgelqt3.f
@@ -159,7 +159,8 @@ RECURSIVE SUBROUTINE CGELQT3( M, N, A, LDA, T, LDT, INFO )
 *
 *        Compute Householder transform when M=1
 *
-         CALL CLARFG( N, A, A( 1, MIN( 2, N ) ), LDA, T )
+         CALL CLARFG( N, A( 1, 1 ), A( 1, MIN( 2, N ) ), LDA,
+     &                T( 1, 1 ) )
          T(1,1)=CONJG(T(1,1))
 *
       ELSE
diff --git a/SRC/dgelqt3.f b/SRC/dgelqt3.f
@@ -173,7 +173,8 @@ RECURSIVE SUBROUTINE DGELQT3( M, N, A, LDA, T, LDT, INFO )
 *
 *        Compute Householder transform when M=1
 *
-         CALL DLARFG( N, A, A( 1, MIN( 2, N ) ), LDA, T )
+         CALL DLARFG( N, A( 1, 1 ), A( 1, MIN( 2, N ) ), LDA,
+     &                T( 1, 1 ) )
 *
       ELSE
 *
diff --git a/SRC/dlatrs.f b/SRC/dlatrs.f
@@ -261,6 +261,9 @@ SUBROUTINE DLATRS( UPLO, TRANS, DIAG, NORMIN, N, A, LDA, X, SCALE,
       DOUBLE PRECISION   BIGNUM, GROW, REC, SMLNUM, SUMJ, TJJ, TJJS,
      $                   TMAX, TSCAL, USCAL, XBND, XJ, XMAX
 *     ..
+*     .. Local Arrays ..
+      DOUBLE PRECISION   WORK(1)
+*     ..
 *     .. External Functions ..
       LOGICAL            LSAME
       INTEGER            IDAMAX
@@ -362,7 +365,7 @@ SUBROUTINE DLATRS( UPLO, TRANS, DIAG, NORMIN, N, A, LDA, X, SCALE,
 *              A is upper triangular.
 *
                DO J = 2, N
-                  TMAX = MAX( DLANGE( 'M', J-1, 1, A( 1, J ), 1, SUMJ ),
+                  TMAX = MAX( DLANGE( 'M', J-1, 1, A( 1, J ), 1, WORK ),
      $                        TMAX )
                END DO
             ELSE
@@ -371,7 +374,7 @@ SUBROUTINE DLATRS( UPLO, TRANS, DIAG, NORMIN, N, A, LDA, X, SCALE,
 *
                DO J = 1, N - 1
                   TMAX = MAX( DLANGE( 'M', N-J, 1, A( J+1, J ), 1,
-     $                        SUMJ ), TMAX )
+     $                        WORK ), TMAX )
                END DO
             END IF
 *
diff --git a/SRC/dtrsyl3.f b/SRC/dtrsyl3.f
@@ -1220,7 +1220,7 @@ SUBROUTINE DTRSYL3( TRANA, TRANB, ISGN, M, N, A, LDA, B, LDB, C,
 *
          SCALOC = MIN( BIGNUM / SCAL, ONE / BUF )
          BUF = BUF * SCALOC
-         CALL DLASCL( 'G', -1, -1, ONE, SCALOC, M, N, C, LDC, IWORK )
+         CALL DLASCL( 'G', -1, -1, ONE, SCALOC, M, N, C, LDC, IWORK(1) )
       END IF
 *
 *     Combine with buffer scaling factor. SCALE will be flushed if
diff --git a/SRC/sgelqt3.f b/SRC/sgelqt3.f
@@ -158,7 +158,8 @@ RECURSIVE SUBROUTINE SGELQT3( M, N, A, LDA, T, LDT, INFO )
 *
 *        Compute Householder transform when M=1
 *
-         CALL SLARFG( N, A, A( 1, MIN( 2, N ) ), LDA, T )
+         CALL SLARFG( N, A( 1, 1 ), A( 1, MIN( 2, N ) ), LDA,
+     &                T( 1, 1 ) )
 *
       ELSE
 *
diff --git a/SRC/slatrs.f b/SRC/slatrs.f
@@ -261,6 +261,9 @@ SUBROUTINE SLATRS( UPLO, TRANS, DIAG, NORMIN, N, A, LDA, X, SCALE,
       REAL               BIGNUM, GROW, REC, SMLNUM, SUMJ, TJJ, TJJS,
      $                   TMAX, TSCAL, USCAL, XBND, XJ, XMAX
 *     ..
+*     .. Local Arrays ..
+      REAL               WORK(1)
+*     ..
 *     .. External Functions ..
       LOGICAL            LSAME
       INTEGER            ISAMAX
@@ -362,7 +365,7 @@ SUBROUTINE SLATRS( UPLO, TRANS, DIAG, NORMIN, N, A, LDA, X, SCALE,
 *              A is upper triangular.
 *
                DO J = 2, N
-                  TMAX = MAX( SLANGE( 'M', J-1, 1, A( 1, J ), 1, SUMJ ),
+                  TMAX = MAX( SLANGE( 'M', J-1, 1, A( 1, J ), 1, WORK ),
      $                        TMAX )
                END DO
             ELSE
@@ -371,7 +374,7 @@ SUBROUTINE SLATRS( UPLO, TRANS, DIAG, NORMIN, N, A, LDA, X, SCALE,
 *
                DO J = 1, N - 1
                   TMAX = MAX( SLANGE( 'M', N-J, 1, A( J+1, J ), 1,
-     $                        SUMJ ), TMAX )
+     $                        WORK ), TMAX )
                END DO
             END IF
 *
diff --git a/SRC/strsyl3.f b/SRC/strsyl3.f
@@ -1223,7 +1223,7 @@ SUBROUTINE STRSYL3( TRANA, TRANB, ISGN, M, N, A, LDA, B, LDB, C,
 *
          SCALOC = MIN( BIGNUM / SCAL, ONE / BUF )
          BUF = BUF * SCALOC
-         CALL SLASCL( 'G', -1, -1, ONE, SCALOC, M, N, C, LDC, IWORK )
+         CALL SLASCL( 'G', -1, -1, ONE, SCALOC, M, N, C, LDC, IWORK(1) )
       END IF
 *
 *     Combine with buffer scaling factor. SCALE will be flushed if
diff --git a/SRC/zgelqt3.f b/SRC/zgelqt3.f
@@ -174,7 +174,8 @@ RECURSIVE SUBROUTINE ZGELQT3( M, N, A, LDA, T, LDT, INFO )
 *
 *        Compute Householder transform when M=1
 *
-         CALL ZLARFG( N, A, A( 1, MIN( 2, N ) ), LDA, T )
+         CALL ZLARFG( N, A( 1, 1 ), A( 1, MIN( 2, N ) ), LDA,
+     &                T( 1, 1 ) )
          T(1,1)=CONJG(T(1,1))
 *
       ELSE
diff --git a/TESTING/EIG/cerrst.f b/TESTING/EIG/cerrst.f
@@ -160,13 +160,13 @@ SUBROUTINE CERRST( PATH, NUNIT )
 *
          SRNAMT = 'CHETD2'
          INFOT = 1
-         CALL CHETD2( '/', 0, A, 1, D, E, TAU, W, 1, INFO )
+         CALL CHETD2( '/', 0, A, 1, D, E, TAU, INFO )
          CALL CHKXER( 'CHETD2', INFOT, NOUT, LERR, OK )
          INFOT = 2
-         CALL CHETD2( 'U', -1, A, 1, D, E, TAU, W, 1, INFO )
+         CALL CHETD2( 'U', -1, A, 1, D, E, TAU, INFO )
          CALL CHKXER( 'CHETD2', INFOT, NOUT, LERR, OK )
          INFOT = 4
-         CALL CHETD2( 'U', 2, A, 1, D, E, TAU, W, 1, INFO )
+         CALL CHETD2( 'U', 2, A, 1, D, E, TAU, INFO )
          CALL CHKXER( 'CHETD2', INFOT, NOUT, LERR, OK )
          NT = NT + 3
 *
diff --git a/TESTING/EIG/derrst.f b/TESTING/EIG/derrst.f
@@ -161,13 +161,13 @@ SUBROUTINE DERRST( PATH, NUNIT )
 *
          SRNAMT = 'DSYTD2'
          INFOT = 1
-         CALL DSYTD2( '/', 0, A, 1, D, E, TAU, W, 1, INFO )
+         CALL DSYTD2( '/', 0, A, 1, D, E, TAU, INFO )
          CALL CHKXER( 'DSYTD2', INFOT, NOUT, LERR, OK )
          INFOT = 2
-         CALL DSYTD2( 'U', -1, A, 1, D, E, TAU, W, 1, INFO )
+         CALL DSYTD2( 'U', -1, A, 1, D, E, TAU, INFO )
          CALL CHKXER( 'DSYTD2', INFOT, NOUT, LERR, OK )
          INFOT = 4
-         CALL DSYTD2( 'U', 2, A, 1, D, E, TAU, W, 1, INFO )
+         CALL DSYTD2( 'U', 2, A, 1, D, E, TAU, INFO )
          CALL CHKXER( 'DSYTD2', INFOT, NOUT, LERR, OK )
          NT = NT + 3
 *
diff --git a/TESTING/EIG/serrst.f b/TESTING/EIG/serrst.f
@@ -161,13 +161,13 @@ SUBROUTINE SERRST( PATH, NUNIT )
 *
          SRNAMT = 'SSYTD2'
          INFOT = 1
-         CALL SSYTD2( '/', 0, A, 1, D, E, TAU, W, 1, INFO )
+         CALL SSYTD2( '/', 0, A, 1, D, E, TAU, INFO )
          CALL CHKXER( 'SSYTD2', INFOT, NOUT, LERR, OK )
          INFOT = 2
-         CALL SSYTD2( 'U', -1, A, 1, D, E, TAU, W, 1, INFO )
+         CALL SSYTD2( 'U', -1, A, 1, D, E, TAU, INFO )
          CALL CHKXER( 'SSYTD2', INFOT, NOUT, LERR, OK )
          INFOT = 4
-         CALL SSYTD2( 'U', 2, A, 1, D, E, TAU, W, 1, INFO )
+         CALL SSYTD2( 'U', 2, A, 1, D, E, TAU, INFO )
          CALL CHKXER( 'SSYTD2', INFOT, NOUT, LERR, OK )
          NT = NT + 3
 *
diff --git a/TESTING/EIG/zerrst.f b/TESTING/EIG/zerrst.f
@@ -160,13 +160,13 @@ SUBROUTINE ZERRST( PATH, NUNIT )
 *
          SRNAMT = 'ZHETD2'
          INFOT = 1
-         CALL ZHETD2( '/', 0, A, 1, D, E, TAU, W, 1, INFO )
+         CALL ZHETD2( '/', 0, A, 1, D, E, TAU, INFO )
          CALL CHKXER( 'ZHETD2', INFOT, NOUT, LERR, OK )
          INFOT = 2
-         CALL ZHETD2( 'U', -1, A, 1, D, E, TAU, W, 1, INFO )
+         CALL ZHETD2( 'U', -1, A, 1, D, E, TAU, INFO )
          CALL CHKXER( 'ZHETD2', INFOT, NOUT, LERR, OK )
          INFOT = 4
-         CALL ZHETD2( 'U', 2, A, 1, D, E, TAU, W, 1, INFO )
+         CALL ZHETD2( 'U', 2, A, 1, D, E, TAU, INFO )
          CALL CHKXER( 'ZHETD2', INFOT, NOUT, LERR, OK )
          NT = NT + 3
 *
diff --git a/TESTING/LIN/cchktr.f b/TESTING/LIN/cchktr.f
@@ -201,7 +201,8 @@ SUBROUTINE CCHKTR( DOTYPE, NN, NVAL, NNB, NBVAL, NNS, NSVAL,
 *     .. Local Arrays ..
       CHARACTER          TRANSS( NTRAN ), UPLOS( 2 )
       INTEGER            ISEED( 4 ), ISEEDY( 4 )
-      REAL               RESULT( NTESTS ), SCALE3( 2 )
+      REAL               RESULT( NTESTS ), RWORK2( 2*NMAX ),
+     $                   SCALE3( 2 )
 *     ..
 *     .. External Functions ..
       LOGICAL            LSAME
@@ -542,10 +543,10 @@ SUBROUTINE CCHKTR( DOTYPE, NN, NVAL, NNB, NBVAL, NNS, NSVAL,
                   SRNAMT = 'CLATRS3'
                   CALL CCOPY( N, X, 1, B, 1 )
                   CALL CCOPY( N, X, 1, B( N+1 ), 1 )
-                  CALL CSCAL( N, BIGNUM, B( N+1 ), 1 )
+                  CALL CSSCAL( N, BIGNUM, B( N+1 ), 1 )
                   CALL CLATRS3( UPLO, TRANS, DIAG, 'N', N, 2, A, LDA,
-     $                          B, MAX(1, N), SCALE3, RWORK, WORK, NMAX,
-     $                          INFO )
+     $                          B, MAX(1, N), SCALE3, RWORK, RWORK2,
+     $                          2*NMAX, INFO )
 *
 *                 Check error code from CLATRS3.
 *
diff --git a/TESTING/LIN/cerrtr.f b/TESTING/LIN/cerrtr.f
@@ -70,7 +70,7 @@ SUBROUTINE CERRTR( PATH, NUNIT )
 *     .. Local Scalars ..
       CHARACTER*2        C2
       INTEGER            INFO
-      REAL               RCOND, SCALE
+      REAL               RCOND, SCALE, SCALES(0)
 *     ..
 *     .. Local Arrays ..
       REAL               R1( NMAX ), R2( NMAX ), RW( NMAX )
@@ -245,40 +245,40 @@ SUBROUTINE CERRTR( PATH, NUNIT )
 *
          SRNAMT = 'CLATRS3'
          INFOT = 1
-         CALL CLATRS3( '/', 'N', 'N', 'N', 0, 0, A, 1, X, 1, SCALE, RW,
-     $                 RW( 2 ), 1, INFO )
+         CALL CLATRS3( '/', 'N', 'N', 'N', 0, 0, A, 1, X, 1, SCALES,
+     $                 RW, RW( 2 ), 1, INFO )
          CALL CHKXER( 'CLATRS3', INFOT, NOUT, LERR, OK )
          INFOT = 2
-         CALL CLATRS3( 'U', '/', 'N', 'N', 0, 0, A, 1, X, 1, SCALE, RW,
-     $                 RW( 2 ), 1, INFO )
+         CALL CLATRS3( 'U', '/', 'N', 'N', 0, 0, A, 1, X, 1, SCALES,
+     $                 RW, RW( 2 ), 1, INFO )
          CALL CHKXER( 'CLATRS3', INFOT, NOUT, LERR, OK )
          INFOT = 3
-         CALL CLATRS3( 'U', 'N', '/', 'N', 0, 0, A, 1, X, 1, SCALE, RW,
-     $                 RW( 2 ), 1, INFO )
+         CALL CLATRS3( 'U', 'N', '/', 'N', 0, 0, A, 1, X, 1, SCALES,
+     $                 RW, RW( 2 ), 1, INFO )
          CALL CHKXER( 'CLATRS3', INFOT, NOUT, LERR, OK )
          INFOT = 4
-         CALL CLATRS3( 'U', 'N', 'N', '/', 0, 0, A, 1, X, 1, SCALE, RW,
-     $                 RW( 2 ), 1, INFO )
+         CALL CLATRS3( 'U', 'N', 'N', '/', 0, 0, A, 1, X, 1, SCALES,
+     $                 RW, RW( 2 ), 1, INFO )
          CALL CHKXER( 'CLATRS3', INFOT, NOUT, LERR, OK )
          INFOT = 5
-         CALL CLATRS3( 'U', 'N', 'N', 'N', -1, 0, A, 1, X, 1, SCALE, RW,
-     $                 RW( 2 ), 1, INFO )
+         CALL CLATRS3( 'U', 'N', 'N', 'N', -1, 0, A, 1, X, 1, SCALES,
+     $                 RW, RW( 2 ), 1, INFO )
          CALL CHKXER( 'CLATRS3', INFOT, NOUT, LERR, OK )
          INFOT = 6
-         CALL CLATRS3( 'U', 'N', 'N', 'N', 0, -1, A, 1, X, 1, SCALE, RW,
-     $                 RW( 2 ), 1, INFO )
+         CALL CLATRS3( 'U', 'N', 'N', 'N', 0, -1, A, 1, X, 1, SCALES,
+     $                 RW, RW( 2 ), 1, INFO )
          CALL CHKXER( 'CLATRS3', INFOT, NOUT, LERR, OK )
          INFOT = 8
-         CALL CLATRS3( 'U', 'N', 'N', 'N', 2, 0, A, 1, X, 1, SCALE, RW,
-     $                 RW( 2 ), 1, INFO )
+         CALL CLATRS3( 'U', 'N', 'N', 'N', 2, 0, A, 1, X, 1, SCALES,
+     $                 RW, RW( 2 ), 1, INFO )
          CALL CHKXER( 'CLATRS3', INFOT, NOUT, LERR, OK )
          INFOT = 10
-         CALL CLATRS3( 'U', 'N', 'N', 'N', 2, 0, A, 2, X, 1, SCALE, RW,
-     $                 RW( 2 ), 1, INFO )
+         CALL CLATRS3( 'U', 'N', 'N', 'N', 2, 0, A, 2, X, 1, SCALES,
+     $                 RW, RW( 2 ), 1, INFO )
          CALL CHKXER( 'CLATRS3', INFOT, NOUT, LERR, OK )
          INFOT = 14
-         CALL CLATRS3( 'U', 'N', 'N', 'N', 1, 0, A, 1, X, 1, SCALE, RW,
-     $                 RW( 2 ), 0, INFO )
+         CALL CLATRS3( 'U', 'N', 'N', 'N', 1, 0, A, 1, X, 1, SCALES,
+     $                 RW, RW( 2 ), 0, INFO )
          CALL CHKXER( 'CLATRS3', INFOT, NOUT, LERR, OK )
 *
 *     Test error exits for the packed triangular routines.
diff --git a/TESTING/LIN/derrtr.f b/TESTING/LIN/derrtr.f
@@ -71,7 +71,7 @@ SUBROUTINE DERRTR( PATH, NUNIT )
 *     .. Local Scalars ..
       CHARACTER*2        C2
       INTEGER            INFO
-      DOUBLE PRECISION   RCOND, SCALE
+      DOUBLE PRECISION   RCOND, SCALE, SCALES(0)
 *     ..
 *     .. Local Arrays ..
       INTEGER            IW( NMAX )
@@ -250,40 +250,40 @@ SUBROUTINE DERRTR( PATH, NUNIT )
 *
          SRNAMT = 'DLATRS3'
          INFOT = 1
-         CALL DLATRS3( '/', 'N', 'N', 'N', 0, 0, A, 1, X, 1, SCALE, W,
-     $                 W( 2 ), 1, INFO )
+         CALL DLATRS3( '/', 'N', 'N', 'N', 0, 0, A, 1, X, 1, SCALES,
+     $                 W, W( 2 ), 1, INFO )
          CALL CHKXER( 'DLATRS3', INFOT, NOUT, LERR, OK )
          INFOT = 2
-         CALL DLATRS3( 'U', '/', 'N', 'N', 0, 0, A, 1, X, 1, SCALE, W,
-     $                 W( 2 ), 1, INFO )
+         CALL DLATRS3( 'U', '/', 'N', 'N', 0, 0, A, 1, X, 1, SCALES,
+     $                 W, W( 2 ), 1, INFO )
          CALL CHKXER( 'DLATRS3', INFOT, NOUT, LERR, OK )
          INFOT = 3
-         CALL DLATRS3( 'U', 'N', '/', 'N', 0, 0, A, 1, X, 1, SCALE, W,
-     $                 W( 2 ), 1, INFO )
+         CALL DLATRS3( 'U', 'N', '/', 'N', 0, 0, A, 1, X, 1, SCALES,
+     $                 W, W( 2 ), 1, INFO )
          CALL CHKXER( 'DLATRS3', INFOT, NOUT, LERR, OK )
          INFOT = 4
-         CALL DLATRS3( 'U', 'N', 'N', '/', 0, 0, A, 1, X, 1, SCALE, W,
-     $                 W( 2 ), 1, INFO )
+         CALL DLATRS3( 'U', 'N', 'N', '/', 0, 0, A, 1, X, 1, SCALES,
+     $                 W, W( 2 ), 1, INFO )
          CALL CHKXER( 'DLATRS3', INFOT, NOUT, LERR, OK )
          INFOT = 5
-         CALL DLATRS3( 'U', 'N', 'N', 'N', -1, 0, A, 1, X, 1, SCALE, W,
-     $                 W( 2 ), 1, INFO )
+         CALL DLATRS3( 'U', 'N', 'N', 'N', -1, 0, A, 1, X, 1, SCALES,
+     $                 W, W( 2 ), 1, INFO )
          CALL CHKXER( 'DLATRS3', INFOT, NOUT, LERR, OK )
          INFOT = 6
-         CALL DLATRS3( 'U', 'N', 'N', 'N', 0, -1, A, 1, X, 1, SCALE, W,
-     $                 W( 2 ), 1, INFO )
+         CALL DLATRS3( 'U', 'N', 'N', 'N', 0, -1, A, 1, X, 1, SCALES,
+     $                 W, W( 2 ), 1, INFO )
          CALL CHKXER( 'DLATRS3', INFOT, NOUT, LERR, OK )
          INFOT = 8
-         CALL DLATRS3( 'U', 'N', 'N', 'N', 2, 0, A, 1, X, 1, SCALE, W,
-     $                 W( 2 ), 1, INFO )
+         CALL DLATRS3( 'U', 'N', 'N', 'N', 2, 0, A, 1, X, 1, SCALES,
+     $                 W, W( 2 ), 1, INFO )
          CALL CHKXER( 'DLATRS3', INFOT, NOUT, LERR, OK )
          INFOT = 10
-         CALL DLATRS3( 'U', 'N', 'N', 'N', 2, 0, A, 2, X, 1, SCALE, W,
-     $                 W( 2 ), 1, INFO )
+         CALL DLATRS3( 'U', 'N', 'N', 'N', 2, 0, A, 2, X, 1, SCALES,
+     $                 W, W( 2 ), 1, INFO )
          CALL CHKXER( 'DLATRS3', INFOT, NOUT, LERR, OK )
          INFOT = 14
-         CALL DLATRS3( 'U', 'N', 'N', 'N', 1, 0, A, 1, X, 1, SCALE, W,
-     $                 W( 2 ), 0, INFO )
+         CALL DLATRS3( 'U', 'N', 'N', 'N', 1, 0, A, 1, X, 1, SCALES,
+     $                 W, W( 2 ), 0, INFO )
          CALL CHKXER( 'DLATRS3', INFOT, NOUT, LERR, OK )
 *
       ELSE IF( LSAMEN( 2, C2, 'TP' ) ) THEN
diff --git a/TESTING/LIN/serrtr.f b/TESTING/LIN/serrtr.f
diff --git a/TESTING/LIN/zchktr.f b/TESTING/LIN/zchktr.f
diff --git a/TESTING/LIN/zerrtr.f b/TESTING/LIN/zerrtr.f

Original file line number	Diff line number	Diff line change
`@@ -159,7 +159,8 @@ RECURSIVE SUBROUTINE CGELQT3( M, N, A, LDA, T, LDT, INFO )`
`159`	`159`	`*`
`160`	`160`	`* Compute Householder transform when M=1`
`161`	`161`	`*`
`162`		`- CALL CLARFG( N, A, A( 1, MIN( 2, N ) ), LDA, T )`
	`162`	`+ CALL CLARFG( N, A( 1, 1 ), A( 1, MIN( 2, N ) ), LDA,`
	`163`	`+ & T( 1, 1 ) )`
`163`	`164`	`T(1,1)=CONJG(T(1,1))`
`164`	`165`	`*`
`165`	`166`	`ELSE`
Original file line number	Diff line number	Diff line change
`@@ -173,7 +173,8 @@ RECURSIVE SUBROUTINE DGELQT3( M, N, A, LDA, T, LDT, INFO )`
`173`	`173`	`*`
`174`	`174`	`* Compute Householder transform when M=1`
`175`	`175`	`*`
`176`		`- CALL DLARFG( N, A, A( 1, MIN( 2, N ) ), LDA, T )`
	`176`	`+ CALL DLARFG( N, A( 1, 1 ), A( 1, MIN( 2, N ) ), LDA,`
	`177`	`+ & T( 1, 1 ) )`
`177`	`178`	`*`
`178`	`179`	`ELSE`
`179`	`180`	`*`
Original file line number	Diff line number	Diff line change
`@@ -1220,7 +1220,7 @@ SUBROUTINE DTRSYL3( TRANA, TRANB, ISGN, M, N, A, LDA, B, LDB, C,`
`1220`	`1220`	`*`
`1221`	`1221`	`SCALOC = MIN( BIGNUM / SCAL, ONE / BUF )`
`1222`	`1222`	`BUF = BUF * SCALOC`
`1223`		`- CALL DLASCL( 'G', -1, -1, ONE, SCALOC, M, N, C, LDC, IWORK )`
	`1223`	`+ CALL DLASCL( 'G', -1, -1, ONE, SCALOC, M, N, C, LDC, IWORK(1) )`
`1224`	`1224`	`END IF`
`1225`	`1225`	`*`
`1226`	`1226`	`* Combine with buffer scaling factor. SCALE will be flushed if`
Original file line number	Diff line number	Diff line change
`@@ -158,7 +158,8 @@ RECURSIVE SUBROUTINE SGELQT3( M, N, A, LDA, T, LDT, INFO )`
`158`	`158`	`*`
`159`	`159`	`* Compute Householder transform when M=1`
`160`	`160`	`*`
`161`		`- CALL SLARFG( N, A, A( 1, MIN( 2, N ) ), LDA, T )`
	`161`	`+ CALL SLARFG( N, A( 1, 1 ), A( 1, MIN( 2, N ) ), LDA,`
	`162`	`+ & T( 1, 1 ) )`
`162`	`163`	`*`
`163`	`164`	`ELSE`
`164`	`165`	`*`
Original file line number	Diff line number	Diff line change
`@@ -1223,7 +1223,7 @@ SUBROUTINE STRSYL3( TRANA, TRANB, ISGN, M, N, A, LDA, B, LDB, C,`
`1223`	`1223`	`*`
`1224`	`1224`	`SCALOC = MIN( BIGNUM / SCAL, ONE / BUF )`
`1225`	`1225`	`BUF = BUF * SCALOC`
`1226`		`- CALL SLASCL( 'G', -1, -1, ONE, SCALOC, M, N, C, LDC, IWORK )`
	`1226`	`+ CALL SLASCL( 'G', -1, -1, ONE, SCALOC, M, N, C, LDC, IWORK(1) )`
`1227`	`1227`	`END IF`
`1228`	`1228`	`*`
`1229`	`1229`	`* Combine with buffer scaling factor. SCALE will be flushed if`
Original file line number	Diff line number	Diff line change
`@@ -174,7 +174,8 @@ RECURSIVE SUBROUTINE ZGELQT3( M, N, A, LDA, T, LDT, INFO )`
`174`	`174`	`*`
`175`	`175`	`* Compute Householder transform when M=1`
`176`	`176`	`*`
`177`		`- CALL ZLARFG( N, A, A( 1, MIN( 2, N ) ), LDA, T )`
	`177`	`+ CALL ZLARFG( N, A( 1, 1 ), A( 1, MIN( 2, N ) ), LDA,`
	`178`	`+ & T( 1, 1 ) )`
`178`	`179`	`T(1,1)=CONJG(T(1,1))`
`179`	`180`	`*`
`180`	`181`	`ELSE`