Merge pull request Reference-LAPACK#728 from angsch/master

langou · web-flow · commit 28f7e8309608 · 2022-10-12T18:19:11.000-06:00
Fix lower triangular band matrix generation
diff --git a/TESTING/LIN/clatb4.f b/TESTING/LIN/clatb4.f
@@ -517,17 +517,18 @@ SUBROUTINE CLATB4( PATH, IMAT, M, N, TYPE, KL, KU, ANORM, MODE,
 *
 *        Set the norm and condition number.
 *
-         IF( IMAT.EQ.2 .OR. IMAT.EQ.8 ) THEN
+         MAT = ABS( IMAT )
+         IF( MAT.EQ.2 .OR. MAT.EQ.8 ) THEN
             CNDNUM = BADC1
-         ELSE IF( IMAT.EQ.3 .OR. IMAT.EQ.9 ) THEN
+         ELSE IF( MAT.EQ.3 .OR. MAT.EQ.9 ) THEN
             CNDNUM = BADC2
          ELSE
             CNDNUM = TWO
          END IF
 *
-         IF( IMAT.EQ.4 ) THEN
+         IF( MAT.EQ.4 ) THEN
             ANORM = SMALL
-         ELSE IF( IMAT.EQ.5 ) THEN
+         ELSE IF( MAT.EQ.5 ) THEN
             ANORM = LARGE
          ELSE
             ANORM = ONE
diff --git a/TESTING/LIN/dlatb4.f b/TESTING/LIN/dlatb4.f
@@ -518,17 +518,18 @@ SUBROUTINE DLATB4( PATH, IMAT, M, N, TYPE, KL, KU, ANORM, MODE,
 *
 *        Set the norm and condition number.
 *
-         IF( IMAT.EQ.2 .OR. IMAT.EQ.8 ) THEN
+         MAT = ABS( IMAT )
+         IF( MAT.EQ.2 .OR. MAT.EQ.8 ) THEN
             CNDNUM = BADC1
-         ELSE IF( IMAT.EQ.3 .OR. IMAT.EQ.9 ) THEN
+         ELSE IF( MAT.EQ.3 .OR. MAT.EQ.9 ) THEN
             CNDNUM = BADC2
          ELSE
             CNDNUM = TWO
          END IF
 *
-         IF( IMAT.EQ.4 ) THEN
+         IF( MAT.EQ.4 ) THEN
             ANORM = SMALL
-         ELSE IF( IMAT.EQ.5 ) THEN
+         ELSE IF( MAT.EQ.5 ) THEN
             ANORM = LARGE
          ELSE
             ANORM = ONE
diff --git a/TESTING/LIN/dlattb.f b/TESTING/LIN/dlattb.f
@@ -195,7 +195,7 @@ SUBROUTINE DLATTB( IMAT, UPLO, TRANS, DIAG, ISEED, N, KD, AB,
       IF( N.LE.0 )
      $   RETURN
 *
-*     Call DLATB4 to set parameters for SLATMS.
+*     Call DLATB4 to set parameters for DLATMS.
 *
       UPPER = LSAME( UPLO, 'U' )
       IF( UPPER ) THEN
diff --git a/TESTING/LIN/dlattp.f b/TESTING/LIN/dlattp.f
@@ -187,7 +187,7 @@ SUBROUTINE DLATTP( IMAT, UPLO, TRANS, DIAG, ISEED, N, A, B, WORK,
       IF( N.LE.0 )
      $   RETURN
 *
-*     Call DLATB4 to set parameters for SLATMS.
+*     Call DLATB4 to set parameters for DLATMS.
 *
       UPPER = LSAME( UPLO, 'U' )
       IF( UPPER ) THEN
diff --git a/TESTING/LIN/dlattr.f b/TESTING/LIN/dlattr.f
@@ -193,7 +193,7 @@ SUBROUTINE DLATTR( IMAT, UPLO, TRANS, DIAG, ISEED, N, A, LDA, B,
       IF( N.LE.0 )
      $   RETURN
 *
-*     Call DLATB4 to set parameters for SLATMS.
+*     Call DLATB4 to set parameters for DLATMS.
 *
       UPPER = LSAME( UPLO, 'U' )
       IF( UPPER ) THEN
diff --git a/TESTING/LIN/slatb4.f b/TESTING/LIN/slatb4.f
@@ -518,17 +518,18 @@ SUBROUTINE SLATB4( PATH, IMAT, M, N, TYPE, KL, KU, ANORM, MODE,
 *
 *        Set the norm and condition number.
 *
-         IF( IMAT.EQ.2 .OR. IMAT.EQ.8 ) THEN
+         MAT = ABS( IMAT )
+         IF( MAT.EQ.2 .OR. MAT.EQ.8 ) THEN
             CNDNUM = BADC1
-         ELSE IF( IMAT.EQ.3 .OR. IMAT.EQ.9 ) THEN
+         ELSE IF( MAT.EQ.3 .OR. MAT.EQ.9 ) THEN
             CNDNUM = BADC2
          ELSE
             CNDNUM = TWO
          END IF
 *
-         IF( IMAT.EQ.4 ) THEN
+         IF( MAT.EQ.4 ) THEN
             ANORM = SMALL
-         ELSE IF( IMAT.EQ.5 ) THEN
+         ELSE IF( MAT.EQ.5 ) THEN
             ANORM = LARGE
          ELSE
             ANORM = ONE
diff --git a/TESTING/LIN/zlatb4.f b/TESTING/LIN/zlatb4.f
@@ -517,17 +517,18 @@ SUBROUTINE ZLATB4( PATH, IMAT, M, N, TYPE, KL, KU, ANORM, MODE,
 *
 *        Set the norm and condition number.
 *
-         IF( IMAT.EQ.2 .OR. IMAT.EQ.8 ) THEN
+         MAT = ABS( IMAT )
+         IF( MAT.EQ.2 .OR. MAT.EQ.8 ) THEN
             CNDNUM = BADC1
-         ELSE IF( IMAT.EQ.3 .OR. IMAT.EQ.9 ) THEN
+         ELSE IF( MAT.EQ.3 .OR. MAT.EQ.9 ) THEN
             CNDNUM = BADC2
          ELSE
             CNDNUM = TWO
          END IF
 *
-         IF( IMAT.EQ.4 ) THEN
+         IF( MAT.EQ.4 ) THEN
             ANORM = SMALL
-         ELSE IF( IMAT.EQ.5 ) THEN
+         ELSE IF( MAT.EQ.5 ) THEN
             ANORM = LARGE
          ELSE
             ANORM = ONE
diff --git a/TESTING/LIN/zlattb.f b/TESTING/LIN/zlattb.f
@@ -204,7 +204,7 @@ SUBROUTINE ZLATTB( IMAT, UPLO, TRANS, DIAG, ISEED, N, KD, AB,
       IF( N.LE.0 )
      $   RETURN
 *
-*     Call ZLATB4 to set parameters for CLATMS.
+*     Call ZLATB4 to set parameters for ZLATMS.
 *
       UPPER = LSAME( UPLO, 'U' )
       IF( UPPER ) THEN
diff --git a/TESTING/LIN/zlattr.f b/TESTING/LIN/zlattr.f
@@ -201,7 +201,7 @@ SUBROUTINE ZLATTR( IMAT, UPLO, TRANS, DIAG, ISEED, N, A, LDA, B,
       IF( N.LE.0 )
      $   RETURN
 *
-*     Call ZLATB4 to set parameters for CLATMS.
+*     Call ZLATB4 to set parameters for ZLATMS.
 *
       UPPER = LSAME( UPLO, 'U' )
       IF( UPPER ) THEN

Original file line number	Diff line number	Diff line change
`@@ -195,7 +195,7 @@ SUBROUTINE DLATTB( IMAT, UPLO, TRANS, DIAG, ISEED, N, KD, AB,`
`195`	`195`	`IF( N.LE.0 )`
`196`	`196`	`$ RETURN`
`197`	`197`	`*`
`198`		`-* Call DLATB4 to set parameters for SLATMS.`
	`198`	`+* Call DLATB4 to set parameters for DLATMS.`
`199`	`199`	`*`
`200`	`200`	`UPPER = LSAME( UPLO, 'U' )`
`201`	`201`	`IF( UPPER ) THEN`
Original file line number	Diff line number	Diff line change
`@@ -187,7 +187,7 @@ SUBROUTINE DLATTP( IMAT, UPLO, TRANS, DIAG, ISEED, N, A, B, WORK,`
`187`	`187`	`IF( N.LE.0 )`
`188`	`188`	`$ RETURN`
`189`	`189`	`*`
`190`		`-* Call DLATB4 to set parameters for SLATMS.`
	`190`	`+* Call DLATB4 to set parameters for DLATMS.`
`191`	`191`	`*`
`192`	`192`	`UPPER = LSAME( UPLO, 'U' )`
`193`	`193`	`IF( UPPER ) THEN`
Original file line number	Diff line number	Diff line change
`@@ -193,7 +193,7 @@ SUBROUTINE DLATTR( IMAT, UPLO, TRANS, DIAG, ISEED, N, A, LDA, B,`
`193`	`193`	`IF( N.LE.0 )`
`194`	`194`	`$ RETURN`
`195`	`195`	`*`
`196`		`-* Call DLATB4 to set parameters for SLATMS.`
	`196`	`+* Call DLATB4 to set parameters for DLATMS.`
`197`	`197`	`*`
`198`	`198`	`UPPER = LSAME( UPLO, 'U' )`
`199`	`199`	`IF( UPPER ) THEN`
Original file line number	Diff line number	Diff line change
`@@ -204,7 +204,7 @@ SUBROUTINE ZLATTB( IMAT, UPLO, TRANS, DIAG, ISEED, N, KD, AB,`
`204`	`204`	`IF( N.LE.0 )`
`205`	`205`	`$ RETURN`
`206`	`206`	`*`
`207`		`-* Call ZLATB4 to set parameters for CLATMS.`
	`207`	`+* Call ZLATB4 to set parameters for ZLATMS.`
`208`	`208`	`*`
`209`	`209`	`UPPER = LSAME( UPLO, 'U' )`
`210`	`210`	`IF( UPPER ) THEN`
Original file line number	Diff line number	Diff line change
`@@ -201,7 +201,7 @@ SUBROUTINE ZLATTR( IMAT, UPLO, TRANS, DIAG, ISEED, N, A, LDA, B,`
`201`	`201`	`IF( N.LE.0 )`
`202`	`202`	`$ RETURN`
`203`	`203`	`*`
`204`		`-* Call ZLATB4 to set parameters for CLATMS.`
	`204`	`+* Call ZLATB4 to set parameters for ZLATMS.`
`205`	`205`	`*`
`206`	`206`	`UPPER = LSAME( UPLO, 'U' )`
`207`	`207`	`IF( UPPER ) THEN`