diff options
Diffstat (limited to 'TESTING/LIN/dlatb4.f')
| -rw-r--r-- | TESTING/LIN/dlatb4.f | 44 |
1 files changed, 40 insertions, 4 deletions
diff --git a/TESTING/LIN/dlatb4.f b/TESTING/LIN/dlatb4.f index e1638b08..7eb1fa43 100644 --- a/TESTING/LIN/dlatb4.f +++ b/TESTING/LIN/dlatb4.f @@ -339,11 +339,10 @@ 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. * @@ -376,6 +375,43 @@ 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. |