aboutsummaryrefslogtreecommitdiff
path: root/TESTING/EIG/schkst2stg.f
diff options
context:
space:
mode:
authorJulie <julie@cs.utk.edu>2016-11-19 17:09:13 -0800
committerJulie <julie@cs.utk.edu>2016-11-19 17:09:13 -0800
commited4e95e4d25f5f995db75445cbc2ba1db5736d23 (patch)
tree0b68f343e62958a8e71dc34c084e9e899eec990d /TESTING/EIG/schkst2stg.f
parent01cdfedf1d726a003f7a4e7331f32a7e434f1707 (diff)
First minor change to 2 stage contribution
Update Makefile to handle .F Update Cmake build Note: We have duplicate loop numbers due to Define statement that stops the compilation with Cmake default option. Fix coming very shortly
Diffstat (limited to 'TESTING/EIG/schkst2stg.f')
-rw-r--r--TESTING/EIG/schkst2stg.f58
1 files changed, 3 insertions, 55 deletions
diff --git a/TESTING/EIG/schkst2stg.f b/TESTING/EIG/schkst2stg.f
index 8db1cf73..4bf9f107 100644
--- a/TESTING/EIG/schkst2stg.f
+++ b/TESTING/EIG/schkst2stg.f
@@ -1,7 +1,5 @@
*> \brief \b SCHKST2STG
*
-* @generated from dchkst2stg.f, fortran d -> s, Sat Nov 5 22:51:30 2016
-*
* =========== DOCUMENTATION ===========
*
* Online html documentation available at
@@ -604,7 +602,7 @@
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date November 2016
*
*> \ingroup single_eig
*
@@ -614,10 +612,10 @@
$ WA1, WA2, WA3, WR, U, LDU, V, VP, TAU, Z, WORK,
$ LWORK, 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, LWORK, NOUNIT, NSIZES,
@@ -1080,56 +1078,6 @@
RESULT( 3 ) = TEMP2 / MAX( UNFL, ULP*MAX( TEMP1, TEMP2 ) )
RESULT( 4 ) = TEMP4 / MAX( UNFL, ULP*MAX( TEMP3, TEMP4 ) )
*
-* Skip the SSYTRD for lower that since we replaced its testing
-* 3 and 4 by the 2-stage one.
- GOTO 101
-*
-* Call SSYTRD and SORGTR to compute S and U from
-* lower triangle, do tests.
-*
- CALL SLACPY( 'L', N, N, A, LDA, V, LDU )
-*
- NTEST = 3
- CALL SSYTRD( 'L', N, V, LDU, SD, SE, TAU, WORK, LWORK,
- $ IINFO )
-*
- IF( IINFO.NE.0 ) THEN
- WRITE( NOUNIT, FMT = 9999 )'SSYTRD(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 SLACPY( 'L', N, N, V, LDU, U, LDU )
-*
- NTEST = 4
- CALL SORGTR( 'L', N, U, LDU, TAU, WORK, LWORK, IINFO )
- IF( IINFO.NE.0 ) THEN
- WRITE( NOUNIT, FMT = 9999 )'SORGTR(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 SSYT21( 2, 'Lower', N, 1, A, LDA, SD, SE, U, LDU, V,
- $ LDU, TAU, WORK, RESULT( 3 ) )
- CALL SSYT21( 3, 'Lower', N, 1, A, LDA, SD, SE, U, LDU, V,
- $ LDU, TAU, WORK, RESULT( 4 ) )
-*
-*after skipping old tests 3 4 back to the normal
-*
- 101 CONTINUE
-*
* Store the upper triangle of A in AP
*
I = 0