diff options
Diffstat (limited to 'TESTING/EIG/zchkst2stg.f')
| -rw-r--r-- | TESTING/EIG/zchkst2stg.f | 58 |
1 files changed, 3 insertions, 55 deletions
diff --git a/TESTING/EIG/zchkst2stg.f b/TESTING/EIG/zchkst2stg.f index a1aaffbc..d93191a6 100644 --- a/TESTING/EIG/zchkst2stg.f +++ b/TESTING/EIG/zchkst2stg.f @@ -1,7 +1,5 @@ *> \brief \b ZCHKST2STG * -* @precisions fortran z -> c -* * =========== DOCUMENTATION =========== * * Online html documentation available at @@ -616,7 +614,7 @@ *> \author Univ. of Colorado Denver *> \author NAG Ltd. * -*> \date November 2011 +*> \date November 2016 * *> \ingroup complex16_eig * @@ -627,10 +625,10 @@ $ LWORK, RWORK, LRWORK, IWORK, LIWORK, RESULT, $ INFO ) * -* -- LAPACK test routine (version 3.4.0) -- +* -- LAPACK test routine (version 3.7.0) -- * -- LAPACK is a software package provided by Univ. of Tennessee, -- * -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- -* November 2011 +* November 2016 * * .. Scalar Arguments .. INTEGER INFO, LDA, LDU, LIWORK, LRWORK, LWORK, NOUNIT, @@ -1097,56 +1095,6 @@ RESULT( 3 ) = TEMP2 / MAX( UNFL, ULP*MAX( TEMP1, TEMP2 ) ) RESULT( 4 ) = TEMP4 / MAX( UNFL, ULP*MAX( TEMP3, TEMP4 ) ) * -* Skip the DSYTRD for lower that since we replaced its testing -* 3 and 4 by the 2-stage one. - GOTO 101 -* -* Call ZHETRD and ZUNGTR to compute S and U from -* lower triangle, do tests. -* - CALL ZLACPY( 'L', N, N, A, LDA, V, LDU ) -* - NTEST = 3 - CALL ZHETRD( 'L', N, V, LDU, SD, SE, TAU, WORK, LWORK, - $ IINFO ) -* - IF( IINFO.NE.0 ) THEN - WRITE( NOUNIT, FMT = 9999 )'ZHETRD(L)', IINFO, N, JTYPE, - $ IOLDSD - INFO = ABS( IINFO ) - IF( IINFO.LT.0 ) THEN - RETURN - ELSE - RESULT( 3 ) = ULPINV - GO TO 280 - END IF - END IF -* - CALL ZLACPY( 'L', N, N, V, LDU, U, LDU ) -* - NTEST = 4 - CALL ZUNGTR( 'L', N, U, LDU, TAU, WORK, LWORK, IINFO ) - IF( IINFO.NE.0 ) THEN - WRITE( NOUNIT, FMT = 9999 )'ZUNGTR(L)', IINFO, N, JTYPE, - $ IOLDSD - INFO = ABS( IINFO ) - IF( IINFO.LT.0 ) THEN - RETURN - ELSE - RESULT( 4 ) = ULPINV - GO TO 280 - END IF - END IF -* - CALL ZHET21( 2, 'Lower', N, 1, A, LDA, SD, SE, U, LDU, V, - $ LDU, TAU, WORK, RWORK, RESULT( 3 ) ) - CALL ZHET21( 3, 'Lower', N, 1, A, LDA, SD, SE, U, LDU, V, - $ LDU, TAU, WORK, RWORK, RESULT( 4 ) ) -* -*after skipping old tests 3 4 back to the normal -* - 101 CONTINUE -* * Store the upper triangle of A in AP * I = 0 |