make SY/HE imat match LAWN 41 (different than PO)
authorMark Gates <mgates3@utk.edu>
Mon, 17 Oct 2016 18:10:48 +0000 (14:10 -0400)
committerMark Gates <mgates3@utk.edu>
Mon, 17 Oct 2016 18:10:48 +0000 (14:10 -0400)
TESTING/LIN/clatb4.f
TESTING/LIN/dlatb4.f
TESTING/LIN/slatb4.f
TESTING/LIN/zlatb4.f

index 274179b..c936f90 100644 (file)
             ANORM = ONE
          END IF
 *
-      ELSE IF( LSAMEN( 2, C2, 'PO' ) .OR. LSAMEN( 2, C2, 'PP' ) .OR.
-     $         LSAMEN( 2, C2, 'HE' ) .OR. LSAMEN( 2, C2, 'HP' ) .OR.
-     $         LSAMEN( 2, C2, 'SY' ) .OR. LSAMEN( 2, C2, 'SP' ) ) THEN
+      ELSE IF( LSAMEN( 2, C2, 'PO' ) .OR. LSAMEN( 2, C2, 'PP' ) ) THEN
 *
-*        xPO, xPP, xHE, xHP, xSY, xSP: Set parameters to generate a
-*        symmetric or Hermitian matrix.
+*        xPO, xPP: Set parameters to generate a
+*        symmetric or Hermitian positive definite matrix.
 *
 *        Set TYPE, the type of matrix to be generated.
 *
             ANORM = ONE
          END IF
 *
+      ELSE IF( LSAMEN( 2, C2, 'HE' ) .OR. LSAMEN( 2, C2, 'HP' ) .OR.
+     $         LSAMEN( 2, C2, 'SY' ) .OR. LSAMEN( 2, C2, 'SP' ) ) THEN
+*
+*        xHE, xHP, xSY, xSP: Set parameters to generate a
+*        symmetric or Hermitian matrix.
+*
+*        Set TYPE, the type of matrix to be generated.
+*
+         TYPE = C2( 1: 1 )
+*
+*        Set the lower and upper bandwidths.
+*
+         IF( IMAT.EQ.1 ) THEN
+            KL = 0
+         ELSE
+            KL = MAX( N-1, 0 )
+         END IF
+         KU = KL
+*
+*        Set the condition number and norm.
+*
+         IF( IMAT.EQ.7 ) THEN
+            CNDNUM = BADC1
+         ELSE IF( IMAT.EQ.8 ) THEN
+            CNDNUM = BADC2
+         ELSE
+            CNDNUM = TWO
+         END IF
+*
+         IF( IMAT.EQ.9 ) THEN
+            ANORM = SMALL
+         ELSE IF( IMAT.EQ.10 ) THEN
+            ANORM = LARGE
+         ELSE
+            ANORM = ONE
+         END IF
+*
       ELSE IF( LSAMEN( 2, C2, 'PB' ) ) THEN
 *
 *        xPB:  Set parameters to generate a symmetric band matrix.
index e1638b0..7eb1fa4 100644 (file)
             ANORM = ONE
          END IF
 *
-      ELSE IF( LSAMEN( 2, C2, 'PO' ) .OR. LSAMEN( 2, C2, 'PP' ) .OR.
-     $         LSAMEN( 2, C2, 'SY' ) .OR. LSAMEN( 2, C2, 'SP' ) ) THEN
+      ELSE IF( LSAMEN( 2, C2, 'PO' ) .OR. LSAMEN( 2, C2, 'PP' ) ) THEN
 *
-*        xPO, xPP, xSY, xSP: Set parameters to generate a
-*        symmetric matrix.
+*        xPO, xPP: Set parameters to generate a
+*        symmetric positive definite matrix.
 *
 *        Set TYPE, the type of matrix to be generated.
 *
             ANORM = ONE
          END IF
 *
+*
+      ELSE IF( LSAMEN( 2, C2, 'SY' ) .OR. LSAMEN( 2, C2, 'SP' ) ) THEN
+*
+*        xSY, xSP: Set parameters to generate a
+*        symmetric matrix.
+*
+*        Set TYPE, the type of matrix to be generated.
+*
+         TYPE = C2( 1: 1 )
+*
+*        Set the lower and upper bandwidths.
+*
+         IF( IMAT.EQ.1 ) THEN
+            KL = 0
+         ELSE
+            KL = MAX( N-1, 0 )
+         END IF
+         KU = KL
+*
+*        Set the condition number and norm.
+*
+         IF( IMAT.EQ.7 ) THEN
+            CNDNUM = BADC1
+         ELSE IF( IMAT.EQ.8 ) THEN
+            CNDNUM = BADC2
+         ELSE
+            CNDNUM = TWO
+         END IF
+*
+         IF( IMAT.EQ.9 ) THEN
+            ANORM = SMALL
+         ELSE IF( IMAT.EQ.10 ) THEN
+            ANORM = LARGE
+         ELSE
+            ANORM = ONE
+         END IF
+*
       ELSE IF( LSAMEN( 2, C2, 'PB' ) ) THEN
 *
 *        xPB:  Set parameters to generate a symmetric band matrix.
index 4ead31d..145c65d 100644 (file)
             ANORM = ONE
          END IF
 *
-      ELSE IF( LSAMEN( 2, C2, 'PO' ) .OR. LSAMEN( 2, C2, 'PP' ) .OR.
-     $         LSAMEN( 2, C2, 'SY' ) .OR. LSAMEN( 2, C2, 'SP' ) ) THEN
+      ELSE IF( LSAMEN( 2, C2, 'PO' ) .OR. LSAMEN( 2, C2, 'PP' ) ) THEN
 *
 *        xPO, xPP, xSY, xSP: Set parameters to generate a
-*        symmetric matrix.
+*        symmetric positive definite matrix.
 *
 *        Set TYPE, the type of matrix to be generated.
 *
             ANORM = ONE
          END IF
 *
+*
+      ELSE IF( LSAMEN( 2, C2, 'SY' ) .OR. LSAMEN( 2, C2, 'SP' ) ) THEN
+*
+*        xSY, xSP: Set parameters to generate a
+*        symmetric matrix.
+*
+*        Set TYPE, the type of matrix to be generated.
+*
+         TYPE = C2( 1: 1 )
+*
+*        Set the lower and upper bandwidths.
+*
+         IF( IMAT.EQ.1 ) THEN
+            KL = 0
+         ELSE
+            KL = MAX( N-1, 0 )
+         END IF
+         KU = KL
+*
+*        Set the condition number and norm.
+*
+         IF( IMAT.EQ.7 ) THEN
+            CNDNUM = BADC1
+         ELSE IF( IMAT.EQ.8 ) THEN
+            CNDNUM = BADC2
+         ELSE
+            CNDNUM = TWO
+         END IF
+*
+         IF( IMAT.EQ.9 ) THEN
+            ANORM = SMALL
+         ELSE IF( IMAT.EQ.10 ) THEN
+            ANORM = LARGE
+         ELSE
+            ANORM = ONE
+         END IF
+*
       ELSE IF( LSAMEN( 2, C2, 'PB' ) ) THEN
 *
 *        xPB:  Set parameters to generate a symmetric band matrix.
index bcfa234..de26637 100644 (file)
             ANORM = ONE
          END IF
 *
-      ELSE IF( LSAMEN( 2, C2, 'PO' ) .OR. LSAMEN( 2, C2, 'PP' ) .OR.
-     $         LSAMEN( 2, C2, 'HE' ) .OR. LSAMEN( 2, C2, 'HP' ) .OR.
-     $         LSAMEN( 2, C2, 'SY' ) .OR. LSAMEN( 2, C2, 'SP' ) ) THEN
+      ELSE IF( LSAMEN( 2, C2, 'PO' ) .OR. LSAMEN( 2, C2, 'PP' ) ) THEN
 *
-*        xPO, xPP, xHE, xHP, xSY, xSP: Set parameters to generate a
-*        symmetric or Hermitian matrix.
+*        xPO, xPP: Set parameters to generate a
+*        symmetric or Hermitian positive definite matrix.
 *
 *        Set TYPE, the type of matrix to be generated.
 *
             ANORM = ONE
          END IF
 *
+      ELSE IF( LSAMEN( 2, C2, 'HE' ) .OR. LSAMEN( 2, C2, 'HP' ) .OR.
+     $         LSAMEN( 2, C2, 'SY' ) .OR. LSAMEN( 2, C2, 'SP' ) ) THEN
+*
+*        xHE, xHP, xSY, xSP: Set parameters to generate a
+*        symmetric or Hermitian matrix.
+*
+*        Set TYPE, the type of matrix to be generated.
+*
+         TYPE = C2( 1: 1 )
+*
+*        Set the lower and upper bandwidths.
+*
+         IF( IMAT.EQ.1 ) THEN
+            KL = 0
+         ELSE
+            KL = MAX( N-1, 0 )
+         END IF
+         KU = KL
+*
+*        Set the condition number and norm.
+*
+         IF( IMAT.EQ.7 ) THEN
+            CNDNUM = BADC1
+         ELSE IF( IMAT.EQ.8 ) THEN
+            CNDNUM = BADC2
+         ELSE
+            CNDNUM = TWO
+         END IF
+*
+         IF( IMAT.EQ.9 ) THEN
+            ANORM = SMALL
+         ELSE IF( IMAT.EQ.10 ) THEN
+            ANORM = LARGE
+         ELSE
+            ANORM = ONE
+         END IF
+*
       ELSE IF( LSAMEN( 2, C2, 'PB' ) ) THEN
 *
 *        xPB:  Set parameters to generate a symmetric band matrix.