derrec.f

Go to the documentation of this file.

*> \brief \b DERREC
*
*  =========== DOCUMENTATION ===========
*
* Online html documentation available at
*            http://www.netlib.org/lapack/explore-html/
*
*  Definition:
*  ===========
*
*       SUBROUTINE DERREC( PATH, NUNIT )
*
*       .. Scalar Arguments ..
*       CHARACTER*3        PATH
*       INTEGER            NUNIT
*       ..
*
*
*> \par Purpose:
*  =============
*>
*> \verbatim
*>
*> DERREC tests the error exits for the routines for eigen- condition
*> estimation for DOUBLE PRECISION matrices:
*>    DTRSYL, DTREXC, DTRSNA and DTRSEN.
*> \endverbatim
*
*  Arguments:
*  ==========
*
*> \param[in] PATH
*> \verbatim
*>          PATH is CHARACTER*3
*>          The LAPACK path name for the routines to be tested.
*> \endverbatim
*>
*> \param[in] NUNIT
*> \verbatim
*>          NUNIT is INTEGER
*>          The unit number for output.
*> \endverbatim
*
*  Authors:
*  ========
*
*> \author Univ. of Tennessee
*> \author Univ. of California Berkeley
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
*> \ingroup double_eig
*
*  =====================================================================
      SUBROUTINE derrec( PATH, NUNIT )
*
*  -- LAPACK test routine --
*  -- LAPACK is a software package provided by Univ. of Tennessee,    --
*  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
*
*     .. Scalar Arguments ..
      CHARACTER*3        PATH
      INTEGER            NUNIT
*     ..
*
*  =====================================================================
*
*     .. Parameters ..
      INTEGER            NMAX
      DOUBLE PRECISION   ONE, ZERO
      parameter( nmax = 4, one = 1.0d0, zero = 0.0d0 )
*     ..
*     .. Local Scalars ..
      INTEGER            I, IFST, ILST, INFO, J, M, NT
      DOUBLE PRECISION   SCALE
*     ..
*     .. Local Arrays ..
      LOGICAL            SEL( NMAX )
      INTEGER            IWORK( NMAX )
      DOUBLE PRECISION   A( NMAX, NMAX ), B( NMAX, NMAX ),
     $                   C( NMAX, NMAX ), S( NMAX ), SEP( NMAX ),
     $                   WI( NMAX ), WORK( NMAX ), WR( NMAX )
*     ..
*     .. External Subroutines ..
      EXTERNAL           chkxer, dtrexc, dtrsen, dtrsna, dtrsyl
*     ..
*     .. Scalars in Common ..
      LOGICAL            LERR, OK
      CHARACTER*32       SRNAMT
      INTEGER            INFOT, NOUT
*     ..
*     .. Common blocks ..
      COMMON             / infoc / infot, nout, ok, lerr
      COMMON             / srnamc / srnamt
*     ..
*     .. Executable Statements ..
*
      nout = nunit
      ok = .true.
      nt = 0
*
*     Initialize A, B and SEL
*
      DO 20 j = 1, nmax
         DO 10 i = 1, nmax
            a( i, j ) = zero
            b( i, j ) = zero
   10    CONTINUE
   20 CONTINUE
      DO 30 i = 1, nmax
         a( i, i ) = one
         sel( i ) = .true.
   30 CONTINUE
*
*     Test DTRSYL
*
      srnamt = 'DTRSYL'
      infot = 1
      CALL dtrsyl( 'X', 'N', 1, 0, 0, a, 1, b, 1, c, 1, scale, info )
      CALL chkxer( 'DTRSYL', infot, nout, lerr, ok )
      infot = 2
      CALL dtrsyl( 'N', 'X', 1, 0, 0, a, 1, b, 1, c, 1, scale, info )
      CALL chkxer( 'DTRSYL', infot, nout, lerr, ok )
      infot = 3
      CALL dtrsyl( 'N', 'N', 0, 0, 0, a, 1, b, 1, c, 1, scale, info )
      CALL chkxer( 'DTRSYL', infot, nout, lerr, ok )
      infot = 4
      CALL dtrsyl( 'N', 'N', 1, -1, 0, a, 1, b, 1, c, 1, scale, info )
      CALL chkxer( 'DTRSYL', infot, nout, lerr, ok )
      infot = 5
      CALL dtrsyl( 'N', 'N', 1, 0, -1, a, 1, b, 1, c, 1, scale, info )
      CALL chkxer( 'DTRSYL', infot, nout, lerr, ok )
      infot = 7
      CALL dtrsyl( 'N', 'N', 1, 2, 0, a, 1, b, 1, c, 2, scale, info )
      CALL chkxer( 'DTRSYL', infot, nout, lerr, ok )
      infot = 9
      CALL dtrsyl( 'N', 'N', 1, 0, 2, a, 1, b, 1, c, 1, scale, info )
      CALL chkxer( 'DTRSYL', infot, nout, lerr, ok )
      infot = 11
      CALL dtrsyl( 'N', 'N', 1, 2, 0, a, 2, b, 1, c, 1, scale, info )
      CALL chkxer( 'DTRSYL', infot, nout, lerr, ok )
      nt = nt + 8
*
*     Test DTREXC
*
      srnamt = 'DTREXC'
      ifst = 1
      ilst = 1
      infot = 1
      CALL dtrexc( 'X', 1, a, 1, b, 1, ifst, ilst, work, info )
      CALL chkxer( 'DTREXC', infot, nout, lerr, ok )
      infot = 2
      CALL dtrexc( 'N', -1, a, 1, b, 1, ifst, ilst, work, info )
      CALL chkxer( 'DTREXC', infot, nout, lerr, ok )
      infot = 4
      ilst = 2
      CALL dtrexc( 'N', 2, a, 1, b, 1, ifst, ilst, work, info )
      CALL chkxer( 'DTREXC', infot, nout, lerr, ok )
      infot = 6
      CALL dtrexc( 'V', 2, a, 2, b, 1, ifst, ilst, work, info )
      CALL chkxer( 'DTREXC', infot, nout, lerr, ok )
      infot = 7
      ifst = 0
      ilst = 1
      CALL dtrexc( 'V', 1, a, 1, b, 1, ifst, ilst, work, info )
      CALL chkxer( 'DTREXC', infot, nout, lerr, ok )
      infot = 7
      ifst = 2
      CALL dtrexc( 'V', 1, a, 1, b, 1, ifst, ilst, work, info )
      CALL chkxer( 'DTREXC', infot, nout, lerr, ok )
      infot = 8
      ifst = 1
      ilst = 0
      CALL dtrexc( 'V', 1, a, 1, b, 1, ifst, ilst, work, info )
      CALL chkxer( 'DTREXC', infot, nout, lerr, ok )
      infot = 8
      ilst = 2
      CALL dtrexc( 'V', 1, a, 1, b, 1, ifst, ilst, work, info )
      CALL chkxer( 'DTREXC', infot, nout, lerr, ok )
      nt = nt + 8
*
*     Test DTRSNA
*
      srnamt = 'DTRSNA'
      infot = 1
      CALL dtrsna( 'X', 'A', sel, 0, a, 1, b, 1, c, 1, s, sep, 1, m,
     $             work, 1, iwork, info )
      CALL chkxer( 'DTRSNA', infot, nout, lerr, ok )
      infot = 2
      CALL dtrsna( 'B', 'X', sel, 0, a, 1, b, 1, c, 1, s, sep, 1, m,
     $             work, 1, iwork, info )
      CALL chkxer( 'DTRSNA', infot, nout, lerr, ok )
      infot = 4
      CALL dtrsna( 'B', 'A', sel, -1, a, 1, b, 1, c, 1, s, sep, 1, m,
     $             work, 1, iwork, info )
      CALL chkxer( 'DTRSNA', infot, nout, lerr, ok )
      infot = 6
      CALL dtrsna( 'V', 'A', sel, 2, a, 1, b, 1, c, 1, s, sep, 2, m,
     $             work, 2, iwork, info )
      CALL chkxer( 'DTRSNA', infot, nout, lerr, ok )
      infot = 8
      CALL dtrsna( 'B', 'A', sel, 2, a, 2, b, 1, c, 2, s, sep, 2, m,
     $             work, 2, iwork, info )
      CALL chkxer( 'DTRSNA', infot, nout, lerr, ok )
      infot = 10
      CALL dtrsna( 'B', 'A', sel, 2, a, 2, b, 2, c, 1, s, sep, 2, m,
     $             work, 2, iwork, info )
      CALL chkxer( 'DTRSNA', infot, nout, lerr, ok )
      infot = 13
      CALL dtrsna( 'B', 'A', sel, 1, a, 1, b, 1, c, 1, s, sep, 0, m,
     $             work, 1, iwork, info )
      CALL chkxer( 'DTRSNA', infot, nout, lerr, ok )
      infot = 13
      CALL dtrsna( 'B', 'S', sel, 2, a, 2, b, 2, c, 2, s, sep, 1, m,
     $             work, 2, iwork, info )
      CALL chkxer( 'DTRSNA', infot, nout, lerr, ok )
      infot = 16
      CALL dtrsna( 'B', 'A', sel, 2, a, 2, b, 2, c, 2, s, sep, 2, m,
     $             work, 1, iwork, info )
      CALL chkxer( 'DTRSNA', infot, nout, lerr, ok )
      nt = nt + 9
*
*     Test DTRSEN
*
      sel( 1 ) = .false.
      srnamt = 'DTRSEN'
      infot = 1
      CALL dtrsen( 'X', 'N', sel, 0, a, 1, b, 1, wr, wi, m, s( 1 ),
     $             sep( 1 ), work, 1, iwork, 1, info )
      CALL chkxer( 'DTRSEN', infot, nout, lerr, ok )
      infot = 2
      CALL dtrsen( 'N', 'X', sel, 0, a, 1, b, 1, wr, wi, m, s( 1 ),
     $             sep( 1 ), work, 1, iwork, 1, info )
      CALL chkxer( 'DTRSEN', infot, nout, lerr, ok )
      infot = 4
      CALL dtrsen( 'N', 'N', sel, -1, a, 1, b, 1, wr, wi, m, s( 1 ),
     $             sep( 1 ), work, 1, iwork, 1, info )
      CALL chkxer( 'DTRSEN', infot, nout, lerr, ok )
      infot = 6
      CALL dtrsen( 'N', 'N', sel, 2, a, 1, b, 1, wr, wi, m, s( 1 ),
     $             sep( 1 ), work, 2, iwork, 1, info )
      CALL chkxer( 'DTRSEN', infot, nout, lerr, ok )
      infot = 8
      CALL dtrsen( 'N', 'V', sel, 2, a, 2, b, 1, wr, wi, m, s( 1 ),
     $             sep( 1 ), work, 1, iwork, 1, info )
      CALL chkxer( 'DTRSEN', infot, nout, lerr, ok )
      infot = 15
      CALL dtrsen( 'N', 'V', sel, 2, a, 2, b, 2, wr, wi, m, s( 1 ),
     $             sep( 1 ), work, 0, iwork, 1, info )
      CALL chkxer( 'DTRSEN', infot, nout, lerr, ok )
      infot = 15
      CALL dtrsen( 'E', 'V', sel, 3, a, 3, b, 3, wr, wi, m, s( 1 ),
     $             sep( 1 ), work, 1, iwork, 1, info )
      CALL chkxer( 'DTRSEN', infot, nout, lerr, ok )
      infot = 15
      CALL dtrsen( 'V', 'V', sel, 3, a, 3, b, 3, wr, wi, m, s( 1 ),
     $             sep( 1 ), work, 3, iwork, 2, info )
      CALL chkxer( 'DTRSEN', infot, nout, lerr, ok )
      infot = 17
      CALL dtrsen( 'E', 'V', sel, 2, a, 2, b, 2, wr, wi, m, s( 1 ),
     $             sep( 1 ), work, 1, iwork, 0, info )
      CALL chkxer( 'DTRSEN', infot, nout, lerr, ok )
      infot = 17
      CALL dtrsen( 'V', 'V', sel, 3, a, 3, b, 3, wr, wi, m, s( 1 ),
     $             sep( 1 ), work, 4, iwork, 1, info )
      CALL chkxer( 'DTRSEN', infot, nout, lerr, ok )
      nt = nt + 10
*
*     Print a summary line.
*
      IF( ok ) THEN
         WRITE( nout, fmt = 9999 )path, nt
      ELSE
         WRITE( nout, fmt = 9998 )path
      END IF
*
      RETURN
 9999 FORMAT( 1x, a3, ' routines passed the tests of the error exits (',
     $      i3, ' tests done)' )
 9998 FORMAT( ' *** ', a3, ' routines failed the tests of the error ex',
     $      'its ***' )
*
*     End of DERREC
*
      END