Just C precision for testing

Author: julielangou

SRC/zhetrf_aa_2stage.f

@@ -53,6 +53,7 @@
 *> CHR   10               List types on next line if 0 < NTYPES < 10
 *> CHK   10               List types on next line if 0 < NTYPES < 10
 *> CHA   10               List types on next line if 0 < NTYPES < 10
+*> CH2   10               List types on next line if 0 < NTYPES < 10
 *> CHP   10               List types on next line if 0 < NTYPES < 10
 *> CSY   11               List types on next line if 0 < NTYPES < 11
 *> CSK   11               List types on next line if 0 < NTYPES < 11

TESTING/LIN/cchkaa.f

@@ -14,15 +14,16 @@
 *                             IWORK, NOUT )
 *
 *       .. Scalar Arguments ..
-*       LOGICAL            TSTERR
-*       INTEGER            NMAX, NN, NNB, NNS, NOUT
-*       COMPLEX            THRESH
+*       LOGICAL    TSTERR
+*       INTEGER    NN, NNB, NNS, NOUT
+*       REAL       THRESH
 *       ..
 *       .. Array Arguments ..
-*       LOGICAL            DOTYPE( * )
-*       INTEGER            IWORK( * ), NBVAL( * ), NSVAL( * ), NVAL( * )
-*       COMPLEX            A( * ), AFAC( * ), AINV( * ), B( * ),
-*      $                   RWORK( * ), WORK( * ), X( * ), XACT( * )
+*       LOGICAL    DOTYPE( * )
+*       INTEGER    IWORK( * ), NBVAL( * ), NSVAL( * ), NVAL( * )
+*       REAL       RWORK( * )
+*       COMPLEX    A( * ), AFAC( * ), AINV( * ), B( * ),
+*      $           WORK( * ), X( * ), XACT( * )
 *       ..
 *
 *
@@ -31,7 +32,7 @@
 *>
 *> \verbatim
 *>
-*> DCHKSY_AA_2STAGE tests CHETRF_AA_2STAGE, -TRS_AA_2STAGE.
+*> CCHKSY_AA_2STAGE tests CHETRF_AA_2STAGE, -TRS_AA_2STAGE.
 *> \endverbatim
 *
 *  Arguments:
@@ -83,7 +84,7 @@
 *>
 *> \param[in] THRESH
 *> \verbatim
-*>          THRESH is COMPLEX
+*>          THRESH is REAL
 *>          The threshold value for the test ratios.  A result is
 *>          included in the output file if RESULT >= THRESH.  To have
 *>          every test ratio printed, use THRESH = 0.
@@ -140,7 +141,7 @@
 *>
 *> \param[out] RWORK
 *> \verbatim
-*>          RWORK is COMPLEX array, dimension (max(NMAX,2*NSMAX))
+*>          RWORK is REAL array, dimension (max(NMAX,2*NSMAX))
 *> \endverbatim
 *>
 *> \param[out] IWORK
@@ -162,7 +163,8 @@
 *> \author Univ. of Colorado Denver
 *> \author NAG Ltd.
 *
-*> \date June 2017
+*> \date November 2017
+*
 *
 *> \ingroup complex_lin
 *
@@ -184,10 +186,11 @@ SUBROUTINE CCHKHE_AA_2STAGE( DOTYPE, NN, NVAL, NNB, NBVAL, NNS,
       REAL               THRESH
 *     ..
 *     .. Array Arguments ..
-      LOGICAL            DOTYPE( * )
-      INTEGER            IWORK( * ), NBVAL( * ), NSVAL( * ), NVAL( * )
-      COMPLEX            A( * ), AFAC( * ), AINV( * ), B( * ),
-     $                   RWORK( * ), WORK( * ), X( * ), XACT( * )
+      LOGICAL      DOTYPE( * )
+      INTEGER      IWORK( * ), NBVAL( * ), NSVAL( * ), NVAL( * )
+      REAL         RWORK( * )
+      COMPLEX      A( * ), AFAC( * ), AINV( * ), B( * ),
+     $             WORK( * ), X( * ), XACT( * )
 *     ..
 *
 *  =====================================================================
@@ -208,18 +211,18 @@ SUBROUTINE CCHKHE_AA_2STAGE( DOTYPE, NN, NVAL, NNB, NBVAL, NNS,
       INTEGER            I, I1, I2, IMAT, IN, INB, INFO, IOFF, IRHS,
      $                   IUPLO, IZERO, J, K, KL, KU, LDA, LWORK, MODE,
      $                   N, NB, NERRS, NFAIL, NIMAT, NRHS, NRUN, NT
-      COMPLEX            ANORM, CNDNUM
+      REAL               ANORM, CNDNUM
 *     ..
 *     .. Local Arrays ..
       CHARACTER          UPLOS( 2 )
       INTEGER            ISEED( 4 ), ISEEDY( 4 )
       REAL               RESULT( NTESTS )
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           ALAERH, ALAHD, ALASUM, CERRHE, CLACPY, CLARHS,
-     $                   CLATB4, CLATMS, CPOT02, DSYT01, 
-     $                   CHETRF_AA_2STAGE, CHETRS_AA_2STAGE,
-     $                   XLAENV
+      EXTERNAL           ALAERH, ALAHD, ALASUM, CERRHE, CLACPY, 
+     $                   CLARHS, CLATB4, CLATMS, CPOT02, 
+     $                   CHETRF_AA_2STAGE, 
+     $                   CHETRS_AA_2STAGE, XLAENV
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          MAX, MIN
@@ -241,6 +244,7 @@ SUBROUTINE CCHKHE_AA_2STAGE( DOTYPE, NN, NVAL, NNB, NBVAL, NNS,
 *
 *     Initialize constants and the random number seed.
 *
+*
 *     Test path
 *
       PATH( 1: 1 ) = 'Complex precision'
@@ -464,10 +468,12 @@ SUBROUTINE CCHKHE_AA_2STAGE( DOTYPE, NN, NVAL, NNB, NBVAL, NNS,
 *+    TEST 1
 *                 Reconstruct matrix from factors and compute residual.
 *
-c                  CALL DSYT01_AA( UPLO, N, A, LDA, AFAC, LDA, IWORK,
+*                   
+c                 NEED TO WRITE CHET01_AA_2STAGE
+c                  CALL CHET01_AA( UPLO, N, A, LDA, AFAC, LDA, IWORK,
 c     $                            AINV, LDA, RWORK, RESULT( 1 ) )
 c                  NT = 1
-                  NT = 0
+                   NT = 0
 *
 *
 *                 Print information about the tests that did not pass
@@ -530,7 +536,6 @@ SUBROUTINE CCHKHE_AA_2STAGE( DOTYPE, NN, NVAL, NNB, NBVAL, NNS,
                         CALL CPOT02( UPLO, N, NRHS, A, LDA, X, LDA,
      $                               WORK, LDA, RWORK, RESULT( 2 ) )
 *
-*
 *                       Print information about the tests that did not pass
 *                       the threshold.
 *
@@ -567,6 +572,6 @@ SUBROUTINE CCHKHE_AA_2STAGE( DOTYPE, NN, NVAL, NNB, NBVAL, NNS,
      $      I6 )
       RETURN
 *
-*     End of DCHKSY_AA_2STAGE
+*     End of CCHKSY_AA_2STAGE
 *
       END

TESTING/LIN/cchkhe_aa_2stage.f

@@ -204,7 +204,7 @@ SUBROUTINE CDRVHE_AA_2STAGE(
 *     .. External Subroutines ..
       EXTERNAL           ALADHD, ALAERH, ALASVM, XLAENV, CERRVX,
      $                   CGET04, CLACPY, CLARHS, CLATB4, CLATMS,
-     $                   CHESV_AA_2STAGE, CHET01_AA, CPOT02,
+     $                   CHESV_AA_2STAGE, CPOT02,
      $                   CHETRF_AA_2STAGE
 *     ..
 *     .. Scalars in Common ..
@@ -453,7 +453,7 @@ SUBROUTINE CDRVHE_AA_2STAGE(
 c     $                                  IWORK, AINV, LDA, RWORK,
 c     $                                  RESULT( 2 ) )
 c                     NT = 2
-                     NT = 1
+					  NT = 1
 *
 *                    Print information about the tests that did not pass
 *                    the threshold.
@@ -485,6 +485,6 @@ SUBROUTINE CDRVHE_AA_2STAGE(
      $      ', test ', I2, ', ratio =', G12.5 )
       RETURN
 *
-*     End of CDRVHE_AA
+*     End of CDRVHE_AA_2STAGE
 *
       END

TESTING/LIN/cdrvhe_aa_2stage.f

@@ -53,6 +53,7 @@
 *> DSR   10               List types on next line if 0 < NTYPES < 10
 *> DSK   10               List types on next line if 0 < NTYPES < 10
 *> DSA   10               List types on next line if 0 < NTYPES < 10
+*> DS2   10               List types on next line if 0 < NTYPES < 10
 *> DSP   10               List types on next line if 0 < NTYPES < 10
 *> DTR   18               List types on next line if 0 < NTYPES < 18
 *> DTP   18               List types on next line if 0 < NTYPES < 18

TESTING/LIN/dchkaa.f

@@ -53,6 +53,7 @@
 *> SSR   10               List types on next line if 0 < NTYPES < 10
 *> SSK   10               List types on next line if 0 < NTYPES < 10
 *> SSA   10               List types on next line if 0 < NTYPES < 10
+*> SS2   10               List types on next line if 0 < NTYPES < 10
 *> SSP   10               List types on next line if 0 < NTYPES < 10
 *> STR   18               List types on next line if 0 < NTYPES < 18
 *> STP   18               List types on next line if 0 < NTYPES < 18

TESTING/LIN/schkaa.f

@@ -218,7 +218,7 @@ SUBROUTINE SCHKSY_AA_2STAGE( DOTYPE, NN, NVAL, NNB, NBVAL, NNS,
 *     ..
 *     .. External Subroutines ..
       EXTERNAL           ALAERH, ALAHD, ALASUM, SERRSY, SLACPY, SLARHS,
-     $                   SLATB4, SLATMS, SPOT02, SSYT01, 
+     $                   SLATB4, SLATMS, SPOT02, SSYT01_AA, 
      $                   SSYTRF_AA_2STAGE, SSYTRS_AA_2STAGE,
      $                   XLAENV
 *     ..
@@ -465,10 +465,9 @@ SUBROUTINE SCHKSY_AA_2STAGE( DOTYPE, NN, NVAL, NNB, NBVAL, NNS,
 *+    TEST 1
 *                 Reconstruct matrix from factors and compute residual.
 *
-c                  CALL SSYT01_AA( UPLO, N, A, LDA, AFAC, LDA, IWORK,
-c     $                            AINV, LDA, RWORK, RESULT( 1 ) )
-c                  NT = 1
-                  NT = 0
+                  CALL SSYT01_AA( UPLO, N, A, LDA, AFAC, LDA, IWORK,
+     $                            AINV, LDA, RWORK, RESULT( 1 ) )
+                  NT = 1
 *
 *
 *                 Print information about the tests that did not pass