Skip to content

Commit 0c40422

Browse files
committed
Testers for CSY two-stage Aasen's
1 parent 852827d commit 0c40422

7 files changed

Lines changed: 1185 additions & 4 deletions

File tree

TESTING/LIN/CMakeLists.txt

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -57,11 +57,11 @@ set(CLINTST cchkaa.f
5757
cchkhp.f cchklq.f cchkpb.f
5858
cchkpo.f cchkps.f cchkpp.f cchkpt.f cchkq3.f cchkql.f
5959
cchkqr.f cchkrq.f cchksp.f cchksy.f cchksy_rook.f cchksy_rk.f
60-
cchksy_aa.f
60+
cchksy_aa.f cchksy_aa_2stage.f
6161
cchktb.f
6262
cchktp.f cchktr.f cchktz.f
6363
cdrvgt.f cdrvhe_rook.f cdrvhe_rk.f
64-
cdrvhe_aa.f cdrvhe_aa_2stage.f
64+
cdrvhe_aa.f cdrvhe_aa_2stage.f cdrvsy_aa_2stage.f
6565
cdrvhp.f
6666
cdrvls.f cdrvpb.f cdrvpp.f cdrvpt.f
6767
cdrvsp.f cdrvsy_rook.f cdrvsy_rk.f

TESTING/LIN/Makefile

Lines changed: 3 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -94,11 +94,12 @@ CLINTST = cchkaa.o \
9494
cchkhe_aa.o cchkhe_aa_2stage.o cchkhp.o cchklq.o cchkpb.o \
9595
cchkpo.o cchkps.o cchkpp.o cchkpt.o cchkq3.o cchkql.o \
9696
cchkqr.o cchkrq.o cchksp.o cchksy.o cchksy_rook.o cchksy_rk.o \
97-
cchksy_aa.o cdrvhe_aa_2stage.o cchktb.o \
97+
cchksy_aa.o cchksy_aa_2stage.o cchktb.o \
9898
cchktp.o cchktr.o cchktz.o \
9999
cdrvgt.o cdrvhe_rook.o cdrvhe_rk.o cdrvhe_aa.o cdrvhp.o \
100+
cdrvhe_aa_2stage.o \
100101
cdrvls.o cdrvpb.o cdrvpp.o cdrvpt.o \
101-
cdrvsp.o cdrvsy_rook.o cdrvsy_rk.o cdrvsy_aa.o \
102+
cdrvsp.o cdrvsy_rook.o cdrvsy_rk.o cdrvsy_aa.o cdrvsy_aa_2stage.o \
102103
cerrgt.o cerrlq.o \
103104
cerrls.o cerrps.o cerrql.o cerrqp.o \
104105
cerrqr.o cerrrq.o cerrtr.o cerrtz.o \

TESTING/LIN/cchkaa.f

Lines changed: 27 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -882,6 +882,33 @@ PROGRAM CCHKAA
882882
ELSE
883883
WRITE( NOUT, FMT = 9988 )PATH
884884
END IF
885+
*
886+
ELSE IF( LSAMEN( 2, C2, 'S2' ) ) THEN
887+
*
888+
* SA: symmetric indefinite matrices with Aasen's algorithm,
889+
*
890+
NTYPES = 11
891+
CALL ALAREQ( PATH, NMATS, DOTYPE, NTYPES, NIN, NOUT )
892+
*
893+
IF( TSTCHK ) THEN
894+
CALL CCHKSY_AA_2STAGE( DOTYPE, NN, NVAL, NNB2, NBVAL2, NNS,
895+
$ NSVAL, THRESH, TSTERR, LDA,
896+
$ A( 1, 1 ), A( 1, 2 ), A( 1, 3 ),
897+
$ B( 1, 1 ), B( 1, 2 ), B( 1, 3 ),
898+
$ WORK, RWORK, IWORK, NOUT )
899+
ELSE
900+
WRITE( NOUT, FMT = 9989 )PATH
901+
END IF
902+
*
903+
IF( TSTDRV ) THEN
904+
CALL CDRVSY_AA_2STAGE(
905+
$ DOTYPE, NN, NVAL, NRHS, THRESH, TSTERR,
906+
$ LDA, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ),
907+
$ B( 1, 1 ), B( 1, 2 ), B( 1, 3 ), WORK,
908+
$ RWORK, IWORK, NOUT )
909+
ELSE
910+
WRITE( NOUT, FMT = 9988 )PATH
911+
END IF
885912
*
886913
ELSE IF( LSAMEN( 2, C2, 'SP' ) ) THEN
887914
*

0 commit comments

Comments
 (0)