Merge pull request #86 from iyamazaki/lapack-aasen

Lapack aasen

12 files changed, 196 insertions, 46 deletions

diff --git a/TESTING/LIN/cchkaa.f b/TESTING/LIN/cchkaa.f
index c5dd3330..cffaa1d6 100644
--- a/TESTING/LIN/cchkaa.f
+++ b/TESTING/LIN/cchkaa.f
@@ -160,13 +160,13 @@
 *     ..
 *     .. External Subroutines ..
       EXTERNAL           ALAREQ, CCHKEQ, CCHKGB, CCHKGE, CCHKGT, CCHKHE,
-     $                   CCHKHE_ROOK, CCHKHP, CCHKLQ, CCHKPB, CCHKPO,
-     $                   CCHKPS, CCHKPP, CCHKPT, CCHKQ3, CCHKQL,
+     $                   CCHKHE_ROOK, CCHKHE_AA, CCHKHP, CCHKLQ, CCHKPB,
+     $                   CCHKPO, CCHKPS, CCHKPP, CCHKPT, CCHKQ3, CCHKQL,
      $                   CCHKQR, CCHKRQ, CCHKSP, CCHKSY, CCHKSY_ROOK,
      $                   CCHKTB, CCHKTP, CCHKTR, CCHKTZ, CDRVGB, CDRVGE,
-     $                   CDRVGT, CDRVHE, CDRVHE_ROOK, CDRVHP, CDRVLS,
-     $                   CDRVPB, CDRVPO, CDRVPP, CDRVPT, CDRVSP, CDRVSY,
-     $                   CDRVSY_ROOK, ILAVER, CCHKQRT, CCHKQRTP
+     $                   CDRVGT, CDRVHE, CDRVHE_ROOK, CDRVHE_AA, CDRVHP,
+     $                   CDRVLS, CDRVPB, CDRVPO, CDRVPP, CDRVPT, CDRVSP,
+     $                   CDRVSY, CDRVSY_ROOK, ILAVER, CCHKQRT, CCHKQRTP
 
 *     ..
 *     .. Scalars in Common ..
diff --git a/TESTING/LIN/cchkhe_aa.f b/TESTING/LIN/cchkhe_aa.f
index 9666f8c2..702677bd 100644
--- a/TESTING/LIN/cchkhe_aa.f
+++ b/TESTING/LIN/cchkhe_aa.f
@@ -9,8 +9,8 @@
 *  ===========
 *
 *       SUBROUTINE CCHKHE_AA( DOTYPE, NN, NVAL, NNB, NBVAL, NNS, NSVAL,
-*                                THRESH, TSTERR, NMAX, A, AFAC, AINV, B, X,
-*                                XACT, WORK, RWORK, IWORK, NOUT )
+*                             THRESH, TSTERR, NMAX, A, AFAC, AINV, B, X,
+*                             XACT, WORK, RWORK, IWORK, NOUT )
 *
 *       .. Scalar Arguments ..
 *       LOGICAL    TSTERR
@@ -169,8 +169,8 @@
 *
 *  =====================================================================
       SUBROUTINE CCHKHE_AA( DOTYPE, NN, NVAL, NNB, NBVAL, NNS, NSVAL,
-     $                         THRESH, TSTERR, NMAX, A, AFAC, AINV, B,
-     $                         X, XACT, WORK, RWORK, IWORK, NOUT )
+     $                      THRESH, TSTERR, NMAX, A, AFAC, AINV, B,
+     $                      X, XACT, WORK, RWORK, IWORK, NOUT )
 *
 *  -- LAPACK test routine (version 3.7.0) --
 *  -- LAPACK is a software package provided by Univ. of Tennessee,    --
diff --git a/TESTING/LIN/cerrhe.f b/TESTING/LIN/cerrhe.f
index 19a3182f..22defe6e 100644
--- a/TESTING/LIN/cerrhe.f
+++ b/TESTING/LIN/cerrhe.f
@@ -89,9 +89,10 @@
 *     ..
 *     .. External Subroutines ..
       EXTERNAL           ALAESM, CHECON, CHECON_ROOK, CHERFS, CHETF2,
-     $                   CHETF2_ROOK, CHETRF, CHETRF_ROOK, CHETRI,
-     $                   CHETRI_ROOK, CHETRI2, CHETRS, CHETRS_ROOK,
-     $                   CHKXER, CHPCON, CHPRFS, CHPTRF, CHPTRI, CHPTRS
+     $                   CHETF2_ROOK, CHETRF, CHETRF_ROOK, CHETRF_AA, 
+     $                   CHETRI, CHETRI_ROOK, CHETRI2, CHETRS, 
+     $                   CHETRS_ROOK, CHETRS_AA, CHKXER, CHPCON, CHPRFS,
+     $                   CHPTRF, CHPTRI, CHPTRS
 *     ..
 *     .. Scalars in Common ..
       LOGICAL            LERR, OK
@@ -334,6 +335,43 @@
          CALL CHKXER( 'CHECON_ROOK', INFOT, NOUT, LERR, OK )
 *
 *        Test error exits of the routines that use factorization
+*        of a Hermitian indefinite matrix with Aasen's algorithm.
+*
+      ELSE IF( LSAMEN( 2, C2, 'HA' ) ) THEN
+*
+*        CHETRF_AA
+*
+         SRNAMT = 'CHETRF_AA'
+         INFOT = 1
+         CALL CHETRF_AA( '/', 0, A, 1, IP, W, 1, INFO )
+         CALL CHKXER( 'CHETRF_AA', INFOT, NOUT, LERR, OK )
+         INFOT = 2
+         CALL CHETRF_AA( 'U', -1, A, 1, IP, W, 1, INFO )
+         CALL CHKXER( 'CHETRF_AA', INFOT, NOUT, LERR, OK )
+         INFOT = 4
+         CALL CHETRF_AA( 'U', 2, A, 1, IP, W, 4, INFO )
+         CALL CHKXER( 'CHETRF_AA', INFOT, NOUT, LERR, OK )
+*
+*        CHETRS_AA
+*
+         SRNAMT = 'CHETRS_AA'
+         INFOT = 1
+         CALL CHETRS_AA( '/', 0, 0, A, 1, IP, B, 1, W, 1, INFO )
+         CALL CHKXER( 'CHETRS_AA', INFOT, NOUT, LERR, OK )
+         INFOT = 2
+         CALL CHETRS_AA( 'U', -1, 0, A, 1, IP, B, 1, W, 1, INFO )
+         CALL CHKXER( 'CHETRS_AA', INFOT, NOUT, LERR, OK )
+         INFOT = 3
+         CALL CHETRS_AA( 'U', 0, -1, A, 1, IP, B, 1, W, 1, INFO )
+         CALL CHKXER( 'CHETRS_AA', INFOT, NOUT, LERR, OK )
+         INFOT = 5
+         CALL CHETRS_AA( 'U', 2, 1, A, 1, IP, B, 2, W, 1, INFO )
+         CALL CHKXER( 'CHETRS_AA', INFOT, NOUT, LERR, OK )
+         INFOT = 8
+         CALL CHETRS_AA( 'U', 2, 1, A, 2, IP, B, 1, W, 1, INFO )
+         CALL CHKXER( 'CHETRS_AA', INFOT, NOUT, LERR, OK )
+*
+*        Test error exits of the routines that use factorization
 *        of a Hermitian indefinite packed matrix with patrial
 *        (Bunch-Kaufman) diagonal pivoting method.
 *
diff --git a/TESTING/LIN/dchkaa.f b/TESTING/LIN/dchkaa.f
index 85d2b2a4..8bcb8217 100644
--- a/TESTING/LIN/dchkaa.f
+++ b/TESTING/LIN/dchkaa.f
@@ -159,11 +159,11 @@
       EXTERNAL           ALAREQ, DCHKEQ, DCHKGB, DCHKGE, DCHKGT, DCHKLQ,
      $                   DCHKPB, DCHKPO, DCHKPS, DCHKPP, DCHKPT, DCHKQ3,
      $                   DCHKQL, DCHKQR, DCHKRQ, DCHKSP, DCHKSY,
-     $                   DCHKSY_ROOK, DCHKTB, DCHKTP, DCHKTR, DCHKTZ,
-     $                   DDRVGB, DDRVGE, DDRVGT, DDRVLS, DDRVPB, DDRVPO,
-     $                   DDRVPP, DDRVPT, DDRVSP, DDRVSY, DDRVSY_ROOK,
-     $                   ILAVER, DCHKQRT, DCHKQRTP, DCHKLQTP, DCHKTSQR,
-     $                   DCHKLQT
+     $                   DCHKSY_ROOK, DCHKSY_AA, DCHKTB, DCHKTP, DCHKTR,
+     $                   DCHKTZ, DDRVGB, DDRVGE, DDRVGT, DDRVLS, DDRVPB,
+     $                   DDRVPO, DDRVPP, DDRVPT, DDRVSP, DDRVSY, 
+     $                   DDRVSY_ROOK, DDRVSY_AA, ILAVER, DCHKQRT, 
+     $                   DCHKQRTP, DCHKLQTP, DCHKTSQR, DCHKLQT
 
 *     ..
 *     .. Scalars in Common ..
diff --git a/TESTING/LIN/dchksy_aa.f b/TESTING/LIN/dchksy_aa.f
index 5bfe71c6..a596bcb9 100644
--- a/TESTING/LIN/dchksy_aa.f
+++ b/TESTING/LIN/dchksy_aa.f
@@ -9,8 +9,8 @@
 *  ===========
 *
 *       SUBROUTINE DCHKSY_AA( DOTYPE, NN, NVAL, NNB, NBVAL, NNS, NSVAL,
-*                                THRESH, TSTERR, NMAX, A, AFAC, AINV, B, X,
-*                                XACT, WORK, RWORK, IWORK, NOUT )
+*                             THRESH, TSTERR, NMAX, A, AFAC, AINV, B, X,
+*                             XACT, WORK, RWORK, IWORK, NOUT )
 *
 *       .. Scalar Arguments ..
 *       LOGICAL            TSTERR
@@ -168,8 +168,8 @@
 *
 *  =====================================================================
       SUBROUTINE DCHKSY_AA( DOTYPE, NN, NVAL, NNB, NBVAL, NNS, NSVAL,
-     $                         THRESH, TSTERR, NMAX, A, AFAC, AINV, B,
-     $                         X, XACT, WORK, RWORK, IWORK, NOUT )
+     $                      THRESH, TSTERR, NMAX, A, AFAC, AINV, B,
+     $                      X, XACT, WORK, RWORK, IWORK, NOUT )
 *
 *  -- LAPACK test routine (version 3.7.0) --
 *  -- LAPACK is a software package provided by Univ. of Tennessee,    --
diff --git a/TESTING/LIN/derrsy.f b/TESTING/LIN/derrsy.f
index 8523fc8c..a453ab19 100644
--- a/TESTING/LIN/derrsy.f
+++ b/TESTING/LIN/derrsy.f
@@ -88,8 +88,9 @@
 *     .. External Subroutines ..
       EXTERNAL           ALAESM, CHKXER, DSPCON, DSPRFS, DSPTRF, DSPTRI,
      $                   DSPTRS, DSYCON, DSYCON_ROOK, DSYRFS, DSYTF2,
-     $                   DSYTF2_ROOK, DSYTRF, DSYTRF_ROOK, DSYTRI,
-     $                   DSYTRI_ROOK, DSYTRI2, DSYTRS, DSYTRS_ROOK
+     $                   DSYTF2_ROOK, DSYTRF, DSYTRF_ROOK, DSYTRF_AA,
+     $                   DSYTRI, DSYTRI_ROOK, DSYTRI2, DSYTRS, 
+     $                   DSYTRS_ROOK, DSYTRS_AA
 *     ..
 *     .. Scalars in Common ..
       LOGICAL            LERR, OK
@@ -333,6 +334,42 @@
          CALL DSYCON_ROOK( 'U', 1, A, 1, IP, -1.0D0, RCOND, W, IW, INFO)
          CALL CHKXER( 'DSYCON_ROOK', INFOT, NOUT, LERR, OK )
 *
+      ELSE IF( LSAMEN( 2, C2, 'SA' ) ) THEN
+*
+*        Test error exits of the routines that use factorization
+*        of a symmetric indefinite matrix with Aasen's algorithm.
+*
+*        DSYTRF_AA
+*
+         SRNAMT = 'DSYTRF_AA'
+         INFOT = 1
+         CALL DSYTRF_AA( '/', 0, A, 1, IP, W, 1, INFO )
+         CALL CHKXER( 'DSYTRF_AA', INFOT, NOUT, LERR, OK )
+         INFOT = 2
+         CALL DSYTRF_AA( 'U', -1, A, 1, IP, W, 1, INFO )
+         CALL CHKXER( 'DSYTRF_AA', INFOT, NOUT, LERR, OK )
+         INFOT = 4
+         CALL DSYTRF_AA( 'U', 2, A, 1, IP, W, 4, INFO )
+         CALL CHKXER( 'DSYTRF_AA', INFOT, NOUT, LERR, OK )
+*
+*        DSYTRS_AA
+*
+         SRNAMT = 'DSYTRS_AA'
+         INFOT = 1
+         CALL DSYTRS_AA( '/', 0, 0, A, 1, IP, B, 1, W, 1, INFO )
+         CALL CHKXER( 'DSYTRS_AA', INFOT, NOUT, LERR, OK )
+         INFOT = 2
+         CALL DSYTRS_AA( 'U', -1, 0, A, 1, IP, B, 1, W, 1, INFO )
+         CALL CHKXER( 'DSYTRS_AA', INFOT, NOUT, LERR, OK )
+         INFOT = 3
+         CALL DSYTRS_AA( 'U', 0, -1, A, 1, IP, B, 1, W, 1, INFO )
+         CALL CHKXER( 'DSYTRS_AA', INFOT, NOUT, LERR, OK )
+         INFOT = 5
+         CALL DSYTRS_AA( 'U', 2, 1, A, 1, IP, B, 2, W, 1, INFO )
+         CALL CHKXER( 'DSYTRS_AA', INFOT, NOUT, LERR, OK )
+         INFOT = 8
+         CALL DSYTRS_AA( 'U', 2, 1, A, 2, IP, B, 1, W, 1, INFO )
+         CALL CHKXER( 'DSYTRS_AA', INFOT, NOUT, LERR, OK )
       ELSE IF( LSAMEN( 2, C2, 'SP' ) ) THEN
 *
 *        Test error exits of the routines that use factorization
diff --git a/TESTING/LIN/schkaa.f b/TESTING/LIN/schkaa.f
index 0d53dedc..37984e14 100644
--- a/TESTING/LIN/schkaa.f
+++ b/TESTING/LIN/schkaa.f
@@ -158,10 +158,11 @@
       EXTERNAL           ALAREQ, SCHKEQ, SCHKGB, SCHKGE, SCHKGT, SCHKLQ,
      $                   SCHKPB, SCHKPO, SCHKPS, SCHKPP, SCHKPT, SCHKQ3,
      $                   SCHKQL, SCHKQR, SCHKRQ, SCHKSP, SCHKSY,
-     $                   SCHKSY_ROOK, SCHKTB, SCHKTP, SCHKTR, SCHKTZ,
-     $                   SDRVGB, SDRVGE, SDRVGT, SDRVLS, SDRVPB, SDRVPO,
-     $                   SDRVPP, SDRVPT, SDRVSP, SDRVSY, SDRVSY_ROOK,
-     $                   ILAVER, SCHKLQTP, SCHKQRT, SCHKQRTP
+     $                   SCHKSY_ROOK, SCHKSY_AA, SCHKTB, SCHKTP, SCHKTR,
+     $                   SCHKTZ, SDRVGB, SDRVGE, SDRVGT, SDRVLS, SDRVPB,
+     $                   SDRVPO, SDRVPP, SDRVPT, SDRVSP, SDRVSY,
+     $                   SDRVSY_ROOK, SDRVSY_AA, ILAVER, SCHKLQTP,
+     $                   SCHKQRT, SCHKQRTP
 *     ..
 *     .. Scalars in Common ..
       LOGICAL            LERR, OK
diff --git a/TESTING/LIN/schksy_aa.f b/TESTING/LIN/schksy_aa.f
index 75a9ebaa..0f668723 100644
--- a/TESTING/LIN/schksy_aa.f
+++ b/TESTING/LIN/schksy_aa.f
@@ -8,9 +8,9 @@
 *  Definition:
 *  ===========
 *
-*       SUBROUTINE SCHKSY_AAEN( DOTYPE, NN, NVAL, NNB, NBVAL, NNS, NSVAL,
-*                               THRESH, TSTERR, NMAX, A, AFAC, AINV, B, X,
-*                               XACT, WORK, RWORK, IWORK, NOUT )
+*       SUBROUTINE SCHKSY_AA( DOTYPE, NN, NVAL, NNB, NBVAL, NNS, NSVAL,
+*                             THRESH, TSTERR, NMAX, A, AFAC, AINV, B, X,
+*                             XACT, WORK, RWORK, IWORK, NOUT )
 *
 *       .. Scalar Arguments ..
 *       LOGICAL    TSTERR
@@ -168,8 +168,8 @@
 *
 *  =====================================================================
       SUBROUTINE SCHKSY_AA( DOTYPE, NN, NVAL, NNB, NBVAL, NNS, NSVAL,
-     $                         THRESH, TSTERR, NMAX, A, AFAC, AINV, B,
-     $                         X, XACT, WORK, RWORK, IWORK, NOUT )
+     $                      THRESH, TSTERR, NMAX, A, AFAC, AINV, B,
+     $                      X, XACT, WORK, RWORK, IWORK, NOUT )
 *
 *  -- LAPACK test routine (version 3.7.0) --
 *  -- LAPACK is a software package provided by Univ. of Tennessee,    --
diff --git a/TESTING/LIN/serrsy.f b/TESTING/LIN/serrsy.f
index 5326612a..8fd38687 100644
--- a/TESTING/LIN/serrsy.f
+++ b/TESTING/LIN/serrsy.f
@@ -88,8 +88,9 @@
 *     .. External Subroutines ..
       EXTERNAL           ALAESM, CHKXER, SSPCON, SSYCON_ROOK, SSPRFS,
      $                   SSPTRF, SSPTRI, SSPTRS, SSYCON, SSYRFS, SSYTF2,
-     $                   SSYTF2_ROOK, SSYTRF, SSYTRF_ROOK, SSYTRI,
-     $                   SSYTRI_ROOK, SSYTRI2, SSYTRS, SSYTRS_ROOK
+     $                   SSYTF2_ROOK, SSYTRF, SSYTRF_ROOK, SSYTRF_AA,
+     $                   SSYTRI, SSYTRI_ROOK, SSYTRI2, SSYTRS, 
+     $                   SSYTRS_ROOK, SSYTRS_AA
 *     ..
 *     .. Scalars in Common ..
       LOGICAL            LERR, OK
@@ -337,6 +338,42 @@
 *        of a symmetric indefinite packed matrix with patrial
 *        (Bunch-Kaufman) pivoting.
 *
+      ELSE IF( LSAMEN( 2, C2, 'SA' ) ) THEN
+*
+*        Test error exits of the routines that use factorization
+*        of a symmetric indefinite matrix with Aasen's algorithm.
+*
+*        SSYTRF_AA
+*
+         SRNAMT = 'SSYTRF_AA'
+         INFOT = 1
+         CALL SSYTRF_AA( '/', 0, A, 1, IP, W, 1, INFO )
+         CALL CHKXER( 'SSYTRF_AA', INFOT, NOUT, LERR, OK )
+         INFOT = 2
+         CALL SSYTRF_AA( 'U', -1, A, 1, IP, W, 1, INFO )
+         CALL CHKXER( 'SSYTRF_AA', INFOT, NOUT, LERR, OK )
+         INFOT = 4
+         CALL SSYTRF_AA( 'U', 2, A, 1, IP, W, 4, INFO )
+         CALL CHKXER( 'SSYTRF_AA', INFOT, NOUT, LERR, OK )
+*
+*        SSYTRS_AA
+*
+         SRNAMT = 'SSYTRS_AA'
+         INFOT = 1
+         CALL SSYTRS_AA( '/', 0, 0, A, 1, IP, B, 1, W, 1, INFO )
+         CALL CHKXER( 'SSYTRS_AA', INFOT, NOUT, LERR, OK )
+         INFOT = 2
+         CALL SSYTRS_AA( 'U', -1, 0, A, 1, IP, B, 1, W, 1, INFO )
+         CALL CHKXER( 'SSYTRS_AA', INFOT, NOUT, LERR, OK )
+         INFOT = 3
+         CALL SSYTRS_AA( 'U', 0, -1, A, 1, IP, B, 1, W, 1, INFO )
+         CALL CHKXER( 'SSYTRS_AA', INFOT, NOUT, LERR, OK )
+         INFOT = 5
+         CALL SSYTRS_AA( 'U', 2, 1, A, 1, IP, B, 2, W, 1, INFO )
+         CALL CHKXER( 'SSYTRS_AA', INFOT, NOUT, LERR, OK )
+         INFOT = 8
+         CALL SSYTRS_AA( 'U', 2, 1, A, 2, IP, B, 1, W, 1, INFO )
+         CALL CHKXER( 'SSYTRS_AA', INFOT, NOUT, LERR, OK )
       ELSE IF( LSAMEN( 2, C2, 'SP' ) ) THEN
 *
 *        SSPTRF
diff --git a/TESTING/LIN/zchkaa.f b/TESTING/LIN/zchkaa.f
index 283e3162..766f873f 100644
--- a/TESTING/LIN/zchkaa.f
+++ b/TESTING/LIN/zchkaa.f
@@ -160,13 +160,13 @@
 *     ..
 *     .. External Subroutines ..
       EXTERNAL           ALAREQ, ZCHKEQ, ZCHKGB, ZCHKGE, ZCHKGT, ZCHKHE,
-     $                   ZCHKHE_ROOK, ZCHKHP, ZCHKLQ, ZCHKPB, ZCHKPO,
-     $                   ZCHKPS, ZCHKPP, ZCHKPT, ZCHKQ3, ZCHKQL,
+     $                   ZCHKHE_ROOK, ZCHKHE_AA, ZCHKHP, ZCHKLQ, ZCHKPB,
+     $                   ZCHKPO, ZCHKPS, ZCHKPP, ZCHKPT, ZCHKQ3, ZCHKQL,
      $                   ZCHKQR, ZCHKRQ, ZCHKSP, ZCHKSY, ZCHKSY_ROOK,
      $                   ZCHKTB, ZCHKTP, ZCHKTR, ZCHKTZ, ZDRVGB, ZDRVGE,
-     $                   ZDRVGT, ZDRVHE, ZDRVHE_ROOK, ZDRVHP, ZDRVLS,
-     $                   ZDRVPB, ZDRVPO, ZDRVPP, ZDRVPT, ZDRVSP, ZDRVSY,
-     $                   ZDRVSY_ROOK, ILAVER, ZCHKQRT, ZCHKQRTP,
+     $                   ZDRVGT, ZDRVHE, ZDRVHE_ROOK, ZDRVHE_AA, ZDRVHP,
+     $                   ZDRVLS, ZDRVPB, ZDRVPO, ZDRVPP, ZDRVPT, ZDRVSP,
+     $                   ZDRVSY, ZDRVSY_ROOK, ILAVER, ZCHKQRT, ZCHKQRTP,
      $                   ZCHKLQT, ZCHKLQTP, ZCHKTSQR
 *     ..
 *     .. Scalars in Common ..
diff --git a/TESTING/LIN/zchkhe_aa.f b/TESTING/LIN/zchkhe_aa.f
index 26baae13..20d595be 100644
--- a/TESTING/LIN/zchkhe_aa.f
+++ b/TESTING/LIN/zchkhe_aa.f
@@ -9,8 +9,8 @@
 *  ===========
 *
 *       SUBROUTINE ZCHKHE_AA( DOTYPE, NN, NVAL, NNB, NBVAL, NNS, NSVAL,
-*                                THRESH, TSTERR, NMAX, A, AFAC, AINV, B, X,
-*                                XACT, WORK, RWORK, IWORK, NOUT )
+*                             THRESH, TSTERR, NMAX, A, AFAC, AINV, B, X,
+*                             XACT, WORK, RWORK, IWORK, NOUT )
 *
 *       .. Scalar Arguments ..
 *       LOGICAL            TSTERR
@@ -169,8 +169,8 @@
 *
 *  =====================================================================
       SUBROUTINE ZCHKHE_AA( DOTYPE, NN, NVAL, NNB, NBVAL, NNS, NSVAL,
-     $                         THRESH, TSTERR, NMAX, A, AFAC, AINV, B,
-     $                         X, XACT, WORK, RWORK, IWORK, NOUT )
+     $                      THRESH, TSTERR, NMAX, A, AFAC, AINV, B,
+     $                      X, XACT, WORK, RWORK, IWORK, NOUT )
 *
 *  -- LAPACK test routine (version 3.7.0) --
 *  -- LAPACK is a software package provided by Univ. of Tennessee,    --
diff --git a/TESTING/LIN/zerrhe.f b/TESTING/LIN/zerrhe.f
index 60066288..47b64ae0 100644
--- a/TESTING/LIN/zerrhe.f
+++ b/TESTING/LIN/zerrhe.f
@@ -90,9 +90,9 @@
 *     .. External Subroutines ..
       EXTERNAL           ALAESM, CHKXER, ZHECON, ZHECON_ROOK, ZHERFS,
      $                   ZHETF2, ZHETF2_ROOK, ZHETRF, ZHETRF_ROOK,
-     $                   ZHETRI, ZHETRI_ROOK, ZHETRI2, ZHETRS,
-     $                   ZHETRS_ROOK, ZHPCON, ZHPRFS, ZHPTRF, ZHPTRI,
-     $                   ZHPTRS
+     $                   ZHETRF_AA, ZHETRI, ZHETRI_ROOK, ZHETRI2, 
+     $                   ZHETRS, ZHETRS_ROOK, ZHETRS_AA, ZHPCON, ZHPRFS,
+     $                   ZHPTRF, ZHPTRI, ZHPTRS
 *     ..
 *     .. Scalars in Common ..
       LOGICAL            LERR, OK
@@ -337,6 +337,43 @@
          CALL CHKXER( 'ZHECON_ROOK', INFOT, NOUT, LERR, OK )
 *
 *        Test error exits of the routines that use factorization
+*        of a Hermitian indefinite matrix with Aasen's algorithm.
+*
+      ELSE IF( LSAMEN( 2, C2, 'HA' ) ) THEN
+*
+*        ZHETRF_AA
+*
+         SRNAMT = 'ZHETRF_AA'
+         INFOT = 1
+         CALL ZHETRF_AA( '/', 0, A, 1, IP, W, 1, INFO )
+         CALL CHKXER( 'ZHETRF_AA', INFOT, NOUT, LERR, OK )
+         INFOT = 2
+         CALL ZHETRF_AA( 'U', -1, A, 1, IP, W, 1, INFO )
+         CALL CHKXER( 'ZHETRF_AA', INFOT, NOUT, LERR, OK )
+         INFOT = 4
+         CALL ZHETRF_AA( 'U', 2, A, 1, IP, W, 4, INFO )
+         CALL CHKXER( 'ZHETRF_AA', INFOT, NOUT, LERR, OK )
+*
+*        ZHETRS_AA
+*
+         SRNAMT = 'ZHETRS_AA'
+         INFOT = 1
+         CALL ZHETRS_AA( '/', 0, 0, A, 1, IP, B, 1, W, 1, INFO )
+         CALL CHKXER( 'ZHETRS_AA', INFOT, NOUT, LERR, OK )
+         INFOT = 2
+         CALL ZHETRS_AA( 'U', -1, 0, A, 1, IP, B, 1, W, 1, INFO )
+         CALL CHKXER( 'ZHETRS_AA', INFOT, NOUT, LERR, OK )
+         INFOT = 3
+         CALL ZHETRS_AA( 'U', 0, -1, A, 1, IP, B, 1, W, 1, INFO )
+         CALL CHKXER( 'ZHETRS_AA', INFOT, NOUT, LERR, OK )
+         INFOT = 5
+         CALL ZHETRS_AA( 'U', 2, 1, A, 1, IP, B, 2, W, 1, INFO )
+         CALL CHKXER( 'ZHETRS_AA', INFOT, NOUT, LERR, OK )
+         INFOT = 8
+         CALL ZHETRS_AA( 'U', 2, 1, A, 2, IP, B, 1, W, 1, INFO )
+         CALL CHKXER( 'ZHETRS_AA', INFOT, NOUT, LERR, OK )
+*
+*        Test error exits of the routines that use factorization
 *        of a Hermitian indefinite packed matrix with patrial
 *        (Bunch-Kaufman) diagonal pivoting method.
 *